*Heading Verification reference for s4 elements with distributed general section stiffness *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=s4, 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, elset=shell1, density=0.8E+01 0.32967E+07, 0.98901E+06, 0.32967E+07, 0.00000E+00, 0.00000E+00, 0.11539E+07, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.27473E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.82418E+03, 0.27473E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.96154E+03 *shell general section, elset=shell2, density=0.8E+01 0.49451E+07, 0.14835E+07, 0.49451E+07, 0.00000E+00, 0.00000E+00, 0.17308E+07, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.92720E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.27816E+04, 0.92720E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.32452E+04 *shell general section, elset=shell3, density=0.8E+01 0.65934E+07, 0.19780E+07, 0.65934E+07, 0.00000E+00, 0.00000E+00, 0.23077E+07, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.21978E+05, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.65934E+04, 0.21978E+05, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.76923E+04 ** *Shell general section, elset=shellSecL2, density=0.8E+01 0.32967E+07, 0.98901E+06, 0.32967E+07, 0.00000E+00, 0.00000E+00, 0.11539E+07, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.27473E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.82418E+03, 0.27473E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.96154E+03 ** *shell general section, elset=shellSecL3, density=0.8E+01 0.49451E+07, 0.14835E+07, 0.49451E+07, 0.00000E+00, 0.00000E+00, 0.17308E+07, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.92720E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.27816E+04, 0.92720E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.32452E+04 ** *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=shellSecL1 3, tshr13, *end step ** ** -------------------------------------------------- **