*Heading Verification test for s4r 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=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=distTableStiff SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF2,SHELLSTIFF2, SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF3,SHELLSTIFF2, SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF3,SHELLSTIFF3 ** *Distribution, name=distStiff, table=distTableStiff, location=element , 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 shell2, 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 shell3, 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=shellSecL1, density=0.8E+01 distStiff ** *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 ** ** -------------------------------------------------- **