/COM,ANSYS RELEASE  5.6    UP19991022       10:17:11    02/09/2000
!
FINISH
/CLEAR,START
! 
/dev,font,1,charter,bold,r,26,,,,,
/title,AAAV Modal Analysis / Aluminum
/units,bin
/prep7
!
!   TITLE:  aaaval56.mac (AAAV)
!
!**************************** 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
!
!------------------------------- AlBe -----------------------------
!
UIMP,1,EX, , ,27.0E6,
UIMP,1,DENS, , ,.076/386.4,
UIMP,1,NUXY, , ,0.17,
!
!-------------------------------- Al ------------------------------
!
UIMP,2,EX, , ,10.4E6,
UIMP,2,DENS, , ,.097/386.4,
UIMP,2,NUXY, , ,0.33,
!
!************************ VARIABLE/CONSTANTS LIST ******************
!
pi=3.1415927
rconv=pi/180        ! radian conversion factor from degrees
!
w1=2.000        ! Rib width (rib 1)
w2=1.00         ! Rib width (ribs 2 thru 7)
wf=2.1298/2     ! front VEXT surface width (rib 1 only)
wb=4.2612/2     ! back VEXT surface width (rib 1 only)
hb=.8468        ! back VEXT surface heigth (rib 1 only)
!
hf=.4125        ! front VEXT surface heigth (rib 1 only)
f2b=.4064       ! distance from front to back of rib areas (rib 1 only)
r1=.25          ! Rib fillet radius
ribh=.636       ! Rib height
extd=1.1736     ! extrusion distance for VEXT command
!
r2=.125         ! Front rib fillet radius
fribh=ribh-r1       ! Height of front surface ribs
topthk=.50      ! Top surface thickness
ftopthk=.65     ! Front top surface thickness
cbr=1.750       ! Center Boss Radius
!
bthk=.050       ! Boss 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
waldistl=6.065      ! left wall inner surface distance from global c.s. center
!
netop=2         ! # elements through top surface thickness
netopw=4        ! # elements across top surface of rib 1 (1/2 rib)
neribw=5        ! # elements across rib 1 (1/2 rib)
!
neribh=netop+1      ! # elements through rib 1 height
neftop=4        ! # elements through front top surface thickness
nefill=3        ! # elements in fillet radius (must be even #)
neffill=nefill      ! # elements in front rib fillet
neL1=1          ! # elements in rib 1 length
!
neL1ext=5       ! # elements along rib 1 VEXT
neL2=7          ! # elements in rib 2 length (must be odd #)
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=7          ! # elements in rib 7 length
!
L1=.2384        ! Length of rib 1 (normal extrusion segment)
L2=3.2662-.0711-.0868   ! 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.1651       ! Point 1 x
r1p1y=-1.136        ! Point 1 y
r1p1z=1.00      ! Point 1 z
!
r1p2x=-2.1651       ! Point 2 x
r1p2y=-1.136        ! Point 2 y
r1p2z=-1.00     ! Point 2 z
!
r1p3x=-2.1651       ! Point 3 x
r1p3y=-.750     ! Point 3 y
r1p3z=1.00      ! Point 3 z
!
!   RIB 2 Point Coordinates
!
r2p1x=-1.7032       ! Point 1 x
r2p1y=-1.136        ! Point 1 y
r2p1z=-1.8386       ! Point 1 z
!
r2p2x=-.8488        ! Point 2 x
r2p2y=-1.136        ! Point 2 y
r2p2z=-2.3582       ! Point 2 z
!
r2p3x=-1.7032       ! Point 3 x
r2p3y=-.750     ! Point 3 y
r2p3z=-1.8386       ! Point 3 z
!
r2p4x=-.750     ! Point 4 x Intersection point with rib 3
r2p4y=-.500     ! Point 4 y
r2p4z=-2.6768       ! Point 4 z
!
!   RIB 7 Point Coordinates
!
r7p1x=-.8248        ! Point 1 x
r7p1y=-1.136        ! Point 1 y
r7p1z=2.3667        ! Point 1 z
!
r7p2x=-1.6845       ! Point 2 x
r7p2y=-1.136        ! Point 2 y
r7p2z=1.8558        ! Point 2 z
!
r7p3x=-.8248        ! Point 3 x
r7p3y=-.750     ! Point 3 y
r7p3z=2.3667        ! Point 3 z
!
r7p4x=-.750     ! Point 4 x
r7p4y=-.500     ! Point 4 y
r7p4z=2.7302        ! Point 4 z
!
!*********************** 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,,
!
!***************************** Top Surface Layout ********************** 
!
csys,0
k,,-3.8770,0,-4.8195
k,,-1.4052,0,-6.585
k,,3.1218,0,-6.585
k,,3.1218,0,6.115
k,,-1.4052,0,6.115
k,,-3.8770,0,4.3494
*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
!
!******************************* End Top Surface ***********************   
!
!**************************** 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,,,neribw,1,
lesize,lm-2,,,neribh,1,
lesize,lm-3,,,netop+1,1,
lesize,lm-4,,,netopw,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,neribh,1,
lesize,lm-7,,,neribw,1,
!
*get,am,area,,num,max
!
aatt,,,2,0          ! set element attributes                
mshkey,0            ! select free meshing
amesh,am-1,am,1         ! mesh areas
!
!
ESIZE,0,neL1,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,2               ! 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
!
/VIEW, 1 ,1,1,1
/ANG, 1
/DIST, 1 ,1.371742,1
/REP,FAST
/DIST, 1 ,1.371742,1
/REP,FAST
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
k,,0,-.0265,-f2b
k,,wf,-.0265,-f2b
k,,wf,ribh-r1,-f2b
k,,r1+wf,ribh,-f2b
k,,r1+wf,ribh+topthk,-f2b
k,,0,ribh+topthk,-f2b
k,,0,ribh-r1,-f2b
!
*get,km,kp,,num,max
!
l,km-6,km-5
l,km-5,km-4
l,km-3,km-2
l,km-2,km-1
l,km-1,km
!
l,km,km-6
l,km,km-4
!
k,,r1+wf,hf,-f2b
!
csn=csn+1
cskp,csn,1,km+1,km-4,km-3,,
!
l,km-4,km-3
!
*get,lm,line,,num,max
!
lesize,lm-7,,,neribw,1,
lesize,lm-6,,,neribh,1,
lesize,lm-5,,,netop,1,
lesize,lm-4,,,netopw,1,
lesize,lm-3,,,netop+1,1,
!
lesize,lm-2,,,neribh,1,
lesize,lm-1,,,neribw,1,
lesize,lm,,,nefill,1,
!
a,km,km-4,km-5,km-6
a,km,km-1,km-2,km-3,km-4
!
*get,am,area,,num,max
!
aatt,,,2,0              ! set element attributes                
mshkey,0                ! select free meshing
amesh,am-1,am,1             ! mesh areas
!
ESIZE,,neL1ext,         ! # element divisions along extrusion
!
csn=csn+1
vextcs1=csn
cskp,vextcs1,0,km,km-4,km-6,,
!
vext,am-1,,,,,extd,wb/wf,hb/hf,
!
csn=csn+1
vextcs2=csn
cskp,vextcs2,0,km-1,km-2,km,,
!
vext,am,,,,,extd,wb/wf,,
!
csys,rib1cs
!
lstr,kp(w1/2+r1,ribh+topthk,-L1),kp(wf+r1,ribh+topthk,-f2b)
lstr,kp(w1/2+r1,ribh,-L1),kp(wf+r1,ribh,-f2b)
lstr,kp(w1/2,ribh-r1,-L1),kp(wf,ribh-r1,-f2b)
lstr,kp(w1/2,0,-L1),kp(wf,-.0265,-f2b)
lstr,kp(0,0,-L1),kp(0,-.0265,-f2b)
!
lstr,kp(0,ribh-r1,-L1),kp(0,ribh-r1,-f2b)
lstr,kp(0,ribh+topthk,-L1),kp(0,ribh+topthk,-f2b)
!
*get,lm,line,,num,max
lsel,s,line,,lm-6,lm,1,
lesize,all,,,1,1,
!
lsel,all
!
a,kp(w1/2+r1,ribh+topthk,-L1),kp(wf+r1,ribh+topthk,-f2b),kp(wf+r1,ribh,-f2b),kp(w1/2+r1,ribh,-L1)
a,kp(w1/2+r1,ribh,-L1),kp(wf+r1,ribh,-f2b),kp(wf,ribh-r1,-f2b),kp(w1/2,ribh-r1,-L1)
a,kp(w1/2,ribh-r1,-L1),kp(wf,ribh-r1,-f2b),kp(wf,-.0265,-f2b),kp(w1/2,0,-L1)
a,kp(w1/2,0,-L1),kp(wf,-.0265,-f2b),kp(0,-.0265,-f2b),kp(0,0,-L1)
a,kp(0,0,-L1),kp(0,-.0265,-f2b),kp(0,ribh-r1,-f2b),kp(0,ribh-r1,-L1)
!
a,kp(0,ribh-r1,-L1),kp(0,ribh-r1,-f2b),kp(0,ribh+topthk,-f2b),kp(0,ribh+topthk,-L1)
a,kp(0,ribh+topthk,-L1),kp(0,ribh+topthk,-f2b),kp(wf+r1,ribh+topthk,-f2b),kp(w1/2+r1,ribh+topthk,-L1)
!
asel,s,loc,z,-L1,-f2b,,
va,all
!
asel,all
!
mopt,pyra,on            ! set mesh option to transitional
mshape,1,3d         ! set element shape to 3d
shpp,warn
!
*get,vm,volu,,num,max
!
vatt,2,,1,0
mshkey,0
vmesh,vm,vm,
!
vsymm,x,all,,,,0,0
!
nummrg,all,.001,.001
numcmp,all
!
LPLOT
!
!*************************** End First (Mid) Rib Layout ***************
!
!******************************** 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
!
x4=r2p4x
z4=r2p4z
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< MATH SECTION >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
topw=w2+2*r1                ! top of rib width
!
m=(r2p2x-r2p1x)/(r2p2z-r2p1z)       ! slope in XZ plane
!
C=sqrt(1+(1/m**2))          ! defining a constant
!
DELx45=topw/C
x5=x4-DELx45
!
DELz45=(1/m)*(x4-x5)
z5=z4-DELz45
!
DELx46=r1/C
x6=x4-DELx46
z6=z4-(1/m)*DELx46
!
DELx57=r1/C
x7=x5+DELx57
z7=z5+(1/m)*DELx57
!
d1=sqrt(((z6-r2p2z)**2)+((x6-r2p2x)**2))
!
d1prime=sqrt(((z7-r2p1z)**2)+((x7-r2p1x)**2))
!
check1=d1/d1prime
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< END MATH SECTION >>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
*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
!
kdele,km+1          ! delete center keypoint
!
k,km+1,0,0,-d1          ! rib center point in cross-section plane
k,,0,ribh-r1,-d1        ! remaining rib cross-section keypoints
k,,w2/2,0,-d1
k,,w2/2,ribh-r1,-d1
k,,w2/2+r1,ribh,-d1
! 
k,,w2/2+r1,ribh+topthk,-d1
k,,0,ribh+topthk,-d1
!
*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,-d1      ! 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,,,neribw-1,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,netop+1,1,
lesize,lm-4,,,netopw-1,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,netop,1,
lesize,lm-7,,,neribw-1,1,
!
*get,am,area,,num,max
!
aatt,,,2,0          ! set element attributes                
mshkey,0            ! select free meshing
amesh,am-1,am,1         ! mesh areas
!
ESIZE,0,neL2,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,2               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,L2-(d1-.0868),,   ! extrude areas into 1st volumes
VOFFST,am-1,L2-(d1-.0868),,
!
ESIZE,0,1,          ! # element divisions 2nd along extrusion
!
VOFFST,am,-(d1-.0868),,     ! extrude areas into 2nd volumes
VOFFST,am-1,-(d1-.0868),,
!
csys,rib2cs
!
asel,s,loc,z,-.0868,-.0868,,    ! select rear areas for extrusion
*get,a1,area,,num,max
*get,a2,area,,num,min
!
ESIZE,0,1,          ! # element divisions 3rd extrusion
!
VOFFST,a1,-.0868,,      ! extrude areas into 3rd volumes
VOFFST,a2,-.0868,,
!
asel,all
!
aclear,am-1,am,1        ! clear area meshes
!
csys,rib2cs         ! mirror volumes to create complete rib
!
vsel,s,loc,x,0,w2/2+r1,,    ! select volumes to be mirrored
vsel,r,loc,z,0,-(L2+.0868),,
!
vsymm,x,all,,,,0,0
!
vsel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
csys,0
!
k,,-3.3634,-2.2829,-4.5686  ! create keypoint for rib triangular area
*get,km,kp,,num,max
!
csys,rib2cs
!
l,kp(-w2/2,0,-.0868),km     ! lines forming rib
l,km,kp(-w2/2,0,-(L2+.0868))
!
*get,lm,line,,num,max
!
lesize,lm,,,neL2+1,1,       ! size element lines
lesize,lm-1,,,neL2+1,1,
!
a,km,kp(-w2/2,0,-(L2+.0868)),kp(-w2/2,0,-d1),kp(-w2/2,0,-.0868)
!                 area containing rib planform
!
lsel,s,loc,x,-w2/2,-w2/2,,
lsel,r,loc,z,-.0868,-(L2+.0868),,
lsel,r,loc,y,0,0,,
!
lccat,all           ! concatenate lines
!
*get,am,area,,num,max
!
aatt,,,2,0          ! set element attributes                
mshkey,1            ! select mapped meshing
amesh,am,am,            ! mesh areas
!
*get,lm,line,,num,max       ! delete concatenated lines
ldele,lm
lsel,all
!
ESIZE,0,2*(neribw-1),       ! # element divisions along extrusion
!
type,1              ! set element type
mat,2               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,w2,,          ! extrude areas into volume
!
aclear,am,am,           ! clear area meshes
!
nummrg,all,.001,.001
numcmp,all
!
vsel,s,volu,,13,16,1
vsel,a,volu,,19,22,1
vclear,all
vsel,all
!
LPLOT
!
!**************************** End 2nd Rib Layout  ************************
!
!**************************** 3rd/4th Rib Layout *************************
!
csys,0
!
k,,.625,-.625,-4.7473       ! common points of ribs 3 & 4
k,,.500,-.625,-4.7473
k,,.625,-.500,-4.7473
k,,.7457,-.625,-4.7150
!
p1=kp(.625,-.625,-4.7473)   ! assigning keypoint numbers to parameters
p2=kp(.500,-.625,-4.7473)
p3=kp(.625,-.500,-4.7473)
p8=kp(.7457,-.625,-4.7150)
!
csn=csn+1           ! create spherical c.s. between ribs 3 & 4
spher34=csn
cskp,spher34,2,p1,p2,p3,,
!
csn=csn+1           ! create rib 3 c.s.
rib3cs=csn
cskp,rib3cs,0,p2,p1,p3,,
!
csn=csn+1           ! create rib 4 c.s.
rib4cs=csn
cskp,rib4cs,0,p1,p8,p3,,
!
csys,rib3cs         ! rib 3 keypoints
!
k,,0,-(ribh-r2),0
k,,-w2,-(ribh-r2),0
k,,-w2,-r2,0
k,,-(w2+r1),r2,0
k,,-(w2+r1),(r2+topthk),0
k,,r2,(r2+topthk),0
!
*get,p7,kp,,num,max     ! create rib 3 lines
!
k,,-(w2+r1),-r2,0       ! front fillet center keypoint
!
l,p2,p7-5
l,p7-5,p7-4
l,p7-4,p7-3
l,p7-1,p7-2
l,p7,p7-1
l,p7,p3
!
csys,spher34
!
l,p2,p3             ! fillet line for rib 3/common to rib 4
!
csn=csn+1           ! front fillet c.s.         
cskp,csn,1,p7+1,p7-3,p7-2,,
!
l,p7-3,p7-2         ! front fillet line
!
a,p3,p2,p7-5,p7-4,p7-3,p7-2,p7-1,p7
!
*get,rib3area,area,,num,max ! assign rib 3 area parameter
!
csys,rib4cs         ! rib 4 keypoints
!
k,,r2,-(ribh-r2),0
k,,(w2+r1+r2),-(ribh-r2),0
k,,(w2+r1+r2),(r2+topthk),0
!
*get,km,kp,,num,max     ! create rib 4 lines
!
l,p8,km-2
l,km-2,km-1
l,km-1,km
l,km,p7
!
csn=csn+1
cskp,csn,1,p1,p8,p3,,
!
l,p8,p3             ! fillet line for rib 4/common to rib 3
!
a,p3,p8,km-2,km-1,km,p7     ! creation of rib 6 area
!
*get,rib4area,area,,num,max ! assign rib 6 area parameter
!
csys,0              ! find extrude distance for rib 5 front
!
x1=-.500
z1=-2.4559
x2=1.3602
z2=-2.1786
!
x3=kx(p7-4)
z3=kz(p7-4)
x4=kx(km-1)
z4=kz(km-1)
!
dfext3=sqrt(((x3-x1)**2)+((z3-z1)**2))  ! calculation of front extrusion
dfext4=sqrt(((x4-x2)**2)+((z4-z2)**2))  ! distances
!
drext3=2.500*cos(20*rconv)
drext4=drext3
!
VOFFST,rib3area,drext3,,        ! extrude area into rear volume
VOFFST,rib3area,-(dfext3-.0868),,   ! extrude area into front volume
!
VOFFST,rib4area,-drext4,,       ! extrude area into rear volume
VOFFST,rib4area,(dfext4-.0868),,    ! extrude area into front volume
!
nummrg,all,.001,.001            ! merge all coincident entities
numcmp,all              ! renumber & compress all entities
!
LPLOT
!
!<<<<<<<<<<<<<<<<<<<<<<< create rib 3 triangular section >>>>>>>>>>>>>>>>>>>
!
csys,rib3cs
!
dext3=(dfext3-.0868)+drext3         ! total rib 3 length
!
k,,-1.00,-(.511+dext3*tan(20*rconv)),-drext3    ! section keypoint
!
*get,pt5,kp,,num,max
!
l,pt5,kp(-1.00,-.511,-drext3)           ! section lines
l,pt5,kp(-1.00,-.511,dfext3-.0868)
!
a,pt5,kp(-1.00,-.511,-drext3),p7-4,kp(-1.00,-.511,dfext3-.0868)
!
*get,am,area,,num,max
!
VOFFST,am,1.00,,            ! extrude area into rear volume
!
!<<<<<<<<<<<<<<<<<<<<<<< create rib 4 triangular section >>>>>>>>>>>>>>>>>>>
!
csys,rib4cs
!
dext4=(dfext4-.0868)+drext4         ! total rib 6 length
!
k,,.125,-(.511+dext4*tan(20*rconv)),-drext4     ! section keypoint
!
*get,pt6,kp,,num,max
!
l,pt6,kp(.125,-.511,-drext4)            ! section lines
l,pt6,kp(.125,-.511,dfext4-.0868)
!
a,pt6,kp(.125,-.511,-drext4),km-2,kp(.125,-.511,dfext4-.0868)
!
*get,am,area,,num,max
!
VOFFST,am,1.25,,            ! extrude area into rear volume
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< END TRIANGULAR SECTIONS >>>>>>>>>>>>>>>>>>>>>>>>>
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< ADD INTERSECTING VOLUMES >>>>>>>>>>>>>>>>>>>>>>>>>
!
*get,vm,volu,,num,max
!
vadd,vm-2,vm-4
!
btol,.01
vadd,vm,vm-1,vm+1
btol,defa
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
csys,rib3cs         ! partition rib triangular sections
!
k,,0,-3.00,0            ! keypoints of dividing area
k,,-1.00,-3.00,0
!
*get,km,kp,,num,max
!
a,km,km-1,kp(0,-.511,0),kp(-1.00,-.511,0)   ! create dividing area
!
*get,am,area,,num,max
!
vsel,s,loc,z,0,(dfext3-.0868),,
vsel,r,loc,x,-(1.00+r1),(2*r2+1.00+r1),
!
vsba,all,am,,delete,delete  ! divide volumes
!
numcmp,all
!
csys,rib4cs         ! partition rib triangular section 6
!
k,,.125,-3.00,0         ! keypoints of dividing area
k,,1.375,-3.00,0
!
*get,km,kp,,num,max
!
a,km,km-1,kp(.125,-.511,0),kp(1.375,-.511,0)    ! create dividing area
!
*get,am,area,,num,max
!
vsel,s,loc,z,0,(dfext3-.0868),,
vsel,r,loc,x,-(1.00+r1+r2),(r2+1.25),
!
vsba,all,am,,delete,delete      ! partition volume
!
vsel,all
!
*get,km,kp,,num,max
lfillt,259,261,.125,km+1        ! create line fillet
r34filcp=km+1
!
csys,0                  ! rib 3/4 fillet c.p. x,y,z corrdinates
r34cpx=kx(r34filcp)
r34cpy=ky(r34filcp)
r34cpz=kz(r34filcp)
kdele,r34filcp
!
csys,rib4cs
!
nummrg,all,.001,.001            ! merge entities
numcmp,all              ! renumber & compress all entities
!
a,145,104,106               ! create filleted area
!
VOFFST,168,.511,,           ! extrude area into lower volumes
VOFFST,169,3.00,,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
csn=csn+1           ! create partitioning c.s.
cskp,csn,1,126,149,150,,
!
k,,3.0,330,0            ! partition area keypoints
k,,3.0,90,0
k,,3.0,210,0
!
*get,km,kp,,num,max
!
l,km-2,km-1         ! partition area lines
l,km-1,km
l,km,km-2
!
a,km-2,km-1,km          ! partition area lines
!
*get,am,area,,num,max
*get,vm,volu,,num,max
!
vsba,vm,am,,delete,delete   ! partition volume
!
vdele,32,,,1            ! delete unused volume
!
nummrg,all,.001,.001
numcmp,all
!
vsel,s,volu,,28,        ! select volume for modification
aslv,s
cm,aset1,area
vdele,28
vsel,all
asel,all
adele,151,153,1,1
!
btol,.001
asbl,172,273,,delete,keep   ! partition area
btol,defa
!
nummrg,all,.001,.001
!
csys,spher34            ! spherical c.s.
!
a,104,105,106           ! create spherical areas
!
vdele,29,30,1,          ! delete unused volumes
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
adele,148,150,1,1       ! delete unnecessary areas
adele,158,,,1
adele,133,,,1
adele,165,,,1
!
adele,168,,,1
adele,167,,,1
!
a,111,112,116,137,145,130   ! replace top areas
!
asel,s,area,,aset1,     ! select areas for volume recreation 
asel,a,area,,133
asel,a,area,,151,152
asel,a,area,,169,171
cm,aset34,area          ! area set for rib 3/4
!
va,all              ! recreate volumes
!
vsel,all
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
LPLOT
!
!************************** End 3rd/4th Rib Layout ***********************
!
!
!**************************** 5th/6th Rib Layout *************************
!
csys,0
!
k,,.625,-.625,4.7473        ! common points of ribs 5 & 6
k,,.7457,-.625,4.7150
k,,.625,-.500,4.7473
k,,.500,-.625,4.7473
!
p1=kp(.625,-.625,4.7473)    ! assigning keypoint numbers to parameters
p2=kp(.7457,-.625,4.7150)
p3=kp(.625,-.500,4.7473)
p8=kp(.500,-.625,4.7473)
!
csn=csn+1           ! create spherical c.s. between ribs 5 & 6
spher56=csn
cskp,spher56,2,p1,p2,p3,,
!
csn=csn+1           ! create rib 5 c.s.
rib5cs=csn
cskp,rib5cs,0,p2,p1,p3,,
!
csn=csn+1           ! create rib 6 c.s.
rib6cs=csn
cskp,rib6cs,0,p1,p8,p3,,
!
csys,rib5cs         ! rib 5 keypoints
!
k,,0,-(ribh-r2),0
k,,-(w2+r1),-(ribh-r2),0
k,,-(w2+r1),(r2+topthk),0
k,,r2,(r2+topthk),0
!
*get,p7,kp,,num,max     ! create rib 5 lines
!
l,p2,p7-3
l,p7-3,p7-2
l,p7-2,p7-1
l,p7-1,p7
l,p7,p3
!
csys,spher56
!
l,p2,p3             ! fillet line for rib 5/common to rib 6
!
a,p3,p2,p7-3,p7-2,p7-1,p7
!
*get,rib5area,area,,num,max ! assign rib 5 area parameter
!
csys,rib6cs         ! rib 6 keypoints
!
k,,r2,-(ribh-r2),0
k,,w2+r2,-(ribh-r2),0
k,,w2+r2,-r2,0
k,,w2+r2+r1,r2,0
k,,w2+r2+r1,(r2+topthk),0
!
k,,w2+r2+r1,-r2,0
!
*get,km,kp,,num,max     ! create rib 6 lines
!
l,p8,km-5
l,km-5,km-4
l,km-4,km-3
l,km-2,km-1
l,km-1,p7
!
csys,spher56
!
l,p8,p3                 ! fillet line for rib 6/common to rib 5
!
csn=csn+1
cskp,csn,1,km,km-3,km-2,,       ! fillet c.s. for rib 6
!
l,km-3,km-2             ! fillet line
!
a,p3,p7,km-1,km-2,km-3,km-4,km-5,p8 ! creation of rib 6 area
!
*get,rib6area,area,,num,max     ! assign rib 6 area parameter
!
csys,0                  ! find extrude distance for rib 5 front
!
x1=1.3643
z1=2.1938
x2=-.500
z2=2.4559
!
x3=kx(p7-2)
z3=kz(p7-2)
x4=kx(km-4)
z4=kz(km-4)
!
dfext5=sqrt(((x3-x1)**2)+((z3-z1)**2))  ! calculation of front extrusion
dfext6=sqrt(((x4-x2)**2)+((z4-z2)**2))  ! distances
!
drext5=2.000*cos(20*rconv)
drext6=drext5
!
VOFFST,rib6area,drext5,,        ! extrude area into rear volume
VOFFST,rib6area,-(dfext6-.0868),,   ! extrude area into front volume
!
VOFFST,rib5area,drext6,,        ! extrude area into rear volume
VOFFST,rib5area,-(dfext5-.0868),,   ! extrude area into front volume
!
nummrg,all,.001,.001            ! merge all coincident entities
numcmp,all              ! renumber & compress all entities
!
LPLOT
!
!<<<<<<<<<<<<<<<<<<<<<<< create rib 5 triangular section >>>>>>>>>>>>>>>>>>>
!
csys,rib5cs
!
dext5=(dfext5-.0868)+drext5 ! total rib 5 length
!
k,,-1.25,-(.511+dext5*tan(20*rconv)),-drext5  ! section keypoint
!
*get,pt5,kp,,num,max
!
l,pt5,kp(-1.25,-.511,-drext5)       ! section lines
l,pt5,kp(-1.25,-.511,dfext5-.0868)
!
a,pt5,kp(-1.25,-.511,-drext5),p7-2,kp(-1.25,-.511,dfext5-.0868)
!
*get,am,area,,num,max
!
VOFFST,am,1.25,,        ! extrude area into rear volume
!
!<<<<<<<<<<<<<<<<<<<<<<< create rib 6 triangular section >>>>>>>>>>>>>>>>>>>
!
csys,rib6cs
!
dext6=(dfext6-.0868)+drext6         ! total rib 6 length
!
k,,1.125,-(.511+dext6*tan(20*rconv)),-drext6    ! section keypoint
!
*get,pt6,kp,,num,max
!
l,pt6,kp(1.125,-.511,-drext6)           ! section lines
l,pt6,kp(1.125,-.511,dfext6-.0868)
!
a,pt6,kp(1.125,-.511,-drext6),km-4,kp(1.125,-.511,dfext6-.0868)
!
*get,am,area,,num,max
!
VOFFST,am,-1.00,,       ! extrude area into rear volume
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< END TRIANGULAR SECTIONS >>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< ADD INTERSECTING VOLUMES >>>>>>>>>>>>>>>>>>>>>>>>>
!
*get,vm,volu,,num,max
!
vadd,vm-2,vm-4
!
btol,.001
vadd,vm,vm-1,35
btol,defa
!
nummrg,all,.001,.001                ! merge all coincident entities
numcmp,all                  ! renumber & compress all entities
!
csys,rib5cs                 ! partition rib triangular section 5
!
k,,0,-3.00,0                    ! keypoints of dividing area
k,,-1.25,-3.00,0
!
*get,km,kp,,num,max
!
a,km,km-1,kp(0,-.511,0),kp(-1.25,-.511,0)   ! create dividing area
!
vsba,30,204,,delete,delete          ! divide volumes
!
numcmp,all
!
csys,rib6cs                 ! partition rib triangular section 6
!
k,,.100,-3.00,0                 ! keypoints of dividing area
k,,1.2,-3.00,0
!
*get,km,kp,,num,max
!
a,km,km-1,kp(.125,-.511,0),kp(1.125,-.511,0)    ! create dividing area
!
vsba,32,208,,delete,delete
!
*get,km,kp,,num,max
lfillt,343,345,.125,km+1        ! create line fillet
r32filcp=km+1
!
csys,0                  ! rib 3/4 fillet c.p. x,y,z corrdinates
r32cpx=kx(r32filcp)
r32cpy=ky(r32filcp)
r32cpz=kz(r32filcp)
kdele,r32filcp
!
csys,rib6cs
!               ! create line fillet
nummrg,all,.001,.001                ! merge entities
a,193,150,152                   ! create filleted area
!
VOFFST,182,.511,,               ! extrude area into lower volumes
VOFFST,191,3.00,,
!
nummrg,all,.001,.001                ! merge all coincident entities
numcmp,all                  ! renumber & compress all entities
!
csn=csn+1
cskp,csn,1,184,195,197,,
!
k,,3.0,330,0
k,,3.0,90,0
k,,3.0,210,0
!
*get,km,kp,,num,max
!
lstr,km-2,km-1
lstr,km-1,km
lstr,km,km-2
!
a,km-2,km-1,km
!
*get,am,area,,num,max
*get,vm,volu,,num,max
!
vsba,vm,am,,delete,delete
!
vdele,37,,,1
!
nummrg,all,.001,.001            ! merge all coincident entities
numcmp,all              ! renumber & compress all entities
!
vsel,s,volu,,34,
aslv,s
cm,aset2,area
vdele,34
vsel,all
asel,all
adele,197,199,1,1
!
asbl,213,367,,delete,keep
!
nummrg,all,.001,.001
!
csys,spher56
!
a,150,151,152
!
vdele,32
vdele,35
!
nummrg,all,.001,.001            ! merge all coincident entities
numcmp,all              ! renumber & compress all entities
!
adele,177,178,1,1           ! delete unnecessary areas
adele,182,191,191-182,1
adele,195,196,1,1
adele,198,205,205-198,1
!
lstr,156,192
!
a,155,156,192,186
a,192,156,161,172
!
asel,s,area,,aset2,
asel,a,area,,177,178,1,
asel,a,area,,197,199,2,
asel,a,area,,207,215,215-207
asel,a,area,,214
cm,aset56,area              ! area set for rib 5/6
!
va,all
!
vsel,all
asel,all
!
nummrg,all,.001,.001            ! merge all coincident entities
numcmp,all              ! renumber & compress all entities
!
LPLOT
!
!******************************** 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
!
x4=abs(r7p4x)
z4=abs(r7p4z)
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< MATH SECTION >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
topw=w2+2*r1                ! top of rib width
!
x1=abs(r7p1x)
z1=abs(r7p1z)
!
x2=abs(r7p2x)
z2=abs(r7p2z)
!
m=(z2-z1)/(x2-x1)           ! slope in XZ plane
!
C=sqrt(1+(m**2))            ! defining a constant
!
DELx54=topw/C
x5=x4+DELx54    !<<<<<<<<<<<<<<<<<<<<<<<< 1st coordinate
!
DELz54=m*DELx54
z5=z4+DELz54    !<<<<<<<<<<<<<<<<<<<<<<<< 2nd coordinate
!
DELx56=r1/C
x6=x5-DELx56    !<<<<<<<<<<<<<<<<<<<<<<<< 3rd coordinate
z6=z5-m*DELx56  !<<<<<<<<<<<<<<<<<<<<<<<< 4th coordinate
!
DELx74=r1/C
x7=x4+DELx74    !<<<<<<<<<<<<<<<<<<<<<<<< 5th coordinate
z7=z4+m*DELx74  !<<<<<<<<<<<<<<<<<<<<<<<< 6th coordinate
!
d2=sqrt(((z7-z1)**2)+((x7-x1)**2))
!
d2prime=sqrt(((z6-z2)**2)+((x6-x2)**2))
!
check2=d2/d2prime
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< END MATH SECTION >>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
*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
!
kdele,km+1          ! delete center keypoint
!
k,km+1,0,0,-d2          ! rib center point in cross-section plane
k,,0,ribh-r1,-d2        ! remaining rib cross-section keypoints
k,,w2/2,0,-d2
k,,w2/2,ribh-r1,-d2
k,,w2/2+r1,ribh,-d2
! 
k,,w2/2+r1,ribh+topthk,-d2
k,,0,ribh+topthk,-d2
!
*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,-d2      ! 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,,,neribw-1,1,
lesize,lm-2,,,netop,1,
lesize,lm-3,,,netop+1,1,
lesize,lm-4,,,netopw-1,1,
!
lesize,lm-5,,,netop,1,
lesize,lm-6,,,netop,1,
lesize,lm-7,,,neribw-1,1,
!
*get,am,area,,num,max
!
aatt,,,2,0          ! set element attributes                
mshkey,0            ! select free meshing
amesh,am-1,am,1         ! mesh areas
!
ESIZE,0,neL7,           ! # element divisions along extrusion
!
type,1              ! set element type
mat,2               ! set material #
esys,0              ! set element c.s.
!
dfext7=d2-.0868         ! calculation of front extrusion
!
VOFFST,am,L2-dfext7,,       ! extrude area into rear volume
VOFFST,am-1,L2-dfext7,,
!
ESIZE,0,1,          ! # element divisions 2nd along extrusion
!
VOFFST,am,-dfext7,,     ! extrude area into front volume
VOFFST,am-1,-dfext7,,
!
csys,rib7cs
!
asel,s,loc,z,-.0868,-.0868,,    ! select rear areas for extrusion
*get,a1,area,,num,max
*get,a2,area,,num,min
!
ESIZE,0,1,          ! # element divisions 3rd extrusion
!
VOFFST,a1,-.0868,,      ! extrude areas into 3rd volumes
VOFFST,a2,-.0868,,
!
asel,all
!
aclear,am-1,am,1        ! clear area meshes
!
csys,rib7cs         ! mirror volumes to create complete rib
!
vsel,s,loc,x,0,w2/2+r1,,    ! select volumes to be mirrored
vsel,r,loc,z,0,-(L2+.0868),,
!
vsymm,x,all,,,,0,0
!
vsel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
LPLOT
!
!<<<<<<<<<<<<<<<<<<<<<<< create rib 7 triangular section >>>>>>>>>>>>>>>>>>>
!
csys,rib7cs
!
k,,-w2/2,-L2*tan(20*rconv),-(L2+.0868)          ! section keypoint
!
*get,pt5,kp,,num,max
!
l,pt5,kp(-w2/2,0,-(L2-dfext7))          ! section lines
l,pt5,kp(-w2/2,0,dfext7)
!
*get,lm,line,,num,max
!
lesize,lm,,,neL2+1,1,               ! size element lines
lesize,lm-1,,,neL2+1,1,
!
a,pt5,kp(-w2/2,0,-L2),kp(-w2/2,0,-dfext7),kp(-w2/2,0,-.0868)
!
lsel,s,loc,x,-w2/2,-w2/2,,
lsel,r,loc,z,-.0868,-(L2+.0868),,
lsel,r,loc,y,0,0,,
!
lccat,all           ! concatenate lines
!
*get,am,area,,num,max
!
aatt,,,2,0          ! set element attributes                
mshkey,1            ! select mapped meshing
amesh,am,am,            ! mesh areas
!
*get,lm,line,,num,max       ! delete concatenated lines
ldele,lm
lsel,all
!
ESIZE,0,2*(neribw-1),       ! # element divisions along extrusion
!
type,1              ! set element type
mat,2               ! set material #
esys,0              ! set element c.s.
!
VOFFST,am,w2,,          ! extrude areas into volume
!
aclear,am,am,           ! clear area meshes
!
!****************************** End 7th Rib Layout  **************************
!
!
!******************* Prepare Ribs for Center Boss Construction ***************
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Ribs 2,3 & 4 >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
vdele,13,16,1,1         ! delete unused volumes
vdele,19,22,1,1
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
LPLOT
!
L=distkp(96,127)
w=distkp(61,110)        ! set area parameters for partitioning
dw1=distkp(123,110)
dw2=distkp(123,117)
!
csn=csn+1
cs1=csn             ! 1st c.s. for keypoints of area
cskp,cs1,0,61,110,60,,
!
k,,0,-1.0,dw1+0.5*dw1
k,,0,-1.0,-0.1
!
*get,km,kp,,num,max     ! create 1st line section
l,km-1,km
!
csn=csn+1
cs2=csn             ! 2nd c.s. for keypoints of area
cskp,cs2,0,117,96,123,,
!
k,,w,dw2+0.5*dw2,1.00
k,,w,-0.1,1.00
!
*get,km,kp,,num,max     ! create 2nd line section
l,km-1,km
!
*get,lm,line,,num,max
!
l,km,km-2
lsbl,lm-1,lm,,delete,delete ! create line intersection
!
l,km+1,km           ! recreate line
ldele,lm+1,,,1          ! delete unused lines
ldele,411,,,1
!
numcmp,all          ! renumber & compress all entities
!
a,107,88,87,94,95,122       ! create partition area
!
*get,am,area,,num,max
!
asel,s,area,,aset34
vsla,s,1
asel,all
!
vsba,all,am,,delete,delete  ! partition volume
!
vsel,all
!
numcmp,all          ! renumber & compress all entities
!
csys,cs2            ! remaining keypoints for partitioning volume
k,,-1.0,0,1.0
csys,cs1
k,,1.0,-1.0,0
!
*get,km,kp,,num,max
a,km,km-1,km-3,km-2,km-4    ! create area
!
*get,am,area,,num,max       ! extrude area into partitioning volume
VOFFST,am,-(L+1.1),,
!
*get,vm,volu,,num,max
vsbv,vm-1,vm            ! subtract volume
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
LPLOT
!
/auto,1
/rep
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Ribs 5,6 & 7 >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
ldele,389,,,1
vclear,27,30,1
vclear,33,36,1
!
vdele,27,30,1,1         ! delete unused volumes
vdele,33,36,1,1
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
L=distkp(159,160)       ! set area parameters for partitioning
w=distkp(190,147)
dw1=distkp(147,165)
dw2=distkp(165,160)
!
csn=csn+1
cs1=csn             ! 1st c.s. for keypoints of area
cskp,cs1,0,190,147,191,,
!
k,,0,-1.0,-(dw1+0.5*dw1)
k,,0,-1.0,0.1
!
*get,km,kp,,num,max     ! create 1st line section
l,km-1,km
!
csn=csn+1
cs2=csn             ! 2nd c.s. for keypoints of area
cskp,cs2,0,160,130,165,,
!
k,,w,dw2+0.5*dw2,-1.00
k,,w,-0.1,-1.00
!
*get,km,kp,,num,max     ! create 2nd line section
l,km-1,km
!
*get,lm,line,,num,max
!
l,km,km-2
lsbl,lm-1,lm,,delete,delete ! create line intersection
!
l,km+1,km           ! recreate line
ldele,lm+1,,,1          ! delete unused lines
ldele,365,,,1
!
numcmp,all          ! renumber & compress all entities
!
a,133,132,128,129,159,150   ! create partition area
!
*get,am,area,,num,max
!
asel,s,area,,aset56
vsla,s,1
!
asel,all
!
btol,.001
vsba,all,am,,delete,delete  ! partition volume
btol,defa
!
vsel,all
!
numcmp,all          ! renumber & compress all entities
!
csys,cs2            ! remaining keypoints for partitioning volume
k,,-1.0,0,-1.0
csys,cs1
k,,1.0,-1.0,0
!
*get,km,kp,,num,max
a,km,km-1,km-3,km-2,km-4    ! create area
!
*get,am,area,,num,max       ! extrude area into partitioning volume
VOFFST,am,(L+1.1),,
!
btol,.001
vsbv,32,33          ! subtract volume
btol,defa
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
LPLOT
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< Fillet / Ribs 2 & 3 >>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
csys,gcs
!
k,,kx(60),ky(60),kz(60)-.25
!
*get,filcp,kp,,num,max
!
csn=csn+1
filcs1=csn
cskp,filcs1,1,filcp,59,60
!
csn=csn+1
filcs2=csn
cskp,filcs2,1,filcp,193,60
!
csys,filcs1
!
k,,r1,0,0
k,,r1,90,0
!
*get,km,kp,,num,max
!
l,km,km-1
*get,fillin1,line,,num,max
!
csys,filcs2
!
k,,r1,0,0
!
l,km,km+1
*get,fillin2,line,,num,max
!
csys,filcs1
!
k,,r1,0,-2.00
l,km+2,km-1
!
csys,filcs2
!
k,,r1,0,2.00
l,km+3,km+1
!
*get,lm,line,,num,max
!
adrag,fillin1,,,,,,lm-1,        ! line drag into areas
adrag,fillin2,,,,,,lm,
!
*get,am,area,,num,max
!
aptn,am-1,am
!
adele,210,211,1,1
ldele,368,369,,1
!
angfil=anglek(215,km-1,km+1)/2      ! angle of fillet (/2)
!
d=(r1-.1465)/tan(angfil)        ! distance of cutting plane from fillet
!
numcmp,all              ! renumber entities
!
csys,filcs1
!
k,,1.0,0,-d
k,,1.0,90,-d
k,,1.0,180,-d
k,,1.0,270,-d
!
*get,km,kp,,num,max
!
csys,0
a,km,km-1,km-2,km-3
!
*get,am,area,,num,max
!
asba,am-2,am,,delete,delete
!
numcmp,all              ! renumber entities
!
csys,filcs2
!
k,,1.0,0,d
k,,1.0,90,d
k,,1.0,180,d
k,,1.0,270,d
!
*get,km,kp,,num,max
!
csys,0
a,km,km-1,km-2,km-3
!
*get,am,area,,num,max
!
asba,am-3,am,,delete,delete
!
adele,207,210,3,1
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,0
!
k,,kx(206),ky(207),kz(206)  ! create center keypoint for rotation
*get,filcp,kp,,num,max
!
rotang=(anglek(filcp,205,207))/rconv    ! rotation angle
arotat,366,,,,,,filcp,206,rotang,2  ! rotate line into area
!
LPLOT
!
/auto,1
/rep
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< Fillet / Ribs 6 & 7 >>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
csys,gcs
!
k,,kx(186),ky(186),kz(186)-.25
!
*get,filcp,kp,,num,max
!
csn=csn+1
filcs1=csn
cskp,filcs1,1,filcp,201,186
!
csn=csn+1
filcs2=csn
cskp,filcs2,1,filcp,181,186
!
csys,filcs1
!
k,,r1,0,0
k,,r1,90,0
!
*get,km,kp,,num,max
!
l,km,km-1
*get,fillin1,line,,num,max
!
csys,filcs2
!
k,,r1,0,0
!
l,km,km+1
*get,fillin2,line,,num,max
!
csys,filcs1
!
k,,r1,0,-2.00
l,km+2,km-1
!
csys,filcs2
!
k,,r1,0,2.00
l,km+3,km+1
!
*get,lm,line,,num,max
!
adrag,fillin1,,,,,,lm-1,        ! line drag into areas
adrag,fillin2,,,,,,lm,
!
*get,am,area,,num,max
!
aptn,am-1,am
!
adele,214,215,1,1
ldele,377,378,,1
!
angfil=anglek(221,km-1,km+1)/2      ! angle of fillet (/2)
d=(r1-.1306)/tan(angfil)        ! distance of cutting plane from fillet
!
numcmp,all              ! renumber entities
!
csys,filcs1
!
k,,1.0,0,-d
k,,1.0,90,-d
k,,1.0,180,-d
k,,1.0,270,-d
!
*get,km,kp,,num,max
!
csys,0
a,km,km-1,km-2,km-3
!
*get,am,area,,num,max
!
asba,am-2,am,,delete,delete
!
numcmp,all              ! renumber entities
!
csys,filcs2
!
k,,1.0,0,d
k,,1.0,90,d
k,,1.0,180,d
k,,1.0,270,d
!
*get,km,kp,,num,max
!
csys,0
a,km,km-1,km-2,km-3
!
*get,am,area,,num,max
!
asba,am-3,am,,delete,delete
!
adele,211,214,3,1
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,0
!
k,,kx(212),ky(211),kz(212)  ! create center keypoint for rotation
*get,filcp,kp,,num,max
!
rotang=(anglek(filcp,211,213))/rconv    ! rotation angle
arotat,374,,,,,,filcp,212,rotang,2  ! rotate line into area
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< Fillet Flat Areas / 2 & 3 >>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
csys,0
!
k,,kx(205),ky(71),kz(205)
!
*get,km,kp,,num,max
!
a,58,71,km,205,59
!
k,,kx(207),ky(107),kz(207)
!
a,192,107,km+1,207,193
!
adrag,372,,,,,,383,
adrag,373,,,,,,386,
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<< Fillet Flat Areas / 6 & 7 >>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
csys,0
!
k,,kx(211),ky(147),kz(211)
!
*get,km,kp,,num,max
!
a,200,147,km,211,201
!
k,,kx(213),ky(189),kz(213)
!
a,182,189,km+1,213
a,213,181,182
!
adrag,380,,,,,,395,
adrag,381,,,,,,397,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!********************************** Front Rib Fillets *********************************
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Rib 3/4 Fillet >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
csys,rib4cs         ! create fillet profile
!
k,,1.250+.125,-(.511+r1),3.00   ! profile keypoints
k,,1.00+.125,-(.511+r1),3.00
k,,1.250+.125,-(.511),3.00
k,,1.250+.125+r1,-(.511+r1),3.00
k,,1.250+.125+r1,-(.511+r1+3.00),3.00
!
k,,1.00+.125,-(.511+r1+3.00),3.00
!
*get,km,kp,,num,max
!
l,km,km-4           ! profile 'straight' lines
l,km,km-1
l,km-1,km-2
!
csn=csn+1
FLfilcs=csn         ! fillet profile c.s.
cskp,FLfilcs,1,km-5,km-4,km-3,,
!
l,km-4,km-3         ! fillet 'round' lines
l,km-3,km-2
!
a,km-4,km-3,km-2,km-1,km    ! profile area
!
*get,am,area,,num,max
!
VOFFST,am,6.00,         ! create rib fillet profile volume
!
*get,vm,volu,,num,max
!
vsbv,29,vm,,delete,keep     ! subtract profile volume from rib 4 volumes
vsbv,19,vm,,delete,delete
!
vsel,s,volu,,34,  !------------------- Rib 3/4 Fillet Modification ---------------------
aslv,s
lsla,s
!
cm,asetv34,area     ! create area set for volume
cm,lsetv34,line     ! create area set for lines
vdele,34
!
adele,120,,,1       ! delete areas for modifying volume
adele,105,,,1
adele,121,,,1
!
numcmp,all
!
csys,0
!
k,,r34cpx,ky(118),r34cpz
*get,r34filcp,kp,,num,max
!
csn=csn+1
r34filcs=csn
cskp,r34filcs,1,r34filcp,118,106
csys,r34filcs
!
l,118,120
!
a,94,87,118,120
a,106,120,118
!
va,all
!
!------------------------------ End Volume modification -----------------------------
!
allsel,all
vadd,32,29          ! add rib 4 top & bottom volumes
!
numcmp,all          ! renumber entities
!
LPLOT
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Rib 5/6 Fillet >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!< Clear volume mesh to clear line mesh attached to volume to be modified
!
vsel,s,volu,,26,        ! clear volume mesh
vclear,26
vdele,26,,,1
vsel,all
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< End volume clear
!
btol,.001
vadd,29,30
btol,defa
!
numcmp,all          ! renumber entities
!
csys,rib5cs         ! create fillet profile
!
k,,-1.250,-(.511+r1),-2.00  ! profile keypoints
k,,-1.00,-(.511+r1),-2.00
k,,-1.250,-(.511),-2.00
k,,-(1.250+r1),-(.511+r1),-2.00
k,,-(1.250+r1),-(.511+r1+3.00),-2.00
!
k,,-1.00,-(.511+r1+3.00),-2.00
!
*get,km,kp,,num,max
!
l,km,km-4           ! profile 'straight' lines
l,km,km-1
l,km-1,km-2
!
csn=csn+1
FRfilcs=csn             ! fillet profile c.s.
cskp,FRfilcs,1,km-5,km-4,km-3,,
!
l,km-4,km-3         ! fillet 'round' lines
l,km-3,km-2
!
a,km-4,km-3,km-2,km-1,km    ! profile area
!
*get,am,area,,num,max
!
VOFFST,am,6.00,         ! create rib fillet profile volume
!
*get,vm,volu,,num,max
!
vsbv,26,vm,,delete,keep     ! subtract profile volume from rib 4 volumes
vsbv,21,vm,,delete,delete
!
csys,rib7cs
vsymm,x,24,,,,0,0
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
vsel,s,volu,,30,  !------------------- Rib 5/6 Fillet Modification ---------------------
aslv,s
lsla,s
!
cm,asetv30,area     ! create area set for volume
cm,lsetv30,line     ! create area set for lines
vdele,30
!
adele,148,,,1       ! delete areas for modifying volume
adele,147,,,1
adele,134,,,1
adele,146,,,1
!
numcmp,all
!
csys,0
!
k,,r32cpx,ky(158),r32cpz
*get,r32filcp,kp,,num,max
!
csn=csn+1
r32filcs=csn
cskp,r32filcs,1,r32filcp,158,151
csys,r32filcs
!
l,160,158
lstr,151,160
lstr,151,144
!
a,129,125,158,160
a,158,160,151
a,160,151,144,161
a,151,152,143,144
!
cm,r56aset,area
va,all
*get,r56volu,volu,,num,max
!
!---------------------------- End Volume modification -----------------------------
!
allsel,all
!
lplot
!
!************************** Center Section Construction ************************
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Section Areas >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Rib 1, Area 1 >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!****************************** 1st 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)
theta=theta1
!
k,,rcb1,theta,0
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,rcb1,theta,-(ribh+topthk)
k,,rbout,theta,-(ribh+topthk)
!   
k,,rbout,theta,-(topthk+r1)
k,,p1r,theta,-topthk
k,,p1r,theta,0
k,,rcb1,theta,-(topthk+r1)
k,,rcb1,theta,-cbore2
!
k,,p1r,theta,-(topthk+r1)   ! center point of rectangular & spherical cs
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,cbr,theta,-(topthk+ribh)
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn
cskp,plncs,1,km-3,km-8,km-7
!
lstr,km-6,km-18         ! create straight lines     
lstr,km-18,km-17
lstr,km-17,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
lstr,km,km-10
!
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-7,km-6
lstr,km-8,km-18
lstr,km-17,km-5
!
lstr,km-8,km-5
lstr,km-15,km-5
lstr,km-8,km-4
lstr,km-4,km-5
lstr,km-4,km-14
!
lstr,km-10,km-4
lstr,km-11,km-14
lstr,km-1,km-12
lstr,km-1,km-2
lstr,km-2,km
!
l,km-8,km-7         ! creation of fillet line 1
!
a,km-12,km-11,km-14,km-13   ! creation of areas
a,km-11,km,km-10,km-4,km-14
a,km-14,km-4,km-5,km-15
a,km-15,km-5,km-17,km-16
a,km-10,km-9,km-8,km-4
!       
a,km-4,km-8,km-5
a,km-5,km-8,km-18,km-17
a,km-8,km-7,km-6,km-18
a,km-1,km-2,km,km-11,km-12
!
csn=csn+1           ! creation of cylindrical cs
CS1=csn
cskp,CS1,1,km-3,km-8,p3
!
csn=csn+1           ! creation of spherical cs
CSr2=csn
cskp,CSr2,2,km-3,km-8,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-25,,,2,1,
lesize,lm-24,,,1,1,
lesize,lm-23,,,2,1,
lesize,lm-22,,,2,1,
lesize,lm-21,,,2,1,
!
lesize,lm-20,,,2,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,1,1,
lesize,lm-16,,,1,1,
!
lesize,lm-15,,,netop,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,netop,1,
lesize,lm-12,,,netop,1,
lesize,lm-11,,,2,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,,,1,1,
lesize,lm-2,,,3,1,
lesize,lm-1,,,1,1,
!
lesize,lm,,,3,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-8           ! creation of cylindrical lines
l,p4,km-9
!
a,p3,km-8,km-9,p4       ! creation of cylindrical area
!
lesize,lm+1,,,3,1,      ! sizing of cylindrical lines
lesize,lm+2,,,3,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,plncs
!
sp1=kp(r1,0,0)          ! spherical point 1
sp2=kp(r1,90,0)         ! spherical point 2
!
csys,rib1cs
!
sp3=kp(-w1/2,ribh-r1,0)     ! spherical point 3
!
csys,CSr2           ! spherical coordinate system
!
a,sp1,sp2,sp3           ! creation of spherical area
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!****************************** 1st 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)
theta=theta2
!
k,,rcb1,theta,0
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,rcb1,theta,-(ribh+topthk)
k,,rbout,theta,-(ribh+topthk)
!   
k,,rbout,theta,-(topthk+r1)
k,,p1r,theta,-topthk
k,,p1r,theta,0
k,,rcb1,theta,-(topthk+r1)
k,,rcb1,theta,-cbore2
!
k,,p1r,theta,-(topthk+r1)   ! center point of rectangular & spherical cs
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,cbr,theta,-(topthk+ribh)
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn
cskp,plncs,1,km-3,km-8,km-7
!
lstr,km-6,km-18         ! create straight lines     
lstr,km-18,km-17
lstr,km-17,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
lstr,km,km-10
!
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-7,km-6
lstr,km-8,km-18
lstr,km-17,km-5
!
lstr,km-8,km-5
lstr,km-15,km-5
lstr,km-8,km-4
lstr,km-4,km-5
lstr,km-4,km-14
!
lstr,km-10,km-4
lstr,km-11,km-14
lstr,km-1,km-12
lstr,km-1,km-2
lstr,km-2,km
!
l,km-8,km-7         ! creation of fillet line 1
!
a,km-12,km-11,km-14,km-13   ! creation of areas
a,km-11,km,km-10,km-4,km-14
a,km-14,km-4,km-5,km-15
a,km-15,km-5,km-17,km-16
a,km-10,km-9,km-8,km-4
!       
a,km-4,km-8,km-5
a,km-5,km-8,km-18,km-17
a,km-8,km-7,km-6,km-18
a,km-1,km-2,km,km-11,km-12
!
csn=csn+1           ! creation of cylindrical cs
CS1=csn
cskp,CS1,1,km-3,km-8,p3
!
csn=csn+1           ! creation of spherical cs
CSr2=csn
cskp,CSr2,2,km-3,km-8,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-25,,,2,1,
lesize,lm-24,,,1,1,
lesize,lm-23,,,2,1,
lesize,lm-22,,,2,1,
lesize,lm-21,,,2,1,
!
lesize,lm-20,,,2,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,1,1,
lesize,lm-16,,,1,1,
!
lesize,lm-15,,,netop,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,netop,1,
lesize,lm-12,,,netop,1,
lesize,lm-11,,,2,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,,,1,1,
lesize,lm-2,,,3,1,
lesize,lm-1,,,1,1,
!
lesize,lm,,,3,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-8           ! creation of cylindrical lines
l,p4,km-9
!
a,p3,km-8,km-9,p4       ! creation of cylindrical area
!
lesize,lm+1,,,3,1,      ! sizing of cylindrical lines
lesize,lm+2,,,3,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,plncs
!
sp1=kp(r1,0,0)          ! spherical point 1
sp2=kp(r1,90,0)         ! spherical point 2
!
csys,rib1cs
!
sp3=kp(w1/2,ribh-r1,0)      ! spherical point 3
!
csys,CSr2           ! spherical coordinate system
!
a,sp1,sp2,sp3           ! creation of spherical area
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!<<<<<<<<<<<<<<<<<<< Extrude rib 2 lines inward to intersect center section
!
csys,rib2cs
! 
k,,-(w2/2+r1),(topthk+ribh),2.00
*get,km,kp,,num,max
l,km,kp(-(w2/2+r1),(topthk+ribh),-d1)
! 
csys,gcs
!
k,,p1r,45,0
l,km+1,kp(p1r,theta2,0)
*get,lm,line,,num,max
lsbl,lm-1,lm,,delete,delete
ldele,462,,,1
!
*get,Ltot,line,461,leng
*get,Lext1,line,136,leng
!
Lext2=Ltot-Lext1
adrag,121,130,129,,,,136,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib2cs
!
k,,-(w2/2),0,kz(265)
!
*get,km,kp,,num,max
!
lstr,80,km
!
*get,lm,line,,num,max
!
adrag,462,464,459,,,,lm,
ldele,460,,,1
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< End Extrusion >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!
!****************************** 2nd Rib, AREA 3 ****************************
!
csys,rib2cs         ! rib 1 coordinate system
!
p1=kp(-(w2/2+r1),ribh+topthk,kz(269))   ! keypoints from position on rib 
p2=kp(-(w2/2+r1),ribh,kz(269))
p3=kp(-w2/2,ribh-r1,kz(269))
p4=kp(-w2/2,0,kz(269))
!
csys,gcs            ! global coordinate system
!
p1r=kx(p1)
theta3=ky(p1)
theta=theta3
!
k,,rcb1,theta,0
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,rcb1,theta,-(ribh+topthk)
k,,rbout,theta,-(ribh+topthk)
!   
k,,rbout,theta,-(topthk+r1)
k,,p1r,theta,-topthk
k,,p1r,theta,0
k,,rcb1,theta,-(topthk+r1)
k,,rcb1,theta,-cbore2
!
k,,p1r,theta,-(topthk+r1)   ! center point of rectangular & spherical cs
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,cbr,theta,-(topthk+ribh)
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn
cskp,plncs,1,km-3,km-8,km-7
!
lstr,km-6,km-18         ! create straight lines     
lstr,km-18,km-17
lstr,km-17,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
lstr,km,km-10
!
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-7,km-6
lstr,km-8,km-18
lstr,km-17,km-5
!
lstr,km-8,km-5
lstr,km-15,km-5
lstr,km-8,km-4
lstr,km-4,km-5
lstr,km-4,km-14
!
lstr,km-10,km-4
lstr,km-11,km-14
lstr,km-1,km-12
lstr,km-1,km-2
lstr,km-2,km
!
l,km-8,km-7         ! creation of fillet line 1
!
a,km-12,km-11,km-14,km-13   ! creation of areas
a,km-11,km,km-10,km-4,km-14
a,km-14,km-4,km-5,km-15
a,km-15,km-5,km-17,km-16
a,km-10,km-9,km-8,km-4
!       
a,km-4,km-8,km-5
a,km-5,km-8,km-18,km-17
a,km-8,km-7,km-6,km-18
a,km-1,km-2,km,km-11,km-12
!
csn=csn+1           ! creation of cylindrical cs
CS1=csn
cskp,CS1,1,km-3,km-8,p3
!
csn=csn+1           ! creation of spherical cs
CSr2=csn
cskp,CSr2,2,km-3,km-8,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-25,,,2,1,
lesize,lm-24,,,1,1,
lesize,lm-23,,,2,1,
lesize,lm-22,,,2,1,
lesize,lm-21,,,2,1,
!
lesize,lm-20,,,2,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,1,1,
lesize,lm-16,,,1,1,
!
lesize,lm-15,,,netop,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,netop,1,
lesize,lm-12,,,netop,1,
lesize,lm-11,,,2,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,,,1,1,
lesize,lm-2,,,3,1,
lesize,lm-1,,,1,1,
!
lesize,lm,,,3,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-8           ! creation of cylindrical lines
l,p4,km-9
!
a,p3,km-8,km-9,p4       ! creation of cylindrical area
!
lesize,lm+1,,,3,1,      ! sizing of cylindrical lines
lesize,lm+2,,,3,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,plncs
!
sp1=kp(r1,0,0)          ! spherical point 1
sp2=kp(r1,90,0)         ! spherical point 2
!
csys,rib2cs
!
sp3=kp(-w2/2,ribh-r1,0)     ! spherical point 3
!
csys,CSr2           ! spherical coordinate system
!
a,sp1,sp2,sp3           ! creation of spherical area
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!<<<<<<<<<<<<<<<<<<< Extrude rib 7 lines inward to intersect center section
!
csys,rib7cs
! 
k,,(w2/2+r1),(topthk+ribh),2.00
*get,km,kp,,num,max
l,km,kp((w2/2+r1),(topthk+ribh),-d2)
!
!---------------------- Create point p1 --------------------
csys,rib1cs
p1=kp(-(w1/2+r1),ribh+topthk,0)
csys,gcs
p1r=kx(p1)
!-----------------------------------------------------------
! 
csys,gcs
!
k,,p1r,150,0
l,km+1,kp(p1r,theta1,0)
*get,lm,line,,num,max
lsbl,lm-1,lm,,delete,delete
ldele,502,,,1
!
*get,Ltot,line,501,leng
*get,Lext1,line,294,leng
!
Lext2=Ltot-Lext1
adrag,273,274,278,,,,294,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,rib7cs
!
k,,(w2/2),0,kz(289)
!
*get,km,kp,,num,max
!
l,176,km
!
*get,lm,line,,num,max
!
adrag,502,504,499,,,,lm,
ldele,500,,,1
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<< End Extrusion >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
!****************************** 7th Rib, AREA 10 ******************************
!
csys,rib7cs         ! 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            ! model global coordinate system
!
p1r=kx(p1)
theta10=ky(p1)
theta=theta10
!
k,,rcb1,theta,0
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,rcb1,theta,-(ribh+topthk)
k,,rbout,theta,-(ribh+topthk)
!   
k,,rbout,theta,-(topthk+r1)
k,,p1r,theta,-topthk
k,,p1r,theta,0
k,,rcb1,theta,-(topthk+r1)
k,,rcb1,theta,-cbore2
!
k,,p1r,theta,-(topthk+r1)   ! center point of rectangular & spherical cs
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,cbr,theta,-(topthk+ribh)
!
*get,km,kp,,num,max     ! creation of planar cylindrical cs 
csn=csn+1
plncs=csn
cskp,plncs,1,km-3,km-8,km-7
!
lstr,km-6,km-18         ! create straight lines     
lstr,km-18,km-17
lstr,km-17,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
lstr,km,km-10
!
lstr,km-10,km-9
lstr,km-9,km-8
lstr,km-7,km-6
lstr,km-8,km-18
lstr,km-17,km-5
!
lstr,km-8,km-5
lstr,km-15,km-5
lstr,km-8,km-4
lstr,km-4,km-5
lstr,km-4,km-14
!
lstr,km-10,km-4
lstr,km-11,km-14
lstr,km-1,km-12
lstr,km-1,km-2
lstr,km-2,km
!
l,km-8,km-7         ! creation of fillet line 1
!
a,km-12,km-11,km-14,km-13   ! creation of areas
a,km-11,km,km-10,km-4,km-14
a,km-14,km-4,km-5,km-15
a,km-15,km-5,km-17,km-16
a,km-10,km-9,km-8,km-4
!       
a,km-4,km-8,km-5
a,km-5,km-8,km-18,km-17
a,km-8,km-7,km-6,km-18
a,km-1,km-2,km,km-11,km-12
!
csn=csn+1           ! creation of cylindrical cs
CS1=csn
cskp,CS1,1,km-3,km-8,p3
!
csn=csn+1           ! creation of spherical cs
CSr2=csn
cskp,CSr2,2,km-3,km-8,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-25,,,2,1,
lesize,lm-24,,,1,1,
lesize,lm-23,,,2,1,
lesize,lm-22,,,2,1,
lesize,lm-21,,,2,1,
!
lesize,lm-20,,,2,1,
lesize,lm-19,,,2,1,
lesize,lm-18,,,2,1,
lesize,lm-17,,,1,1,
lesize,lm-16,,,1,1,
!
lesize,lm-15,,,netop,1,
lesize,lm-14,,,2,1,
lesize,lm-13,,,netop,1,
lesize,lm-12,,,netop,1,
lesize,lm-11,,,2,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,,,1,1,
lesize,lm-2,,,3,1,
lesize,lm-1,,,1,1,
!
lesize,lm,,,3,1,
!
csys,CS1            ! cylindrical cs at center point
!
l,p3,km-8           ! creation of cylindrical lines
l,p4,km-9
!
a,p3,km-8,km-9,p4       ! creation of cylindrical area
!
lesize,lm+1,,,3,1,      ! sizing of cylindrical lines
lesize,lm+2,,,3,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,plncs
!
sp1=kp(r1,0,0)          ! spherical point 1
sp2=kp(r1,90,0)         ! spherical point 2
!
csys,rib7cs
!
sp3=kp(w2/2,ribh-r1,0)      ! spherical point 3
!
csys,CSr2           ! spherical coordinate system
!
a,sp1,sp2,sp3           ! creation of spherical area
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
!************************** 3rd / 4th Ribs, AREA 4 ****************************
!
csys,gcs            ! model global coordinate system
!
phi1=ky(61)
phi2=ky(60)
phi3=ky(199)
phi4=ky(202)
phi5=ky(211)
!
p1=61
p2=60
p3=199
p4=202
p5=211
!
theta4=(phi1+phi2+phi3+phi4+phi5)/5
!
theta=theta4
!
k,,rcb1,theta,0         ! area profile keypoints 
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,rcb1,theta,-(ribh+topthk)
k,,rcb1,theta,-(topthk+r1)
!
k,,rcb1,theta,-cbore2
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,cbr,theta,-(topthk+ribh)
!
*get,km,kp,,num,max
!
lstr,p1,km-13           ! create straight lines
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-6,km
lstr,km,km-5
!
lstr,km-5,p5
lstr,km-13,p3
lstr,p4,km-12
lstr,km-4,km-12
lstr,p4,km-4
!
lstr,km-10,km-4
lstr,p4,km-3
lstr,km-4,km-3
lstr,km-3,km-9
lstr,km-3,km-5
!
lstr,km-6,km-9
lstr,km-7,km-1
lstr,km-1,km-2
lstr,km-2,km
!
a,km-7,km-6,km-9,km-8       ! creation of areas
a,km-6,km,km-5,km-3,km-9
a,km-9,km-3,km-4,km-10
a,km-10,km-4,km-12,km-11
a,km-5,p5,p4,km-3
!       
a,km-3,p4,km-4
a,km-4,p4,km-12
a,p4,p3,km-13,km-12
a,p3,p2,p1,km-13
a,km-7,km-1,km-2,km,km-6
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-23,,,3,1,
lesize,lm-22,,,1,1,
lesize,lm-21,,,2,1,
lesize,lm-20,,,2,1,
lesize,lm-19,,,2,1,
!
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
lesize,lm-16,,,2,1,
lesize,lm-15,,,1,1,
lesize,lm-14,,,1,1,
!
lesize,lm-13,,,netop,1,
lesize,lm-12,,,netop,1,
lesize,lm-11,,,netop,1,
lesize,lm-10,,,2,1,
lesize,lm-9,,,netop,1,
!
lesize,lm-8,,,2,1,
lesize,lm-7,,,netop,1,
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
!
lesize,lm-3,,,2,1,
lesize,lm-2,,,1,1,
lesize,lm-1,,,3,1,
lesize,lm,,,1,1,
!
lesize,344,,,3,1,       ! sizing of pre-built fillet lines
lesize,347,,,2,1,
lesize,364,,,2,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
!************************** 5th / 6th Ribs, AREA 9 ****************************
!
csys,gcs            ! model global coordinate system
!
phi1=ky(180)
phi2=ky(181)
phi3=ky(205)
phi4=ky(208)
phi5=ky(214)
!
p1=180
p2=181
p3=205
p4=208
p5=214
!
theta9=(phi1+phi2+phi3+phi4+phi5)/5
!
theta=theta9
!
k,,rcb1,theta,0         ! area profile keypoints 
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,rcb1,theta,-(ribh+topthk)
k,,rcb1,theta,-(topthk+r1)
!
k,,rcb1,theta,-cbore2
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,cbr,theta,-(topthk+ribh)
!
*get,km,kp,,num,max
!
lstr,p1,km-13           ! create straight lines
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-6,km
lstr,km,km-5
!
lstr,km-5,p5
lstr,km-13,p3
lstr,p4,km-12
lstr,km-4,km-12
lstr,p4,km-4
!
lstr,km-10,km-4
lstr,p4,km-3
lstr,km-4,km-3
lstr,km-3,km-9
lstr,km-3,km-5
!
lstr,km-6,km-9
lstr,km-7,km-1
lstr,km-1,km-2
lstr,km-2,km
!
a,km-7,km-6,km-9,km-8       ! creation of areas
a,km-6,km,km-5,km-3,km-9
a,km-9,km-3,km-4,km-10
a,km-10,km-4,km-12,km-11
a,km-5,p5,p4,km-3
!       
a,km-3,p4,km-4
a,km-4,p4,km-12
a,p4,p3,km-13,km-12
a,p3,p2,p1,km-13
a,km-7,km-1,km-2,km,km-6
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-23,,,3,1,
lesize,lm-22,,,1,1,
lesize,lm-21,,,2,1,
lesize,lm-20,,,2,1,
lesize,lm-19,,,2,1,
!
lesize,lm-18,,,2,1,
lesize,lm-17,,,2,1,
lesize,lm-16,,,2,1,
lesize,lm-15,,,1,1,
lesize,lm-14,,,1,1,
!
lesize,lm-13,,,netop,1,
lesize,lm-12,,,netop,1,
lesize,lm-11,,,netop,1,
lesize,lm-10,,,2,1,
lesize,lm-9,,,netop,1,
!
lesize,lm-8,,,2,1,
lesize,lm-7,,,netop,1,
lesize,lm-6,,,2,1,
lesize,lm-5,,,2,1,
lesize,lm-4,,,2,1,
!
lesize,lm-3,,,2,1,
lesize,lm-2,,,1,1,
lesize,lm-1,,,3,1,
lesize,lm,,,1,1,
!
lesize,352,,,3,1,       ! sizing of pre-built fillet lines
lesize,355,,,2,1,
lesize,373,,,2,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
!****************************** 4th Rib, AREA 6 ****************************
!
csys,gcs            ! model global coordinate system
!
p1=188
p2=187
p3=117
!
theta6=ky(p3)
theta=theta6
!
k,,rcb1,theta,0         ! area profile keypoints 
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,cbr,theta,-(topthk+ribh)
k,,rcb1,theta,-(ribh+topthk)
!
k,,rcb1,theta,-cbore2
k,,rcb1,theta,-(topthk+r1)
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,kx(p3),theta,0
!
*get,km,kp,,num,max
!
lstr,km,km-14           ! create straight lines
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-6,km-5
!
lstr,km-5,km-4
lstr,km-4,km-3
lstr,km-3,km-13
lstr,km-3,km-11
lstr,km-4,km-10
!
lstr,km-10,km-7
lstr,p3,km-5
lstr,km,p3
lstr,km-8,km-1
lstr,km-1,km-2
!
lstr,km-2,km-6
lstr,p1,km
!
a,km-9,km-8,km-7,km-10      ! creation of areas
a,km-10,km-7,km-6,km-5,km-4
a,km-10,km-4,km-3,km-11
a,km-11,km-3,km-13,km-12
a,km-14,km-13,km-3,km-4,km-5,p3,km
!       
a,km,p3,p2,p1
a,km-8,km-1,km-2,km-6,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,3,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,,,1,1,
lesize,lm-12,,,1,1,
!
lesize,lm-11,,,2,1,
lesize,lm-10,,,2,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,2,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,3,1,
lesize,lm-4,,,7,1,
lesize,lm-3,,,1,1,
lesize,lm-2,,,3,1,
!
lesize,lm-1,,,1,1,
lesize,lm,,,1,1,
!
lesize,320,,,7,1,       ! sizing of pre-built fillet lines
lesize,386,,,1,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
!****************************** 5th Rib, AREA 7 ****************************
!
csys,gcs            ! model global coordinate system
!
p1=192
p2=191
p3=153
!
theta7=ky(p3)
theta=theta7
!
k,,rcb1,theta,0         ! area profile keypoints 
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,cbr,theta,-(topthk+ribh)
k,,rcb1,theta,-(ribh+topthk)
!
k,,rcb1,theta,-cbore2
k,,rcb1,theta,-(topthk+r1)
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,kx(p3),theta,0
!
*get,km,kp,,num,max
!
lstr,km,km-14           ! create straight lines
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-6,km-5
!
lstr,km-5,km-4
lstr,km-4,km-3
lstr,km-3,km-13
lstr,km-3,km-11
lstr,km-4,km-10
!
lstr,km-10,km-7
lstr,p3,km-5
lstr,km,p3
lstr,km-8,km-1
lstr,km-1,km-2
!
lstr,km-2,km-6
lstr,p1,km
!
a,km-9,km-8,km-7,km-10      ! creation of areas
a,km-10,km-7,km-6,km-5,km-4
a,km-10,km-4,km-3,km-11
a,km-11,km-3,km-13,km-12
a,km-14,km-13,km-3,km-4,km-5,p3,km
!       
a,km,p3,p2,p1
a,km-8,km-1,km-2,km-6,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,3,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,,,1,1,
lesize,lm-12,,,1,1,
!
lesize,lm-11,,,2,1,
lesize,lm-10,,,2,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,2,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,3,1,
lesize,lm-4,,,7,1,
lesize,lm-3,,,1,1,
lesize,lm-2,,,3,1,
!
lesize,lm-1,,,1,1,
lesize,lm,,,1,1,
!
lesize,328,,,7,1,       ! sizing of pre-built fillet lines
lesize,302,,,1,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!------------- Condition areas for creation of area profiles 5 & 8 -----------------
!
csys,0
!
l,105,189
l,196,143
!
*get,lm,line,,num,max
!
asbl,210,lm-1,,delete,delete
asbl,162,lm,,delete,delete
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
!----------------------------------------------------------------------------------
!
!
!****************************** 3rd/4th Rib, AREA 5 ****************************
!
csys,gcs            ! model global coordinate system
!
p1=190
p2=189
p3=105
!
theta5=ky(p3)
theta=theta5
!
k,,rcb1,theta,0         ! area profile keypoints 
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,cbr,theta,-(topthk+ribh)
k,,rcb1,theta,-(ribh+topthk)
!
k,,rcb1,theta,-cbore2
k,,rcb1,theta,-(topthk+r1)
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,kx(p3),theta,0
!
*get,km,kp,,num,max
!
lstr,km,km-14           ! create straight lines
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-6,km-5
!
lstr,km-5,km-4
lstr,km-4,km-3
lstr,km-3,km-13
lstr,km-3,km-11
lstr,km-4,km-10
!
lstr,km-10,km-7
lstr,p3,km-5
lstr,km,p3
lstr,km-8,km-1
lstr,km-1,km-2
!
lstr,km-2,km-6
lstr,p1,km
!
a,km-9,km-8,km-7,km-10      ! creation of areas
a,km-10,km-7,km-6,km-5,km-4
a,km-10,km-4,km-3,km-11
a,km-11,km-3,km-13,km-12
a,km-14,km-13,km-3,km-4,km-5,p3,km
!       
a,km,p3,p2,p1
a,km-8,km-1,km-2,km-6,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,3,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,,,1,1,
lesize,lm-12,,,1,1,
!
lesize,lm-11,,,2,1,
lesize,lm-10,,,2,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,2,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,3,1,
lesize,lm-4,,,7,1,
lesize,lm-3,,,1,1,
lesize,lm-2,,,3,1,
!
lesize,lm-1,,,1,1,
lesize,lm,,,1,1,
!
lesize,323,,,7,1,       ! sizing of pre-built fillet lines
lesize,631,,,1,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!
!****************************** 5th Rib, AREA 8 ****************************
!
csys,gcs            ! model global coordinate system
!
p1=197
p2=196
p3=143
!
theta8=ky(p3)
theta=theta8
!
k,,rcb1,theta,0         ! area profile keypoints 
k,,rcb1,theta,-cbore1
k,,rcb2,theta,-cbore1
k,,rcb2,theta,-(topthk+r1)
k,,rcb2,theta,-cbore2
! 
k,,rbore,theta,-cbore2
k,,rbore,theta,-(ribh+topthk)
k,,rcb2,theta,-(ribh+topthk)
k,,cbr,theta,-(topthk+ribh)
k,,rcb1,theta,-(ribh+topthk)
!
k,,rcb1,theta,-cbore2
k,,rcb1,theta,-(topthk+r1)
k,,cbr,theta,-(topthk+ribh+bthk)
k,,rbore,theta,-(topthk+ribh+bthk)
k,,kx(p3),theta,0
!
*get,km,kp,,num,max
!
lstr,km,km-14           ! create straight lines
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-6,km-5
!
lstr,km-5,km-4
lstr,km-4,km-3
lstr,km-3,km-13
lstr,km-3,km-11
lstr,km-4,km-10
!
lstr,km-10,km-7
lstr,p3,km-5
lstr,km,p3
lstr,km-8,km-1
lstr,km-1,km-2
!
lstr,km-2,km-6
lstr,p1,km
!
a,km-9,km-8,km-7,km-10      ! creation of areas
a,km-10,km-7,km-6,km-5,km-4
a,km-10,km-4,km-3,km-11
a,km-11,km-3,km-13,km-12
a,km-14,km-13,km-3,km-4,km-5,p3,km
!       
a,km,p3,p2,p1
a,km-8,km-1,km-2,km-6,km-7
!
*get,lm,line,,num,max       ! element sizing of lines
!
lesize,lm-21,,,3,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,,,1,1,
lesize,lm-12,,,1,1,
!
lesize,lm-11,,,2,1,
lesize,lm-10,,,2,1,
lesize,lm-9,,,2,1,
lesize,lm-8,,,2,1,
lesize,lm-7,,,2,1,
!
lesize,lm-6,,,2,1,
lesize,lm-5,,,3,1,
lesize,lm-4,,,7,1,
lesize,lm-3,,,1,1,
lesize,lm-2,,,3,1,
!
lesize,lm-1,,,1,1,
lesize,lm,,,1,1,
!
lesize,339,,,7,1,       ! sizing of pre-built fillet lines
lesize,632,,,1,1,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT               ! plot lines
!
!----------------------- Create Center-Axis Keypoints-----------------------
!
csys,gcs
!
k,,0,0,-1.00
*get,axpt1,kp,,num,min
*get,axpt2,kp,,num,max
lstr,axpt1,axpt2
!
!---------------------------------------------------------------------------
!
!
!<<<<<<<<<<<<<<<<<<<<<<< Mesh & Rotate Area 2 Profile >>>>>>>>>>>>>>>>>>>>>>>
!
csys,gcs
!
theta1st=theta2
theta2nd=theta3
!
asel,s,loc,y,theta1st,theta1st,,   ! select profile @ specified angle
lsla,s
!
aatt,,,2,0          ! set area element attributes
mshkey,1            ! select mapped meshing
!
*get,am,area,,num,max
*get,lm,line,,num,max
!
lsel,r,line,,lm-17,lm-16,,  ! select lines to be concatenated
lccat,all           ! concatenate lines
!
lsla,s              ! reselect line/area set
!
lsel,r,line,,lm-16,lm-15,,  ! select lines to be concatenated
lccat,all           ! concatenate lines
!
lsla,s
!
asel,u,area,,am-2,am-1,,    ! select areas for mapped meshing
!
amesh,all           ! map mesh areas
!
asel,s,loc,y,theta1st,theta1st,,   ! reselect areas for free meshing
asel,u,area,,am-8,am-3,,
asel,u,area,,am,am,,
!
aatt,,,2,0          ! set area element attributes
mshkey,0            ! select free meshing
!
amesh,all           ! free mesh areas
!
asel,all
lsel,all
!
*get,lm,line,,num,max
ldele,lm-1,lm,
!
!
!----------------- Rotate Area Meshes into Volume Meshes ---------------------
!
nerot=2             ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta1st,theta1st,,    ! select area profile @ specified angle
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt1,axpt2,phi,1   ! create volumes
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
aclear,all
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!<<<<<<<<<<<<<<<<<<<<<<< Mesh & Rotate Area 10 Profile >>>>>>>>>>>>>>>>>>>>>>>
!
csys,rib7cs         ! rib 1 coordinate system
!
p1=kp((w2/2+r1),ribh+topthk,0)  ! keypoints from position on rib 
!
csys,gcs
!
theta1st=theta10
theta2nd=theta1
!
asel,s,loc,y,theta1st,theta1st,,   ! select profile @ specified angle
asel,r,loc,x,rbore,kx(p1),,
lsla,s
!
aatt,,,2,0          ! set area element attributes
mshkey,1            ! select mapped meshing
!
*get,am,area,,num,max
*get,lm,line,,num,max
!
lsel,r,line,,lm-17,lm-16,,  ! select lines to be concatenated
lccat,all           ! concatenate lines
!
lsla,s              ! reselect line/area set
!
lsel,r,line,,lm-16,lm-15,,  ! select lines to be concatenated
lccat,all           ! concatenate lines
!
lsla,s
!
asel,u,area,,am-2,am-1,,    ! select areas for mapped meshing
!
amesh,all           ! map mesh areas
!
asel,s,loc,y,theta1st,theta1st,,   ! reselect areas for free meshing
asel,r,loc,x,rbore,kx(p1),,
asel,u,area,,am-8,am-3,,
asel,u,area,,am,am,,
!
aatt,,,2,0          ! set area element attributes
mshkey,0            ! select free meshing
!
amesh,all           ! free mesh areas
!
asel,all
lsel,all
!
*get,lm,line,,num,max
ldele,lm-1,lm,
!
!
!----------------- Rotate Area Meshes into Volume Meshes ---------------------
!
nerot=2             ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
asel,r,loc,x,rbore,kx(p1),,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt1,axpt2,phi,1   ! create volumes
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
asel,r,loc,x,rbore,kx(p1),,
aclear,all
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 1 Meshes into Area 2 ---------------------
!
csys,gcs
!
theta1st=theta1
theta2nd=theta2
!
nerot=5             ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
!
*get,am,area,,num,max
!
asel,u,area,,am-4,am-1,,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt1,axpt2,phi,2   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 3 Meshes into Area 4 ---------------------
!
csys,gcs
!
theta1st=theta3
theta2nd=theta4
!
nerot=5             ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
!
*get,am,area,,num,max
!
asel,u,area,,am-4,am-1,,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt1,axpt2,phi,1   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 4 Meshes into Area 5 ---------------------
!
csys,gcs
!
theta1st=theta4
theta2nd=theta5
!
nerot=4             ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
!
*get,am,area,,num,max
!
asel,u,area,,am-5,am-1,,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt1,axpt2,phi,1   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 10 Meshes into Area 9 ---------------------
!
csys,rib7cs         ! rib 7 coordinate system
!
p1=kp((w2/2+r1),ribh+topthk,0)  ! keypoints from position on rib 
!
csys,gcs
!
theta1st=theta9
theta2nd=theta10
!
nerot=5             ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta2nd,theta2nd,,   ! select area profile @ specified angle
asel,r,loc,x,rbore,kx(p1),,
!
*get,am,area,,num,max
!
asel,u,area,,am-4,am-1,,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt2,axpt1,phi,1   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 9 Meshes into Area 8 ---------------------
!
csys,gcs
!
theta1st=theta8
theta2nd=theta9
!
nerot=4             ! # elements along rotation
theta=-180-theta8
theta1st=theta
phi=abs(180-theta2nd-theta1st)  ! calculate rotation angle
!
asel,s,loc,y,theta2nd,theta2nd,,   ! select area profile @ specified angle
!
*get,am,area,,num,max
!
asel,u,area,,am-5,am-1,,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt2,axpt1,phi,1   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 8 Meshes into Area 7 ---------------------
!
csys,gcs
!
ksel,s,loc,x,rcb1,,   ! select keypoint of outer radius bounding areas
ksel,r,loc,y,theta8,
ksel,r,loc,z,0,,
*get,p1,kp,,num,max
ksel,all
!
theta1st=theta7
theta2nd=theta8
!
nerot=4             ! # elements along rotation
phi=abs(theta2nd-theta1st)  ! calculate rotation angle
!
asel,s,loc,y,theta2nd,theta2nd,,   ! select area profile @ specified angle
asel,r,loc,x,rbore,kx(p1),,
!
*get,am,area,,num,max
!
asel,u,area,,am-2,am-1,,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt2,axpt1,phi,1   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 5 Meshes into Area 6 ---------------------
!
csys,gcs
!
ksel,s,loc,x,rcb1,,   ! select keypoint of outer radius bounding areas
ksel,r,loc,y,theta5,
ksel,r,loc,z,0,,
*get,p1,kp,,num,max
ksel,all
!
theta1st=theta5
theta2nd=theta6
!
nerot=4             ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
asel,r,loc,x,rbore,kx(p1),,
!
*get,am,area,,num,max
!
asel,u,area,,am-2,am-2,,
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt1,axpt2,phi,1   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Rotate Area 6 Meshes into Area 7 ---------------------
!
csys,gcs
!
ksel,s,loc,x,rcb1,,   ! select keypoint of outer radius bounding areas
ksel,r,loc,y,theta6,
ksel,r,loc,z,0,,
*get,p1,kp,,num,max
ksel,all
!
theta1st=theta6
theta2nd=theta7
!
nerot=20            ! # elements along rotation
phi=abs(theta1st-theta2nd)  ! calculate rotation angle
!
asel,s,loc,y,theta1st,theta1st,,   ! select area profile @ specified angle
asel,r,loc,x,rbore,kx(p1),,
!
*get,am,area,,num,max
!
ESIZE,0,nerot,
vrotat,all,,,,,,axpt1,axpt2,phi,1   ! create volumes
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
!
!------------------------------------------------------------------------------------------
!  All points referenced by 'pxy', where the x and y denote a digit, indicates a KP:
!  
!               x:  arbitrarily selected  reference number
!               y:  angular position i.e., y means 'thetay' variable
!
!           Therefore, all points with the same x value on a given area at
!           different y angles are the same corresponding points on the 
!           different areas.
!------------------------------------------------------------------------------------------
!   
csys,gcs        ! Create lines forming volume
!
ksel,s,loc,x,rcb1,,   !----------- BEGIN line creation algoritm by select kp's ---------------
ksel,r,loc,y,theta6,
ksel,r,loc,z,0,,
*get,p16,kp,,num,max
!
ksel,s,loc,x,rcb1,,
ksel,r,loc,y,theta7,
ksel,r,loc,z,0,,
*get,p17,kp,,num,max
!
ksel,all
l,p16,p17   !----------------- END line creation algorithm -------------------------------
!
!
ksel,s,loc,x,kx(354)-.020,kx(354)+.020,   !--- BEGIN line creation algoritm by select kp's ---
ksel,r,loc,y,theta6,
ksel,r,loc,z,0,,
*get,p26,kp,,num,max
!
ksel,s,loc,x,kx(354)-.020,kx(354)+.020,
ksel,r,loc,y,theta7,
ksel,r,loc,z,0,,
*get,p27,kp,,num,max
!
ksel,all
l,p26,p27   !----------------- END line creation algorithm -------------------------------
!
!
ksel,s,loc,x,kx(354)-.020,kx(354)+.020,   !--- BEGIN line creation algoritm by select kp's ---
ksel,r,loc,y,theta6,
ksel,r,loc,z,-(topthk+ribh),,
*get,p36,kp,,num,max
!
ksel,s,loc,x,kx(354)-.020,kx(354)+.020,
ksel,r,loc,y,theta7,
ksel,r,loc,z,-(topthk+ribh),,
*get,p37,kp,,num,max
!
ksel,all
l,p36,p37   !----------------- END line creation algorithm -------------------------------
!
!
*get,lm,line,,num,max
!
lsel,s,line,,lm-2,lm,,
!
lesize,all,,,nerot,1,
lsel,all
!
asel,s,loc,x,kx(p16),kx(p16),,  ! create area set for creating volume
asel,r,loc,y,theta7,theta6,,
!
cm,aset1,area
!
asel,all
!
ksel,s,loc,x,rcb1,,         ! create KP's for creation of areas
ksel,r,loc,y,theta6,
ksel,r,loc,z,-cbore1,,
*get,p46,kp,,num,max
!
ksel,s,loc,x,rcb1,,
ksel,r,loc,y,theta7,
ksel,r,loc,z,-cbore1,,
*get,p47,kp,,num,max
!
ksel,s,loc,x,rcb1,,         ! create KP's for creation of area
ksel,r,loc,y,theta6,
ksel,r,loc,z,-(topthk+ribh),,
*get,p56,kp,,num,max
!
ksel,s,loc,x,rcb1,,
ksel,r,loc,y,theta7,
ksel,r,loc,z,-(topthk+ribh),,
*get,p57,kp,,num,max
!
ksel,all
!
a,p16,p46,p47,p17       ! forming remaining areas of volume
a,p16,p26,p27,p17
a,p26,p36,p37,p27
a,p36,p56,p57,p37
!
asel,s,loc,y,theta6,theta7,,
asel,r,loc,x,rcb1-.050,kx(p26)+.020,
!
va,all
!
asel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-------------------------------------------------------------------------------------
!
!
!---------------------- Transition Mesh Volume Between Areas 1 & 2 -------------------
!
csys,gcs
!
ksel,s,loc,x,rcb1,,     ! point p1
ksel,r,loc,y,theta1
ksel,r,loc,z,0,,
*get,p1,kp,,num,max
ksel,all
!
ksel,s,loc,x,p1r,,      ! point p5
ksel,r,loc,y,theta1
ksel,r,loc,z,0,,
*get,p5,kp,,num,max
ksel,all
!
csys,0
!
rp5=abs(kx(p5))
!
csys,gcs
!
ksel,s,loc,x,rp5,,      ! point p2
ksel,r,loc,y,90
ksel,r,loc,z,0,,
*get,p2,kp,,num,max
ksel,all
!
k,,kx(p1),ky(p2),kz(p1)
!
*get,km,kp,,num,max
!
ksel,s,loc,x,kx(km),,       ! point p3
ksel,r,loc,y,90
ksel,r,loc,z,-cbore1,,
*get,p3,kp,,num,max
ksel,all
!
ksel,s,loc,x,kx(km),,       ! point p4
ksel,r,loc,y,theta2
ksel,r,loc,z,0,,
*get,p4,kp,,num,max
ksel,all
!
l,p1,km
lstr,km,p2
lstr,km,p3
l,km,p4
lstr,13,406
!
ksel,s,loc,x,p1r,,      ! point p6
ksel,r,loc,y,theta2
ksel,r,loc,z,0,,
*get,p6,kp,,num,max
ksel,all
!
ksel,s,loc,x,rcb1,,     ! point p7
ksel,r,loc,y,theta1
ksel,r,loc,z,-cbore1,,
*get,p7,kp,,num,max
ksel,all
!
ksel,s,loc,x,rcb1,,     ! point p8
ksel,r,loc,y,theta2
ksel,r,loc,z,-cbore1,,
*get,p8,kp,,num,max
ksel,all
!
a,p1,p5,p2,km
a,km,p2,p6,p4
a,p1,km,p3,p7
a,km,p3,p8,p4
a,13,406,239,240,10
!
a,13,406,256,257,11
!
asel,s,loc,x,kx(p1),kx(p5)
asel,r,loc,y,theta2,theta1
!
csys,rib1cs
!
asel,u,loc,z,0,-1.00,,
asel,a,loc,z,0,0,,
!
va,all
!
asel,all
!
*get,vm,volu,,num,max
*get,lm,line,,num,max
!
lesize,lm-4,,,4,1,
lesize,lm-3,,,1,1,
lesize,lm-2,,,2,1,
lesize,lm-1,,,4,1,
lesize,lm,,,2,1,
!
vatt,,,1,
mopt,pyra,on
mshape,1,3d
mshkey,0
shpp,warn
!
vmesh,vm
!
allsel,all
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!---------------------- Map Mesh Volume Between Areas 7 & 8 ------------------
!
csys,gcs
!
l,355,385
lstr,369,399
!
a,369,399,385,355
a,153,143,394,364
a,192,369,399,197
a,355,385,386,356
!
*get,lm,line,,num,max
!
lesize,lm,,,4,1,
lesize,lm-1,,,4,1,
lesize,337,,,4,1,
lesize,338,,,4,1,
lesize,235,,,4,1,
!
asel,s,loc,x,kx(p16),kx(197)
asel,r,loc,y,ky(197),theta7
!
va,all
!
*get,vm,volu,,num,max
!
asel,all
!
asel,s,loc,x,kx(p16),kx(p16)
asel,r,loc,y,theta8,theta7
!
accat,all
!
asel,all
!
lsel,s,loc,y,theta8,theta8,,
lsel,r,loc,x,kx(p16),kx(p16),,
!
lccat,all
!
lsel,s,loc,y,theta7,theta7,,
lsel,r,loc,x,kx(p16),kx(p16),,
!
lccat,all
!
lsel,all
!
asel,s,loc,x,kx(p16),kx(197)
asel,r,loc,z,0,0,,
asel,r,loc,y,ky(197),theta7
!
accat,all
!
asel,s,loc,x,kx(p16),kx(197)
asel,r,loc,z,kz(153),kz(153),,
asel,r,loc,y,ky(197),theta7
!
accat,all
!
asel,s,loc,x,kx(p16),kx(197)
asel,r,loc,y,theta7-5,theta7
!
accat,all
!
asel,s,loc,x,kx(p16),kx(197)
asel,r,loc,y,theta8-5,theta8
!
accat,all
!
asel,all
!
vatt,,,1,
mshape,0,3d
mshkey,1
!
vmesh,vm
!
*get,am,area,,num,max
*get,lm,line,,num,max
!
ldele,lm-1,lm
adele,am-4,am
!
allsel,all
!
LDELE,858
LDELE,859
LDELE,860
LDELE,861
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------  
!
!
!---------------------- Map Mesh Volume Between Areas 5 & 6 -------------------
!
csys,gcs
!
l,370,340
lstr,384,354
!
a,354,384,190,188
a,340,370,384,354
a,340,370,371,341
a,105,117,349,379
!
*get,lm,line,,num,max
!
lesize,lm,,,4,1,
lesize,lm-1,,,4,1,
lesize,322,,,4,1,
lesize,321,,,4,1,
lesize,174,,,4,1,
!
asel,s,loc,x,kx(p16),kx(188)
asel,r,loc,y,theta6-5,theta5
!
va,all
!
*get,vm,volu,,num,max
!
asel,all
!
asel,s,loc,x,kx(p16),kx(p16)
asel,r,loc,y,theta6,theta5
!
accat,all
!
asel,all
!
lsel,s,loc,y,theta5,theta5,,
lsel,r,loc,x,kx(p16),kx(p16),,
!
lccat,all
!
lsel,s,loc,y,theta6,theta6,,
lsel,r,loc,x,kx(p16),kx(p16),,
!
lccat,all
!
lsel,all
!
asel,s,loc,x,kx(p16),kx(188)
asel,r,loc,z,0,0,,
asel,r,loc,y,theta6,theta5
!
accat,all
!
asel,s,loc,x,kx(p16),kx(188)
asel,r,loc,z,kz(117),kz(117),,
asel,r,loc,y,theta6,theta5
!
accat,all
!
asel,s,loc,x,kx(p16),kx(188)
asel,r,loc,y,theta6,theta6+5
!
accat,all
!
asel,s,loc,x,kx(p16),kx(188)
asel,r,loc,y,theta5,theta5
!
accat,all
!
asel,all
!
vatt,,,1,
mshape,0,3d
mshkey,1
!
vmesh,vm
!
*get,am,area,,num,max
*get,lm,line,,num,max
!
ldele,lm-1,lm
adele,am-4,am
!
LDELE,860
LDELE,861
LDELE,862
LDELE,863
!
LDELE,860
LDELE,861
LDELE,862
LDELE,863
!
nummrg,all,.003,.003        ! merge all coincident entities
numcmp,all          ! renumber entities
!
LPLOT
!
!-----------------------------------------------------------------------------
!
!
!----------------------- Create Volumes / Areas 8 to 10 ----------------------
!
vsel,s,volu,,27,    ! select triangular volume to modify
aslv,s
lsla,s
asel,u,area,,162,
!
vclear,all
vdele,all
adele,all,,,1
!
allsel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
lesize,301,,,neL7,1,
lesize,303,,,1,1,
lesize,305,,,neL7+1,1,
!
lccat,301,303       ! concatenate lines
!
esize,0,neribw-1,
aatt,,,2,0      ! set element attributes
mshkey,1        ! select mapped meshing
amesh,162       ! mesh areas
!
*get,lm,line,,num,max
ldele,lm
!
*get,dext,line,299,leng
!
voffst,162,dext,
!
voffst,553,dext,
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
csys,gcs
!
asel,s,loc,y,theta10,360+theta8
asel,r,loc,x,kx(384),kx(166)
!
a,194,397,383,324,180
a,180,324,293,290,288,166,167
a,383,384,325,324
a,324,325,294,293
a,143,392,332,212,210,144
!
a,182,211,212,332,301,302,289,176,414
!
adele,288,292,1,1
!
va,all
!
aslv,s
lsla,s
*get,lm,line,,num,max
!
lesize,lm-627,,,4,1,
lesize,lm-538,,,2,1,
lesize,lm-534,,,6,1,
lesize,lm-533,,,3,1,
lesize,lm-529,,,4,1,
!
lesize,lm-520,,,2,1,
lesize,lm-519,,,3,1,
lesize,lm-517,,,2,1,
lesize,lm-516,,,3,1,
lesize,lm-514,,,2,1,
!
lesize,lm-513,,,2,1,
lesize,lm-504,,,1,1,
lesize,lm-503,,,2,1,
lesize,lm-502,,,3,1,
lesize,lm-501,,,2,1,
!
lesize,lm-500,,,3,1,
lesize,lm-499,,,4,1,
lesize,lm-498,,,2,1,
lesize,lm-496,,,2,1,
lesize,lm-371,,,2,1,
!
lesize,lm-370,,,1,1,
lesize,lm-369,,,2,1,
lesize,lm-368,,,1,1,
lesize,lm-367,,,2,1,
lesize,lm-366,,,1,1,
!
lesize,lm-365,,,1,1,
lesize,lm-364,,,2,1,
lesize,lm-363,,,1,1,
lesize,lm-362,,,1,1,
lesize,lm-361,,,3,1,
!
lesize,lm-360,,,1,1,
lesize,lm-359,,,2,1,
lesize,lm-6,,,4,1,
lesize,lm-2,,,4,1,
lesize,lm-1,,,4,1,
lesize,lm,,,5,1,
!
allsel,all
!
lplot
!
!--------------------------------------------------------------------------------
!
!
!------------------------- Create Volumes / Areas 3 to 5 ------------------------
!
vsel,s,volu,,15,    ! select triangular volume to modify
aslv,s
lsla,s
asel,u,area,,79,
!
vclear,all
vdele,all
adele,all,,,1
!
allsel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
lesize,127,,,neL2,1,
lesize,136,,,1,1,
lesize,137,,,neL2+1,1,
!
lccat,127,136       ! concatenate lines
!
esize,0,neribw-1,
aatt,,,2,0      ! set element attributes
mshkey,1        ! select mapped meshing
amesh,79        ! mesh areas
!
*get,lm,line,,num,max
ldele,lm
!
*get,dext,line,122,leng
!
voffst,79,dext,
!
voffst,557,dext,
!
allsel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
*get,km,kp,,num,max
!
csys,gcs
!
asel,s,loc,y,theta5-5,theta3
asel,r,loc,x,rcb1,kx(km-356)
!
lsla,s
ksll,s
!
*get,km,kp,,num,max
*get,lm,line,,num,max
!
a,km-340,km-153,km-151,km-148,km-107,km-355,km-354
a,km-107,km-355,km-228,km-35,km-49
a,km-148,km-107,km-106,km-147
a,km-107,km-106,km-48,km-49
a,km-336,km-152,km-139,km-140,km-99,km-208,km-210,km-345,km
!
a,km-310,km-209,km-208,km-99,km-40,km-312
!
*get,am,area,,num,max
!
adele,am-297,am-293,1,1
!
va,all
!
*get,vm,volu,,num,max
vsel,s,volu,,vm
aslv,s
lsla,s
ksll,s
*get,lm,line,,num,max
!
lesize,lm-681,,,5,1,
lesize,lm-562,,,3,1,
lesize,lm-559,,,3,1,
lesize,lm-552,,,5,1,
lesize,lm-551,,,5,1,
!
lesize,lm-535,,,2,1,
lesize,lm-534,,,3,1,
lesize,lm-532,,,2,1,
lesize,lm-531,,,3,1,
lesize,lm-529,,,2,1,
!
lesize,lm-528,,,2,1,
lesize,lm-519,,,2,1,
lesize,lm-518,,,3,1,
lesize,lm-517,,,1,1,
lesize,lm-516,,,2,1,
!
lesize,lm-515,,,3,1,
lesize,lm-514,,,2,1,
lesize,lm-512,,,2,1,
lesize,lm-418,,,2,1,
lesize,lm-417,,,1,1,
!
lesize,lm-416,,,3,1,
lesize,lm-415,,,1,1,
lesize,lm-414,,,2,1,
lesize,lm-413,,,1,1,
lesize,lm-412,,,1,1,
!
lesize,lm-410,,,1,1,
lesize,lm-409,,,1,1,
lesize,lm-408,,,3,1,
lesize,lm-407,,,1,1,
lesize,lm-406,,,2,1,
!
lesize,lm-6,,,4,1,
lesize,lm-2,,,4,1,
lesize,lm-1,,,5,1,
lesize,lm,,,4,1,
!
allsel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber entities
!
lplot
!
!************************* End Center Section Construction ***********************
!
!
!************************ Modify Ribs For Wall Construction **********************
!
!
!----------------------------------- Ribs 5 & 6 ----------------------------------
!
csys,0
!
plndistr=5.8976   ! distance of cutting plane through ribs 5 & 6
          ! from global c.s. center
!
k,,5,1,plndistr
k,,-5,1,plndistr
k,,-5,-5,plndistr
k,,5,-5,plndistr
!
*get,km,kp,,num,max
!
a,km,km-1,km-2,km-3
!
*get,am,area,,num,max
!
csys,gcs
!
vsel,s,loc,y,theta9,360+theta7
vsel,r,loc,x,kx(127),15.00
!
*get,vhigh,volu,,num,max
*get,vlow,volu,,num,min
!
*do,i,vlow,vhigh,1
!
  *if,i,ne,vhigh,then
    vsba,i,am,,delete,keep
  *else
    vsba,i,am,,delete,delete
  *endif
!
*enddo
!
vsel,s,loc,y,theta9,360+theta7
vsel,r,loc,x,kx(127),15.00
!
csys,0
vsel,r,loc,z,plndistr,15.00
!
vdele,all,,,1
!
allsel,all
!
numcmp,all          ! renumber entities
!
LPLOT
!
!----------------------------------- Ribs 3 & 4 ----------------------------------
!
csys,0
!
plndistl=6.3667   ! distance of cutting plane through ribs 5 & 6
          ! from global c.s. center
!
k,,5,1,-plndistl
k,,-5,1,-plndistl
k,,-5,-5,-plndistl
k,,5,-5,-plndistl
!
*get,km,kp,,num,max
!
a,km,km-1,km-2,km-3
!
*get,am,area,,num,max
!
csys,gcs
!
vsel,s,loc,y,theta6,theta4
vsel,r,loc,x,kx(91),15.00
!
*get,vhigh,volu,,num,max
*get,vlow,volu,,num,min
!
*do,i,vlow,vhigh,1
!
  *if,i,ne,vhigh,then
    vsba,i,am,,delete,keep
  *else
    vsba,i,am,,delete,delete
  *endif
!
*enddo
!
vsel,s,loc,y,theta6,theta4
vsel,r,loc,x,kx(91),15.00
!
csys,0
vsel,r,loc,z,-plndistl,-15.00
!
vdele,all,,,1
!
allsel,all
!
numcmp,all          ! renumber entities
!
LPLOT
!
!---------------------- Front Housing Fillet / Right --------------------------
!
csys,0
!
k,,2.8471,-.500,.04358      ! coordinates creating front fillet surface
k,,2.5981,-.750,.02179
k,,2.5981,-1.136,.02179
k,,2.8471,0,.04358
!
*get,km,kp,,num,max
!
csn=csn+1
cskp,csn,0,km-2,km-1,km-3,,
!
k,,0,r1,0
!
csn=csn+1
RFfilcs=csn
cskp,RFfilcs,1,km+1,km-2,km-3,,
!
l,km-2,km-3
lstr,km,km-3
lstr,km-2,km-1
!
k,,r1,0,-10.00
!
lstr,km-2,km+2
!
*get,lm,line,,num,max
!
csys,rib5cs         ! select right/front rib area
!
asel,s,loc,z,0,-15.00
asel,r,loc,x,-(w2+r1),
asel,r,loc,y,-(ribh-r2),r2+topthk
!
adrag,lm-1,lm-2,lm-3,,,,lm,
!
*get,am,area,,num,max
!
aptn,all
!
adele,am+4,am+6,1,1
!
allsel,all
lsla,u
lsel,u,line,,1,6
!
ldele,all
!
allsel,all
numcmp,all          ! renumber entities
!
LPLOT
!
!----------------------- Front Housing Fillet / Left --------------------------
!
csys,0
!
k,,2.8471,-.500,-.04358
k,,2.5981,-.750,-.02179
k,,2.5981,-1.136,-.02179
k,,2.8471,0,-.04358
!
*get,km,kp,,num,max
!
csn=csn+1
cskp,csn,0,km-2,km-1,km-3,,
!
k,,0,r1,0
!
csn=csn+1
LFfilcs=csn
cskp,LFfilcs,1,km+1,km-2,km-3,,
!
l,km-2,km-3
lstr,km,km-3
lstr,km-2,km-1
!
k,,r1,0,10.00
!
lstr,km-2,km+2
!
*get,lm,line,,num,max
!
csys,rib4cs
!
asel,s,loc,z,0,-15.00
asel,r,loc,x,(w2+r1+r2),
asel,r,loc,y,-(ribh-r2),r2+topthk
!
adrag,lm-1,lm-2,lm-3,,,,lm,
!
*get,am,area,,num,max
!
aptn,all
!
adele,am+4,am+6,1,1
!
allsel,all
lsla,u
lsel,u,line,,1,6
!
ldele,all
!
allsel,all
numcmp,all          ! renumber entities
!
LPLOT
!
!----------------------------- Modify Front Volume & Mesh ---------------------------
!
csys,0
!
lstr,kp(2.8471,-.500,.04358),kp(2.8471,-.500,-.04358)
lstr,kp(2.5981,-.750,.02179),kp(2.5981,-.750,-.02179)
lstr,kp(2.5981,-1.136,.02179),kp(2.5981,-1.136,-.02179)
lstr,kp(2.8471,0,.04358),kp(2.8471,0,-.04358)
!
*get,km,kp,,num,max
!
a,kp(2.8471,0,.04358),kp(2.8471,0,-.04358),kp(2.8471,-.500,-.04358),kp(2.8471,-.500,.04358)
a,kp(2.8471,-.500,.04358),kp(2.8471,-.500,-.04358),kp(2.5981,-.750,-.02179),kp(2.5981,-.750,.02179)
a,kp(2.5981,-.750,.02179),kp(2.5981,-.750,-.02179),kp(2.5981,-1.136,-.02179),kp(2.5981,-1.136,.02179)
!
vsel,s,volu,,88,
aslv,s
cm,asetv88,area
vdele,88
asel,all
!
numcmp,all          ! renumber & compress all entities
!
k,,2.330,1,-3
k,,2.330,1,3
k,,2.330,-5,3
k,,2.330,-5,-3
!
*get,km,kp,,num,max
!
asel,s,area,,asetv88
!
a,km,km-1,km-2,km-3
!
*get,am,area,,num,max
!
aptn,all
adele,580,581,1,1
adele,584,588,588-584,1
!
va,all
!
numcmp,all          ! renumber & compress all entities
!
csys,0
!
asel,s,loc,x,2.330,
csys,RFfilcs-1
!
k,,kx(439),ky(421)+1,kz(439)
k,,kx(441),ky(421)+1,kz(441)
!
*get,km,kp,,num,max
!
a,439,km-1,km,441
!
*get,am,area,,num,max
!
asel,s,loc,y,0,r1
asel,r,loc,z,0,15.00
asel,a,area,,am
!
aptn,all
!
adele,589,,,1
!
numcmp,all          ! renumber & compress all entities
!
csys,LFfilcs-1
!
k,,kx(438),ky(431)+1,kz(438)
k,,kx(440),ky(431)+1,kz(440)
!
*get,km,kp,,num,max
!
a,438,km-1,km,440
!
*get,am,area,,num,max
!
asel,s,loc,y,0,r1
asel,r,loc,z,0,-15.00
asel,a,area,,am
!
aptn,all
!
adele,593,,,1
!
allsel,all
!
lplot
!
!---------------------------- Meshing of Volumes -------------------------------
!
vsel,s,volu,,101        ! size lines for meshing of volume
aslv,s
lsla,s
!
csys,gcs
!
lsel,r,loc,x,kx(299)
lesize,all,,,20,1,
lsla,s
!
lsel,r,loc,x,kx(328)
lsel,r,loc,z,0,0
lesize,all,,,5,1,
lsla,s
!
lsel,r,loc,x,kx(328)
lsel,r,loc,z,kz(440)
lesize,all,,,5,1,
lsla,s
!
lsel,r,loc,x,kx(439)
lsel,r,loc,y,ky(439)
lsel,r,loc,z,0,kz(441)
lesize,all,,,7,1,
lsla,s
!
lsel,r,loc,x,kx(438)
lsel,r,loc,y,ky(438)
lsel,r,loc,z,0,kz(440)
lesize,all,,,7,1,
lsla,s
!
csys,0
!
lsel,r,loc,x,2.330
lsel,r,loc,y,0,0,,
lesize,all,,,10,1,
lsla,s
!
lsel,r,loc,x,2.330
lsel,r,loc,y,ky(440)
lesize,all,,,10,1,
lsla,s
!
csys,0          !<<<<<<<<<<<<<<<<< concatenation of areas
asel,r,loc,x,2.330
asel,r,loc,z,kz(439),kz(438)
cm,a1,area
aslv,s
!
csys,gcs
asel,r,loc,x,kx(328)
asel,r,loc,y,360+ky(328),360+ky(313)
cm,a2,area
aslv,s
!
numcmp,area
!
asel,s,area,,a1
asel,a,area,,a2
accat,all       ! concatenation area 1
aslv,s
!
csys,gcs
asel,r,loc,x,kx(299)
accat,all       ! concatenation area 2
aslv,s
!
numcmp,line
!
csys,gcs        !<<<<<<<<<<<<<<<<< concatenation of lines
lsel,r,loc,x,kx(328)
lsel,r,loc,z,0
cm,l1,line
lsla,s
!
csys,0
lsel,r,loc,x,2.330
lsel,r,loc,y,0,0
cm,l2,line
lsla,s
!
lsel,s,line,,l1
lsel,a,line,,l2
lccat,all       ! concatenation line 1
lsla,s
!
csys,gcs
lsel,r,loc,x,kx(328)
lsel,r,loc,z,kz(441)
cm,l1,line
lsla,s
!
csys,0
lsel,r,loc,x,2.330
lsel,r,loc,y,ky(441)
cm,l2,line
lsla,s
!
lsel,s,line,,l1
lsel,a,line,,l2
lccat,all       ! concatenation line 2
lsla,s
!
csys,gcs
lsel,r,loc,x,kx(314)
lsel,r,loc,y,ky(314)
lccat,all       ! concatenation line 3
lsla,s
!
lsel,r,loc,x,kx(299)
lsel,r,loc,y,ky(299)
lccat,all       ! concatenation line 4
lsla,s
!
mshape,0,3d
vatt,,,1,0
mshkey,1
!
vmesh,all
!
aslv,s
!
*get,am,area,,num,max
*get,lm,line,,num,max
adele,am-1,am,1,1
ldele,lm-3,lm,1,1
!
allsel,all
!
nummrg,all,.001,.001        ! merge all coincident entities
numcmp,all          ! renumber & compress all entities
!
lplot
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Volume Modification of Right Front Rib
! 
csys,rib5cs         ! select right/front rib volume
!
vsel,s,loc,z,0,-15.00
vsel,r,loc,x,0,-(w2+r1)
vsel,r,loc,y,-(ribh-r2),r2+topthk
!
*get,vm,volu,,num,max
!
vsel,s,volu,,vm,
aslv,s
lsla,s
cm,aset1,area
vdele,vm
!
asel,r,loc,z,0,-15.00
asel,r,loc,x,-(w2+r1),
asel,r,loc,y,-(ribh-r2),r2+topthk
adele,all
asel,s,area,,aset1
!
asel,s,loc,z,0,-15.00
asel,r,loc,x,-(w2+r1),
asel,r,loc,y,-(ribh-r2),r2+topthk
cm,a1,area
asel,s,area,,aset1
!
asel,r,loc,y,r2+topthk
adele,all
asel,s,area,,aset1
!
asel,r,loc,y,-(ribh-r2)
vsla,s,0
aslv,s
cm,aset2,area
vdele,all
!
asel,s,area,,aset1
asel,r,loc,y,-(ribh-r2)
adele,all
!
asel,s,area,,aset1
asel,a,area,,a1
a,kp(-1.25,.625,0),kp(.125,.625,0),kp(.125,.625,-1.19),kp(-1.25,.625,-1.56),kp(-1.25,.625,-1.50)
a,kp(-1.25,-.511,0),kp(0,-.511,0),kp(0,-.511,-1.22),kp(-1.25,-.511,-1.56),kp(-1.25,-.511,-.765)
!
numcmp,volu
va,all
*get,vm,volu,,num,max
!
asel,s,area,,aset2
asel,r,loc,x,-w2,-(w2+r1)
asel,r,loc,y,-(ribh-r2+r1),-(ribh-r2)
adele,all,,,1
!
asel,s,area,,aset2
csys,FRfilcs
a,kp(.250,90,.441),kp(.250,90,1.23),kp(.250,90,2.00),kp(.2500,0,2.00),kp(.2500,0,.508)
cm,aset2,area
!
csys,rib5cs
asel,s,loc,z,0,-15.00
asel,r,loc,x,0,-(w2+r1)
asel,r,loc,y,-(ribh-r2)
asel,a,area,,aset2
!
va,all
!
allsel.all
numcmp,all
!
lplot
!