| 查看: 1553 | 回复: 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 |
回复此楼» 猜你喜欢
5-氯茚酮
已经有6人回复
-
换导师学院里面不同意咋办，原导师现导师都同意了
已经有18人回复
-
投稿返修中，增加了作者，论文被拒了，还有办法补救嘛？
已经有18人回复
-
请教大家一个问题
已经有5人回复
-
求个博导把我捞走，一篇JMCA，另有一篇AFM共一
已经有7人回复
-
搞科研是不是很锻炼身体？高血压，痛风是小病
已经有8人回复
-
76批博后面上中标，小满贯！
已经有17人回复
-
咨询一下，去江西高校工作和发展前景怎么样
已经有13人回复
-
哪里招科研助理
已经有3人回复
-
长江学者评审开始了，有评审专家爆料被地毯式打招呼
已经有14人回复
风起天阑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
10