| 查看: 1233 | 回复: 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人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
物理学I论文润色/翻译怎么收费?
已经有291人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有31人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有65人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
2楼2015-12-06 08:03:03
