*HEADING
   CLOAD test problem.  
**  
**  There are two beams made from 6-node, 2-D solids.  Each one has 25
**  elements along the length and 2 elements through the thickness.
**  The nodes on the 
**  right boundary are pinned.  The nodes on the left boundary are 
**  constrained to remain in a straight line without sliding - the 
**  LINEAR MPC.
**
**  Both beams are loaded with a CLOAD couple in the X-direction so 
**  that a pure moment is applied to the left end of the beam.  The 
**  CLOADs on the top beam are follower forces, while the bottom are 
**  constant in direction. All loads are ramped linearly through time 
**  and have the same magnitude.  On the right boundary the nodes are
**  pinned.  The material is elastic.
**
** Nodes at corners of the two beams
**
*NODE
   101,0.,0.
   151,.4,0.
   501,0.,.02
   551,.4,.02
   1101,0.,0.3
   1151,.4,0.3
   1501,0.,.32
   1551,.4,.32
**
** Bottom of lower beam
**
*NGEN,NSET=B1
   101,151
**
** Top of lower beam
**
*NGEN,NSET=T1
   501,551
**
** Bottom of upper beam
**
*NGEN,NSET=B2
   1101,1151
**
** Top of upper beam
**
*NGEN,NSET=T2
   1501,1551
**
** Fill both beams
**
*NFILL
   B1,T1,4,100
   B2,T2,4,100
**
** Constrain left end of beams.
**
*MPC
   BEAM,101,501
   LINEAR,201,101,501
   LINEAR,301,101,501
   LINEAR,401,101,501
   BEAM,1101,1501
   LINEAR,1201,1101,1501
   LINEAR,1301,1101,1501
   LINEAR,1401,1101,1501
**
** Elements for both beams.
**
*ELEMENT,TYPE=CPE6M,ELSET=SOLID
   101,101,103,303,102,203,202
   102,303,301,101,302,201,202	
   1101,1101,1103,1303,1102,1203,1202
   1102,1303,1301,1101,1302,1201,1202
*ELGEN,ELSET=SOLID1
   101,25,2,2,2,200,100
   102,25,2,2,2,200,100   
*ELGEN,ELSET=SOLID2
   1101,25,2,2,2,200,100
   1102,25,2,2,2,200,100
*ELSET,ELSET=EOUTPUT
110, 111
*ELSET,ELSET=SOLID
SOLID1,SOLID2
*SECTION CONTROLS, NAME=RATE
*SOLID SECTION,ELSET=SOLID,MATERIAL=SAMPLE, CONTROLS=RATE
**
** Material properties are purely fictional.
**
*MATERIAL,NAME=SAMPLE
*ELASTIC
   1.E9,0.,
*DENSITY
   10000., 
***
*ELEMENT, TYPE=ROTARYI,ELSET=ROTARY
 50000,101
 60000,501
 70000,1101
 80000,1501
*ROTARY INERTIA,ELSET=ROTARY
1.E-10,1.E-10,1.E-10
**
** Right side boundary conditions.
**
*BOUNDARY
   151,PINNED
   251,PINNED
   351,PINNED
   451,PINNED
   551,PINNED
   1151,PINNED
   1251,PINNED
   1351,PINNED
   1451,PINNED
   1551,PINNED
**
** Ramp function for loads
**
*AMPLITUDE,NAME=FUNC
   0.,0.,.2,1.
*RESTART,TIMEMARKS=YES,WRITE,NUMBER INTERVAL=1
**
*NSET,NSET=NOUT
101,1101
**
*STEP
*DYNAMIC,EXPLICIT
   ,.2
**
** Constant direction couple on lower beam.
**
*CLOAD,AMPLITUDE=FUNC
   101,1,5.1E5
   501,1,-5.1E5
**
** Follower couple on upper beam.
**
*CLOAD,FOLLOWER,AMPLITUDE=FUNC
   1101,1,5.1E5
   1501,1,-5.1E5
**
*MONITOR,NODE=1101,DOF=6
*FILE OUTPUT,NUMBER=4,TIMEMARKS=YES
*NODE FILE, NSET=NOUT
U,
*NSET, NSET=QA_TEST
NOUT
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT,HISTORY,TIME INTERVAL=0.2
*ENERGY OUTPUT, VAR=ALL
***OUTPUT,HISTORY,TIME INTERVAL=0.0002
***NODE OUTPUT,NSET=NOUT
**UR3,
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLWK,ETOTAL
*END STEP