*Heading Multistep verification test for stri65 elements using distributions *Node 1, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00 2, 0.8333333E-01, 0.0000000E+00, 0.0000000E+00 3, 0.1666667E+00, 0.0000000E+00, 0.0000000E+00 4, 0.2500000E+00, 0.0000000E+00, 0.0000000E+00 5, 0.3333333E+00, 0.0000000E+00, 0.0000000E+00 6, 0.4166667E+00, 0.0000000E+00, 0.0000000E+00 7, 0.5000000E+00, 0.0000000E+00, 0.0000000E+00 8, 0.0000000E+00, 0.8333333E-01, 0.0000000E+00 9, 0.8333333E-01, 0.8333333E-01, 0.0000000E+00 10, 0.1666667E+00, 0.8333333E-01, 0.0000000E+00 11, 0.2500000E+00, 0.8333333E-01, 0.0000000E+00 12, 0.3333333E+00, 0.8333333E-01, 0.0000000E+00 13, 0.4166667E+00, 0.8333333E-01, 0.0000000E+00 14, 0.5000000E+00, 0.8333333E-01, 0.0000000E+00 15, 0.0000000E+00, 0.1666667E+00, 0.0000000E+00 16, 0.8333333E-01, 0.1666667E+00, 0.0000000E+00 17, 0.1666667E+00, 0.1666667E+00, 0.0000000E+00 18, 0.2500000E+00, 0.1666667E+00, 0.0000000E+00 19, 0.3333333E+00, 0.1666667E+00, 0.0000000E+00 20, 0.4166667E+00, 0.1666667E+00, 0.0000000E+00 21, 0.5000000E+00, 0.1666667E+00, 0.0000000E+00 22, 0.0000000E+00, 0.2500000E+00, 0.0000000E+00 23, 0.8333333E-01, 0.2500000E+00, 0.0000000E+00 24, 0.1666667E+00, 0.2500000E+00, 0.0000000E+00 25, 0.2500000E+00, 0.2500000E+00, 0.0000000E+00 26, 0.3333333E+00, 0.2500000E+00, 0.0000000E+00 27, 0.4166667E+00, 0.2500000E+00, 0.0000000E+00 28, 0.5000000E+00, 0.2500000E+00, 0.0000000E+00 29, 0.0000000E+00, 0.3333333E+00, 0.0000000E+00 30, 0.8333333E-01, 0.3333333E+00, 0.0000000E+00 31, 0.1666667E+00, 0.3333333E+00, 0.0000000E+00 32, 0.2500000E+00, 0.3333333E+00, 0.0000000E+00 33, 0.3333333E+00, 0.3333333E+00, 0.0000000E+00 34, 0.4166667E+00, 0.3333333E+00, 0.0000000E+00 35, 0.5000000E+00, 0.3333333E+00, 0.0000000E+00 36, 0.0000000E+00, 0.4166667E+00, 0.0000000E+00 37, 0.8333333E-01, 0.4166667E+00, 0.0000000E+00 38, 0.1666667E+00, 0.4166667E+00, 0.0000000E+00 39, 0.2500000E+00, 0.4166667E+00, 0.0000000E+00 40, 0.3333333E+00, 0.4166667E+00, 0.0000000E+00 41, 0.4166667E+00, 0.4166667E+00, 0.0000000E+00 42, 0.5000000E+00, 0.4166667E+00, 0.0000000E+00 43, 0.0000000E+00, 0.5000000E+00, 0.0000000E+00 44, 0.8333333E-01, 0.5000000E+00, 0.0000000E+00 45, 0.1666667E+00, 0.5000000E+00, 0.0000000E+00 46, 0.2500000E+00, 0.5000000E+00, 0.0000000E+00 47, 0.3333333E+00, 0.5000000E+00, 0.0000000E+00 48, 0.4166667E+00, 0.5000000E+00, 0.0000000E+00 49, 0.5000000E+00, 0.5000000E+00, 0.0000000E+00 *Nset, nset=side1, generate 1, 7, 1, *Nset, nset=side3, generate 43, 49, 1, *Elset, Elset=eall, generate 1, 18, 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=stri65, elset=eall 1, 1, 3, 17, 2, 10, 9, 2, 17, 15, 1, 16, 8, 9, 3, 3, 5, 19, 4, 12, 11, 4, 19, 17, 3, 18, 10, 11, 5, 5, 7, 21, 6, 14, 13, 6, 21, 19, 5, 20, 12, 13, 7, 15, 17, 31, 16, 24, 23, 8, 31, 29, 15, 30, 22, 23, 9, 17, 19, 33, 18, 26, 25, 10, 33, 31, 17, 32, 24, 25, 11, 19, 21, 35, 20, 28, 27, 12, 35, 33, 19, 34, 26, 27, 13, 29, 31, 45, 30, 38, 37, 14, 45, 43, 29, 44, 36, 37, 15, 31, 33, 47, 32, 40, 39, 16, 47, 45, 31, 46, 38, 39, 17, 33, 35, 49, 34, 42, 41, 18, 49, 47, 33, 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, 1, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, 2, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, 3, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00, 4, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00, 5, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00, 6, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00, 7, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00, 8, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00, 9, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00, 10, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00, 11, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00, 12, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00, 13, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00, 14, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00, 15, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00, 16, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00, 17, 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00, 18, 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.1000000E-01, 1, 0.1000000E-01, 2, 0.1000000E-01, 3, 0.2000000E-01, 4, 0.2000000E-01, 5, 0.3000000E-01, 6, 0.3000000E-01, 7, 0.4000000E-01, 8, 0.4000000E-01, 9, 0.5000000E-01, 10, 0.5000000E-01, 11, 0.6000000E-01, 12, 0.6000000E-01, 13, 0.7000000E-01, 14, 0.7000000E-01, 15, 0.8000000E-01, 16, 0.8000000E-01, 17, 0.9000000E-01, 18, 0.9000000E-01, *Distribution Table, name=offsetTable ratio *Distribution, name=distOffset, location=element, table=offsetTable , -0.5000000E+00, 1, -0.5000000E+00, 2, -0.5000000E+00, 3, -0.4750000E+00, 4, -0.4750000E+00, 5, -0.4500000E+00, 6, -0.4500000E+00, 7, -0.4250000E+00, 8, -0.4250000E+00, 9, -0.4000000E+00, 10, -0.4000000E+00, 11, -0.3750000E+00, 12, -0.3750000E+00, 13, -0.3500000E+00, 14, -0.3500000E+00, 15, -0.3250000E+00, 16, -0.3250000E+00, 17, -0.3000000E+00, 18, -0.3000000E+00, *Distribution Table, name=OriAngleTable angle *Distribution, name=distOriAngle, location=element, table=OriAngleTable , 0.0, 1, -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, 10, -2.5, 11, -3.0, 12, -3.0, 13, -3.5, 14, -3.5, 15, -4.0, 16, -4.0, 17, -4.5, 18, -4.5, *orientation, name=ori distOrient 3, distOriAngle *Shell General Section, composite, elset=eall, orientation=ori, shell thickness=distThick, offset=distOffset 0.1, , mat1, 45 0.1, , mat1, 30 0.1, , mat1, -30 0.1, , mat1, -45 *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, *Node print u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *end step ** ** -------------------------------------------------- ** *step, name=step-2 SSD mode-based *steady state dynamics, interval=eigenfrequency 60, *El print, freq=999 s, sf, ssavg, se, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, *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, 3, 3, 0.050 *El print, freq=999 s, sf, ssavg, se, sth, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, sth, *end step ** ** -------------------------------------------------- ** *step, perturbation, name=step-9 static *static *El print, freq=999 s, sf, ssavg, se, sth, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sf, ssavg, se, 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