*Heading
Multistep verification test for cgax4h elements using distributions
*Node
1, 0.1000000E-01, 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.1000000E-01, 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.1000000E-01, 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.1000000E-01, 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,
*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=cgax4h, 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
, 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