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