*HEADING
   SHELL BENDING PROBLEM
**
** Analytical result: vertical static 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
**
*STEP, NLGEOM=YES
*STATIC
 0.1,1.,0.1,0.1
*DSLOAD, CONSTANT RESULTANT=YES, FOLLOWER=NO
 TIP1,EDLD,12500.,0.,-1.,0.
*DSLOAD, CONSTANT RESULTANT=YES, FOLLOWER=YES
 TIP2,EDLD,12500.,0.,-1.,0.
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=NOUT
 U2,
*EL PRINT, FREQ=0
*NODE PRINT, FREQ=999, NSET=NOUT
 U,
*END STEP