*HEADING
 LARGE STRAIN VISCO - RIGID BODY MOTION - CPS4
 APPLY LOAD INSTANTANEOUSLY AND LET RELAX .DO A STATIC ROTATION
 TO MAKE SURE STATE VARIABLES REMAIN ENFORCED AND ROTATE
*RESTART,WRITE
*NODE
1,
2,10.
3,10.,5.,
4,0.,5.,
99,-10.,0.
*ELEMENT,TYPE=CPS4,ELSET=ONE
1,1,2,3,4
*SOLID SECTION,ELSET=ONE,MATERIAL=BIDERMAN
1., 
*MATERIAL,NAME=BIDERMAN
*HYPERELASTIC,N=1,MODULI=INSTANTANEOUS
27.02,1.42,0.000001
*VISCOELASTIC,TIME=PRONY
0.25,0.25,5.
0.25,0.25,10.
*ELSET,ELSET=EFILE
1, 
*surface,NAME=SLVC
ONE,S2
*surface,NAME=SLVD
ONE,S4
*rigid body,analytical surface=MSTC,REF NODE=99
*surface,name=MSTC,TYPE=SEGMENTS
START,10.,16.
LINE,10.,-11.
*rigid body,analytical surface=MSTD,REF NODE=999
*surface,TYPE=SEGMENTS,name=MSTD
START,0.,-11.
LINE,0.,16.
*CONTACT PAIR,INTERACTION=INTER1
SLVC,MSTC
SLVD,MSTD
*SURFACE INTERACTION,NAME=INTER1
*ELEMENT,TYPE=SPRING1,ELSET=SPRINGX
11,1
*SPRING,ELSET=SPRINGX
1, 
1., 
*BOUNDARY
1,2,2
2,2,2
999,1,6
99,6,6
99,2,2
**NODE 1 WITH FIXED BC IN DIR. 2 AND IN CONTACT WITH MAIN SURFACE MSTD
**CONFLICTS WITH THE NORMAL CONTACT CONSTRAINT IN STEP 3 WHEN THE CONTACT
**SURFACES ARE ROTATED. OVERCONSTRAINT CHECKS GIVE A WARNING AND CHANGE 
**THE CONTACT STATUS TO OPEN AUTOMATICALLY TO AVOID OVERCONSTRAINT. THIS
**CAUSES CONVERGENCE PROBLEM. THE FOLLOWING KEYWORD IS BEING ADDED TO AVOID
**THIS PROBLEM. IN THIS CASE NODE 1 AND REF. NODE 999 FOR MAIN SURFACE
**MSTD ARE COINCIDENT AND THIS OVERCONSTRAINT SHOULD NOT CAUSE ANY PROBLEM.
**OVERCONSTRAINT CHECKS WILL BE IMPROVED LATER ON TO DETECT AND IGNORE
**THESE CASES (WHEN SECONDARY AND MAIN NODES ARE COINCIDENT).
*CONSTRAINT CONTROLS,NO CHECKS
*STEP,NLGEOM,INC=25
*STATIC
.1,1.
*BOUNDARY
99,1,1,2.
*END STEP
*STEP,NLGEOM,INC=25
*VISCO
5.,50.,2.,10.,
*BOUNDARY
99,1,1,10.
*END STEP
*STEP,NLGEOM,INC=200
*STATIC
30.,90.,1.,30.
*BOUNDARY,OP=NEW
1,2,2
999,1,2
999,6,6,1.5707963
99,2,2
99,1,1,10.
99,6,6,1.5707963
*END STEP
*STEP,NLGEOM
*VISCO
10.,100.
*BOUNDARY,OP=NEW
1,1,1
999,1,2
999,6,6,1.5707963
99,2,2
99,1,1,10.
99,6,6,1.5707963
*EL FILE,F=1,ELSET=EFILE
S,E,CE,CEP
*CONTACT FILE,F=1
CSTRESS,
CDISP, 
*NODE FILE,F=1
U,RF
*END STEP