24小时热门版块排行榜

>论坛更新日志 (3710)
>虫友互识 (515)
>文献求助 (450)
>导师招生 (338)
>考博 (181)
>休闲灌水 (179)
>论文投稿 (169)
>基金申请 (115)
>硕博家园 (106)
>招聘信息布告栏 (82)
>博后之家 (79)
>教师之家 (63)
>论文道贺祈福 (58)
>找工作 (38)
>公派出国 (37)
>考研 (37)

返回列表

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

Dr_Kevin

木虫 (小有名气)

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

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

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

回复此楼

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

摩擦疲劳振动与可靠性

滑动轴承及铸造专题

» 猜你喜欢

解决渗碳体倾斜角度问题已经有0人回复
北京科技大学/李亚庚教授招收2026年秋季入学博士生已经有11人回复
机械工程论文润色/翻译怎么收费? 已经有114人回复
省级国家级人才计划申报已经有0人回复
探究TGF-β在癌症免疫调控中的作用机制|肿瘤已经有2人回复
探究TGF-β在癌症免疫调控中的作用机制|肿瘤已经有10人回复
所以当初说在欧美生活&#144308;用讲人情世故的人怎么生活的已经有12人回复
求助SambVca最新网址或等效工具已经有0人回复
求助Journal of Alloys and Compounds剘刊2014年612卷的封皮、目录、版权页或封底已经有5人回复
英国伯明翰大学招收计算流体力学博士研究生已经有0人回复
英国伯明翰大学招收计算流体力学和气动声学博士研究生已经有0人回复

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

已阅回复此楼关注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的回帖

查看全部 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的回帖