应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 不可压稳态线接触滑块程序改编成MATLAB程序
51/1返回列表
查看: 1552  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者风起天阑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
回复此楼

» 猜你喜欢
知行合一，心外无物
1楼2018-09-17 14:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
4楼: Originally posted by 风起天阑1 at 2018-09-18 14:10:21
下面图片就是最后数据图的结果，楼主现在只是初学者，只能拾人牙慧，简单的编程，对利用矩阵等方法还不是太过熟悉，还希望大佬不吝赐教，能否帮忙利用矩阵的方法修改一下 。谢谢


...

改了改部分代码，供参考：
CODE:
clear

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);
if KG == 1 % 直线滑块膜厚计算
    X = X1 : ( X2 - X1 ) / ( N - 1 ) : X2;
    H = H1 / H2 - ( H1 / H2 - 1.0 ) * X;
else     % 曲线滑块膜厚计算
    X = X1 : ( X2 - X1 ) /  ( N - 1 ) : X2;
    H = 1.0 + ( H1 / H2 - 1.0 ) * X .^ 2;
end
P( 1 ) = 0.0; % 压力起点初始化
P( 2 : N -1 ) = 0.5; % 压力各点赋初值
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 = sum( P .* X' ) * AL;
ALOAD=ALOAD*DX*AL*6.0*U*EDA*AL/(H2^2);% 承载能力计算
P0 = zeros(121,1);
if KG==1
     for i=1:N
       P0( i )=-(-1.0/(H(i)*H2)+H1*H2/(H1+H2)/(H2*H(i))^2+1.0/(H1+H2))/(H1/H2-1.0)*H2;
     end
end
figure
subplot( 2, 2, 1 )
plot( X, P0, '--b' )
subplot( 2, 2, 3 )
plot( X, H, 'r' )
subplot( 2, 2, [ 2, 4 ] )
plot( X, P0, '--b' )
hold on
plot( X, H, 'r' )

结果和书上有出入，压力分布的值有数量级的差异，还需仔细检查代码。
不可压稳态线接触滑块程序改编成MATLAB程序
计算结果.png
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

风起天阑1
明德厚学、求是创新
5楼2018-09-18 15:27:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

风起天阑1

木虫 (正式写手)
自己顶一下，大佬帮帮忙

发自小木虫Android客户端
回复此楼
知行合一，心外无物
2楼2018-09-18 07:58:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主
正确的结果是什么？
运行了一下程序，可以运行完。没数学模型，也没有书，不知道楼主的具体问题是什么。

不过程序有优化的空间，建议尽量采用矩阵运算，减少for循环，一方面代码简洁，可读性高，再就是运行效率略高。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

风起天阑1
明德厚学、求是创新
3楼2018-09-18 13:54:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风起天阑1

木虫 (正式写手)
送红花一朵
引用回帖:
3楼: Originally posted by xiegangmai at 2018-09-18 13:54:32
正确的结果是什么？
运行了一下程序，可以运行完。没数学模型，也没有书，不知道楼主的具体问题是什么。

不过程序有优化的空间，建议尽量采用矩阵运算，减少for循环，一方面代码简洁，可读性高，再就是运行效率 ...

下面图片就是最后数据图的结果，楼主现在只是初学者，只能拾人牙慧，简单的编程，对利用矩阵等方法还不是太过熟悉，还希望大佬不吝赐教，能否帮忙利用矩阵的方法修改一下 。谢谢
不可压稳态线接触滑块程序改编成MATLAB程序-1


不可压稳态线接触滑块程序改编成MATLAB程序-2



发自小木虫Android客户端
赞一下 回复此楼
知行合一，心外无物
4楼2018-09-18 14:10:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[工艺技术] 5-氯茚酮 +4 红色阡陌 2024-11-20 6/300 2024-11-23 23:47 by linweikang
[论文投稿] 投稿返修中，增加了作者，论文被拒了，还有办法补救嘛？ +9 mollyzhang_2003 2024-11-22 18/900 2024-11-23 22:02 by mollyzhang_2003
[基金申请] 请教大家一个问题 +4 愿每一个人都平� 2024-11-22 5/250 2024-11-23 21:04 by RRRKKK
[考博] 求个博导把我捞走，一篇JMCA，另有一篇AFM共一 +5 骑萝卜钓鱼 2024-11-21 7/350 2024-11-23 21:03 by Zixiao_Wang
[博后之家] 76批博后面上中标，小满贯！ (金币+5) +17 UnUnc_Wang 2024-11-19 17/850 2024-11-23 19:36 by Bob_robert
[教师之家] 咨询一下，去江西高校工作和发展前景怎么样 +11 akslis2024 2024-11-17 13/650 2024-11-23 18:37 by luokereng
[论文投稿] 这是必须要润色吗 +6 hfuucjh 2024-11-20 9/450 2024-11-23 14:10 by 4399喜助
[教师之家] 55岁女副教授，博士，该退休还是不退？ +14 aizheteng 2024-11-22 14/700 2024-11-23 11:35 by Quakerbird
[有机交流] 氨基酸乙酯盐酸盐结晶 +3 Wen-1901243 2024-11-22 3/150 2024-11-23 08:08 by 591950582
[考博] 读博求助 10+3 姜大胖仙女 2024-11-18 5/250 2024-11-23 03:53 by dosupy
[教师之家] 坐标西安，找个搭子 +10 hamihappy 2024-11-19 12/600 2024-11-22 23:14 by yiran909
[考博] 招申请考核博士研究生 +4 zhangqi1012 2024-11-22 4/200 2024-11-22 22:24 by 求上岸的
[考博] 24级双9硕钙钛矿材料方向求博导收留 +6 yfdsefx 2024-11-19 8/400 2024-11-22 14:54 by 鱼儿梦水深
[硕博家园] 预答辩被狠批 +10 呵呵orc 2024-11-21 10/500 2024-11-22 14:43 by 120890215
[有机交流] 实验室减压蒸馏二甲苯 20+4 lijiajun09 2024-11-19 8/400 2024-11-21 22:50 by s66230911
[论文投稿] assigning reviewer一周就awaiting ECI decision是什么意思 +3 mollyzhang_2003 2024-11-21 5/250 2024-11-21 20:08 by mollyzhang_2003
[论文投稿] under review 50+7 星辰2022 2024-11-19 13/650 2024-11-21 14:39 by TopEdit
[论文投稿] 返修意见：有关于XRD精修R值过高的问题 5+4 zychiang 2024-11-18 10/500 2024-11-20 19:33 by zychiang
[论文投稿] 请问JCR分区与WOS分区是一样的吗? 5+4 zhongyuan764 2024-11-19 7/350 2024-11-20 14:19 by cqupenghao
[基金申请] 海优出来了吗 +4 fjshida 2024-11-17 5/250 2024-11-20 08:12 by kudofaye
信息提示
关闭
请填处理意见
关闭