*HEADING
STEADY STATE DYNAMICS,DIRECT VERIFICATION OF CIN3D18R USING DASHPOTS
*NODE,NSET=NALL
1,0.,0.
2,0.,-1.
3,1.,0.
4,1.,-1.
5,0.,0.,1
6,0.,-1.,1.
7,1.,0.,1.
8,1.,-1.,1.
16,0.5,-1.,1.
17,1.,-1.,0.5
18,0.5,-1.,0.
19,0.,-1.,0.5
20,0.5,0.,1.
21,1.,0.,0.5
22,0.5,0.,0.,
23,0.,0.,0.5
24,0.,-0.5,1.
25,1.,-0.5,1.
26,1.,-0.5,0.
27,0.,-0.5,0.
28,0.5,-0.5,0.5
29,0.5,-1.,0.5
30,0.5,0.,0.5
31,0.5,-0.5,1.
32,1.,-0.5,0.5
33,0.5,-0.5,0.
34,0.,-0.5,0.5
1000,0.5,-1.5,0.5
*NSET,NSET=TOP,GENERATE
1,7,2
20,23,1
*ELEMENT,TYPE=C3D27,ELSET=FE
1,6,8,4,2,5,7,3,1,16,17,18,19,
20,21,22,23,24,25,26,27,28,29,
30,31,32,33,34
*SOLID SECTION,ELSET=FE,MATERIAL=ONE
*MATERIAL,NAME=ONE
*ELASTIC
1.,.1
*DENSITY
.01, 
*ELEMENT,TYPE=DASHPOT1,ELSET=DPOT2
101,1000
*DASHPOT,ELSET=DPOT2
2, 
0.1011299,
*ELEMENT,TYPE=DASHPOT1,ELSET=DPOT1
201,1000
*DASHPOT,ELSET=DPOT1
1, 
0.0674199,
*NSET,NSET=BOT
2,4,6,8,16,17,18,19,29
*MPC
TIE,BOT,1000
*RESTART,WRITE,FREQUENCY=1
*STEP
*STEADYSTATEDYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
0.01,0.04,4,1
*BOUNDARY
NALL,1,1
NALL,3,3
TOP,2,2,1.E-4
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*ENERGY PRINT
*EL PRINT,FREQUENCY=0
*NODE PRINT
U, 
PU, 
*OUTPUT,FIELD
*NODE OUTPUT
U, 
*END STEP
*STEP
*STEADYSTATEDYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
0.01,0.04,4,1
*BOUNDARY,OP=NEW
NALL,2,3
TOP,1,1,1.E-4
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*ENERGY PRINT
*EL PRINT,FREQUENCY=0
*NODE PRINT
U, 
PU, 
*OUTPUT,FIELD
*NODE OUTPUT
U, 
*END STEP