资源描述：
《非线性分析作业(偏压构件非线性分析)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、非线性分析作业（偏压构件非线性分析）作业2：对一压弯钢筋混凝土或预应力混凝土构件的截面进行全过程分析答：钢筋混凝土压弯构件的全过程分析1.计算假定1）平截面假定梁正截面变形后仍保持平面，截面应变为直线分布，不考虑钢筋与混凝土之间的相对位移。从理论上来讲平截面假定仅适用于跨高比较大的连续均质弹性材料的构件。对由混凝土及钢筋组成的构件，由于材料的非均质性，以及混凝土开裂，特别是在纵筋屈服，受压区高度减小而临近破坏的阶段，在开裂截面上的平截面假定已不能适用。但是，考虑到构件破坏是产生在某一区段长度内的，而且试验结果表明，只要应变量测标距有一定长
2、度，量测的截烦平均应变值从加荷开始直到构件破坏，都能较好地符合平截面假定。2）钢筋的应力■应变关系应采用单向加载下，反复加载下的不同情况而定。3）混凝土应力■应变关系选用单轴作用下的混凝土本构关系模型，且在混凝土的受压区带有下降段。4）不考虑剪切变形的影响。2.编程计算举例1）计算构件及计算简图钢筋混凝土梁截面配筋及计算简图如下所示，混凝土等级为C30,抗压强度为fck?20.1MPa,抗拉强度为ftk?2.01MPao钢筋等级为HRB335,屈服强度fyk?335MPa,轴力N=150kNo2）计算程序程序釆用C语言编写。釆用分级加应变
3、，先迭代求受压区高度，再求对应的外弯矩、曲率，得到弯矩、曲率数据通过excel处理得到弯矩一曲率关系图。具休程序如下：include"stdio.h"#include"conio.h"#include<math.h>main（）{FILE*fsl;FILE*fs2;FILE*fs3;FILE*fs4;intn=10000,n0,nl,n2,ij;（将截面分成10000份，nO为受压区条带数，小1、n2为钢筋跟数。）floath,bdd：Ld241,t2,t3,t4xy$sl,s2,s3,s4,s5Jl，As：LAs2,c,M/
4、Ml/IVI2,r/sy,fyJ=3000?w,N=150000,s0=0.00025,st,cl,a,e,al=0/el=0.0033/c2,t5,t6,t7,s6,s7;(a、e、al>el均为二分法迭代区间值。)printf("inputh?b,fy,dl,d2,nl,n2:");scanf("%f%f%f%f%f%d%d",&h,&b,&fy,&dl,&d2,&nl,&n2);d=h/n;x=10000;(初始受压区高度迭代值。)sy=fy/200000;Asl=pa*n
5、l*dl*dl/4;As2=pa*n2*d2*d2/4;fsl=fopen("r.txt","w");fs2=fopen("M.txt","w”)；fs3=fopen("w.txt","w");fs4=fopen("nO.txt"z"w");cl=10000;j=l；while(((cl*cl)>100000)&&(j<1000))(由于受压区在程序屮的离散性，为防止循环无法终止，控制循环次数j。){cl=20.1*h*b*(l-(l-s0/0.002)*(l-s0/0.002))+200000*(Asl+As
6、2)*s0-N;sO=sO+0.000001;if(cl>=0)el=s0;elseal=s0;s0=al+(el-al)/2;(用二分法迭代确定初始状态截面应变。)j=j+l；}printf(,,s0=%10.8f:,,,s0);j=l；c2=10000;al=0;el=0.0033;st=0.0001;while(((c2*c2)>100000)&&(j<1000)){tl=0;for(i=l;i<n+l;i++){s5=(h-i*d)*st/h;tl=tl+20.1*d*b*(l-(l-s
7、5/0.002)*(l-s5/0.002));s6=(h-as2)*st/h;t3=200000*s6*As2;s7=asl*st/h;t4=200000*s7*Asl;c2=tl+t3+t4-N;if(c2>=O)el=st;elseal=st;st=al+(el-al)/2;j=j+l;}printf(ust=%10.8f:H,st);s=sO+0.000001;while(s<0.0033){c=1000000;j=l；x=5000;a=0;（迭代受压区高e=20000;while(((c*c)>50000)&
8、amp;&(j<5000))度循环。）｛nO=x*n/h;if(nO>n)nO=n;elsenO=nO;tl=O;for(i=l;i<nO;i++){sl=(x-i*
