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

Fotran90版—平面刚架有限元分源程序代码 program main real,allocatable::ks(:,:) allocatable lnd(:,:) allocatable crd(:,:) allocatable ea(:) allocatable ei(:) allocatable jcs(:,:) allocatable pj(:,:) allocatable bl(:) allocatable p(:) open (5,file=inputdates.in) read (5,*) ne,nj,ns,npj nj3=3*nj allocate (lnd(ne,2)) allocate (crd(nj,2)) allocate (ea(ne)) allocate (ei(ne)) allocate (jcs(ns,4)) allocate (pj(npj,4)) allocate (bl(ne)) allocate (ks(nj3,nj3)) allocate (p(nj3)) 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 readin(ne,nj,ns,npj,lnd,crd,ea,ei,jcs,pj,bl) call formf(nj3,npj,pj,p) call cks(nj3,ne,ea,ei,bl,lnd,crd,nj,ks) call dealbc(ns,jcs,nj3,ks,p) call solve(nj,nj3,ks,p) close (5) stop end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! subroutine readin(ne,nj,ns,npj,lnd,crd,ea,ei,jcs,pj,bl) dimension lnd(ne,2),crd(nj,2),ea(ne),ei(ne),jcs(ns,4),pj(npj,4),bl(ne) read (5,*) ((lnd(i,j),j=1,2),i=1,ne) write (*,(/1x,element dates,/1x,element,4x,conection,8x)) write (*,(1x,i5,2x,i5,3x,to,i5,3x)) (i,(lnd(i,j),j=1,2),i=1,ne) read (5,*) ((crd(i,j),j=1,2),i=1,nj) write (*,(1x,nodel coordinates/3x,node,6x,x-coordinates,7x,y-coordinates)) write (*,(1x,i5,5x,f10.4,10x,f10.4)) (i,(crd(i,j),j=1,2),i=1,nj) read (5,*) (ea(i),i=1,ne) read (5,*) (ei(i),i=1,ne) write (*,(/1x,material parameters,/1x,element,11x,ea,13x,ei)) write (*,(1x,i5,5x,2e15.6)) (i,ea(i),ei(i),i=1,ne) read (5,*) ((jcs(i,j),j=1,4),i=1,ns) write (*,(/1x,constrained nodes,/3x,nodes,1x,X,4x,Y,4x,R)) write (*,(4i5)) ((jcs(i,j),j=1,4),i=1,ns) read (5,*) ((pj(i,j),j=1,4),i=1,npj) write (*,(/1x,joints of londs,/1x,joint,5x,PX,8x,PY,8x,Mxy)) write (*,(1x,f5.0,3f10.4)) ((pj(i,j),j=1,4),i=1,npj) do ie=1,ne i=lnd(ie,1) j=lnd(ie,2) dx=crd(j,1)-crd(i,1) dy=crd(j,2)-crd(i,2) bl(ie)=sqrt(dx**2+dy**2) end do write (*,(/1x,the length of the e