*Heading
Multistep verification test for s4 elements with nodal 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,
*Elset, Elset=shellSecL1, generate
1, 7, 3,
*Elset, Elset=shellSecL2, generate
2, 8, 3,
*Elset, Elset=shellSecL3, generate
3, 9, 3,
*elset, elset=sideEl
shellSecL1, shellSecL3
*elset, elset=midEl
shellSecL2,
*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,
** 0.750E-03,
**------------------------------------------------------------------------------------------
*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,
*Distribution Table, name=orientationTable
coord3d, coord3d
*Distribution, name=distOrient, location=element, table=orientationTable
, 1., 0., 0., 0., 1.0, 0.0
shellSecL1, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00
shellSecL2, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00, 0.0000000E+00
shellSecL3, -0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00
*orientation, name=ori
distOrient
3, 0.0
*Distribution Table, name=ThicknessTable
length
*Distribution, name=distThick, location=element, table=thicknessTable
, 1.0
shellSecL2, 0.01
*Shell Section, material=mat1, elset=midEl, offset=-0.50000E+00, orient=ori, shell thickness=distThick
1.0,
*shell section, material=mat1, elset=sideEl, nodal thickness, offset=-0.5, orientation=ori
1.,
*Boundary
side1, 1, 6
**
** --------------------------------------------------
**
**
*Elset, Elset =QA_TEST
eall,
*Nset, Nset =QA_TEST
Side3,
**
*step, nlgeom=yes, name=step-8
dynamic
*Dynamic, Explicit
, 1.0
*Amplitude, Name=Smooth, Def=Smooth Step
0., 0., 0.5, 1.0, 1., 1.
*Boundary, Amp=Smooth
side3, 3, 3, 0.050
*Output, field, op=new, Time Int = 0.25, Time Marks =Yes
*Node output, Nset =QA_TEST
RF,
*Node output
u,
*Element output, Elset =QA_Test
s,
sf,
ssavg,
se,
sth,
*Element output, Elset =QA_Test
3,
tshr,
*Output, hist, filter=antialiasing, Time Int=0.1
*Element Output, Elset=QA_TEST
s11,
sf1,
ssavg1,
se1,
sk1,
sth,
*Element output, Elset=QA_TEST
3,
tshr13,
*Node Output, Nset=QA_TEST
RF3,
*Energy output, VAR=PRESELECT
*end step
**