*HEADING
 BRAZIER PROBLEM WITH 6 SHELL ELEMENTS          
 KINEMATIC COUPLING CONSTRAINT AND DISTRIBUTING COUPLING ELEMENT
*NODE
 101, 0.,-7.807,0.
 113, 0.,7.807,0.
 201, .5,-7.807,0.
 213, .5,7.807,0.
 301, 1.,-7.807,0.
 313, 1.,7.807,0.
*NODE,NSET=N300
300,1.,0.,0.
*NODE
3000,1.,0.,0.
*NGEN,LINE=C,NSET=SYMM
 101,113,1,0,0.,0.,0.,-1.
*NSET,NSET=NSYMMOUT
 101,113
*NGEN,LINE=C
 201,213,2,0,.5,0.,0.,-1.
 301,313,1,0,1.,0.,0.,-1.
*ELEMENT,TYPE=S8R5
1,101,301,303,103,201,302,203,102
*ELGEN,ELSET=EALL
1,6,2
*ELSET,ELSET=PRIN
 1,6
*SHELL SECTION,MATERIAL=A1,ELSET=EALL
.41,7
*MATERIAL,NAME=A1
*ELASTIC
28.1E6,.2642
**
** uncomment the following option for the 
** elastic-plastic case
**
***PLASTIC
**39440.,0.
**50170.,.00473
**54950.,.01264
**58540.,.02836
**61520.,.0491
**76520.,.105
**
** REMOVE RIGID-BODY MOTION IN Y DIRECTION
** 
*BOUNDARY
3000,2
*SURFACE,NAME=SDCOUP,TYPE=NODE
301,1
302,4
303,2
304,4
305,2
306,4
307,2
308,4
309,2
310,4
311,2
312,4
313,1
*COUPLING,CONSTRAINT NAME=DCOUP,REF NODE=3000,SURFACE=SDCOUP
*DISTRIBUTING
**
** CONSTRAIN THE SECTION TO REMAIN IN THE A3 PLANE
**
*SURFACE,NAME=SKCOUP1,TYPE=NODE
301,
313,
*COUPLING,CONSTRAINT NAME=KC1,REF NODE=300,SURFACE=SKCOUP1
*KINEMATIC
1,
3,6 
*SURFACE,NAME=SKCOUP2,TYPE=NODE
302,
303,
304,
305,
306,
307,
308,
309,
310,
311,
312,
*COUPLING,CONSTRAINT NAME=KC2,REF NODE=300,SURFACE=SKCOUP2
*KINEMATIC
1,
5,6
**
*NSET,NSET=PXY
101,201,301,113,213,313
*BOUNDARY
PXY,3,5
SYMM,1
SYMM,5,6
300,2,3
300,4,5
*NSET,NSET=N300
300, 
3000, 
*RESTART,WRITE,FREQUENCY=1
**
** uncomment the following step
** to add pressurization
**
***STEP,INC=1,NLGEOM
**  PRESSURIZE PIPE
***STATIC
**1.0, 1.0
***DLOAD,OP=NEW
** EALL,P,500.
***CLOAD,OP=NEW,FOLLOWER
** 301,1,1329.71
** 302,1,5318.82
** 303,1,2659.41
** 304,1,5318.82
** 305,1,2659.41
** 306,1,5318.82
** 307,1,2659.41
** 308,1,5318.82
** 309,1,2659.41
** 310,1,5318.82
** 311,1,2659.41
** 312,1,5318.82
** 313,1,1329.71
***EL PRINT
**LOADS, 
**S, 
**SINV, 
**E, 
***END STEP
*STEP,INC=1000,NLGEOM
*STATIC
 .1,1.0
*BOUNDARY
 300,6,6,.01
*OUTPUT,FIELD,FREQUENCY=2
*ELEMENT OUTPUT
S,E,PE,PEEQ,PEMAG
*NODE OUTPUT
U,RF
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=N300
U,RF
*NODE OUTPUT,NSET=NSYMMOUT
U,
*EL FILE, FREQUENCY=100
S,E, 
*NODE FILE,NSET=N300
U, 
RF, 
*MONITOR,NODE=313,DOF=2
*END STEP