衬砌结构配筋程序

.-衬砌结构配筋程序程序解释本程序根据钢筋混凝土结构中的设计原理来编制，首先读入材料各特征系数，然后建立循环，读入各截面的厚度、轴力和弯矩，如果弯矩为负值，则要转为正值，但算出受拉区和受压区钢筋量后，要将它们互换，最终才得到衬砌内侧和外侧的钢筋量Ass（j）感谢阅读Ass1（j）。算出各截面所需的配筋量，取衬砌内侧的钢筋最大值作为整体衬砌结构内侧每截面的配筋量即maxAs，取衬砌外侧的钢筋最大值作为整体衬砌结构外侧每截面的配筋量maxAs1，式中各主要符号的含义如下：Ak——安全系数；谢谢阅读Rw——混凝土弯曲抗压极限强度标准值；Ra——混凝土弯曲抗压极限强度；N、NN（j）——轴力、轴力数组；M、MM（j）——弯矩、弯矩数组；b——截面宽度（沿隧道走向取单位长度1m）；x——混凝土受压区高度；Rg——钢筋的抗压计算强度标准值；As——受拉区钢筋的截面面积；As1——受压区钢筋的截面面积；精品文档放心下载y1——形心轴到受拉区边缘的距离；e——轴力作用点到到受拉钢筋重心的距离；e1——轴力作用点到受压钢筋重心的距离；精品文档放心下载h0——受压区边缘到受拉钢筋重心的距离；a——受拉钢筋重心到受拉区混凝土边缘的距离；a1——受压钢筋重心到受压区混凝土边缘的距离；e0——偏心距；h——截面高度（即衬砌厚度）；Ec——混凝土的受压弹性模量；Es——钢筋的弹性模量；Wmax——最大裂缝宽度；W——裂缝宽度允许值；afai——构件受力特征系数；csa——裂缝间纵向受拉钢筋应变不均匀系数；fctk——混凝土抗拉强度标准值；rte——纵向受力钢筋配筋率；Ace——有效受拉混凝土截面面积；R——纵向受拉钢筋表面特征系数；c——钢筋保护层厚度；感谢阅读.-sigma——纵向受拉钢筋的应力；d——受拉区钢筋直径；d1——受压区钢筋直径；fai——纵向弯曲系数；px——偏心情况：1，为大偏心；px=2，为小偏心。建立输入文件（1.dat）格式如下：xinghao，biaohao精品文档放心下载Ak，Ra，Rw，Rg，Ec，Es，fai，fc，fy1，fctk，W，r，afaic，b，d，d1谢谢阅读hh(0)，NN(0)，MM(0)hh(1)，NN(1)，MM(1)…hh(i)，NN(i)，MM(i)注：i为隧道衬砌截面个数。dimensionhh(60),Ass(60),Ass1(60),TT(60)realMM(60),NN(60),M,N,maxAs,maxAs1,As,As1integerpx感谢阅读characterxinghao*10,biaohao*10精品文档放心下载open(1,file="1.txt",status='unknown')谢谢阅读open(2,file="2.txt",status='unknown')感谢阅读read(1,*)xinghao,biaohaowrite(2,*)'钢筋型号:'write(2,*)xinghaowrite(2,*)'混凝土标号:'write(2,*)biaohaoread(1,*)AK,Ra,Rw,Rg,Ec,Es,fai,fc,fy1,fctk,W,r,afai感谢阅读read(1,*)c,b,d,d1write(2,*)'配筋参数:'write(2,100)1000*c,1000*b,1000*d,1000*d1谢谢阅读100format(1x,'c=',f8.3,4x,'b=',f8.3,4x,'d=',f8.3,4x,'d1=',f8.3)write(2,*)'******************************************************谢谢阅读*1********'pause10read(1,*)imaxAs=0.-maxAs1=0do10,j=0,iread(1,*)hh(j),NN(j),MM(j)精品文档放心下载if(MM(j).lt.0)thenTT(j)=MM(j)MM(j)=-MM(j)endifh=hh(j)M=MM(j)N=NN(j)e0=M/Na=c+d/2a1=c+d1/2h0=h-ae=e0+(h/2-c-d/2)e1=e0-(h/2-c-d1/2)As1=-0.00000420 As1=As1+0.000004As=-0.000004As=As+0.000004if(As.gt.0.04)goto20精品文档放心下载!/*计算受压区混凝土高度x*/if(e0.le.(h/2-a1))thenx=sqrt((e-h0)**2+2*Rg*(As*e+As1*e1)/(Rw*b))-(e-h0)谢谢阅读elsex=sqrt(ABS((e-h0)**2+2*Rg*(As*e-As1*e1)/(Rw*b)))-(e-h0)精品文档放心下载endif!/*判断大小偏心情况*/if(x.le.(0.55*h0))thenpx=1else.-px=2endif! /*检算截面强度*/if(px.eq.1)thenif(x.ge.(2*a1))thenif((AK*N).gt.(Rw*b*x+Rg*(As1-As)))goto30else感谢阅读if((AK*N*e1).gt.(Rg*As*(h0-a1)))goto30精品文档放心下载endifelseif((AK*N*e).gt.(0.5*Ra*b*h0**2+Rg*As1*(h0-a1)))goto20精品文档放心下载if(e0.le.(h/2-a1))thenif((AK*N*e1).gt.(0.5*Ra*b*h0**2+Rg*As*(h0-a)))goto30endif感谢阅读endif! /*检算按轴心受压构件考虑时的截面强度*/AA=b*hif(N.gt.(0.9*fai*(fc*AA+fy1*(As+As1))))goto30精品文档放心下载! /*检算受拉边缘混凝土裂缝宽度*/Ace=0.5*b*hrta=As/Aceif(rta.lt.0.01)rta=0.01z=(0.87-0.12*(h0/e)**2)*h0感谢阅读if(z.gt.(0.87*h0))z=0.87*h0精品文档放心下载if(As.eq.0)thencsa=0.4elsesigma=N*(e-z)/(As*z)csa=1.1-0.65*fctk*(rta*sigma)感谢阅读endifif(csa.gt.1.0)csa=1.0if(csa.lt.0.4)csa=0.4Wmax=afai*csa*r*(2.7*c+0.1*d/rta)*sigma/Es谢谢阅读.-if(Wmax.le.W.or.e0.le.(0.55*h0))then谢谢阅读Ass(j)=As*1000000Ass1(j)=As1*1000000if(TT(j).lt.0)thensss=Ass(j)Ass(j)=Ass1(j)Ass1(j)=sssendifwrite(2,200)j,1000*h,Ass(j),Ass1(j)精品文档放心下载200format(1x,'j=',i2,4x,'h=',f8.3,4x,'Ass(j)=',f8.1,5x,'Ass1(j)=',f8.1)if(Ass(j).le.0.and.Ass1(j).le.0)then谢谢阅读write(2,*)'按构造要求配筋'elsewrite(2,*)'按受力要求配筋,同时得满足构造要求所以：'endif谢谢阅读if(Ass(j).lt.0.002*b*h*1000000)Ass(j)=0.002*b*h*1000000if(Ass1(j).lt.0.002*b*h*1000000)Ass1(j)=0.002*b*h*1000000精品文档放心下载write(2,400)Ass(j),Ass1(j)精品文档放心下载400 format(1x,'Ass(j)=',f8.1,5x,'Ass1(j)=',f8.1)write(2,*)'******************************************************谢谢阅读**1********'if(Ass(j).gt.maxAs)maxAs=Ass(j)感谢阅读if(Ass1(j).gt.maxAs1)maxAs1=Ass1(j)精品文档放心下载elsegoto30endifcontinuewrite(2,*)'******************************************************感谢阅读*1********'write(2,300)maxAs,maxAs1format(1x,'maxAs=',f8.1,5x,'maxAs1=',f8.1)感谢阅读close(1)close(2).-end根据第3章中素混凝土衬砌结构的内力检算结果，如果检算后内力不符合要求，为了保证衬砌结构安全，必须按受力要求进行配筋。通过检算知道，各类围岩素混凝土衬砌结构均有部分截面不符合要求，但除了Ⅴ级浅埋外大部分只有个别点不通过，且安全系数或容许值已大于1，在实际工程中，这种检算不通过可以忽略。因此按构造要求对Ⅴ级浅埋段进行配筋。谢谢阅读（1）配筋原理隧道衬砌结构属于偏心受压矩形构件，根据偏心受压钢筋混凝土结构设计原理进行配筋。采用试算法，假设受拉区和受压区一定钢筋分别为AS和AS′，检算钢筋混凝土结构的强度和裂缝是否符合要求，直到衬砌每个截面都符合要求为至，最后得到各截面的钢筋情况。这部分可编写程序进行计算。感谢阅读1）计算偏心距e0eM/N（5-4-1）0式中：M——弯矩；N——轴力。2）确定截面受压区高度x先假设衬砌截面受拉钢筋和受压钢筋面积分别为A和A，按下列公式计算谢谢阅读s s出受压区高度x，即可以确定截面中性轴位置。（5-4-2）R(AeAe)Rbx(ehx/2)gssw0解出上式x：x(eh)（5-4-3）0gssw0式中：e——轴向力到受拉钢筋重心的距离，ee(h/2a)；0e——轴向力到受压钢筋重心的距离，ee(h/2a)；谢谢阅读0a——钢筋A的重心到截面受拉边缘的距离；sa——钢筋A的重心到截面受压边缘的距离；s——衬砌厚度；——钢筋的抗拉计算强度标准值/MPa；g——混凝土弯曲抗压极限强度标准值/MPa；w——混凝土受压区边缘到受拉钢筋重心的距离；0——衬砌计算截面宽度，取1m。当轴向力N位于作用于钢筋A和A重心之间时，式（5-4-2）和式（5-4-3）谢谢阅读s s.-中取正值；当N作用于两重心以外时，则取负号。3）截面大小偏心类型的确定如果x055h，则矩形截面为大偏心受压构件，否则为小偏心受压构件。谢谢阅读04）强度检算a、大偏心受压构件如果x2a，按下式进行强度检算：KNRbxR(AA)（5-4-4）wgss如果x2a，按下式进行强度检算：RA(h（5-4-5）gs0b、小偏心受压构件强度检算公式如下：KNe0.5Rbh（5-4-6）RA(ha)2a0gs0如果轴向力N位于钢筋A的重心与钢筋A的重心之间，强度检算还应符合ss下式：KNe0.5Rbh2RA(ha)（5-4-7）a0gs0除了以上按偏心受压构件进行检算外，还应按轴心受压构件验算截面的轴心抗压强度：精品文档放心下载(A)]（5-4-8）N0.9[fAfAcyss式中：——安全系数；——混凝土弯曲抗压极限强度/Mpa；a——钢筋混凝土构件的纵向弯曲系数；c——混凝土轴心抗压强度设计值；——受力钢筋的抗压强度设计值；y——构件截面面积。其它符号意义同前。5）裂缝检算当e<0.55h,可以不进行裂缝检算,否则需要按下式检算衬砌截面的最大精品文档放心下载0 0裂缝宽度 是否符合要求：max (2.7c0.1d/) /E （5-4-8）谢谢阅读max te s s式中：——构件受力特征系数，对偏心受压构件取=2.7；谢谢阅读——裂缝间纵向受拉钢筋应变不均匀系数，1.10.65f /()，当精品文档放心下载ctk te s.-0.4时，取0.4，当1.0时，取1.0，其中，为按有效受拉混凝精品文档放心下载te土面积计算的纵向受拉钢筋配筋率， A/A，当 0.01时，取0.01；精品文档放心下载te s ce te——有效受拉混凝土截面面积，对于偏心受压矩形构件，取A0.5bh，感谢阅读ce ce为截面高度，即衬砌厚度；——纵向受拉钢筋表面特征系数，螺旋钢筋取0.7，光面钢筋取1.0；感谢阅读——混凝土保护层厚度；——纵向受拉钢筋的应力，N(ez)/(Az)，其中z为纵向受拉钢筋谢谢阅读s s s合力点至受压区合力之间的距离，则取z0.87h。

z[0.87

0.12(h0

/e)

2]h0

，且z

0.87h0

，否0采用以上的钢筋混凝土结构配筋方法，因为衬砌各截面内力不同，按受力要求各截面所需配筋量不同，为了施工时铺设钢筋方便，以及整个衬砌结构受拉（或受压）边并不是都是位于衬砌内侧或外侧，拱顶内侧受拉，而边墙等其他截面可能是外侧受拉，完全按受力要求配筋会给施工带来很大麻烦，因此为了满足衬砌结构每个截面受力配筋要求，按所需钢筋量最多的截面对整体衬砌的每个截面进行配筋。精品文档放心下载（3）编写配筋程序隧道衬砌结构截面有几十个，如果采用手工计算，显然很麻烦，因此采用编制程序来配置钢筋。根据以上配筋原理，程序假设衬砌结构内侧和外侧钢筋分别为A和A，整个衬砌结构受拉（或受压）边并不是都是位于衬砌内侧或外侧，精品文档放心下载s s拱顶内侧受拉，而边墙等其他截面可能是外侧受拉，也就是说并不是所有截面试算中的A都为受拉钢筋，A都为受压钢筋，因此为了与上面钢筋配置原理一致，精品文档放心下载s s衬砌结构外侧受拉时，要将A

和A倒过来，即受拉钢筋为A，受压钢筋为A精品文档放心下载

，s

s

s

s而其他的参数受拉钢筋直径d与受压钢筋直径d也要倒过去。而最后配筋结果恢复原来假设的钢筋。其实主要是考虑到偏心受压构件的等效轴力N的作用点的问题，衬砌内侧受拉与外侧受拉的等效轴力N的作用点位置分别位于衬砌轴线的内侧和外侧，因此需要作出这样的处理。精品文档放心下载衬砌配筋程序中假设受拉钢筋直径为d，受压钢筋直径为d，保护层厚度为感谢阅读，具体程序及其解释详见附录C。设计原理为，首先读入材