一维非稳态导热通用程序.docVIP

!*************************************************************** ! 一维非稳态导热通用程序(不变部分) ! This is a general purpose program to solve 1-D diffusion ! problem in the form of: ! ρcdt/dz=1/a(x)d/dx(a(x)Γdt/dx)+s !******************2003.7 revised******************************** !...................Define Variables......................... MODULE VARIABLES INTEGER,PARAMETER::L1=130 REAL,DIMENSION(L1):: X,XF,XM,XP,R,RF,AP REAL,DIMENSION(L1)::AE,AW,CN,T,TA REAL,DIMENSION(L1)::TG,GM,RC INTEGER:: K=1,KM=1,KP=1,OM=1 INTEGER:: JB,JE,KE,KI,KF,KN,KR,KT,LS,MD,M1,M2,NF REAL:: AEC,AI,BE,BI,DF,DS,DT,EP,EX,PW,TU,TM,XE,XI REAL:: A1,A2,T2,TC,SC,SP,RO,TE,DN,LM END !...............................................Main Program........................................... PROGRAM MAIN USE VARIABLES IMPLICIT NONE INTEGER I OPEN(1,FILE=q.dat) OPEN(2,FILE=temp.dat) NF=1！（求解变量指标） KN=1！（非稳态问题输出局部变量，输出一次加1） TU=0！（当前时间） 150 KT=1！（非线性问题迭代次数） CALL Speci !First to specify the problem CALL Grid !Set up grid points 200 CALL Difsor !Specify the diff-coeff and source term 220 CALL InterOutput !Output intermediate results CALL Coeff !Set up coefficients of discretization equation CALL TDMA !Solve the algebraic equation by TDMA IF(LS.EQ.2.OR.LS.EQ.4) THEN ！（对非线性问题） IF(DF.GT.EP) THEN！（如果最大偏差大于允许值） DO I=1,M1 TA(I)=TA(I)+OM*(T(I)-TA(I))！（采用亚松弛方式将当前解付给上一次迭代值） END DO DF=0！（最大偏差置零） KT=KT+1！（非线性问题迭代次数加1） GOTO 200！（转到DIFSOR模块重新计算扩散系数与源项） END IF END IF CALL GPRINT !(四类问题均要经过一般输出) IF(LS.EQ.3.OR.LS.EQ.4) THEN！（对非稳态问题） IF(TU.LT.TM) THEN！（时间小于设定的最大值） DO I=1,M1 TG(I)=T(I)！（当前计算结果付给上一时层） END DO KT=1 IF(LS.EQ.3) THEN！（非稳态线性问题） GOTO 220！（转到中间输出模块） ELSE（非稳态非线性问题） GOTO 200！（转到DIFSOR模块重新计算扩散系数与源项） END IF END IF END IF !.............special results print out, if not, just leave it open.......... CALL SPRINT