!   TITLE:  c552002670_stat_lam.mac (EMAG)
!
!---------------------------------------------------------------------------
!       Rev -   07/06/2000      J. Rinaldi
!---------------------------------------------------------------------------
!
!   This macro creates the stator lamination profile, P/N C552002670
!
!   CSsl0:      Cartesian Stator Lamination Coordinate System
!   CSsl1:  Cylindrical Stator Lamination Coordinate System
!
!
!******************************* TABLE OF CONTENTS **************************
!
!       Model Assembly Preparation
!
!   1.  Parameter List
!   2.  Model Coordinate System
!   3.  Material and Element Type
!   4.  Solid Model
!       4.1  Circular Area Layout
!       4.2  Lamination Spacer Layout
!
!
!***************************** END Table of Contents ************************
!
!
!
!************************** Model Assembly Preparation **********************
!
/prep7
!
/dev,font,1,charter,bold,r,24,,,,,
!
/title,Stator Lamination C552002670
!
/VIEW, 1 ,,,1
/ANG, 1
!
/units,bin      ! British/inch units
!
!**************************** End Model Preparation *************************
!
!
!
!<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>>

                    d_out=1.011     ! lamination outer diameter
                    d_mid=0.892     ! lamination mid-circle diameter
                    d_in=0.4        ! lamination inner diameter  0.369
                    w_spacer=0.078/3    ! lamination spacer width

                    ang_intr=10     ! lamination spacer angle interval
                    lam_fil=.015    ! lamination fillet  0.015

!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!
!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate Systems >>>>>>>>>>>>>>>>>>>>
!
acs=0       ! active coordinate system (c.s.) locating model c.s.
!
!------------- Specify model c.s. center location --------------
!
xc=0.00     ! x (r) coordinate center
yc=0.00     ! y (theta) coordinate center
zc=0.00     ! z (z or phi) coordinate center
!
!------------- Specify model c.s. orientation ------------------
!
rotxy=0     ! 1st rotation angle (+x to +y)
rotyz=0     ! 2nd rotation angle (+y to +z)
rotzx=0     ! 3rd rotation angle (+z to +x)
!
!---------------------------------------------------------------
!
csn=11
!
csn=csn+1
CSsl0=csn       ! cartesian model c.s. number
!
csn=csn+1
CSsl1=csn       ! cylindrical model c.s. number
!
csys,acs
clocal,CSsl0,0,xc,yc,zc,rotxy,rotyz,rotzx
!
csys,acs
clocal,CSsl1,1,xc,yc,zc,rotxy,rotyz,rotzx
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!
!<<<<<<<<<<<<<<<<<<<< 3.  Material and Element Type >>>>>>>>>>>>>>>>>>>
!
!
!---------------------- Material Properties --------------------
!
UIMP,1,EX, , ,10.3E6,       ! Material:  XXXXX
UIMP,1,DENS, , ,.097/386.4,
UIMP,1,NUXY, , ,0.33,
!UIMP,XX,XXXX
!
!---------------------------------------------------------------
!
!
!------------------------- Element Type ------------------------
!
!ET,XX,XXXXX
!ET,XX,XXXXX
!
!---------------------------------------------------------------
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 4.  Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!
!---------------------- 4.1  Circular Area Layout ---------------------
!
csys,CSsl1
    !
    k,,0,0,0
    k,,d_out/2,0,0
    k,,d_out/2,ang_intr,0
    !
    *get,km,kp,,num,max
    !
    a,km-2,km-1,km
    !
    *get,lm,line,,num,max
    *get,am,area,,num,max
    !
    !
    k,,0,0,0
    k,,d_mid/2,0,0
    k,,d_mid/2,ang_intr,0
    !
    *get,km,kp,,num,max
    !
    a,km-2,km-1,km
    !
    *get,lm,line,,num,max
    *get,am,area,,num,max
    !
    !
    asba,am-1,am,,delete,delete
    !
*get,km,kp,,num,max
numstr,kp,km
!
*get,lm,line,,num,max
numstr,line,lm
!
*get,am,area,,num,max
numstr,area,am
!
lplot
!
!-------------------- 4.2  Lamination Spacer Layout -------------------
!
pi=acos(-1)
phi=(asin((w_spacer/2)/(d_in/2)))*(180/pi)
!
csys,CSsl1
    !
    k,,d_in/2,0,0
    k,,d_in/2,phi,0
    !
    *get,km,kp,,num,max
    !
    l,km-1,km
    !
    *get,lm,line,,num,max
    !
    !
csys,CSsl0
    !
    k,,(d_mid/2)+0.5*(d_out/2-d_mid/2),w_spacer/2,0
    k,,(d_mid/2)+0.5*(d_out/2-d_mid/2),0,0
    !
    *get,km,kp,,num,max
    !
    lstr,km-2,km-1
    lstr,km-1,km
    lstr,km,km-3
    !
    *get,lm,line,,num,max
    !
    !
    al,lm-3,lm-2,lm-1,lm
    !
*get,km,kp,,num,max
numstr,kp,km
!
*get,lm,line,,num,max
numstr,line,lm
!
*get,am,area,,num,max
numstr,area,am
!
lplot
!
csys,CSsl1
    !
    k,,d_in/2,ang_intr-phi,0
    k,,d_in/2,ang_intr,0
    !
    *get,km,kp,,num,max
    !
    l,km-1,km
    !
    *get,lm,line,,num,max
    !
csys,CSsl0
    !
    csn=csn+1
    symcs=csn
    clocal,symcs,0,0,0,0,ang_intr,0,0
    !
    k,,(d_mid/2)+0.5*(d_out/2-d_mid/2),0,0
    k,,(d_mid/2)+0.5*(d_out/2-d_mid/2),-w_spacer/2,0
    !
    *get,km,kp,,num,max
    !
    lstr,km-2,km-1
    lstr,km-1,km
    lstr,km,km-3
    !
    *get,lm,line,,num,max
    !
    al,lm-3,lm-2,lm-1,lm
    !
    *get,am,area,,num,max
    !
    aadd,am-2,am-1,am
    !
*get,km,kp,,num,max
numstr,kp,km
!
*get,lm,line,,num,max
numstr,line,lm
!
*get,am,area,,num,max
numstr,area,am
!
asel,s,area,,am
lsla,s
cm,lset1,line
!
adele,am
!
lfillt,lm,lm-1,lam_fil
lfillt,lm,lm-2,lam_fil
!
*get,km,kp,,num,max
*get,lm,line,,num,max
!
allsel,all
!
lsel,s,line,,lset1
lsel,a,line,,lm-1,lm
!
al,all
!
*get,km,kp,,num,max
numstr,kp,km
!
*get,lm,line,,num,max
numstr,line,lm
!
*get,am,area,,num,max
numstr,area,am
!
lplot
!
!
n_intr=90/ang_intr
!
*do,i,1,n_intr-1,1
    !
    csys,symcs
    !
    arsym,y,am,,,,0,0
    !
    *get,km,kp,,num,max
    numstr,kp,km
    !
    *get,lm,line,,num,max
    numstr,line,lm
    !
    *get,am,area,,num,max
    numstr,area,am
    !
    clocal,symcs,0,0,0,0,ang_intr,0,0
*enddo
!
csys,CSsl0
!
*do,j,1,n_intr,1
    *if,j,eq,1,then
        asel,s,area,,am-(n_intr-j)
    *else
        asel,a,area,,am-(n_intr-j)
    *endif
*enddo
!
arsym,x,all,,,,0,0
!
*get,km,kp,,num,max
numstr,kp,km
!
*get,lm,line,,num,max
numstr,line,lm
!
*get,am,area,,num,max
numstr,area,am
!
*do,j,1,2*n_intr,1
    *if,j,eq,1,then
        asel,s,area,,am-(2*n_intr-j)
    *else
        asel,a,area,,am-(2*n_intr-j)
    *endif
*enddo
!
arsym,y,all,,,,0,0
!
*get,km,kp,,num,max
numstr,kp,km
!
*get,lm,line,,num,max
numstr,line,lm
!
*get,am,area,,num,max
numstr,area,am
!
csys,CSsl1
asel,s,loc,x,d_in/2,d_out/2
!
aadd,all
!
*get,km,kp,,num,max
numstr,kp,km
!
*get,lm,line,,num,max
numstr,line,lm
!
*get,am,area,,num,max
numstr,area,am
!
lplot