| 查看: 1682 | 回复: 5 | ||||
Dr_Kevin木虫 (小有名气)
|
[求助]
小弟求滑动轴承非线性油膜力matlab程序，最好是差分法的，谢谢 已有2人参与
|
| 小弟求滑动轴承非线性油膜力matlab程序，最好是差分法的，谢谢 |
回复此楼» 收录本帖的淘帖专辑推荐
 摩擦疲劳振动与可靠性 | 滑动轴承及铸造专题 |
» 猜你喜欢
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”学硕与专硕有缺额
已经有0人回复
-
重庆长江师范学院材料工程接收材料类和机械类考生调剂硕士研究生
已经有26人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有179人回复
-
重庆长江师范学院材料工程接收材料类考生调剂硕士研究生
已经有0人回复
-
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”方向有调剂名额
已经有0人回复
-
天津大学化学工程专业专业课
已经有0人回复
-
二氧化铂的XRD和XPS 检测结果完全对不上，这是为什么？？？？
已经有0人回复
-
辽宁工程技术大学材料科学与工程、储能科学与工程、材料与化工专业接收调剂，名额充足
已经有0人回复
-
332，085601求调剂
已经有14人回复
-
英语一来了就有
已经有0人回复
-
福建理工大学招收英语一 数学一或者二
已经有0人回复
【答案】应助回帖
|
% 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楼2016-01-14 15:28:03
Dr_Kevin
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1080.9
- 散金: 118
- 红花: 1
- 帖子: 163
- 在线: 57.7小时
- 虫号: 2458958
- 注册: 2013-05-11
- 专业: 机械动力学
3楼2016-01-15 09:54:04
Dr_Kevin
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1080.9
- 散金: 118
- 红花: 1
- 帖子: 163
- 在线: 57.7小时
- 虫号: 2458958
- 注册: 2013-05-11
- 专业: 机械动力学
4楼2016-01-15 10:23:49
Fairyflower
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1068.7
- 散金: 38
- 红花: 1
- 帖子: 271
- 在线: 8小时
- 虫号: 5188909
- 注册: 2016-11-01
- 专业: 机械工程
5楼2016-11-01 14:31:31
6楼2018-07-19 23:10:31
5