介绍两个例子:
1 点对面
!Radiation modeling using surface effect element method for 2d model.
finish
/clear
/PREP7
ET,1,PLANE55
ET,2,SURF151
KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,9,1
R,1,.6,5.67e-8, , ,
RMORE, , , , , ,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,50
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,500
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,1500
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EMIS,1,,0.85
k,1,-0.75,0.15,0
k,2,-0.25,0.25,0
k,3,0.3,0.1,0
k,4,0.8,0.25,0
k,5,1.45,0.2,0
k,6,-0.75,-0.55,0
k,7,1.45,-0.55,0
k,8,0.5,0.5,0
bsplin,1,2,3,4,5
a,1,5,7,6
lsel,s,,,1
lesi,all,,,25
lsel,s,,,2
lesi,all,,,8
TYPE, 1
MAT, 1
REAL, 1
MSHAPE,0,2D
MSHKEY,1
AMESH,all
NKPT,0, 8
*get,nmax,node,,num,max
lsel,s,,,1
nsll,s,1
nsel,a,,,nmax
TYPE, 2
MAT, 1
REAL, 1
TSHAP,LINE
ESURF,nmax
ALLSEL,ALL
D,nmax,TEMP,500
lsel,s,,,3
nsll,s,1
d,all,temp,300
alls
FINISH
/SOLU
ANTYPE,0
tunif,300
NROPT,FULL, ,
OUTRES,ALL,ALL,
TIME,1
AUTOTS,1
NSUBST,10,1000,5,1
KBC,1
SOLVE
FINISH
2 面对面
finish
/clear
/filname,surf_radi1,1
/PREP7
ET,1,SOLID70
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,190
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,2700
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,900
rd1=1.0
rd2=0.9
rd3=0.3
rd4=0.25
z0=0
z1=0.2
dt=0.0001
CYLIND,rd1,rd2,z0,z1,0,180,
wpof,-0.30000
CYLIND,rd3,rd4,z0,z1,0,180,
VSYMM,Y,all, , , ,0,0
NUMMRG,KP, , , ,LOW
CSWPLA,11,1,1,1
lsel,s,loc,z,z0+dt,z1-dt
lesi,all,,,6,1
lsel,s,loc,x,rd3
lsel,r,loc,z,z0
lesi,all,,,20,1
lsel,s,loc,y,0
nsel,r,loc,z,z0
lesi,all,,,4
csys,1
lsel,s,loc,x,rd1
lsel,r,loc,z,z0
lesi,all,,,25
alls
MSHAPE,0,3D
MSHKEY,1
VMESH,all
asel,s,loc,x,rd2
sfa,all,,rdsf,0.85,1
csys,11
asel,s,loc,x,rd3
sfa,all,,rdsf,0.85,1
alls
nsel,s,loc,x,rd4
d,all,temp,1000
csys,1
nsel,s,loc,x,rd1
d,all,temp,300
alls
finish
/solu
tunif,300
STEF,0.567E-07
TOFFST,0
RADOPT,0.1,0.1E-03,0,1000,0.1,0.1
SPCTEMP,1,300
HEMIOPT,10
V2DOPT,0,20,1,200
VFOPT,OFF, , , ,BINA
alls
time,1
kbc,1
nsubst,10,100,5
save
solve |