*Heading Multistep reference solution test for s4 elements *preprint,model=yes *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, *Elset, Elset=eall, generate 1, 9, 1, **------------------------------------------------------------------------------------------ *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, *Shell General Section, orient=orientL1, material=mat1, elset=shellSecL1, offset=-0.50000E+00 0.0100, *Orientation, name=orientL1 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00 3,0.5 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL2 2, 2, 3, 7, 6, *Shell General Section, orient=orientL2, material=mat1, elset=shellSecL2, offset=-0.47500E+00 0.0200, *Orientation, name=orientL2 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00 3,1.0 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL3 3, 3, 4, 8, 7, *Shell General Section, orient=orientL3, material=mat1, elset=shellSecL3, offset=-0.45000E+00 0.0300, *Orientation, name=orientL3 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00 3,2.0 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL4 4, 5, 6, 10, 9, *Shell General Section, orient=orientL4, material=mat1, elset=shellSecL4, offset=-0.42500E+00 0.0400, *Orientation, name=orientL4 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00 3,4.0 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL5 5, 6, 7, 11, 10, *Shell General Section, orient=orientL5, material=mat1, elset=shellSecL5, offset=-0.40000E+00 0.0500, *Orientation, name=orientL5 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00 3,8.0 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL6 6, 7, 8, 12, 11, *Shell General Section, orient=orientL6, material=mat1, elset=shellSecL6, offset=-0.37500E+00 0.0600, *Orientation, name=orientL6 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00 3,9.0 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL7 7, 9, 10, 14, 13, *Shell General Section, orient=orientL7, material=mat1, elset=shellSecL7, offset=-0.35000E+00 0.0700, *Orientation, name=orientL7 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00 3,10.0 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL8 8, 10, 11, 15, 14, *Shell General Section, orient=orientL8, material=mat1, elset=shellSecL8, offset=-0.32500E+00 0.0800, *Orientation, name=orientL8 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00 3,-8.0 **------------------------------------------------------------------------------------------ *Element, type=s4, elset=shellSecL9 9, 11, 12, 16, 15, *Shell General Section, orient=orientL9, material=mat1, elset=shellSecL9, offset=-0.30000E+00 0.0900, *Orientation, name=orientL9 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00 3,-10.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 2, tshr, *Node print u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 2, 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 2, tshr, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 2, 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, 2, 2, 0.025 side3, 3, 3, 0.050 *El print, freq=999 s, sf, ssavg, se, sth, *El print, freq=999 2, tshr, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 2, tshr, *Output, history *Element output,elset=shellSecL1 s11, sf1, ssavg1, se1, sk1, sth, *Element output,elset=shellSecL1 2, tshr13, *end step ** ** -------------------------------------------------- ** *step, perturbation, name=step-9 static *static *El print, freq=999 s, sf, ssavg, se, sth, *El print, freq=999 2, tshr, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 2, 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