| 1 | 1/1 | 返回列表 |
| 查看: 1556 | 回復(fù): 0 | ||
hanchun5085金蟲 (小有名氣)
|
[求助]
一個(gè)循環(huán)判斷問題
|
|
%剪力-變形骨架曲線繪制 clc clear format short %截面尺寸 Nc=140000; %柱頂軸向力N Nb=0; %梁端軸向力N hc=250; %柱橫截面高度mm hb=300; %梁橫截面高度mm hc1=230; %柱橫截面有效高度mm hb1=280; %梁橫截面有效高度mm angle_c=atan(hb1/hc1); %混凝土壓桿角度 angle_s=atan(hb/hc); b=250; %節(jié)點(diǎn)區(qū)寬度 %混凝土壓縮材料屬性 fcu=62; %N/mm^2 試驗(yàn)值或者規(guī)范取值 fc=0.776*fcu; %根據(jù)規(guī)范得到抗壓值,系數(shù)為0.76--0.82 fc2_max=fc; %FRC最大抗壓強(qiáng)度 epsilon_0=-0.0065; %混凝土峰值壓應(yīng)變 A=1.101; %FRC上升段受壓參數(shù) A1=(0.0234+0.1577*0.02)*(0.0016*fc*fc+0.2315*fc+8.4633); %FRC下降段受壓參數(shù) %混凝土拉伸材料屬性 E0=11702; %根據(jù)李艷公式計(jì)算 sigma_tu=4.47; %極限拉應(yīng)力 sigma_ss=0.896*sigma_tu; %最大橋接應(yīng)力 epsilon_10=sigma_ss/E0; Eie=1.618*sigma_tu^2.104; %硬化段斜率 epsilon_tu=(sigma_tu-sigma_ss)/Eie; %受拉極限拉應(yīng)變 sigma_i=sigma_ss*(1-Eie/E0); %鋼筋材料屬性 %梁方向縱筋 fyh=427; %梁縱向屈服強(qiáng)度MPa ph=0.01608; %梁縱向鋼筋面積配箍率 Esv=200000; %梁縱向鋼筋彈性模量 ETv=2000; %梁縱向鋼筋塑性模量 epsilon_yv=0.002135; %梁縱向鋼筋屈服應(yīng)變 %柱方向縱筋 fyv=427; %y向縱筋屈服強(qiáng)度MPa pv=0.013; %節(jié)點(diǎn)豎向鋼筋面積配箍率 Esv=200000; %豎向鋼筋彈性模量 ETv=2000; %豎向鋼筋塑性模量 epsilon_yv=0.002135; %豎向鋼筋屈服應(yīng)變 %節(jié)點(diǎn)區(qū)箍筋 pjv=0.0015; %節(jié)點(diǎn)區(qū)箍筋面積配筋率 fshj=273; %箍筋屈服強(qiáng)度MPa Esh=210000; %水平鋼筋彈性模量 ETh=2100; %水平鋼筋塑性模量 epsilon_yh=0.0013; %水平鋼筋屈服應(yīng)變 co=0.33; %縱筋參與系數(shù) % 第一步:計(jì)算水平鋼筋應(yīng)力d=ph*fsh d=pjv*fshj+co*ph*fyh; %水平鋼筋應(yīng)力ph*fsh % 第二步:水平FRC應(yīng)力fch fch=-d; %水平FRC應(yīng)力fch c=1.0;v=0.9; if abs(c-v)>0.0001 c=c-0.0001; % 第三步:受壓本構(gòu)關(guān)系求水平應(yīng)變epsilon_h fc2_max=c*fc; a=fch/fc2_max; epsilon_h=((-A*a+2*a+A)-sqrt((A*a-2*a-A)*(A*a-2*a-A)-4*a))*epsilon_0/2; % 第四步:剪應(yīng)變Yhv b=epsilon_h*hc1/(hc*(sin(angle_s))*(cos(angle_s))); Yhv=0.5*(b+sqrt(b*b+4*b/125)); cj=hc/(2+250*Yhv); % 第五步:軟化系數(shù)v epsilon_1=(sin(angle_s)*sin(angle_s)*(hc-2*cj)/hb1+0.5*tan(angle_c))*Yhv; v=5.8/sqrt(0.86*fc*(1+400*epsilon_1)); %if abs(c-v)>0.01 % v=v; %end %c=v; end % 第六步:FRC主壓應(yīng)變epsilon_2 epsilon_2=epsilon_h-Yhv*(tan(angle_c))/2; % 第七步:FRC抗壓強(qiáng)度sigma_fc2 c=epsilon_2/epsilon_0; if epsilon_2<epsilon_0 sigma_fc2=-(A*c-c*c)*fc2_max/(1+(A-2)*c); else sigma_fc2=-A1*c*fc2_max/(1+(A1-2)*c+c*c); end % 第八步:FRC主拉應(yīng)變epsilon_1 epsilon_1=epsilon_2+Yhv/2/sin(angle_c)/cos(angle_c); fc2_max=fc/(0.8-0.34*epsilon_1/epsilon_0); % 第九步:主拉應(yīng)力 if epsilon_1>epsilon_tu fc1=0; elseif epsilon_1<=epsilon_tu if epsilon_1<epsilon_10 fc1=E0*epsilon_1; else fc1=sigma_i+Eie*epsilon_1; end end % 第十步:剪應(yīng)力 Thv=(fc1-sigma_fc2)*sin(angle_c)*cos(angle_c); '成功' |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +4 | 零八# 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 4/200 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +6 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +8 | 北道巷 2026-03-24 | 8/400 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|