*HEADING
: PRESSURE LOADING Of A QUARTER CYLINDER; SINGLE LAYER OF GASKET
          ELEMENTS CONNECTED TO C3D27R ELEMENTS USING SHARED NODES.  [GK3D18]
**
** Add stick out element and nodes - extension of element 139
**-----------------------------------------------
*NODE,NSET=EXTENSION1
*RESTART,WRITE,FREQ=999
*NODE
*NSET,NSET=EXTENSION
EXTENSION1,439,440,441,502,503,504
**The above set corresponds to the nodes that make up the overhanging
**part of the gasket (on one side).
*NCOPY,CHANGE NUMBER=500,OLD SET=EXTENSION,SHIFT,NEW SET=EXTENSION2
**-----------------------------------------------
*ELEMENT, TYPE=C3D27R, ELSET=SOLID
**
*ELEMENT, TYPE=GK3D18, ELSET=GASKET
**
*ELEMENT, TYPE=GK3D18, ELSET=OVERHANG1
*ELEMENT, TYPE=GK3D18, ELSET=OVERHANG2
**
**-----------------------------------------------
*SOLID SECTION, ELSET=SOLID, MATERIAL=STEEL,ORIENTATION=CYL
*ORIENTATION,SYSTEM=C,NAME=CYL
*ORIENTATION,SYSTEM=C,NAME=CYL1
*MATERIAL, NAME=STEEL
*ELASTIC, TYPE=ISO
2.E5,0.3
**-----
*GASKET SECTION,ELSET=GASKET,BEHAVIOR=IS,ORIENTATION=CYL1
*GASKET BEHAVIOR,NAME=IS
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING
**-----
*GASKET SECTION,ELSET=OVERHANG1,BEHAVIOR=IS2,ORIENTATION=CYL1
*GASKET BEHAVIOR,NAME=IS2
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING
*GASKET ELASTICITY,COMPONENT=MEMBRANE
**-----
*GASKET SECTION,ELSET=OVERHANG2,BEHAVIOR=IS3,ORIENTATION=CYL1,
STABILIZATION STIFFNESS=0.0
*GASKET BEHAVIOR,NAME=IS3
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING
**
**-----------------------------------------------
*NSET,NSET=YSYMM0,GEN
*NSET,NSET=YSYMM0_T
*NSET,NSET=YSYMM0_B
*NSET,NSET=YSYMM0_M
*NSET,NSET=YSYMM1,GEN
*NSET,NSET=YSYMM1_T
*NSET,NSET=YSYMM1_B
*NSET,NSET=YSYMM1_M
*NSET,NSET=TOP,GEN
*NSET,NSET=BOTTOM,GEN
*NSET,NSET=IN_CORN,GEN
*NSET,NSET=IN_MID,GEN
*NSET,NSET=IN_MID_T
*NSET,NSET=IN_MID_B
*NSET,NSET=IN_MID_M
*TRANSFORM,TYPE=R,NSET=IN_MID_T
*TRANSFORM,TYPE=R,NSET=IN_MID_B
*TRANSFORM,TYPE=R,NSET=IN_MID_M
*TRANSFORM,TYPE=R,NSET=YSYMM1_T
*TRANSFORM,TYPE=R,NSET=YSYMM1_B
*TRANSFORM,TYPE=R,NSET=YSYMM1_M
*BOUNDARY
IN_CORN,1
IN_MID,1
BOTTOM,3
TOP,3
YSYMM0,2,3
YSYMM1,2,3
**
*ELSET,ELSET=EOUT
**-----------------------------------------------
*STEP
*STATIC
**
** apply outside load
**
*DLOAD, OP=NEW
109,p2,100.
*NODE PRINT
U,RF
*EL PRINT,FREQ=0
S,
E,
*EL PRINT,FREQ=999,ELSET=EOUT
E,
S,
NFORC,
CS11,
*NODE FILE,FREQ=1
U,RF
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT
U,RF
*EL FILE,FREQ=1
S,E
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
S,E
*END STEP