*HEADING
  THERMAL (*TEMPERATURE) LOAD VERIFICATION FOR COH3D8, COH3D6, COH2D4 & COHAX4
*PREPRINT, MODEL=YES
*NODE,NSET=COH3D8
*SYSTEM
*NODE,NSET=COH3D6
*SYSTEM
*NODE,NSET=COH2D4
*SYSTEM
*NODE,NSET=COHAX4
*ELEMENT,TYPE=COH3D8, ELSET=COH3D8
*ELEMENT,TYPE=COH3D6, ELSET=COH3D6
*ELEMENT,TYPE=COH2D4, ELSET=COH2D4
*ELEMENT,TYPE=COHAX4, ELSET=COHAX4
*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
*NSET, NSET=TOP_2D
*NSET, NSET=BOTTOM_3D
*NSET, NSET=BOTTOM_2D
*NSET, NSET=RIGHT_3D
*NSET, NSET=RIGHT_2D
*NSET, NSET=LEFT_3D
*NSET, NSET=LEFT_2D
*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