!   TITLE:  mirror_support.mac
!
!---------------------------------------------------------------------------
!       Rev -   03/14/2000  J. Rinaldi
!---------------------------------------------------------------------------
!
!   This is a macro that builds the AAAV mirror support.  It is current
!   for ANSYS 5.6.  P/N: 400A171
!
!   CSmirsp0:       Cartesian Mirror Support Coordinate System
!   CSmirsp1:   Cylindrical Mirror Support Coordinate System
!
!
!******************************* TABLE OF CONTENTS **************************
!
!      Model Assembly Preparation
!
!   1. Model Coordinate System
!   2. Material and Element Type
!   3. Parameter List
!   4. Solid Model
!       4.1  Create Initial Profile Volume
!       4.2  Create Top Surface Cut
!       4.3  Create Rib Pockets
!       4.4  Modify Bottom Surface
!           4.4.1  Beveled Rectanglular Cut
!           4.4.2  Beveled Surface Cut
!           4.4.3  Left Beveled Surface Material Cut
!       4.5  Center Hole Bores
!           4.5.1  Center Hole
!           4.5.2  Counterbore
!       4.6  Elevation Axis Journals
!           4.6.1  Left Journal
!           4.6.2  Right Journal
!   5.  FEA Model
!       5.1  Prepare Volumes
!           5.1.1  Long Axis Cut
!           5.1.2  Left Journal Volumes
!           5.1.3  Left Side Rib Volumes
!           5.1.4  Right Side Rib Volumes
!           5.1.5  Short Axis Cut
!           5.1.6  Right Journal Volumes
!           5.1.7  Back Surface Cuts
!           5.1.8  Back Contour Cuts
!           5.1.9  Rib Pocket Contours
!
!***************************** END Table of Contents ************************
!
!
!
!************************** Model Assembly Preparation **********************
!
!
mod_assemb_prep,1               ! prepare database for assembly
!
!
/title,AAAV Mirror Support Construction
!
!
/units,bin                  ! British/inch units
!
!
!**************************** End Model Preparation *************************
!
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>
!
!
mirror_support.parm
!
!
!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate Systems >>>>>>>>>>>>>>>>>>>>
!
!
!mod_coord_system,0,1.00,2.00,0.00,180,0,0
mod_coord_system,acs,xc,yc,zc,rotxyc,rotyzc,rotzxc
!
create_cs,0,'CSms_cg0',0,0,0,h11,0,0,0      ! Want CG relative to here
!
!<<<<<<<<<<<<<<<<<<<< 3.  Material and Element Type >>>>>>>>>>>>>>>>>>>
!
!
!---------------------- Material Properties --------------------
!
UIMP,1,EX, , ,10.3E6,       ! Aluminum A356  2
UIMP,1,DENS, , ,.097/386.4,
UIMP,1,NUXY, , ,0.33,

!UIMP,1,EX, , ,30.0E6,      ! Beralcast 191 or 363
!UIMP,1,DENS, , ,.078/386.4,
!UIMP,1,NUXY, , ,0.2,
!
!
!---------------------------------------------------------------
!
!
!------------------------- Element Type ------------------------
!
ET,1,SOLID95
ET,2,SHELL93
!
!---------------------------------------------------------------
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 4. Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!----------------- 4.1  Create Initial Profile Volume ----------------
!
csys,CSmirsp0
!
k,,l4,h1,0
k,,-l3,h1,0
k,,-l3,-h2,0
k,,l4,-h2,0
!
init_entity_num
!
lstr,km-3,km-2
lstr,km,km-1
!
k,,l2,h1,0
k,,l2,-h2,0
k,,-l1,h1,0
k,,-l1,-h2,0
!
init_entity_num
!
lstr,km-3,km-2
lstr,km,km-1
!
create_cs,1,,1,km-7,km-3,km-4           ! create 1st beveled line
k,,l2,phi1,0
lstr,km+1,km-7
!
create_cs,1,,1,km-6,km-1,km-5           ! create 2nd beveled line
k,,l1,phi1,0
lstr,km+2,km-6
!
create_cs,1,,1,km-5,km,km-6             ! create 3rd beveled line
k,,l1,phi1,0
lstr,km+3,km-5
!
create_cs,1,,1,km-4,km-2,km-7           ! create 4th beveled line
k,,l2,phi1,0
lstr,km+4,km-4
!
init_entity_num
!
lsbl,lm-3,lm-5,,delete,keep
lsbl,lm,lm-5,,delete,delete
!
lstr,km+1,km+2
!
lsbl,lm-2,lm-4,,delete,keep
lsbl,lm-1,lm-4,,delete,delete
!
lstr,km+3,km+4
!
init_entity_num
!
ldele,lm-7,lm-1,2,1
!
init_entity_num
!
a,km-15,km-14,km-1,km,km-13,km-12,km-2,km-3
!
init_entity_num
!
voffst,am,extd1
!
init_entity_num
!
lplot
!
!----------------------- 4.2  Create Top Surface Cut -------------------
!
csys,CSmirsp0
!
create_cs,0,'CStprof',0,0,-(1.2*h2),0,0,90,0
!
k,,-l1,h4,0             ! create horizontal surface line
k,,l2,h4,0
!
create_lines,,1
!
create_cs,0,,1,-l1,h3,0,0,0,0
!
k,,l1,-phi2,0           ! create left beveled surface line
k,,l1,180-phi2,0
!
create_lines,,1
!
csys,CStprof                ! create right beveled surface line
!
create_cs,0,,1,l2,h3,0,0,0,0
!
k,,l2,180+phi2,0
k,,l2,phi2,0
!
create_lines,,1
!
lsbl,lm-1,lm-2,,delete,keep
    csys,CSmirsp0       ! left partition distance, for meshing purposes
    pdl=kx(km+1)
    csys,csn
lsbl,lm,lm-2,,delete,delete
    csys,CSmirsp0       ! right partition distance, for meshing purposes
    pdr=kx(km+2)
    csys,csn
lstr,km+1,km+2
ldele,lm+1,lm+3,2,1
!
csys,CStprof
!
k,,kx(km-2),2*extd1,0
k,,kx(km),2*extd1,0
!
init_entity_num
!
lstr,km,km-1
lstr,km-1,km-6
lstr,km,km-4
!
init_entity_num
!
a,km,km-1,km-6,km-3,km-2,km-4
!
init_entity_num
!
voffst,am,-(1.5*(h1+h2))
!
init_entity_num
!
vsbv,vm-1,vm,,delete,delete
!
init_entity_num
!
lplot
!
!------------------------ 4.3  Create Rib Pockets ---------------------
!
csys,CSmirsp0
!
k,,kx(km-1),ky(km-1),2*extd1
k,,kx(km-7),ky(km-7),2*extd1
k,,0,h6,2*extd1
k,,1.2*l2,h6,2*extd1
!
init_entity_num
!
lstr,km,km-1
lstr,km-2,km-3
!
create_cs,1,,0,km-3,km-2,km-1
!
k,,-2.00,h5,0
k,,l2,h5,0
!
init_entity_num
!
lstr,km-1,km
!
csys,CSmirsp0
!
k,,l5,0,2*extd1
k,,l5,h1,2*extd1
!
init_entity_num
!
lstr,km,km-1
!
init_entity_num
!
lsbl,lm-1,lm,,delete,keep
lsbl,lm-3,lm,,delete,delete
!
init_entity_num
!
lstr,km-1,km
!
create_cs,1,,1,km,km-7,km-1
!
k,,l2,phi3,0
!
lstr,km,km+1
!
csys,CSmirsp0           ! formation of rib lines
!
k,,l6,h1,2*extd1
k,,l6,-h1,2*extd1
!
init_entity_num
!
lstr,km-1,km
!
k,,l7,h1,2*extd1
k,,l7,-h1,2*extd1
!
init_entity_num
!
lstr,km-1,km
!
k,,l8,h1,2*extd1
k,,l8,-h1,2*extd1
!
init_entity_num
!
lstr,km-1,km
!
k,,l9,h1,2*extd1
k,,l9,-h1,2*extd1
!
init_entity_num
!
lstr,km-1,km
!
k,,l10,h1,2*extd1
k,,l10,-h1,2*extd1
!
init_entity_num
!
lstr,km-1,km
!
init_entity_num
!
lsel,selv,line,,lm-4,lm,1
lsel,a,line,,lm-10
!
ln_part_ln,lm-10
!
init_entity_num
!
lstr,km,km-1
lstr,km-2,km-3
lstr,km-4,km-17
!
init_entity_num
!
lsel,selv,line,,lm-11,lm-3,2
lsel,a,line,,lm-18
!
ln_part_ln,lm-18
!
init_entity_num
!
lstr,km-4,km-21
lstr,km-2,km-3
lstr,km,km-1
!
init_entity_num
!
lfillt,lm-15,lm-3,rfil      ! fillet rib lines
lfillt,lm-3,lm,rfil
lfillt,lm,lm-5,rfil
lfillt,lm-5,lm-15,rfil
!
lfillt,lm-14,lm-7,rfil
lfillt,lm-7,lm-1,rfil
lfillt,lm-1,lm-9,rfil
lfillt,lm-9,lm-14,rfil
!
lfillt,lm-13,lm-11,rfil
lfillt,lm-11,lm-2,rfil
lfillt,lm-2,lm-32,rfil
lfillt,lm-32,lm-13,rfil
!
init_entity_num
!
a,km-23,km-22,km-21,km-20,km-19,km-18,km-17,km-16
a,km-15,km-14,km-13,km-12,km-11,km-10,km-9,km-8
a,km-7,km-6,km-5,km-4,km-3,km-2,km-1,km
!
init_entity_num
!
asel,selv,area,,am-2,am,1
lsla,s
cm,lset1,line
lsel,all
lsel,u,line,,lset
!
csys,CSmirsp0
!
lsel,selv,loc,z,2*extd1
lsel,u,line,,lset1
!
ldele,all
!
lsel,all
lsel,u,line,,lset
!
arsym,y,am-2,am,1,,1,0
!
create_cs,0,'CSsymm0',0,-(abs(llext-lrext)/2),0,0,0,0,0
!
arsym,x,am-2,am+3,1,,1,0        ! create areas by 2nd reflection
!
init_entity_num
!
*do,i,am-11,am,1
    voffst,i,-((2*extd1)-h7),,
*enddo
!
init_entity_num
!
vsel,selv,volu,,vm-11,vm
cm,vset1,volu
vsel,all
vsel,u,volu,,vset
!
init_entity_num
!
vsbv,vm-12,vset1,,delete,delete
!
init_entity_num
!
lplot
!
!------------------------ 4.4  Modify Bottom Surface ---------------------
!
!
!--------------------- 4.4.1  Beveled Rectanglular Cut -------------------
!
csys,CSmirsp0
!
rec_bore,7,0,0,0,0,0,2*rtangle,l12,l11,abs(h2-h9),abs(h2-h8),dcut,phi4
!
init_entity_num
!
lplot
!
!------------------------ 4.4.2  Beveled Surface Cut ---------------------
!
csys,CSmirsp0
!
create_cs,0,'CSbevsf0',0,-(1.5*l1),-(1.5*h2),(h10+chmf1),90,-90,0
!
k,,0,0,0
k,,(1.5*h2)-h2,0,0
k,,(1.5*h2)-h2+chmf1,chmf1,0
k,,(1.5*h2)-h2+l13,chmf1,0
k,,(1.5*h2)-h2+l13,chmf1+h10,0
!
k,,(1.5*h2)-h2+l14,chmf1+h10,0
k,,(1.5*h2)-h2+l14,chmf1,0
k,,(1.5*h2)+h1-chmf2,chmf1,0
!
init_entity_num
!
lstr,km-7,km-6
lstr,km-6,km-5
lstr,km-5,km-4
lstr,km-3,km-2
lstr,km-1,km
!
create_cs,1,,1,km-4,km-1,km-3
!
k,,3*h10,phi5,0
lstr,km-4,km+1
!
create_cs,1,,1,km-1,km,km-2             ! create 2nd beveled line
!
k,,3*h10,180-phi5,0
lstr,km-1,km+2
!
init_entity_num
!
lsbl,lm-1,lm-3,,delete,keep
lsbl,lm,lm-3,,delete,delete
!
init_entity_num
!
lstr,km-1,km
!
ldele,lm-2,lm,2,1
!
init_entity_num
!
csys,CSbevsf0
!
k,,(1.5*h2)+h1,chmf1-chmf2,0
k,,(1.5*h2)+(1.5*h1),chmf1-chmf2,0
k,,(1.5*h2)+(1.5*h1),2*(h10+chmf1),0
k,,0,2*(h10+chmf1),0
!
init_entity_num
!
lstr,km,km-1
lstr,km-1,km-2
lstr,km-2,km-3
lstr,km-3,km-8
lstr,km,km-15
!
init_entity_num
!
lsel,selv,line,,lm-5,lm,1
lsel,a,line,,lm-9,lm-7,2
lsel,a,line,,lm-16,lm-14,1
lsel,a,line,,lm-12
!
al,all
!
lsel,all
lsel,u,line,,lset
!
init_entity_num
!
voffst,am,-(1.75*(l1+l2))
!
init_entity_num
!
vsbv,vm-1,vm,,delete,delete
!
init_entity_num
!
lplot
!
!--------------- 4.4.3  Left Beveled Surface Material Cut ------------
!
csys,CSmirsp0
!
k,,-l15,-h2+l16,-0.5
k,,-l15,-(1.2*h2),-0.5
k,,-(1.2*l1),-(1.2*h2),-0.5
k,,-(1.2*l1),-h2+l16,-0.5
!
init_entity_num
!
lstr,km,km-1
lstr,km-1,km-2
lstr,km-2,km-3
lstr,km-3,km
!
init_entity_num
!
lfillt,lm,lm-1,rrad
!
init_entity_num
!
a,km,km-1,km-2,km-3,km-4
!
init_entity_num
!
voffst,am,0.5+h10
!
init_entity_num
!
vsbv,vm-1,vm,,delete,delete
!
init_entity_num
!
lplot
!
!------------------------ 4.5  Center Hole Bores ---------------------
!
!
!
!-------------------------- 4.5.1  Center Hole -----------------------
!
csys,CSmirsp1
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)
    k,,dhole/2,i,-extd1
*enddo
!
create_area,1
!
voffst,am,2*extd1
!
init_entity_num
!
vsbv,vm-1,vm,,delete,delete
!
init_entity_num
!
!------------------------ 4.5.2  Counterbore ------------------------
!
/out,ms,debug

csys,CSmirsp1
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)
    k,,dbore1/2,i,2*h4
*enddo
!
/out

create_area,1
!
/out,ms,debug,,append

!voffst,am,-(((2*h4)-h11)+boredpth)
voffst,am,-2*h4+boredpth

/out

!
init_entity_num
!
vsbv,vm-1,vm,,delete,delete
!
init_entity_num
!
lplot
!
!-------------------- 4.6  Elevation Axis Journals -------------------
!
!
!
!---------------- 4.6.1  Left (Resolver) Journal ----------------------
!
csys,CSmirsp0
!
create_cs,0,'CSm_res1',1,-l1,0,h11,0,0,-90
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)
    k,,dlj/2,i,0
*enddo
!
create_area,1
!
voffst,am,extd2
!
init_entity_num
!
vadd,vm-1,vm
!
init_entity_num
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)     ! machine inner bore
    k,,dbore2/2,i,1.2*extd2
*enddo
!
create_area,1
!
voffst,am,-(0.2*extd2+extd3)
!
init_entity_num
!
vsbv,vm-1,vm,,delete,delete
!
init_entity_num
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)     ! machine outer surface
    k,,1.2*dlj/2,i,1.2*extd2
*enddo
!
create_area,1
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)
    k,,dos/2,i,1.2*extd2
*enddo
!
create_area,1
!
asba,am-1,am,,delete,delete
!
init_entity_num
!
voffst,am,-(0.2*extd2+extd4)
!
init_entity_num
!
vsbv,vm-1,vm,,delete,delete
!
init_entity_num
!
lplot
!
!----------- 3.6.1.1  Resolver Journal Bending modal ----------------------
!
*if,f_bend,eq,1,then
!
csys,CSm_res1
lstr,kp(dbore2/2,0,extd2-extd3),kp(dbore2/2,180,extd2-extd3)
lstr,kp(dbore2/2,90,extd2-extd3),kp(dbore2/2,270,extd2-extd3)
!
asel,s,loc,z,extd2-extd3
!
asbl,all,lm+1,,delete,delete
!
asel,s,loc,z,extd2-extd3
!
asbl,all,lm+2,,delete,delete
!
init_entity_num
!
lplot
!
*endif
!
!------------------- 4.6.2  Right (Torquer) Journal ----------------------
!
csys,CSmirsp0
!
create_cs,0,'CSm_trq1',1,l2,0,h11,0,0,90
!
create_circle,drj/2                         ! 1st volume
!
create_vol-lines,extd5
!
create_circle,d1/2,extd5                    ! 2nd volume
!
*if,f_bend,eq,1,then
!
!------------- 3.6.2.1 Modifyied model for bending modal analysis ---------
!
vadd,vm-1,vm
!
init_entity_num
!
create_vol-lines,l18/2                      ! Bearing Journal extrude
!
vsel,s,volu,,vm
!
ang_cut_pln,0,0
!
vsel,s,volu,,vm-1,vm
!
ang_cut_pln,90,0
!
vadd,all
!
init_entity_num
!
*else
!
!------------ 3.6.2.2 Full model for tosional modal analysis --------------
!
create_vol-lines,l18                        ! Bearing Journal extrude
!
create_circle,d2/2,l18+extd5                ! 3rd volume
!
create_vol-lines,l19
!
create_circle,d3/2,l19+l18+extd5            ! 4th volume
!
create_vol-lines,l20
!
create_circle,d4/2,l20+l19+l18+extd5        ! 5th volume
!
create_vol-lines,l21
!
vadd,vm-5,vm-4,vm-3,vm-2,vm-1,vm
!
init_entity_num
!
*endif
!
lplot
!
*go,:skip
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 5.  FEA Model >>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!----------------------- 5.1  Prepare Volumes ----------------------
!
!
!----------------------- 5.1.1  Long Axis Cut ----------------------
!
csys,CSm_res1
!
k,,1.5*dlj/2,0,(1.2*extd2)
k,,1.5*dlj/2,0,-(1.5*(l1+l2))
k,,1.5*dlj/2,2*rtangle,-(1.5*(l1+l2))
k,,1.5*dlj/2,2*rtangle,(1.2*extd2)
!
init_entity_num
!
lstr,km-3,km-2
lstr,km-2,km-1
lstr,km-1,km
lstr,km,km-3
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vol_part_area,am
!
init_entity_num
!
lplot
!
!-------------------- 5.1.2  Left Journal Volumes ------------------
!
csys,CSm_res1
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)  ! create 1st cutting plane
    k,,1.5*dlj/2,i,l22-l1
*enddo
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vol_part_area,am
!
init_entity_num
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)  ! create volume partitions
    k,,1.5*h1,i,-l23
*enddo
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vol_part_area,am
!
init_entity_num
!
vsel,selv,loc,z,-l23,(l22-l1)
!
vadd,all
!
vsel,all
vsel,u,volu,,vset
!
init_entity_num
!
csys,CSm_res1           ! quarter-cut journal volumes
!
k,,1.5*dlj/2,90,0
k,,1.5*dlj/2,90,1.5*extd2
k,,1.5*dlj/2,270,1.5*extd2
k,,1.5*dlj/2,270,0
!
init_entity_num
!
lstr,km-3,km-2
lstr,km-2,km-1
lstr,km-1,km
lstr,km,km-3
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vsel,selv,loc,z,(extd2-extd3),5*(extd2-extd3)
!
vol_part_area,am
!
init_entity_num
!
lplot
!
!-------------------- 5.1.3  Left Side Rib Volumes ------------------
!
csys,CSmirsp0
!
k,,pdl,1.2*h1,-1.5*extd1
k,,pdl,1.2*h1,1.5*extd1
k,,pdl,-1.2*h2,1.5*extd1
k,,pdl,-1.2*h2,-1.5*extd1
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vsel,selv,loc,x,-(l1-l23),l2
!
vol_part_area,am
!
init_entity_num
!
lplot
!
!-------------------- 5.1.4  Right Side Rib Volumes ------------------
!
csys,CSmirsp0
!
k,,pdr,1.2*h1,-1.5*extd1
k,,pdr,1.2*h1,1.5*extd1
k,,pdr,-1.2*h2,1.5*extd1
k,,pdr,-1.2*h2,-1.5*extd1
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vsel,selv,loc,x,pdl,l2
!
vol_part_area,am
!
init_entity_num
!
lplot
!
!----------------------- 5.1.5  Short Axis Cut ---------------------
!
csys,CSmirsp0
!
k,,0,1.2*h1,-1.5*extd1
k,,0,1.2*h1,1.5*extd1
k,,0,-1.2*h2,1.5*extd1
k,,0,-1.2*h2,-1.5*extd1
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vsel,selv,loc,x,pdl,pdr
!
vol_part_area,am
!
init_entity_num
!
lplot
!
!-------------------- 5.1.6  Right Journal Volumes --------------------
!
csys,CSm_trq1
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)  ! create 1st cutting plane
    k,,1.5*h1,i,-l24
*enddo
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
csys,CSmirsp0
vsel,selv,loc,x,pdr,1.5*l2
csys,CSm_trq1
!
vol_part_area,am
!
init_entity_num
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)
    k,,1.5*drj/2,i,extd5
*enddo
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vsel,selv,loc,z,-l24,5.00
!
vol_part_area,am
!
init_entity_num
!
vsel,selv,loc,z,-l24,extd5
!
vadd,all
!
init_entity_num
!
k,,1.5*drj/2,90,5.00
k,,1.5*drj/2,90,-l24
k,,1.5*drj/2,270,-l24
k,,1.5*drj/2,270,5.00
!
init_entity_num
!
lstr,km-3,km-2
lstr,km-2,km-1
lstr,km-1,km
lstr,km,km-3
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vsel,selv,loc,z,extd5,5.00
!
vol_part_area,am
!
init_entity_num
!
lplot
!
!---------------------- 5.1.7  Back Surface Cuts ---------------------
!
csys,CSmirsp1
!
*do,i,0,3*((pi/2)*dconv),((pi/2)*dconv)   ! create cutting plane
    k,,1.5*l1,i,h10
*enddo
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
csys,CSmirsp0
vsel,selv,loc,x,-(l1-l23),(l2-l24)
!
vol_part_area,am
!
init_entity_num
!
lplot
!
!----------------------- 5.1.8  Back Contour Cuts ----------------------
!
csys,CSmirsp0
!
k,,1.5*l2,1.5*h1,h7
k,,-(1.5*l1),1.5*h1,h7
k,,-(1.5*l1),-(1.5*h2),h7
k,,1.5*l2,-(1.5*h2),h7
!
init_entity_num
!
a,km-3,km-2,km-1,km
!
init_entity_num
!
vsel,selv,loc,x,-(l1-l23),(l2-l24)
!
vol_part_area,am
!
init_entity_num
!
lplot
!






!---------------------- 5.1.9  Rib Pocket Contours ---------------------
!
csys,CSmirsp0
!
!asel,selv,loc,z,h7
!

!-(abs(llext-lrext)/2)
! Create CS's for the extreme positions
csys,CSm_res1
create_cs,0,'CSms_up1',1,0,0,0,-el_up,0,0
!
csys,CSm_res1
create_cs,0,'CSms_dn1',1,0,0,0,el_down,0,0


:skip


finalize_model

save,mirror_support,db