*Heading
 Multistep verification test for m3d4r elements using distributions
*Preprint, model=yes
*Node
*Nset, nset=side1, generate
*Nset, nset=side2
*Nset, nset=side3, generate
*Nset, nset=side4
*Nset, nset=load
*Elset, Elset=eall, generate
*Elset, elset=MembraneSecL1
**------------------------------------------------------------------------------------------
*Material, name=mat
*Elastic
 1.e6, 0.3
*Density
  0.1E-03,
**------------------------------------------------------------------------------------------
*Element, type=m3d4r, elset=eall
*Membrane Section, elset=eall, material=mat,  membrane thickness=distThickSect, orientation=ori
 0.1
***Orientation, name=ori
** 0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00
** 3, 0.
*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 Table, name=orientationTable
 coord3d, coord3d
*Distribution Table, name=oriAngleTable
 angle
*Distribution Table, name=thicknessTable
 length
**------------------------------------------------------------------------------------------
*Initial Conditions, Type=Stress
 1, 101., 101.
 2, 102., 102.
 3, 103., 103.
 4, 104., 104.
 5, 105., 105.
 6, 106., 106.
 7, 107., 107.
 8, 108., 108.
*Boundary
 side1, 1, 3
 side2, 1, 3
 side3, 1, 3
 side4, 1, 3
*psd-definition, name=white, type=force
*spectrum, name=spec,type=acceleration
**
** --------------------------------------------------
**
*step, name=step-0, nlgeom
 Resolve initial stresses
*static
*El print
 s,
*Node print
 u,
*Output, field
*Node output
 u,
*Element output
 s,
 sth,
*end step
**
** --------------------------------------------------
**
*step, name=step-1
 frequency analysis
*frequency
 12,
*Node print
 u,
*Output, field
*Node output
 u,
*Element output
 sth
*end step
**
** --------------------------------------------------
**
*step, name=step-2
 SSD mode-based
*steady state dynamics, interval=eigenfrequency
 60.
*El print, freq=999
 s,
 sth,
*Node print, freq=999
 u,
*Output, field
*Node output
 u,
*Element output
 s,
 sth,
*Cload
 load, 1,1.
 load, 3,1.
*Cload, imaginary
 load, 2,1.
*select eigenmodes, generate
*modal damping, rayleigh
*end step
**
** --------------------------------------------------
**
*step, name=step-3
 modal dynamic
*modal dynamic
 0.1,1.
*Cload
 load, 1,1.
 load, 2,1.
 load, 3,1.
*select eigenmodes, generate
*modal damping, rayleigh
*end step
**
** --------------------------------------------------
**
*step, name=step-4
 SSD direct
*steady state dynamics, direct, interval=eigenfrequency
 60.,
*Cload
 load, 1,1.
 load, 3,1.
*Cload, imaginary
 load, 2,1.
*end step
**
** --------------------------------------------------
**
*step, name=step-5
 SSD subspace
*steady state dynamics, subspace projection, interval=eigenfrequency
 60.,
*Cload
 load, 1,1.
 load, 3,1.
*Cload, imaginary
 load, 2,1.
*end step
**
** --------------------------------------------------
**
*step, name=step-6
 random response
*random response
*select eigenmodes, generate
*cload, load case=1
 load, 1, 1.
 load, 3, 1.
*cload, load case=2
 load, 2, 1.
*correlation, psd=white, type=uncorrelated, complex=no
*Node print, freq=999
 u,
*Output, field
*Node output
 u,
*end step
**
** --------------------------------------------------
**
*step, name=step-7
 response spectrum
*response spectrum, sum=abs, comp=alg
 spec,1.,1.,1.,1.
*select eigenmodes, generate
*end step
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-8
 static
*static
 1.0, 1.0
*Restart, write
*Boundary, op=new
 side1, 3, 3
 side2, 3, 3
 side3,3,3
 side4,3,3
 side1, 2, 2
 side4, 1, 1
 side2, 1,1,.025
 side3, 2,2,.025
*El print, freq=999
 s,
 sth,
*Node print, freq=999
 u,
*Output, field
*Node output
 u,
*Element output
 s,
 sth,
*Output, history
*Element output,elset=MembraneSecL1
 s11,
 sth,
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-9
 static
*static
*El print, freq=999
 s,
 sth,
*Node print, freq=999
 u,
*Output, field
*Node output
 u,
*Element output
 s,
 sth,
*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