*HEADING
 PURE BENDING OF CURVED PIPE WITH SHELL ELEMENTS     
 UNITS: LB, IN
*NODE
 101, 0.,-7.807,0.
 113, 0., 7.807,0.
 201, .3373297637,-7.810513979,0.
 213, .6625978999,7.800097689,0.
 301,.6745131222,-7.82105439,0.
 313, 1.324908224,7.779393751,0.
*NODE,NSET=N300
 300,.9997106733,-.0208303195,0.
 3000,.9997106733,-.0208303195,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,.4999638318,-.005208145,0.,-.999782994,.0208318263,0.
 301,313,1,0,.9997106733,-.0208303195,0.,-.99913207,.0416546114,0.
*ELEMENT,TYPE=S8R5,ELSET=ALL
 1, 101,301,303,103,201,302,203,102
*ELGEN,ELSET=ALL
 1,6,2
*SHELL SECTION,MATERIAL=A1,ELSET=ALL
 .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
**
*orientation,name=curv,def=coord,sys=rect
0.99913207,-0.0416546114,0.,0.,1.,0.
*SURFACE,NAME=SKCOUP1,TYPE=NODE
301,
313,
*COUPLING,CONSTRAINT NAME=KC1,REF NODE=300,SURFACE=SKCOUP1,orientation=curv
*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,orientation=curv
*KINEMATIC
1,
5,6
**
*NSET,NSET=PXY
 101,201,301,113,213,313,300
*BOUNDARY
 PXY,3,5
 SYMM,1
 SYMM,5,6
 300,2,3
 300,4,5
*RESTART,WRITE,FREQUENCY=2
**
** uncomment the following step
** to add pressurization
**
***STEP,INC=1,NLGEOM
**  PRESSURIZE PIPE
***STATIC
**1.0, 1.0
***DLOAD,OP=NEW
** ALL,P,500.
***CLOAD,OP=NEW,FOLLOWER
** 301,1,1328.54
** 301,2,-55.38
** 302,1,5314.2
** 302,2,-221.56
** 303,1,2657.12
** 303,2,-110.78
** 304,1,5314.2
** 304,2,-221.56
** 305,1,2657.12
** 305,2,-110.78
** 306,1,5314.2
** 306,2,-221.56
** 307,1,2657.12
** 307,2,-110.78
** 308,1,5314.2
** 308,2,-221.56
** 309,1,2657.12
** 309,2,-110.78
** 310,1,5314.2
** 310,2,-221.56
** 311,1,2657.12
** 311,2,-110.78
** 312,1,5314.2
** 312,2,-221.56
** 313,1,1328.54
** 313,2,-55.38
***EL PRINT
**LOADS, 
**S, 
**SINV, 
**E, 
***END STEP
*STEP,NLGEOM,INC=100
*STATIC
 .1,1
*BOUNDARY
300,6,6,-.010
*OUTPUT,FIELD,FREQUENCY=2
*ELEMENT OUTPUT
S,E,PE,PEEQ,PEMAG
*NODE OUTPUT
U,
*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