《加藤模型——刚劲混凝土非线性分析Matlab程序.pdf》由会员分享,可在线阅读,更多相关《加藤模型——刚劲混凝土非线性分析Matlab程序.pdf(4页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、钢筋混凝土非线性分析加藤模型solution:(1)钢材弹塑性应力强化模型:steel_sigma_epsilon_double.mfunction sig_s=steel_sigma_epsilon_double(eps_s,fy,Es,Es2)%根据钢筋双线模型计算计算钢筋应力if nargin=3Es2=0.01*Es;endey=fy/Es;if abs(eps_s)es(i-1)%判断是加载还是卸载,s=1表示加载,s=-1表示卸载s=1;else s=-1;endwhile 2%画出该次加载的曲线e=e+de*s;f=steel_sigma_epsilon_double(e,fy,
2、Es);plot(e,f)if abs(e)=abs(es(i)breakendendi=i+1;if abs(es(i-1)ey%如果前一次加载已经发生塑性变形,跳出循环breakendendeq=abs(e);%记录加载结束后应力峰值对应于骨架曲线中的应变%塑性滞回过程while 1f0=f;e0=e;%记录加载结束后的应力应变if es(i)es(i-1)s=1;else s=-1;endwhile 2%卸载过程f=f+df*s;e=e+df/Es*s;if f*s=0%如果应力变为0,则卸载完成,跳出循环breakendplot(e,f)endf1=f;e1=e;x=0;dx=0.00
3、2;%软化过程e_res=abs(e0)-ey;Eb=-Es/6*log(10*e_res);%软化段割线刚度a=Es/(Es-Eb);%软化段曲线系数while 2x=x+dx;y=a*x/(x+a-1);e=e1+abs(f0)/Eb*x*s;f=f1+abs(f0)*y*s;if abs(f)=abs(f0)%若达到前一次加载应力峰值,则软化段结束,跳出循环breakendplot(e,f)endwhile 2%强化阶段e=e+de*s;eq=eq+de;%等效于骨架曲线中的应变f=steel_sigma_epsilon_double(eq,fy,Es)*s;if abs(e)=abs(es(i)breakendplot(e,f)endi=i+1;if i=length(es)+1breakend%pauseendif es(i-1)=0%如果应力变为0,则卸载完成,跳出循环breakendplot(e,f)endxlabel(应变);ylabel(应力(MPa));title(钢材在拉压反复循环加载下的应力应变曲线)grid on-5-4-3-2-1012345x 10-3-500-400-300-200-1000100200300400500 (MPa)