*Heading
Multistep verification test for cgax8r elements using distributions
*Node
1, 0.1000000E-01, 0.0000000E+00, 0.0000000E+00
2, 0.1666667E+00, 0.0000000E+00, 0.0000000E+00
3, 0.3333333E+00, 0.0000000E+00, 0.0000000E+00
4, 0.5000000E+00, 0.0000000E+00, 0.0000000E+00
5, 0.6666667E+00, 0.0000000E+00, 0.0000000E+00
6, 0.8333333E+00, 0.0000000E+00, 0.0000000E+00
7, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00
8, 0.1000000E-01, 0.1666667E+00, 0.0000000E+00
9, 0.1666667E+00, 0.1666667E+00, 0.0000000E+00
10, 0.3333333E+00, 0.1666667E+00, 0.0000000E+00
11, 0.5000000E+00, 0.1666667E+00, 0.0000000E+00
12, 0.6666667E+00, 0.1666667E+00, 0.0000000E+00
13, 0.8333333E+00, 0.1666667E+00, 0.0000000E+00
14, 0.1000000E+01, 0.1666667E+00, 0.0000000E+00
15, 0.1000000E-01, 0.3333333E+00, 0.0000000E+00
16, 0.1666667E+00, 0.3333333E+00, 0.0000000E+00
17, 0.3333333E+00, 0.3333333E+00, 0.0000000E+00
18, 0.5000000E+00, 0.3333333E+00, 0.0000000E+00
19, 0.6666667E+00, 0.3333333E+00, 0.0000000E+00
20, 0.8333333E+00, 0.3333333E+00, 0.0000000E+00
21, 0.1000000E+01, 0.3333333E+00, 0.0000000E+00
22, 0.1000000E-01, 0.5000000E+00, 0.0000000E+00
23, 0.1666667E+00, 0.5000000E+00, 0.0000000E+00
24, 0.3333333E+00, 0.5000000E+00, 0.0000000E+00
25, 0.5000000E+00, 0.5000000E+00, 0.0000000E+00
26, 0.6666667E+00, 0.5000000E+00, 0.0000000E+00
27, 0.8333333E+00, 0.5000000E+00, 0.0000000E+00
28, 0.1000000E+01, 0.5000000E+00, 0.0000000E+00
29, 0.1000000E-01, 0.6666667E+00, 0.0000000E+00
30, 0.1666667E+00, 0.6666667E+00, 0.0000000E+00
31, 0.3333333E+00, 0.6666667E+00, 0.0000000E+00
32, 0.5000000E+00, 0.6666667E+00, 0.0000000E+00
33, 0.6666667E+00, 0.6666667E+00, 0.0000000E+00
34, 0.8333333E+00, 0.6666667E+00, 0.0000000E+00
35, 0.1000000E+01, 0.6666667E+00, 0.0000000E+00
36, 0.1000000E-01, 0.8333333E+00, 0.0000000E+00
37, 0.1666667E+00, 0.8333333E+00, 0.0000000E+00
38, 0.3333333E+00, 0.8333333E+00, 0.0000000E+00
39, 0.5000000E+00, 0.8333333E+00, 0.0000000E+00
40, 0.6666667E+00, 0.8333333E+00, 0.0000000E+00
41, 0.8333333E+00, 0.8333333E+00, 0.0000000E+00
42, 0.1000000E+01, 0.8333333E+00, 0.0000000E+00
43, 0.1000000E-01, 0.1000000E+01, 0.0000000E+00
44, 0.1666667E+00, 0.1000000E+01, 0.0000000E+00
45, 0.3333333E+00, 0.1000000E+01, 0.0000000E+00
46, 0.5000000E+00, 0.1000000E+01, 0.0000000E+00
47, 0.6666667E+00, 0.1000000E+01, 0.0000000E+00
48, 0.8333333E+00, 0.1000000E+01, 0.0000000E+00
49, 0.1000000E+01, 0.1000000E+01, 0.0000000E+00
*Nset, nset=side1, generate
1, 7, 1,
*Nset, nset=side3, generate
43, 49, 1,
*Elset, Elset=eall, generate
1, 9, 1,
**------------------------------------------------------------------------------------------
*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-03,
**------------------------------------------------------------------------------------------
*Element, type=cgax8r, elset=eall
1, 1, 3, 17, 15, 2, 10, 16, 8,
2, 3, 5, 19, 17, 4, 12, 18, 10,
3, 5, 7, 21, 19, 6, 14, 20, 12,
4, 15, 17, 31, 29, 16, 24, 30, 22,
5, 17, 19, 33, 31, 18, 26, 32, 24,
6, 19, 21, 35, 33, 20, 28, 34, 26,
7, 29, 31, 45, 43, 30, 38, 44, 36,
8, 31, 33, 47, 45, 32, 40, 46, 38,
9, 33, 35, 49, 47, 34, 42, 48, 40,
*Distribution table, name=OrientationTable
coord3d, coord3d
*Distribution, name=distOrient, location=element, table=orientationTable
, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00,
1, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00,
2, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00,
3, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00,
4, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00,
5, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00,
6, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00,
7, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00,
8, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00,
9, 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00,
*Distribution Table, name=AngleTable
angle
*Distribution, name=distAngle, table=AngleTable, location=element
, 0.0
1, 1.0
2, -1.0
3, 2.0
4, -2.0
5, 3.0
6, -3.0
7, 4.0
8, -4.0
9, 5.0
*Solid Section, material=mat1, elset=eall, orientation=ori
1.,
*Orientation, name=ori
distOrient
3, distAngle
*Boundary
side1, 1, 2
side1, 5, 5
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-1
static
*static
.25, 1.0
*Boundary
side3, 1, 1, 0.01
side3, 5, 5, 0.01
*El print, freq=999
s,
*Node print, freq=999
u,
*Output, field, op=new
*Node output
u,
*Element output
s,
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-2
static
*static
*El print, freq=999
s,
*Node print, freq=0
*Output, field, op=new
*Element output
s,
*load case, name=load1
*dload
eall, br, 1.e3
*end load case
*load case, name=load2
*dload
eall, grav, 1.e6, 0.,1.,0.
*end load case
*load case, name=load3
*dload
eall, cent, 100.,0.,0.,0.,1.,0.
*end load case
*load case, name=load4
*dload
eall, centrif, 1.e6,0.,0.,0.,1.,0.
*end load case
*end step