*HEADING
 MULTISTEP VERIFICATION TEST FOR S3R ELEMENTS WITH NODAL THICKNESS
*NODE
*NSET, NSET=SIDE1, GENERATE
*NSET, NSET=SIDE3, GENERATE
*NSET, NSET=NODETHICK1, GENERATE
*NSET, NSET=NODETHICK2, GENERATE
*ELSET, ELSET=SHELLSECL1, GENERATE
*ELSET, ELSET=SHELLSECL2, GENERATE
*ELSET, ELSET=SHELLSECL3, GENERATE
*ELSET, ELSET=SIDEEL
 SHELLSECL1, SHELLSECL3
*ELSET, ELSET=MIDEL
 SHELLSECL2,
*NODAL THICKNESS
 NODETHICK1, 0.01
 NODETHICK2, 0.02
**------------------------------------------------------------------------------------------
*MATERIAL, NAME=MAT1
*ELASTIC, TYPE=ENGINEERING CONSTANTS
  0.300E+07,  0.400E+07,  0.500E+07,  0.100E+00,  0.200E+00,  0.300E+00,  0.300E+06,  0.400E+06
  0.500E+06,
*DENSITY
  0.750E-03,
*MATERIAL, NAME=MAT2
*ELASTIC, TYPE=ENGINEERING CONSTANTS
  0.600E+07,  0.200E+07,  0.200E+07,  0.100E+00,  0.200E+00,  0.300E+00,  0.300E+06,  0.400E+06
  0.200E+06,
*DENSITY
  0.700E-03,
**------------------------------------------------------------------------------------------
*ELEMENT, TYPE=S3R, ELSET=EALL
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
 COORD3D, COORD3D
*DISTRIBUTION, NAME=DISTORIENT, LOCATION=ELEMENT, TABLE=ORIENTATIONTABLE
 SHELLSECL1,  0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00
 SHELLSECL2,  0.0000000E+00,   0.1000000E+01,   0.0000000E+00,  -0.1000000E+01,   0.0000000E+00,   0.0000000E+00
 SHELLSECL3, -0.1000000E+01,   0.0000000E+00,   0.0000000E+00,  0.0000000E+00,   -0.1000000E+01,   0.0000000E+00
*DISTRIBUTION TABLE, NAME=ORIANGLETABLE
 ANGLE
*DISTRIBUTION, NAME=DISTORIANGLE, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
 SHELLSECL1, 10
 SHELLSECL2, 15
 SHELLSECL3, 20
*DISTRIBUTION TABLE, NAME=THICKNESSTABLE
 LENGTH
*DISTRIBUTION, NAME=DISTTHICK, LOCATION=ELEMENT, TABLE=THICKNESSTABLE
 SHELLSECL2, 0.01
*ORIENTATION, NAME=ORI
 DISTORIENT
 3, DISTORIANGLE
*SHELL GENERAL SECTION, COMPOSITE, ELSET=MIDEL, OFFSET=-0.50000E+00, SHELL THICKNESS=DISTTHICK, ORIENTATION=ORI
 1.0, , MAT1, -30, PLY-1
 1.0, , MAT1, +45, PLY-2
 1.0, , MAT2, -45, PLY-3
 1.0, , MAT2, +45, PLY-4
*SHELL GENERAL SECTION, COMPOSITE, ELSET=SIDEEL, NODAL THICKNESS, OFFSET=-0.5, ORIENTATION=ORI
 1.0, , MAT1, -30, PLY-1
 1.0, , MAT1, +45, PLY-2
 1.0, , MAT2, -45, PLY-3
 1.0, , MAT2, +45, PLY-4
*BOUNDARY
 SIDE1, 1, 6
*PSD-DEFINITION, NAME=WHITE, TYPE=FORCE
*SPECTRUM, NAME=SPEC,TYPE=ACCELERATION
**
** --------------------------------------------------
**
*STEP, NAME=STEP-1
 FREQUENCY ANALYSIS
*FREQUENCY,SIM,normalization=mass
 5,
*EL PRINT
 S,
 SF,
 SSAVG,
 SE,
*EL PRINT
 TSHR,
*NODE PRINT
 U,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 S,
 SF,
 SSAVG,
 SE,
 STH,
*ELEMENT OUTPUT
 TSHR,
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-2
 SSD MODE-BASED
*STEADY STATE DYNAMICS, INTERVAL=EIGENFREQUENCY
 60,
*EL PRINT, FREQ=999
 S,
 SF,
 SSAVG,
 SE,
 STH,
*EL PRINT, FREQ=999
 TSHR,
*NODE PRINT, FREQ=999
 U,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 S,
 SF,
 SSAVG,
 SE,
 STH,
*ELEMENT OUTPUT
 TSHR,
*CLOAD
 SIDE3, 4, 0.1
 SIDE3, 3, 2.0
 SIDE3, 2, 0.5
*SELECT EIGENMODES, GENERATE
*MODAL DAMPING, RAYLEIGH
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-3
 MODAL DYNAMIC
*MODAL DYNAMIC
 0.1,1.
*CLOAD
 SIDE3, 4, 0.1
 SIDE3, 3, 0.2
 SIDE3, 2, 0.5
*SELECT EIGENMODES, GENERATE
*MODAL DAMPING, RAYLEIGH
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-4
 SSD DIRECT
*STEADY STATE DYNAMICS, DIRECT, INTERVAL=EIGENFREQUENCY
 60.,
*CLOAD
 SIDE3, 4, 0.1
 SIDE3, 3, 0.2
 SIDE3, 2, 0.5
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-5
 SSD SUBSPACE
*STEADY STATE DYNAMICS, SUBSPACE PROJECTION, INTERVAL=EIGENFREQUENCY
 60.,
*CLOAD
 SIDE3, 4, 0.1
 SIDE3, 3, 0.2
 SIDE3, 2, 0.5
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-6
 RANDOM RESPONSE
*RANDOM RESPONSE
*SELECT EIGENMODES, GENERATE
*CLOAD, LOAD CASE=1
 SIDE3, 4, 0.1
 SIDE3, 3, 0.2
 SIDE3, 2, 0.5
*CLOAD, LOAD CASE=2
 SIDE3, 4, 0.1
 SIDE3, 3, 0.2
 SIDE3, 2, 0.5
*CORRELATION, PSD=WHITE, TYPE=UNCORRELATED, COMPLEX=NO
*EL PRINT, FREQ=999
 SF,
 SE,
*NODE PRINT, FREQ=999
 U,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 SF,
 SE,
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-7
 RESPONSE SPECTRUM
*RESPONSE SPECTRUM, SUM=ABS, COMP=ALG
 SPEC,0.,0.,1.,1.
*SELECT EIGENMODES, GENERATE
*END STEP
**
** --------------------------------------------------
**
*STEP, NLGEOM=YES, NAME=STEP-8
 STATIC
*STATIC
 1.0, 1.0
*BOUNDARY
 SIDE3, 3, 3, 0.050
*EL PRINT, FREQ=999
 S,
 SF,
 SSAVG,
 SE,
 STH,
*EL PRINT, FREQ=999
 TSHR,
*NODE PRINT, FREQ=999
 U,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 S,
 SF,
 SSAVG,
 SE,
 STH,
*ELEMENT OUTPUT
 TSHR,
*END STEP
**
** --------------------------------------------------
**
*STEP, PERTURBATION, NAME=STEP-9
 STATIC
*STATIC
*EL PRINT, FREQ=999
 S,
 SF,
 SSAVG,
 SE,
 STH,
*EL PRINT, FREQ=999
 TSHR,
*NODE PRINT, FREQ=999
 U,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 S,
 SF,
 SSAVG,
 SE,
 STH,
*ELEMENT OUTPUT
 TSHR,
*LOAD CASE, NAME=LOAD1
*DLOAD
 EALL, BZ, 1.E3
*END LOAD CASE
*LOAD CASE, NAME=LOAD2
*DLOAD
 EALL, GRAV, 1.E6, 0.,0.,1.
*END LOAD CASE
*LOAD CASE, NAME=LOAD3
*DLOAD
 EALL, CENT, 100.,0.,0.,0.,1.,0.
*END LOAD CASE
*LOAD CASE, NAME=LOAD4
*DLOAD
 EALL, CENTRIF, 1.E6,0.,0.,0.,1.,0.
*END LOAD CASE
*END STEP