*HEADING
Surface-based Pressure Penetration Analysis of
Air Duct Kiss Seal
*NODE,NSET=NSEAL,
INPUT=presspenairductseal_node.inp
*ELEMENT,TYPE=CAX4H,ELSET=STEEL,
INPUT=presspenairductseal_elem_metal.inp
*ELEMENT,TYPE=CAX4H,ELSET=RUBFIB,
INPUT=presspenairductseal_elem_rub.inp
*ELSET,ELSET=FACE1,GENERATE
437,439,1
443,452,1
*ELSET,ELSET=FACE3,GENERATE
28,36,1
67,76,1
227,276,1
307,316,1
440,442,1
453,462,1
473,475,1
*ELSET,ELSET=FACE4,GENERATE
437,440,3
*ELSET,ELSET=PPRES
476,477,478,479,480,481,482,503,504,
505,506,507,508,509,510,511,512,513,
514,515,516,517,518,519,520,521,522,
533,534,535,536
*NSET,NSET=EDGE
1,2,3,4,5
*NODE,NSET=BCYL
5001,20.5,2.9,0.
5002,20.5,0.,0.
*RIGID BODY,ANALYTICAL SURFACE=CFACE,REFNODE=5001
*SURFACE,NAME=CFACE,TYPE=SEGMENTS
START,21.5,2.9,
LINE,19.5,2.9
*SURFACE,NAME=PPRES
PPRES,S3
*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=5002
*SURFACE,NAME=CFACE2,TYPE=SEGMENTS
START,19.5,0.,
LINE,21.5,0.
*ELSET,ELSET=ELBOT,GENERATE
2,9,1
37,40,1
*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
568,582,1
*ELSET,ELSET=SL1,GENERATE
77,122,1
*SURFACE,NAME=SECONDARY1
MA1, 
*SURFACE,NAME=MAIN1
SL1, 
*CONTACT PAIR, INTERACTION=SELF1
SECONDARY1,MAIN1
*SURFACE INTERACTION, NAME=SELF1
*FRICTION,SLIPTOLERANCE=.2
.1, 
*SURFACE BEHAVIOR, 
PRESSURE-OVERCLOSURE=EXPONENTIAL
.005,1000.
*ELSET,ELSET=MA2,GENERATE
307,315,1
*ELSET,ELSET=SL2,GENERATE
443,451,1
*SURFACE,NAME=MAIN2
MA2, 
*SURFACE,NAME=SECONDARY2
SL2, 
*CONTACT PAIR, INTERACTION=SELF2
SECONDARY2,MAIN2
*SURFACE INTERACTION, NAME=SELF2
*FRICTION,SLIPTOLERANCE=.2
.1, 
*SURFACE BEHAVIOR, 
PRESSURE-OVERCLOSURE=EXPONENTIAL
.0002,1000.
*ELSET,ELSET=MA3,GENERATE
28,36,1
*ELSET,ELSET=SL3,GENERATE
227,270,1
*SURFACE,NAME=MAIN3
MA3, 
*SURFACE,NAME=SECONDARY3
SL3, 
*CONTACT PAIR, INTERACTION=SELF3
SECONDARY3,MAIN3
*SURFACE INTERACTION, NAME=SELF3
*FRICTION
.1, 
*SURFACE BEHAVIOR, 
PRESSURE-OVERCLOSURE=EXPONENTIAL
.0002,1000.
*ELSET,ELSET=MA4,GENERATE
612,622,1
*ELSET,ELSET=SL4,GENERATE
550,562,1
*SURFACE,NAME=MAIN4,TRIM=YES
MA4, 
*SURFACE,NAME=SECONDARY4
SL4, 
*CONTACT PAIR, INTERACTION=SELF4
SECONDARY4,MAIN4
*SURFACE INTERACTION, NAME=SELF4
*FRICTION,SLIPTOLERANCE=.2
.1, 
*SURFACE BEHAVIOR, 
PRESSURE-OVERCLOSURE=EXPONENTIAL
.001,1000.
*ELSET,ELSET=MA5,GENERATE
603,610,1
*ELSET,ELSET=SL5,GENERATE
116,124,1
*SURFACE,NAME=MAIN5
MA5, 
*SURFACE,NAME=SECONDARY5
SL5, 
*CONTACT PAIR, INTERACTION=SELF5
SECONDARY5,MAIN5
*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
-211.,-.4
-172.,-.35
-145.,-.30
-122.,-.25
-106.,-.20
-89.,-.15
-66.,-.1
-36.,-.05
0.,0.
132.,.05
231.,.1
297.,.15
363.,.2
594.,.25
924.,.3
1188.,.4
*SOLID SECTION,ELSET=STEEL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
3.0e+7,.3
*BOUNDARY
EDGE,1,1
1,2,2
*RESTART,WRITE,FREQUENCY=2
*STEP,UNSYMM=YES,NLGEOM
*STATIC
.0125,1.,1.E-8,.05
*BOUNDARY,OP=MOD
5001,1,,0.
5001,2,,-1.4
5001,6,,0.
5002,1,,0.
5002,2,,0.
5002,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,SECONDARY=PPRES,MAIN=CFACE,FREQ=100
CSTRESS,
*ENERGY PRINT
*CONTACT FILE,SECONDARY=PPRES,MAIN=CFACE,FREQ=100
CSTRESS,
*OUTPUT,FIELD,FREQ=100
*CONTACT OUTPUT,SECONDARY=PPRES,MAIN=CFACE
CSTRESS,
*NODE output
U,
*END STEP
*STEP,UNSYMM=YES,NLGEOM,inc=100
*STATIC,stabilize,ALLSDTOL=0.06
.0125,1.,1.E-8,.05
*DLOAD,OP=NEW
FACE1,P1,30.
FACE3,P3,30.
FACE4,P4,30.
*PRESSURE PENETRATION,MAIN=CFACE,SECONDARY=PPRES
597,,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,SECONDARY=PPRES,MAIN=CFACE,FREQ=2
CSTRESS,PPRESS
*ENERGY PRINT
*CONTACT FILE,SECONDARY=PPRES,MAIN=CFACE,FREQ=2
CSTRESS,PPRESS
*OUTPUT,FIELD,FREQ=1
*CONTACT OUTPUT,SECONDARY=PPRES,MAIN=CFACE
PPRESS
*element output
else
elsd
*NODE output
U, 
*Output, history, variable=PRESELECT,freq=1
*END STEP