*Heading Multistep reference solution test for s4 elements with nodal thickness *Node 1, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00 2, 0.3333333E+00, 0.0000000E+00, 0.0000000E+00 3, 0.6666667E+00, 0.0000000E+00, 0.0000000E+00 4, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00 5, 0.0000000E+00, 0.3333333E+00, 0.0000000E+00 6, 0.3333333E+00, 0.3333333E+00, 0.0000000E+00 7, 0.6666667E+00, 0.3333333E+00, 0.0000000E+00 8, 0.1000000E+01, 0.3333333E+00, 0.0000000E+00 9, 0.0000000E+00, 0.6666667E+00, 0.0000000E+00 10, 0.3333333E+00, 0.6666667E+00, 0.0000000E+00 11, 0.6666667E+00, 0.6666667E+00, 0.0000000E+00 12, 0.1000000E+01, 0.6666667E+00, 0.0000000E+00 13, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00 14, 0.3333333E+00, 0.1000000E+01, 0.0000000E+00 15, 0.6666667E+00, 0.1000000E+01, 0.0000000E+00 16, 0.1000000E+01, 0.1000000E+01, 0.0000000E+00 *Nset, nset=side1, generate 1, 4, 1, *Nset, nset=side3, generate 13, 16, 1, *Nset, nset=nodeThick1, generate 1, 13, 4, 4, 16, 4, *Nset, nset=nodeThick2, generate 2, 14, 4, 3, 15, 4, *Elset, Elset=eall, generate 1, 9, 1, *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=shellSecL1 1, 1, 2, 6, 5, 4, 5, 6, 10, 9, 7, 9, 10, 14, 13, *Shell Section, orient=orientL1, material=mat1, elset=shellSecL1, offset=-0.50000E+00, nodal thickness 1., *Orientation, name=orientL1 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00 3,0. **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL2 2, 2, 3, 7, 6, 5, 6, 7, 11, 10, 8, 10, 11, 15, 14, *Shell Section, orient=orientL2, material=mat1, elset=shellSecL2, offset=-0.50000E+00 0.01, *Orientation, name=orientL2 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00, 0.0000000E+00 3,0. **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL3 3, 3, 4, 8, 7, 6, 7, 8, 12, 11, 9, 11, 12, 16, 15, *Shell Section, orient=orientL3, material=mat1, elset=shellSecL3, offset=-0.50000E+00, nodal thickness 1., *Orientation, name=orientL3 -0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00 3,0. *Boundary side1, 1, 6 *psd-definition, name=white, type=force 1.0,0.1,0.0 1.0,0.1,60.0 *spectrum, name=spec,type=acceleration 0., 0., 0. 1.e5, 10., 0. 1.e5, 100., 0. ** ** -------------------------------------------------- ** *step, name=step-1 frequency analysis *frequency, SIM=NO 5, *El print s, sf, ssavg, se, *El print 3, tshr, *Node print u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 3, 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 3, tshr, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 3, tshr, *cload side3, 4, 0.1 side3, 3, 2.0 side3, 2, 0.5 *select eigenmodes, generate 1, 5 *modal damping, rayleigh 1, 5, 0.5 *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 1, 5 *modal damping, rayleigh 1, 5, 0.5 *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 60., *select eigenmodes, generate 1, 5 *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 1,1 2,1 *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 1, 5 *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 3, tshr, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 3, tshr, *Output, history *Element output,elset=shellSecL1 s11, sf1, ssavg1, se1, sk1, sth, *Element output,elset=shellSecL1 3, tshr13, *end step ** ** -------------------------------------------------- ** *step, perturbation, name=step-9 static *static *El print, freq=999 s, sf, ssavg, se, sth, *El print, freq=999 3, tshr, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 3, 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