FEA SOL3: STATIC ANALYSIS OF A SANDWICH PANEL

! covert everything to a consistent unit system
! [N],[NM],[N/m^2=Pa],[N/m]
P = 7e6
! facesheet
EXFS = 66.05E9
EYFS = 66.05E9
EZFS = 11.45E9
NUXYFS = 0.033
NUYZFS = 0.352
NUXZFS = 0.352
GXYFS = 4.88E9
GYZFS = 5.04E9
GXZFS = 5.04E9
TF = 0.00025
! core
EC = 127E6
NUC = 0.3
TC = 0.0125
! parameters
W = 0.15
H = 0.15
D = 0.025
DIV1 = 10
DIV2 = 7
DIV3 = 5
DIV4 = DIV1/2
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,100,100,100,0
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,0,0,0,15
/PREP7/PREP7 Enters the model creation preprocessor.
! define elements
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,1,SHELL281
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,2,SOLID186
! create material properties
! orthotropic (facesheet)
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EX,1,,EXFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EY,1,,EYFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EZ,1,,EZFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,1,,NUXYFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRYZ,1,,NUYZFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXZ,1,,NUXZFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,GXY,1,,GXYFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,GYZ,1,,GYZFS
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,GXZ,1,,GXZFS
! isotropic (core)
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EX,2,,EC
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,2,,NUC
! define shell layup
SECTSECTYPE, SECID, Type, Subtype, Name, REFINEKEY Associates section type information with a section ID number.,1,shell,,
SECDATASECDATA, VAL1, VAL2, VAL3, VAL4, VAL5, VAL6, VAL7, VAL8, VAL9, VAL10, VAL11, VAL12 Describes the geometry of a section.,TFS,1,0.0,3
SECOFFSETSECOFFSET, Location, OFFSET1, OFFSET2, CG-Y, CG-Z, SH-Y, SH-Z Defines the section offset for cross sections.,MID
!
! create local coordinatesystem
LOCALLOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2 Defines a local coordinate system by a location and orientation.,11,0,0,0,0,,,,1,1
! set element properties
TYPETYPE, ITYPE Sets the element type attribute pointer.,1
MATMAT, MAT Sets the element material attribute pointer.,1
REALREAL, NSET Sets the element real constant set attribute pointer.,
ESYSESYS, KCN Sets the element coordinate system attribute pointer.,11
SECNUMSECNUM, SECID Sets the element section attribute pointer.,1
TSHAPTSHAP, Shape Defines simple 2-D and 3-D geometric surfaces for target segment elements.,LINE
! create rectangles
RECTNGRECTNG, X1, X2, Y1, Y2 Creates a rectangular area anywhere on the working plane.,0,W/2,0,H/2,
RECTNGRECTNG, X1, X2, Y1, Y2 Creates a rectangular area anywhere on the working plane.,0,D,0,H/2,
RECTNGRECTNG, X1, X2, Y1, Y2 Creates a rectangular area anywhere on the working plane.,0,W/2,0,D,
RECTNGRECTNG, X1, X2, Y1, Y2 Creates a rectangular area anywhere on the working plane.,0,2*D,0,H/2,
RECTNGRECTNG, X1, X2, Y1, Y2 Creates a rectangular area anywhere on the working plane.,0,W/2,0,2*D
! overlap areas
AOVLAPAOVLAP, NA1, NA2, NA3, NA4, NA5, NA6, NA7, NA8, NA9 Overlaps areas.,ALL
! set element size
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,21,,,DIV1
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,22,,,DIV1
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,37,,,DIV2
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,43,,,DIV2
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,26,,,DIV3
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,34,,,DIV3
! mesh areas as mapped
MSHKEYMSHKEY, KEY Specifies whether free meshing or mapped meshing should be used to mesh a model.,1
AMESHAMESH, NA1, NA2, NINC Generates nodes and area elements within areas.,ALL
! copy areas
AGENAGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE Generates additional areas from a pattern of areas.,2,ALL,,,,,-TF/2
! change elementy type and material properties
TYPETYPE, ITYPE Sets the element type attribute pointer.,2
MATMAT, MAT Sets the element material attribute pointer.,2
! set extrude options
TYPETYPE, ITYPE Sets the element type attribute pointer.,2
EXTOPTEXTOPT, Lab, Val1, Val2, Val3, Val4 Controls options relating to the generation of volume elements from area elements.,ESIZE,6
EXTOPTEXTOPT, Lab, Val1, Val2, Val3, Val4 Controls options relating to the generation of volume elements from area elements.,ACLEAR
!
EXTOPTEXTOPT, Lab, Val1, Val2, Val3, Val4 Controls options relating to the generation of volume elements from area elements.,ATTR,
MATMAT, MAT Sets the element material attribute pointer.,2
REALREAL, NSET Sets the element real constant set attribute pointer.
ESYSESYS, KCN Sets the element coordinate system attribute pointer.,11
! select areas to be extruded
ASELASEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of areas.,S,LOC,Z,-TF/2
! extrude elements
VEXTVEXT, NA1, NA2, NINC, DX, DY, DZ, RX, RY, RZ Generates additional volumes by extruding areas.,ALL,,,0,0,-TC
!
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
!
! create bottom facesheet
TYPETYPE, ITYPE Sets the element type attribute pointer.,1
MATMAT, MAT Sets the element material attribute pointer.,1
ASELASEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of areas.,S,LOC,Z,-TC-TF/2
AGENAGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE Generates additional areas from a pattern of areas.,2,ALL,,,,,-TF/2
ASELASEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of areas.,S,LOC,Z,-TC-TF
MSHKEYMSHKEY, KEY Specifies whether free meshing or mapped meshing should be used to mesh a model.,1
AMESHAMESH, NA1, NA2, NINC Generates nodes and area elements within areas.,ALL
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
! macro to create rigid connection beteen the top facesheet and the core
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Z,-TF/2
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NUM1,NODE,0,COUNT
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,MIN1,NODE,0,NUM,MIN
CURNOD = MIN1
*DO*DO, Par, IVAL, FVAL, INC Defines the beginning of a do-loop.,I,1,NUM1
  ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
SLAVE = NODE(NX(CURNOD),NY(CURNOD),0)
  CERIGCERIG, MASTE, SLAVE, Ldof, Ldof2, Ldof3, Ldof4, Ldof5 Defines a rigid region.,CURNOD,SLAVE,UXYZ
  NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Z,-TF/2
CURNOD = NDNEXT(CURNOD)
*ENDDO*ENDDO Ends a do-loop and starts the looping action.
! macro to create rigid connection between the bottom facesheet and the core
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Z,-TC-TF/2
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NUM1,NODE,0,COUNT
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,MIN1,NODE,0,NUM,MIN
CURNOD = MIN1
*DO*DO, Par, IVAL, FVAL, INC Defines the beginning of a do-loop.,I,1,NUM1
  ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
SLAVE = NODE(NX(CURNOD),NY(CURNOD),-TC-TF)
  CERIGCERIG, MASTE, SLAVE, Ldof, Ldof2, Ldof3, Ldof4, Ldof5 Defines a rigid region.,CURNOD,SLAVE,UXYZ
  NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Z,-TC-TF/2
CURNOD = NDNEXT(CURNOD)
*ENDDO*ENDDO Ends a do-loop and starts the looping action.
! define kinematic constraints
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,X,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UX,0
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Y,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UY,0
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Y,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UY,0
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,X,W/2
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UX,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UZ,0
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Y,H/2
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UY,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UZ,0
! apply loads
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
ASELASEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of areas.,S,LOC,Z,0
SFASFA, Area, LKEY, Lab, VALUE, VALUE2 Specifies surface loads on the selected areas.,ALL,1,PRES,-P
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
FINISHFINISH Exits normally from a processor.
/NERR,5,100000,,0,5,
! solve the finite element equation
/SOL/SOLU Enters the solution processor.
SOLVESOLVE, Action Starts a solution.