資源描述：
《非線性分析作業(yè)(偏壓構(gòu)件非線性分析)》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、非線性分析作業(yè)（偏壓構(gòu)件非線性分析）作業(yè)2：對一壓彎鋼筋混凝土或預(yù)應(yīng)力混凝土構(gòu)件的截面進行全過程分析答：鋼筋混凝土壓彎構(gòu)件的全過程分析1.計算假定1）平截面假定梁正截面變形后仍保持平面，截面應(yīng)變?yōu)橹本€分布，不考慮鋼筋與混凝土之間的相對位移。從理論上來講平截面假定僅適用于跨高比較大的連續(xù)均質(zhì)彈性材料的構(gòu)件。對由混凝土及鋼筋組成的構(gòu)件，由于材料的非均質(zhì)性，以及混凝土開裂，特別是在縱筋屈服，受壓區(qū)高度減小而臨近破壞的階段，在開裂截面上的平截面假定已不能適用。但是，考慮到構(gòu)件破壞是產(chǎn)生在某一區(qū)段長度內(nèi)的，而且試驗結(jié)果表明，只要應(yīng)變量測標距有一定長
2、度，量測的截煩平均應(yīng)變值從加荷開始直到構(gòu)件破壞，都能較好地符合平截面假定。2）鋼筋的應(yīng)力■應(yīng)變關(guān)系應(yīng)采用單向加載下，反復(fù)加載下的不同情況而定。3）混凝土應(yīng)力■應(yīng)變關(guān)系選用單軸作用下的混凝土本構(gòu)關(guān)系模型，且在混凝土的受壓區(qū)帶有下降段。4）不考慮剪切變形的影響。2.編程計算舉例1）計算構(gòu)件及計算簡圖鋼筋混凝土梁截面配筋及計算簡圖如下所示，混凝土等級為C30,抗壓強度為fck?20.1MPa,抗拉強度為ftk?2.01MPao鋼筋等級為HRB335,屈服強度fyk?335MPa,軸力N=150kNo2）計算程序程序釆用C語言編寫。釆用分級加應(yīng)變
3、，先迭代求受壓區(qū)高度，再求對應(yīng)的外彎矩、曲率，得到彎矩、曲率數(shù)據(jù)通過excel處理得到彎矩一曲率關(guān)系圖。具休程序如下：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為受壓區(qū)條帶數(shù)，小1、n2為鋼筋跟數(shù)。）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均為二分法迭代區(qū)間值。)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;(初始受壓區(qū)高度迭代值。)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))(由于受壓區(qū)在程序屮的離散性，為防止循環(huán)無法終止，控制循環(huán)次數(shù)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;(用二分法迭代確定初始狀態(tài)截面應(yīng)變。)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;（迭代受壓區(qū)高e=20000;while(((c*c)>50000)&
8、amp;&(j<5000))度循環(huán)。）｛nO=x*n/h;if(nO>n)nO=n;elsenO=nO;tl=O;for(i=l;i<nO;i++){sl=(x-i*
