诚信·公平·开放·共赢LoyaltyFairOpeningWin-win元计算有限元自动生成系统所开发源代码系列各向同性线性二维稳态传热有限元计算程序1.简介元计算所开发的并行有限元程序自动生成系统（pFEPG）可根据用户需要开发出各种有限元计算程序源代码。该源代码系列即为pFEPG所开发出来的求解各学科典型问题的有限元计算程序。该组程序为各向同性线性二维稳态传热有限元计算程序。2.starta.for，对温度场的数据进行初始化；implicitreal*8(a-h,o-z)character*12fname,filename(20)common/aa/ia(250000000)common/bb/ib(125000000)c....opendisp0filetogetthenumbersofnodesanddegreeoffreedomc....knode....numberofnodes,kdgof....numberofd.o.f.open(1,file='',form='unformatted')read(1)knode,kdgofclose(1)kvar=knode*kdgofwrite(*,*)'knode,kdgof,kvar='write(*,'(1x,4i7)')knode,kdgof,kvarkvar1=kvar+1kcoor=3kelem=31250000knb1=kdgof*knode*1if(knb1/2*2.lt.knb1)knb1=knb1+1kna4=kcoor*knode*2kna1=kdgof*knode*2kna2=kdgof*knode*2kna3=kdgof*knode*2kna5=knode*1if(kna5/2*2.lt.kna5)kna5=kna5+1knb4=kelem*1if(knb4/2*2.lt.knb4)knb4=knb4+1knb2=kvar1*1if(knb2/2*2.lt.knb2)knb2=knb2+1knb3=kvar1*1if(knb3/2*2.lt.knb3)knb3=knb3+1kna0=1kna1=kna1+kna0kna2=kna2+kna1kna3=kna3+kna2地址：北京海淀上地三街9号嘉华大厦A座811室邮编：100085中新天津生态城汉北路7号邮编：300467电话：（8610）82345565传真：（8610）82345545网址：www.ectec.asia邮箱：ectec＠ectec.asia诚信·公平·开放·共赢LoyaltyFairOpeningWin-winkna4=kna4+kna3kna5=kna5+kna4if(kna5-1.gt.250000000)thenwrite(*,*)'exceedmemoryofarrayia'write(*,*)'memoryofia=250000000'write(*,*)'memoryneeded=',kna5,'inprgramstart'stop55555endifknb0=1knb1=knb1+knb0knb2=knb2+knb1knb3=knb3+knb2knb4=knb4+knb3if(knb4-1.gt.125000000)thenwrite(*,*)'exceedmemoryofarrayib'write(*,*)'memoryofib=125000000'write(*,*)'memoryneeded=',knb4,'inprgramstart'stop55555endifcallstart(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)endsubroutinestart(knode,kdgof,kcoor,kvar,*kelem,maxt,kvar1,u0,u1,u2,*coor,inodvar,nodvar,numcol,lm,node,*filename)implicitreal*8(a-h,o-z)character*12filename(20)DIMENSIONNODVAR(KDGOF,KNODE),COOR(KCOOR,KNODE),R(3),*U0(KDGOF,KNODE),U1(KDGOF,KNODE),U2(KDGOF,KNODE),*INODVAR(KN