*HEADING
  DLOAD, DSLOAD  AND *TEMPERATURE FOR COHAX4
*PREPRINT, MODEL=YES
*NODE,NSET=ALL
1, 10., 0.
2, 20., 0.
3, 20., 10.
4, 10., 10.
*NSET,NSET=TOP
3, 4
*NSET,NSET=BOTTOM
1, 2
*NSET,NSET=RIGHT
2, 3
*NSET,NSET=LEFT
1, 4
*ELEMENT,TYPE=COHAX4, ELSET=EALL
1, 1, 2, 3, 4
*SURFACE, NAME=SIDE1
1, S1
*SURFACE, NAME=SIDE2
1, S2
*SURFACE, NAME=SIDE3
1, S3
*SURFACE, NAME=SIDE4
1, S4
*COHESIVE SECTION, THICKNESS=GEOMETRY, MATERIAL=COH, ELSET=EALL, RESPONSE=CONTINUUM
*MATERIAL,NAME=COH
*ELASTIC
3.E6, 0.3
*DENSITY
 10.0, 
*EXPANSION
 .0001, 
*INITIAL CONDITIONS, TYPE=TEMPERATURE
ALL, 0.0
*BOUNDARY
 ALL,1,2
*NSET, NSET=QA_TEST_ALL
ALL
*ELSET, ELSET=QA_TEST_EALL
EALL
***************************
** 	STEP-1
***************************
*STEP,PERTURBATION, NAME=BR
*STATIC
*DLOAD
 1,BR,3000.
*OUTPUT, FIELD
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
 S,BF
*NODE OUTPUT, NSET=QA_TEST_ALL
 RF
*OUTPUT, HISTORY, FREQUENCY=0
*EL PRINT
 LOADS,
*NODE PRINT
 RF, 
*END STEP
***************************
** 	STEP-2
***************************
*STEP,PERTURBATION, NAME=BZ
*STATIC
*DLOAD
 1,BZ,3000.
*END STEP
***************************
** 	STEP-3
***************************
*STEP,PERTURBATION, NAME=GRAV
*STATIC
*DLOAD
 1,GRAV,300.0,0.0,1.0,0.0
*END STEP
***************************
** 	STEP-4
***************************
*STEP,PERTURBATION, NAME=CENT
*STATIC
*DLOAD
1, CENT, 300., 0.0, 0.0, 0.0, 0., 1., 0.
*END STEP
***************************
** 	STEP-5
***************************
*STEP,PERTURBATION, NAME=CENTRIF
*STATIC
*DLOAD
1, CENTRIF, 30., 0.0, 0.0, 0.0, 0., 1., 0.
*END STEP
***************************
** 	STEP-6
***************************
*STEP,PERTURBATION, NAME=P1
*STATIC
*DLOAD
 1,P1,3000.
*END STEP
***************************
** 	STEP-7
***************************
*STEP,PERTURBATION, NAME=P2
*STATIC
*DLOAD
 1,P2,3000.
*END STEP
***************************
** 	STEP-8
***************************
*STEP,PERTURBATION, NAME=P3
*STATIC
*DLOAD
 1,P3,3000.
*END STEP
***************************
** 	STEP-9
***************************
*STEP,PERTURBATION, NAME=P4
*STATIC
*DLOAD
 1,P4,3000.
*END STEP
***************************
** 	STEP-10
***************************
*STEP,PERTURBATION, NAME=S1
*STATIC
*DSLOAD
 SIDE1,P,3000.
*END STEP
***************************
** 	STEP-11
***************************
*STEP,PERTURBATION, NAME=S2
*STATIC
*DSLOAD
 SIDE2,P,3000.
*END STEP
***************************
** 	STEP-12
***************************
*STEP,PERTURBATION, NAME=S3
*STATIC
*DSLOAD
 SIDE3,P,3000.
*END STEP
***************************
** 	STEP-13
***************************
*STEP,PERTURBATION, NAME=S4
*STATIC
*DSLOAD
 SIDE4,P,3000.
*END STEP
***************************
** 	STEP-14
***************************
*STEP,PERTURBATION, NAME=TEMP_ALL
*STATIC
*TEMPERATURE
 ALL,50.
*OUTPUT, FIELD
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
TEMP, S, E
*NODE OUTPUT, NSET=QA_TEST_ALL
U, RF
*EL PRINT
 S, TEMP
*NODE PRINT
 RF, 
*END STEP
***************************
** 	STEP-15
***************************
*STEP, NAME=TEMP_TOP
*STATIC
*BOUNDARY,OP=NEW
 BOTTOM, 1, 2, 0.0
 TOP, 1, 1, 0.0
*TEMPERATURE,OP=MOD
 TOP, 50.0
*OUTPUT, FIELD
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
TEMP, S, E
*NODE OUTPUT, NSET=QA_TEST_ALL
U, RF
*EL PRINT
 S, TEMP
*NODE PRINT
 RF, 
*END STEP
***************************
** 	STEP-16
***************************
*STEP, NAME=TEMP_-TOP
*STATIC
*TEMPERATURE,OP=MOD
 TOP, 0.0
*END STEP
***************************
** 	STEP-17
***************************
*STEP, NAME=TEMP_RIGHT
*STATIC
*TEMPERATURE,OP=MOD
 RIGHT, 50.0
*END STEP
***************************
** 	STEP-18
***************************
*STEP, NAME=TEMP_-RIGHT
*STATIC
*TEMPERATURE,OP=MOD
 RIGHT, 0.0
*END STEP