24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Dr_Kevin

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1080.9
散金: 118
红花: 1
帖子: 163
在线: 57.7小时
虫号: 2458958
注册: 2013-05-11
专业: 机械动力学

[求助] 小弟求滑动轴承非线性油膜力matlab程序，最好是差分法的，谢谢已有2人参与

小弟求滑动轴承非线性油膜力matlab程序，最好是差分法的，谢谢

回复此楼

» 收录本帖的淘帖专辑推荐

摩擦疲劳振动与可靠性

滑动轴承及铸造专题

» 猜你喜欢

1楼 2015-11-02 11:13:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Fairyflower

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1068.7
散金: 38
红花: 1
帖子: 271
在线: 8小时
虫号: 5188909
注册: 2016-11-01
专业: 机械工程

【答案】应助回帖

引用回帖:

2楼: Originally posted by 雪之星愿 at 2016-01-14 15:28:03
% This function is used to calculate the pressure distribution of journal bearing without taking the deformation into account
delta_m=2*pi/m; %沿圆周方向均匀划分m格 ...

你好，能不能给一下这个程序模型的出处，谢谢！

赞一下

回复此楼

5楼2016-11-01 14:31:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 6 个回答

雪之星愿

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 23.1
红花: 1
帖子: 18
在线: 9.2小时
虫号: 4035217
注册: 2015-08-24
性别: MM
专业: 无机化学

【答案】应助回帖

% This function is used to calculate the pressure distribution of journal bearing without taking the deformation into account
delta_m=2*pi/m;                                  %沿圆周方向均匀划分m格
delta_n=2/n;                                     %沿轴向均匀划分n格
for j=1:n+1                                        %对各节点赋压力初值为0
for i=1:m+1
      P(i,j)=0;
      Q(i,j)=0;
end
end
S=0;
T=0;
for j=2:n                                  %第一次计算各内部节点的压力值
for i=2:m
   A=(1+e*cos((i+1/2-1)*delta_m))^3;
B=(1+e*cos((i-1/2-1)*delta_m))^3;
C=((1/ratio)*delta_m/delta_n)^2*(1+e*cos((i-1)*delta_m))^3;
D=C;
E=A+B+C+D;
F=delta_m*e*(cos((i+1/2-1)*delta_m)-cos((i-1/2-1)*delta_m));
P(i,j)=(A*P(i+1,j)+B*P(i-1,j)+C*P(i,j+1)+D*P(i,j-1)-F)/E;
if (P(i,j)<=0)                                  %引入雷诺边界条件
         for k=i:m
            P(k,j)=0;
            S=S+abs(Q(k,j));
            T=T+abs(P(i,j));
            Q(k,j)=P(k,j);
         end
         break;
end
S=S+abs(P(i,j)-Q(i,j));
T=T+abs(P(i,j));
Q(i,j)=P(i,j);
end
end
while(S/T>0.001)                         %循环计算各内部节点的压力值
S=0;
T=0;
for j=2:n
for i=2:m
A=(1+e*cos((i+1/2-1)*delta_m))^3;
B=(1+e*cos((i-1/2-1)*delta_m))^3;
C=((1/ratio)*delta_m/delta_n)^2*(1+e*cos((i-1)*delta_m))^3;
D=C;
E=A+B+C+D;
F=delta_m*e*(cos((i+1/2-1)*delta_m)-cos((i-1/2-1)*delta_m));
P(i,j)=(A*P(i+1,j)+B*P(i-1,j)+C*P(i,j+1)+D*P(i,j-1)-F)/E;
if (P(i,j)<=0)
            for k=i:m
                  P(k,j)=0;
                  S=S+abs(Q(k,j));
                  T=T+abs(P(k,j));
                  Q(k,j)=P(k,j);
            end
            break;
end
         S=S+abs(P(i,j)-Q(i,j));
         T=T+abs(P(i,j));
         Q(i,j)=P(i,j);
      end
end
end
%P;                   %此部分是当本函数单独使用时，用于绘制离散压力分布
%for i=1:m+1
% x(i)=(i-1)*delta_m;
%end
%for j=1:n+1
% y(j)=(j-1)*delta_n;
%end
%figure

赞一下(2人)

回复此楼

2楼2016-01-14 15:28:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖