*HEADING
 block-4.inp: A 10x10x100 inch concrete column modelled as a stack of
 10 inch cubes with interface elements in between.  It has horizontal
 support at the top and is subjected to the ground motion QUAKE.AMP from
 the ABAQUS example set. There is a 1 psi vertical surcharge at the top.
 Units: pounds, inches, seconds.
**
********************************** Nodes
**
*NODE,NSET=INTER0A
1,0.,0.,0.0
2,10.0,0.,0.0
3,10.0,10.0,0.0
4,0.,10.0,0.0
*NODE,NSET=INTER0B
5,0.,0.,0.0
6,10.0,0.,0.0
7,10.0,10.0,0.0
8,0.,10.0,0.0
*NODE,NSET=INTER1A
9,0.,0.,10.0
10,10.0,0.,10.0
11,10.0,10.0,10.0
12,0.,10.0,10.0
*NODE,NSET=INTER1B
13,0.,0.,10.0
14,10.0,0.,10.0
15,10.0,10.0,10.0
16,0.,10.0,10.0
*NODE,NSET=INTER2A
17,0.,0.,20.0
18,10.0,0.,20.0
19,10.0,10.0,20.0
20,0.,10.0,20.0
*NODE,NSET=INTER2B
21,0.,0.,20.0
22,10.0,0.,20.0
23,10.0,10.0,20.0
24,0.,10.0,20.0
*NODE,NSET=INTER3A
25,0.,0.,30.0
26,10.0,0.,30.0
27,10.0,10.0,30.0
28,0.,10.0,30.0
*NODE,NSET=INTER3B
29,0.,0.,30.0
30,10.0,0.,30.0
31,10.0,10.0,30.0
32,0.,10.0,30.0
*NODE,NSET=INTER4A
33,0.,0.,40.0
34,10.0,0.,40.0
35,10.0,10.0,40.0
36,0.,10.0,40.0
*NODE,NSET=INTER4B
37,0.,0.,40.0
38,10.0,0.,40.0
39,10.0,10.0,40.0
40,0.,10.0,40.0
*NODE,NSET=INTER5A
41,0.,0.,50.0
42,10.0,0.,50.0
43,10.0,10.0,50.0
44,0.,10.0,50.0
*NODE,NSET=INTER5B
45,0.,0.,50.0
46,10.0,0.,50.0
47,10.0,10.0,50.0
48,0.,10.0,50.0
*NODE,NSET=INTER6A
49,0.,0.,60.0
50,10.0,0.,60.0
51,10.0,10.0,60.0
52,0.,10.0,60.0
*NODE,NSET=INTER6B
53,0.,0.,60.0
54,10.0,0.,60.0
55,10.0,10.0,60.0
56,0.,10.0,60.0
*NODE,NSET=INTER7A
57,0.,0.,70.0
58,10.0,0.,70.0
59,10.0,10.0,70.0
60,0.,10.0,70.0
*NODE,NSET=INTER7B
61,0.,0.,70.0
62,10.0,0.,70.0
63,10.0,10.0,70.0
64,0.,10.0,70.0
*NODE,NSET=INTER8A
65,0.,0.,80.0
66,10.0,0.,80.0
67,10.0,10.0,80.0
68,0.,10.0,80.0
*NODE,NSET=INTER8B
69,0.,0.,80.0
70,10.0,0.,80.0
71,10.0,10.0,80.0
72,0.,10.0,80.0
*NODE,NSET=INTER9A
73,0.,0.,90.0
74,10.0,0.,90.0
75,10.0,10.0,90.0
76,0.,10.0,90.0
*NODE,NSET=INTER9B
77,0.,0.,90.0
78,10.0,0.,90.0
79,10.0,10.0,90.0
80,0.,10.0,90.0
*NODE,NSET=INTER10A
81,0.,0.,100.0
82,10.0,0.,100.0
83,10.0,10.0,100.0
84,0.,10.0,100.0
*NODE,NSET=INTER10B
85,0.,0.,100.0
86,10.0,0.,100.0
87,10.0,10.0,100.0
88,0.,10.0,100.0
** Use an 8-node brick element, reduced integration with hourglass control. 
*ELEMENT, TYPE=C3D8RH, ELSET=BLOCK 
1,5,6,7,8,9,10,11,12
2,13,14,15,16,17,18,19,20
3,21,22,23,24,25,26,27,28
4,29,30,31,32,33,34,35,36
5,37,38,39,40,41,42,43,44
6,45,46,47,48,49,50,51,52
7,53,54,55,56,57,58,59,60
8,61,62,63,64,65,66,67,68
9,69,70,71,72,73,74,75,76
10,77,78,79,80,81,82,83,84
** Use an 8-node interface element 
*ELEMENT, TYPE=INTER4, ELSET=INTER 
1001,1,2,3,4,5,6,7,8
1002,9,10,11,12,13,14,15,16
1003,17,18,19,20,21,22,23,24
1004,25,26,27,28,29,30,31,32
1005,33,34,35,36,37,38,39,40
1006,41,42,43,44,45,46,47,48
1007,49,50,51,52,53,54,55,56
1008,57,58,59,60,61,62,63,64
1009,65,66,67,68,69,70,71,72
1010,73,74,75,76,77,78,79,80
********************************** Materials
**
** Concrete block.
*SOLID SECTION,ELSET=BLOCK,MATERIAL=MASONRY
*MATERIAL,NAME=MASONRY
*ELASTIC
3.0E6,.20
*DENSITY
 2.26e-4
**
** Contact interface.
*INTERFACE,ELSET=INTER
*FRICTION
 1.0,
*SURFACE CONTACT, SOFTENED
0.001, 1.
**A  B
**A: Clearance at which contact pressure is zero.
**B: Pressure at zero clearance (see 6.5.17 of manual)
**
**
*AMPLITUDE,TIME=TOTAL TIME,VALUE=RELATIVE,INPUT=QUAKE.AMP,NAME=EQ
*RESTART,WRITE,FREQUENCY=2
********************************** Loads
** Apply Gravity load in four increments.
*STEP,NLGEOM,INC=100
*STATIC
0.25, 1.
********************************** Boundary Conditions
**
*BOUNDARY 
** Restrain base interface element 
INTER0A,ENCASTRE 
** Restrain X displacements at two top nodes, for roller support 
84,1 
81,1 
*DLOAD
** Block self weight.  g = 386.09 in/s^2
 BLOCK, GRAV, 386.09, 0.0, 0.0, -1.0
** A 1 psi surcharge on the top face of the top block. 0. 1K total.
 10,P2,1.
********************************** OUTPUT
**
** Print dispalcements for all nodes.
*NODE PRINT
 U
**
** Print support reactions and their total
*NODE PRINT,FREQUENCY=1,NSET=INTER0A,TOTAL=YES
 RF
**
** Supress all element printout.
*EL PRINT, FREQUENCY=0
*END STEP
**
********************************** Loads
** Dynamic Analysis for ground motion.
*STEP,NLGEOM,INC=1000
*DYNAMIC,ALPHA=-0.1
.01 , 5.
**A   B
**A: Suggested time increment.
**B: Time period
*BOUNDARY,OP=NEW,TYPE=ACCELERATION,AMPLITUDE=EQ
INTER0B,  1,1,  386.09
INTER0A,  1,1,  386.09
84,  1,1,  386.09
81,  1,1,  386.09
*BOUNDARY,OP=NEW
** Restrain base interface element in all DOFs except X translation 
INTER0A, 2,3
********************************** OUTPUT
**
** Print dispalcements for all nodes.
*NODE PRINT
 U
*NODE PRINT,FREQUENCY=1,NSET=INTER0A,TOTAL=YES
 RF
** Supress element printout for block elements.
*EL PRINT, ELSET=BLOCK, FREQUENCY=0
** Print output for interface elments every second increment.
*EL PRINT, ELSET=INTER, FREQUENCY=2
*END STEP