| 查看: 1285 | 回复: 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次。 |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有237人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
2楼2015-12-06 08:03:03
50