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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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)小點如1.360354400372624e-04，為什么只循環(huán)2次。

回復(fù)此樓

» 猜你喜歡

重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有200人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有21人回復(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ǔ)工程

備注太少了點

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

回復(fù)此樓

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