*Heading
 Multistep verification test for s4 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
*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,
**------------------------------------------------------------------------------------------
*Element, type=s4, elset=eall
*Shell General Section, material=mat1, elset=eall, offset=-0.50000E+00, nodal thickness
  1.,
*Distribution, name=distOrient, location=element, type=orientation
 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, name=distThick, location=element, type=scalar
 shellSecL2, 0.01
*Element Properties
 orientation, distOrient
 thickness, distThick
*Boundary
 side1, 1, 6
*psd-definition, name=white, type=force
*spectrum, name=spec,type=acceleration
**
** --------------------------------------------------
**
*step, name=step-1
 frequency analysis
*frequency
 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
*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
*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
*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.
*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,
*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,
** 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