*HEADING
: CONTACT PAIR WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL,
          HIGH INITIAL STIFFNESS, NORMAL ONLY BEHAVIOR [GK3D2N]
*NODE,NSET=NALL
*ELEMENT,TYPE=GK3D2N,ELSET=EINT
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
*GASKET BEHAVIOR,NAME=IS
*GASKET THICKNESS BEHAVIOR,DIRECTION=LOADING,VARIABLE=STRESS
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING,VARIABLE=STRESS
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
0.0,0.0, 1.0,-30., 2.0,-35., 3.0,-0.1,
4.0,-35., 5.0,-10., 6.0,-0.1, 7.0,-37.,
8.0,-40., 9.0,-0.0, 10.0,-40., 11.0,-42.,
12.0,-0.1, 13.0,-42., 14.0,-43., 15.0,-0.1,
16.0,-43., 17.0,-45, 18.0,-0.1, 19.0,-45,
20.0,-60, 21.0,-0.1, 22.0,-60, 23.0,-75.,
24.0,-0.1, 25.0,-75., 26.0,-83., 27.0,-0.1,
28.0,-83., 29.0,-90., 30.0,-0.1
*BOUNDARY
1,1,1
*******************************************
*NODE
*ELEMENT,TYPE=R3D4,ELSET=RBODY
*RIGID BODY,ELSET=RBODY,REFNODE=10000
*surface,type=node,NAME=CNSET
*surface,NAME=RSURF
RBODY,SNEG
*CONTACT PAIR,INTERACTION=INT1,TIED,ADJUST=0.01
CNSET,RSURF
*SURFACE INTERACTION,NAME=INT1
1.,
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,30.,0.0001,0.01
**CONTROLS,PARAMETER=TIME INCREMENT
**
**,0.25
*BOUNDARY
10000,1,1
10000,3,6
*CLOAD,AMP=HIST
10000,2,1.
*NODE PRINT,FREQ=0
U,
RF,
*NODE FILE,FREQ=148
U,
RF,
*EL PRINT,FREQ=0,ELSET=EINT
S,
E,
NE,
EE,
PE,
EVOL,
THE,
TEMP,
*EL FILE,FREQ=148,ELSET=EINT
E,
EE,
NE,
PE,
S,
EVOL,
THE,
*OUTPUT,FIELD,FREQ=148,VAR=PRE
*CONTACT PRINT,MASTER=RSURF,SLAVE=CNSET,FREQ=0
*CONTACT FILE,MASTER=RSURF,SLAVE=CNSET,FREQ=148
*END STEP
*******************************************
**No Test of Membrane Behavior of Element
*******************************************
**No Test Transverse Shear Stiffness of Element
*******************************************
*STEP
Thermal Expansion of Element
*STATIC
*CONTROLS,RESET
*CLOAD,OP=NEW
*BOUNDARY,OP=NEW
10000,1,6
1,1
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************