*Heading Verification test for s3r elements with distributed layer angle and thickness *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, *Nset, nset=nodeThick1, generate 1, 13, 4, 4, 16, 4, *Nset, nset=nodeThick2, generate 2, 14, 4, 3, 15, 4, *Nodal Thickness nodeThick1, 0.01 nodeThick2, 0.02 **------------------------------------------------------------------------------------------ *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+01, **------------------------------------------------------------------------------------------ *Element, type=s3r, elset=eall 1, 1, 2, 6 2, 2, 3, 7 3, 3, 4, 8 4, 5, 6, 10 5, 6, 7, 11 6, 7, 8, 12 7, 9, 10, 14 8, 10, 11, 15 9, 11, 12, 16 11, 6, 5, 1 12, 7, 6, 2 13, 8, 7, 3 14, 10, 9, 5 15, 11, 10, 6 16, 12, 11, 7 17, 14, 13, 9 18, 15, 14, 10 19, 16, 15, 11 ** *Elset, Elset=shell1 1, 11 *Elset, Elset=shell2 2, 12 *Elset, Elset=shell3 3, 13 *Elset, Elset=shell4 4, 14 *Elset, Elset=shell5 5, 15 *Elset, Elset=shell6 6, 16 *Elset, Elset=shell7 7, 17 *Elset, Elset=shell8 8, 18 *Elset, Elset=shell9 9, 19 ** *Elset, Elset=shellSecL1 shell1, shell2, shell3 *Elset, Elset=shellSecL2 shell4, shell5, shell6 *Elset, Elset=shellSecL3 shell7, shell8, shell9 ** *orientation, name=ori -1., 0., 0., 0., -1., 0. 3, 15. ** *Distribution Table, name=thicknessTable length *Distribution, name=distThickPly1, location=element, table=thicknessTable , 0.150 shell4, 0.175 shell6, 0.165 *Distribution, name=distThickPly2, location=element, table=thicknessTable , 1.400 shell7, 1.250 shell8, 1.100 *Distribution, name=distThickPly4, location=element, table=thicknessTable , 0.130 shell5, 0.095 shell6, 0.075 ** *Distribution Table, name=orientAngTable angle *Distribution, name=distAngPly1, location=element, table=orientAngTable , 10. shell2, 40. shell3, 70. *Distribution, name=distAngPly2, location=element, table=orientAngTable , 30. shell7, 60. shell9, 90. *Distribution, name=distAngPly3, location=element, table=orientAngTable , -20. shell4, -80. shell5, -70. *Distribution, name=distAngPly4, location=element, table=orientAngTable , -40. shell7, -60. shell8, -90. ** *Shell general Section, composite, elset=shellSecL1, nodal thickness, offset=-0.5, orientation=ori 1.000, 3, mat1, distAngPly1, ply-1 1.000, 3, mat1, 78., ply-2 1.000, 3, mat1, 18., ply-3 1.000, 3, mat1, 28., ply-4 ** *Shell general Section, composite, elset=shellSecL2, offset=-0.5, orient=ori distThickPly1, 3, mat1, 38., ply-1 1.250, 3, mat1, 78., ply-2 1.250, 3, mat1, distAngPly3, ply-3 distThickPly4, 3, mat1, 58., ply-4 ** *Shell general Section, composite, elset=shellSecL3, nodal thickness, orientation=ori 1.000, 3, mat1, 68., ply-1 distThickPly2, 3, mat1, distAngPly2, ply-2 1.000, 3, mat1, 48., ply-3 1.000, 3, mat1, distAngPly4, ply-4 *Boundary side1, 1, 6 ** -------------------------------------------------- ** *amplitude, def=smooth step, name=smooth 0., 0., 1., 1.0 *step, nlgeom=yes, name=step-8 explicit *dynamic, explicit , 1.0 *Boundary, amplitude=smooth side3, 3, 3, 0.050 *Output, field *Node output u, *Element output 1, 3, 4, 6, 7, 9, 10, 12 s, sf, ssavg, se, sth, *Element output 5, tshr, *Output, history *Element output,elset=shellSecL1 s11, sf1, ssavg1, se1, sk1, sth, *Element output,elset=shellSecL2 3, tshr13, *end step ** ** -------------------------------------------------- **