*Heading Multistep verification test for m3d6 elements using distributions *Preprint, model=yes *Node 1, 0., 0., 0. 2, 0.333333343, 0., 0. 3, 0.666666687, 0., 0. 4, 1., 0., 0. 5, 0., 0.333333343, 0. 6, 0.333333343, 0.333333343, 0. 7, 0.666666687, 0.333333343, 0. 8, 1., 0.333333343, 0. 9, 0., 0.666666687, 0. 10, 0.333333343, 0.666666687, 0. 11, 0.666666687, 0.666666687, 0. 12, 1., 0.666666687, 0. 13, 0., 1., 0. 14, 0.333333343, 1., 0. 15, 0.666666687, 1., 0. 16, 1., 1., 0. 17, 0.166666672, 0., 0. 18, 0.333333343, 0.166666672, 0. 19, 0.166666672, 0.166666672, 0. 20, 0.166666672, 0.333333343, 0. 21, 0., 0.166666672, 0. 22, 0.5, 0., 0. 23, 0.666666687, 0.166666672, 0. 24, 0.5, 0.166666672, 0. 25, 0.5, 0.333333343, 0. 26, 0.833333373, 0., 0. 27, 1., 0.166666672, 0. 28, 0.833333373, 0.166666672, 0. 29, 0.833333373, 0.333333343, 0. 30, 0.333333343, 0.5, 0. 31, 0.166666672, 0.5, 0. 32, 0.166666672, 0.666666687, 0. 33, 0., 0.5, 0. 34, 0.666666687, 0.5, 0. 35, 0.5, 0.5, 0. 36, 0.5, 0.666666687, 0. 37, 1., 0.5, 0. 38, 0.833333373, 0.5, 0. 39, 0.833333373, 0.666666687, 0. 40, 0.333333343, 0.833333373, 0. 41, 0.166666672, 0.833333373, 0. 42, 0.166666672, 1., 0. 43, 0., 0.833333373, 0. 44, 0.666666687, 0.833333373, 0. 45, 0.5, 0.833333373, 0. 46, 0.5, 1., 0. 47, 1., 0.833333373, 0. 48, 0.833333373, 0.833333373, 0. 49, 0.833333373, 1., 0. *Nset, nset=side1 1, 2, 3,4 ,17,22,26 *Nset, nset=side2 4, 8,12,16,27,37,47 *Nset, nset=side3 13,14,15,16,42,46,49 *Nset, nset=side4 1, 5, 9,13,21,33,43 *Nset, nset=load 6,7,10,11 *Elset, Elset=eall, generate 1, 18, 1, *Elset, elset=MembraneSecL1 1, **------------------------------------------------------------------------------------------ *Material, name=mat *Elastic 1.e6, 0.3 *Density 0.1E-03, **------------------------------------------------------------------------------------------ *Element, type=m3d6, elset=eall 1, 1, 2, 6, 17, 18, 19 2, 6, 5, 1, 20, 21, 19 3, 2, 3, 7, 22, 23, 24 4, 7, 6, 2, 25, 18, 24 5, 3, 4, 8, 26, 27, 28 6, 8, 7, 3, 29, 23, 28 7, 5, 6, 10, 20, 30, 31 8, 10, 9, 5, 32, 33, 31 9, 6, 7, 11, 25, 34, 35 10, 11, 10, 6, 36, 30, 35 11, 7, 8, 12, 29, 37, 38 12, 12, 11, 7, 39, 34, 38 13, 9, 10, 14, 32, 40, 41 14, 14, 13, 9, 42, 43, 41 15, 10, 11, 15, 36, 44, 45 16, 15, 14, 10, 46, 40, 45 17, 11, 12, 16, 39, 47, 48 18, 16, 15, 11, 49, 44, 48 *Membrane Section, elset=eall, material=mat, membrane thickness=distThickSect, orientation=ori 0.1 *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, 1, 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, 5, 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, 9, 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, 13, 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, 17, 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00, 2, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, 4, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00, 6, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00, 8, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00, 10, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00, 12, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00, 14, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00, 16, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00, 18, 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., 1, 1., 3, 2., 5, 3., 7, 4., 9, 5., 11, 6., 13, 7., 15, 8., 17, 9., 2, 1., 4, 2., 6, 3., 8, 4., 10, 5., 12, 6., 14, 7., 16, 8., 18, 9., *Distribution, name=distThickSect, location=element, table=thicknessTable , 1. 1, .101 3, .102 5, .103 7, .104 9, .105 11, .106 13, .107 15, .108 17, .109 2, .101 4, .102 6, .103 8, .104 10, .105 12, .106 14, .107 16, .108 18, .109 *Distribution Table, name=orientationTable coord3d, coord3d *Distribution Table, name=oriAngleTable angle *Distribution Table, name=thicknessTable length **------------------------------------------------------------------------------------------ *Initial Conditions, Type=Stress 1, 101., 101. 3, 102., 102. 5, 103., 103. 7, 104., 104. 9, 105., 105. 11, 106., 106. 13, 107., 107. 15, 108., 108. 17, 109., 109. 2, 101., 101. 4, 102., 102. 6, 103., 103. 8, 104., 104. 10, 105., 105. 12, 106., 106. 14, 107., 107. 16, 108., 108. 18, 109., 109. *Boundary side1, 1, 3 side2, 1, 3 side3, 1, 3 side4, 1, 3 *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-0, nlgeom Resolve initial stresses *static *El print s, *Node print u, *Output, field *Node output u, *Element output s, sth, *end step ** ** -------------------------------------------------- ** *step, name=step-1 frequency analysis *frequency 75, *Node print u, *Output, field *Node output u, *Element output sth *end step ** ** -------------------------------------------------- ** *step, name=step-2 SSD mode-based *steady state dynamics, interval=eigenfrequency 60. *El print, freq=999 s, sth, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sth, *Cload load, 1,1. load, 3,1. *Cload, imaginary load, 2,1. *select eigenmodes, generate 1, 75 *modal damping, rayleigh 1, 5, 0.5 *end step ** ** -------------------------------------------------- ** *step, name=step-3 modal dynamic *modal dynamic 0.1,1. *Cload load, 1,1. load, 2,1. load, 3,1. *select eigenmodes, generate 1, 75 *modal damping, rayleigh 1, 5, 0.5 *end step ** ** -------------------------------------------------- ** *step, name=step-4 SSD direct *steady state dynamics, direct, interval=eigenfrequency 60., *Cload load, 1,1. load, 3,1. *Cload, imaginary load, 2,1. *end step ** ** -------------------------------------------------- ** *step, name=step-5 SSD subspace *steady state dynamics, subspace projection, interval=eigenfrequency 60., *Cload load, 1,1. load, 3,1. *Cload, imaginary load, 2,1. *end step ** ** -------------------------------------------------- ** *step, name=step-6 random response *random response 60., *select eigenmodes, generate 1, 75 *cload, load case=1 load, 1, 1. load, 3, 1. *cload, load case=2 load, 2, 1. *correlation, psd=white, type=uncorrelated, complex=no 1,1 2,1 *Node print, freq=999 u, *Output, field *Node output u, *end step ** ** -------------------------------------------------- ** *step, name=step-7 response spectrum *response spectrum, sum=abs, comp=alg spec,1.,1.,1.,1. *select eigenmodes, generate 1, 75 *end step ** ** -------------------------------------------------- ** *step, nlgeom=yes, name=step-8 static *static 1.0, 1.0 *Restart, write *Boundary, op=new side1, 3, 3 side2, 3, 3 side3,3,3 side4,3,3 side1, 2, 2 side4, 1, 1 side2, 1,1,.025 side3, 2,2,.025 *El print, freq=999 s, sth, *Node print, freq=999 u, *Output, field *Node output u, *Element output s, sth, *Output, history *Element output,elset=MembraneSecL1 s11, sth, *end step ** ** -------------------------------------------------- ** *step, perturbation, name=step-9 static *static *El print, freq=999 s, sth, *Node print, freq=999 u, *load case, name=load1 *dload eall, bz, 1.e3 *Output, field *Node output u, *Element output s, sth, *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 *Output, field *Element output s, sth, *end step