桥梁荷载横向分布的计算讲.ppt

《桥梁荷载横向分布的计算讲.ppt》由会员分享，可在线阅读，更多相关《桥梁荷载横向分布的计算讲.ppt（112页珍藏版）》请在三一办公上搜索。

1、桥梁结构分析的计算机方法,桥梁荷载横向分布的计算,横向分布定义,整体桥梁结构应采用影响面加载计算最不利荷载,横向分布定义,为简化计算，采用近似影响面来加载近似影响面纵横方向分别相似,轴重,轴重与轮重的关系,各纵向影响线比例关系,横向分布定义,横向分布定义,概念的正确性内力分布、荷载分布在那种情况下是正确横向分布系数沿纵向的变化弯矩剪力,刚性横梁法,刚性横梁法,偏心受压法、修正偏心受压法基本假设横梁刚度无穷大,ej:到纵轴的偏心距,刚性横梁法,由桥梁工程课程知，横向影响线（直线）的计算公式如下，其中n为主梁数：为考虑抗扭刚度影响的修正系数：,刚性横梁法,当各主梁截面相同，且间距相等时：式中：a,

2、I,It为每个主梁的间距、抗弯惯矩和抗扭惯矩。,刚性横梁法,横向分布系数计算公式 2004公路桥梁规范的荷载如何考虑？,刚性横梁法,计算程序程序功能：计算第i_beam主梁在itype荷载下的横向分布系数co。,刚性横梁法,变量说明,控制变量Ib计算类型指数，等于0时各梁截面相同，间距相等；等于1时为一般情况。Itype荷载类型指数，等于1、2、3时分别计算汽车荷载、挂车荷载和人群荷载。N_main主梁片数i_beam要计算横向分布系数的主梁号（从左到右）N_lane车道数Wid_cro人行道宽度,变量说明,wid_lan车道宽Al简支梁的跨度Ag材料的剪切模量Ae材料的弹性模量控制数组Xin

3、f(n_main)-影响线的横坐标Aii(n_main)各主梁抗弯惯矩数组Aiti(n_main)各主梁抗扭惯矩数组,程序源代码,ccc Main program ccdimension aai(30),aii(30),aiti(30)dimension xinf(30),yinf(30)character dat*12,out*12,line*80,程序源代码,write(*,(/)write(*,*)write(*,(ls,a)Enter the name of datafile:read(*,(a)dat 文件名write(*,(/)call cname(dat)去后缀,程序源代码,ou

4、t=datcall fname(dat,.dat)加后缀.datcall fname(out,.out)加后缀.outopen(10,file=dat,status=unknown)open(11,file=out,status=unknown)read(10,*)line 标题read(10,*)itype,n_main,I_beam,n_lane,wid_cro,wid_lan,程序源代码,read(10,*)(xinf(i)=1,n_main)影响线横坐标read(10,*)ib 控制变量，=0各梁同，=1不同write(*,10)10format(/lx,=input data=:)w

5、rite(*,(4i7,2f10.4)itype,n_main,i_beam,n_lane,wid_cro,wid_lanwrite(*,(5f10.4)(xinf(i),i=1,n_main)write(*,(i7)ib,程序源代码,if(ib.eq.1)thenread(10,*)al,ag,ae 跨度l,G,Eread(10,*)(aii(i),i=1,n_main)Iiread(10,*)(aiti(i),i=1,n_main)Itido i=1,n_mainaai(i)=abs(xinf(i)aiend do,程序源代码,else if(ib.eq.0)thendo i=1,n_ma

6、inaii(i)=1.0aiti(i)=1.0end dodx=abs(xinf(n_main)-xinf(1)/float(n_main-1)相邻梁距离ado i=1,n_mainaai(i)=abs(xinf(i)/dx ej/aend doend if,程序源代码,call inf(i_beam,ib,n_main,aai,aii,aiti,al,ag,ae,xinf,yinf)求影响线yinfcall co_eff(itype,n_lane,n_main,wid_lan,wid_cro,xinf,yinf,co2,co,ilane)求分布系数cowrite(11,(10 x,a)-,程

7、序源代码,write(11,(10 x,a)|write(11,(10 x,a)|C61GXHL print filewrite(11,(10 x,a)|write(11,(10 x,a)-write(11,*)write(11,100)100format(/10 x,=rigid crossbeam influence line:)write(11,200)200format(10 x,=xinf(1-n_main):)write(11,300)(xinf(i),i=1,n_main),程序源代码,300format(10 x,5f12.6)write(11,400)400format(10

8、 x,=yinf(1-n_main):)write(11,500)(yinf(i),i=1,n_main)500format(10 x,5f12.6)write(11,600)600format(10 x,=crosswise distribution coefficient:)write(11,700)co700format(10 x,f12.6),程序源代码,close(10)close(11)write(*,(/)write(*,a)The program C61GXHL over!endC-,程序源代码,subroutine cname(c)c-character*12 cdo 10

9、i=12,1,-1if(c(i:i).eq.)thenip=igoto 30endif,程序源代码,10continue30if(ip.gt.1)thendo 40 i=ip,12c(i:i)=endifendc-,程序源代码,subroutine fname(c,ce)c-character c*12,ce*4k=0do 10 i=1,8if(c(i:i).eq.)thenk=Igoto 15endif 10 continue,程序源代码,if(k.eq.0)k=9do 20i=k,12 20c(i:i)=do 25 i=k,k+3 j=i-k+125 c(i:i)=ce(j:j)end,程

10、序源代码,c-subroutine lin_ins(xl,y1,x2,y2,x,y)线性内插c-！linear insertion!y=y1*(x-x2)/(x1-x2)+y2*(x-x1)/(x2-x1)y=(y1*(x-x2)-y2*(x-x1)/(x1-x2)end,程序源代码,c-subroutine eta(ib,n_main,i_beam,j,aai,aii,aiti,al,ag,ae,xinf,etaij)求ijc-!ETAijdimension aai(n_main),aii(n_main),aiti(n_main),xinf(n_main),程序源代码,if(xinf(i_b

11、eam)*xinf(j).ge.0.0)thensign=1.0elsesign=-1.0end ifq1=0.0q2=0.0do i=1,n_main q1=q1+aii(i)q2=q2+aai(i)*aai(i)*aii(i)end do,I,J梁是否同侧判断,程序源代码,call pass(0，177,q1,q2,aai(1),aii(1)if(ib.eq.0)then b=1.0else q3=0.0 do i=1,n_main q3=q3+aiti(i)Iti end do b=1.0/(1.0+ag*al*al*q3/(12.0*ae*q2)end ifetaij=aii(i_be

12、am)*q1+sign*aai(j)*aai(i_beam)*aii(i_beam)*b/q2 vjend,是否考虑抗扭刚度的影响,程序源代码,c-subroutine inf(i_beam,ib,n_main,aai,aii,aiti,al,ag,ae,xinf,yinf)求横向影响线c-!influence line by rigid crossbeamdimension aai(n_main),aii(n_main),aiti(n_main)dimension xinf(n_main),yinf(n_main)call eta(ib,n_main,i_beam,1,aai,aii,ait

13、i,al,ag,ae,xinf,yinf1)求ij yinf(1)=yinf1,1号梁,程序源代码,call eta(ib,n_main,i_beam,n_main,aai,aii,aiti,al,ag,ae,xinf,yinfn)yinf(n_main)=yinfndo i=2,n_main-1 call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),xinf(i),yinf(i)线性内插求中间各值 end do end,n号梁,程序源代码,c-subroutine co_eff(itype,n_lane,n_main,wid_lan,w

14、id_cro,xinf,yinf,co2,co,ilane)求横向分布系数c-!lateral(or crosswise)distribution coefficientdimension xinf(n_main),yinf(n_main),程序源代码,if(yinf(n_main).gt.yinf(1)then sign=1.0else sign=-1.0end ifco=0.0if(itype.eq.1)then 汽车 qq=0.0,程序源代码,x=(wid_lan*0.5-0.5)*sign do i=1,n_lane call lin_ins(xinf(1),yinf(1),xinf(

15、n_main),yinf(n_main),x,v)求最靠边轮位处影响线值 qq=qq+v x=x-1.8*sign call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)求另一轮位处影响线值 qq=qq+v if(n_lane.gt.2.and.i.eq.2)then co2=co 2列车的分布系数 end if,程序源代码,if(co.lt.qq)then co=qq*0.5 求co ilane=I end if x=x-1.3*sign 车列间距 end doend ifif(itype.eq.2)then x=(wid_la

16、n*0.5-1.0)*sign 挂车 call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v),程序源代码,co=co+v do i=1,3 x=x-0.9*sign call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)co=co+v end do co=co*0.25end if,程序源代码,if(itype.eq.3)then if(yinf(1).ge.0.0)then 人群 x=-wid_lan*0.5-0.5*wid_cro call lin_ins(xin

17、f(1).yinf(1),xinf(n_main),yinf(n_main),x,v)按人行道中心线位置计算 co=co+v end if if(yinf(n_main).gt.0.0)then x=wid_lan*0.5+0.5*wid_cro call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)co=co+v end if,程序源代码,end if endsubroutine pass(iunit,line,v1,v2,v3,v4)输出character form*20form=(4(lx,f12.6)if(iunit.gt

18、.0)then write iunit,(a,i6)=pass line=,line,程序源代码,write(iunit,form)v1,v2,v3,v4else write(*,(a,i6)=pass line=,line write(*,form)v1,v2,v3,v4end ifend,问题思考,对于整体箱梁（型梁）如果计算其横向分布系数大跨度桥梁中横向分布的考虑偏载系数多梁式T型梁桥按单根梁计算和按整体计算各自的优缺点,刚（铰）接梁（板）法,将多梁式桥梁简化为数根并列而相互间横向刚（铰）接的狭长板（梁）各主梁接缝间传递剪力、弯矩、水平压力、水平剪力用半波正弦荷载作用在某一板上，计算各板

19、（梁）间的力分配关系,刚（铰）接梁（板）法,刚（铰）接梁（板）法,以各结点处未知力Mi，Qi为未知数，列力法方程求解 设:则i 结点的力法方程为（协调方程）：,对于i=2,3,m-2来说,刚（铰）接梁（板）法,i=1时，无 项；i=m-1时，无 项。以上各式中下标意义：第一个为所考虑的未知数序号第二个为产生影响因素的未知数序号对每个结点列方程，得力法正则方程：,刚（铰）接梁（板）法,柔度矩阵A 的计算注意下标，偶数号为剪力，奇数号为弯矩。,刚（铰）接梁（板）法,结点处的基本位移模式注意图中位移符号,刚（铰）接梁（板）法,刚（铰）接梁（板）法,载变位矩阵B的计算,P=1作用在k梁中线 P=1作用

20、在k梁右端,P=1作用在1#梁左端,刚（铰）接梁（板）法,横向影响线 令P=1分别作用于1#梁左端、中间各结点、各梁中线、m#梁右端，共2m+1个位置，则可求得横向影响线。,刚（铰）接梁（板）法,主梁扭矩影响线,i=1,i=2,3,m-1,i=m,计算机程序,图 3-11,图 3-12,主程序框图 调用关系框图,变量说明,控制变量N_main主梁根数Nod节点数，nod=n_main-1Ndinf影响线点数，ndinf=2*n_main+1Inform1为0时各主梁截面相同，为1时各主梁截面不同Inform2为0时各主梁间同时有铰结和刚接，为1时各主梁间全部铰结，为2时各主梁间全部刚结Span

21、主梁跨度Gg主梁材料剪切模量Ee-主梁材料弹性模量,变量说明,控制数组Inform(nod)主梁联结方式数组Bi(n_main)各主梁抗弯惯矩数组Ti(n_main)各主梁抗扭惯矩数组Bmain(n_main)主梁宽度（两个节点间的距离）数组Rmain(n_main)主梁悬臂板长度数组Bic(n_main)主梁悬臂板（沿梁长方向单位长度）抗弯惯矩数组Xinf(ndinf)影响线横向坐标数组Eta(n_main,ndinf)存放主梁横向分布影响线，数组每一行为一片主梁的横向分布影响线Tor(n_main,ndinf)存放主梁扭矩影响线Tm(nod,ndif)存放节点处横向玩具影响线,刚(铰)接梁

22、(板)实例,刚(铰)接梁(板)实例,程序解读,程序源代码：program c62gxhl!-transverse distribution influence linedimension inform(39),ik(78)dimension bmain(40),rmain(40),bi(40),ti(40),bic(40),xinf(81)dimension gk(78,4),b(78),eta(40,81),tm(39,81)character dat*12,out*12,line*80write(*,(/)write(*,*)write(*,(lx,a)Enter the name of

23、datafile:read(*,(a)dat,程序解读,write(*,(/)call cname(dat)out=datcall fname(dat,.dat)call fname(out,.out)open(10,file=dat,status=unknown)read(10,*)line open(11,file=dat,status=unknown)read(10,*)n_main,inform1,inform2read(10,*)span,gg,eenod=n_main-1 结点数ndinf=2*n_main+1=1作用点数,程序解读,m=n_mainif(inform1.eq.0)

24、m=1 各主梁同do i=1,m read(10,*)bi(i),ti(i),bmain(i),rmain(i),bic(i)end doif(inform1.eq.0)then do i=2,n_main bi(i)=bi(i-1)ti(i)=ti(i-1)bmain(i)=bmain(i-1)bic(i)=bic(i-1)rmain(i)=rmain(i-1),程序解读,end do end if if(inform2.eq.0)thenread(10,*)(inform(i),i=1,nod)刚、铰接混合 end if if(inform2.eq.1)thendo i=1,nod inf

25、orm(i)=1 铰接end do end if if(inform2.eq.2)thendo i=1,nod,程序解读,inform(i)=2end do end if call inf(n_main,nod,ndinf,span,ee,gg,bi,ti,bmain,inform),+ik,rmain,bic,xinf,gk,b,eta,tor,tm)call output(11,n_main,nod,ndinf,xinf,ta,tor,tm)close(10)close(11)write(*,(/)write(*,(a)Thats all for the program C62GJJBL!

26、end,程序解读,ccc Subroutines for C62GJJBLcccccc-bubroutine cnamec-,程序解读,character*12cdo 10i=12,1,-1 if(c(i:i).eq.)then ip=i goto 30 endif10 continue30 if(ip.gt.1)thendo 40i=ip,1240 c(i:i)=endifend,程序解读,c-subroutine fname(c,ce)c-character c*12,ce*4k=0do 10 i=1,8 if(c(i:i).eq.)then k=i goto 15endif10 cont

27、inue15 if(k.eq.0)k=9,程序解读,do 20 i=k,1220 c(i:i)=do 25 i=k,k+325 c(i:i)=ce(j:j)endc-subroutine decom(neq0,neq,nband,gk)c-!Decomposing K=LDLby CHOLESCKY inpoved square root method.Dimension gk(neq0,4),t(4)do i=1,neq,程序解读,1=i-nband+1 if(l.lt.1)l=1 do j=1,i ji=nband-i+j qq=0.0 if(j.ne.1)then j1=j-1 do k

28、=1,j1ki=nband+k-ikj=nband+k-jqq=qq+t(ki)*gk(j,kj)end do,程序解读,end if if(i.eq.j)then gk(i,nband)=gk(i,nband)-qq else if(i.gt.j)then t(ji)=gk(i,ji)-qq gk(i,ji)=t(ji)/gk(j,nband)end ifend doend do end do endc-,程序解读,subroutine back(neq0,neq,nband,gk,b)c-!Return substitution and forward substitutiondimens

29、ion gk(neq0,4),b(neq)do i=2,neq qq=0.0 l=i-nband+1 if(l.lt.1)1=1 il=i-1 do j=1,il ji=nband-i+j qq=qq+gk(i,ji)*b(j),程序解读,end do b(i)=b(i)-qq end do do i=1,neq b(i)=b(i)/gk(i,nband)end do do ii=1,neq-1 i=neq-ii qq=0.0 kn=neq if(neq.gt.i+nband-1)kn=I+nband-1 il=i+1,程序解读,do j=il,kn ij=nband+i-j qq=qq+gk

30、(j,ij)*b(j)end do b(i)=b(i)-qq end do endc-subroutine inf(n_main,nod,ndinf,span,ee,gg,bi,ti,bmain,inform,+ik,rmain,bic,xinf,gk,b,eta,tor,tm)c-dimension inform(nod),ik(2*nod),程序解读,dimension bi(n_main),ti(n_main),bmain(n_main),rmain(n_main)dimension bic(n_main),eta(n_main,ndinf),tor(n_main,ndinf)dimen

31、sion tm(nod,ndinf),xinf(ndinf),gk(2*nod,4),b(2*nod)dimension c(2,4),d(5)nband=4neq0=2*noddo ii=1,ned do k=1,2 i=(ii-1)*2+k ik(i)=1 if(k.eq.1.and.inform(ii).eq.1)ik(i)=0 end do,程序解读,end dok=0do i=1,neq0 k=k+ik(i)if(ik(i).ne.0)ik(i)=kend doneq=kpi1=3.1415926pi2=pi1*pi1pi4=pi2*pi2s2=span*spando ii=1,ne

32、d,程序解读,do i=1,2 do j=1,4 c(i,j)=0.0 end do end do fi=s1/gg/ti(ii)/pi2 fi1=s2/gg/ti(ii+1)/pi2 wi=s2*s2/ee/bi(ii)/pi4 wi1=s2*s2/ee/bi(ii+1)/pi4 cti=rmain(ii)/ee/bic(ii)ui=cti*rmain(ii)*0.5,程序解读,f=ui*rmain(ii)/1.5 cti1=rmain(ii+1)/ee/bic(ii+1)uil=ctil*rmain(ii+1)*0.5 f1=ui1*rmain(ii+1)/1.5 c(1,3)=fi+fi

33、l+ci1+cti1c(2,3)=fi*bmain(ii)*0.5+ui-fil*bmain(ii+1)*0.5-uilc(2,4)=wi+wi1+f+f1+fi*bmain(ii)*bmain(ii)*0.25+fi1*bmain(ii+1)*bmain(ii+1)*0.25 if(ii.gt.1)then c(1,1)=-fi c(1,2)=fi*bmain(ii)*0.5 c(2,2)=-wi+fi*bmain(ii)*bmain(ii)*0.25,程序解读,if(inform(ii).ne.1)then c(2,1)=-c(1,2)i=ik(2*ii-1)j=nband-i+ik(2*

34、ii-3)if(j.ge.1)then gk(i+1,j-1)=c(2,1)end if gk(i,j)=c(1,1)j=nband-i+ik(2*ii-2)gk(i,j)=c(1,2)end if i=ik(2*ii),程序解读,j=nband-i+ik(2*ii-1)gk(i,j)=c(2,2)end if if(inform(ii).ne.1)then i=ik(2*ii-1)j=nband gk(i,j)=c(1,3)gk(i+1,j-1)=c(2,3)end if i=ik(2*ii)j=nband gk(i,j)=c(2,2),程序解读,end do call decom(neq0

35、,neq,nband,gk)qq=0.0 do i=1,n_main qq=qq+bmain(i)end do xinf(1)=-qq*0.5do i=1,n_main xinf(2*i)=xinf(2*i-1)+bmain(i)*0.5 xinf(2*i+1)=xinf(2*i)+bmain(i)0.5end do!-circulation in load-added point,程序解读,do k=1,ndinf kb=k/2 if(k.eq.1)kb=1 do i=1,neq b(i)=0.0 end do fi=s2/gg/ti(kb)/pi2 w=s2*s2/ee/bi(kb)/pi

36、4 ct=rmain(kb)/ee/bic(kb)u=ct*rmain(kb)*0.5 f=u*rmain(kb)/1.5 d(1)=w,程序解读,d(2)=-fi*bmain(kb)*0.5 d(3)-w+fi*bmain(kb)*20.25 d(4)=fi*bmain(kb)*0.5+u d(5)w+fi*bmain(kb)*bmain(kb)*0.25+f if(k.eq.1)then i=ik(2)b(i)=-d(3)i=ik(1)if(i.gt.0)b(i)=d(2)else if(kb*2.ne.k)then if(k.ne.3)then,程序解读,i=ik(2*kb-3)if(i

37、.gt.0)b(i)=d(2)i=ik(2*kb-2)b(i)=d(3)end if i=ik(2*bk)b(i)=d(5)i=ik(2*kb-1)if(i.gt.0)b(i)=d(4)else i=ik(2*kb)b(i)=d(1),程序解读,if(kb.ne.1)then i=ik(2*kb-2)b(i)=-d(1)end if end if end if call back(neq0,neq,nband,gk,b)!transverse distribution influence line;!torsion moment influence line for each beam.Do

38、i=1,n_main if(i.eq.1)then,程序解读,i2=ik(2*i)b2=b(i2)eta(i,k)=1.0-b2if(kb.ne.1)eta(i,k)=-b2i1=ik(1)b1=0.0if(i1.gt.0)b1=b(i1)tor(i,k)=-b1-b2*bmain(i)*0.5if(k.eq.1)tor(i,k)=-b1-(b2+1.0)*bmain(i)*0.5if(k.eq.3)tor(i,k)=-b1-(b2+1.0)*bmain(i)*0.5 else i1=ik(2*I),程序解读,b0=b(i1)i1=ik(2*I-1)b1=0.0 if(i1.gt.0)b1=b

39、(i1)il=ik(2*i-2)b2=b(il)i1=ik(2*i-3)b3=0.0 if(il.gt.0)b3=b(il)if(i.eq.n_main)then eta(i,k)=b2 if(kb.eq.i)eta(i,k)=1.0+b2,程序解读,tor(i,k)=b3-b2*bmain(i)*0.5 if(k.eq.ndinf)tor(i,k)=tor(i,k)+bmain(i)*0.5 else eta(i,k)=b2-b0 if(kb.eq.i)eta(i,k)=1.0+b2-b0 tor(i,k)=b3-b1-(b2+b0)*bmain(i)*0.5 if(k.eq.2*i+1)t

40、or(i,k)=tor(i,k)+bmain(i)*0.5 end if end if end do!transverse joit bend moment influence line do i=1,nod,程序解读,i1=ik(2*i-1)b2=0.0 if(i1.gt.0)b2=b(i1)tm(i,k)=b2 end doend doendc-subroutine output(iunit,n_main,nod,ndinf,xinf,eta,tor,tm)c-dimension xinf(ndinf),eta(n_main,ndinf),tor(n_main,ndinf)dimensio

41、n tm(nod,ndinf),程序解读,character form*20form=(10(f7.4,lx)write(iunit,(10 x,a)+-write(iunt,(10 x,a)+|write(iunit,(10 x,a)+|C62GJJBL print file|write(iunt,(10 x,a)+|write(iunt,(10 x,a)+-,程序解读,write(iunit,*)write(iunit,(a,i4)+abscissa for load-added node(1 to n_main*2+1):write(iunit,form)(xinf(i),i=1,ndi

42、nf)write(iunit,(/)do i=1,n_main write(iunit,(a,a,i4)+transverse distribution influence line for main,write(iunit,form)(eta(i,j),j=1,ndinf)end do write(iunit,(/)do i=1,n_main,程序解读,write(iunit,(a,i4)+torsion moment influence line for main girder No.=,i write(iunit,form)(tor(i,j),j=1,ndinf)end do write

43、(iunit,(/)do i=1,nod write(iunit,(a,a,i4)+transverse jont(1 to n_main-1)moment influence line for,+No.=,i write(iunit,form)(tm(i,j),j=1,ndinf)end do end,比拟正交异性板法（G-M法）,前面介绍的几种计算荷载横向分布系数的方法（刚性横梁法和刚（铰）接梁法），都有一个共同的特点：就是把全桥视为一系列并排放置的主梁所构成的梁格或梁系结构来进行力学分析。他们的不同之处在于根据各种不同桥梁结构的具体特点，对横向结构的连接刚性作了不同程度的假设。对于由主梁

44、、连续的桥面板和多横隔梁所组成的梁桥，当其宽度与跨度之比较大时，将其比拟为一块矩形平板，作为弹性薄板按古典弹性理论进行分析，即所谓的“比拟正交异性板法”或称“G-M法”,比拟正交异性板法（G-M法）,比拟正交异性板法（G-M法）,基本思路：将主梁和桥面板、横梁体系比拟为一块矩形正交异性平板，按弹性薄板理论求解。该方法适用于宽度较大的由主梁、连续桥面板和多横梁组成的梁式桥。,比拟正交异性板法（G-M法）,比拟正交异性板法（G-M法）,基本理论比拟正交异性板挠曲面方程平衡方程：取混凝土泊桑比0，可得出挠曲面微分方程：式中,比拟正交异性板法（G-M法）,如图3-13所示，将板分为n条宽度为b的板条，

45、在k处单位正弦荷载作用下（图3-14），各板条跨中挠度为w1,w2,wi,wn，则各板条分担的荷载与挠度成正比：根据荷载、挠度、内力的关系,根据内、外力的平衡,比拟正交异性板法（G-M法）,比拟正交异性板法（G-M法）,程序源代码（BASIC）,1000 AY=-1020 INPUTinput filename:;F1Y 输入文件名1030 INPUToutput filename:;F2y 输出文件名1040 DIM YE(350),Y(15),WIJ(15,350),ETAIJ(15,350),WB(350),MC(4),MCQ(15),MCG(15),MCL(15),MC R(15)10

46、50 OPENI,#1,F1Y 打开文件1060 OPENO,#2,F2Y 打开文件1070 INPUT#1,N,DX 读入n,dx1080 INPUT#1,L,LB,RB,JX,JY,JTXY 读入：跨度，单主梁宽，行车道宽，人行道宽，Jk,Jy,（JTX+JTY）1085 MM=7:DE=1:G=0.425:DB=LB*N/2 取7项级数（1,3,5,7）1120 CLOSE#1 关闭文件1130 NN=INT(2*DB/DX+0.1)+1 插值系数nn1140 NCD=INT(XB+0.3)/3.1 车道数1150 IF NCD4 THEN NCD=41160 PRINT#2,AY117

47、0 PRINT#2,N=;N1190 PRINT#2,DX=;DX 输出1200 PRINT#2,L=;L1210 PRINT#2,LB=;LB1220 PRINT#2,RB=RB1230 PRINT#2,XB=;XB;(;NCD;EANE)1260 PRINT#2,JX=;JX,程序源代码（BASIC）,1270 PRINT#2,JY=;JY1280 PRINT#2,JTX+JTY=;JTXY1290 PRINT#2,Coordinates of beam location:PRINT#2,No.,Y1300 FOR I=1 TO N1305 Y(I)=-DB+LB*(I-1/2)各主梁Y坐

48、标1310 PRINT#2,I,Y(I)1320 NEXT I1330 PRINT#2,AY1350 PI=3.1415926#1370 YE(1)=-DB1380 FOR I=2 TO NN1390 YE(I)=YE(1)+DX*(I-1)插值点坐标1400 NEXT I1420 ALPHA=G*JTXY/（2*DE*SQR（JX*JY）1425 A=SQR（0.5*(1-ALPHA)）1430 B=SQR(0.5(1+AL;PHA)1440 THETA=DB*(JX/JY)0.25/L:OMEGA=THETA*PI/DB 1450 PHI=THETA*PI*A:ETA=THETA*PI*B

49、,1460 GAMMA=THETA*A:DELTA=THETA*B,1480 DEF FNSH(Z)=(EXP(Z)+EXP(-Z)/2 ch(z)1500 DM=(2*A*A-2*B*B+1)*B*FNSH(PHI)*FNCH(PHI)-(2*A*A-2*B*B-1)*A*SIN(ETA)*COS(ETA)M,程序源代码（BASIC）,1510 DN=(2*A*A-2*B*B+1)*B*FNSH(PHI)*FNCH(PHI)+(2*A*A-2*B*B-1)*A*SIN(ETA)*COS(ETA)N1530 C1=A*SIN(ETA)-B*COS(ETA)1540 C2=A*SIN(ETA)+

50、B*COS(ETA)1550 C3=A*A*SIN(ETA)+2*A*B*COS(ETA)-B*B*SIN(ETA)1560 C4=A*A*COS(ETA)+2*A*B*SIN(ETA-B*B*COS(ETA)1580 FOR K=1 TO NN 对插值点循环1600 E=YE(K):PSI=PI*E/DB ek,1610 D1=GAMMA*PSI:D2=DELTA*PSI d1,d21630 CC=FNCH(D1)*COS(D2):SS=FNSH(D1)*SIN(D2)s1,s21640 SC=FNSH(D1)*COS(D2):CS=FNCH(D1)*SIN(D2)s3,s41650 DP=