!   TITLE:  nozzle.mac (COMMON_MISSILE)
!
!-----------------------------------------------------------------------
!   Rev -   08/06/2002      J. Rinaldi
!-----------------------------------------------------------------------
!
!   This macro creates the variable exhaust nozzle for the common missile
!   project.
!
!   CSnoz0:     Cartesian Nozzle Coordinate System
!   CSnoz1:     Cylindrical Nozzle Coordinate System
!   CSnoz2:     Spherical Nozzle Coordinate System
!
!
!************************** TABLE OF CONTENTS **************************
!
!       Model Assembly Preparation
!
!   1.  Parameter List
!   2.  Model Coordinate System
!   3.  Solid Model
!
!********************* Model Assembly Preparation **********************
!
mod_assemb_prep,arg1,arg2,arg3
!
!<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>>
!
nozzle.parm
!
!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate System >>>>>>>>>>>>>>>>>>>>>
!
mod_coord_system,acs,xc,yc,zc,rotxyc,rotyzc,rotzxc
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 3.  Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!---------------------- 3.1  Nozzle Outer Profile ---------------------
!
csys,CSnoz0
!
x_max=sqrt(l_noz/A1)
x_inc=x_max/n_pts
x=0
!
*do,i,0,n_pts,1
    z=A1*x**2
    !
    k,,x,0,z
    !
    x=x+x_inc
*enddo
!
init_entity_num
!
ksel,s,kp,,km-n_pts,km
bsplin,all,,,,,,0,0
init_entity_num
!
ksel,s,kp,,km-n_pts,km
ksel,u,kp,,km
ksel,u,kp,,km-n_pts
!
kdele,all
!
init_entity_num
!
k,,0,0,-1
k,,0,0,l_noz
k,,1.1*x_max,0,l_noz
!
create_lines,,1
isolate_lines
profile_close
!
create_area-lines
!
w_cch=cch_rat1*x_max
!
k,,w_cch/2,0,w_cch/2
k,,0,0,w_cch/2
k,,0,0,-w_cch/2
k,,w_cch/2,0,-w_cch/2
!
create_lines
!
ksel,s,kp,,km
ln_fillet,,R01
!
create_area-lines
!
aadd,am-1,am
!
init_entity_num
!
isolate_areas
revolve_aset,CSnoz0,360,2
!
cm,v_out,volu
!
lplot
!
!---------------------- 3.2  Nozzle Journals ---------------------
!
csys,CSnoz0
!
create_cs,0,,0,0,0,0,0,0,90
!
create_circle,D01/2
!
create_area-lines
a_journ=am
!
asel,s,area,,a_journ
extrude_aset,L01/2
!
asel,s,area,,a_journ
extrude_aset,-L01/2
!
cmsel,s,v_out
vsel,a,volu,,vm-1,vm,1
vadd,all
!
init_entity_num
v_out=vm
!
lplot
!
!---------------------- 3.3  Nozzle Inner Profile ---------------------
!
csys,CSnoz0
!
x_max=sqrt((l_noz-C2)/A2)
x_inc=x_max/n_pts
x=0
!
*do,i,0,n_pts,1
    z=A2*x**2+C2
    !
    k,,x,0,z
    !
    x=x+x_inc
*enddo
!
init_entity_num
!
ksel,s,kp,,km-n_pts,km
bsplin,all,,,,,,0,0
init_entity_num
!
ksel,s,kp,,km-n_pts,km
ksel,u,kp,,km
ksel,u,kp,,km-n_pts
!
kdele,all
!
init_entity_num
!
k,,0,0,-1
k,,0,0,l_noz
k,,1.1*x_max,0,l_noz
!
create_lines,,1
isolate_lines
profile_close
!
create_area-lines
!
w_cch=cch_rat2*x_max
!
k,,w_cch/2,0,w_cch/2
k,,0,0,w_cch/2
k,,0,0,-w_cch/2
k,,w_cch/2,0,-w_cch/2
!
create_lines
!
ksel,s,kp,,km
ln_fillet,,R02
!
create_area-lines
!
aadd,am-1,am
!
init_entity_num
!
isolate_areas
revolve_aset,CSnoz0,360,2
!
vsel,u,volu,,v_out
cm,v_in,volu
!
init_entity_num
!
vsbv,v_out,v_in,,delete,delete
!
init_entity_num
!
lplot
!
finalize_model