*HEADING
 Surface-based Pressure Penetration Analysis of Air Duct Kiss Seal
 Using C3D8H elements -- Modelling a 5-degree fraction
**PREPRINT,MODEL=YES,HISTORY=YES
*NODE,NSET=NSEAL
*ELEMENT, TYPE=C3D8H,ELSET=RUBFIB
*ELEMENT, TYPE=C3D8H,ELSET=STEEL
**
*ELSET,ELSET=FACE3,GENERATE
*ELSET,ELSET=FACE5,GENERATE
*ELSET,ELSET=FACE6,GENERATE
*ELSET,ELSET=PPRES
*NSET,NSET=EDGE
*NODE,NSET=BCYL
*RIGID BODY,ANALYTICAL SURFACE=CFACE,REFNODE=500001
*SURFACE,NAME=CFACE,TYPE=REVOLUTION
START,21.5,2.9
LINE,19.5,2.9
*SURFACE,NAME=PPRES
PPRES,S5
*CONTACT PAIR,INTERACTION=PPEN
PPRES,CFACE
*SURFACE INTERACTION,NAME=PPEN
*FRICTION
.00001,
*SURFACE BEHAVIOR,
PRESSURE-OVERCLOSURE=EXPONENTIAL
.0002,1000.
*RIGID BODY,ANALYTICAL SURFACE=CFACE2,REF NODE=500002
*SURFACE,NAME=CFACE2,TYPE=REVOLUTION
START,19.5,0.
LINE,21.5,0.
*ELSET,ELSET=ELBOT,GENERATE
*SURFACE,NAME=BOTTOM
ELBOT,
*CONTACT PAIR, INTERACTION=FLANGE
BOTTOM,CFACE2
*SURFACE INTERACTION, NAME=FLANGE
*FRICTION
.1,
*SURFACE BEHAVIOR,
PRESSURE-OVERCLOSURE=EXPONENTIAL
.0005,1000.
*ELSET,ELSET=MA1,GENERATE
*ELSET,ELSET=SL1,GENERATE
*SURFACE,NAME=SLAVE1
MA1,
*SURFACE,NAME=MASTER1
SL1,
*CONTACT PAIR, INTERACTION=SELF1
SLAVE1,MASTER1
*SURFACE INTERACTION, NAME=SELF1
*FRICTION,SLIPTOLERANCE=.2
.1,
*SURFACE BEHAVIOR,
PRESSURE-OVERCLOSURE=EXPONENTIAL
.005,1000.
*ELSET,ELSET=MA2,GENERATE
*ELSET,ELSET=SL2,GENERATE
*SURFACE,NAME=MASTER2
MA2,
*SURFACE,NAME=SLAVE2
SL2,
*CONTACT PAIR, INTERACTION=SELF2
SLAVE2,MASTER2
*SURFACE INTERACTION, NAME=SELF2
*FRICTION,SLIPTOLERANCE=.2
.1,
*SURFACE BEHAVIOR,
PRESSURE-OVERCLOSURE=EXPONENTIAL
.0002,1000.
*ELSET,ELSET=MA3,GENERATE
*ELSET,ELSET=SL3,GENERATE
*SURFACE,NAME=MASTER3
MA3,
*SURFACE,NAME=SLAVE3
SL3,
*CONTACT PAIR, INTERACTION=SELF3
SLAVE3,MASTER3
*SURFACE INTERACTION, NAME=SELF3
*FRICTION
.1,
*SURFACE BEHAVIOR,
PRESSURE-OVERCLOSURE=EXPONENTIAL
.0002,1000.
*ELSET,ELSET=MA4,GENERATE
*ELSET,ELSET=SL4,GENERATE
*SURFACE,NAME=MASTER4,TRIM=YES
MA4,
*SURFACE,NAME=SLAVE4
SL4,
*CONTACT PAIR, INTERACTION=SELF4
SLAVE4,MASTER4
*SURFACE INTERACTION, NAME=SELF4
*FRICTION,SLIPTOLERANCE=.2
.1,
*SURFACE BEHAVIOR,
PRESSURE-OVERCLOSURE=EXPONENTIAL
.001,1000.
*ELSET,ELSET=MA5,GENERATE
*ELSET,ELSET=SL5,GENERATE
*SURFACE,NAME=MASTER5
MA5,
*SURFACE,NAME=SLAVE5
SL5,
*CONTACT PAIR, INTERACTION=SELF5
SLAVE5,MASTER5
*SURFACE INTERACTION, NAME=SELF5
*FRICTION
.1,
*SURFACE BEHAVIOR,
PRESSURE-OVERCLOSURE=EXPONENTIAL
.001,1000.
*ELSET,ELSET=EALL
STEEL,RUBFIB
*SOLID SECTION,ELSET=RUBFIB,MATERIAL=RUBFIB
*MATERIAL,NAME=RUBFIB
*HYPERELASTIC,OGDEN,N=4,TEST DATA INPUT
*UNIAXIAL TEST DATA
*SOLID SECTION,ELSET=STEEL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
3.0e+7,.3
*nset,nset=n1,generate
*nset,nset=n597,generate
**
*nset,nset=bound1,generate
*nset,nset=bound3,generate
**
*BOUNDARY
EDGE,1
EDGE,3
n1,2,2
bound3,3
**
*EQUATION
 bound1,3,0.996194698091746,bound1,1,-0.08715574274765808
**
*STEP,UNSYMM=YES,NLGEOM
*STATIC
.0125,1.,1.E-8,.05
*BOUNDARY,OP=MOD
500001,1,,0.
500001,2,,-1.4
500001,3,6,0.
500002,1,6,0.
*EL PRINT,FREQUENCY=100
S,
*NODE PRINT,FREQUENCY=100
U,
*NODE PRINT,FREQUENCY=1,NSET=BCYL
U1,U2,RF1,RF2
*EL FILE,FREQUENCY=0
*NODE FILE,FREQUENCY=100
U,
*CONTACT PRINT,SLAVE=PPRES,MASTER=CFACE,FREQ=100
CSTRESS,
*CONTACT FILE,SLAVE=PPRES,MASTER=CFACE,FREQ=100
CSTRESS,
*OUTPUT,FIELD,FREQ=100
*CONTACT OUTPUT,SLAVE=PPRES,MASTER=CFACE
CSTRESS,
*NODE output
U,
*END STEP
*STEP,UNSYMM=YES,NLGEOM
*STATIC
.0125,1.,1.E-8,.05
*DLOAD,OP=NEW
FACE3,P3,30.
FACE5,P5,30.
FACE6,P6,30.
*PRESSURE PENETRATION,MASTER=CFACE,SLAVE=PPRES
n597,,30.0
*EL PRINT,FREQUENCY=100
S,
*NODE PRINT,FREQUENCY=100
U,
*NODE PRINT,FREQUENCY=1,NSET=BCYL
U1,U2,RF1,RF2
*EL FILE,FREQUENCY=100
S,
*NODE FILE,FREQUENCY=100
U,
*CONTACT PRINT,SLAVE=PPRES,MASTER=CFACE,FREQ=2
CSTRESS,PPRESS
*CONTACT FILE,SLAVE=PPRES,MASTER=CFACE,FREQ=2
CSTRESS,PPRESS
*OUTPUT,FIELD
*CONTACT OUTPUT,SLAVE=PPRES,MASTER=CFACE
CSTRESS,PPRESS
*NODE output
U,
*OUTPUT,history
*NODE OUTPUT,NSET=BCYL
RF2
*END STEP