/BATCH
/COM,ANSYS RELEASE  5.4    UP19970828         10:56:52    07/26/1999
/input,menust,tmp     ,,,,,,,,,,,,,,,,,1
/GRA,POWER
/GST,ON
/dev,font,1,charter,bold,r,24,,,,,
/title,AAAV Modal Analysis / AlBe
/units,bin
/prep7
!
!**************************** TABLE OF CONTENTS ***********************
!   
!       1.  Material Properties/Element Types
!       2.      VARIABLE/CONSTANTS LIST
!       3.      Create Global Coordinate System
!   
!       Top Surface Layout
!   
!       4.      First (Mid) Rib Layout
!       5.      2nd Rib Layout
!       6.      3rd Rib Layout
!       7.      4th Rib Layout
!       8.      5th Rib Layout
!       9.  6th Rib Layout
!       10. 7th Rib Layout
!   
!       Begin Center Boss Construction
!   
!       11. Mid (First) Rib, AREA 1
!       12. Mid (First) Rib, AREA 2
!       13. 2nd Rib, AREA 3
!       14. 2nd Rib, AREA 4
!       15. 3rd Rib, AREA 5
!       16. 3rd Rib, AREA 6
!       17. 4th Rib, AREA 7
!       18. 5th Rib, AREA 10
!       19. 6th Rib, AREA 11
!       20. 6th Rib, AREA 12
!       21. 7th Rib, AREA 13
!       22. 7th Rib, AREA 14
!       23. 4th Rib, AREA 8
!       24. 5th Rib, AREA 9
!   
!       Center Boss Volume Lines
!   
!       25. Areas 2 & 3
!       26. Areas 4 & 5
!       27. Areas 6 & 7
!       28. Areas 10 & 11
!       29. Areas 12 & 13
!       30. Areas 14 & 1
!       31. Areas 8 & 9
!   
!       Fill In Rib Gaps
!   
!       32. Rib 1
!       33. Rib 2
!       34. Rib 3
!       35. Rib 6
!       36. Rib 7
!       37. Rib 4
!       38. Rib 5
!           
!********************* Material Properties/Element Types **************
!
ET,1,SOLID95
ET,2,SHELL93
!
UIMP,1,EX, , ,30E6,
UIMP,1,DENS, , ,.101/386.4,
UIMP,1,NUXY, , ,0.3,
!
!************************ VARIABLE/CONSTANTS LIST ******************
!
pi=3.1415927
rconv=pi/180     ! radian conversion factor from degrees
!
w1=.25       ! Rib width (ribs 1,2,4,5,7)
w2=.227      ! Rib width (ribs,3,6)
ribh=.736    ! Rib height
r1=.25       ! Rib fillet radius
r2=.125      ! Front rib fillet radius
fribh=ribh-r1    ! Height of front surface ribs
topthk=.40   ! Top surface thickness
ftopthk=.65  ! Front top surface thickness
cbore1=.25   ! Center counterbore depth 1
cbore2=.852  ! Center counterbore depth 2
rcb1=2.0005  ! Radius of counterbore 1
rcb2=1.5595  ! Radius of counterbore 2
rbore=1.1565     ! Center bore hole radius
rbout=2.250  ! Center boss outer radius
!   
netop=2      ! # elements through top surface thickness
neftop=4     ! # elements through front top surface thickness
nefill=4     ! # elements in fillet radius (must be even #)
neffill=nefill-1 ! # elements in front rib fillet
neL1=4       ! # elements in rib 1 length
neL2=10      ! # elements in rib 2 length
neL3=11      ! # elements in rib 3 length
neL4=8       ! # elements in rib 4 length
neL5=7       ! # elements in rib 5 length
neL6=10      ! # elements in rib 6 length
neL7=9       ! # elements in rib 7 length
!
L1=1.405     ! Length of rib 1
L2=3.5725    ! Length of rib 2
L3=3.8099    ! Length of rib 3
L4=2.926     ! Length of rib 4
L5=2.462     ! Length of rib 5
L6=3.647     ! Length of rib 6
L7=3.188     ! Length of rib 7
!
!   RIB 1 Point Coordinates
!
r1p1x=-2.4717    ! Point 1 x
r1p1y=-1.136     ! Point 1 y
r1p1z=.125   ! Point 1 z
!
r1p2x=-2.4717    ! Point 2 x
r1p2y=-1.136     ! Point 2 y
r1p2z=-.125  ! Point 2 z
!
r1p3x=-2.4717    ! Point 3 x
r1p3y=-.650  ! Point 3 y
r1p3z=.125   ! Point 3 z
!
!   RIB 2 Point Coordinates
!
r2p1x=-1.417     ! Point 1 x
r2p1y=-1.136     ! Point 1 y
r2p1z=-2.029     ! Point 1 z
!
r2p2x=-1.2051    ! Point 2 x
r2p2y=-1.136     ! Point 2 y
r2p2z=-2.1616    ! Point 2 z
!
r2p3x=-1.417     ! Point 3 x
r2p3y=-.650  ! Point 3 y
r2p3z=-2.029     ! Point 3 z
!
!   RIB 3 Point Coordinates
!
r3p1x=-.102  ! Point 1 x
r3p1y=-1.136     ! Point 1 y
r3p1z=-2.472     ! Point 1 z
!
r3p2x=.125   ! Point 2 x
r3p2y=-1.136     ! Point 2 y
r3p2z=-2.472     ! Point 2 z
!
r3p3x=-.102  ! Point 3 x
r3p3y=-.650  ! Point 3 y
r3p3z=-2.472     ! Point 3 z
!
!   RIB 4 Point Coordinates
!
r4p1x=.8809  ! Point 1 x
r4p1y=-1.136     ! Point 1 y
r4p1z=-2.326     ! Point 1 z
!
r4p2x=1.11   ! Point 2 x
r4p2y=-1.136     ! Point 2 y
r4p2z=-2.226     ! Point 2 z
!
r4p3x=.8809  ! Point 3 x
r4p3y=-.650  ! Point 3 y
r4p3z=-2.326     ! Point 3 z
!
!   RIB 5 Point Coordinates
!
r5p1x=.9337  ! Point 1 x
r5p1y=-1.136     ! Point 1 y
r5p1z=2.306  ! Point 1 z
!
r5p2x=1.1604     ! Point 2 x
r5p2y=-1.136     ! Point 2 y
r5p2z=2.2007     ! Point 2 z
!
r5p3x=.9337  ! Point 3 x
r5p3y=-.650  ! Point 3 y
r5p3z=2.306  ! Point 3 z
!
!   RIB 6 Point Coordinates
!
r6p1x=.148   ! Point 1 x
r6p1y=-1.136     ! Point 1 y
r6p1z=2.475  ! Point 1 z
!
r6p2x=-.102  ! Point 2 x
r6p2y=-1.136     ! Point 2 y
r6p2z=2.475  ! Point 2 z
!
r6p3x=.148   ! Point 3 x
r6p3y=-.650  ! Point 3 y
r6p3z=2.475  ! Point 3 z
!
!   RIB 7 Point Coordinates
!
r7p1x=-1.2707    ! Point 1 x
r7p1y=-1.136     ! Point 1 y
r7p1z=2.1237     ! Point 1 z
!
r7p2x=-1.4785    ! Point 2 x
r7p2y=-1.136     ! Point 2 y
r7p2z=1.9847     ! Point 2 z
!
r7p3x=-1.2707    ! Point 3 x
r7p3y=-.650  ! Point 3 y
r7p3z=2.1237     ! Point 3 z
!
/com, ****************** Create Global Coordinate System ***************
!
csys,0
k,,0,0,0
k,,0,0,-.5
k,,-.5,0,0
*get,km,kp,,num,max
csn=10
csn=csn+1
gcs=csn
cskp,gcs,1,km-2,km-1,km,,
!
/com, ************************ Top Surface Layout ********************** 
!
csys,0
k,,-3.877,0,-4.786
k,,-.999,0,-6.585
k,,3.207,0,-6.585
k,,3.207,0,6.115
k,,-.999,0,6.115
k,,-3.877,0,4.316
*get,km,kp,,num,max
l,km-5,km-4
l,km-4,km-3
l,km-3,km-2
l,km-2,km-1
l,km-1,km
l,km,km-5
!
/com, ************************** End Top Surface ***********************   
!
/com, ******************** First (Mid) Rib Layout ******************* 
!
csys,0
!
k,,r1p1x,r1p1y,r1p1z        ! 3 keypoints forming planar surface 
k,,r1p2x,r1p2y,r1p2z        ! to build rib cross-sections 
k,,r1p3x,r1p3y,r1p3z
!
*get,km,kp,,num,max     ! creation of c.s. for rib in planar
csn=csn+1           ! surface
rib1cs=csn
cskp,rib1cs,0,km-2,km-1,km,,
!
kbetw,km-2,km-1,km+1,ratio,0.5  ! center keypoint of rib cross-section
!
cskp,rib1cs,0,km+1,km-1,km  ! move planar surface c.s. to rib center
!
k,,0,ribh-r1,0
k,,w1/2,0,0
k,,w1/2,ribh-r1,0       ! remaining rib cross-section keypoints 
k,,w1/2+r1,ribh,0
k,,w1/2+r1,ribh+topthk,0
k,,0,ribh+topthk,0
!
*get,km,kp,,num,max     ! form straight lines of rib
!
l,km-6,km-4
l,km-4,km-3
l,km-2,km-1
l,km-1,km
l,km,km-5
l,km-5,km-6
l,km-5,km-3
!
k,,w1/2+r1,ribh-r1,0        ! cylindrical c.s. center point
!               ! for creating fillet line
csn=csn+1           ! creating cylindrical c.s.
cskp,csn,1,km+1,km-3,km-2,,
!
l,km-2,km-3         ! fillet line
! 
a,km-3,km-5,km,km-1,km-2    ! creating areas
a,km-6,km-5,km-3,km-4
!
*get,lm,line,,num,max       ! sizing elements on lines
!
lesize,lm,,,nefill,1,
lesize,lm-1,,,netop,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,nefill/2,1,
lesize,lm-4,,,nefill/2,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,netop,1,
lesize,lm-7,,,netop,1,
!
lccat,lm-3,lm-4         ! concatentation of lines
!
*get,am,area,,num,max
!
aatt,1,,2,0         ! set element attributes                
mshkey,1            ! select mapped meshing
amesh,am-1,am,1         ! mesh areas
!
*get,lm,line,,num,max       ! delete concatenated lines
ldele,lm
!
ESIZE,0,neL1,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,1               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,L1,,          ! extrude areas
VOFFST,am-1,L1,,
!
aclear,am-1,am,1        ! clear area meshes
!
csys,rib1cs         ! mirror volumes to create complete rib
*get,vm,volu,,num,max
vsymm,x,vm-1,vm,,,0,0
!
/view,1,1,1,1
/auto,1
/replot
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
! 
/view,1,1,1,1
/auto,1
lplot
!
/com, ********************** End First (Mid) Rib Layout ***************
!
/com, *************************** 2nd Rib Layout **********************
!
csys,0
!
k,,r2p1x,r2p1y,r2p1z        ! 3 keypoints forming planar surface 
k,,r2p2x,r2p2y,r2p2z        ! to build rib cross-sections 
k,,r2p3x,r2p3y,r2p3z
!
*get,km,kp,,num,max     ! creation of c.s. for rib in planar
csn=csn+1           ! surface
rib2cs=csn
cskp,rib2cs,0,km-2,km-1,km,,
!
kbetw,km-2,km-1,km+1,ratio,0.5  ! center keypoint of rib cross-section
!
cskp,rib2cs,0,km+1,km-1,km  ! move planar surface c.s. to rib center
!
k,,0,ribh-r1,0
k,,w1/2,0,0
k,,w1/2,ribh-r1,0       ! remaining rib cross-section keypoints 
k,,w1/2+r1,ribh,0
k,,w1/2+r1,ribh+topthk,0
k,,0,ribh+topthk,0
!
*get,km,kp,,num,max     ! form straight lines of rib
!
l,km-6,km-4
l,km-4,km-3
l,km-2,km-1
l,km-1,km
l,km,km-5
l,km-5,km-6
l,km-5,km-3
!
k,,w1/2+r1,ribh-r1,0        ! cylindrical c.s. center point
!               ! for creating fillet line
csn=csn+1           ! creating cylindrical c.s.
cskp,csn,1,km+1,km-3,km-2,,
!
l,km-2,km-3         ! fillet line
! 
a,km-3,km-5,km,km-1,km-2    ! creating areas
a,km-6,km-5,km-3,km-4
!
*get,lm,line,,num,max       ! sizing elements on lines
!
lesize,lm,,,nefill,1,
lesize,lm-1,,,netop,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,nefill/2,1,
lesize,lm-4,,,nefill/2,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,netop,1,
lesize,lm-7,,,netop,1,
!
lccat,lm-3,lm-4         ! concatentation of lines
!
*get,am,area,,num,max
!
aatt,1,,2,0         ! set element attributes                
mshkey,1            ! select mapped meshing
amesh,am-1,am,1         ! mesh areas
!
*get,lm,line,,num,max       ! delete concatenated lines
ldele,lm
!
ESIZE,0,neL2,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,1               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,L2,,          ! extrude areas
VOFFST,am-1,L2,,
!
aclear,am-1,am,1        ! clear area meshes
!
csys,rib2cs         ! mirror volumes to create complete rib
*get,vm,volu,,num,max
vsymm,x,vm-1,vm,,,0,0
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
! 
LPLOT
!
/com, ***************** End 2nd Rib Layout  ********************
!
/com, ***************** 3rd Rib Layout *************************
!
csys,0
!  
k,,r3p1x,r3p1y,r3p1z        ! 3 keypoints forming planar surface 
k,,r3p2x,r3p2y,r3p2z        ! to build rib cross-sections 
k,,r3p3x,r3p3y,r3p3z
!
*get,km,kp,,num,max     ! creation of c.s. for rib in planar
csn=csn+1           ! surface
rib3cs=csn
cskp,rib3cs,0,km-2,km-1,km,,
!
kbetw,km-2,km-1,km+1,ratio,0.5  ! center keypoint of rib cross-section
!
cskp,rib3cs,0,km+1,km-1,km  ! move planar surface c.s. to rib center
!
k,,0,ribh-r1,0
k,,w2/2,0,0
k,,w2/2,ribh-r1,0       ! remaining rib cross-section keypoints 
k,,w2/2+r1,ribh,0
k,,w2/2+r1,ribh+topthk,0
k,,0,ribh+topthk,0
!
*get,km,kp,,num,max     ! form straight lines of rib
!
l,km-6,km-4
l,km-4,km-3
l,km-2,km-1
l,km-1,km
l,km,km-5
l,km-5,km-6
l,km-5,km-3
!
k,,w2/2+r1,ribh-r1,0        ! cylindrical c.s. center point
!               ! for creating fillet line
csn=csn+1           ! creating cylindrical c.s.
cskp,csn,1,km+1,km-3,km-2,,
!
l,km-2,km-3         ! fillet line
! 
a,km-3,km-5,km,km-1,km-2    ! creating areas
a,km-6,km-5,km-3,km-4
!
*get,lm,line,,num,max       ! sizing elements on lines
!
lesize,lm,,,nefill,1,
lesize,lm-1,,,netop,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,nefill/2,1,
lesize,lm-4,,,nefill/2,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,netop,1,
lesize,lm-7,,,netop,1,
!
lccat,lm-3,lm-4         ! concatentation of lines
!
*get,am,area,,num,max
!
aatt,1,,2,0         ! set element attributes                
mshkey,1            ! select mapped meshing
amesh,am-1,am,1         ! mesh areas
!
*get,lm,line,,num,max       ! delete concatenated lines
ldele,lm
!
ESIZE,0,neL3,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,1               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,L3,,          ! extrude areas
VOFFST,am-1,L3,,
!
aclear,am-1,am,1        ! clear area meshes
!
csys,rib3cs         ! mirror volumes to create complete rib
*get,vm,volu,,num,max
vsymm,x,vm-1,vm,,,0,0
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
! 
LPLOT
!
/com, ***************** End 3rd Rib Layout *************************
!
/com, ******************* 4th Rib Layout ***************************
!
csys,0
!
k,,r4p1x,r4p1y,r4p1z        ! 3 keypoints forming planar surface 
k,,r4p2x,r4p2y,r4p2z        ! to build rib cross-sections 
k,,r4p3x,r4p3y,r4p3z
!
*get,km,kp,,num,max     ! creation of rib local cartesian c.s.
csn=csn+1
rib4cs=csn
cskp,rib4cs,0,km-2,km-1,km,,
!
kbetw,km-2,km-1,km+1,ratio,0.5  ! center keypoint of rib cross-section
!
cskp,rib4cs,0,km+1,km-1,km  ! move planar surface c.s. to rib center
!
k,,0,ribh-r1,0          ! creation of keypoints   
k,,0,fribh-r2,0
k,,-w1/2,0,0
k,,-w1/2,ribh-r1,0
k,,-(w1/2+r1),ribh,0
!   
k,,-(w1/2+r1),ribh+topthk,0
k,,0,ribh+topthk,0
k,,w1/2+r2,ribh+topthk,0
k,,w1/2+r2,fribh,0
k,,w1/2,fribh-r2,0
!   
k,,w1/2,0,0
k,,-(w1/2+r1),ribh-r1,0
k,,w1/2+r2,fribh-r2,0
!
*get,km,kp,,num,max     ! creation of straight lines
!  
l,km-10,km-9
l,km-8,km-7
l,km-7,km-6
l,km-6,km-12
l,km-12,km-11
!   
l,km-11,km-13
l,km-13,km-10
l,km-6,km-5
l,km-5,km-4
l,km-11,km-3
!   
l,km-3,km-2
l,km-2,km-13
l,km-9,km-12
! 
csn=csn+1           ! creation of rear fillet c.s.
cskp,csn,1,km-1,km-9,km-8,,
!
l,km-8,km-9         ! creation of rear fillet line
!
csn=csn+1           ! creation of front fillet c.s.
cskp,csn,1,km,km-3,km-4
!   
l,km-3,km-4         ! creation of front fillet line
!
a,km-9,km-8,km-7,km-6,km-12 ! creation of areas 
a,km-10,km-9,km-12,km-11,km-13
a,km-13,km-11,km-3,km-2
a,km-11,km-12,km-6,km-5,km-4,km-3
!
*get,lm,line,,num,max       ! sizing line elements
!
lesize,lm,,,neffill,1,
lesize,lm-1,,,nefill,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,netop,1,
lesize,lm-4,,,netop,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,neftop,1,
lesize,lm-7,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-9,,,netop,1,
!
lesize,lm-10,,,1,1,
lesize,lm-11,,,nefill/2,1,
lesize,lm-12,,,nefill/2,1
lesize,lm-13,,,netop,1,
lesize,lm-14,,,netop+1,1,
!
lccat,lm-12,lm-11       ! concatenate selected lines
lccat,lm-10,lm-9
!
*get,am,area,,num,max
!
aatt,1,,2,0         ! set area attributes
!
mshkey,1            ! select mapped meshing
amesh,am-3,am-1,1       ! map mesh selected areas
!
mshkey,0            ! select free meshing
amesh,am,am,            ! free mesh selected areas
!
ldele,lm+1          ! delete concatenated lines
ldele,lm+2
!
ESIZE,0,neL4,           ! # element divisions along extrusion
!
type,1              ! select element type for extrusion
mat,1               ! select material #
esys,0              ! select element c.s.
!
VOFFST,am,L4,,          ! extrude areas
VOFFST,am-1,L4,,
VOFFST,am-2,L4,,
VOFFST,am-3,L4,,
!
ESIZE,0,1,
!
VOFFST,am,-.1423,,          ! extrude areas
VOFFST,am-1,-.1423,,
!
aclear,am-3,am,1        ! clear area meshes
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
LPLOT
!
!
/com, ************************ End 4th Rib Layout *********************
!
!
/com, ************************ 5th Rib Layout *************************
!
csys,0
!  
k,,r5p1x,r5p1y,r5p1z        ! 3 keypoints forming planar surface 
k,,r5p2x,r5p2y,r5p2z        ! to build rib cross-sections 
k,,r5p3x,r5p3y,r5p3z
!
*get,km,kp,,num,max     ! creation of rib local cartesian c.s.
csn=csn+1
rib5cs=csn
cskp,rib5cs,0,km-2,km-1,km,,
!
kbetw,km-2,km-1,km+1,ratio,0.5  ! center keypoint of rib cross-section
!
cskp,rib5cs,0,km+1,km-1,km  ! move planar surface c.s. to rib center
!
k,,0,ribh-r1,0          ! creation of keypoints   
k,,0,fribh-r2,0
k,,-w1/2,0,0
k,,-w1/2,ribh-r1,0
k,,-(w1/2+r1),ribh,0
!   
k,,-(w1/2+r1),ribh+topthk,0
k,,0,ribh+topthk,0
k,,w1/2+r2,ribh+topthk,0
k,,w1/2+r2,fribh,0
k,,w1/2,fribh-r2,0
!   
k,,w1/2,0,0
k,,-(w1/2+r1),ribh-r1,0
k,,w1/2+r2,fribh-r2,0
!
*get,km,kp,,num,max     ! creation of straight lines
!  
l,km-10,km-9
l,km-8,km-7
l,km-7,km-6
l,km-6,km-12
l,km-12,km-11
!   
l,km-11,km-13
l,km-13,km-10
l,km-6,km-5
l,km-5,km-4
l,km-11,km-3
!   
l,km-3,km-2
l,km-2,km-13
l,km-9,km-12
! 
csn=csn+1           ! creation of rear fillet c.s.
cskp,csn,1,km-1,km-9,km-8,,
!
l,km-8,km-9         ! creation of rear fillet line
!
csn=csn+1           ! creation of front fillet c.s.
cskp,csn,1,km,km-3,km-4
!   
l,km-3,km-4         ! creation of front fillet line
!
a,km-9,km-8,km-7,km-6,km-12 ! creation of areas 
a,km-10,km-9,km-12,km-11,km-13
a,km-13,km-11,km-3,km-2
a,km-11,km-12,km-6,km-5,km-4,km-3
!
*get,lm,line,,num,max       ! sizing line elements
!
lesize,lm,,,neffill,1,
lesize,lm-1,,,nefill,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,netop,1,
lesize,lm-4,,,netop,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,neftop,1,
lesize,lm-7,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-9,,,netop,1,
!
lesize,lm-10,,,1,1,
lesize,lm-11,,,nefill/2,1,
lesize,lm-12,,,nefill/2,1
lesize,lm-13,,,netop,1,
lesize,lm-14,,,netop+1,1,
!
lccat,lm-12,lm-11       ! concatenate selected lines
lccat,lm-10,lm-9
!
*get,am,area,,num,max
!
aatt,1,,2,0         ! set area attributes
!
mshkey,1            ! select mapped meshing
amesh,am-3,am-1,1       ! map mesh selected areas
!
mshkey,0            ! select free meshing
amesh,am,am,            ! free mesh selected areas
!
ldele,lm+1          ! delete concatenated lines
ldele,lm+2
!
ESIZE,0,neL5,           ! # element divisions along extrusion
!
type,1              ! select element type for extrusion
mat,1               ! select material #
esys,0              ! select element c.s.
!
VOFFST,am,-L5,,         ! extrude areas
VOFFST,am-1,-L5,,
VOFFST,am-2,-L5,,
VOFFST,am-3,-L5,,
!
ESIZE,0,1,
!
VOFFST,am,.1423,,           ! extrude areas
VOFFST,am-1,.1423,,
!
aclear,am-3,am,1        ! clear area meshes
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
LPLOT
!
!
/com, ************************** End 5th Rib Layout ***********************
!
!
/com, ************************** 6th Rib Layout ***************************
!
csys,0
!  
k,,r6p1x,r6p1y,r6p1z        ! 3 keypoints forming planar surface 
k,,r6p2x,r6p2y,r6p2z        ! to build rib cross-sections 
k,,r6p3x,r6p3y,r6p3z
!
*get,km,kp,,num,max     ! creation of c.s. for rib in planar
csn=csn+1           ! surface
rib6cs=csn
cskp,rib6cs,0,km-2,km-1,km,,
!
kbetw,km-2,km-1,km+1,ratio,0.5  ! center keypoint of rib cross-section
!
cskp,rib6cs,0,km+1,km-1,km  ! move planar surface c.s. to rib center
!
k,,0,ribh-r1,0
k,,w2/2,0,0
k,,w2/2,ribh-r1,0       ! remaining rib cross-section keypoints 
k,,w2/2+r1,ribh,0
k,,w2/2+r1,ribh+topthk,0
k,,0,ribh+topthk,0
!
*get,km,kp,,num,max     ! form straight lines of rib
!
l,km-6,km-4
l,km-4,km-3
l,km-2,km-1
l,km-1,km
l,km,km-5
l,km-5,km-6
l,km-5,km-3
!
k,,w2/2+r1,ribh-r1,0        ! cylindrical c.s. center point
!               ! for creating fillet line
csn=csn+1           ! creating cylindrical c.s.
cskp,csn,1,km+1,km-3,km-2,,
!
l,km-2,km-3         ! fillet line
! 
a,km-3,km-5,km,km-1,km-2    ! creating areas
a,km-6,km-5,km-3,km-4
!
*get,lm,line,,num,max       ! sizing elements on lines
!
lesize,lm,,,nefill,1,
lesize,lm-1,,,netop,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,nefill/2,1,
lesize,lm-4,,,nefill/2,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,netop,1,
lesize,lm-7,,,netop,1,
!
lccat,lm-3,lm-4         ! concatentation of lines
!
*get,am,area,,num,max
!
aatt,1,,2,0         ! set element attributes                
mshkey,1            ! select mapped meshing
amesh,am-1,am,1         ! mesh areas
!
*get,lm,line,,num,max       ! delete concatenated lines
ldele,lm
!
ESIZE,0,neL6,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,1               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,L6,,          ! extrude areas
VOFFST,am-1,L6,,
!
aclear,am-1,am,1        ! clear area meshes
!
csys,rib6cs         ! mirror volumes to create complete rib
*get,vm,volu,,num,max
vsymm,x,vm-1,vm,,,0,0
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
! 
LPLOT
!
/com, ************************ End 6th Rib Layout ********************
!
!
/com ************************** 7th Rib Layout ***********************
!
csys,0
!
k,,r7p1x,r7p1y,r7p1z        ! 3 keypoints forming planar surface 
k,,r7p2x,r7p2y,r7p2z        ! to build rib cross-sections 
k,,r7p3x,r7p3y,r7p3z
!
*get,km,kp,,num,max     ! creation of c.s. for rib in planar
csn=csn+1           ! surface
rib7cs=csn
cskp,rib7cs,0,km-2,km-1,km,,
!
kbetw,km-2,km-1,km+1,ratio,0.5  ! center keypoint of rib cross-section
!
cskp,rib7cs,0,km+1,km-1,km  ! move planar surface c.s. to rib center
!
k,,0,ribh-r1,0
k,,w1/2,0,0
k,,w1/2,ribh-r1,0       ! remaining rib cross-section keypoints 
k,,w1/2+r1,ribh,0
k,,w1/2+r1,ribh+topthk,0
k,,0,ribh+topthk,0
!
*get,km,kp,,num,max     ! form straight lines of rib
!
l,km-6,km-4
l,km-4,km-3
l,km-2,km-1
l,km-1,km
l,km,km-5
l,km-5,km-6
l,km-5,km-3
!
k,,w1/2+r1,ribh-r1,0        ! cylindrical c.s. center point
!               ! for creating fillet line
csn=csn+1           ! creating cylindrical c.s.
cskp,csn,1,km+1,km-3,km-2,,
!
l,km-2,km-3         ! fillet line
! 
a,km-3,km-5,km,km-1,km-2    ! creating areas
a,km-6,km-5,km-3,km-4
!
*get,lm,line,,num,max       ! sizing elements on lines
!
lesize,lm,,,nefill,1,
lesize,lm-1,,,netop,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,nefill/2,1,
lesize,lm-4,,,nefill/2,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,netop,1,
lesize,lm-7,,,netop,1,
!
lccat,lm-3,lm-4         ! concatentation of lines
!
*get,am,area,,num,max
!
aatt,1,,2,0         ! set element attributes                
mshkey,1            ! select mapped meshing
amesh,am-1,am,1         ! mesh areas
!
*get,lm,line,,num,max       ! delete concatenated lines
ldele,lm
!
ESIZE,0,neL7,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,1               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,L7,,          ! extrude areas
VOFFST,am-1,L7,,
!
aclear,am-1,am,1        ! clear area meshes
!
csys,rib7cs         ! mirror volumes to create complete rib
*get,vm,volu,,num,max
vsymm,x,vm-1,vm,,,0,0
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
! 
LPLOT
!
/com, ********************** End 7th Rib Layout  ********************
!
!
/com, ***************** Begin Center Boss Construction **************
!
/com  ******************** Mid (First) Rib, AREA 1 ******************
!
!
csys,rib1cs         ! rib 1 coordinate system
!
p1=kp(-(w1/2+r1),ribh+topthk,0) ! keypoints from position on rib 
p2=kp(-(w1/2+r1),ribh,0)
p3=kp(-w1/2,ribh-r1,0)
p4=kp(-w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta1=ky(p1)
!
k,,rcb1,theta1,0
k,,rcb1,theta1,-cbore1
k,,rcb2,theta1,-cbore1
k,,rcb2,theta1,-(topthk+r1)
k,,rcb2,theta1,-cbore2
! 
k,,rbore,theta1,-cbore2
k,,rbore,theta1,-(ribh+topthk)
k,,rcb2,theta1,-(ribh+topthk)
k,,rcb1,theta1,-(ribh+topthk)
k,,rbout,theta1,-(ribh+topthk)
!   
k,,rbout,theta1,-(topthk+r1)
k,,p1r,theta1,-topthk
k,,p1r,theta1,0
k,,rcb1,theta1,-(topthk+r1)
k,,rcb1,theta1,-cbore2
!
k,,p1r,theta1,-(topthk+r1)  ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib1cs         ! rib coordinate system
!
sphrp1=kp(-w1/2,ribh-r1,0)  ! derive spherical surface keypoints
sphrp2=kp(-(w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* Mid (First) Rib, AREA 2 ********************
!
!
csys,rib1cs         ! rib 1 coordinate system
!
p1=kp((w1/2+r1),ribh+topthk,0)  ! keypoints from position on rib 
p2=kp((w1/2+r1),ribh,0)
p3=kp(w1/2,ribh-r1,0)
p4=kp(w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta2=ky(p1)
!
k,,rcb1,theta2,0
k,,rcb1,theta2,-cbore1
k,,rcb2,theta2,-cbore1
k,,rcb2,theta2,-(topthk+r1)
k,,rcb2,theta2,-cbore2
! 
k,,rbore,theta2,-cbore2
k,,rbore,theta2,-(ribh+topthk)
k,,rcb2,theta2,-(ribh+topthk)
k,,rcb1,theta2,-(ribh+topthk)
k,,rbout,theta2,-(ribh+topthk)
!   
k,,rbout,theta2,-(topthk+r1)
k,,p1r,theta2,-topthk
k,,p1r,theta2,0
k,,rcb1,theta2,-(topthk+r1)
k,,rcb1,theta2,-cbore2
!
k,,p1r,theta2,-(topthk+r1)  ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib1cs         ! rib coordinate system
!
sphrp1=kp(w1/2,ribh-r1,0)   ! derive spherical surface keypoints
sphrp2=kp((w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ******************** 2nd Rib, AREA 3 ******************
!
!
csys,rib2cs         ! rib coordinate system
!
p1=kp(-(w1/2+r1),ribh+topthk,0) ! keypoints from position on rib 
p2=kp(-(w1/2+r1),ribh,0)
p3=kp(-w1/2,ribh-r1,0)
p4=kp(-w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta3=ky(p1)
!
k,,rcb1,theta3,0
k,,rcb1,theta3,-cbore1
k,,rcb2,theta3,-cbore1
k,,rcb2,theta3,-(topthk+r1)
k,,rcb2,theta3,-cbore2
! 
k,,rbore,theta3,-cbore2
k,,rbore,theta3,-(ribh+topthk)
k,,rcb2,theta3,-(ribh+topthk)
k,,rcb1,theta3,-(ribh+topthk)
k,,rbout,theta3,-(ribh+topthk)
!   
k,,rbout,theta3,-(topthk+r1)
k,,p1r,theta3,-topthk
k,,p1r,theta3,0
k,,rcb1,theta3,-(topthk+r1)
k,,rcb1,theta3,-cbore2
!
k,,p1r,theta3,-(topthk+r1)  ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib2cs         ! rib coordinate system
!
sphrp1=kp(-w1/2,ribh-r1,0)  ! derive spherical surface keypoints
sphrp2=kp(-(w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 2nd Rib, AREA 4 ********************
!
!
csys,rib2cs         ! rib 1 coordinate system
!
p1=kp((w1/2+r1),ribh+topthk,0)  ! keypoints from position on rib 
p2=kp((w1/2+r1),ribh,0)
p3=kp(w1/2,ribh-r1,0)
p4=kp(w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta4=ky(p1)
!
k,,rcb1,theta4,0
k,,rcb1,theta4,-cbore1
k,,rcb2,theta4,-cbore1
k,,rcb2,theta4,-(topthk+r1)
k,,rcb2,theta4,-cbore2
! 
k,,rbore,theta4,-cbore2
k,,rbore,theta4,-(ribh+topthk)
k,,rcb2,theta4,-(ribh+topthk)
k,,rcb1,theta4,-(ribh+topthk)
k,,rbout,theta4,-(ribh+topthk)
!   
k,,rbout,theta4,-(topthk+r1)
k,,p1r,theta4,-topthk
k,,p1r,theta4,0
k,,rcb1,theta4,-(topthk+r1)
k,,rcb1,theta4,-cbore2
!
k,,p1r,theta4,-(topthk+r1)  ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib2cs         ! rib coordinate system
!
sphrp1=kp(w1/2,ribh-r1,0)   ! derive spherical surface keypoints
sphrp2=kp((w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 3rd Rib, AREA 5 ********************
!
!
csys,rib3cs         ! rib coordinate system
!
p1=kp(-(w2/2+r1),ribh+topthk,0) ! keypoints from position on rib 
p2=kp(-(w2/2+r1),ribh,0)
p3=kp(-w2/2,ribh-r1,0)
p4=kp(-w2/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta5=ky(p1)
!
k,,rcb1,theta5,0
k,,rcb1,theta5,-cbore1
k,,rcb2,theta5,-cbore1
k,,rcb2,theta5,-(topthk+r1)
k,,rcb2,theta5,-cbore2
! 
k,,rbore,theta5,-cbore2
k,,rbore,theta5,-(ribh+topthk)
k,,rcb2,theta5,-(ribh+topthk)
k,,rcb1,theta5,-(ribh+topthk)
k,,rbout,theta5,-(ribh+topthk)
!   
k,,rbout,theta5,-(topthk+r1)
k,,p1r,theta5,-topthk
k,,p1r,theta5,0
k,,rcb1,theta5,-(topthk+r1)
k,,rcb1,theta5,-cbore2
!
k,,p1r,theta5,-(topthk+r1)  ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib3cs         ! rib coordinate system
!
sphrp1=kp(-w2/2,ribh-r1,0)  ! derive spherical surface keypoints
sphrp2=kp(-(w2/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 3rd Rib, AREA 6 ********************
!
!
csys,rib3cs         ! rib 1 coordinate system
!
p1=kp((w2/2+r1),ribh+topthk,0)  ! keypoints from position on rib 
p2=kp((w2/2+r1),ribh,0)
p3=kp(w2/2,ribh-r1,0)
p4=kp(w2/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta6=ky(p1)
!
k,,rcb1,theta6,0
k,,rcb1,theta6,-cbore1
k,,rcb2,theta6,-cbore1
k,,rcb2,theta6,-(topthk+r1)
k,,rcb2,theta6,-cbore2
! 
k,,rbore,theta6,-cbore2
k,,rbore,theta6,-(ribh+topthk)
k,,rcb2,theta6,-(ribh+topthk)
k,,rcb1,theta6,-(ribh+topthk)
k,,rbout,theta6,-(ribh+topthk)
!   
k,,rbout,theta6,-(topthk+r1)
k,,p1r,theta6,-topthk
k,,p1r,theta6,0
k,,rcb1,theta6,-(topthk+r1)
k,,rcb1,theta6,-cbore2
!
k,,p1r,theta6,-(topthk+r1)  ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib3cs         ! rib coordinate system
!
sphrp1=kp(w2/2,ribh-r1,0)   ! derive spherical surface keypoints
sphrp2=kp((w2/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 4th Rib, AREA 7 ********************
!
!
csys,rib4cs         ! rib coordinate system
!
p1=kp(-(w1/2+r1),ribh+topthk,0) ! keypoints from position on rib 
p2=kp(-(w1/2+r1),ribh,0)
p3=kp(-w1/2,ribh-r1,0)
p4=kp(-w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta7=ky(p1)
!
k,,rcb1,theta7,0
k,,rcb1,theta7,-cbore1
k,,rcb2,theta7,-cbore1
k,,rcb2,theta7,-(topthk+r1)
k,,rcb2,theta7,-cbore2
! 
k,,rbore,theta7,-cbore2
k,,rbore,theta7,-(ribh+topthk)
k,,rcb2,theta7,-(ribh+topthk)
k,,rcb1,theta7,-(ribh+topthk)
k,,rbout,theta7,-(ribh+topthk)
!   
k,,rbout,theta7,-(topthk+r1)
k,,p1r,theta7,-topthk
k,,p1r,theta7,0
k,,rcb1,theta7,-(topthk+r1)
k,,rcb1,theta7,-cbore2
!
k,,p1r,theta7,-(topthk+r1)  ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib4cs         ! rib coordinate system
!
sphrp1=kp(-w1/2,ribh-r1,0)  ! derive spherical surface keypoints
sphrp2=kp(-(w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 5th Rib, AREA 10 ********************
!
!
csys,rib5cs         ! rib coordinate system
!
p1=kp(-(w1/2+r1),ribh+topthk,0) ! keypoints from position on rib 
p2=kp(-(w1/2+r1),ribh,0)
p3=kp(-w1/2,ribh-r1,0)
p4=kp(-w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta10=ky(p1)
!
k,,rcb1,theta10,0
k,,rcb1,theta10,-cbore1
k,,rcb2,theta10,-cbore1
k,,rcb2,theta10,-(topthk+r1)
k,,rcb2,theta10,-cbore2
! 
k,,rbore,theta10,-cbore2
k,,rbore,theta10,-(ribh+topthk)
k,,rcb2,theta10,-(ribh+topthk)
k,,rcb1,theta10,-(ribh+topthk)
k,,rbout,theta10,-(ribh+topthk)
!   
k,,rbout,theta10,-(topthk+r1)
k,,p1r,theta10,-topthk
k,,p1r,theta10,0
k,,rcb1,theta10,-(topthk+r1)
k,,rcb1,theta10,-cbore2
!
k,,p1r,theta10,-(topthk+r1) ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib5cs         ! rib coordinate system
!
sphrp1=kp(-w1/2,ribh-r1,0)  ! derive spherical surface keypoints
sphrp2=kp(-(w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ******************** 6th Rib, AREA 11 ******************
!
!
csys,rib6cs         ! rib 1 coordinate system
!
p1=kp(-(w2/2+r1),ribh+topthk,0) ! keypoints from position on rib 
p2=kp(-(w2/2+r1),ribh,0)
p3=kp(-w2/2,ribh-r1,0)
p4=kp(-w2/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta11=ky(p1)
!
k,,rcb1,theta11,0
k,,rcb1,theta11,-cbore1
k,,rcb2,theta11,-cbore1
k,,rcb2,theta11,-(topthk+r1)
k,,rcb2,theta11,-cbore2
! 
k,,rbore,theta11,-cbore2
k,,rbore,theta11,-(ribh+topthk)
k,,rcb2,theta11,-(ribh+topthk)
k,,rcb1,theta11,-(ribh+topthk)
k,,rbout,theta11,-(ribh+topthk)
!   
k,,rbout,theta11,-(topthk+r1)
k,,p1r,theta11,-topthk
k,,p1r,theta11,0
k,,rcb1,theta11,-(topthk+r1)
k,,rcb1,theta11,-cbore2
!
k,,p1r,theta11,-(topthk+r1) ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib6cs         ! rib coordinate system
!
sphrp1=kp(-w2/2,ribh-r1,0)  ! derive spherical surface keypoints
sphrp2=kp(-(w2/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 6th Rib, AREA 12 ********************
!
!
csys,rib6cs         ! rib 1 coordinate system
!
p1=kp((w2/2+r1),ribh+topthk,0)  ! keypoints from position on rib 
p2=kp((w2/2+r1),ribh,0)
p3=kp(w2/2,ribh-r1,0)
p4=kp(w2/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta12=ky(p1)
!
k,,rcb1,theta12,0
k,,rcb1,theta12,-cbore1
k,,rcb2,theta12,-cbore1
k,,rcb2,theta12,-(topthk+r1)
k,,rcb2,theta12,-cbore2
! 
k,,rbore,theta12,-cbore2
k,,rbore,theta12,-(ribh+topthk)
k,,rcb2,theta12,-(ribh+topthk)
k,,rcb1,theta12,-(ribh+topthk)
k,,rbout,theta12,-(ribh+topthk)
!   
k,,rbout,theta12,-(topthk+r1)
k,,p1r,theta12,-topthk
k,,p1r,theta12,0
k,,rcb1,theta12,-(topthk+r1)
k,,rcb1,theta12,-cbore2
!
k,,p1r,theta12,-(topthk+r1) ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib6cs         ! rib coordinate system
!
sphrp1=kp(w2/2,ribh-r1,0)   ! derive spherical surface keypoints
sphrp2=kp((w2/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  *********************** 7th Rib, AREA 13 *********************
!
!
csys,rib7cs         ! rib coordinate system
!
p1=kp(-(w1/2+r1),ribh+topthk,0) ! keypoints from position on rib 
p2=kp(-(w1/2+r1),ribh,0)
p3=kp(-w1/2,ribh-r1,0)
p4=kp(-w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta13=ky(p1)
!
k,,rcb1,theta13,0
k,,rcb1,theta13,-cbore1
k,,rcb2,theta13,-cbore1
k,,rcb2,theta13,-(topthk+r1)
k,,rcb2,theta13,-cbore2
! 
k,,rbore,theta13,-cbore2
k,,rbore,theta13,-(ribh+topthk)
k,,rcb2,theta13,-(ribh+topthk)
k,,rcb1,theta13,-(ribh+topthk)
k,,rbout,theta13,-(ribh+topthk)
!   
k,,rbout,theta13,-(topthk+r1)
k,,p1r,theta13,-topthk
k,,p1r,theta13,0
k,,rcb1,theta13,-(topthk+r1)
k,,rcb1,theta13,-cbore2
!
k,,p1r,theta13,-(topthk+r1) ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib7cs         ! rib coordinate system
!
sphrp1=kp(-w1/2,ribh-r1,0)  ! derive spherical surface keypoints
sphrp2=kp(-(w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 7th Rib, AREA 14 ********************
!
!
csys,rib7cs         ! rib 1 coordinate system
!
p1=kp((w1/2+r1),ribh+topthk,0)  ! keypoints from position on rib 
p2=kp((w1/2+r1),ribh,0)
p3=kp(w1/2,ribh-r1,0)
p4=kp(w1/2,0,0)
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta14=ky(p1)
!
k,,rcb1,theta14,0
k,,rcb1,theta14,-cbore1
k,,rcb2,theta14,-cbore1
k,,rcb2,theta14,-(topthk+r1)
k,,rcb2,theta14,-cbore2
! 
k,,rbore,theta14,-cbore2
k,,rbore,theta14,-(ribh+topthk)
k,,rcb2,theta14,-(ribh+topthk)
k,,rcb1,theta14,-(ribh+topthk)
k,,rbout,theta14,-(ribh+topthk)
!   
k,,rbout,theta14,-(topthk+r1)
k,,p1r,theta14,-topthk
k,,p1r,theta14,0
k,,rcb1,theta14,-(topthk+r1)
k,,rcb1,theta14,-cbore2
!
k,,p1r,theta14,-(topthk+r1) ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn           ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!       
k,,r1,45,0          ! creation of fillet keypoint
!
*get,km,kp,,num,max     ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-12,km-2
lstr,km-8,km-2
lstr,km-9,km-12
lstr,km-16,km
!
l,km-6,km           ! creation of fillet line 1
l,km,km-5           ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9    ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-8,km-7,km-6,km-3,km-2
!       
a,km-3,km-6,km,km-16,km-15
a,km,km-5,km-4,km-16
!
csn=csn+1
CS1=csn         ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn            ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,netop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,netop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,netop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,netop,1,
!
lesize,lm-1,,,3,1,
lesize,lm,,,2,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-6           ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4       ! creation of cylindrical area
!
lesize,lm+1,,,4,1,      ! sizing of cylindrical lines
lesize,lm+2,,,4,1,
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib7cs         ! rib coordinate system
!
sphrp1=kp(w1/2,ribh-r1,0)   ! derive spherical surface keypoints
sphrp2=kp((w1/2+r1),ribh,0)
!
csys,plncs          ! planar cylindrical cs
!
sphrp3=kp(r1,45,0)      ! derive spherical surface keypoints
sphrp4=kp(r1,0,0)
!
csys,csn            ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4   ! creation of spherical area
!
nummrg,all,,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
/com  ********************* 4th Rib, AREA 8 ********************
!
!
csys,rib4cs             ! rib 1 coordinate system
!
p1=kp((w1/2+r2),fribh+ftopthk,.1423)    ! keypoints from position on rib 
p2=kp((w1/2+r2),fribh,.1423)
p3=kp(w1/2,fribh-r2,.1423)
p4=kp(w1/2,0,.1423)
!
csys,gcs                ! global coordinate system
!
p1r=kx(p1)
theta8=ky(p1)
!
k,,rcb1,theta8,0
k,,rcb1,theta8,-cbore1
k,,rcb2,theta8,-cbore1
k,,rcb2,theta8,-ftopthk
k,,rcb2,theta8,-cbore2
! 
k,,rbore,theta8,-cbore2
k,,rbore,theta8,-(ribh+topthk)
k,,rcb2,theta8,-(ribh+topthk)
k,,rcb1,theta8,-(ribh+topthk)
k,,p1r-r2,theta8,-(ribh+topthk)
!   
k,,p1r-r2,theta8,-(ftopthk+r2)
k,,p1r,theta8,-ftopthk
k,,p1r,theta8,0
k,,rcb1,theta8,-ftopthk
k,,rcb1,theta8,-cbore2
!
k,,p1r,theta8,-(ftopthk+r2)     ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max         ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn               ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!
k,,r2,45,0
!       
*get,km,kp,,num,max         ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-2,km-12
lstr,km-8,km-2
lstr,km-12,km-9
lstr,km,km-15
!
l,km-6,km               ! creation of fillet line 1
l,km,km-5               ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9        ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-3,km-2,km-8,km-7,km-6
!       
a,km-15,km-3,km-6,km
a,km,km-5,km-4,km-16,km-15
!
csn=csn+1
CS1=csn             ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn                ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max           ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,neftop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,neftop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,neftop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,neftop,1,
!
lesize,lm-1,,,2,1,
lesize,lm,,,3,1,
!
csys,CS1                ! cylindrical cs at center point
!
l,p3,km-6               ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4           ! creation of cylindrical area
!
lesize,lm+1,,,3,1,          ! sizing of cylindrical lines
lesize,lm+2,,,3,1,
!
nummrg,all,,.001            ! merge all coincident entities
numcmp,all              ! renumber entities
!
csys,rib4cs             ! rib coordinate system
!
sphrp1=kp(w1/2,fribh-r2,.1423)      ! derive spherical surface keypoints
sphrp2=kp((w1/2+r2),fribh,.1423)
!
csys,plncs              ! planar cylindrical cs
!
sphrp3=kp(r2,45,0)          ! derive spherical surface keypoints
sphrp4=kp(r2,0,0)
!
csys,csn                ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4       ! creation of spherical area
!
nummrg,all,,.001            ! merge all coincident entities
numcmp,all              ! renumber entities
!
LPLOT                   ! plot lines
!
!
!/com  ********************* 5th Rib, AREA 9 ********************
!
!
csys,rib5cs             ! rib 1 coordinate system
!
p1=kp((w1/2+r2),fribh+ftopthk,-.1423)   ! keypoints from position on rib 
p2=kp((w1/2+r2),fribh,-.1423)
p3=kp(w1/2,fribh-r2,-.1423)
p4=kp(w1/2,0,-.1423)
!
csys,gcs                ! global coordinate system
!
p1r=kx(p1)
theta9=ky(p1)
!
k,,rcb1,theta9,0
k,,rcb1,theta9,-cbore1
k,,rcb2,theta9,-cbore1
k,,rcb2,theta9,-ftopthk
k,,rcb2,theta9,-cbore2
! 
k,,rbore,theta9,-cbore2
k,,rbore,theta9,-(ribh+topthk)
k,,rcb2,theta9,-(ribh+topthk)
k,,rcb1,theta9,-(ribh+topthk)
k,,p1r-r2,theta9,-(ribh+topthk)
!   
k,,p1r-r2,theta9,-(ftopthk+r2)
k,,p1r,theta9,-ftopthk
k,,p1r,theta9,0
k,,rcb1,theta9,-ftopthk
k,,rcb1,theta9,-cbore2
!
k,,p1r,theta9,-(ftopthk+r2)     ! center point of rectangular & spherical cs
!
*get,km,kp,,num,max         ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn               ! for fillet keypoint and fillet lines
cskp,plncs,1,km,km-5,km-4
!
k,,r2,45,0
!       
*get,km,kp,,num,max         ! creation of all lines
!
lstr,km-4,km-16
lstr,km-16,km-15
lstr,km-15,km-14
lstr,km-14,km-13
lstr,km-13,km-12
!
lstr,km-12,km-11
lstr,km-11,km-10
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-8,km-7
!
lstr,km-7,km-6
lstr,km-5,km-4
lstr,km-15,km-3
lstr,km-6,km-3
lstr,km-13,km-3
!
lstr,km-2,km-3
lstr,km-2,km-12
lstr,km-8,km-2
lstr,km-12,km-9
lstr,km,km-15
!
l,km-6,km               ! creation of fillet line 1
l,km,km-5               ! creation of fillet line 2
!
a,km-12,km-11,km-10,km-9        ! creation of areas
a,km-12,km-9,km-8,km-2
a,km-12,km-2,km-3,km-13
a,km-13,km-3,km-15,km-14
a,km-3,km-2,km-8,km-7,km-6
!       
a,km-15,km-3,km-6,km
a,km,km-5,km-4,km-16,km-15
!
csn=csn+1
CS1=csn             ! creation of cylindrical cs
cskp,CS1,1,km-1,km-6,p3
!
csn=csn+1
CSr2=csn                ! creation of spherical cs
cskp,CSr2,2,km-1,km-6,km-5
!
*get,lm,line,,num,max           ! element sizing of lines
!
lesize,lm-21,,,2,1,
lesize,lm-20,,,1,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
!
lesize,lm-16,,,2,1,
lesize,lm-15,,,2,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,2,1,
lesize,lm-12,,,neftop,1,
!
lesize,lm-11,,,4,1,
lesize,lm-10,,,neftop,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,neftop,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
lesize,lm-3,,,2,1,
lesize,lm-2,,,neftop,1,
!
lesize,lm-1,,,2,1,
lesize,lm,,,3,1,
!
csys,CS1                ! cylindrical cs at center point
!
l,p3,km-6               ! creation of cylindrical lines
l,p4,km-7
!
a,p3,km-6,km-7,p4           ! creation of cylindrical area
!
lesize,lm+1,,,3,1,          ! sizing of cylindrical lines
lesize,lm+2,,,3,1,
!
nummrg,all,,.001            ! merge all coincident entities
numcmp,all              ! renumber entities
!
csys,rib5cs             ! rib coordinate system
!
sphrp1=kp(w1/2,fribh-r2,-.1423)     ! derive spherical surface keypoints
sphrp2=kp((w1/2+r2),fribh,-.1423)
!
csys,plncs              ! planar cylindrical cs
!
sphrp3=kp(r2,45,0)          ! derive spherical surface keypoints
sphrp4=kp(r2,0,0)
!
csys,csn                ! spherical coordinate system
!
a,sphrp1,sphrp2,sphrp3,sphrp4       ! creation of spherical area
!
nummrg,all,,.001            ! merge all coincident entities
numcmp,all              ! renumber entities
!
LPLOT                   ! plot lines
!
!
/com, ******************* Center Boss Volume Lines *********************
! 
!
/com, ************************* Areas 2 & 3 ****************************
!
rib1st=rib1cs               ! first rib defining volume
rib2nd=rib2cs               ! second rib defining volume
!
theta1st=theta2             ! first angle defining volume
theta2nd=theta3             ! second angle defining volume
! 
csys,rib1st
!
p11=kp((w1/2+r1),ribh+topthk,0)     ! define keypoints forming volume lines 
p12=kp((w1/2+r1),ribh,0)
!
csys,gcs
!
p11r=kx(p11)
A=p11r-r1*cos(45*rconv)
B=r1-r1*sin(45*rconv)
p13=kp(A,theta1st,-(topthk+B))
!
csys,rib2nd
!
p21=kp(-(w1/2+r1),ribh+topthk,0)
p22=kp(-(w1/2+r1),ribh,0)
!
csys,gcs
!
p21r=kx(p21)
C=p21r-r1*cos(45*rconv)
D=r1-r1*sin(45*rconv)
p23=kp(C,theta2nd,-(topthk+D))
!
p14=kp(rbout,theta1st,-(topthk+r1))
p15=kp(rbout,theta1st,-(topthk+ribh))
p16=kp(rcb1,theta1st,-(topthk+ribh))
p17=kp(rcb1,theta1st,-cbore2)
p18=kp(rcb1,theta1st,-(topthk+r1))
p19=kp(rcb1,theta1st,-cbore1)
p110=kp(rcb1,theta1st,0)
p111=kp(rcb2,theta1st,-(topthk+ribh))
p112=kp(rcb2,theta1st,-cbore2)
p113=kp(rcb2,theta1st,-(topthk+r1))
p114=kp(rcb2,theta1st,-cbore1)
p115=kp(rbore,theta1st,-(topthk+ribh))
p116=kp(rbore,theta1st,-cbore2)
!
p24=kp(rbout,theta2nd,-(topthk+r1))
p25=kp(rbout,theta2nd,-(topthk+ribh))
p26=kp(rcb1,theta2nd,-(topthk+ribh))
p27=kp(rcb1,theta2nd,-cbore2)
p28=kp(rcb1,theta2nd,-(topthk+r1))
p29=kp(rcb1,theta2nd,-cbore1)
p210=kp(rcb1,theta2nd,0)
p211=kp(rcb2,theta2nd,-(topthk+ribh))
p212=kp(rcb2,theta2nd,-cbore2)
p213=kp(rcb2,theta2nd,-(topthk+r1))
p214=kp(rcb2,theta2nd,-cbore1)
p215=kp(rbore,theta2nd,-(topthk+ribh))
p216=kp(rbore,theta2nd,-cbore2)
!
l,p11,p21               ! create lines to form volumes
l,p12,p22
l,p13,p23
l,p14,p24
l,p15,p25
!
l,p16,p26
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
l,p113,p213
l,p114,p214
l,p115,p215
!
l,p116,p216
!
*get,lm,line,,num,max
!
lsel,s,line,,lm-15,lm,1,0
lesize,all,,,6,1,
!
lsel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Create Areas <<<<<<<<<<<<<<<<<<<<<<<
!
csys,gcs
!
asel,s,loc,y,theta1st,theta1st,,
!
*get,am11,area,,num,min
!
am12=am11+1
am13=am11+2
am14=am11+3
am15=am11+4
am16=am11+5
am17=am11+6
!
asel,s,loc,y,theta2nd,theta2nd,,
!
*get,am21,area,,num,min
!
am22=am21+1
am23=am21+2
am24=am21+3
am25=am21+4
am26=am21+5
am27=am21+6
!
a,p11,p12,p22,p21
a,p11,p110,p210,p21
a,p13,p110,p210,p23
a,p12,p13,p23,p22
a,p13,p14,p24,p23
!
a,p14,p18,p28,p24
a,p19,p18,p28,p29
a,p110,p19,p29,p210
a,p14,p15,p25,p24
a,p15,p16,p26,p25
!
a,p17,p16,p26,p27
a,p18,p17,p27,p28
a,p16,p111,p211,p26
a,p112,p111,p211,p212
a,p17,p112,p212,p27
!
a,p113,p112,p212,p213
a,p18,p113,p213,p28
a,p114,p113,p213,p214
a,p19,p114,p214,p29
a,p111,p115,p215,p211
!
a,p116,p115,p215,p216
a,p112,p116,p216,p212
!
asel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<< Create Volumes <<<<<<<<<<<<<<<<<<<<<<<
!
*get,am,area,,num,max
!
va,am11,am21,am,am-1,am-2,am-8          ! create volumes
va,am12,am22,am-7,am-8,am-9,am-11
va,am13,am23,am-5,am-6,am-7,am-10
va,am14,am24,am-3,am-4,am-5,am-15
va,am15,am25,am-16,am-10,am-11,am-12,am-13
!
va,am16,am26,am-19,am-17,am-16,am-15,am-14
va,am17,am27,am-21,am-20,am-19,am-18
!
*get,vm,volu,,num,max
!
vsel,s,volu,,vm-6,vm,,
!
accat,am-15,am-14
accat,am-10,am-11
!
vatt,1,,1,0             ! set element attributes                
mshkey,1                ! select mapped meshing
vmesh,all               ! mesh areas
!
vsel,all
!
adele,am+1              ! delete concatenated areas
adele,am+2
!
*get,lm,line,,num,max
!
ldele,lm
ldele,lm-1
ldele,lm-2
ldele,lm-3
!
LPLOT
!
!
/com, ************************* Areas 4 & 5 ****************************
!
rib1st=rib2cs               ! first rib defining volume
rib2nd=rib3cs               ! second rib defining volume
!
theta1st=theta4             ! first angle defining volume
theta2nd=theta5             ! second angle defining volume
! 
csys,rib1st
!
p11=kp((w1/2+r1),ribh+topthk,0)     ! define keypoints forming volume lines 
p12=kp((w1/2+r1),ribh,0)
!
csys,gcs
!
p11r=kx(p11)
A=p11r-r1*cos(45*rconv)
B=r1-r1*sin(45*rconv)
p13=kp(A,theta1st,-(topthk+B))
!
csys,rib2nd
!
p21=kp(-(w1/2+r1),ribh+topthk,0)
p22=kp(-(w1/2+r1),ribh,0)
!
csys,gcs
!
p21r=kx(p21)
C=p21r-r1*cos(45*rconv)
D=r1-r1*sin(45*rconv)
p23=kp(C,theta2nd,-(topthk+D))
!
p14=kp(rbout,theta1st,-(topthk+r1))
p15=kp(rbout,theta1st,-(topthk+ribh))
p16=kp(rcb1,theta1st,-(topthk+ribh))
p17=kp(rcb1,theta1st,-cbore2)
p18=kp(rcb1,theta1st,-(topthk+r1))
p19=kp(rcb1,theta1st,-cbore1)
p110=kp(rcb1,theta1st,0)
p111=kp(rcb2,theta1st,-(topthk+ribh))
p112=kp(rcb2,theta1st,-cbore2)
p113=kp(rcb2,theta1st,-(topthk+r1))
p114=kp(rcb2,theta1st,-cbore1)
p115=kp(rbore,theta1st,-(topthk+ribh))
p116=kp(rbore,theta1st,-cbore2)
!
p24=kp(rbout,theta2nd,-(topthk+r1))
p25=kp(rbout,theta2nd,-(topthk+ribh))
p26=kp(rcb1,theta2nd,-(topthk+ribh))
p27=kp(rcb1,theta2nd,-cbore2)
p28=kp(rcb1,theta2nd,-(topthk+r1))
p29=kp(rcb1,theta2nd,-cbore1)
p210=kp(rcb1,theta2nd,0)
p211=kp(rcb2,theta2nd,-(topthk+ribh))
p212=kp(rcb2,theta2nd,-cbore2)
p213=kp(rcb2,theta2nd,-(topthk+r1))
p214=kp(rcb2,theta2nd,-cbore1)
p215=kp(rbore,theta2nd,-(topthk+ribh))
p216=kp(rbore,theta2nd,-cbore2)
!
l,p11,p21               ! create lines to form volumes
l,p12,p22
l,p13,p23
l,p14,p24
l,p15,p25
!
l,p16,p26
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
l,p113,p213
l,p114,p214
l,p115,p215
!
l,p116,p216
!
*get,lm,line,,num,max
!
lsel,s,line,,lm-15,lm,1,0
lesize,all,,,3,1,
!
lsel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Create Areas <<<<<<<<<<<<<<<<<<<<<<<
!
!
csys,gcs
!
asel,s,loc,y,theta1st,theta1st,,
!
*get,am11,area,,num,min
!
am12=am11+1
am13=am11+2
am14=am11+3
am15=am11+4
am16=am11+5
am17=am11+6
!
asel,s,loc,y,theta2nd,theta2nd,,
!
*get,am21,area,,num,min
!
am22=am21+1
am23=am21+2
am24=am21+3
am25=am21+4
am26=am21+5
am27=am21+6
!
a,p11,p12,p22,p21
a,p11,p110,p210,p21
a,p13,p110,p210,p23
a,p12,p13,p23,p22
a,p13,p14,p24,p23
!
a,p14,p18,p28,p24
a,p19,p18,p28,p29
a,p110,p19,p29,p210
a,p14,p15,p25,p24
a,p15,p16,p26,p25
!
a,p17,p16,p26,p27
a,p18,p17,p27,p28
a,p16,p111,p211,p26
a,p112,p111,p211,p212
a,p17,p112,p212,p27
!
a,p113,p112,p212,p213
a,p18,p113,p213,p28
a,p114,p113,p213,p214
a,p19,p114,p214,p29
a,p111,p115,p215,p211
!
a,p116,p115,p215,p216
a,p112,p116,p216,p212
!
asel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<< Create Volumes <<<<<<<<<<<<<<<<<<<<<<<
!
*get,am,area,,num,max
!
va,am11,am21,am,am-1,am-2,am-8          ! create volumes
va,am12,am22,am-7,am-8,am-9,am-11
va,am13,am23,am-5,am-6,am-7,am-10
va,am14,am24,am-3,am-4,am-5,am-15
va,am15,am25,am-16,am-10,am-11,am-12,am-13
!
va,am16,am26,am-19,am-17,am-16,am-15,am-14
va,am17,am27,am-21,am-20,am-19,am-18
!
*get,vm,volu,,num,max
!
vsel,s,volu,,vm-6,vm,,
!
accat,am-15,am-14
accat,am-10,am-11
!
vatt,1,,1,0             ! set element attributes                
mshkey,1                ! select mapped meshing
vmesh,all               ! mesh areas
!
vsel,all
!
adele,am+1              ! delete concatenated areas
adele,am+2
!
*get,lm,line,,num,max
!
ldele,lm
ldele,lm-1
ldele,lm-2
ldele,lm-3
!
LPLOT
!
!
/com, ************************* Areas 6 & 7 ****************************
!
rib1st=rib3cs               ! first rib defining volume
rib2nd=rib4cs               ! second rib defining volume
!
theta1st=theta6             ! first angle defining volume
theta2nd=theta7             ! second angle defining volume
! 
csys,rib1st
!
p11=kp((w1/2+r1),ribh+topthk,0)     ! define keypoints forming volume lines 
p12=kp((w1/2+r1),ribh,0)
!
csys,gcs
!
p11r=kx(p11)
A=p11r-r1*cos(45*rconv)
B=r1-r1*sin(45*rconv)
p13=kp(A,theta1st,-(topthk+B))
!
csys,rib2nd
!
p21=kp(-(w1/2+r1),ribh+topthk,0)
p22=kp(-(w1/2+r1),ribh,0)
!
csys,gcs
!
p21r=kx(p21)
C=p21r-r1*cos(45*rconv)
D=r1-r1*sin(45*rconv)
p23=kp(C,theta2nd,-(topthk+D))
!
p14=kp(rbout,theta1st,-(topthk+r1))
p15=kp(rbout,theta1st,-(topthk+ribh))
p16=kp(rcb1,theta1st,-(topthk+ribh))
p17=kp(rcb1,theta1st,-cbore2)
p18=kp(rcb1,theta1st,-(topthk+r1))
p19=kp(rcb1,theta1st,-cbore1)
p110=kp(rcb1,theta1st,0)
p111=kp(rcb2,theta1st,-(topthk+ribh))
p112=kp(rcb2,theta1st,-cbore2)
p113=kp(rcb2,theta1st,-(topthk+r1))
p114=kp(rcb2,theta1st,-cbore1)
p115=kp(rbore,theta1st,-(topthk+ribh))
p116=kp(rbore,theta1st,-cbore2)
!
p24=kp(rbout,theta2nd,-(topthk+r1))
p25=kp(rbout,theta2nd,-(topthk+ribh))
p26=kp(rcb1,theta2nd,-(topthk+ribh))
p27=kp(rcb1,theta2nd,-cbore2)
p28=kp(rcb1,theta2nd,-(topthk+r1))
p29=kp(rcb1,theta2nd,-cbore1)
p210=kp(rcb1,theta2nd,0)
p211=kp(rcb2,theta2nd,-(topthk+ribh))
p212=kp(rcb2,theta2nd,-cbore2)
p213=kp(rcb2,theta2nd,-(topthk+r1))
p214=kp(rcb2,theta2nd,-cbore1)
p215=kp(rbore,theta2nd,-(topthk+ribh))
p216=kp(rbore,theta2nd,-cbore2)
!
l,p11,p21               ! create lines to form volumes
l,p12,p22
l,p13,p23
l,p14,p24
l,p15,p25
!
l,p16,p26
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
l,p113,p213
l,p114,p214
l,p115,p215
!
l,p116,p216
!
*get,lm,line,,num,max
!
lsel,s,line,,lm-15,lm,1,0
lesize,all,,,1,1,
!
lsel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Create Areas <<<<<<<<<<<<<<<<<<<<<<<
!
!
csys,gcs
!
asel,s,loc,y,theta1st,theta1st,,
!
*get,am11,area,,num,min
!
am12=am11+1
am13=am11+2
am14=am11+3
am15=am11+4
am16=am11+5
am17=am11+6
!
asel,s,loc,y,theta2nd,theta2nd,,
!
*get,am21,area,,num,min
!
am22=am21+1
am23=am21+2
am24=am21+3
am25=am21+4
am26=am21+5
am27=am21+6
!
a,p11,p12,p22,p21
a,p11,p110,p210,p21
a,p13,p110,p210,p23
a,p12,p13,p23,p22
a,p13,p14,p24,p23
!
a,p14,p18,p28,p24
a,p19,p18,p28,p29
a,p110,p19,p29,p210
a,p14,p15,p25,p24
a,p15,p16,p26,p25
!
a,p17,p16,p26,p27
a,p18,p17,p27,p28
a,p16,p111,p211,p26
a,p112,p111,p211,p212
a,p17,p112,p212,p27
!
a,p113,p112,p212,p213
a,p18,p113,p213,p28
a,p114,p113,p213,p214
a,p19,p114,p214,p29
a,p111,p115,p215,p211
!
a,p116,p115,p215,p216
a,p112,p116,p216,p212
!
asel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<< Create Volumes <<<<<<<<<<<<<<<<<<<<<<<
!
*get,am,area,,num,max
!
va,am11,am21,am,am-1,am-2,am-8          ! create volumes
va,am12,am22,am-7,am-8,am-9,am-11
va,am13,am23,am-5,am-6,am-7,am-10
va,am14,am24,am-3,am-4,am-5,am-15
va,am15,am25,am-16,am-10,am-11,am-12,am-13
!
va,am16,am26,am-19,am-17,am-16,am-15,am-14
va,am17,am27,am-21,am-20,am-19,am-18
!
*get,vm,volu,,num,max
!
vsel,s,volu,,vm-6,vm,,
!
accat,am-15,am-14
accat,am-10,am-11
!
vatt,1,,1,0             ! set element attributes                
mshkey,1                ! select mapped meshing
vmesh,all               ! mesh areas
!
vsel,all
!
adele,am+1              ! delete concatenated areas
adele,am+2
!
*get,lm,line,,num,max
!
ldele,lm
ldele,lm-1
ldele,lm-2
ldele,lm-3
!
LPLOT
!
!
/com, ************************* Areas 10 & 11 *************************
!
rib1st=rib5cs               ! first rib defining volume
rib2nd=rib6cs               ! second rib defining volume
!
theta1st=theta10            ! first angle defining volume
theta2nd=theta11            ! second angle defining volume
! 
csys,rib1st
!
p11=kp(-(w1/2+r1),ribh+topthk,0)    ! define keypoints forming volume lines 
p12=kp(-(w1/2+r1),ribh,0)
!
csys,gcs
!
p11r=kx(p11)
A=p11r-r1*cos(45*rconv)
B=r1-r1*sin(45*rconv)
p13=kp(A,theta1st,-(topthk+B))
!
csys,rib2nd
!
p21=kp(-(w1/2+r1),ribh+topthk,0)
p22=kp(-(w1/2+r1),ribh,0)
!
csys,gcs
!
p21r=kx(p21)
C=p21r-r1*cos(45*rconv)
D=r1-r1*sin(45*rconv)
p23=kp(C,theta2nd,-(topthk+D))
!
p14=kp(rbout,theta1st,-(topthk+r1))
p15=kp(rbout,theta1st,-(topthk+ribh))
p16=kp(rcb1,theta1st,-(topthk+ribh))
p17=kp(rcb1,theta1st,-cbore2)
p18=kp(rcb1,theta1st,-(topthk+r1))
p19=kp(rcb1,theta1st,-cbore1)
p110=kp(rcb1,theta1st,0)
p111=kp(rcb2,theta1st,-(topthk+ribh))
p112=kp(rcb2,theta1st,-cbore2)
p113=kp(rcb2,theta1st,-(topthk+r1))
p114=kp(rcb2,theta1st,-cbore1)
p115=kp(rbore,theta1st,-(topthk+ribh))
p116=kp(rbore,theta1st,-cbore2)
!
p24=kp(rbout,theta2nd,-(topthk+r1))
p25=kp(rbout,theta2nd,-(topthk+ribh))
p26=kp(rcb1,theta2nd,-(topthk+ribh))
p27=kp(rcb1,theta2nd,-cbore2)
p28=kp(rcb1,theta2nd,-(topthk+r1))
p29=kp(rcb1,theta2nd,-cbore1)
p210=kp(rcb1,theta2nd,0)
p211=kp(rcb2,theta2nd,-(topthk+ribh))
p212=kp(rcb2,theta2nd,-cbore2)
p213=kp(rcb2,theta2nd,-(topthk+r1))
p214=kp(rcb2,theta2nd,-cbore1)
p215=kp(rbore,theta2nd,-(topthk+ribh))
p216=kp(rbore,theta2nd,-cbore2)
!
l,p11,p21               ! create lines to form volumes
l,p12,p22
l,p13,p23
l,p14,p24
l,p15,p25
!
l,p16,p26
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
l,p113,p213
l,p114,p214
l,p115,p215
!
l,p116,p216
!
*get,lm,line,,num,max
!
lsel,s,line,,lm-15,lm,1,0
lesize,all,,,1,1,
!
lsel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Create Areas <<<<<<<<<<<<<<<<<<<<<<<
!
!
csys,gcs
!
asel,s,loc,y,theta1st,theta1st,,
!
*get,am11,area,,num,min
!
am12=am11+1
am13=am11+2
am14=am11+3
am15=am11+4
am16=am11+5
am17=am11+6
!
asel,s,loc,y,theta2nd,theta2nd,,
!
*get,am21,area,,num,min
!
am22=am21+1
am23=am21+2
am24=am21+3
am25=am21+4
am26=am21+5
am27=am21+6
!
a,p11,p12,p22,p21
a,p11,p110,p210,p21
a,p13,p110,p210,p23
a,p12,p13,p23,p22
a,p13,p14,p24,p23
!
a,p14,p18,p28,p24
a,p19,p18,p28,p29
a,p110,p19,p29,p210
a,p14,p15,p25,p24
a,p15,p16,p26,p25
!
a,p17,p16,p26,p27
a,p18,p17,p27,p28
a,p16,p111,p211,p26
a,p112,p111,p211,p212
a,p17,p112,p212,p27
!
a,p113,p112,p212,p213
a,p18,p113,p213,p28
a,p114,p113,p213,p214
a,p19,p114,p214,p29
a,p111,p115,p215,p211
!
a,p116,p115,p215,p216
a,p112,p116,p216,p212
!
asel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<< Create Volumes <<<<<<<<<<<<<<<<<<<<<<<
!
*get,am,area,,num,max
!
va,am11,am21,am,am-1,am-2,am-8          ! create volumes
va,am12,am22,am-7,am-8,am-9,am-11
va,am13,am23,am-5,am-6,am-7,am-10
va,am14,am24,am-3,am-4,am-5,am-15
va,am15,am25,am-16,am-10,am-11,am-12,am-13
!
va,am16,am26,am-19,am-17,am-16,am-15,am-14
va,am17,am27,am-21,am-20,am-19,am-18
!
*get,vm,volu,,num,max
!
vsel,s,volu,,vm-6,vm,,
!
accat,am-15,am-14
accat,am-10,am-11
!
vatt,1,,1,0             ! set element attributes                
mshkey,1                ! select mapped meshing
vmesh,all               ! mesh areas
!
vsel,all
!
adele,am+1              ! delete concatenated areas
adele,am+2
!
*get,lm,line,,num,max
!
ldele,lm
ldele,lm-1
ldele,lm-2
ldele,lm-3
!
LPLOT
!
!
/com, ************************* Areas 12 & 13 *************************
!
rib1st=rib6cs               ! first rib defining volume
rib2nd=rib7cs               ! second rib defining volume
!
theta1st=theta12            ! first angle defining volume
theta2nd=theta13            ! second angle defining volume
! 
csys,rib1st
!
p11=kp((w1/2+r1),ribh+topthk,0)     ! define keypoints forming volume lines 
p12=kp((w1/2+r1),ribh,0)
!
csys,gcs
!
p11r=kx(p11)
A=p11r-r1*cos(45*rconv)
B=r1-r1*sin(45*rconv)
p13=kp(A,theta1st,-(topthk+B))
!
csys,rib2nd
!
p21=kp(-(w1/2+r1),ribh+topthk,0)
p22=kp(-(w1/2+r1),ribh,0)
!
csys,gcs
!
p21r=kx(p21)
C=p21r-r1*cos(45*rconv)
D=r1-r1*sin(45*rconv)
p23=kp(C,theta2nd,-(topthk+D))
!
p14=kp(rbout,theta1st,-(topthk+r1))
p15=kp(rbout,theta1st,-(topthk+ribh))
p16=kp(rcb1,theta1st,-(topthk+ribh))
p17=kp(rcb1,theta1st,-cbore2)
p18=kp(rcb1,theta1st,-(topthk+r1))
p19=kp(rcb1,theta1st,-cbore1)
p110=kp(rcb1,theta1st,0)
p111=kp(rcb2,theta1st,-(topthk+ribh))
p112=kp(rcb2,theta1st,-cbore2)
p113=kp(rcb2,theta1st,-(topthk+r1))
p114=kp(rcb2,theta1st,-cbore1)
p115=kp(rbore,theta1st,-(topthk+ribh))
p116=kp(rbore,theta1st,-cbore2)
!
p24=kp(rbout,theta2nd,-(topthk+r1))
p25=kp(rbout,theta2nd,-(topthk+ribh))
p26=kp(rcb1,theta2nd,-(topthk+ribh))
p27=kp(rcb1,theta2nd,-cbore2)
p28=kp(rcb1,theta2nd,-(topthk+r1))
p29=kp(rcb1,theta2nd,-cbore1)
p210=kp(rcb1,theta2nd,0)
p211=kp(rcb2,theta2nd,-(topthk+ribh))
p212=kp(rcb2,theta2nd,-cbore2)
p213=kp(rcb2,theta2nd,-(topthk+r1))
p214=kp(rcb2,theta2nd,-cbore1)
p215=kp(rbore,theta2nd,-(topthk+ribh))
p216=kp(rbore,theta2nd,-cbore2)
!
l,p11,p21               ! create lines to form volumes
l,p12,p22
l,p13,p23
l,p14,p24
l,p15,p25
!
l,p16,p26
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
l,p113,p213
l,p114,p214
l,p115,p215
!
l,p116,p216
!
*get,lm,line,,num,max
!
lsel,s,line,,lm-15,lm,1,0
lesize,all,,,3,1,
!
lsel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Create Areas <<<<<<<<<<<<<<<<<<<<<<<
!
csys,gcs
!
asel,s,loc,y,theta1st,theta1st,,
!
*get,am11,area,,num,min
!
am12=am11+1
am13=am11+2
am14=am11+3
am15=am11+4
am16=am11+5
am17=am11+6
!
asel,s,loc,y,theta2nd,theta2nd,,
!
*get,am21,area,,num,min
!
am22=am21+1
am23=am21+2
am24=am21+3
am25=am21+4
am26=am21+5
am27=am21+6
!
a,p11,p12,p22,p21
a,p11,p110,p210,p21
a,p13,p110,p210,p23
a,p12,p13,p23,p22
a,p13,p14,p24,p23
!
a,p14,p18,p28,p24
a,p19,p18,p28,p29
a,p110,p19,p29,p210
a,p14,p15,p25,p24
a,p15,p16,p26,p25
!
a,p17,p16,p26,p27
a,p18,p17,p27,p28
a,p16,p111,p211,p26
a,p112,p111,p211,p212
a,p17,p112,p212,p27
!
a,p113,p112,p212,p213
a,p18,p113,p213,p28
a,p114,p113,p213,p214
a,p19,p114,p214,p29
a,p111,p115,p215,p211
!
a,p116,p115,p215,p216
a,p112,p116,p216,p212
!
asel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<< Create Volumes <<<<<<<<<<<<<<<<<<<<<<<
!
*get,am,area,,num,max
!
va,am11,am21,am,am-1,am-2,am-8          ! create volumes
va,am12,am22,am-7,am-8,am-9,am-11
va,am13,am23,am-5,am-6,am-7,am-10
va,am14,am24,am-3,am-4,am-5,am-15
va,am15,am25,am-16,am-10,am-11,am-12,am-13
!
va,am16,am26,am-19,am-17,am-16,am-15,am-14
va,am17,am27,am-21,am-20,am-19,am-18
!
*get,vm,volu,,num,max
!
vsel,s,volu,,vm-6,vm,,
!
accat,am-15,am-14
accat,am-10,am-11
!
vatt,1,,1,0             ! set element attributes                
mshkey,1                ! select mapped meshing
vmesh,all               ! mesh areas
!
vsel,all
!
adele,am+1              ! delete concatenated areas
adele,am+2
!
*get,lm,line,,num,max
!
ldele,lm
ldele,lm-1
ldele,lm-2
ldele,lm-3
!
LPLOT
!
!
/com, ************************* Areas 14 & 1 *************************
!
rib1st=rib7cs               ! first rib defining volume
rib2nd=rib1cs               ! second rib defining volume
!
theta1st=theta14            ! first angle defining volume
theta2nd=theta1             ! second angle defining volume
! 
csys,rib1st
!
p11=kp((w1/2+r1),ribh+topthk,0)     ! define keypoints forming volume lines 
p12=kp((w1/2+r1),ribh,0)
!
csys,gcs
!
p11r=kx(p11)
A=p11r-r1*cos(45*rconv)
B=r1-r1*sin(45*rconv)
p13=kp(A,theta1st,-(topthk+B))
!
csys,rib2nd
!
p21=kp(-(w1/2+r1),ribh+topthk,0)
p22=kp(-(w1/2+r1),ribh,0)
!
csys,gcs
!
p21r=kx(p21)
C=p21r-r1*cos(45*rconv)
D=r1-r1*sin(45*rconv)
p23=kp(C,theta2nd,-(topthk+D))
!
p14=kp(rbout,theta1st,-(topthk+r1))
p15=kp(rbout,theta1st,-(topthk+ribh))
p16=kp(rcb1,theta1st,-(topthk+ribh))
p17=kp(rcb1,theta1st,-cbore2)
p18=kp(rcb1,theta1st,-(topthk+r1))
p19=kp(rcb1,theta1st,-cbore1)
p110=kp(rcb1,theta1st,0)
p111=kp(rcb2,theta1st,-(topthk+ribh))
p112=kp(rcb2,theta1st,-cbore2)
p113=kp(rcb2,theta1st,-(topthk+r1))
p114=kp(rcb2,theta1st,-cbore1)
p115=kp(rbore,theta1st,-(topthk+ribh))
p116=kp(rbore,theta1st,-cbore2)
!
p24=kp(rbout,theta2nd,-(topthk+r1))
p25=kp(rbout,theta2nd,-(topthk+ribh))
p26=kp(rcb1,theta2nd,-(topthk+ribh))
p27=kp(rcb1,theta2nd,-cbore2)
p28=kp(rcb1,theta2nd,-(topthk+r1))
p29=kp(rcb1,theta2nd,-cbore1)
p210=kp(rcb1,theta2nd,0)
p211=kp(rcb2,theta2nd,-(topthk+ribh))
p212=kp(rcb2,theta2nd,-cbore2)
p213=kp(rcb2,theta2nd,-(topthk+r1))
p214=kp(rcb2,theta2nd,-cbore1)
p215=kp(rbore,theta2nd,-(topthk+ribh))
p216=kp(rbore,theta2nd,-cbore2)
!
l,p11,p21               ! create lines to form volumes
l,p12,p22
l,p13,p23
l,p14,p24
l,p15,p25
!
l,p16,p26
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
l,p113,p213
l,p114,p214
l,p115,p215
!
l,p116,p216
!
*get,lm,line,,num,max
!
lsel,s,line,,lm-15,lm,1,0
lesize,all,,,6,1,
!
lsel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Create Areas <<<<<<<<<<<<<<<<<<<<<<<
!
!
csys,gcs
!
asel,s,loc,y,theta1st,theta1st,,
!
*get,am11,area,,num,min
!
am12=am11+1
am13=am11+2
am14=am11+3
am15=am11+4
am16=am11+5
am17=am11+6
!
asel,s,loc,y,theta2nd,theta2nd,,
!
*get,am21,area,,num,min
!
am22=am21+1
am23=am21+2
am24=am21+3
am25=am21+4
am26=am21+5
am27=am21+6
!
a,p11,p12,p22,p21
a,p11,p110,p210,p21
a,p13,p110,p210,p23
a,p12,p13,p23,p22
a,p13,p14,p24,p23
!
a,p14,p18,p28,p24
a,p19,p18,p28,p29
a,p110,p19,p29,p210
a,p14,p15,p25,p24
a,p15,p16,p26,p25
!
a,p17,p16,p26,p27
a,p18,p17,p27,p28
a,p16,p111,p211,p26
a,p112,p111,p211,p212
a,p17,p112,p212,p27
!
a,p113,p112,p212,p213
a,p18,p113,p213,p28
a,p114,p113,p213,p214
a,p19,p114,p214,p29
a,p111,p115,p215,p211
!
a,p116,p115,p215,p216
a,p112,p116,p216,p212
!
asel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<< Create Volumes <<<<<<<<<<<<<<<<<<<<<<<
!
*get,am,area,,num,max
!
va,am11,am21,am,am-1,am-2,am-8          ! create volumes
va,am12,am22,am-7,am-8,am-9,am-11
va,am13,am23,am-5,am-6,am-7,am-10
va,am14,am24,am-3,am-4,am-5,am-15
va,am15,am25,am-16,am-10,am-11,am-12,am-13
!
va,am16,am26,am-19,am-17,am-16,am-15,am-14
va,am17,am27,am-21,am-20,am-19,am-18
!
*get,vm,volu,,num,max
!
vsel,s,volu,,vm-6,vm,,
!
accat,am-15,am-14
accat,am-10,am-11
!
vatt,1,,1,0             ! set element attributes                
mshkey,1                ! select mapped meshing
vmesh,all               ! mesh areas
!
vsel,all
!
adele,am+1              ! delete concatenated areas
adele,am+2
!
*get,lm,line,,num,max
!
ldele,lm
ldele,lm-1
ldele,lm-2
ldele,lm-3
!
LPLOT
!
!
/com, ************************* Areas 8 & 9 *************************
!
rib1st=rib4cs               ! first rib defining volume
rib2nd=rib5cs               ! second rib defining volume
!
theta1st=theta8             ! first angle defining volume
theta2nd=theta9             ! second angle defining volume
! 
csys,rib1st
!
p11=kp((w1/2+r2),fribh+ftopthk,.1423)   ! define keypoints forming volume lines 
p12=kp((w1/2+r2),fribh,.1423)
!
csys,gcs
!
p11r=kx(p11)
A=p11r-r2*cos(45*rconv)
B=r2-r2*sin(45*rconv)
p13=kp(A,theta1st,-(ftopthk+B))
!
csys,rib2nd
!
p21=kp((w1/2+r2),fribh+ftopthk,-.1423)
p22=kp((w1/2+r2),fribh,-.1423)
!
csys,gcs
!
p21r=kx(p21)
C=p21r-r2*cos(45*rconv)
D=r2-r2*sin(45*rconv)
p23=kp(C,theta2nd,-(ftopthk+D))
!
p14=kp(p11r-r2,theta1st,-(ftopthk+r2))
p15=kp(p11r-r2,theta1st,-(ftopthk+fribh))
p16=kp(rcb1,theta1st,-(ftopthk+fribh))
p17=kp(rcb1,theta1st,-cbore2)
p18=kp(rcb1,theta1st,-ftopthk)
p19=kp(rcb1,theta1st,-cbore1)
p110=kp(rcb1,theta1st,0)
p111=kp(rcb2,theta1st,-(ftopthk+fribh))
p112=kp(rcb2,theta1st,-cbore2)
p113=kp(rcb2,theta1st,-ftopthk)
p114=kp(rcb2,theta1st,-cbore1)
p115=kp(rbore,theta1st,-(ftopthk+fribh))
p116=kp(rbore,theta1st,-cbore2)
!
p24=kp(p21r-r2,theta2nd,-(ftopthk+r2))
p25=kp(p21r-r2,theta2nd,-(ftopthk+fribh))
p26=kp(rcb1,theta2nd,-(ftopthk+fribh))
p27=kp(rcb1,theta2nd,-cbore2)
p28=kp(rcb1,theta2nd,-ftopthk)
p29=kp(rcb1,theta2nd,-cbore1)
p210=kp(rcb1,theta2nd,0)
p211=kp(rcb2,theta2nd,-(ftopthk+fribh))
p212=kp(rcb2,theta2nd,-cbore2)
p213=kp(rcb2,theta2nd,-ftopthk)
p214=kp(rcb2,theta2nd,-cbore1)
p215=kp(rbore,theta2nd,-(ftopthk+fribh))
p216=kp(rbore,theta2nd,-cbore2)
!
l,p11,p21               ! create lines to form volumes
l,p12,p22
l,p13,p23
l,p14,p24
l,p15,p25
!
l,p16,p26
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
l,p113,p213
l,p114,p214
l,p115,p215
!
l,p116,p216
!
csys,gcs

lsel,s,loc,y,theta8,theta8,,        ! select lines at specified angle
*get,lm,line,,num,max
!
lccat,lm-20,lm-19           ! concatenate lines for mapped meshing
lccat,lm-6,lm-4
!
lsel,all
!
asel,s,loc,y,theta8,theta8,,        ! select areas at specified angle
*get,am,area,,num,max
!
aatt,1,,2,0             ! set element attributes
!
mshkey,1                ! select mapped meshing
amesh,am-6,am,1             ! mesh  selected areas
!
asel,all
!
*get,lm,line,,num,max
!
ldele,lm                ! delete concatenated lines
ldele,lm-1
!
*get,lm,line,,num,max
!
esize,0,14              ! set element divisions along drag
!
vdrag,am-6,,,,,,lm,         ! drag area meshes into volumes
vdrag,am-5,,,,,,lm-4,
vdrag,am-4,,,,,,lm-4,
vdrag,am-3,,,,,,lm-2,
vdrag,am-2,,,,,,lm-10,
vdrag,am-1,,,,,,lm-7,
vdrag,am,,,,,,lm-15,
!
asel,s,loc,y,theta8,theta8,,        ! select areas and clear area meshes
aclear,all
asel,all
!
nummrg,all,.001,.001            ! merge coincident entities
numcmp,all              ! renumber & compress all entities
!
LPLOT
!
/com, **************************** Fill In Rib Gaps **************************
!
!
/com, ********************************* Rib 1 ********************************
!
!
ribcs=rib1cs                ! set active rib c.s.
theta1st=theta1             ! boundary angles of rib
theta2nd=theta2
!
csys,ribcs  !<<<<<<<<<<<<<<<<<<<<<< begin creation of endpoint keypoints
!
cpbot=kp(0,0,0)             ! top & bottom rib center points
cptop=kp(0,topthk+ribh,0)
!
p1=kp(-(w1/2+r1),ribh+topthk,0)     ! top area enpoints
p2=kp((w1/2+r1),ribh+topthk,0)
!
p11=kp(-w1/2,0,0)           ! bottom area enpoints
p21=kp(w1/2,0,0)
!
csys,gcs
!
p1r=kx(p1)              ! radius of outermost points
p2r=kx(p2)
!
p3=kp(rcb1,theta1st,0)          ! remaining top area enpoints
p4=kp(rcb1,theta2nd,0)
!
p12=kp(rbout,theta1st,-(topthk+ribh))   ! remaining bottom area endpoints   
p13=kp(rcb1,theta1st,-(topthk+ribh))
p22=kp(rbout,theta2nd,-(topthk+ribh))
p23=kp(rcb1,theta2nd,-(topthk+ribh))
!
p14=kp(rcb1,theta1st,-cbore2)       ! remaining endpoints for vdrag
p15=kp(rcb1,theta1st,-(topthk+r1))
p16=kp(rcb1,theta1st,-cbore1)
p17=kp(rcb2,theta1st,-(topthk+ribh))
p18=kp(rcb2,theta1st,-cbore2)
p19=kp(rcb2,theta1st,-(topthk+r1))
p110=kp(rcb2,theta1st,-cbore1)
p111=kp(rbore,theta1st,-(topthk+ribh))
p112=kp(rbore,theta1st,-cbore2)
!
p24=kp(rcb1,theta2nd,-cbore2)
p25=kp(rcb1,theta2nd,-(topthk+r1))
p26=kp(rcb1,theta2nd,-cbore1)
p27=kp(rcb2,theta2nd,-(topthk+ribh))
p28=kp(rcb2,theta2nd,-cbore2)
p29=kp(rcb2,theta2nd,-(topthk+r1))
p210=kp(rcb2,theta2nd,-cbore1)
p211=kp(rbore,theta2nd,-(topthk+ribh))
p212=kp(rbore,theta2nd,-cbore2)
!
l,p12,p22               ! creation of boundary lines
l,p13,p23
l,p14,p24
l,p15,p25
l,p16,p26
!
l,p3,p4
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
!
*get,lm,line,,num,max
!
asel,s,loc,y,theta1st,theta1st,,
asel,r,loc,x,rbore,rcb1,,
!
*get,am,area,,num,max
!
vatt,1,,1,0             ! set element attributes
esize,0,4               ! set element divisions along drag
!
a,cpbot,p11,p12,p22,p21         ! create top & bottom areas
a,p12,p13,p23,p22
a,p1,p3,p4,p2,cptop
a,p3,p4,p26,p16
!
vdrag,am-3,,,,,,lm,         ! drag areas into volumes
vdrag,am-2,,,,,,lm-4,
vdrag,am-1,,,,,,lm-4,
vdrag,am,,,,,,lm-3,
!
nummrg,all,,.001
!
asel,s,loc,y,theta2nd,theta1st,,    ! select out areas to create volume
asel,r,loc,x,rcb1,p1r,,
!
va,all
!
mopt,pyra,on                ! set mesh option to transitional
mshape,1,3d
shpp,warn               ! set element shape to 3d
!
*get,vm,volu,,num,max
!
vatt,1,,1,0             ! set element attributes
mshkey,0                ! select free meshing
vmesh,vm,vm,                ! mesh volume
!
LPLOT
!
!
/com, ********************************* Rib 2 ********************************
!
!
ribcs=rib2cs                ! set active rib c.s.
theta1st=theta3             ! boundary angles of rib
theta2nd=theta4
!
csys,ribcs  !<<<<<<<<<<<<<<<<<<<<<< begin creation of endpoint keypoints
!
cpbot=kp(0,0,0)             ! top & bottom rib center points
cptop=kp(0,topthk+ribh,0)
!
p1=kp(-(w1/2+r1),ribh+topthk,0)     ! top area enpoints
p2=kp((w1/2+r1),ribh+topthk,0)
!
p11=kp(-w1/2,0,0)           ! bottom area enpoints
p21=kp(w1/2,0,0)
!
csys,gcs
!
p1r=kx(p1)              ! radius of outermost points
p2r=kx(p2)
!
p3=kp(rcb1,theta1st,0)          ! remaining top area enpoints
p4=kp(rcb1,theta2nd,0)
!
p12=kp(rbout,theta1st,-(topthk+ribh))   ! remaining bottom area endpoints   
p13=kp(rcb1,theta1st,-(topthk+ribh))
p22=kp(rbout,theta2nd,-(topthk+ribh))
p23=kp(rcb1,theta2nd,-(topthk+ribh))
!
p14=kp(rcb1,theta1st,-cbore2)       ! remaining endpoints for vdrag
p15=kp(rcb1,theta1st,-(topthk+r1))
p16=kp(rcb1,theta1st,-cbore1)
p17=kp(rcb2,theta1st,-(topthk+ribh))
p18=kp(rcb2,theta1st,-cbore2)
p19=kp(rcb2,theta1st,-(topthk+r1))
p110=kp(rcb2,theta1st,-cbore1)
p111=kp(rbore,theta1st,-(topthk+ribh))
p112=kp(rbore,theta1st,-cbore2)
!
p24=kp(rcb1,theta2nd,-cbore2)
p25=kp(rcb1,theta2nd,-(topthk+r1))
p26=kp(rcb1,theta2nd,-cbore1)
p27=kp(rcb2,theta2nd,-(topthk+ribh))
p28=kp(rcb2,theta2nd,-cbore2)
p29=kp(rcb2,theta2nd,-(topthk+r1))
p210=kp(rcb2,theta2nd,-cbore1)
p211=kp(rbore,theta2nd,-(topthk+ribh))
p212=kp(rbore,theta2nd,-cbore2)
!
l,p12,p22               ! creation of boundary lines
l,p13,p23
l,p14,p24
l,p15,p25
l,p16,p26
!
l,p3,p4
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
!
*get,lm,line,,num,max
!
asel,s,loc,y,theta1st,theta1st,,
asel,r,loc,x,rbore,rcb1,,
!
*get,am,area,,num,max
!
vatt,1,,1,0             ! set element attributes
esize,0,4               ! set element divisions along drag
!
a,cpbot,p11,p12,p22,p21         ! create top & bottom areas
a,p12,p13,p23,p22
a,p1,p3,p4,p2,cptop
a,p3,p4,p26,p16
!
vdrag,am-3,,,,,,lm,         ! drag areas into volumes
vdrag,am-2,,,,,,lm-4,
vdrag,am-1,,,,,,lm-4,
vdrag,am,,,,,,lm-3,
!
nummrg,all,,.001
!
asel,s,loc,y,theta2nd,theta1st,,    ! select out areas to create volume
asel,r,loc,x,rcb1,p1r,,
!
va,all
!
mopt,pyra,on                ! set mesh option to transitional
mshape,1,3d
shpp,warn               ! set element shape to 3d
!
*get,vm,volu,,num,max
!
vatt,1,,1,0             ! set element attributes
mshkey,0                ! select free meshing
vmesh,vm,vm,                ! mesh volume
!
LPLOT
!
!
/com, ********************************* Rib 3 ********************************
!
!
ribcs=rib3cs                ! set active rib c.s.
theta1st=theta5             ! boundary angles of rib
theta2nd=theta6
!
csys,ribcs  !<<<<<<<<<<<<<<<<<<<<<< begin creation of endpoint keypoints
!
cpbot=kp(0,0,0)             ! top & bottom rib center points
cptop=kp(0,topthk+ribh,0)
!
p1=kp(-(w2/2+r1),ribh+topthk,0)     ! top area enpoints
p2=kp((w2/2+r1),ribh+topthk,0)
!
p11=kp(-w2/2,0,0)           ! bottom area enpoints
p21=kp(w2/2,0,0)
!
csys,gcs
!
p1r=kx(p1)              ! radius of outermost points
p2r=kx(p2)
!
p3=kp(rcb1,theta1st,0)          ! remaining top area enpoints
p4=kp(rcb1,theta2nd,0)
!
p12=kp(rbout,theta1st,-(topthk+ribh))   ! remaining bottom area endpoints   
p13=kp(rcb1,theta1st,-(topthk+ribh))
p22=kp(rbout,theta2nd,-(topthk+ribh))
p23=kp(rcb1,theta2nd,-(topthk+ribh))
!
p14=kp(rcb1,theta1st,-cbore2)       ! remaining endpoints for vdrag
p15=kp(rcb1,theta1st,-(topthk+r1))
p16=kp(rcb1,theta1st,-cbore1)
p17=kp(rcb2,theta1st,-(topthk+ribh))
p18=kp(rcb2,theta1st,-cbore2)
p19=kp(rcb2,theta1st,-(topthk+r1))
p110=kp(rcb2,theta1st,-cbore1)
p111=kp(rbore,theta1st,-(topthk+ribh))
p112=kp(rbore,theta1st,-cbore2)
!
p24=kp(rcb1,theta2nd,-cbore2)
p25=kp(rcb1,theta2nd,-(topthk+r1))
p26=kp(rcb1,theta2nd,-cbore1)
p27=kp(rcb2,theta2nd,-(topthk+ribh))
p28=kp(rcb2,theta2nd,-cbore2)
p29=kp(rcb2,theta2nd,-(topthk+r1))
p210=kp(rcb2,theta2nd,-cbore1)
p211=kp(rbore,theta2nd,-(topthk+ribh))
p212=kp(rbore,theta2nd,-cbore2)
!
l,p12,p22               ! creation of boundary lines
l,p13,p23
l,p14,p24
l,p15,p25
l,p16,p26
!
l,p3,p4
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
!
*get,lm,line,,num,max
!
asel,s,loc,y,theta1st,theta1st,,
asel,r,loc,x,rbore,rcb1,,
!
*get,am,area,,num,max
!
vatt,1,,1,0             ! set element attributes
esize,0,4               ! set element divisions along drag
!
a,cpbot,p11,p12,p22,p21         ! create top & bottom areas
a,p12,p13,p23,p22
a,p1,p3,p4,p2,cptop
a,p3,p4,p26,p16
!
vdrag,am-3,,,,,,lm,         ! drag areas into volumes
vdrag,am-2,,,,,,lm-4,
vdrag,am-1,,,,,,lm-4,
vdrag,am,,,,,,lm-3,
!
nummrg,all,,.001
!
asel,s,loc,y,theta2nd,theta1st,,    ! select out areas to create volume
asel,r,loc,x,rcb1,p1r,,
!
va,all
!
mopt,pyra,on                ! set mesh option to transitional
mshape,1,3d
shpp,warn               ! set element shape to 3d
!
*get,vm,volu,,num,max
!
vatt,1,,1,0             ! set element attributes
mshkey,0                ! select free meshing
vmesh,vm,vm,                ! mesh volume
!
LPLOT
!
!
/com, ********************************* Rib 6 ********************************
!
!
ribcs=rib6cs                ! set active rib c.s.
theta1st=theta11            ! boundary angles of rib
theta2nd=theta12
!
csys,ribcs  !<<<<<<<<<<<<<<<<<<<<<< begin creation of endpoint keypoints
!
cpbot=kp(0,0,0)             ! top & bottom rib center points
cptop=kp(0,topthk+ribh,0)
!
p1=kp(-(w2/2+r1),ribh+topthk,0)     ! top area enpoints
p2=kp((w2/2+r1),ribh+topthk,0)
!
p11=kp(-w2/2,0,0)           ! bottom area enpoints
p21=kp(w2/2,0,0)
!
csys,gcs
!
p1r=kx(p1)              ! radius of outermost points
p2r=kx(p2)
!
p3=kp(rcb1,theta1st,0)          ! remaining top area enpoints
p4=kp(rcb1,theta2nd,0)
!
p12=kp(rbout,theta1st,-(topthk+ribh))   ! remaining bottom area endpoints   
p13=kp(rcb1,theta1st,-(topthk+ribh))
p22=kp(rbout,theta2nd,-(topthk+ribh))
p23=kp(rcb1,theta2nd,-(topthk+ribh))
!
p14=kp(rcb1,theta1st,-cbore2)       ! remaining endpoints for vdrag
p15=kp(rcb1,theta1st,-(topthk+r1))
p16=kp(rcb1,theta1st,-cbore1)
p17=kp(rcb2,theta1st,-(topthk+ribh))
p18=kp(rcb2,theta1st,-cbore2)
p19=kp(rcb2,theta1st,-(topthk+r1))
p110=kp(rcb2,theta1st,-cbore1)
p111=kp(rbore,theta1st,-(topthk+ribh))
p112=kp(rbore,theta1st,-cbore2)
!
p24=kp(rcb1,theta2nd,-cbore2)
p25=kp(rcb1,theta2nd,-(topthk+r1))
p26=kp(rcb1,theta2nd,-cbore1)
p27=kp(rcb2,theta2nd,-(topthk+ribh))
p28=kp(rcb2,theta2nd,-cbore2)
p29=kp(rcb2,theta2nd,-(topthk+r1))
p210=kp(rcb2,theta2nd,-cbore1)
p211=kp(rbore,theta2nd,-(topthk+ribh))
p212=kp(rbore,theta2nd,-cbore2)
!
l,p12,p22               ! creation of boundary lines
l,p13,p23
l,p14,p24
l,p15,p25
l,p16,p26
!
l,p3,p4
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
!
*get,lm,line,,num,max
!
asel,s,loc,y,theta1st,theta1st,,
asel,r,loc,x,rbore,rcb1,,
!
*get,am,area,,num,max
!
vatt,1,,1,0             ! set element attributes
esize,0,4               ! set element divisions along drag
!
a,cpbot,p11,p12,p22,p21         ! create top & bottom areas
a,p12,p13,p23,p22
a,p1,p3,p4,p2,cptop
a,p3,p4,p26,p16
!
vdrag,am-3,,,,,,lm,         ! drag areas into volumes
vdrag,am-2,,,,,,lm-4,
vdrag,am-1,,,,,,lm-4,
vdrag,am,,,,,,lm-3,
!
nummrg,all,,.001
!
asel,s,loc,y,360+theta1st,theta2nd,,    ! select out areas to create volume
asel,r,loc,x,rcb1,p1r,,
!
va,all
!
mopt,pyra,on                ! set mesh option to transitional
mshape,1,3d
shpp,warn               ! set element shape to 3d
!
*get,vm,volu,,num,max
!
vatt,1,,1,0             ! set element attributes
mshkey,0                ! select free meshing
vmesh,vm,vm,                ! mesh volume
!
LPLOT
!
!
/com, ********************************* Rib 7 ********************************
!
!
ribcs=rib7cs                ! set active rib c.s.
theta1st=theta13            ! boundary angles of rib
theta2nd=theta14
!
csys,ribcs  !<<<<<<<<<<<<<<<<<<<<<< begin creation of endpoint keypoints
!
cpbot=kp(0,0,0)             ! top & bottom rib center points
cptop=kp(0,topthk+ribh,0)
!
p1=kp(-(w1/2+r1),ribh+topthk,0)     ! top area enpoints
p2=kp((w1/2+r1),ribh+topthk,0)
!
p11=kp(-w1/2,0,0)           ! bottom area enpoints
p21=kp(w1/2,0,0)
!
csys,gcs
!
p1r=kx(p1)              ! radius of outermost points
p2r=kx(p2)
!
p3=kp(rcb1,theta1st,0)          ! remaining top area enpoints
p4=kp(rcb1,theta2nd,0)
!
p12=kp(rbout,theta1st,-(topthk+ribh))   ! remaining bottom area endpoints   
p13=kp(rcb1,theta1st,-(topthk+ribh))
p22=kp(rbout,theta2nd,-(topthk+ribh))
p23=kp(rcb1,theta2nd,-(topthk+ribh))
!
p14=kp(rcb1,theta1st,-cbore2)       ! remaining endpoints for vdrag
p15=kp(rcb1,theta1st,-(topthk+r1))
p16=kp(rcb1,theta1st,-cbore1)
p17=kp(rcb2,theta1st,-(topthk+ribh))
p18=kp(rcb2,theta1st,-cbore2)
p19=kp(rcb2,theta1st,-(topthk+r1))
p110=kp(rcb2,theta1st,-cbore1)
p111=kp(rbore,theta1st,-(topthk+ribh))
p112=kp(rbore,theta1st,-cbore2)
!
p24=kp(rcb1,theta2nd,-cbore2)
p25=kp(rcb1,theta2nd,-(topthk+r1))
p26=kp(rcb1,theta2nd,-cbore1)
p27=kp(rcb2,theta2nd,-(topthk+ribh))
p28=kp(rcb2,theta2nd,-cbore2)
p29=kp(rcb2,theta2nd,-(topthk+r1))
p210=kp(rcb2,theta2nd,-cbore1)
p211=kp(rbore,theta2nd,-(topthk+ribh))
p212=kp(rbore,theta2nd,-cbore2)
!
l,p12,p22               ! creation of boundary lines
l,p13,p23
l,p14,p24
l,p15,p25
l,p16,p26
!
l,p3,p4
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
!
*get,lm,line,,num,max
!
asel,s,loc,y,theta1st,theta1st,,
asel,r,loc,x,rbore,rcb1,,
!
*get,am,area,,num,max
!
vatt,1,,1,0             ! set element attributes
esize,0,4               ! set element divisions along drag
!
a,cpbot,p11,p12,p22,p21         ! create top & bottom areas
a,p12,p13,p23,p22
a,p1,p3,p4,p2,cptop
a,p3,p4,p26,p16
!
vdrag,am-3,,,,,,lm,         ! drag areas into volumes
vdrag,am-2,,,,,,lm-4,
vdrag,am-1,,,,,,lm-4,
vdrag,am,,,,,,lm-3,
!
nummrg,all,,.001
!
asel,s,loc,y,theta2nd,theta1st,,    ! select out areas to create volume
asel,r,loc,x,rcb1,p1r,,
!
va,all
!
mopt,pyra,on                ! set mesh option to transitional
mshape,1,3d
shpp,warn               ! set element shape to 3d
!
*get,vm,volu,,num,max
!
vatt,1,,1,0             ! set element attributes
mshkey,0                ! select free meshing
vmesh,vm,vm,                ! mesh volume
!
LPLOT
!
!
/com, ********************************* Rib 4 ********************************
!
!
ribcs=rib4cs                ! set active rib c.s.
theta1st=theta7             ! boundary angles of rib
theta2nd=theta8
!
csys,ribcs  !<<<<<<<<<<<<<<<<<<<<<<<<< begin creation of endpoint keypoints
!
cpbot=kp(0,0,0)             ! front, top & bottom rib center points
cpbotfnt=kp(0,0,.1423)
cptop=kp(0,topthk+ribh,0)
cptopfnt=kp(0,topthk+ribh,.1423)
!
p1=kp(-(w1/2+r1),ribh+topthk,0)     ! top area enpoints
p2=kp((w1/2+r1),ribh+topthk,.1423)
!
p11=kp(-w1/2,0,0)           ! bottom area enpoints
p21=kp(w1/2,0,.1423)
!
csys,gcs
!
p1r=kx(p1)              ! radius of outermost points
p2r=kx(p2)
thetacp=ky(cpbot)
!
p3=kp(rcb1,theta1st,0)          ! remaining top area enpoints
p4=kp(rcb1,theta2nd,0)
!
p12=kp(rbout,theta1st,-(topthk+ribh))   ! remaining bottom area endpoints   
p13=kp(rcb1,theta1st,-(topthk+ribh))
p22=kp(rbout,theta2nd,-(topthk+ribh))
p23=kp(rcb1,theta2nd,-(topthk+ribh))
!
p14=kp(rcb1,theta1st,-cbore2)       ! remaining endpoints for vdrag
p15=kp(rcb1,theta1st,-(topthk+r1))
p16=kp(rcb1,theta1st,-cbore1)
p17=kp(rcb2,theta1st,-(topthk+ribh))
p18=kp(rcb2,theta1st,-cbore2)
p19=kp(rcb2,theta1st,-(topthk+r1))
p110=kp(rcb2,theta1st,-cbore1)
p111=kp(rbore,theta1st,-(topthk+ribh))
p112=kp(rbore,theta1st,-cbore2)
!
p24=kp(rcb1,theta2nd,-cbore2)
p25=kp(rcb1,theta2nd,-(topthk+r1))
p26=kp(rcb1,theta2nd,-cbore1)
p27=kp(rcb2,theta2nd,-(topthk+ribh))
p28=kp(rcb2,theta2nd,-cbore2)
p29=kp(rcb2,theta2nd,-(topthk+r1))
p210=kp(rcb2,theta2nd,-cbore1)
p211=kp(rbore,theta2nd,-(topthk+ribh))
p212=kp(rbore,theta2nd,-cbore2)
!
l,p12,p22               ! creation of boundary lines
l,p13,p23
l,p14,p24
l,p15,p25
l,p16,p26
!
l,p3,p4
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
!
*get,lm,line,,num,max
!
asel,s,loc,y,theta1st,theta1st,,
asel,r,loc,x,rbore,rcb1,,
!
*get,am,area,,num,max
!
vatt,1,,1,0             ! set element attributes
esize,0,4               ! set element divisions along drag
!
a,cpbot,p11,p12,p22,p21,cpbotfnt            ! create top & bottom areas
a,p12,p13,p23,p22
a,p1,p3,p4,p2,cptopfnt,cptop
a,p3,p4,p26,p16
!
vdrag,am-3,,,,,,lm,         ! drag areas into volumes
vdrag,am-2,,,,,,lm-4,
vdrag,am-1,,,,,,lm-4,
vdrag,am,,,,,,lm-3,
!
nummrg,all,,.001
!
csys,ribcs              ! select out offset areas on rib
!
asel,s,loc,x,0,w1/2+r2,,
asel,r,loc,z,0,0,,
asel,r,loc,y,0,fribh+ftopthk,,
cm,aoff1,area
!
asel,s,loc,x,w1/2,w1/2+r2,,
asel,r,loc,z,0,.1423,,
asel,r,loc,y,0,fribh+ftopthk,,
cm,aoff2,area
!
asel,s,loc,y,fribh-r2,fribh-r2,,
asel,r,loc,z,0,.1423,,
asel,r,loc,x,0,w1/2,,
cm,aoff3,area
!
asel,s,loc,y,0,0,,
asel,r,loc,z,0,.1423,,
asel,r,loc,x,0,w1/2,,
cm,aoff4,area
!
csys,gcs
!
asel,s,loc,y,theta2nd,theta1st,,    ! select out areas to create volume
asel,r,loc,x,rcb1,p1r,,
asel,u,area,,aoff1,,,
asel,u,area,,aoff2,,,
asel,u,area,,aoff3,,,
asel,u,area,,aoff4,,,
asel,u,area,,93,,,
!
va,all
!
mopt,pyra,on                ! set mesh option to transitional
mshape,1,3d
shpp,warn               ! set element shape to 3d
!
*get,vm,volu,,num,max
!
vatt,1,,1,0             ! set element attributes
mshkey,0                ! select free meshing
vmesh,vm,vm,                ! mesh volume
!
LPLOT
!
!
/com, ********************************* Rib 5 ********************************
!
!
ribcs=rib5cs                ! set active rib c.s.
theta1st=theta9             ! boundary angles of rib
theta2nd=theta10
!
csys,ribcs  !<<<<<<<<<<<<<<<<<<<<<<<<< begin creation of endpoint keypoints
!
cpbot=kp(0,0,0)             ! front, top & bottom rib center points
cpbotfnt=kp(0,0,-.1423)
cptop=kp(0,topthk+ribh,0)
cptopfnt=kp(0,topthk+ribh,-.1423)
!
p1=kp(-(w1/2+r1),ribh+topthk,0)     ! top area enpoints
p2=kp((w1/2+r1),ribh+topthk,-.1423)
!
p11=kp(-w1/2,0,0)           ! bottom area enpoints
p21=kp(w1/2,0,-.1423)
!
csys,gcs
!
p1r=kx(p1)              ! radius of outermost points
p2r=kx(p2)
thetacp=ky(cpbot)
!
p3=kp(rcb1,theta2nd,0)          ! remaining top area enpoints
p4=kp(rcb1,theta1st,0)
!
p12=kp(rbout,theta2nd,-(topthk+ribh))   ! remaining bottom area endpoints   
p13=kp(rcb1,theta2nd,-(topthk+ribh))
p22=kp(rbout,theta1st,-(topthk+ribh))
p23=kp(rcb1,theta1st,-(topthk+ribh))
!
p14=kp(rcb1,theta2nd,-cbore2)       ! remaining endpoints for vdrag
p15=kp(rcb1,theta2nd,-(topthk+r1))
p16=kp(rcb1,theta2nd,-cbore1)
p17=kp(rcb2,theta2nd,-(topthk+ribh))
p18=kp(rcb2,theta2nd,-cbore2)
p19=kp(rcb2,theta2nd,-(topthk+r1))
p110=kp(rcb2,theta2nd,-cbore1)
p111=kp(rbore,theta2nd,-(topthk+ribh))
p112=kp(rbore,theta2nd,-cbore2)
!
p24=kp(rcb1,theta1st,-cbore2)
p25=kp(rcb1,theta1st,-(topthk+r1))
p26=kp(rcb1,theta1st,-cbore1)
p27=kp(rcb2,theta1st,-(topthk+ribh))
p28=kp(rcb2,theta1st,-cbore2)
p29=kp(rcb2,theta1st,-(topthk+r1))
p210=kp(rcb2,theta1st,-cbore1)
p211=kp(rbore,theta1st,-(topthk+ribh))
p212=kp(rbore,theta1st,-cbore2)
!
l,p12,p22               ! creation of boundary lines
l,p13,p23
l,p14,p24
l,p15,p25
l,p16,p26
!
l,p3,p4
l,p17,p27
l,p18,p28
l,p19,p29
l,p110,p210
!
l,p111,p211
l,p112,p212
!
*get,lm,line,,num,max
!
asel,s,loc,y,theta1st,theta1st,,
asel,r,loc,x,rbore,rcb1,,
!
*get,am,area,,num,max
!
vatt,1,,1,0             ! set element attributes
esize,0,4               ! set element divisions along drag
!
a,cpbot,p11,p12,p22,p21,cpbotfnt    ! bottom outer area
a,p12,p13,p23,p22           ! bottom inner area
a,p1,p3,p4,p2,cptopfnt,cptop        ! top area
a,p3,p4,p26,p16             ! top vertical area
!
vdrag,am-3,,,,,,lm,         ! drag areas into volumes
vdrag,am-2,,,,,,lm-4,
vdrag,am-1,,,,,,lm-4,
vdrag,am,,,,,,lm-3,
!
nummrg,all,,.001
!
csys,ribcs              ! select out offset areas on rib
!
asel,s,loc,x,0,w1/2+r2,,
asel,r,loc,z,0,0,,
asel,r,loc,y,0,fribh+ftopthk,,
cm,aoff1,area
!
asel,s,loc,x,w1/2,w1/2+r2,,
asel,r,loc,z,0,-.1423,,
asel,r,loc,y,0,fribh+ftopthk,,
cm,aoff2,area
!
asel,s,loc,y,fribh-r2,fribh-r2,,
asel,r,loc,z,0,-.1423,,
asel,r,loc,x,0,w1/2,,
cm,aoff3,area
!
asel,s,loc,y,0,0,,
asel,r,loc,z,0,-.1423,,
asel,r,loc,x,0,w1/2,,
cm,aoff4,area
!
csys,gcs
!
asel,s,loc,y,theta2nd,theta1st,,    ! select out areas to create volume
asel,r,loc,x,rcb1,p1r,,
asel,u,area,,aoff1,,,
asel,u,area,,aoff2,,,
asel,u,area,,aoff3,,,
asel,u,area,,aoff4,,,
asel,u,area,,127,,,
!
va,all
!
mopt,pyra,on                ! set mesh option to transitional
mshape,1,3d             ! set element shape to 3d
shpp,warn               ! set element shape checks to 'warn'                
!
*get,vm,volu,,num,max
!
vatt,1,,1,0             ! set element attributes
mshkey,0                ! select free meshing
vmesh,vm,vm,                ! mesh volume
!
asel,all
lsel,all
!
LPLOT
!
!