*HEADING
FREEZING OF A SQUARE PLATE  (C3D8RT)
**
** The target is temperature history at FOUT nodes
** see ch. 1.6.2 of  ABAQUS/Benchmarks Manual
**
*NODE
101,0.,0.
109,4.,0.
901,0.,4.
909,4.,4.
1101,0.,0. ,0.1
1109,4.,0. ,0.1
1901,0.,4. ,0.1
1909,4.,4. ,0.1
*NGEN,NSET=Y0
101,109
1101,1109
*NGEN,NSET=Y4
901,909
1901,1909
*NGEN,NSET=X0
101,901,100
1101,1901,100
*NGEN,NSET=X4
109,909,100
1109,1909,100
*NFILL
Y0,Y4,8,100
*ELEMENT,TYPE=C3D8RT
101,101,102,202,201 ,1101,1102,1202,1201
*ELGEN,ELSET=ALL
101,8,1,1,8,100,100
*SOLID SECTION,ELSET=ALL,MATERIAL=A1
*MATERIAL,NAME=A1
*ELASTIC
1.e6,0.3
*EXPANSION
0.0,
*DENSITY
1.,
*USER MATERIAL, TYPE=THERMAL, CONSTANTS=5
** conductivity, specific heat, latent heat, solidus temp., liquidus temp.
 1.08, 1., 70.26, -.25, -.15
*AMPLITUDE,NAME=RAMP,TIME=TOTAL TIME,VALUE=ABSOLUTE
0.,0.,.05,-45.,5.,-45.
*NSET,NSET=FOUT
 303,305
*NSET,NSET=QA_TEST
 303,305
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,5.
*BOUNDARY,AMPLITUDE=RAMP
 X0,11,11,-45.
 Y0,11,11,-45.
*BOUNDARY,TYPE=VELOCITY
 Y0, 1,3, 0.0
***************************
*FILE OUTPUT, NUM=4
*NODE FILE,NSET=FOUT
NT, 
***************************
*OUTPUT,FIELD,NUM=100
*NODE OUTPUT,NSET=QA_TEST
NT, 
**************************
***OUTPUT,HISTORY,FREQ=1000
***NODE OUTPUT,NSET=FOUT
***NT,
*******
*OUTPUT,HISTORY,TIME INTERVAL=1.25 
*ENERGY OUTPUT, VARIABLE=ALL
***************************
*END STEP