*Heading Multistep verification test for s3 elements using distributions *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, 18, 1, *Elset, elset=shellSecL1 1, **------------------------------------------------------------------------------------------ *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, *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, *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, *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, **------------------------------------------------------------------------------------------ *Element, type=s3, elset=e1 1, 1, 2, 6, 2, 6, 5, 1, *Element, type=s3, elset=e2 3, 2, 3, 7, 4, 7, 6, 2, *Element, type=s3, elset=e3 5, 3, 4, 8, 6, 8, 7, 3, *Element, type=s3, elset=e4 7, 5, 6, 10, 8, 10, 9, 5, *Element, type=s3, elset=e5 9, 6, 7, 11, 10, 11, 10, 6, *Element, type=s3, elset=e6 11, 7, 8, 12, 12, 12, 11, 7, *Element, type=s3, elset=e7 13, 9, 10, 14, 14, 14, 13, 9, *Element, type=s3, elset=e8 15, 10, 11, 15, 16, 15, 14, 10, *Element, type=s3, elset=e9 17, 11, 12, 16, 18, 16, 15, 11, *Shell Section, composite, elset=eall, offset=distOffset, shell thickness=distThickSect, orientation=ori,layup=layup-1 distThickPly1, 3, mat1, distAnglePly1, ply-1 distThickPly2, 3, mat2, -60, ply-2 .01253, 3, mat3, distAnglePly3, ply-3 distThickPly4, 3, mat4, distAnglePly4, ply-4 *Orientation, name=ori distOrient 3, distOrientAngle *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, e1, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, e2, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00, e3, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00, e4, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00, e5, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00, e6, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00, e7, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00, e8, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00, e9, 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00, *Distribution, name=distOrientAngle, location=element, table=oriAngleTable , 0., e1, 1., e2, 2., e3, 3., e4, 4., e5, 5., e6, 6., e7, 7., e8, 8., e9, 9., *Distribution, name=distThickSect, location=element, table=thicknessTable , 1. e1, 2.2110E-02 e2, 3.1053E-02 e3, 3.8610E-02 e4, 5.4483E-02 e5, 5.5110E-02 e6, 6.4053E-02 e7, 7.1610E-02 e8, 8.7483E-02 e9, 8.8110E-02 *Distribution, name=distOffset, location=element, table=offsetTable , 0. e1, -0.5000000E+00, e2, -0.4750000E+00, e3, -0.4500000E+00, e4, -0.4250000E+00, e5, -0.4000000E+00, e6, -0.3750000E+00, e7, -0.3500000E+00, e8, -0.3250000E+00, e9, -0.3000000E+00, *Distribution, name=distThickPly1, location=element, table=thicknessTable , 0.00251 e2, 0.0051 e3, 0.00751 e4, 0.011 e5, 0.01251 e6, 0.0151 e7, 0.01751 e8, 0.021 e9, 0.02251 *Distribution, name=distThickPly2, location=element, table=thicknessTable , .02252 e1, 0.00252 e2, 0.0052 e3, 0.00752 e4, 0.012 e5, 0.01252 e6, 0.0152 e7, 0.01752 e8, 0.022 *Distribution, name=distThickPly4, location=element, table=thicknessTable , .02254 e1, 0.00254 e2, 0.0054 e3, 0.00754 e4, 0.014 e5, 0.01254 e6, 0.0154 e7, 0.01754 e8, 0.024 e9, 0.02254 *Distribution, name=distAnglePly1, location=element, table=oriAngleTable , 60. e2, 60.2 e3, 60.3 e4, 60.4 e5, 60.5 e6, 60.6 e7, 60.7 e8, 60.8 e9, 0. *Distribution, name=distAnglePly3, location=element, table=oriAngleTable , 0. *Distribution, name=distAnglePly4, location=element, table=oriAngleTable , 0. e1, 90. e2, 90.2 e3, 90.3 e4, 90.4 e5, 90.5 e6, 90.6 e7, 90.7 e8, 90.8 e9, 0. *Distribution Table, name=orientationTable coord3d, coord3d *Distribution Table, name=oriAngleTable angle *Distribution Table, name=thicknessTable length *Distribution Table, name=offsetTable ratio *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 5, *El print s, sf, ssavg, se, *El print 3, tshr, *Node print u, *Output, field *Node output u, *Element output s, sf, ssavg, 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, 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, *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, 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, 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