!   TITLE:  stator_lam.mac (PART_LIB)
!
!
!   This macro creates a stator lamination stack / coil wire assembly:
!
!               P/N XXXXXXXXXXXX
!
!       UPPER ASSEMBLY:     XXXXXXXXXXXXX
!
!   ctslcs:         Cartesian Stator Lamination Coordinate System
!   cylslcs:        Cylindrical Stator Lamination Coordinate System
!   sphslcs:        Spherical Stator Lamination Coordinate System
!
!
!******************************* TABLE OF CONTENTS **************************
!
!       Model Assembly Preparation
!
!   1.  Parameter List
!   2.  Model Coordinate System
!   2.  Material and Element Type
!   4.  Solid Model
!   5.  FE model
!
!
!************************** Model Assembly Preparation **********************
!
!
!
mod_assemb_prep                 ! prepare database for assembly
!
!
/title,XXXXXXXX
!
!
/units,bin                      ! British/inch units
!
!
!<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>>


        r1=1.225/2-0.1
        r2=1.800/2+0.1

        r3=1.225/2              ! wire inner radius
        r4=1.800/2              ! wire outer radius

        l1=0.600                ! lam. axial length
        l2=0.350                ! wire bundle width

        w_stat=0.29176          ! stator lamination assembly weight


!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate Systems >>>>>>>>>>>>>>>>>>>>
!
!
mod_coord_system,cthscs,0.00,0.00,0.00,0,0,0
!
!
!<<<<<<<<<<<<<<<<<<<< 3.  Material and Element Type >>>>>>>>>>>>>>>>>>>
!
!
!
UIMP,3,EX, , ,28.5E6,       ! Material:  ASTM A582 TYPE 416 COND A
UIMP,3,DENS, , ,.28/386.4,
UIMP,3,NUXY, , ,0.27,
!alpha=5.9E-6 in/in/F
!
!
ET,1,SOLID95
ET,2,SHELL93
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 4.  Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
csys,ctslcs
!
k,,0,0,0        ! create axis rotation points
k,,0,0,1
!
init_entity_num
!
axp1=km-1
axp2=km
!
k,,0,r1,l1      ! points 1 thru 4
k,,0,r1,0
k,,0,r2,0
k,,0,r2,l1
!
create_area
!
lplot
!
!
csys,ctslcs
!
k,,0,r3,l1      ! points 1 thru 4
k,,0,r4,l1
k,,0,r4,l1+l2
k,,0,r3,l1+l2
!
create_lines
!
lfillt,lm-2,lm-1,0.2*((r2-r3)-(r2-r4))
lfillt,lm-1,lm,0.2*((r2-r3)-(r2-r4))
!
create_area-lines
!
k,,0,r4,0       ! points 5 thru 8
k,,0,r3,0
k,,0,r3,-l2
k,,0,r4,-l2
!
create_lines
!
lfillt,lm-2,lm-1,0.2*((r2-r3)-(r2-r4))
lfillt,lm-1,lm,0.2*((r2-r3)-(r2-r4))
!
create_area-lines
!
aglue,am-2,am-1,am
!
init_entity_num
!
lplot
!
!
revolve_aset,axp1,axp2,4*rtangle,4
!
kdele,axp1,axp2
!
init_entity_num
!
lplot
!
!
vsum
!
*get,v_stat,volu,,volu      ! total volume of assembly
!
rho_stat=w_stat/v_stat      ! density of assembly material
!
!
finalize_model_data