*Heading Multistep verification test for s9r5 elements using distributions *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, **------------------------------------------------------------------------------------------ *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=s9r5, elset=eall 1, 1, 3, 17, 15, 2, 10, 16, 8, 9 2, 3, 5, 19, 17, 4, 12, 18, 10, 11 3, 5, 7, 21, 19, 6, 14, 20, 12, 13 4, 15, 17, 31, 29, 16, 24, 30, 22, 23 5, 17, 19, 33, 31, 18, 26, 32, 24, 25 6, 19, 21, 35, 33, 20, 28, 34, 26, 27 7, 29, 31, 45, 43, 30, 38, 44, 36, 37 8, 31, 33, 47, 45, 32, 40, 46, 38, 39 9, 33, 35, 49, 47, 34, 42, 48, 40, 41 *Distribution Table, name=orientationTable coord3d, coord3d *Distribution, name=distOrient, location=element, table=orientationTable , 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, 2, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00, 3, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00, 4, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00, 5, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00, 6, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00, 7, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00, 8, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00, 9, 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00, *Distribution Table, name=thicknessTable length *Distribution, name=distThick, location=element, table=thicknessTable , 0.2000000E-01, 1, 0.1000000E-01, 3, 0.3000000E-01, 4, 0.4000000E-01, 5, 0.5000000E-01, 6, 0.6000000E-01, 7, 0.7000000E-01, 8, 0.8000000E-01, 9, 0.9000000E-01, *Distribution Table, name=offsetTable ratio *Distribution, name=distOffset, location=element, table=offsetTable , -0.4250000E+00, 1, -0.5000000E+00, 2, -0.4750000E+00, 3, -0.4500000E+00, 5, -0.4000000E+00, 6, -0.3750000E+00, 7, -0.3500000E+00, 8, -0.3250000E+00, 9, -0.3000000E+00, *Distribution Table, name=OriAngleTable angle *Distribution, name=distOriAngle, location=element, table=OriAngleTable , -0.5 2, -0.5 3, -1.0 4, -1.0 5, -1.5 6, -1.5, 7, -2.0, 8, -2.0, 9, -2.5, *orientation, name=ori distOrient 3, distOriAngle *Shell Section, composite, elset=eall, orientation=ori, offset=distOffset, shell thickness=distThick distThickPly1, 3, mat1, 90., ply-1 1., 3, mat1, distAnglePly2, ply-2 *Distribution, name=distThickPly1, location=element, table=thicknessTable , 0.1000000E-01 eall, 1.0 *Distribution, name=distAnglePly2, location=element, table=OriAngleTable , 0.1000000E-01 eall, 0. *Boundary side1, 1, 6 *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, 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-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 3, tshr, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *Element output 3, tshr, *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