| 查看: 1341 | 回复: 1 | ||
[求助]
求助matlab循环迭代问题
|
|
clear clc format long S1=0.01; %桩身第n单元顶位移 deL=[2 2 1.5]; %分段长度 dln=deL(1) E=3.5*10^4; %桩的弹性模量 A=0.2922; %桩横截面面积 U=1.9164; %桩周长 Pn0=1.10861617800175; %桩身n-1单元底轴力 for X=1:length(deL) qs=0.9*S1/(1/(0.24095*55.25)+S1/(60*10^-3)); %桩侧摩阻力 dPn=qs*U*dln; %轴力增量 Pn=Pn0+dPn; %n-1单元底轴力 Pm=(Pn+Pn0)/2; %n-1单元中心截面处轴力 dSn=Pm*dln/(E*A); %n-1单元平均变形 Sn=S1+dSn; %n-1单元顶位移 qs1=0.9*Sn/(1/(0.24095*55.25)+Sn/(60*10^-3)); dPn1=U*E*A*qs1*dSn/Pn; Pn1=Pn0+dPn1; Pm1=(Pn0+Pn1)/2; dSn1=Pm1*dln/(E*A); Sn1=S1+dSn1; Wucha1=abs(dPn1-dPn)/dPn1; if Wucha1<0.001 disp('不用循环') end qs2=0.9*Sn1/(1/(0.24095*55.25)+Sn1/(60*10^-3)); dPn2=U*E*A*qs2*dSn1/Pn1; Pn2=Pn0+dPn2; Pm2=(Pn0+Pn2)/2; dSn2=Pm2*dln/(E*A); Sn2=S1+dSn2; Wucha2=abs(dPn2-dPn1)/dPn2; disp('Wucha2') disp(Wucha2) if Wucha2<0.001 disp(n) disp(Pn1) disp(Sn1) end count=1; while Wucha2>0.002 % if count>10000 % disp('初始位移S0设置不合理') % break %end Sn1=Sn2; dSn1=dSn2; Pn1=Pn2; qs3=0.9*Sn1/(1/(0.24095*55.25)+Sn1/(60*10^-3)); dPn3=U*E*A*qs3*dSn1/Pn1;disp(dPn3) Pn3=Pn0+dPn3; Pm3=(Pn0+Pn3)/2; dSn3=Pm3*dln/(E*A); Sn3=S1+dSn3; Wucha2=abs(dPn3-dPn2)/dPn3; count=count+1; disp('count') disp(Wucha2) disp(count) disp(Pn1) disp(Sn1) end end 求各位大神帮忙解决wucha2>0.000001时，为什么输出的全是一个值。如果在把精度调小点如1.360354400372624e-04，为什么只循环2次。 |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有158人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
2楼2015-12-06 08:03:03