*INCLUDE,INPUT = B031_work_r1.inp
*INCLUDE,INPUT = B031_lowerPlate_r1.inp
*INCLUDE,INPUT = B031_upperPlate.inp
*Node, NSET=cNode
*ELSET, ELSET=Eall1
EvolumesC3D20R
EvolumesC3D20R_L
EvolumesC3D20R_U
**
*MATERIAL, NAME=SolidMaterial
*ELASTIC
210000, 0.300
*PLASTIC
2100.,0.
2110.,0.03
*DENSITY
7.8e-9
*SOLID SECTION,ELSET=Eall1,MATERIAL=SolidMaterial
*RIGID BODY,ELSET=EvolumesC3D20R_L,REF NODE=21000,ROT NODE=21001
*RIGID BODY,ELSET=EvolumesC3D20R_U,REF NODE=22000,ROT NODE=22001
*CONTACT PAIR, INTERACTION=INTFemConstraintContact,TYPE=SURFACE TO SURFACE
W_L,P_L
*SURFACE INTERACTION, NAME=INTFemConstraintContact
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=LINEAR
 1000000.0
 *FRICTION
0.3,5000.
*CONTACT PAIR, INTERACTION=INTFemConstraintContact,TYPE=SURFACE TO SURFACE
W_U,P_U
*SURFACE INTERACTION, NAME=INTFemConstraintContact
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=LINEAR
 1000000.0
 *FRICTION
0.01,5000.
*AMPLITUDE,NAME=Graph1
0,0,
0.5,1,
1,0,
*BOUNDARY
NodeBC_Group_1,2,2
NodeBC_Group_2,1,1
**NodeBC_Group_2,2,3
**NodeBC_Group_2,1,1,-60
**NodeBC_Group_3,3,3
21000,1,3
21001,1,3
22000,1,1
22000,2,2
22001,1,3
*BOUNDARY,AMPLITUDE=graph1
22000,3,3,-5
*STEP, NLGEOM
*STATIC, SOLVER=SPOOLES,direct
0.05,1.0
*DLOAD
Eall1,GRAV,9.8,0.0,0.0,-0.1
*NODE FILE, OUTPUT=3d
 U
*EL FILE
 S, E
**
*NSET,NSET=NodeBC_Group_4
**
*NODE PRINT,NSET=NodeBC_Group_4
 u
**
*NODE PRINT, NSET=NodeBC_Group_4, TOTALS=ONLY
 RF
*END STEP