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

横向分布定义l l整体桥梁结构应采用影响面加载计算最不利荷载整体桥梁结构应采用影响面加载计算最不利荷载横向分布定义l l为简化计算，采用近为简化计算，采用近似影响面来加载似影响面来加载l l近似影响面纵横方向近似影响面纵横方向分别相似分别相似轴重轴重与轮重的关系各纵向影响线比例关系各纵向影响线比例关系横向分布定义横向分布定义l l概念的正确性概念的正确性l l内力分布、荷载分布内力分布、荷载分布l l在那种情况下是正确在那种情况下是正确l l横向分布系数沿纵向横向分布系数沿纵向的变化的变化l l弯矩弯矩l l剪力剪力刚性横梁法刚性横梁法l l偏心受压法、修正偏心受压法偏心受压法、修正偏心受压法l l基本假设基本假设l l横梁刚度无穷大横梁刚度无穷大ej:到纵轴的偏心距刚性横梁法l l由桥梁工程课程知，横向影响线（直线）的计算公式如下，其中n为主梁数：l l为考虑抗扭刚度影响的修正系数：刚性横梁法l l当各主梁截面相同，且间距相等时：l l式中：a,I,It为每个主梁的间距、抗弯惯矩和抗扭惯矩。刚性横梁法l l横向分布系数计算公式横向分布系数计算公式 l l20042004公路桥梁规范的荷载如何考虑？公路桥梁规范的荷载如何考虑？刚性横梁法l l计算程序计算程序计算程序计算程序l l程序功能：计算第程序功能：计算第i_beami_beam主梁在主梁在itypeitype荷载下的横向分荷载下的横向分布系数布系数coco。刚性横梁法变量说明l l控制变量l lIbIb计算类型指数，等于计算类型指数，等于0 0时各梁截面相同，间时各梁截面相同，间距相等；等于距相等；等于1 1时为一般情况。时为一般情况。l lItypeItype荷载类型指数，等于荷载类型指数，等于1 1、2 2、3 3时分别计时分别计算汽车荷载、挂车荷载和人群荷载。算汽车荷载、挂车荷载和人群荷载。l lN_mainN_main主梁片数主梁片数l li_beami_beam要计算横向分布系数的主梁号（从左要计算横向分布系数的主梁号（从左到右）到右）l lN_laneN_lane车道数车道数l lWid_croWid_cro人行道宽度人行道宽度变量说明l lwid_lanwid_lan车道宽车道宽l lAlAl简支梁的跨度简支梁的跨度l lAgAg材料的剪切模量材料的剪切模量l lAeAe材料的弹性模量材料的弹性模量l l控制数组l lXinf(n_main)-Xinf(n_main)-影响线的横坐标影响线的横坐标l lAii(n_main)各主梁抗弯惯矩数组l lAiti(n_main)各主梁抗扭惯矩数组程序源代码l l l lc cl lc cl lc Main program c Main program l lc cl lc cl ldimension aai(30),aii(30),aiti(30)dimension aai(30),aii(30),aiti(30)l ldimension xinf(30),yinf(30)dimension xinf(30),yinf(30)l lcharacter dat*12,out*12,line*80character dat*12,out*12,line*80程序源代码l lwrite(*,(/)write(*,(/)l lwrite(*,*)write(*,*)l lwrite(*,(ls,a)Enter the name of datafile:write(*,(ls,a)Enter the name of datafile:l lread(*,(a)dat read(*,(a)dat 文件名文件名l lwrite(*,(/)write(*,(/)l lcall cname(dat)call cname(dat)去后缀去后缀程序源代码l lout=datout=datl lcall fname(dat,.dat)call fname(dat,.dat)加后缀加后缀.dat.datl lcall fname(out,.out)call fname(out,.out)加后缀加后缀.out.outl lopen(10,file=dat,status=unknown)open(10,file=dat,status=unknown)l lopen(11,file=open(11,file=outout,status=unknown),status=unknown)l lread(10,*)line read(10,*)line 标题标题l lread(10,*)itype,n_main,I_beam,n_lane,wid_cro,wid_lanread(10,*)itype,n_main,I_beam,n_lane,wid_cro,wid_lan程序源代码l lread(10,*)(xinf(i)=1,n_main)read(10,*)(xinf(i)=1,n_main)影响线横坐标影响线横坐标l lread(10,*)ib read(10,*)ib 控制变量，控制变量，=0=0各梁同，各梁同，=1=1不同不同l lwrite(*,10)write(*,10)l l1010format(/lx,=input data=:)format(/lx,=input data=:)l lwrite(*,(4i7,2f10.4)write(*,(4i7,2f10.4)itype,itype,n_main,n_main,i_beam,i_beam,n_lane,n_lane,wid_cro,wid_cro,wid_lanwid_lanl lwrite(*,(5f10.4)(xinf(i),i=1,n_main)write(*,(5f10.4)(xinf(i),i=1,n_main)l lwrite(*,(i7)ibwrite(*,(i7)ib程序源代码l lif(ib.eq.1)thenl lread(10,*)al,ag,ae read(10,*)al,ag,ae 跨度跨度l,G,El,G,El lread(10,*)(aii(i),i=1,n_main)Iiread(10,*)(aii(i),i=1,n_main)Iil lread(10,*)(aiti(i),i=1,n_main)Itiread(10,*)(aiti(i),i=1,n_main)Itil ldo i=1,n_maindo i=1,n_mainl laai(i)=abs(xinf(i)aiaai(i)=abs(xinf(i)ail lend doend do程序源代码l lelse if(ib.eq.0)thenl ldo i=1,n_maindo i=1,n_mainl laii(i)=1.0aii(i)=1.0l laiait ti(i)=i(i)=1.01.0l lend doend dol ldx=abs(xinf(n_main)dx=abs(xinf(n_main)-xinf(1)/float(n_mainxinf(1)/float(n_main-1)1)相邻梁距离相邻梁距离a al ldo i=1,n_mainl laai(i)=abs(xinf(i)/dx ej/aaai(i)=abs(xinf(i)/dx ej/al lend dol lend if程序源代码l lcall inf(i_beam,ib,n_main,aai,aii,aiti,al,ag,ae,xinf,yinf)求影响线yinfl lcall co_eff(itype,n_lane,n_main,wid_lan,wid_cro,xinf,yinf,co2,co,ilane)求分布系数col lwrite(11,(10 x,a)-程序源代码l lwrite(11,(10 x,a)write(11,(10 x,a)|l lwrite(11,(10 x,a)write(11,(10 x,a)|C61GXHL print file|C61GXHL print filel lwrite(11,(10 x,a)write(11,(10 x,a)|l lwrite(11,(10 x,a)-write(11,(10 x,a)-l lwrite(11,*)write(11,*)l lwrite(11,100)write(11,100)l l100100 format(/10 x,=rigid crossbeam influence line:)format(/10 x,=rigid crossbeam influence line:)l lwrite(11,200)write(11,200)l l200200 format(10 x,=xinf(1-n_main):)format(10 x,=xinf(1-n_main):)l lwrite(11,300)(xinf(i),i=1,n_main)write(11,300)(xinf(i),i=1,n_main)程序源代码l l300300 format(10 x,5f12.6)format(10 x,5f12.6)l lwrite(11,400)write(11,400)l l400400 format(10 x,=yinf(1-n_main):)format(10 x,=yinf(1-n_main):)l lwrite(11,500)(yinf(i),i=1,n_main)write(11,500)(yinf(i),i=1,n_main)l l500500 format(10 x,5format(10 x,5f f12.6)12.6)l lwrite(11,600)write(11,600)l l600600 format(10 x,=crosswise distribution coefficient:)format(10 x,=crosswise distribution coefficient:)l lwrite(11,700)cowrite(11,700)col l700700 format(10 x,f12.6)format(10 x,f12.6)程序源代码l lclose(10)close(10)l lclose(11)close(11)l lwrite(*,(/)write(*,(/)l lwrite(*,a)The program C61GXHL over!write(*,a)The program C61GXHL over!l lendendl lC -C -程序源代码l lsubroutine cname(c)l lc -l lcharacter*12 cl ldo 10 i=12,1,-1l lif(c(i:i).eq.)thenif(c(i:i).eq.)thenl lip=iip=il lgoto 30goto 30l lendifendif程序源代码l l10 continuel l30 if(ip.gt.1)thenl ldo 40 i=ip,12do 40 i=ip,12l lc(i:i)=c(i:i)=l lendifl lendl lc -程序源代码l lsubroutine fname(c,ce)l lc -l lcharacter c*12,ce*4l lk=0l ldo 10 i=1,8l lif(c(i:i).eq.)thenif(c(i:i).eq.)thenl lk=Ik=Il lgoto 15goto 15l lendifendifl l 10 continue程序源代码l lif(k.eq.0)k=9l ldo 20i=k,12l l 20c(i:i)=l l do 25 i=k,k+3l l j=i-k+1l l25 c(i:i)=ce(j:j)l l end程序源代码l lc -l lsubroutine lin_ins(xl,y1,x2,y2,x,y)线性内插l lc -l l！linear insertionl l!y=y1*(x-x2)/(x1-x2)+y2*(x-x1)/(x2-x1)l ly=(y1*(x-x2)-y2*(x-x1)/(x1-x2)l lend程序源代码l lc -c -l lsubroutine eta(ib,n_main,i_beam,j,aai,aii,aiti,al,ag,subroutine eta(ib,n_main,i_beam,j,aai,aii,aiti,al,ag,ae,xinf,etaijae,xinf,etaij)求求ijijl lc -c -l l!ETAij!ETAijl ldimension aai(n_main),aii(n_main),aiti(n_main),dimension aai(n_main),aii(n_main),aiti(n_main),xinf(n_main)xinf(n_main)程序源代码l lif(xinfif(xinf(i_beam)*xinf(j).ge.0.0)theni_beam)*xinf(j).ge.0.0)thenl lsign=1.0sign=1.0l lelseelsel lsign=-1.0sign=-1.0l lend ifend ifl lq1=0.0q1=0.0l lq2=0.0q2=0.0l ldo i=1,n_maindo i=1,n_mainl l q1=q1+aii(i)q1=q1+aii(i)l l q2=q2+aai(i)*aai(i)*aii(i)q2=q2+aai(i)*aai(i)*aii(i)l lend doend doI,J梁是否同侧判断程序源代码l lcall pass(0call pass(0，177,q1,q2,aai(1),aii(1)177,q1,q2,aai(1),aii(1)l lif(ib.eq.0)thenif(ib.eq.0)thenl l b=1.0 b=1.0l lelseelsel l q3=0.0q3=0.0l l do i=1,n_maindo i=1,n_mainl l q3=q3+aiti(i)Itiq3=q3+aiti(i)Itil l end doend dol l b=1.0/(1.0+ag*al*al*q3/(12.0*ae*q2)b=1.0/(1.0+ag*al*al*q3/(12.0*ae*q2)l lend ifend ifl letaij=aii(i_beam)etaij=aii(i_beam)*q1+sign*aai(j)*aai(i_beam)*aii(i_beam)q1+sign*aai(j)*aai(i_beam)*aii(i_beam)*b/q2 vj*b/q2 vjl lendend是否考虑抗扭刚度的影响程序源代码l lc -c -l lsubroutine inf(i_beam,ib,n_main,aai,aii,aiti,al,ag,ae,subroutine inf(i_beam,ib,n_main,aai,aii,aiti,al,ag,ae,xinf,yinf)xinf,yinf)求横向影响线求横向影响线l lc -c -l l!influence line by rigid crossbeam!influence line by rigid crossbeaml ldimension aai(n_main),aii(n_main),aiti(n_main)dimension aai(n_main),aii(n_main),aiti(n_main)l ldimension xinf(n_main),yinf(n_main)dimension xinf(n_main),yinf(n_main)l lcall eta(ib,n_main,i_beam,1,aai,aii,aiti,al,ag,ae,call eta(ib,n_main,i_beam,1,aai,aii,aiti,al,ag,ae,xinf,yinf1)xinf,yinf1)求求ijijl l yinf(1)=yinf1 yinf(1)=yinf11号梁程序源代码l lcall eta(ib,n_main,i_beam,n_main,aai,aii,aiti,al,ag,ae,call eta(ib,n_main,i_beam,n_main,aai,aii,aiti,al,ag,ae,xinf,yinfn)xinf,yinfn)l lyinf(n_main)=yinfnyinf(n_main)=yinfnl ldo i=2,n_main-1do i=2,n_main-1l l call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),xinf(i),yinf(i)xinf(i),yinf(i)线性内插求中间各值线性内插求中间各值l l end doend dol l end endn号梁程序源代码l lc -c -l lsubroutine subroutine co_eff(itype,n_lane,n_main,wid_lan,wid_cro,xinf,yinf,co_eff(itype,n_lane,n_main,wid_lan,wid_cro,xinf,yinf,co2,co,ilane)co2,co,ilane)求横向分布系数求横向分布系数l lc -c -l l!lateral(or crosswise)distribution coefficient!lateral(or crosswise)distribution coefficientl ldimension xinf(n_main),yinf(n_main)dimension xinf(n_main),yinf(n_main)程序源代码l lif(yinf(n_main).gt.yinf(1)thenif(yinf(n_main).gt.yinf(1)thenl l sign=1.0sign=1.0l lelseelsel l sign=-1.0sign=-1.0l lend ifend ifl lco=0.0co=0.0l lif(itype.eq.1)then if(itype.eq.1)then 汽车汽车l l qq=0.0qq=0.0程序源代码l l x=(wid_lan*0.5-0.5)*signx=(wid_lan*0.5-0.5)*signl l do i=1,n_lanedo i=1,n_lanel l call lin_ins(xinf(1),yinf(1),xinf(n_main),call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)yinf(n_main),x,v)求最靠边轮位处影响线值求最靠边轮位处影响线值l l qq=qq+vqq=qq+vl l x=x-1.8*signx=x-1.8*signl l call lin_ins(xinf(1),yinf(1),xinf(n_main),call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)yinf(n_main),x,v)求另一轮位处影响线值求另一轮位处影响线值l l qq=qq+vqq=qq+vl l if(n_lane.gt.2.and.i.eq.2)thenif(n_lane.gt.2.and.i.eq.2)thenl l co2=co 2co2=co 2列车的分布系数列车的分布系数l l end ifend if程序源代码l l if(co.lt.qq)thenif(co.lt.qq)thenl l co=qq*0.5 co=qq*0.5 求求cocol l ilane=Iilane=Il l end ifend ifl l x=x-1.3*sign x=x-1.3*sign 车列间距车列间距l l end doend dol lend ifend ifl lif(itype.eq.2)thenif(itype.eq.2)thenl l x=(wid_lan*0.5-1.0)*sign x=(wid_lan*0.5-1.0)*sign 挂车挂车l l call lin_ins(xinf(1),yinf(1),xinf(n_main),call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)yinf(n_main),x,v)程序源代码l l co=co+vco=co+vl l do i=1,3do i=1,3l l x=x-0.9*signx=x-0.9*signl l call lin_ins(xinf(1),yinf(1),xinf(n_main),call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)yinf(n_main),x,v)l l co=co+vco=co+vl l end doend dol l co=co*0.25co=co*0.25l lend ifend if程序源代码l lif(itype.eq.3)thenif(itype.eq.3)thenl l if(yinf(1).ge.0.0)then if(yinf(1).ge.0.0)then 人群人群l l x=-wid_lan*0.5-0.5*wid_crox=-wid_lan*0.5-0.5*wid_crol l call lin_ins(xinf(1).yinf(1),xinf(n_main),call lin_ins(xinf(1).yinf(1),xinf(n_main),yinf(n_main),x,v)yinf(n_main),x,v)按人行道中心线位置计算按人行道中心线位置计算l l co=co+vco=co+vl l end ifend ifl l if(yinf(n_main).gt.0.0)thenif(yinf(n_main).gt.0.0)thenl l x=wid_lan*0.5+0.5*wid_crox=wid_lan*0.5+0.5*wid_crol l call lin_ins(xinf(1),yinf(1),xinf(n_main),call lin_ins(xinf(1),yinf(1),xinf(n_main),yinf(n_main),x,v)yinf(n_main),x,v)l l co=co+vco=co+vl l end ifend if程序源代码l l end if end ifl l end endl lsubroutine pass(iunit,line,v1,v2,v3,v4)subroutine pass(iunit,line,v1,v2,v3,v4)输出输出l lcharacter form*20character form*20l lform=(4(lx,f12.6)form=(4(lx,f12.6)l lif(iunit.gt.0)thenif(iunit.gt.0)thenl l write iunit,(a,i6)=pass line=,linewrite iunit,(a,i6)=pass line=,line程序源代码l l write(iunit,form)v1,v2,v3,v4write(iunit,form)v1,v2,v3,v4l lelseelsel l write(*,(a,i6)=pass line=,line write(*,(a,i6)=pass line=,linel l write(*,form)v1,v2,v3,v4 write(*,form)v1,v2,v3,v4l lend ifend ifl lendend问题思考l l对于整体箱梁（型梁）如果计算其横向分布系数l l大跨度桥梁中横向分布的考虑大跨度桥梁中横向分布的考虑l l偏载系数偏载系数l l多梁式T型梁桥按单根梁计算和按整体计算各自的优缺点刚（铰）接梁（板）法l l将多梁式桥梁简化为数根并列而相互间横向刚（铰）接的狭长板（梁）l l各主梁接缝间传递剪力、弯矩、水平压力、水平剪力l l用半波正弦荷载作用在某一板上，计算各板（梁）间的力分配关系刚（铰）接梁（板）法 刚（铰）接梁（板）法l l以各结点处未知力Mi，Qi为未知数，列力法方程求解 l l设:l l则i 结点的力法方程为（协调方程）：对于对于i=2,3,m-2i=2,3,m-2来说来说刚（铰）接梁（板）法l li=1时，无 项；i=m-1时，无 项。l l以上各式中下标意义：l l第一个为所考虑的未知数序号第一个为所考虑的未知数序号l l第二个为产生影响因素的未知数序号第二个为产生影响因素的未知数序号l l对每个结点列方程，得力法正则方程：刚（铰）接梁（板）法l l柔度矩阵柔度矩阵A A 的计算的计算l l注意下标，偶数号为剪力，奇数号为弯矩。注意下标，偶数号为剪力，奇数号为弯矩。刚（铰）接梁（板）法l l结点处的基本位移模式结点处的基本位移模式l l注意图中位移符号注意图中位移符号 刚（铰）接梁（板）法刚（铰）接梁（板）法载变位矩阵载变位矩阵B的计算的计算 P=1作用在k梁中线 P=1作用在k梁右端 P=1作用在1#梁左端 刚（铰）接梁（板）法l l横向影响线 l l令令P=1P=1分别作用于分别作用于1#1#梁左端、中间各结点、各梁左端、中间各结点、各梁中线、梁中线、m#m#梁右端，共梁右端，共2m+12m+1个位置，则可求个位置，则可求得横向影响线。得横向影响线。刚（铰）接梁（板）法l l主梁扭矩影响线 i=1 i=2,3,m-1i=m 计算机程序图 3-11图 3-12主程序框图 调用关系框图变量说明l l控制变量控制变量l lN_mainN_main主梁根数主梁根数l lNodNod节点数，节点数，nod=n_main-1nod=n_main-1l lNdinfNdinf影响线点数，影响线点数，ndinf=2*n_main+1ndinf=2*n_main+1l lInform1Inform1为为0 0时各主梁截面相同，为时各主梁截面相同，为1 1时各主梁截面不时各主梁截面不同同l lInform2Inform2为为0 0时各主梁间同时有铰结和刚接，为时各主梁间同时有铰结和刚接，为1 1时各时各主梁间全部铰结，为主梁间全部铰结，为2 2时各主梁间全部刚结时各主梁间全部刚结l lSpanSpan主梁跨度主梁跨度l lGgGg主梁材料剪切模量主梁材料剪切模量l lEe-Ee-主梁材料弹性模量主梁材料弹性模量变量说明l l控制数组控制数组l lInform(nod)Inform(nod)主梁联结方式数组主梁联结方式数组l lBi(n_main)Bi(n_main)各主梁抗弯惯矩数组各主梁抗弯惯矩数组l lTi(n_main)Ti(n_main)各主梁抗扭惯矩数组各主梁抗扭惯矩数组l lBmain(n_main)Bmain(n_main)主梁宽度（两个节点间的距离）数组主梁宽度（两个节点间的距离）数组l lRmain(n_main)Rmain(n_main)主梁悬臂板长度数组主梁悬臂板长度数组l lBic(n_main)Bic(n_main)主梁悬臂板（沿梁长方向单位长度）抗弯惯矩数组主梁悬臂板（沿梁长方向单位长度）抗弯惯矩数组l lXinf(ndinf)Xinf(ndinf)影响线横向坐标数组影响线横向坐标数组l lEta(n_main,ndinf)Eta(n_main,ndinf)存放主梁横向分布影响线，数组每一行为一存放主梁横向分布影响线，数组每一行为一片主梁的横向分布影响线片主梁的横向分布影响线l lTor(n_main,ndinf)Tor(n_main,ndinf)存放主梁扭矩影响线存放主梁扭矩影响线l lTm(nod,ndif)Tm(nod,ndif)存放节点处横向玩具影响线存放节点处横向玩具影响线刚(铰)接梁(板)实例刚(铰)接梁(板)实例程序解读l l程序源代码：程序源代码：l lprogram c62gxhlprogram c62gxhll l!-transverse distribution influence line!-transverse distribution influence linel ldimension inform(39),ik(78)dimension inform(39),ik(78)l ldimension bmain(40),rmain(40),bi(40),ti(40),bic(40),dimension bmain(40),rmain(40),bi(40),ti(40),bic(40),xinf(81)xinf(81)l ldimension gk(78,4),b(78),eta(40,81),tm(39,81)dimension gk(78,4),b(78),eta(40,81),tm(39,81)l lcharacter dat*12,out*12,line*80character dat*12,out*12,line*80l lwrite(*,(/)write(*,(/)l lwrite(*,*)write(*,*)l lwrite(*,(lx,a)Enter the name of datafile:write(*,(lx,a)Enter the name of datafile:l lread(*,(a)dat read(*,(a)dat 程序解读l lwrite(*,(/)write(*,(/)l lcall cname(dat)call cname(dat)l lout=datout=datl lcall fname(dat,.dat)call fname(dat,.dat)l lcall fname(out,.out)call fname(out,.out)l lopen(10,file=dat,status=unknown)open(10,file=dat,status=unknown)l lread(10,*)line read(10,*)line l lopen(11,file=dat,status=unknown)open(11,file=dat,status=unknown)l lread(10,*)n_main,inform1,inform2read(10,*)n_main,inform1,inform2l lread(10,*)span,gg,eeread(10,*)span,gg,eel lnod=n_main-1 nod=n_main-1 结点数结点数l lndinf=2*n_main+1 =1ndinf=2*n_main+1 =1作用点数作用点数程序解读l lm=n_mainm=n_mainl lif(inform1.eq.0)m=1 if(inform1.eq.0)m=1 各主梁同各主梁同l ldo i=1,mdo i=1,ml l read(10,*)bi(i),ti(i),bmain(i),rmain(i),bic(i)read(10,*)bi(i),ti(i),bmain(i),rmain(i),bic(i)l lend doend dol lif(inform1.eq.0)thenif(inform1.eq.0)thenl l do i=2,n_main do i=2,n_mainl l bi(i)=bi(i-1)bi(i)=bi(i-1)l l ti(i)=ti(i-1)ti(i)=ti(i-1)l l bmain(i)=bmain(i-1)bmain(i)=bmain(i-1)l l bic(i)=bic(i-1)bic(i)=bic(i-1)l l rmain(i)=rmain(i-1)rmain(i)=rmain(i-1)程序解读l lend doend dol l end if end ifl l if(inform2.eq.0)then if(inform2.eq.0)thenl lread(10,*)(inform(i),i=1,nod)read(10,*)(inform(i),i=1,nod)刚、铰接混合刚、铰接混合l l end ifend ifl l if(inform2.eq.1)then if(inform2.eq.1)thenl ldo i=1,noddo i=1,nodl l inform(i)=1 inform(i)=1 铰接铰接l lend doend dol l end if end ifl l if(inform2.eq.2)then if(inform2.eq.2)thenl ldo i=1,noddo i=1,nod程序解读l l inform(i)=2 inform(i)=2l lend doend dol l end if end ifl l call inf call inf(n_main,nod,ndinf,span,ee,gg,bi,ti,bmain,inform),(n_main,nod,ndinf,span,ee,gg,bi,ti,bmain,inform),l l+ik,rmain,bic,xinf,gk,b,eta,tor,tm)+ik,rmain,bic,xinf,gk,b,eta,tor,tm)l l call output(11,n_main,nod,ndinf,xinf,ta,tor,tm)call output(11,n_main,nod,ndinf,xinf,ta,tor,tm)l l close(10)close(10)l l close(11)close(11)l l write(*,(/)write(*,(/)l l write(*,(a)Thats all for the program C62GJJBL!write(*,(a)Thats all for the program C62GJJBL!l l end end程序解读l lc cl lc cl lc Subroutines for C62GJJBLc Subroutines for C62GJJBLl lc cl lc cl lc cl lc Computer room/1995.05.17c Computer room/1995.05.17l lc cl lc cl lc -c -l lbubroutine cnamebubroutine cnamel lc -c -程序解读l lcharacter*12ccharacter*12cl ldo 10i=12,1,-1do 10i=12,1,-1l l if(c(i:i).eq.)then if(c(i:i).eq.)thenl l ip=i ip=il l goto 30 goto 30l l endif endifl l10 continue10 continuel l30 if(ip.gt.1)then30 if(ip.gt.1)thenl ldo 40i=ip,12do 40i=ip,12l l40 c(i:i)=40 c(i:i)=l lendifendifl lendend程序解读l lc -c -l lsubroutine fname(c,ce)subroutine fname(c,ce)l lc -c -l lcharacter c*12,ce*4character c*12,ce*4l lk=0k=0l ldo 10 i=1,8do 10 i=1,8l l if(c(i:i).eq