/PREP7
KEYW,PR_THERM,1
ET,1,PLANE55
KEYOPT,1,3,1
MP,KXX,1,0.0635
MP,DENS,1,2000
MP,C,1,443
*DIM,a,,30,144
*DIM,T,,30,144
*DIM,q,,30,144
RECTNG,0,0.15,0,0.02
MSHAPE,0,2D
MSHKEY,1
/REPLOT,RESIZE
/REPLOT,RESIZE
LESIZE,1,,,15
LESIZE,2,,,8
AMESH,1
FINISH
/SOLU
ANTYPE,TRANS
*SET,time,0
*SET,delt,0.5
time=15
TIME,time
*DO,j,1,144
T(1,j)=25
*ENDDO
*DO,j,1,144
a(1,j)=0.03
*ENDDO
LSEL,S,,,3
LSEL,A,,,1
NSLL,S,1
D,ALL,TEMP,120
*DO,i,1,30
*DO,j,1,144
NSEL,S,LOC,X,0,0.15
D,j,TEMP,T(i,j)
*ENDDO
*DO,j,1,144
q(i,j)=1.53E14*EXP(-8054.4/(T(i,j)+273))*a(i,j)**0.6*(1-a(i,j))**1.4
NSEL,S,LOC,X,0,0.15
BF,j,HGEN,q(i,j)
*ENDDO
SOLV
/POST1
i=i+1
*DO,j,1,144
*GET,T(i,j),NODE,j,TEMP
*ENDDO
*DO,j,1,144
a(i,j)=a(i-1,j)+delt*1.53E14*EXP(-8054.4/(T(i,j)+273))*a(i-1,j)**0.6*(1-a(i-1,j))**1.4
*ENDDO
*ENDDO |