*Heading
 Multistep verification test for s4 elements using distributions
*Preprint, model=yes
*Node
*Nset, nset=side1, generate
*Nset, nset=side3, generate
*Elset, Elset=eall, generate
*Elset, elset=shellSecL1
**------------------------------------------------------------------------------------------
*Element, type=s4, elset=eall
*Shell General Section, elset=eall, density=.04
   0.30405E+05, 0.40541E+04, 0.40541E+05, 0.00000E+00, 0.00000E+00, 0.30000E+04, 0.15203E+03, 0.20270E+02
   0.00000E+00, 0.10135E+01, 0.20270E+02, 0.20270E+03, 0.00000E+00, 0.13514E+00, 0.13514E+01, 0.00000E+00
   0.00000E+00, 0.15000E+02, 0.00000E+00, 0.00000E+00, 0.10000E+00
*Distribution, name=distOrient, location=element, type=orientation, input=distrib_s4_orient_in.inp
*Distribution, name=distStiff, location=element, type=shell3d stiffness, input=distrib_s4_stiff_in.inp
*Element Properties
 orientation, distOrient
 shell3dstiffness, distStiff
*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
 sf,
 se,
*Node print
 u,
*Output, field
*Node output
 u,
*Element output
 sf,
*end step
**
** --------------------------------------------------
**
*step, name=step-2
 SSD mode-based
*steady state dynamics, interval=eigenfrequency
 60,
*El print, freq=999
 sf,
 se,
*Node print, freq=999
 u,
*Output, field
*Node output
 u,
*Element output
 sf,
*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
*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, 2, 2, 0.025
 side3, 3, 3, 0.050
*El print, freq=999
 sf,
 se,
*Node print, freq=999
 u,
*Output, field
*Node output
 u,
*Element output
 sf,
*Output, history
*Element output,elset=shellSecL1
 sf1,
 se1,
 sk1,
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-9
 static
*static
*El print, freq=999
 sf,
 se,
*Node print, freq=999
 u,
*Output, field
*Node output
 u,
*Element output
 sf,
*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