橋梁荷載橫向分布的計算講.ppt

橋梁結(jié)構(gòu)分析的計算機(jī)方法,橋梁荷載橫向分布的計算,橫向分布定義,整體橋梁結(jié)構(gòu)應(yīng)采用影響面加載計算最不利荷載,橫向分布定義,為簡化計算，采用近似影響面來加載 近似影響面縱橫方向分別相似,軸重,軸重與輪重的關(guān)系,各縱向影響線比例關(guān)系,橫向分布定義,橫向分布定義,概念的正確性 內(nèi)力分布、荷載分布 在那種情況下是正確 橫向分布系數(shù)沿縱向的變化 彎矩 剪力,剛性橫梁法,剛性橫梁法,偏心受壓法、修正偏心受壓法 基本假設(shè) 橫梁剛度無窮大,ej: 到縱軸的偏心距,剛性橫梁法,由橋梁工程課程知，橫向影響線（直線）的計算公式如下，其中n為主梁數(shù)： 為考慮抗扭剛度影響的修正系數(shù)：,剛性橫梁法,當(dāng)各主梁截面相同，且間距相等時： 式中：a, I, It為每個主梁的間距、抗彎慣矩和抗扭慣矩。,剛性橫梁法,橫向分布系數(shù)計算公式 2004公路橋梁規(guī)范的荷載如何考慮？,剛性橫梁法,計算程序 程序功能：計算第i_beam主梁在itype荷載下的橫向分布系數(shù)co。,剛性橫梁法,變量說明,控制變量 Ib計算類型指數(shù)，等于0時各梁截面相同，間距相等；等于1時為一般情況。 Itype荷載類型指數(shù)，等于1、2、3時分別計算汽車荷載、掛車荷載和人群荷載。 N_main主梁片數(shù) i_beam要計算橫向分布系數(shù)的主梁號（從左到右） N_lane車道數(shù) Wid_cro人行道寬度,變量說明,wid_lan車道寬 Al簡支梁的跨度 Ag材料的剪切模量 Ae材料的彈性模量 控制數(shù)組 Xinf(n_main)-影響線的橫坐標(biāo) Aii(n_main)各主梁抗彎慣矩數(shù)組 Aiti(n_main) 各主梁抗扭慣矩數(shù)組,程序源代碼,c c c Main program c c dimension 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) 去后綴,程序源代碼,out=dat call fname(dat,.dat) 加后綴.dat call fname(out,.out) 加后綴.out open(10,file=dat,status=unknown) open(11,file=out,status=unknown) read(10,*)line 標(biāo)題 read(10,*)itype,n_main,I_beam,n_lane,wid_cro,wid_lan,程序源代碼,read(10,*) (xinf(i)=1,n_main) 影響線橫坐標(biāo) read(10,*)ib 控制變量，=0各梁同，=1不同 write(*,10) 10 format(/lx,=input data=:) write(*,(4i7,2f10.4) itype, n_main, i_beam, n_lane, wid_cro, wid_lan write(*,(5f10.4) (xinf(i),i=1,n_main) write(*,(i7)ib,程序源代碼,if (ib.eq.1)then read(10,*)al,ag,ae 跨度l,G,E read(10,*) (aii(i),i=1,n_main) Ii read(10,*) (aiti(i),i=1,n_main) Iti do i=1,n_main aai(i)=abs(xinf(i) ai end do,程序源代碼,else if (ib.eq.0) then do i=1,n_main aii(i)=1.0 aiti(i)=1.0 end do dx=abs(xinf(n_main)-xinf(1)/float(n_main-1)相鄰梁距離a do i=1,n_main aai(i)=abs(xinf(i)/dx ej/a end do end if,程序源代碼,call inf (i_beam, ib, n_main, aai, aii, aiti, al, ag, ae, xinf, yinf) 求影響線yinf call co_eff(itype, n_lane, n_main, wid_lan, wid_cro, xinf, yinf ,co2 ,co ,ilane) 求分布系數(shù)co write(11,(10x,a) -,程序源代碼,write(11,(10x,a) | write(11,(10x,a) | C61GXHL print file write(11,(10x,a) | write(11,(10x,a) - write(11,*) write(11,100) 100 format(/10x,=rigid crossbeam influence line:) write(11,200) 200 format(10x,=xinf(1-n_main):) write(11,300) (xinf(i),i=1,n_main),程序源代碼,300 format(10x,5f12.6) write(11,400) 400 format(10x,=yinf(1-n_main):) write(11,500) (yinf(i),i=1,n_main) 500 format(10x,5f12.6) write(11,600) 600 format(10x,=crosswise distribution coefficient:) write(11,700)co 700 format(10x,f12.6),程序源代碼,close(10) close(11) write(*,(/) write(*,a) The program C61GXHL over! end C -,程序源代碼,subroutine cname(c) c - character*12 c do 10 i=12,1,-1 if (c(i:i).eq.)then ip=i goto 30 endif,程序源代碼,10 continue 30 if(ip.gt.1)then do 40 i=ip,12 c(i:i)= endif end c -,程序源代碼,subroutine fname(c,ce) c - character c*12,ce*4 k=0 do 10 i=1,8 if(c(i:i).eq.)then k=I goto 15 endif 10 continue,程序源代碼,if(k.eq.0)k=9 do 20i=k,12 20 c(i:i)= do 25 i=k,k+3 j=i-k+1 25 c(i:i)=ce(j:j) end,程序源代碼,c - subroutine lin_ins (xl,y1,x2,y2,x,y) 線性內(nèi)插 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) 求ij c - ! ETAij dimension aai(n_main), aii(n_main), aiti(n_main), xinf(n_main),程序源代碼,if (xinf(i_beam)*xinf(j).ge.0.0)then sign=1.0 else sign=-1.0 end if q1=0.0 q2=0.0 do i=1,n_main q1=q1+aii(i) q2=q2+aai(i)*aai(i)*aii(i) end do,I,J梁是否同側(cè)判斷,程序源代碼,call pass(0，177,q1,q2,aai(1),aii(1) if(ib.eq.0)then b=1.0 else 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 if etaij=aii(i_beam)*q1+sign*aai(j)*aai(i_beam)*aii(i_beam)*b/q2 vj end,是否考慮抗扭剛度的影響,程序源代碼,c - subroutine inf (i_beam,ib,n_main,aai,aii,aiti,al,ag,ae, xinf,yinf) 求橫向影響線 c - ! influence line by rigid crossbeam dimension 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,aiti,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)=yinfn do i=2,n_main-1 call lin_ins (xinf(1),yinf(1),xinf(n_main), yinf(n_main), xinf(i),yinf(i) 線性內(nèi)插求中間各值 end do end,n號梁,程序源代碼,c - subroutine co_eff(itype,n_lane,n_main,wid_lan,wid_cro,xinf, yinf, co2, co,ilane) 求橫向分布系數(shù) c - ! lateral(or crosswise)distribution coefficient dimension xinf(n_main),yinf(n_main),程序源代碼,if(yinf(n_main).gt.yinf(1) then sign=1.0 else sign=-1.0 end if co=0.0 if (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(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列車的分布系數(shù) end if,程序源代碼,if(co.lt.qq)then co=qq*0.5 求co ilane=I end if x=x-1.3*sign 車列間距 end do end if if(itype.eq.2) then x=(wid_lan*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.25 end 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(xinf(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 end subroutine pass(iunit,line,v1,v2,v3,v4) 輸出 character form*20 form=(4(lx,f12.6) if (iunit.gt.0) then write iunit,(a,i6) =pass line=,line,程序源代碼,write (iunit,form)v1,v2,v3,v4 else write(*,(a,i6) =pass line=,line write(*,form)v1,v2,v3,v4 end if end,問題思考,對于整體箱梁（型梁）如果計算其橫向分布系數(shù) 大跨度橋梁中橫向分布的考慮 偏載系數(shù) 多梁式T型梁橋按單根梁計算和按整體計算各自的優(yōu)缺點,剛（鉸）接梁（板）法,將多梁式橋梁簡化為數(shù)根并列而相互間橫向剛（鉸）接的狹長板（梁） 各主梁接縫間傳遞剪力、彎矩、水平壓力、水平剪力 用半波正弦荷載作用在某一板上，計算各板（梁）間的力分配關(guān)系,剛（鉸）接梁（板）法,剛（鉸）接梁（板）法,以各結(jié)點處未知力Mi，Qi為未知數(shù)，列力法方程求解 設(shè): 則i 結(jié)點的力法方程為（協(xié)調(diào)方程）：,對于i=2, 3, , m-2來說,剛（鉸）接梁（板）法,i=1時，無 項；i=m-1時，無 項。 以上各式中下標(biāo)意義： 第一個為所考慮的未知數(shù)序號 第二個為產(chǎn)生影響因素的未知數(shù)序號 對每個結(jié)點列方程，得力法正則方程：,剛（鉸）接梁（板）法,柔度矩陣A 的計算 注意下標(biāo)，偶數(shù)號為剪力，奇數(shù)號為彎矩。,剛（鉸）接梁（板）法,結(jié)點處的基本位移模式 注意圖中位移符號,剛（鉸）接梁（板）法,剛（鉸）接梁（板）法,載變位矩陣B的計算,P=1作用在k梁中線 P=1作用在k梁右端,P=1作用在1#梁左端,剛（鉸）接梁（板）法,橫向影響線 令P=1分別作用于1# 梁左端、中間各結(jié)點、各梁中線、m# 梁右端，共2m+1個位置，則可求得橫向影響線。,剛（鉸）接梁（板）法,主梁扭矩影響線,i=1,i=2, 3, , m-1,i=m,計算機(jī)程序,圖 3-11,圖 3-12,主程序框圖 調(diào)用關(guān)系框圖,變量說明,控制變量 N_main主梁根數(shù) Nod節(jié)點數(shù)，nod=n_main-1 Ndinf影響線點數(shù)，ndinf=2*n_main+1 Inform1為0時各主梁截面相同，為1時各主梁截面不同 Inform2為0時各主梁間同時有鉸結(jié)和剛接，為1時各主梁間全部鉸結(jié)，為2時各主梁間全部剛結(jié) Span主梁跨度 Gg主梁材料剪切模量 Ee-主梁材料彈性模量,變量說明,控制數(shù)組 Inform(nod)主梁聯(lián)結(jié)方式數(shù)組 Bi(n_main)各主梁抗彎慣矩數(shù)組 Ti(n_main)各主梁抗扭慣矩數(shù)組 Bmain(n_main)主梁寬度（兩個節(jié)點間的距離）數(shù)組 Rmain(n_main)主梁懸臂板長度數(shù)組 Bic(n_main)主梁懸臂板（沿梁長方向單位長度）抗彎慣矩數(shù)組 Xinf(ndinf)影響線橫向坐標(biāo)數(shù)組 Eta(n_main,ndinf) 存放主梁橫向分布影響線，數(shù)組每一行為一片主梁的橫向分布影響線 Tor(n_main,ndinf) 存放主梁扭矩影響線 Tm(nod,ndif)存放節(jié)點處橫向玩具影響線,剛(鉸)接梁(板)實例,剛(鉸)接梁(板)實例,程序解讀,程序源代碼： program c62gxhl !-transverse distribution influence line dimension 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*80 write(*,(/) write(*,*) write(*,(lx,a)Enter the name of datafile: read(*,(a)dat,程序解讀,write(*,(/) call cname(dat) out=dat call 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, inform2 read(10,*)span,gg,ee nod=n_main-1 結(jié)點數(shù) ndinf=2*n_main+1 =1作用點數(shù),程序解讀,m=n_main if (inform1.eq.0) m=1 各主梁同 do i=1,m read (10,*)bi(i),ti(i),bmain(i), rmain(i),bic(i) end do if (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) then read(10,*) (inform(i),i=1,nod) 剛、鉸接混合 end if if (inform2.eq.1) then do i=1,nod inform(i)=1 鉸接 end do end if if (inform2.eq.2) then do i=1,nod,程序解讀,inform(i)=2 end 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! end,程序解讀,c c c Subroutines for C62GJJBL c c c c Computer room/1995.05.17 c c c - bubroutine cname c -,程序解讀,character*12c do 10i=12,1,-1 if (c(i:i).eq.) then ip=i goto 30 endif 10 continue 30 if (ip.gt.1) then do 40i=ip,12 40 c(i:i)= endif end,程序解讀,c - subroutine fname(c,ce) c - character c*12,ce*4 k=0 do 10 i=1,8 if(c(i:i).eq.) then k=i goto 15 endif 10 continue 15 if(k.eq.0)k=9,程序解讀,do 20 i=k,12 20 c(i:i)= do 25 i=k,k+3 25 c(i:i)=ce(j:j) end c - 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=1,j1 ki=nband+k-i kj=nband+k-j qq=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 if end do end do end do end c -,程序解讀,subroutine back (neq0,neq,nband,gk,b) c - ! Return substitution and forward substitution dimension 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(j,ij)*b(j) end do b(i)=b(i)-qq end do end c - 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) dimension tm(nod,ndinf),xinf(ndinf),gk(2*nod,4),b(2*nod) dimension c(2,4),d(5) nband=4 neq0=2*nod do ii=1,ned do k=1,2 i=(ii-1)*2+k ik(i)=1 if(rm(ii).eq.1)ik(i)=0 end do,程序解讀,end do k=0 do i=1,neq0 k=k+ik(i) if(ik(i).ne.0)ik(i)=k end do neq=k pi1=3.1415926 pi2=pi1*pi1 pi4=pi2*pi2 s2=span*span do ii=1,ned,程序解讀,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+fil+ci1+cti1 c(2,3)=fi*bmain(ii)*0.5+ui-fil*bmain(ii+1)*0.5-uil c(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*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,neq,nband,gk) qq=0.0 do i=1,n_main qq=qq+bmain(i) end do xinf(1)=-qq*0.5 do 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.5 end 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)/pi4 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.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 i=1,n_main if (i.eq.1) then,程序解讀,i2=ik(2*i) b2=b(i2) eta(i,k)=1.0-b2 if (kb.ne.1) eta(i,k)=-b2 i1=ik(1) b1=0.0 if (i1.gt.0)b1=b(i1) tor(i,k)=-b1-b2*bmain(i)*0.5 if(k.eq.1)tor(i,k)=-b1-(b2+1.0)*bmain(i)*0.5 if(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(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) tor(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 do end do end c - subroutine output (iunit,n_main,nod,ndinf,xinf,eta,tor,tm) c - dimension xinf(ndinf), eta(n_main,ndinf),tor(n_main,ndinf) dimension tm(nod,ndinf),程序解讀,character form*20 form=(10(f7.4,lx) write(iunit,(10x,a) + - write(iunt,(10x,a) + | | write(iunit,(10x,a) + | C62GJJBL print file | write(iunt,(10x,a) + | | write(iunt,(10x,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,ndinf) 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 (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法）,前面介紹的幾種計算荷載橫向分布系數(shù)的方法（剛性橫梁法和剛（鉸）接梁法），都有一個共同的特點： 就是把全橋視為一系列并排放置的主梁所構(gòu)成的梁格或梁系結(jié)構(gòu)來進(jìn)行力學(xué)分析。 他們的不同之處在于根據(jù)各種不同橋梁結(jié)構(gòu)的具體特點，對橫向結(jié)構(gòu)的連接剛性作了不同程度的假設(shè)。 對于由主梁、連續(xù)的橋面板和多橫隔梁所組成的梁橋，當(dāng)其寬度與跨度之比較大時，將其比擬為一塊矩形平板，作為彈性薄板按古典彈性理論進(jìn)行分析，即所謂的“比擬正交異性板法”或稱“G-M法 ”,比擬正交異性板法（G-M法）,比擬正交異性板法（G-M法）,基本思路：將主梁和橋面板、橫梁體系比擬為一塊矩形正交異性平板，按彈性薄板理論求解。 該方法適用于寬度較大的由主梁、連續(xù)橋面板和多橫梁組成的梁式橋。,比擬正交異性板法（G-M法）,比擬正交異性板法（G-M法）,基本理論 比擬正交異性板撓曲面方程 平衡方程： 取混凝土泊桑比0，可得出撓曲面微分方程： 式中,比擬正交異性板法（G-M法）,如圖3-13所示，將板分為n條寬度為b的板條，在k處單位正弦荷載作用下（圖3-14），各板條跨中撓度為w1, w2, , wi, , wn，則各板條分擔(dān)的荷載與撓度成正比：根據(jù)荷載、撓度、內(nèi)力的關(guān)系,根據(jù)內(nèi)、外力的平衡,比擬正交異性板法（G-M法）,比擬正交異性板法（G-M法）,程序源代碼（BASIC）,1000 AY=“-“ 1020 INPUT“input filename: “; F1Y 輸入文件名 1030 INPUT“output 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) 1050 OPEN“I“, #1, F1Y 打開文件 1060 OPEN“O“, #2, F2Y 打開文件 1070 INPUT #1, N, DX 讀入n, dx 1080 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項級數(shù)（1,3,5,7） 1120 CLOSE # 1 關(guān)閉文件 1130 NN=INT(2*DB/DX+0.1)+1 插值系數(shù)nn 1140 NCD=INT(XB+0.3)/3.1 車道數(shù) 1150 IF NCD4 THEN NCD=4 1160 PRINT # 2, AY 1170 PRINT # 2, “N=“; N 1190 PRINT # 2, “DX=“; DX 輸出 1200 PRINT # 2, “L=“; L 1210 PRINT # 2, “LB=“; LB 1220 PRINT # 2, “RB=“RB 1230 PRINT # 2, “XB=“; XB;“ (“; NCD; “EANE)“ 1260 PRINT # 2, “JX=“; JX,程序源代碼（BASIC）,1270 PRINT # 2, “JY=“; JY 1280 PRINT # 2, “JTX+JTY=“; JTXY 1290 PRINT # 2, “Coordinates of beam location“: PRINT # 2, “No.“, “Y“ 1300 FOR I=1 TO N 1305 Y(I)=-DB+LB*(I-1/2) 各主梁Y坐標(biāo) 1310 PRINT # 2, I, Y(I) 1320 NEXT I 1330 PRINT # 2, AY 1350 PI=3.1415926# 1370 YE(1)=-DB 1380 FOR I=2 TO NN 1390 YE(I)=YE(1)+DX*(I-1) 插值點坐標(biāo) 1400 NEXT I 1420 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 , 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) N 1530 C1=A*SIN (ETA)-B*COS (ETA) 1540 C2=A*SIN (ETA)+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 對插值點循環(huán) 1600 E=YE(K): PSI=PI*E/DB ek, 1610 D1=GAMMA*PSI: D2=DELTA*PSI d1, d2 1630 CC=FNCH(D1)*COS(D2): SS=FNSH(D1)*SIN(D2) s1, s2 1640 SC=FNSH(D1)*COS(D2): CS=FNCH(D