*HEADING
 MULTISTEP VERIFICATION TEST FOR S4R ELEMENTS USING DISTRIBUTIONS
*PREPRINT, MODEL=NO
*NODE
*NSET, NSET=SIDE1, GENERATE
*NSET, NSET=SIDE3, GENERATE
*ELSET, ELSET=EALL, GENERATE
*ELSET, ELSET=SHELLSECL1
**------------------------------------------------------------------------------------------
*MATERIAL, NAME=MAT1
*ELASTIC, TYPE=ENGINEERING CONSTANTS
  0.310E+07,  0.410E+07,  0.510E+07,  0.110E+00,  0.210E+00,  0.310E+00,  0.310E+06,  0.410E+06
  0.510E+06,
*DENSITY
  0.751E-03,
*MATERIAL, NAME=MAT2
*ELASTIC, TYPE=ENGINEERING CONSTANTS
  0.320E+07,  0.420E+07,  0.520E+07,  0.120E+00,  0.220E+00,  0.320E+00,  0.320E+06,  0.420E+06
  0.520E+06,
*DENSITY
  0.752E-03,
*MATERIAL, NAME=MAT3
*ELASTIC, TYPE=ENGINEERING CONSTANTS
  0.330E+07,  0.430E+07,  0.530E+07,  0.130E+00,  0.230E+00,  0.330E+00,  0.330E+06,  0.430E+06
  0.530E+06,
*DENSITY
  0.753E-03,
*MATERIAL, NAME=MAT4
*ELASTIC, TYPE=ENGINEERING CONSTANTS
  0.340E+07,  0.440E+07,  0.540E+07,  0.140E+00,  0.240E+00,  0.340E+00,  0.340E+06,  0.440E+06
  0.540E+06,
*DENSITY
  0.754E-03,
**------------------------------------------------------------------------------------------
*ELEMENT, TYPE=S4R, ELSET=EALL
*SHELL GENERAL SECTION, COMPOSITE, ELSET=EALL, OFFSET=DISTOFFSET, SHELL THICKNESS=DISTTHICKSECT, ORIENTATION=ORI,LAYUP=LAYUP-1
 DISTTHICKPLY1, 3, MAT1,  DISTANGLEPLY1, PLY-1
 DISTTHICKPLY2, 3, MAT2,            -60, PLY-2
        .01253, 3, MAT3,  DISTANGLEPLY3, PLY-3
 DISTTHICKPLY4, 3, MAT4,  DISTANGLEPLY4, PLY-4
*ORIENTATION, NAME=ORI
 DISTORIENT
 3, DISTORIENTANGLE
*DISTRIBUTION, NAME=DISTORIENT, LOCATION=ELEMENT, TABLE=ORIENTATIONTABLE
*DISTRIBUTION, NAME=DISTORIENTANGLE, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
*DISTRIBUTION, NAME=DISTTHICKSECT, LOCATION=ELEMENT, TABLE=THICKNESSTABLE
*DISTRIBUTION, NAME=DISTOFFSET, LOCATION=ELEMENT, TABLE=OFFSETTABLE
*DISTRIBUTION, NAME=DISTTHICKPLY1, LOCATION=ELEMENT, TABLE=THICKNESSTABLE
*DISTRIBUTION, NAME=DISTTHICKPLY2, LOCATION=ELEMENT, TABLE=THICKNESSTABLE
*DISTRIBUTION, NAME=DISTTHICKPLY4, LOCATION=ELEMENT, TABLE=THICKNESSTABLE
*DISTRIBUTION, NAME=DISTANGLEPLY1, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
*DISTRIBUTION, NAME=DISTANGLEPLY3, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
*DISTRIBUTION, NAME=DISTANGLEPLY4, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
 COORD3D, COORD3D
*DISTRIBUTION TABLE, NAME=ORIANGLETABLE
 ANGLE
*DISTRIBUTION TABLE, NAME=THICKNESSTABLE
 LENGTH
*DISTRIBUTION TABLE, NAME=OFFSETTABLE
 RATIO
*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,
 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,
 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,
*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
*RESTART, WRITE
*BOUNDARY
 SIDE3, 2, 2, 0.025
 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,
 STH,
*ELEMENT OUTPUT
 TSHR,
*OUTPUT, HISTORY
*ELEMENT OUTPUT,ELSET=SHELLSECL1
 S11,
 SF1,
 SSAVG1,
 SE1,
 SK1,
 STH,
*ELEMENT OUTPUT,ELSET=SHELLSECL1
 TSHR13,
*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,
 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