*Heading Verification test for s4r 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=s4r, elset=eall 1, 1, 2, 6, 5, 2, 2, 3, 7, 6, 3, 3, 4, 8, 7, 4, 5, 6, 10, 9, 5, 6, 7, 11, 10, 6, 7, 8, 12, 11, 7, 9, 10, 14, 13, 8, 10, 11, 15, 14, 9, 11, 12, 16, 15, ** *Elset, Elset=shell1 1, *Elset, Elset=shell2 2, *Elset, Elset=shell3 3, *Elset, Elset=shell4 4, *Elset, Elset=shell5 5, *Elset, Elset=shell6 6, *Elset, Elset=shell7 7, *Elset, Elset=shell8 8, *Elset, Elset=shell9 9, ** *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 ** ** -------------------------------------------------- **