[Fortran] 纯文本查看 复制代码
!Master plane
program Bok1(input,output,tape10=input,tape12=output,tape2,tape4,tape6,tape8)
dimension title(12)
common/Contro/npoin,nelem,nnode,ndofn,ndime,nstre,ntype,ngaus,nprop,nmats,nvfix,nevab,&
icase,ncase,itemp,iprob,nprob
common/Lgdata/coord(80,2),props(10,5),presc(40,2),asdis(160),eload(25,16),strin(4,225),&
nofix(40),ifpre(40,2),lnods(25,8),matno(25)
common/Work/elcod(2,8),shape(8),deriv(2,8),dmatx(3,3),cartd(2,8),dbmat(3,16),bmatx(3,16),&
smatx(3,16,9),posgp(3),weigp(3),gpcod(2,9),neror(24)
!**********************************************************
! program for the solution of plane
! frameworks
!**********************************************************
!
! this segment controls the calling,in order, of all subroutines
!
!call Openfile
read(10,900) nprob
900 format(i5)
write(12,905) nprob
905 format(1ho,5x,23htotal no. of problems =,i5)
do 20 iprob=1,nprob
rewind 2
rewind 4
rewind 6
rewind 8
read(10,910) title
910 format(12a6)
write(12,915) iprob,title
915 format(/////,6x,12hproblem no. ,i3,10x,12a6)
!
! call the subroutine which reads most of
! the problem data
!
call Input
!
! next create the element stiffness file
!
call Stifps
do 10 icase=1,ncase
!
! compute loads, after reading the relevant
! extra data
!
call Loadps
!
! merge and solve the resulting equations
! by the prontal solver
!
call Front
!
! compute the stresses in all the elements
!
call Streps
10 continue
20 continue
stop
end