Fotran90版—平面刚架有限元分源程序代码.doc

Fotran90版—平面刚架有限元分源程序代码 program main real,allocatable::ks:,: allocatable lnd:,: allocatable crd:,: allocatable ea: allocatable ei: allocatable jcs:,: allocatable pj:,: allocatable bl: allocatable p: open 5,fileinputdates.in read 5,* ne,nj,ns,npj nj33*nj allocate lndne,2 allocate crdnj,2 allocate eane allocate eine allocate jcsns,4 allocate pjnpj,4 allocate blne allocate ksnj3,nj3 allocate pnj3 write *,1x,plane fram structural analysis// write *,1x,structural parameters/ write *,/1x,total number of write *,1x,element,i5/1x,joints,i5/1x,constructed joints,i5/1x,loads,i5/ne,nj,ns,npj call readinne,nj,ns,npj,lnd,crd,ea,ei,jcs,pj,blcall formfnj3,npj,pj,p call cksnj3,ne,ea,ei,bl,lnd,crd,nj,ks call dealbcns,jcs,nj3,ks,p call solvenj,nj3,ks,p close 5 stop end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! subroutine readinne,nj,ns,npj,lnd,crd,ea,ei,jcs,pj,bl dimension lndne,2,crdnj,2,eane,eine,jcsns,4,pjnpj,4,blne read 5,* lndi,j,j1,2,i1,ne write *,/1x,element dates,/1x,element,4x,conection,8x write *,1x,i5,2x,i5,3x,to,i5,3x i,lndi,j,j1,2,i1,ne read 5,* crdi,j,j1,2,i1,nj write *,1x,nodel coordinates/3x,node,6x,x-coordinates,7x,y-coordinates write *,1x,i5,5x,f10.4,10x,f10.4 i,crdi,j,j1,2,i1,nj read 5,* eai,i1,ne read 5,* eii,i1,ne write *,/1x,material parameters,/1x,element,11x,ea,13x,ei write *,1x,i5,5x,2e15.6 i,eai,eii,i1,ne read 5,* jcsi,j,j1,4,i1,ns write *,/1x,constrained nodes,/3x,nodes,1x,X,4x,Y,4x,R write *,4i5 jcsi,j,j1,4,i1,ns read 5,* pji,j,j1,4,i1,npj write *,/1x,joints of londs,/1x,joint,5x,PX,8x,PY,8x,Mxy write *,1x,f5.0,3f10.4 pji,j,j1,4,i1,npj do ie1,neilndie,1jlndie,2dxcrdj,1-crdi,1dycrdj,2-crdi,2bliesqrtdx**2+dy**2 end do write *,/1x,the length of the elements,/1x,element numers,5x,numembers length write *,1x,i5,10x,f10.4 i,bli,i1,ne return end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! subroutine turnne,lnd,crd,bl,t,ie,nj dimension lndne,2,crdnj,2,blne,t6,6 i1lndie,1 j1lndie,2 dxcrdj1,1-crdi1,1 dycrdj1,2-crdi1,2 sidy/blie codx/