西南交大混凝土桥作业条带法.doc

《西南交大混凝土桥作业条带法.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、基于MATLAB的纯弯构件M-φ关系分析学院：土木工程学院班级：桥梁二班姓名：日期：2013年11月3日摘要：对于一个确定的钢筋混凝土构件截面，在一个确定的受力条件下，可求出一条M-φ曲线。以下一个具体的算例通过基于matlab的条带法分析程序来分析纯弯构件的正截面全过程受力。并且比对不同配筋率下（少筋、适筋和超筋）构件受弯全过程行为以及分析增加混凝土强度和增加钢筋强度对抗弯行为的影响。1.截面概况钢筋混凝土T梁截面：,,,,,；采用C25混凝土，，，，；HRB335钢筋，，；。2.由规范简化计算2.1判断T形截面类型，故为第一类T形截面。计算极限弯矩设计值：3.条带法分析3.1四个基本假

2、定(1)截面应变保持平面(2)要考虑混凝土的抗拉强度(3)混凝土的应力应变关系曲线如图1，图2(4)钢筋的应力应变曲线如图3图1混凝土受压区应力应变曲线图2混凝土受拉区应力应变曲线图3钢筋的应力应变曲线3.2条带法分析全过程1.首先取；2.取中性轴高度的上下限值和，计算。通过二分法最终迭代出中性轴高度，其中和分别取截面上下缘高度；3.分别将，和代入求轴力的子函数，求出各个中性轴对应的轴力；(1)由，以及截面数据根据变形协调计算每个条带的应变以及钢筋处的应变；(2)由每个条带的应变根据本构关系求每个条带的应力和钢筋的应力，应力乘以条带的面积得条带轴力；（受拉区混凝土应力应变曲线按斜直线考虑，

3、程序中将混凝土受压受拉应力应变曲线考虑在一个图中，以受压为正受拉为负，当超过混凝土的极限拉应变时不考虑混凝土对截面的抗拉贡献）(3)将所有条带的轴力叠加求和，并加上钢筋的拉力即得截面轴力；4.如果进入下一步。否则，如果，令，，否则，令。返回第二步迭代；5.由迭代收敛的中性轴高度计算曲率，并调用计算弯矩的子函数，计算弯矩M。6.如果钢筋的应变超过极限拉应变，则终止计算。否则，每次增加，，返回第二步继续计算。（若程序不输出“钢筋拉断”，则受压区混凝土被压碎）4.分析结论4.1不同配筋率下构件受弯全过程行为（讨论少筋、适筋和超筋）图4不同配筋率下的受弯过程(1)适筋梁采用钢筋面积6836，（非超

4、筋梁），（非少筋梁）由MATLAB计算的极限承载力为,与规范所给的简化计算公式的结果相差，说明了程序计算的可靠性。(2)超筋梁采用钢筋面积11196由判断属于第二类T形截面故为超筋梁。(3)少筋梁采用钢筋面积452故为少筋梁。4.2分析混凝土强度对抗弯行为的影响图5不同混凝土强度下的受弯过程（1）采用适筋梁，通过改变混凝土的标号来分析混凝土强度对承载能力的影响。（2）结论：当混凝土强度较小时增加混凝土强度可以提高梁的抗弯承载能力，但混凝土强度较高时提高混凝土的强度对承载能力的影响不明显。4.3分析钢筋强度对抗弯行为的影响图6不同钢筋强度下的受弯过程(1)采用适筋梁，通过改变钢筋的级别来分析

5、钢筋强度对承载能力的影响。(2)结论：通过改变钢筋的级别来改变钢筋的强度，可以有效提高梁的抗弯承载能力。4.4研究由M-φ曲线求得的刚度与《公路桥规》计算的刚度的差别1.通过MATLAB计算的弯矩和曲率数据，计算出程序求出的刚度2.根据规范给出的公式计算出抗弯刚度本算例分别通过编程计算和手算两种方式计算结构的刚度编程结果为，手算结果为。3.可以比较，规范计算的结果与条带法计算的结果相差在40%左右，如图4所示。图7规范计算与程序计算的刚度差别附录1条带法程序functiondata=Wuyx(fc,fy,Es,ft,Ec,b,h,hf1,bf,As,as,n)%主函数；输入截面特性%输入截

6、面特性data=[];ford=1:2001psc=(d-1)*0.00000165;%混凝土最大压应变循环y1=0;y2=h;axf1=Wuyx2(fc,fy,Es,ft,Ec,psc,y1,b,h,hf1,bf,As,as,n);%调动求轴力函数axf2=Wuyx2(fc,fy,Es,ft,Ec,psc,y2,b,h,hf1,bf,As,as,n);ifaxf1*axf2>0,break,endmax=1+round((log(y2-y1)-log(0.001))/log(2));%计算最大二分法次数forz=1:max%二分法循环求中性轴y0=(y1+y2)/2;axf0=Wuyx2

7、(fc,fy,Es,ft,Ec,psc,y0,b,h,hf1,bf,As,as,n);ifaxf0==0y1=y0;y2=y0;elseifaxf2*axf0>0y2=y0;axf2=axf0;elsey1=y0;axf1=axf0;endify2-y1<0.0001,break,end%控制最小误差终止循环endmaxss=psc/(h-y0)*(y0-as);ifmaxss>0.01,disp('钢筋拉断'),bre