*HEADING
BOLT DOWN AND OPERATIONAL CYCLE OF COOLANT MANIFOLD COVER    
Units: N-mm-sec-deg C
**
*RESTART,WRITE,FREQ=999
**
** read nodes, elements, and surfaces
** from external file
** 4639 nodes in nset NALL
** 1304 C3D8I elements
**  208 C3D6 elements
**  973 GK3D8N elements
**   29 GK3D6N elements
** create element sets: COVER1,COVER2,BASE1,BASE2,BOLT1,MATL1,BEAD1,
**                      COVER,BASE,COVERF1,COVERF2,COVERF3,COVERF4,
**                      COVERF5,COVERF6,BASEF1,BASEF3,BASEF4,BASEF5,
**                      BASEF6
** create node sets:    NALL,BASE_BOT
** create surfaces: SCOVER2,SBOLTL,SBOLTR,SCOV1L,SCOV1R,
**                             SGTOP,SGBOT,SBASE1
**
*INCLUDE,INPUT=manifoldgasket_thick_mesh.inp
**
*****
*GASKET SECTION, ELSET=BEAD1, BEHAVIOR=GBEAD1
*GASKET BEHAVIOR,NAME=GBEAD1
*GASKET THICKNESS BEHAVIOR,TYPE=ELASTIC-PLASTIC,DIRECTION=LOADING
0.0  ,0.0   
3.27 ,0.0254
8.40 ,0.0508
18.52,0.0762
34.26,0.1016
41.69,0.1118
45.50,0.13
47.80,0.15
*GASKET THICKNESS BEHAVIOR,TYPE=ELASTIC-PLASTIC,DIRECTION=UNLOADING
0.0,0.11,0.11
10.0,0.13,0.11
20.0,0.14,0.11
30.0,0.145,0.11
40.0,0.1475,0.11
47.8,0.15,0.11
*EXPANSION
1.2E-4,
*GASKET SECTION, ELSET=MATL1, BEHAVIOR=GMATL1
,0.076
*GASKET BEHAVIOR,NAME=GMATL1
*GASKET THICKNESS BEHAVIOR,TYPE=ELASTIC-PLASTIC,DIRECTION=LOADING
0.0  ,0.0   
4.83 ,0.0254
11.72,0.0508
25.17,0.0762
43.09,0.1016
52.75,0.1118
57.02,0.13
58.85,0.15
*GASKET THICKNESS BEHAVIOR,TYPE=ELASTIC-PLASTIC,DIRECTION=UNLOADING
0.0,0.09,0.09
10.0,0.108,0.09
20.0,0.122,0.09
30.0,0.131,0.09
40.0,0.142,0.09
50.0,0.147,0.09
58.85,0.15,0.09
*EXPANSION
3.E-5,
*****
*SOLID SECTION, ELSET=COVER, MATERIAL=ALUM
*SOLID SECTION, ELSET=BASE, MATERIAL=STEEL
*SOLID SECTION, ELSET=BOLT1, MATERIAL=STEEL
*****
*MATERIAL, NAME=STEEL
*ELASTIC
200000.,.28
*EXPANSION
1.6E-5,
*****
*MATERIAL, NAME=ALUM
*ELASTIC
71000.,.33
*EXPANSION
2.3E-5,
*****
**
** INITIAL TEMPERATURE -- 20 DEGREES CELSIUS
**
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL,20.0
** 
**
** CONTACT DEFINITIONS
**
*CONTACT PAIR,INTERACTION=IBOLT,SMALL SLIDING,ADJUST=.01
SBOLTL,SCOV1L
SBOLTR,SCOV1R
*CONTACT PAIR,INTERACTION=ITOP,SMALL SLIDING,ADJUST=.01
SGTOP,SCOVER2
*CONTACT PAIR,INTERACTION=IBOT,SMALL SLIDING, TIED,ADJUST=.01
SGBOT,SBASE1
*SURFACE INTERACTION,NAME=IBOLT
*SURFACE INTERACTION,NAME=ITOP
*SURFACE INTERACTION,NAME=IBOT
**
** PRE-TENSION SECTION DEFINITIONS
**
*ELSET,ELSET=RCUT
3701,3702,3703,3708
*ELSET,ELSET=LCUT
3709,3710,3711,3716
*surface,NAME=PREBOLTL
LCUT,   S1
*surface,NAME=PREBOLTR
RCUT,   S1
*PRE-TENSION SECTION,SURFACE=PREBOLTL,NODE=1000001
0.,0.,-1.
*PRE-TENSION SECTION,SURFACE=PREBOLTR,NODE=1000002
0.,0.,-1.
**
** SOFT SPRINGS ON COVER, BASE AND GASKET
** TO ELIMINATE RIGID BODY MOTION
**  - Nodes 643 and 2012 are located on COVER
**  - Nodes 5560 and 5689 are located on BASE 

**
*ELEMENT,TYPE=SPRING1,ELSET=SP1
999001,643
999002,2012
999003,5560
999004,5689
*ELEMENT,TYPE=SPRING1,ELSET=SP3
999005,643
999006,2012
*SPRING,ELSET=SP1
1, 
1., 
*SPRING,ELSET=SP3
3, 
1., 
**
** DISPLACEMENT BOUNDARY CONDITIONS 
**
*BOUNDARY
BASE_BOT,3
YSYMM,2
**
** NODE SETS FOR OUTPUT
**
*NSET,NSET=REFNODE
1000001,1000002
**********************************************
**
** STEP 1 -- PRE-TENSION LEFT BOLT
**
*STEP
*STATIC
0.01,1.0
*CLOAD
1000001,1,6000.
*BOUNDARY
1000002,1,,0.0
**
*MONITOR,NODE=15737,DOF=1
*NODE FILE,NSET=REFNODE,FREQ=1
U,RF,CF
*EL FILE,ELSET=IBEAD,FREQ=1
S,E,THE
*FILE FORMAT,ZERO INCREMENT
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*END STEP
**
** STEP 1A -- SUBSTITUTE LEFT BOLT LOAD WITH RESULTANT LENGTH CHANGE 
**
*STEP,INC=1
*STATIC
1.E-10,1.0E-10,1.0E-10,1.0E-10
*BOUNDARY,FIXED
1000001,1
*CLOAD,OP=NEW
*ENDSTEP
**
** STEP 2 -- PRE-TENSION RIGHT BOLT
**
*STEP
*STATIC
0.01,1.0
*CLOAD,OP=NEW
1000002,1,6000.
*BOUNDARY,OP=NEW,FIXED
1000001,1
*BOUNDARY,OP=NEW
BASE_BOT,3
YSYMM,2
*END STEP
**
** STEP 2A -- SUBSTITUTE RIGHT BOLT LOAD WITH RESULTANT LENGTH CHANGE 
**
*STEP,INC=1
*STATIC
1.E-10,1.0E-10,1.0E-10,1.0E-10
*BOUNDARY,FIXED
1000002,1
*CLOAD,OP=NEW
*ENDSTEP
**
** STEP 3 -- INCREASE TEMPERATURE TO 150 DEGREES CELSIUS 
**           and PRESSURIZE INTERIOR CAVITY TO 0.6895 MPa
**
*STEP
*STATIC
0.125,1.0
*TEMPERATURE
NALL,150.
*DLOAD,OP=NEW
coverf1,p1,0.6895
coverf2,p2,0.6895
coverf3,p3,0.6895
coverf4,p4,0.6895
coverf5,p5,0.6895
coverf6,p6,0.6895
basef1 ,p1,0.6895
basef3 ,p3,0.6895
basef4 ,p4,0.6895
basef5 ,p5,0.6895
basef6 ,p6,0.6895
*MONITOR,NODE=15737,DOF=1
*END STEP
**
** STEP 4 -- DECREASE TEMPERATURE TO -40 DEGREES CELSIUS
**
*STEP
*STATIC
0.125,1.0
*TEMPERATURE
NALL,-40.0
*END STEP
**
** STEP 5 -- INCREASE TEMPERATURE TO 20 DEGREES CELSIUS
**           and REMOVE INTERNAL PRESSURE LOAD
**
*STEP
*STATIC
0.05,1.0
*TEMPERATURE
NALL,20.0
*DLOAD,OP=NEW
*END STEP
**
** STEP 6 -- UNLOAD THE GASKET
**
*STEP
*STATIC
0.125,1.0
*BOUNDARY,OP=NEW
BASE_BOT,3
YSYMM,2
1000001,1,,0.005
1000002,1,,0.005
*END STEP