线弹性小变形平面壳静力有限元计算程序.docVIP

元计算有限元自动生成系统所开发源代码系列 线弹性小变形平面板壳静力有限元计算程序 简介 元计算（www.ectec.asia）公司所开发的并行有限元程序自动生成系统（pFEPG）可根据用户需要开发出各种有限元计算程序源代码。该源代码系列即为pFEPG所开发出来的求解各学科典型问题的有限元计算程序。 该组程序为线弹性小变形平面板壳静力有限元计算程序。 starta.for，对位移场的数据进行初始化； implicit real*8 a-h,o-z character*12 fname,filename 20 common /aa/ ia 座机电话号码0 common /bb/ ib 1座机电话号码 c.... open disp0 file to get the numbers of nodes and degree of freedom c.... knode .... number of nodes, kdgof .... number of d.o.f. open 1,file ,form unformatted read 1 knode,kdgof close 1 kvar knode*kdgof write *,* knode,kdgof,kvar write *, 1x,4i7 knode,kdgof,kvar kvar1 kvar+1 kcoor 3 kelem 座机电话号码 knb1 kdgof*knode*1 if knb1/2*2 .lt. knb1 knb1 knb1+1 kna4 kcoor*knode*2 kna1 kdgof*knode*2 kna2 kdgof*knode*2 kna3 kdgof*knode*2 kna5 knode*1 if kna5/2*2 .lt. kna5 kna5 kna5+1 knb4 kelem*1 if knb4/2*2 .lt. knb4 knb4 knb4+1 knb2 kvar1*1 if knb2/2*2 .lt. knb2 knb2 knb2+1 knb3 kvar1*1 if knb3/2*2 .lt. knb3 knb3 knb3+1 kna0 1 kna1 kna1+kna0 kna2 kna2+kna1 kna3 kna3+kna2 kna4 kna4+kna3 kna5 kna5+kna4 if kna5-1.gt.座机电话号码0 then write *,* exceed memory of array ia write *,* memory of ia 座机电话号码0 write *,* memory needed ,kna5, in prgram start stop 55555 endif knb0 1 knb1 knb1+knb0 knb2 knb2+knb1 knb3 knb3+knb2 knb4 knb4+knb3 if knb4-1.gt.1座机电话号码 then write *,* exceed memory of array ib write *,* memory of ib 1座机电话号码 write *,* memory needed ,knb4, in prgram start stop 55555 endif call start knode,kdgof,kcoor,kvar, *kelem,maxt,kvar1,ia kna0 ,ia kna1 ,ia kna2 , *ia kna3 ,ia kna4 ,ib knb0 ,ib knb1 ,ib knb2 , *ib knb3 , *filename end subroutine start knode,kdgof,kcoor,kvar, *kelem,maxt,kvar1,u0,u1,u2, *coor,inodvar,nodvar,numcol,lm,node, *filename implicit real*8 a-h,o-z character*12 filename 20 DIMENSION NODVAR KDGOF,KNODE ,COOR KCOOR,KNODE ,R 3 , * U0 KDGOF,KNODE ,U1 KDGOF,KNODE ,U2 KDGOF,KNODE , * INODVAR KNODE ,node kelem DIMENSION NUMCOL KVAR1 ,LM KVAR1 CHARACTER*1 MATERIAL logical filflg C .....................................................