| 查看: 1176 | 回复: 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次。 |
回复此楼» 猜你喜欢
🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队
已经有0人回复
-
第一性原理计算方向2026级博士申请 PRB*1，四级484
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有174人回复
-
求助VISSIM破解版软件
已经有0人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有13人回复
2楼2015-12-06 08:03:03