*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 GK3D8 elements
**   29 GK3D6 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_mesh.inp
**
*****
*GASKET SECTION, ELSET=BEAD1, BEHAVIOR=GBEAD1
*GASKET BEHAVIOR,NAME=GBEAD1
*GASKET ELASTICITY,COMPONENT=MEMBRANE
75.,
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR
40.0,
*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 ELASTICITY,COMPONENT=MEMBRANE
105.,
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR
55.,
*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,
ADJUST=.01
SGBOT,SBASE1
*SURFACE INTERACTION,NAME=IBOLT
*FRICTION
0.20,
*SURFACE INTERACTION,NAME=ITOP
*FRICTION
0.20,
*SURFACE INTERACTION,NAME=IBOT
*SURFACE BEHAVIOR,NO SEPARATION
*FRICTION
0.20,
**
** 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
**  - Nodes 15827 and 15737 are located on GASKET
**
*ELEMENT,TYPE=SPRING1,ELSET=SP1
999001,643
999002,2012
999003,5560
999004,5689
999005,15827
999006,15737
*ELEMENT,TYPE=SPRING1,ELSET=SP3
999007,643
999008,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,Convert Sdi=No
*STATIC
0.01,1.0
*CLOAD
1000001,1,6000.
*BOUNDARY
1000002,1,,0.0
**
*MONITOR,NODE=15737,DOF=3
*NODE FILE,NSET=REFNODE,FREQ=0
U,TF,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=REFNODE
TF,
*ELEMENT OUTPUT,ELSET=IBEAD
S11,E11,THE11,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
U,
*ELEMENT OUTPUT,ELSET=GASKET
PE,S
*CONTACT OUTPUT
CSTRESS,CSTRESSERI
*EL FILE,ELSET=IBEAD,FREQ=0
S,E,THE
***FILE FORMAT,ZERO INCREMENT
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*END STEP
**
** STEP 2 -- PRE-TENSION RIGHT BOLT
**           SUBSTITUTE LEFT BOLT LOAD
**           WITH RESULTANT LENGTH CHANGE
**
*STEP,Convert Sdi=No
*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 3 -- INCREASE TEMPERATURE TO 150
**           DEGREES CELSIUS and PRESSURIZE
**           INTERIOR CAVITY TO 0.6895 MPa
**           SUBSTITUTE RIGHT BOLT LOAD
**           WITH RESULTANT LENGTH CHANGE
**
*STEP,Convert Sdi=No
*STATIC
0.125,1.0
*BOUNDARY,FIXED
1000002,1
*CLOAD,OP=NEW
*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,Convert Sdi=No
*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,Convert Sdi=No
*STATIC
0.05,1.0
*TEMPERATURE
NALL,20.0
*DLOAD,OP=NEW
*END STEP
**
** STEP 6 -- UNLOAD THE GASKET
**
*STEP,Convert Sdi=No
*STATIC
0.125,1.0
*BOUNDARY,OP=NEW
BASE_BOT,3
YSYMM,2
1000001,1,,0.005
1000002,1,,0.005
*END STEP