 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1228 | 回復(fù): 1 | ||
[求助]
求助matlab循環(huán)迭代問(wèn)題
|
|
clear clc format long S1=0.01; %樁身第n單元頂位移 deL=[2 2 1.5]; %分段長(zhǎng)度 dln=deL(1) E=3.5*10^4; %樁的彈性模量 A=0.2922; %樁橫截面面積 U=1.9164; %樁周長(zhǎng) 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時(shí)，為什么輸出的全是一個(gè)值。如果在把精度調(diào)小點(diǎn)如1.360354400372624e-04，為什么只循環(huán)2次。 |
回復(fù)此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有155人回復(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ù)
-
請(qǐng)問(wèn)還有沒(méi)有用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ì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
2樓2015-12-06 08:03:03
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xuanqi521 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +8 | 何潤(rùn)采123 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +8 | 薛云鵬 2026-03-13 | 8/400 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|