*Heading Multistep reference solution for cpe8 elements *Node 1, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00 2, 0.1666667E+00, 0.0000000E+00, 0.0000000E+00 3, 0.3333333E+00, 0.0000000E+00, 0.0000000E+00 4, 0.5000000E+00, 0.0000000E+00, 0.0000000E+00 5, 0.6666667E+00, 0.0000000E+00, 0.0000000E+00 6, 0.8333333E+00, 0.0000000E+00, 0.0000000E+00 7, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00 8, 0.0000000E+00, 0.1666667E+00, 0.0000000E+00 9, 0.1666667E+00, 0.1666667E+00, 0.0000000E+00 10, 0.3333333E+00, 0.1666667E+00, 0.0000000E+00 11, 0.5000000E+00, 0.1666667E+00, 0.0000000E+00 12, 0.6666667E+00, 0.1666667E+00, 0.0000000E+00 13, 0.8333333E+00, 0.1666667E+00, 0.0000000E+00 14, 0.1000000E+01, 0.1666667E+00, 0.0000000E+00 15, 0.0000000E+00, 0.3333333E+00, 0.0000000E+00 16, 0.1666667E+00, 0.3333333E+00, 0.0000000E+00 17, 0.3333333E+00, 0.3333333E+00, 0.0000000E+00 18, 0.5000000E+00, 0.3333333E+00, 0.0000000E+00 19, 0.6666667E+00, 0.3333333E+00, 0.0000000E+00 20, 0.8333333E+00, 0.3333333E+00, 0.0000000E+00 21, 0.1000000E+01, 0.3333333E+00, 0.0000000E+00 22, 0.0000000E+00, 0.5000000E+00, 0.0000000E+00 23, 0.1666667E+00, 0.5000000E+00, 0.0000000E+00 24, 0.3333333E+00, 0.5000000E+00, 0.0000000E+00 25, 0.5000000E+00, 0.5000000E+00, 0.0000000E+00 26, 0.6666667E+00, 0.5000000E+00, 0.0000000E+00 27, 0.8333333E+00, 0.5000000E+00, 0.0000000E+00 28, 0.1000000E+01, 0.5000000E+00, 0.0000000E+00 29, 0.0000000E+00, 0.6666667E+00, 0.0000000E+00 30, 0.1666667E+00, 0.6666667E+00, 0.0000000E+00 31, 0.3333333E+00, 0.6666667E+00, 0.0000000E+00 32, 0.5000000E+00, 0.6666667E+00, 0.0000000E+00 33, 0.6666667E+00, 0.6666667E+00, 0.0000000E+00 34, 0.8333333E+00, 0.6666667E+00, 0.0000000E+00 35, 0.1000000E+01, 0.6666667E+00, 0.0000000E+00 36, 0.0000000E+00, 0.8333333E+00, 0.0000000E+00 37, 0.1666667E+00, 0.8333333E+00, 0.0000000E+00 38, 0.3333333E+00, 0.8333333E+00, 0.0000000E+00 39, 0.5000000E+00, 0.8333333E+00, 0.0000000E+00 40, 0.6666667E+00, 0.8333333E+00, 0.0000000E+00 41, 0.8333333E+00, 0.8333333E+00, 0.0000000E+00 42, 0.1000000E+01, 0.8333333E+00, 0.0000000E+00 43, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00 44, 0.1666667E+00, 0.1000000E+01, 0.0000000E+00 45, 0.3333333E+00, 0.1000000E+01, 0.0000000E+00 46, 0.5000000E+00, 0.1000000E+01, 0.0000000E+00 47, 0.6666667E+00, 0.1000000E+01, 0.0000000E+00 48, 0.8333333E+00, 0.1000000E+01, 0.0000000E+00 49, 0.1000000E+01, 0.1000000E+01, 0.0000000E+00 *Nset, nset=side1, generate 1, 7, 1, *Nset, nset=side3, generate 43, 49, 1, *Elset, Elset=eall, generate 1, 9, 1, **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL1 1, 1, 3, 17, 15, 2, 10, 16, 8, *Solid Section, orient=orientL1, material=mat1, elset=solidSecL1 1., *Orientation, name=orientL1 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00 3, 1.0 *Material, name=mat1 *Elastic, type=engineering constants 0.310E+07, 0.410E+07, 0.510E+07, 0.110E+00, 0.210E+00, 0.310E+00, 0.310E+06, 0.410E+06 0.510E+06, *Density 0.751E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL2 2, 3, 5, 19, 17, 4, 12, 18, 10, *Solid Section, orient=orientL2, material=mat2, elset=solidSecL2 1., *Orientation, name=orientL2 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00 3, -1.0 *Material, name=mat2 *Elastic, type=engineering constants 0.320E+07, 0.420E+07, 0.520E+07, 0.120E+00, 0.220E+00, 0.320E+00, 0.320E+06, 0.420E+06 0.520E+06, *Density 0.752E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL3 3, 5, 7, 21, 19, 6, 14, 20, 12, *Solid Section, orient=orientL3, material=mat3, elset=solidSecL3 1., *Orientation, name=orientL3 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00 3, 2.0 *Material, name=mat3 *Elastic, type=engineering constants 0.330E+07, 0.430E+07, 0.530E+07, 0.130E+00, 0.230E+00, 0.330E+00, 0.330E+06, 0.430E+06 0.530E+06, *Density 0.753E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL4 4, 15, 17, 31, 29, 16, 24, 30, 22, *Solid Section, orient=orientL4, material=mat4, elset=solidSecL4 1., *Orientation, name=orientL4 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00 3, -2.0 *Material, name=mat4 *Elastic, type=engineering constants 0.340E+07, 0.440E+07, 0.540E+07, 0.140E+00, 0.240E+00, 0.340E+00, 0.340E+06, 0.440E+06 0.540E+06, *Density 0.754E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL5 5, 17, 19, 33, 31, 18, 26, 32, 24, *Solid Section, orient=orientL5, material=mat5, elset=solidSecL5 1., *Orientation, name=orientL5 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00 3, 3.0 *Material, name=mat5 *Elastic, type=engineering constants 0.350E+07, 0.450E+07, 0.550E+07, 0.150E+00, 0.250E+00, 0.350E+00, 0.350E+06, 0.450E+06 0.550E+06, *Density 0.755E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL6 6, 19, 21, 35, 33, 20, 28, 34, 26, *Solid Section, orient=orientL6, material=mat6, elset=solidSecL6 1., *Orientation, name=orientL6 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00 3, -3.0 *Material, name=mat6 *Elastic, type=engineering constants 0.360E+07, 0.460E+07, 0.560E+07, 0.160E+00, 0.260E+00, 0.360E+00, 0.360E+06, 0.460E+06 0.560E+06, *Density 0.756E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL7 7, 29, 31, 45, 43, 30, 38, 44, 36, *Solid Section, orient=orientL7, material=mat7, elset=solidSecL7 1., *Orientation, name=orientL7 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00 3, 4.0 *Material, name=mat7 *Elastic, type=engineering constants 0.370E+07, 0.470E+07, 0.570E+07, 0.170E+00, 0.270E+00, 0.370E+00, 0.370E+06, 0.470E+06 0.570E+06, *Density 0.757E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL8 8, 31, 33, 47, 45, 32, 40, 46, 38, *Solid Section, orient=orientL8, material=mat8, elset=solidSecL8 1., *Orientation, name=orientL8 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00 3, -4.0 *Material, name=mat8 *Elastic, type=engineering constants 0.380E+07, 0.480E+07, 0.580E+07, 0.180E+00, 0.280E+00, 0.380E+00, 0.380E+06, 0.480E+06 0.580E+06, *Density 0.758E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 **------------------------------------------------------------------------------------------ *Element, type=cpe8, elset=solidSecL9 9, 33, 35, 49, 47, 34, 42, 48, 40, *Solid Section, orient=orientL9, material=mat9, elset=solidSecL9 1., *Orientation, name=orientL9 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00 3, 5.0 *Material, name=mat9 *Elastic, type=engineering constants 0.390E+07, 0.490E+07, 0.590E+07, 0.190E+00, 0.290E+00, 0.390E+00, 0.390E+06, 0.490E+06 0.590E+06, *Density 0.759E-03, *Plastic 2.0e5, 0. 2.2e5, .0009 2.2e5, .0029 *Boundary side1, 1, 2 *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 5, *El print s, *Node print u, *Output, field *Node output u, *Element output s, *end step ** ** -------------------------------------------------- ** *step, name=step-2 SSD mode-based *steady state dynamics, interval=eigenfrequency 60, *El print, freq=999 s, *Node print, freq=0 *Output, field, op=new *Element output s, *cload side3, 2, 100. *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, 2, 100. *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, 2, 100. *end step ** ** -------------------------------------------------- ** *step, name=step-5 SSD subspace *steady state dynamics, subspace projection, interval=eigenfrequency 60., *cload side3, 2, 100. *end step ** ** -------------------------------------------------- ** *step, name=step-6 random response *random response 60., *select eigenmodes, generate 1, 5 *cload, load case=1 side3, 2, 100. *cload, load case=2 side3, 2, 100. *correlation, psd=white, type=uncorrelated, complex=no 1,1 2,1 *end step ** ** -------------------------------------------------- ** *step, name=step-7 response spectrum *response spectrum, sum=abs, comp=alg spec,0.,1.,0.,1. *select eigenmodes, generate 1, 5 *end step ** ** -------------------------------------------------- ** *step, nlgeom=yes, name=step-8 static *static 1.0, 1.0 *Boundary side3, 1, 1, 0.025 side3, 2, 2, 0.050 *El print, freq=999 s, peeq, *Node print, freq=999 u, *Output, field, op=new *Node output u, *Element output s,peeq *end step ** ** -------------------------------------------------- ** *step, perturbation, name=step-9 static *static *El print, freq=999 s, *Node print, freq=0 *Output, field, op=new *Element output s, *load case, name=load1 *dload eall, bx, 1.e3 *end load case *load case, name=load2 *dload eall, grav, 1.e6, 0.,1.,0. *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