*************************************************
*Heading
** Job name: Tmeca Model on a piston
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
************ units: mm,N,MPa, J et °C
***** read mesh
*include,input=piston-mesh.dat
*************
** MATERIALS
**
*Material, name=alu
*Conductivity
 0.2,
*Density
 2.7e-06,
*Elastic
70000., 0.3
*Expansion
 3e-05,
*Plastic
150.,0.,20.
170., 0.02,20.
200., 0.05,20.
300.,  0.1,20.
50.,0.,420.
70., 0.02,420.
100., 0.05,420.
200.,  0.1,420.
*Specific Heat
420.
*****************************************************
*initial conditions, type=temperature
*****  Set-4 = nall
Set-4,20.
*****************************************************
*Step, name=Step-1, nlgeom=NO, inc=1000
heating
*Coupled Temperature-displacement, deltmx=5.
0.05, 3600.
**
** BC's
**
** Name: BC-Tfix Type: Temperature
*Boundary
Set-1, 11, 11,20.
** Name: BC-fixnodes
*Boundary
Set-3,1,3,0.
**
** Surface heat flux
*Dsflux
Surf-1, S, 0.2
**
** results
*Output, field, frequency=200
*Node Output
NT
*Element Output
HFL, LE, PE, PEEQ, PEMAG, S
**
*Output, history
*Node Output, nset=Set-1
NT
*Element Output, elset=Set-1
mises,temp,peeq,HFLM
S11,S22,S33
**
*Output, history, variable=PRESELECT, frequency=200
*End Step
*****************************************************
** STEP: Step-2
*Step, name=Step-2, nlgeom=NO, inc=1000
cooling
*Coupled Temperature-displacement,deltmx=5.
0.05, 3600.
**
** Surface heat flux
*Dsflux
Surf-1, S, 0.
**
*End Step