*HEADING
   CLOAD test problem.  File: 
**  
**  There are two beams made from 4-node, 2-D solids.  Each one has 20
**  elements along the length and 4 elements through the thickness for
**  a total of 160 CPE4R elements and 210 nodes.  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.
   121,.4,0.
   501,0.,.02
   521,.4,.02
   1101,0.,0.3
   1121,.4,0.3
   1501,0.,.32
   1521,.4,.32
**
** Bottom of lower beam
**
*NGEN,NSET=B1
   101,121
**
** Top of lower beam
**
*NGEN,NSET=T1
   501,521
**
** Bottom of upper beam
**
*NGEN,NSET=B2
   1101,1121
**
** Top of upper beam
**
*NGEN,NSET=T2
   1501,1521
**
** 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=CPE4R,ELSET=SOLID
   101,101,102,202,201
   1101,1101,1102,1202,1201
*ELGEN,ELSET=SOLID1
   101,20,1,1,4,100,100
*ELGEN,ELSET=SOLID2
   1101,20,1,1,4,100,100
*ELSET,ELSET=SOLID
SOLID1,SOLID2
*SOLID SECTION,ELSET=SOLID,MATERIAL=SAMPLE
**
** 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
   121,PINNED
   221,PINNED
   321,PINNED
   421,PINNED
   521,PINNED
   1121,PINNED
   1221,PINNED
   1321,PINNED
   1421,PINNED
   1521,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,4.90E5
   501,1,-4.90E5
**
** Follower couple on upper beam.
**
*CLOAD,FOLLOWER,AMPLITUDE=FUNC
   1101,1,4.90E5
   1501,1,-4.90E5
**
*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,ELSET=SOLID1
**ALLSE,ALLKE
***ENERGY OUTPUT,ELSET=SOLID2
**ALLSE,ALLKE
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLWK,ETOTAL
**DT,
*PRINT, ETOTAL=YES
*END STEP