** ******************************************
**  Abaqus  6.20
**************************
*HEADING
Intro to ABAQUS heat flow calculation: Al casting
radius:30 mm heigth 150 mm, coarse mesh with linear elements
*SYSTEM
************
***** units:  kg,mm,J,s,K
******** mesh with ca 3 x 3 mm2 elements
******** lower radius: 30 mm
******** upper radius: 60 mm
******** height: 150 mm
**** Axisymetric Mesh (r-z): nodes, elements and sets
*node
**** number,r,z
*** equivalent to 1,0.,0.
*ngen,nset=n0
***** ngen = node generate
*ngen,nset=n50
*nfill,nset=nall
************* 50 layers of 3 mm
n0,n50,50,100
******** element
*element,type=dcax4
*** dcax4 = axi elements with 4 nodes for diffusion pb
*elgen,elset=eall
***** elgen = element generate
******** element and node sets for BCs
*elset,elset=elow,generate
*elset,elset=etop,generate
*nset,nset=nmes
*elset,elset=emes
*******
*SOLID SECTION, ELSET=eall,MATERIAL=alu
*** affect mesh elements to a material
****
*MATERIAL,NAME=alu
*DENSITY
**** en kg/mm3
2.7e-6
*CONDUCTIVITY
*** W/mmK
*** value, temp
50.e-3,20.
55.e-3,600.
*SPECIFIC HEAT
*** J/kgC
*** value, temp
980.,20.
900.,500.
*latent heat
**** in J/kg  400.e3
400.e3,600.,650.
*film property,name=compact
**** htc W/mm2K, temp
********************************************************
*INITIAL CONDITIONS,TYPE=TEMPERATURE
nall,700.
*************
*** cooling to 100°C
*************
*STEP,inc=10000,amp=step
cooling to 100°C in 500 sec
******
*heat transfer,DELTMX=1.
***********   dt0,step-time,min,max
0.01,500.,1e-08
**************************
*film
***** condition mixte: flux = h*(tsurf-tex)
***** h en W/mm2K
** elements,film,tex,h
elow,f1,100.,200.e-6
etop,f2,100.,500.e-6
**********************************
***** results:
**** follow one dof at a node (result in *.sta file)
*monitor,node=11,dof=11
**** print results in *.dat file:
*el print,position=centroidal,elset=emes,freq=100
temp,hfl
*el print,elset=emes,freq=0,position=averaged at nodes
temp,hfl
*el print,position=averaged at nodes,elset=eall,freq=0
temp,hfl
*node print,nset=nmes,freq=10
coord,nt11
*************************************
**** save field results in *.odb file:
*OUTPUT,FIELD,freq=10,var=preselect
*element output,elset=eall
temp
*node output,nset=nall
nt
**** save history results in *.odb file:
*output,history,freq=1
*node output,nset=nmes
nt
**********************************
*END STEP
*************