24小时热门版块排行榜

返回列表

Dr_Kevin

木虫 (小有名气)

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

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

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

回复此楼

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

摩擦疲劳振动与可靠性

滑动轴承及铸造专题

» 猜你喜欢

[2027 考博] 上海・｜激光焊接 / 增材制造 / 机器人智能装备招博士已经有112人回复
英国埃克塞特大学（QS 136）招CSC博士生，3D打印增材制造方向已经有4人回复
金属材料论文润色/翻译怎么收费? 已经有230人回复
【CSC博士招生】英国拉夫堡大学｜机器人+3D视觉/结构光测量+AI｜2027入学已经有1人回复
振动测试台、盐雾试验机、氙灯耐候试验箱怎么选？环境可靠性试验设备厂家推荐指南已经有0人回复
工业高温烘烤箱哪家好？源头厂家设备性价比高已经有0人回复
高低温冲击试验箱哪个厂家的好？十大品牌推荐及质量点评已经有0人回复
无损检测校准机构怎么选？美标报告、超声涡流设备校准服务商已经有0人回复
加拿大Western University招收2027年博士生｜增材制造、复合材料与可持续材料已经有0人回复

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的回帖