!   TITLE:  hyfly_thermal.mac (HYFLY)
!
!-----------------------------------------------------------------------
!       Rev -   03/14/2002      J. Rinaldi
!-----------------------------------------------------------------------
!
!   This macro creates the actuator envelope for the thermal analysis
!   of the hyfly missile.
!
!   CShyfly0:       Cartesian Hyfly Coordinate System
!   CShyfly1:       Cylindrical Hyfly Coordinate System
!   CShyfly2:       Spherical Hyfly Coordinate System
!
!
!************************** TABLE OF CONTENTS **************************
!
!       Model Assembly Preparation
!
!   1.  Parameter List
!   2.  Model Coordinate System
!   3.  Solid Model
!   4.  Meshing
!
!********************* Model Assembly Preparation **********************
!
mod_assemb_prep,arg1,arg2,arg3
!
!<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>>
!
hyfly_thermal.parm                          ! name of parameter file
!
!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate System >>>>>>>>>>>>>>>>>>>>>
!
mod_coord_system,acs,xc,yc,zc,rotxyc,rotyzc,rotzxc
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 3.  Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!------------------------ 3.1  Actuator Section -----------------------
!
csys,CShyfly1
!
init_entity_num
!
kmi=km+1
!
k,,Dr03,90,Dz03
k,,Dr04,90,Dz04
k,,Dr05,90,Dz05
k,,Dr06,90,Dz06
k,,Dr07,90,Dz07
!
k,,Dr08,90,Dz08
k,,Dr09,90,Dz09
k,,Dr10,90,Dz10
k,,Dr11,90,Dz11
k,,Dr12,90,Dz12
!
k,,Dr13,90,Dz13
!
init_entity_num
!
csys,CShyfly0
kscale,0,kmi,km,1,scalgeom,scalgeom,scalgeom
csys,CShyfly1
!
ksel,s,kp,,km-10,km
!bsplin,all,,,,,,0,-1
spline,all,,,,,,0,-1
!
init_entity_num
!
kmi=km+1
!
k,,Dr13,90,Dz13
k,,Dr01,90,Dz01
k,,Dr02,90,Dz02
k,,Dr03,90,Dz03
!
csys,CShyfly0
kscale,0,kmi,kmi+3,1,scalgeom,scalgeom,scalgeom
csys,CShyfly1
!
create_lines,,1
nummrg,kp
!
create_area-lines
view,'f'
!
lplot
!
!-------------------------- 3.2  Skin Section -------------------------
!
csys,CShyfly1
!
kmi=km+1
!
k,,Dr15,90,Dz15
k,,Dr16,90,Dz16
k,,Dr17,90,Dz17
k,,Dr18,90,Dz18
k,,Dr19,90,Dz19
!
k,,Dr20,90,Dz20
k,,Dr21,90,Dz21
k,,Dr22,90,Dz22
k,,Dr23,90,Dz23
k,,Dr24,90,Dz24
!
k,,Dr25,90,Dz25
!
init_entity_num
!
csys,CShyfly0
kscale,0,kmi,km,1,scalgeom,scalgeom,scalgeom
csys,CShyfly1
!
ksel,s,kp,,km-10,km
!bsplin,all,,,,,,0,-1
spline,all,,,,,,0,-1
init_entity_num
!
kmi=km+1
!
k,,Dr25,90,Dz25
k,,Dr13+skn_thk,90,Dz13
k,,Dr01+skn_thk,90,Dz01
k,,Dr02-skn_thk,90,Dz02
k,,Dr03-skn_thk,90,Dz03
!
csys,CShyfly0
kscale,0,kmi,kmi+4,1,scalgeom,scalgeom,scalgeom
csys,CShyfly1
!
create_lines,,1
nummrg,kp
!
create_area-lines
aovlap,all
!
init_entity_num
!
lplot
!
!----------------------- 3.3  Insulation Section ----------------------
!
csys,CShyfly1
!
kmi=km+1
!
k,,Dr27,90,Dz27
k,,Dr28,90,Dz28
k,,Dr29,90,Dz29
k,,Dr30,90,Dz30
k,,Dr31,90,Dz31
!
k,,Dr32,90,Dz32
k,,Dr33,90,Dz33
k,,Dr34,90,Dz34
k,,Dr35,90,Dz35
k,,Dr36,90,Dz36
!
k,,Dr37,90,Dz37
!
init_entity_num
!
csys,CShyfly0
kscale,0,kmi,km,1,scalgeom,scalgeom,scalgeom
csys,CShyfly1
!
ksel,s,kp,,km-10,km
!bsplin,all,,,,,,0,-1
spline,all,,,,,,0,-1
init_entity_num
!
kmi=km+1
!
k,,Dr37,90,Dz37
k,,Dr13+skn_thk,90,Dz13
k,,Dr01+skn_thk,90,Dz01
k,,Dr27,90,Dz02
k,,Dr27,90,Dz03
!
csys,CShyfly0
kscale,0,kmi,kmi+4,1,scalgeom,scalgeom,scalgeom
csys,CShyfly1
!
create_lines,,1
!
isolate_lines
nummrg,kp
!
create_area-lines
aovlap,all
!
init_entity_num
!
lplot
!
!------------------------------ 4 Define Mesh -----------------------------
!
esize=.100*scalgeom
!
asel,s,area,,3,4,1
msh_area,,,1,esize,esize,,ti4900

! Inner area
asel,s,area,,1
!MSH_AREA, mstyp, mshdens, etpa, elmnl, elmxl, lszmsh, mat, real, ecs
msh_area,,,1,2*esize,2*esize,,sand
!
! Nozzle insulation
asel,s,area,,6
msh_area,,,1,esize,esize,,sand
!
eplot

!
!----------------------- 5  Loads -----------------------------------------
!

tunif,Tinit

! Apply load to inner skin
lsel,s,line,,32,41,1
lsel,a,line,,45
nsll,s,1
sf,all,conv,h_air,Tnoz
!d,all,temp,Tnoz

!dl,all,,temp,Tnoz
!ic,all,temp,Tnoz

! Apply load to outer skin
lsel,s,line,,25
nsll,s,1
sf,all,conv,h_air,Tskin
!d,all,temp,Tskin

!dl,25,,temp,Tskin
!ic,all,temp,Tskin

!
!----------------------- 6  Solve -----------------------------------------
!
alls
!
/solu
!
antype,4
trnopt,full
lumpm,0
!
*if,scalgeom,eq,1,then
    ! Units of time: hours
    time,12/60
    autots,-1
    deltim,.005,.0001,.01,1
*else
    ! Units of time: seconds
    time,12*60
    autots,-1
    deltim,18,.36,36,1
*endif
!
kbc,0
!
!tsres,erase
!
outres,all,all
!
/status,solu
solve
!
/post1
set,first
!
eplot
!
/cval
!/cval,1,2643,2750,3000,0,0,0,0,0
/replot
!
plns,temp,
antime,30,0.5,,0,0,0,0

!
! Select elements in actuator envelope
!
csys,0
nsel,s,loc,y,-1.05,10.44
esln,s,1
asel,s,area,,1
esla,r


:skip