帮忙看个辐射的问题
要有个初期的初始温度，然后考虑通过辐射来冷却。在把初始温度加载进去以后，好像是个加热的过程。请各位帮忙看看。/filnam,initial-temperature
r1=0.6e-3
r2=0.8e-3
r3=0.005e-3
r4=0.02e-3
y1=10e-3
y2=12e-3
/prep7
et,1,plane223,,,1
et,2,link32
mptemp,1,100,200,300,400,450,500,
mptemp,7,550,600,800,900
mpdata,ex,1,1,6.762e10,6.2e10,5.8e10,5.1e10,4.7e10,4.0e10
mpdata,ex,1,7,8.2e9,6.9e9,4.8e9,3.2e9
mpdata,prxy,1,1,0.23,0.23,0.23,0.23,0.23,0.23
mpdata,prxy,1,7,0.23,0.23,0.23,0.23
mpdata,kxx,1,1,0.01125e2,0.01175e2,0.0125e2,0.01325e2,0.01375e2,0.0143e2
mpdata,kxx,1,7,0.0148e2,0.0148e2,0.0137e2,0.01285e2
mpdata,c,1,1,897,1005,1070,1123,1147,1223
mpdata,c,1,7,1346,1391,1397,1403
mpdata,dens,1,1,2436,2425,2415,2405,2401,2398
mpdata,dens,1,7,2387,2375,2353,2327
mpdata,alpx,1,1,85.8e-7,88.7e-7,95.1e-7,102.2e-7,108.4e-7,111.8e-7
mpdata,alpx,1,7,198.6e-7,201.2e-7,205.1e-7,208.5e-7
k,1,r1
k,2,r2
k,3,r1,y1
k,4,r2,y1
k,5,r3,y2
k,6,r4,y2
a,1,2,4,3
a,3,4,6,5
aglue,all
lsel,s,,,2
lesize,all,,,40
lsel,s,,,1
lesize,1,,,4
mat,1
mshkey,1
amesh,all
!
nsel,s,loc,y,11.8e-3,12e-3
d,all,temp,960
nsel,s,loc,y,11.6e-3,11.8e-3
d,all,temp,820
nsel,s,loc,y,11.4e-3,11.6e-3
d,all,temp,700
nsel,s,loc,y,11.0e-3,11.4e-3
d,all,temp,580
nsel,s,loc,y,10.3e-3,11.0e-3
d,all,temp,475
nsel,s,loc,y,10e-3,10.3e-3
d,all,temp,350
allsel
!
/solu
antype,trans
time,0.2
tref,25
timint,0,therm
solve
!/post1$plnsol,temp
/solu
nsel,s,loc,y,10e-3,12e-3
ddele,all,temp
!lsel,s,,,1,3,2$sfl,all,conv,3.5,,25$lsel,a,,,6,7,1$sfl,all,conv,3.5,,25
allsel
antype,4
time,1
tref,450
timint,1,therm
solve
finish
save
/filnam,frozening
/prep7
type,2
lsel,s,ext
lplot
nsll,s,1
nplot
n,1000,3*r1,2*y1
esurf
/PSYMB,ESYS,1 ! turn on esys plotting
esel,s,type,,2
eplot
allsel
finish
/AUX12 ! enter radiation matrix utility
EMIS,1,0.7, ! specify emissivity
STEF,5.67e-8, ! specify Stefan-Boltzman constant
GEOM,1,0, ! specify geometry type
SPACE,1000, ! specify space node
VTYPE,0, , ! specify HIDDEN method
MPRINT,0 ! do not print matrix
WRITE,radtest ! write radation matrix file, radtest1
FINISH ! exit the radiation matrix utility
!
/prep7
esel,s,type,,2 ! select all type 3 elements
/PSYMB,ESYS,0 ! turn off esys plotting
eplot
EDELE,ALL ! delete all type 3 elements
ALLSEL
et,3,matrix50
keyopt,3,1,1
type,3
se,radtest,,,0.0001
finish
!
/solu
antype,0
!LDREAD,TEMP,2,1, ,1,'initial-temperature','rst',' '
LDREAD,TEMP,,,0.8,1,'initial-temperature','rst',' '
!nsel,s,loc,y,0$d,all,uy
nlgeom,0
NROPT,AUTO, ,on !定义H法选项
LUMPM,0
EQSLV,FRONT, ,0, !求解器选择
toffst,273
d,1000,temp,200
time,10
nsubst,10$outres,all,1
solve
页:
[1]