*HEADING
360 DEGREES MODEL CORRESPONDING TO 1X2 MESH OF C3D8R ELEMENTS.
*PREPRINT, ECHO=YES, MODEL=YES, HISTORY=YES
**
*NODE, NSET=NALL
1, -5., -5., 10.
2, -5., -5., 2.71932E-29
3, -25., -25., 7.69139E-30
4, -25., -25., 10.
5, 0., -10., 0.
6, 0., -50., 0.
7, 0., -50., 10.
8, 0., -10., 10.
9, 5., 5., 2.75941E-29
10, 5., 5., 10.
11, 25., 25., 10.
12, 25., 25., 1.99023E-29
13, 10., 0., 10.
14, 50., 0., 10.
15, 50., 0., 0.
16, 10., 0., 0.
17, 25., -25., 3.54987E-30
18, 25., -25., 10.
19, 5., -5., 10.
20, 5., -5., 2.09995E-29
21, -5., 5., 2.11967E-29
22, -25., 25., 3.74709E-30
23, -25., 25., 10.
24, -5., 5., 10.
25, -10., 0., 0.
26, -50., 0., 0.
27, -50., 0., 10.
28, -10., 0., 10.
29, 0., 10., 0.
30, 0., 50., 0.
31, 0., 50., 10.
32, 0., 10., 10.
**
*ELEMENT, TYPE=C3D8R, ELSET=EALL
1, 8, 5, 6, 7, 1, 2, 3, 4
2, 16, 13, 14, 15, 9, 10, 11, 12
3, 15, 14, 13, 16, 17, 18, 19, 20
4, 25, 26, 27, 28, 21, 22, 23, 24
5, 21, 22, 23, 24, 29, 30, 31, 32
6, 1, 2, 3, 4, 28, 25, 26, 27
7, 17, 18, 19, 20, 6, 7, 8, 5
8, 9, 10, 11, 12, 29, 32, 31, 30
**
*NSET, NSET=NBC
1, 2, 5, 8, 9, 10, 13, 16, 19, 20, 21, 24, 25, 28, 29, 32
**
*SOLID SECTION, ELSET=EALL,MATERIAL=MAT
**
*MATERIAL, NAME=MAT
*DENSITY
8.22E-9
*ELASTIC
205000, 0.3, 293.149
202000, 0.3, 373.149
196000, 0.3, 473.149
190000, 0.3, 573.149
184000, 0.3, 673.150
177000, 0.3, 773.150
170000, 0.3, 873.150
163000, 0.3, 973.150
*PLASTIC
900.0, 0.000000
950.0, 0.000529
1000.0, 0.001611
1050.0, 0.004433
1100.0, 0.012430
1150.0, 0.035075
1200.0, 0.097323
1217.0, 0.136826
*CREEP,LAW=TIME
4.3514E-14, 4.3286, -0.2442, 1250.0
2.4022E-10, 3.3482, -0.2290, 1365.0
3.2310E-06, 2.0677, -0.3762, 1480.0
*EXPANSION
1.300E-5, 293.149
1.330E-5, 373.149
1.362E-5, 473.149
1.399E-5, 573.150
1.432E-5, 673.150
1.470E-5, 773.150
1.512E-5, 873.150
1.571E-5, 973.150
**
*BOUNDARY
NBC, 1,3
**
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL, 293.149
**
*************************************************************************
** Static step
**
*STEP, INC=100, NLGEOM
*STATIC
0.1, 1.0
**
*TEMPERATURE
NALL, 1350.0
**
*DLOAD
EALL, CENT, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0
**
*OUTPUT, FIELD, FREQ=1
*ELEMENT OUTPUT
S, LE, PE
*NODE OUTPUT
U, RF
**
*OUTPUT, HISTORY, FREQ=1
*ENERGY OUTPUT
ALLAE, ALLCD, ALLIE, ALLPD, ALLSE, ALLVD, ALLWK, ETOTAL
**
*NODE PRINT
U
*EL PRINT, FREQ=0
**
*CONTROLS, PARAMETERS=FIELD
,0.0,
*END STEP
*************************************************************************
**
*STEP, INC=1000, NLGEOM
*VISCO, CETOL=3.0E-3
0.01, 50.0, 0.00001, 10.0
**
*OUTPUT, FIELD, FREQ=1
*ELEMENT OUTPUT
S, LE, PE, CE
*NODE OUTPUT
U, RF
**
*OUTPUT, HISTORY, FREQ=1
*ENERGY OUTPUT
ALLAE, ALLCD, ALLIE, ALLPD, ALLSE, ALLVD, ALLWK, ETOTAL
**
*node print,freq=99
u
*END STEP
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
22
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
*STEP
*STEADY STATE DYNAMICS
30.
*MODAL DAMPING
1,21,0.01
*SELECT EIGENMODES, GENERATE
1,21,1
*DLOAD,LOAD CASE=1
2, GRAV, 10.,0.,0.,1.
8, GRAV, 10.,0.,0.,1.
1, GRAV,-10.,0.,0.,1.
6, GRAV,-10.,0.,0.,1.
*DLOAD,LOAD CASE=2
4, GRAV,-10.,0.,0.,1.
5, GRAV,-10.,0.,0.,1.
3, GRAV, 10.,0.,0.,1.
7, GRAV, 10.,0.,0.,1.
*NODE PRINT
U
*OUTPUT,FIELD
U
*MODAL PRINT
GU
*END STEP