*HEADING
: DISTRIBUTING COUPLING WITH NONLINEAR BEHAVIOR  [DCOUP3D]
*RESTART,W
*NODE
10,0.0,  0.0
1 ,2.0,  0.5
2 ,2.0, -0.5
3 ,1.0,  1.0
*ELEMENT,TYPE=DCOUP3D,ELSET=CONNECTSET
200,10
*DISTRIBUTING COUPLING,ELSET=CONNECTSET,MASS=10.0
2,2.0
3,3.0
1,1.0
*ELEMENT,TYPE=SPRINGA,ELSET=SPRINGS
401,1,2
402,2,3
404,3,1
*SPRING,ELSET=SPRINGS

1.0E8, 
*ELEMENT,TYPE=DASHPOT1,ELSET=EDASHPOT1
110,1
120,2
130,3
*ELEMENT,TYPE=DASHPOT1,ELSET=EDASHPOT2
210,1
220,2
230,3
*ELEMENT,TYPE=DASHPOT1,ELSET=EDASHPOT3
310,1
320,2
330,3
*DASHPOT,ELSET=EDASHPOT1
1, 
01.0E2, 
*DASHPOT,ELSET=EDASHPOT2
2, 
02.0E2, 
*DASHPOT,ELSET=EDASHPOT3
3, 
03.0E2, 
*BOUNDARY
10,2,5
*AMPLITUDE,NAME=SINE,DEFINITION=PERIODIC
1,6.28,0.0,0.0
0.0,1.0
**
**  STEP 1: ROTATION IN THE X-Y PLANE
**
*STEP,NLGEOM,INC=100
*STATIC
0.0125,1.0,0.0125,0.0125
*BOUNDARY
10,6,6,6.283185308
*BOUNDARY,AMPLITUDE=SINE
10,1,1,1.0
*EL FILE,ELSET=CONNECTSET,FREQ=10
NFORC,
*NODE FILE,FREQ=10
U,RF
*END STEP
**
**  STEP 2: ROTATION IN THE Z-X PLANE
**
*STEP,NLGEOM,INC=100
*STATIC
0.0125,1.0,0.0125,0.0125
*BOUNDARY,TYPE=VELOCITY,OP=NEW
10,5,5,6.283185308
*BOUNDARY,AMPLITUDE=SINE,OP=NEW
10,3,3,1.0
10,1,2,0.0
10,4,4,0.0
10,6,6,0.0
*EL FILE,ELSET=CONNECTSET,FREQ=10
NFORC,
*NODE FILE,FREQ=10
U,RF
*END STEP
**
**  STEP 3: ROTATION IN THE Y-Z PLANE
**
*STEP,NLGEOM,INC=100
*STATIC
0.0125,1.0,0.0125,0.0125
*BOUNDARY,TYPE=VELOCITY,OP=NEW
10,4,4,6.283185308
*BOUNDARY,AMPLITUDE=SINE,OP=NEW
10,2,2,1.0
10,1,1,0.0
10,3,3,0.0
10,5,6,0.0
*EL FILE,ELSET=CONNECTSET,FREQ=10
NFORC,
*NODE FILE,FREQ=10
U,RF
*END STEP
**
**  STEP 4: DIRECT INTEGRATION DYNAMIC
**
*STEP,INC=200,NLGEOM
*DYNAMIC,DIRECT
0.09375,1.0,
*BOUNDARY,OP=NEW,TYPE=VELOCITY
10,6,6,6.283185308
*BOUNDARY,OP=NEW,AMPLITUDE=SINE
10,1,1,1.0
10,2,5,0.0
*EL FILE,ELSET=CONNECTSET,FREQ=4
NFORC,
*NODE FILE,FREQ=4
U, 
V, 
A, 
*EL PRINT, ELSET=CONNECTSET, SUMMARY=NO,FREQ=5
 ELKE,
*END STEP