24小时热门版块排行榜

返回列表

Dr_Kevin

木虫 (小有名气)

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

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

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

回复此楼

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

摩擦疲劳振动与可靠性

滑动轴承及铸造专题

» 猜你喜欢

请问45岁后（高校在编教师）还能申请在职博士吗？已经有27人回复
工材面上未中，大佬们帮忙看看这些意见评分如何？给小弟指明方向已经有24人回复
金属材料论文润色/翻译怎么收费? 已经有141人回复
疲劳寿命预测求收留已经有3人回复
SCI论文投稿with editor十天正常吗？已经有2人回复
26年秋季金属方向博士申请-目前研三在读-发表SCI 1区和2区学术论文共两篇已经有1人回复
审稿日期短时间内频繁更换，是被拒的节奏吗？已经有10人回复
如何产生Frank分位错已经有2人回复
请大家帮忙看下，这样的审稿状态是否凶多吉少？已经有12人回复

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

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

雪之星愿

新虫 (初入文坛)

应助: 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的回帖

Dr_Kevin

木虫 (小有名气)

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

引用回帖:

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格 ...

谢谢~~只是。。。您这个程序中的F表达式对么。。。

赞一下

回复此楼

3楼2016-01-15 09:54:04

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