版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

xuanqi521

新蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 719.7
散金: 211
紅花: 2
帖子: 1058
在線: 111.8小時
蟲號: 3529424
注冊: 2014-11-10
性別: GG
專業(yè): 免疫生物學(xué)

[求助] 求助matlab循環(huán)迭代問題

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));          %樁側(cè)摩阻力
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('不用循環(huán)')
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設(shè)置不合理')
      % 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時，為什么輸出的全是一個值。如果在把精度調(diào)小點(diǎn)如1.360354400372624e-04，為什么只循環(huán)2次。

回復(fù)此樓

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有114人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有62人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)

1樓 2015-12-04 17:07:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

GOzyb

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 794.4
紅花: 2
帖子: 70
在線: 47.5小時
蟲號: 3286418
注冊: 2014-06-21
專業(yè): 巖土與基礎(chǔ)工程

備注太少了點(diǎn)

發(fā)自小木蟲Android客戶端

回復(fù)此樓

2樓2015-12-06 08:03:03