*HEADING
  THERMAL (*TEMPERATURE) LOAD VERIFICATION FOR COH3D8, COH3D6, COH2D4 & COHAX4
*PREPRINT, MODEL=YES
*NODE,NSET=COH3D8
11, 0., 0., 0.
12, 0., 0., 7.
13, 7., 0., 7.
14, 7., 0., 0.
15, 0., 7., 0.
16, 0., 7., 7.
17, 7., 7., 7.
18, 7., 7., 0.
*SYSTEM
0., 10., 0.
*NODE,NSET=COH3D6
21, 0., 0., 0.
22, 0., 0., 7.
23, 7., 0., 7.
24, 7., 0., 0.
25, 0., 7., 0.
26, 0., 7., 7.
27, 7., 7., 7.
28, 7., 7., 0.
*SYSTEM
0., 20., 0.
*NODE,NSET=COH2D4
31, 0., 0., 0.
32, 7., 0., 0.
33, 7., 7., 0.
34, 0., 7., 0.
*SYSTEM
0., 30., 0.
*NODE,NSET=COHAX4
41, 10., 0., 0.
42, 20., 0., 0.
43, 20., 10., 0.
44, 10., 10., 0.
*ELEMENT,TYPE=COH3D8, ELSET=COH3D8
11, 11, 12, 13, 14, 15, 16, 17, 18
*ELEMENT,TYPE=COH3D6, ELSET=COH3D6
21, 21, 22, 23, 25, 26, 27
22, 24, 21, 23, 28, 25, 27
*ELEMENT,TYPE=COH2D4, ELSET=COH2D4
31, 31, 32, 33, 34
*ELEMENT,TYPE=COHAX4, ELSET=COHAX4
41, 41, 42, 43, 44
*NSET, NSET=ALL
COH3D8, COH3D6, COH2D4, COHAX4
*ELSET, ELSET=EALL
COH3D8, COH3D6, COH2D4, COHAX4
*COHESIVE SECTION, THICKNESS=GEOMETRY, MATERIAL=COH, ELSET=EALL, RESPONSE=CONTINUUM
*MATERIAL,NAME=COH
*ELASTIC
3.E6, 0.3
*EXPANSION
 .0001, 
*DENSITY
 10.0, 
*BOUNDARY
 COH3D8,1,3
 COH3D6,1,3
 COH2D4,1,2
 COHAX4,1,2
*NSET, NSET=QA_TEST_ALL
ALL
*ELSET, ELSET=QA_TEST_EALL
EALL
*NSET, NSET=TOP_3D
15, 16, 17, 18, 25, 26, 27, 28
*NSET, NSET=TOP_2D
33, 34, 43, 44
*NSET, NSET=BOTTOM_3D
11, 12, 13, 14, 21, 22, 23, 24
*NSET, NSET=BOTTOM_2D
31, 32, 41, 42
*NSET, NSET=RIGHT_3D
13, 14, 17, 18, 23, 24, 27, 28
*NSET, NSET=RIGHT_2D
32, 33, 42, 43
*NSET, NSET=LEFT_3D
11, 12, 15, 16, 21, 22, 25, 26
*NSET, NSET=LEFT_2D
31, 34, 41, 44
*INITIAL CONDITIONS, TYPE=TEMPERATURE
ALL, 0.0
***************************
** 	STEP-1
***************************
*amplitude,name=amp1,definition=smooth
0,0,1,1
*amplitude,name=amp2,definition=smooth
0,1,1,0
*STEP, NAME=TEMP_ALL
*DYNAMIC, EXPLICIT
,1.
*TEMPERATURE,amplitude=amp1
 ALL,50.
*OUTPUT, FIELD, NUMBER INTERVAL=20, TIME MARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
 LE, TEMP
*NODE OUTPUT, NSET=QA_TEST_ALL
 U
*OUTPUT, HISTORY, VARIABLE=PRESELECT,TIME INTERVAL=1.0
*END STEP
***************************
** 	STEP-2
***************************
*STEP, NAME=RELEASE_STRESS
*DYNAMIC, EXPLICIT
,1.
*TEMPERATURE, OP=MOD,amplitude=amp2
 ALL, 50.
*END STEP
***************************
** 	STEP-3
***************************
*STEP, NAME=TEMP_TOP
*DYNAMIC, EXPLICIT
,1.
*BOUNDARY,OP=NEW, TYPE=VELOCITY
 BOTTOM_3D, 1, 3, 0.0
 BOTTOM_2D, 1, 2, 0.0
 TOP_3D, 1, 1, 0.0
 TOP_3D, 3, 3, 0.0
 TOP_2D, 1, 1, 0.0
*TEMPERATURE,OP=MOD,amplitude=amp1
 TOP_3D, 50.0
 TOP_2D, 50.0
*END STEP
***************************
** 	STEP-4
***************************
*STEP, NAME=TEMP_-TOP
*DYNAMIC, EXPLICIT
,1.
*TEMPERATURE,OP=MOD,amplitude=amp2
 TOP_3D, 50.0
 TOP_2D, 50.0
*END STEP
***************************
** 	STEP-5
***************************
*STEP, NAME=TEMP_RIGHT
*DYNAMIC, EXPLICIT
,1.
*TEMPERATURE,OP=MOD,amplitude=amp1
 RIGHT_3D, 50.0
 RIGHT_2D, 50.0
*END STEP
***************************
** 	STEP-6
***************************
*STEP, NAME=TEMP_-RIGHT
*DYNAMIC, EXPLICIT
,1.
*TEMPERATURE,OP=MOD,amplitude=amp2
 RIGHT_3D, 50.0
 RIGHT_2D, 50.0
*END STEP