| 查看: 1612 | 回复: 0 | ||
hanchun5085金虫 (小有名气)
|
[求助]
一个循环判断问题
|
|
%剪力-变形骨架曲线绘制 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; %节点区宽度 %混凝土压缩材料属性 fcu=62; %N/mm^2 试验值或者规范取值 fc=0.776*fcu; %根据规范得到抗压值，系数为0.76--0.82 fc2_max=fc; %FRC最大抗压强度 epsilon_0=-0.0065; %混凝土峰值压应变 A=1.101; %FRC上升段受压参数 A1=(0.0234+0.1577*0.02)*(0.0016*fc*fc+0.2315*fc+8.4633); %FRC下降段受压参数 %混凝土拉伸材料属性 E0=11702; %根据李艳公式计算 sigma_tu=4.47; %极限拉应力 sigma_ss=0.896*sigma_tu; %最大桥接应力 epsilon_10=sigma_ss/E0; Eie=1.618*sigma_tu^2.104; %硬化段斜率 epsilon_tu=(sigma_tu-sigma_ss)/Eie; %受拉极限拉应变 sigma_i=sigma_ss*(1-Eie/E0); %钢筋材料属性 %梁方向纵筋 fyh=427; %梁纵向屈服强度MPa ph=0.01608; %梁纵向钢筋面积配箍率 Esv=200000; %梁纵向钢筋弹性模量 ETv=2000; %梁纵向钢筋塑性模量 epsilon_yv=0.002135; %梁纵向钢筋屈服应变 %柱方向纵筋 fyv=427; %y向纵筋屈服强度MPa pv=0.013; %节点竖向钢筋面积配箍率 Esv=200000; %竖向钢筋弹性模量 ETv=2000; %竖向钢筋塑性模量 epsilon_yv=0.002135; %竖向钢筋屈服应变 %节点区箍筋 pjv=0.0015; %节点区箍筋面积配筋率 fshj=273; %箍筋屈服强度MPa Esh=210000; %水平钢筋弹性模量 ETh=2100; %水平钢筋塑性模量 epsilon_yh=0.0013; %水平钢筋屈服应变 co=0.33; %纵筋参与系数 % 第一步：计算水平钢筋应力d=ph*fsh d=pjv*fshj+co*ph*fyh; %水平钢筋应力ph*fsh % 第二步：水平FRC应力fch fch=-d; %水平FRC应力fch c=1.0;v=0.9; if abs(c-v)>0.0001 c=c-0.0001; % 第三步：受压本构关系求水平应变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; % 第四步：剪应变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); % 第五步：软化系数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主压应变epsilon_2 epsilon_2=epsilon_h-Yhv*(tan(angle_c))/2; % 第七步：FRC抗压强度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主拉应变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); % 第九步：主拉应力 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 % 第十步：剪应力 Thv=(fc1-sigma_fc2)*sin(angle_c)*cos(angle_c); '成功' |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
什么时候能放榜呀？
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
循环伏安过程中出现了一个还原峰而出现了两个氧化峰，为什么出现这种现象？
已经有8人回复
- 请大家帮忙看一个循环伏安曲线~
已经有9人回复
- 请教一个C++读文件的问题，
已经有12人回复
- 编一个fortran循环程序求思路
已经有6人回复
- Fortran 函数中的循环无法进行
已经有4人回复
- 用CHI660D如何测循环阳极极化曲线
已经有12人回复
- 【小白必看】JAVA快速入门—《Java学习路线图：Java必须知道的300个问题》
已经有273人回复
- 怎样由循环伏安图判断可逆反应？
已经有5人回复
- 循环伏安曲线只有还原峰没有氧化峰
已经有19人回复
- 循环伏安法测钴配合物的问题
已经有8人回复
- Python中for循环的一个问题
已经有4人回复
- FORTRAN新手 求助主程序循环问题
已经有10人回复
- 请教循环伏安曲线一个问题
已经有10人回复
- 如何判断一个材料是否具有热塑性？
已经有9人回复
- 请教一个循环伏安曲线
已经有12人回复
- 【求助】循环伏安峰面积的一个问题
已经有4人回复
- 【求助】循环伏安测试的主要目的
已经有28人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50