第3章 桥梁荷载横向分布的计算讲稿-2010.ppt

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

2、相等时：式中：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材料的

3、剪切模量 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 文件名 wr

4、ite(*,(/) 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 控制变量，

5、=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

6、,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

7、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) 2

8、00 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,f

9、12.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 endif10 continue,程序源代码,if(k.eq.0)k=9 do 20i=k,1220 c(i:i)=do 25 i=k,k+3j=i-k+1 25 c(i:i)=ce(j:j)end,