24小时热门版块排行榜

返回列表

hanchun5085

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 673.2
散金: 5
红花: 3
帖子: 209
在线: 82.5小时
虫号: 1957193
注册: 2012-08-26
专业: 结构工程

[求助] 一个循环判断问题

%剪力-变形骨架曲线绘制
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);
'成功'

回复此楼