*HEADING
   SHELL BENDING PROBLEM
**
** Analytical result: vertical steady-state displacement of 240mm 
** at the tip of the beam. This compares with the second beam.
**
** Nodes at corners of the two beams
**
*NODE,NSET=NALL
   101,0.,0.
   141,.4,0.
   201,0.,0.,.02
   241,.4,0.,.02
   1101,0.,0.3
   1141,.4,0.3
   1201,0.,0.3,.02
   1241,.4,0.3,.02
**
** Lower beam
**
*NGEN,NSET=NALL
   101,141
   201,241
**
** Upper beam
**
*NGEN,NSET=NALL
   1101,1141
   1201,1241
**
** Output
**
*NSET,NSET=NOUT
 101,1101
** Elements for both beams.
**
*ELEMENT,TYPE=S4R,ELSET=SHELL
   101,101,201,202,102
   1101,1101,1201,1202,1102
*ELGEN,ELSET=SHELL
   101,40,1,1
   1101,40,1,1
*SHELL SECTION,ELSET=SHELL,MATERIAL=SAMPLE,SECTION INTEGRATION=GAUSS
   .02,
**
** Material properties are purely fictional.
**
*MATERIAL,NAME=SAMPLE
*ELASTIC
   1.E9,0.,
*DENSITY
   10000.,
**
** Right side boundary conditions.
**
*BOUNDARY
   141,ENCASTRE
   241,ENCASTRE
   1141,ENCASTRE
   1241,ENCASTRE
   NALL,3,5
*SURFACE, NAME=TIP1
 101, E1
*SURFACE, NAME=TIP2
 1101, E1
**
** Ramp function for loads
**
*AMPLITUDE,NAME=FUNC
   0.,0.,.2,1.
*RESTART,WRITE,NUMBER INTERVAL=2,TIMEMARKS=NO
**
**
*NSET,NSET=QA_TEST_NOUT
NOUT,
**
*STEP
*DYNAMIC,EXPLICIT
   ,.8
*DSLOAD, CONSTANT RESULTANT=YES, FOLLOWER=NO, AMPLITUDE=FUNC
 TIP1,EDLD,12500.,0.,-1.,0.
*DSLOAD, CONSTANT RESULTANT=YES, FOLLOWER=YES, AMPLITUDE=FUNC
 TIP2,EDLD,12500.,0.,-1.,0.
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT,NSET=QA_TEST_NOUT
 U,
***ELEMENT OUTPUT
** SF,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=NOUT
 U2,
*END STEP