| 查看: 1647 | 回复: 6 | ||
| 【悬赏金币】回答本帖问题，作者风起天阑1将赠送您 50 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
风起天阑1木虫 (正式写手)
|
[求助]
不可压稳态线接触滑块程序改编成MATLAB程序已有1人参与
|
|
|
本人在将黄平老师所著的润滑数值计算方法的第二章不可压稳态线接触滑块程序改编成MATLAB程序时压力的计算出了问题，但是自己找不到解决办法，请大佬解惑，谢谢。程序如下： X = zeros(121,1); % 坐标节点数组 H = zeros(121,1); % 膜厚节点数组 P = zeros(121,1); % 压力节点数组 H1=1.0; % 量纲化最大膜厚 H2=0.5; % 量纲化最小膜厚 U=1; % 量纲化速度 EDA=0.02;% 润滑油粘度 AL=0.01; % 滑块长度 KG=1;% 等于1时为直线滑块，2时为曲线滑块 N=121;% 节点取值 if KG==2 X1=-1.0; X2=1.0; else X1=0.0; X2=1.0; end DX=1/(N-1); for i=1:N if KG==1 % 直线滑块膜厚计算 X(i)=X1-(i-1.0)*DX*(X1-X2); H(i)=H1/H2-(H1/H2-1.0)*X(i); else % 曲线滑块膜厚计算 X(i)=X1-(i-1.0)*DX*(X1-X2); H(i)=1.0+(H1/H2-1.0)*X(i)^2; end end for i=2:N-1 P(i)=0.5;% 压力各点赋初值 end P(1)=0.0;% 压力起点初始化 P(N)=0.0;% 压力终点初始化 IK=0; ERO=1; while ERO>1e-7 C1=0.0; ALOAD=0.0; for i=2:N-1 A1=(0.5*(H(i+1)+H(i)))^3; A2=(0.5*(H(i)+H(i-1)))^3; PD=P(i); P(i)=(-0.5*DX*(H(i+1)-H(i-1))+A1*P(i+1)+A2*P(i-1))/(A1+A2); P2(i)=0.3*PD+0.7*P(i); if P2(i)<0 P2(i)=0; end C1=C1+abs(P2(i)-PD); ALOAD=ALOAD+P2(i); end ERO=C1/ALOAD;% 相对误差，用于判断收敛情况 IK=IK+1; end X0=0; for i=1:N X0=X0+P(i)*X(i); end X0=X0*AL; ALOAD=ALOAD*DX*AL*6.0*U*EDA*AL/(H2^2);% 承载能力计算 if KG==1 for i=1:N P0=-(-1.0/(H(i)*H2)+H1*H2/(H1+H2)/(H2*H(i))^2+1.0/(H1+H2))/(H1/H2-1.0)*H2; end end  @wuming524 |
回复此楼» 猜你喜欢
复试调剂
已经有3人回复
-
欢迎211本科，线上复试，低分同学报考课题组，A区一本，交叉学科
已经有8人回复
-
322调剂
已经有3人回复
-
学科化学352分求调剂
已经有4人回复
-
309求调剂
已经有6人回复
-
286生物学求调剂
已经有3人回复
-
269求调剂
已经有5人回复
-
化学工程 288分求调剂
已经有10人回复
-
调剂
已经有5人回复
-
材料与化工085600求调剂
已经有8人回复
风起天阑1
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2563.1
- 散金: 122
- 红花: 1
- 帖子: 821
- 在线: 79.1小时
- 虫号: 6327111
- 注册: 2017-04-16
- 性别: GG
- 专业: 机械结构强度学
2楼2018-09-18 07:58:26
风起天阑1
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2563.1
- 散金: 122
- 红花: 1
- 帖子: 821
- 在线: 79.1小时
- 虫号: 6327111
- 注册: 2017-04-16
- 性别: GG
- 专业: 机械结构强度学
送红花一朵|
下面图片就是最后数据图的结果，楼主现在只是初学者，只能拾人牙慧，简单的编程，对利用矩阵等方法还不是太过熟悉，还希望大佬不吝赐教，能否帮忙利用矩阵的方法修改一下 。谢谢    发自小木虫Android客户端 |
4楼2018-09-18 14:10:21
风起天阑1
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2563.1
- 散金: 122
- 红花: 1
- 帖子: 821
- 在线: 79.1小时
- 虫号: 6327111
- 注册: 2017-04-16
- 性别: GG
- 专业: 机械结构强度学
6楼2018-09-18 18:59:41
30
