*Heading
Multistep verification test for s9r5 elements using distributions
*Node
1, 0.0000000E+00, 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.0000000E+00, 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.0000000E+00, 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.0000000E+00, 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.0000000E+00, 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.0000000E+00, 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.0000000E+00, 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=s9r5, elset=eall
1, 1, 3, 17, 15, 2, 10, 16, 8, 9
2, 3, 5, 19, 17, 4, 12, 18, 10, 11
3, 5, 7, 21, 19, 6, 14, 20, 12, 13
4, 15, 17, 31, 29, 16, 24, 30, 22, 23
5, 17, 19, 33, 31, 18, 26, 32, 24, 25
6, 19, 21, 35, 33, 20, 28, 34, 26, 27
7, 29, 31, 45, 43, 30, 38, 44, 36, 37
8, 31, 33, 47, 45, 32, 40, 46, 38, 39
9, 33, 35, 49, 47, 34, 42, 48, 40, 41
*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=thicknessTable
length
*Distribution, name=distThick, location=element, table=thicknessTable
, 0.1000000E-01,
2, 0.2000000E-01,
3, 0.3000000E-01,
4, 0.4000000E-01,
5, 0.5000000E-01,
6, 0.6000000E-01,
7, 0.7000000E-01,
8, 0.8000000E-01,
9, 0.9000000E-01,
*Distribution Table, name=offsetTable
ratio
*Distribution, name=distOffset, location=element, table=offsetTable
, -0.5000000E+00,
2, -0.4750000E+00,
3, -0.4500000E+00,
4, -0.4250000E+00,
5, -0.4000000E+00,
6, -0.3750000E+00,
7, -0.3500000E+00,
8, -0.3250000E+00,
9, -0.3000000E+00,
*Distribution Table, name=OriAngleTable
angle
*Distribution, name=distOriAngle, location=element, table=OriAngleTable
, 0.0,
1, -0.5
2, -0.5
3, -1.0
4, -1.0
5, -1.5
6, -1.5,
7, -2.0,
8, -2.0,
9, -2.5,
*orientation, name=ori
distOrient
3, distOriAngle
*Shell General Section, material=mat1, elset=eall, orientation=ori, offset=distOffset, shell thickness=distThick
0.0100,
*Boundary
side1, 1, 6
*spectrum, name=spec,type=acceleration
0., 0., 0.
1.e5, 10., 0.
1.e5, 100., 0.
**
** --------------------------------------------------
**
*step, name=step-1
frequency analysis
*frequency
5,
*El print
s,
sf,
ssavg,
se,
*El print
2,
tshr,
*Node print
u,
*Output, field
*Node output
u,
*Element output
s,
sf,
ssavg,
se,
sth,
*Element output
2,
tshr,
*end step
**
** --------------------------------------------------
**
*step, name=step-2
SSD mode-based
*steady state dynamics, interval=eigenfrequency
60,
*El print, freq=999
s,
sf,
ssavg,
se,
sth,
*El print, freq=999
2,
tshr,
*Node print, freq=999
u,
*Output, field
*Node output
u,
*Element output
s,
sf,
ssavg,
se,
sth,
*Element output
2,
tshr,
*cload
side3, 4, 0.1
side3, 3, 2.0
side3, 2, 0.5
*select eigenmodes, generate
1, 5
*modal damping, rayleigh
1, 5, 0.5
*end step
**
** --------------------------------------------------
**
*step, name=step-3
modal dynamic
*modal dynamic
0.1,1.
*cload
side3, 4, 0.1
side3, 3, 0.2
side3, 2, 0.5
*select eigenmodes, generate
1, 5
*modal damping, rayleigh
1, 5, 0.5
*end step
**
** --------------------------------------------------
**
*step, name=step-4
SSD direct
*steady state dynamics, direct, interval=eigenfrequency
60.,
*cload
side3, 4, 0.1
side3, 3, 0.2
side3, 2, 0.5
*end step
**
** --------------------------------------------------
**
*step, name=step-7
response spectrum
*response spectrum, sum=abs, comp=alg
spec,0.,0.,1.,1.
*select eigenmodes, generate
1, 5
*end step
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-8
static
*static
1.0, 1.0
*Boundary
side3, 2, 2, 0.025
side3, 3, 3, 0.050
*El print, freq=999
s,
sf,
ssavg,
se,
sth,
*El print, freq=999
2,
tshr,
*Node print, freq=999
u,
*Output, field
*Node output
u,
*Element output
s,
sf,
ssavg,
se,
sth,
*Element output
2,
tshr,
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-9
static
*static
*El print, freq=999
s,
sf,
ssavg,
se,
sth,
*El print, freq=999
2,
tshr,
*Node print, freq=999
u,
*Output, field
*Node output
u,
*Element output
s,
sf,
ssavg,
se,
sth,
*Element output
2,
tshr,
*load case, name=load1
*dload
eall, bz, 1.e3
*end load case
*load case, name=load2
*dload
eall, grav, 1.e6, 0.,0.,1.
*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