*HEADING
 block-2.inp: An unreinfoced masonr column modelled as a stack of discrete
 10x10x10 in cube elements, connected by interface elements.  There are
 horizontal concentrated loads at mid-height.
 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,
**
********************************** Boundary Conditions
**
*BOUNDARY 
** Restrain base interface element 
INTER0A,ENCASTRE 
** Restrain X displacements at two top nodes, for roller support 
84,1 
81,1 
**
*RESTART,WRITE,FREQUENCY=1
********************************** Loads
** Apply Gravity load in four increments.
*STEP,NLGEOM,INC=100
*STATIC
0.25, 1.
** g = 384 in/s^2
*DLOAD
 BLOCK, GRAV, 384.0, 0.0, 0.0, -1.0
********************************** OUTPUT
**
** Print dispalcements for all nodes.
*NODE PRINT
 U
**
** Print support reactions and their total
*NODE PRINT,FREQUENCY=1,NSET=INTER0A,TOTAL=YES
 RF
*END STEP
**
**
*STEP,NLGEOM,INC=100
** Apply lateral load in four increments.
*STATIC,RIKS
 .25,1.,1.0e-5,1.0e10,1.0,45,1,20.0
** A B      C   D      E   F   G  H
** A: Suggested time increment.
** B: Time interval
** C: Minimum time increment.
** D: Maximum time increment.
** E: Maximum load proportionality factor.
** F: Number of monitored node.
** G: DOF number for monitored node.
** H: Maximum value for displacement at monitored node.
********************************** OUTPUT
**
** Print dispalcements for all nodes.
*NODE PRINT
 U
**
** Print support reactions and their total
*NODE PRINT,FREQUENCY=1,NSET=INTER0A,TOTAL=YES
 RF
**
** Load two of the mid-span nodes with horizontal loads of 100 pounds each.
*CLOAD
 45,1,100.
 48,1,100.
**A B C
** A: Node number
** B: DOF number
** C: Load magnitude
**
*END STEP