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

整理ppt程序解读inform(i)=2enddoendifcallinf(n_main,nod,ndinf,span,ee,gg,bi,ti,bmain,inform),+ik,rmain,bic,xinf,gk,b,eta,tor,tm)calloutput(11,n_main,nod,ndinf,xinf,ta,tor,tm)close(10)close(11)write(*,′(//)′)write(*,′(a\)′)′That’sallfortheprogramC62GJJBL!′end程序解读cccSubroutinesforC62GJJBLccccComputerroom/1995.05.17ccc----------------------------------------bubroutinecname?c----------------------------------------程序解读character*12cdo10i=12,1,-1if(c(i:i).eq.′.′)thenip=igoto30endif10continue30if(ip.gt.1)thendo40i=ip,1240c(i:i)=′′endifend程序解读c-------------------------------------subroutinefname(c,ce)c-------------------------------------characterc*12,ce*4k=0do10i=1,8if(c(i:i).eq.′′)thenk=igoto15endif10continue15if(k.eq.0)k=9程序解读do20i=k,1220c(i:i)=′′do25i=k,k+325c(i:i)=ce(j:j)endc----------------------------------------------subroutinedecom(neq0,neq,nband,gk)c-----------------------------------------------!Decomposing[K]=LDL′byCHOLESCKYinpovedsquarerootmethod.Dimensiongk(neq0,4),t(4)doi=1,neq程序解读1=i-nband+1if(l.lt.1)l=1doj=1,iji=nband-i+jqq=0.0if(j.ne.1)thenj1=j-1dok=1,j1ki=nband+k-ikj=nband+k-jqq=qq+t(ki)*gk(j,kj)enddo程序解读endifif(i.eq.j)thengk(i,nband)=gk(i,nband)-qqelseif(i.gt.j)thent(ji)=gk(i,ji)-qqgk(i,ji)=t(ji)/gk(j,nband)endifenddoenddoenddoendc----------------------------------------程序解读subroutineback(neq0,neq,nband,gk,b)c------------------------------------------!Returnsubstitutionandforwardsubstitutiondimensiongk(neq0,4),b(neq)doi=2,neqqq=0.0l=i-nband+1if(l.lt.1)1=1il=i-1doj=1,ilji=nband-i+j