/COM, Thermal
/COM, Magnetic-Nodal
!*
/PREP7
!*
ET,1,PLANE13
!*
KEYOPT,1,1,6
KEYOPT,1,2,0
KEYOPT,1,3,0
KEYOPT,1,4,0
KEYOPT,1,5,0
!*
ET,2,PLANE13
!*
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,0
KEYOPT,2,5,0
!*
ET,3,PLANE13
!*
KEYOPT,3,1,6
KEYOPT,3,2,0
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,0
!*
ET,4,SURF151
!*
KEYOPT,4,1,0
KEYOPT,4,2,0
KEYOPT,4,3,0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,6,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
!*
!*
!*
R,4, , , , , ,
RMORE, , , , , ,
RMORE, , ,
!*
!*
!*
EMUNIT,MKS
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,1,,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,1,,1.678e-8
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,2,,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,25.5
MPTEMP,3,160
MPTEMP,4,291.5
MPTEMP,5,477.6
MPTEMP,6,635
MPTEMP,7,698
MPTEMP,8,709
MPTEMP,9,720.3
MPTEMP,10,742
MPTEMP,11,761
MPTEMP,12,1000
MPDATA,MURX,3,,205
MPDATA,MURX,3,,200
MPDATA,MURX,3,,190
MPDATA,MURX,3,,182
MPDATA,MURX,3,,161
MPDATA,MURX,3,,135
MPDATA,MURX,3,,104
MPDATA,MURX,3,,84
MPDATA,MURX,3,,35
MPDATA,MURX,3,,17
MPDATA,MURX,3,,1
MPDATA,MURX,3,,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,125
MPTEMP,3,250
MPTEMP,4,375
MPTEMP,5,500
MPTEMP,6,625
MPTEMP,7,750
MPTEMP,8,875
MPTEMP,9,1000
MPDATA,RSVX,3,,0.184e-6
MPDATA,RSVX,3,,0.272e-6
MPDATA,RSVX,3,,0.384e-6
MPDATA,RSVX,3,,0.512e-6
MPDATA,RSVX,3,,0.656e-6
MPDATA,RSVX,3,,0.824e-6
MPDATA,RSVX,3,,1.032e-6
MPDATA,RSVX,3,,1.152e-6
MPDATA,RSVX,3,,1.2e-6
SAVE
PCIRC,0.003, ,0,360, !cable
PCIRC,0.012, ,0,360, !inner radius of guan
PCIRC,0.018, ,0,360, !outer radius of guan
PCIRC,0.060, ,0,360, !outer radius of model
FLST,2,4,5,ORDE,2
FITEM,2,1
FITEM,2,-4
AOVLAP,P51X
SAVE
!*
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 1, 4, 1, 0, !copper cable
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,2,5,ORDE,2
FITEM,5,5
FITEM,5,7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 2, 4, 2, 0, !air
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,AREA
ASEL, , , , 6
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 3, 4, 3, 0, !guan
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
SMRT,6 !mesh all
SMRT,2
MSHAPE,0,2D
MSHKEY,0
!*
FLST,5,4,5,ORDE,3
FITEM,5,1
FITEM,5,5
FITEM,5,-7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,2,4,5,ORDE,3
FITEM,2,1
FITEM,2,5
FITEM,2,-7
ACLEAR,P51X
APLOT
MSHAPE,1,2D
!*
FLST,5,4,5,ORDE,3
FITEM,5,1
FITEM,5,5
FITEM,5,-7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y '
CMDELE,_Y1
CMDELE,_Y2
!*
APLOT
EPLOT
SAVE
n ! create space node for SURF151
*get,nmax,node,,num,max
FLST,5,4,4,ORDE,2
FITEM,5,5
FITEM,5,-8
LSEL,S, , ,P51X
NSLL,S,1
/UI,MESH,OFF
ALLSEL,ALL
FLST,5,4,4,ORDE,2
FITEM,5,5
FITEM,5,-8
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y
!*
!*
CMSEL,S,_Y1
LATT,3,4,4, , , ,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
ALLSEL,ALL
FLST,2,4,4,ORDE,2
FITEM,2,5
FITEM,2,-8
LMESH,P51X
SAVE
*get,emax,elem,,num,max
emodif,emax,4,nmax ! modify element to add space node for radiation
EPLOT
et,4,0 ! reset type 4 to null element
/UI,MESH,OFF
APLOT
ALLSEL,ALL
ASEL,S, , , 1
NSLA,S,1
FLST,4,49,1,ORDE,6
FITEM,4,1
FITEM,4,-4
FITEM,4,17
FITEM,4,-48
FITEM,4,185
FITEM,4,-197
CP,1,VOLT,P51X !couple volt in cupper cable
FLST,2,49,1,ORDE,6
FITEM,2,1
FITEM,2,-4
FITEM,2,17
FITEM,2,-48
FITEM,2,185
FITEM,2,-197
/GO
!*
F,P51X,AMPS,160 !apply current in copper cable
ALLSEL,ALL
ASEL,S, , , 6
NSLA,S,1
FLST,4,148,1,ORDE,8
FITEM,4,5
FITEM,4,-12
FITEM,4,49
FITEM,4,-116
FITEM,4,481
FITEM,4,-516
FITEM,4,1124
FITEM,4,-1159
CP,2,VOLT,P51X !couple volt in guan
FLST,2,148,1,ORDE,8
FITEM,2,5
FITEM,2,-12
FITEM,2,49
FITEM,2,-116
FITEM,2,481
FITEM,2,-516
FITEM,2,1124
FITEM,2,-1159
/GO
!*
F,P51X,AMPS,-160 !aplly current in guan
ALLSEL,ALL
FINISH
/SOL
!*
ANTYPE,3
HARFRQ,50,0,
NSUBST,0,
KBC,0
!*
FINISH
/PREP7
PHYSICS,WRITE,emage, , ,
FINISH
/prep7
lsclear,all ! clear all b.c.'s and options
et,1,0
et,2,0
et,3,55,,,1
et,4,151,,,1,1,1 ! SURF151 element for radiation
keyopt,4,9,1
r,4,1,5.67e-8 ! form factor, Stefan-Boltzmann constant
SAVE
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,730
MPTEMP,3,930
MPTEMP,4,1000
MPDATA,KXX,3,,60.64
MPDATA,KXX,3,,29.5
MPDATA,KXX,3,,28
MPDATA,KXX,3,,28.
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,27
MPTEMP,3,127
MPTEMP,4,327
MPTEMP,5,527
MPTEMP,6,727
MPTEMP,7,765
MPTEMP,8,765.001
MPTEMP,9,927
MPDATA,ENTH,3,,0
MPDATA,ENTH,3,,91609056
MPDATA,ENTH,3,,453285756
MPDATA,ENTH,3,,1.2748e9
MPDATA,ENTH,3,,2.2519e9
MPDATA,ENTH,3,,3.3396e9
MPDATA,ENTH,3,,3.548547e9
MPDATA,ENTH,3,,3.548556e9
MPDATA,ENTH,3,,4.3520e9
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EMIS,3,,0.68
SAVE
FINISH
/PREP7
FINISH
/SOLU
!*
antype,trans
!*
TOFFST,273,
!*
TUNIF,18,
cnvtol,heat,1 ! convergence tolerance
kbc,1 ! step loads
trnopt,full
autos,on ! auto time-stepping
*SET,ftime,12 ! final time
*SET,tinc,.05 ! time increment for harmonic analysis
*SET,time,0 ! initialize time
*SET,delt,.01 ! maximum delta time step
deltim,1e-5,1e-6,delt,on ! time step control
outres,basic,all ! save all load step information
physics,write,thermal ! write thermal physics file
finish
*do,i,1,ftime/tinc ! solution *do loop
*SET,time,time+tinc ! increment time
SAVE
physics,read,emag ! read emag physics file
/solu
*if,i,eq,1,then
tunif,18
*else
ldread,temp,last,,,,,rth ! read thermal analysis temperatures
*endif
solve ! solve harmonic analysis
finish
physics,read,thermal ! read thermal physics file
/assign,esav,therm,esav ! redirect files for use in thermal restart
/assign,emat,therm,emat
/solu
*if,i,gt,1,then
antype,trans,rest ! thermal restart
*endif
time,time ! time at end of thermal run
esel,s,mat,,3 ! select billet region
ldread,hgen,,,,2,,rmg ! apply coupled joule heating load from emag
esel,all
solve
finish
/assign,esav ! reassign files to default
/assign,emat
*enddo ! end of solution looping
finish
save ! save database
/post26 ! time-history postprocessor
/show
nsol,2,1132,temp,,tempcl ! store temperature at billet centerline
nsol,3,85,temp,,tempsurf ! store temperature at billet outer diameter
plvar,2,3 ! plot temperature rise over time
prvar,2,3 ! print temperature rise over time
finish |