*Heading Verification reference for s4rs 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=s4rs, 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. ** *Shell general Section, composite, elset=shell1, nodal thickness, offset=-0.5, orientation=ori 1.000, 3, mat1, 10., 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=shell2, nodal thickness, offset=-0.5, orientation=ori 1.000, 3, mat1, 40., 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=shell3, nodal thickness, offset=-0.5, orientation=ori 1.000, 3, mat1, 70., 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=shell4, offset=-0.5, orient=ori 0.175, 3, mat1, 38., ply-1 1.250, 3, mat1, 78., ply-2 1.250, 3, mat1, -80., ply-3 0.130, 3, mat1, 58., ply-4 *Shell general Section, composite, elset=shell5, offset=-0.5, orient=ori 0.150, 3, mat1, 38., ply-1 1.250, 3, mat1, 78., ply-2 1.250, 3, mat1, -70., ply-3 0.095, 3, mat1, 58., ply-4 *Shell general Section, composite, elset=shell6, offset=-0.5, orient=ori 0.165, 3, mat1, 38., ply-1 1.250, 3, mat1, 78., ply-2 1.250, 3, mat1, -20., ply-3 0.075, 3, mat1, 58., ply-4 ** *Shell general Section, composite, elset=shell7, nodal thickness, orientation=ori 1.000, 3, mat1, 68., ply-1 1.250, 3, mat1, 60., ply-2 1.000, 3, mat1, 48., ply-3 1.000, 3, mat1, -60., ply-4 *Shell general Section, composite, elset=shell8, nodal thickness, orientation=ori 1.000, 3, mat1, 68., ply-1 1.100, 3, mat1, 30., ply-2 1.000, 3, mat1, 48., ply-3 1.000, 3, mat1, -90., ply-4 *Shell general Section, composite, elset=shell9, nodal thickness, orientation=ori 1.000, 3, mat1, 68., ply-1 1.400, 3, mat1, 90., ply-2 1.000, 3, mat1, 48., ply-3 1.000, 3, mat1, -40., 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 ** ** -------------------------------------------------- **