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

桥梁结构分析的计算机方法,桥梁荷载横向分布的计算,横向分布定义,整体桥梁结构应采用影响面加载计算最不利荷载,横向分布定义,为简化计算，采用近似影响面来加载 近似影响面纵横方向分别相似,轴重,轴重与轮重的关系,各纵向影响线比例关系,横向分布定义,横向分布定义,概念的正确性 内力分布、荷载分布 在那种情况下是正确 横向分布系数沿纵向的变化 弯矩 剪力,刚性横梁法,刚性横梁法,偏心受压法、修正偏心受压法 基本假设 横梁刚度无穷大,ej: 到纵轴的偏心距,刚性横梁法,由桥梁工程课程知，横向影响线（直线）的计算公式如下，其中n为主梁数： 为考虑抗扭刚度影响的修正系数：,刚性横梁法,当各主梁截面相同，且间距相等时： 式中：a, 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材料的弹性模量 控制数组 Xinf(n_main)-影响线的横坐标 Aii(n_main)各主梁抗弯惯矩数组 Aiti(n_main) 各主梁抗扭惯矩数组,程序源代码,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 标题 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) 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) 求分布系数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) 线性内插 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梁是否同侧判断,程序源代码,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) 线性内插求中间各值 end do end,n号梁,程序源代码,c - subroutine co_eff(itype,n_lane,n_main,wid_lan,wid_cro,xinf, yinf, co2, co,ilane) 求横向分布系数 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列车的分布系数 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,问题思考,对于整体箱梁（型梁）如果计算其横向分布系数 大跨度桥梁中横向分布的考虑 偏载系数 多梁式T型梁桥按单根梁计算和按整体计算各自的优缺点,刚（铰）接梁（板）法,将多梁式桥梁简化为数根并列而相互间横向刚（铰）接的狭长板（梁） 各主梁接缝间传递剪力、弯矩、水平压力、水平剪力 用半波正弦荷载作用在某一板上，计算各板（梁）间的力分配关系,刚（铰）接梁（板）法,刚（铰）接梁（板）法,以各结点处未知力Mi，Qi为未知数，列力法方程求解 设: 则i 结点的力法方程为（协调方程）：,对于i=2, 3, , m-2来说,刚（铰）接梁（板）法,i=1时，无 项；i=m-1时，无 项。 以上各式中下标意义： 第一个为所考虑的未知数序号 第二个为产生影响因素的未知数序号 对每个结点列方程，得力法正则方程：,刚（铰）接梁（板）法,柔度矩阵A 的计算 注意下标，偶数号为剪力，奇数号为弯矩。,刚（铰）接梁（板）法,结点处的基本位移模式 注意图中位移符号,刚（铰）接梁（板）法,刚（铰）接梁（板）法,载变位矩阵B的计算,P=1作用在k梁中线 P=1作用在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-1 Ndinf影响线点数，ndinf=2*n_main+1 Inform1为0时各主梁截面相同，为1时各主梁截面不同 Inform2为0时各主梁间同时有铰结和刚接，为1时各主梁间全部铰结，为2时各主梁间全部刚结 Span主梁跨度 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)存放节点处横向玩具影响线,刚(铰)接梁(板)实例,刚(铰)接梁(板)实例,程序解读,程序源代码： 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 结点数 ndinf=2*n_main+1 =1作用点数,程序解读,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法）,前面介绍的几种计算荷载横向分布系数的方法（刚性横梁法和刚（铰）接梁法），都有一个共同的特点： 就是把全桥视为一系列并排放置的主梁所构成的梁格或梁系结构来进行力学分析。 他们的不同之处在于根据各种不同桥梁结构的具体特点，对横向结构的连接刚性作了不同程度的假设。 对于由主梁、连续的桥面板和多横隔梁所组成的梁桥，当其宽度与跨度之比较大时，将其比拟为一块矩形平板，作为弹性薄板按古典弹性理论进行分析，即所谓的“比拟正交异性板法”或称“G-M法 ”,比拟正交异性板法（G-M法）,比拟正交异性板法（G-M法）,基本思路：将主梁和桥面板、横梁体系比拟为一块矩形正交异性平板，按弹性薄板理论求解。 该方法适用于宽度较大的由主梁、连续桥面板和多横梁组成的梁式桥。,比拟正交异性板法（G-M法）,比拟正交异性板法（G-M法）,基本理论 比拟正交异性板挠曲面方程 平衡方程： 取混凝土泊桑比0，可得出挠曲面微分方程： 式中,比拟正交异性板法（G-M法）,如图3-13所示，将板分为n条宽度为b的板条，在k处单位正弦荷载作用下（图3-14），各板条跨中挠度为w1, w2, , wi, , wn，则各板条分担的荷载与挠度成正比：根据荷载、挠度、内力的关系,根据内、外力的平衡,比拟正交异性板法（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项级数（1,3,5,7） 1120 CLOSE # 1 关闭文件 1130 NN=INT(2*DB/DX+0.1)+1 插值系数nn 1140 NCD=INT(XB+0.3)/3.1 车道数 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坐标 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) 插值点坐标 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 对插值点循环 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