** **************************************************************
**  Abaqus  6.14-1
********
*HEADING
Bilame Invar-Steel
beam: section 5x10 mm2, length 100 mm
*SYSTEM
***** units: mm, N, kg, s, MPa
*node
*ngen,nset=nn0
***** dx = 20 mm/8 = 2.25 mm et dy = 10/4 = 2.5 mm
*ngen,nset=nn1
**** section in xy plane
*nfill,nset=n0
nn0,nn1,4,10
****
*ncopy,change number=10000,old set=n0, new set=n100,shift
*nfill,nset=nall,bias=1.
n0,n100,100,100
******************************
*element,type=c3d8t
****
*elgen,elset=eall
********************************
*elset,elset=esteel,generate
*elset,elset=einvar,generate
*elset,elset=etop,generate
*elset,elset=eunder,generate
*Surface, type=ELEMENT, name=top
etop, S5
*Surface, type=ELEMENT, name=under
eunder, S3
*nset,nset=nmes,generate
*elset,elset=emes
*nset,nset=nfix,generate
*******
*SOLID SECTION, ELSET=esteel, MATERIAL=steel
*material, name=steel
*Conductivity
30.e-3
*Density
 7.8e-06
*Elastic
210000., 0.3
*Expansion
 1e-05
*Specific Heat
 1000.
*SOLID SECTION, ELSET=einvar, MATERIAL=invar
*material, name=invar
*Conductivity
**** W/mmK
10.e-3
*Density
***** Kg/mm3
8e-06
*Elastic
******  MPa
141000., 0.3
*Expansion
1e-06
*Specific Heat
***** J/KgK
500.
********************************************
*initial conditions,type=temperature
nall, 20.
**********************************
*STEP, inc=1000
heating during 100 sec
******
*Coupled Temperature-displacement,deltmx=2.
0.1, 100.
**************************
*boundary
**** fixed nodes
n0,encastre
********* no gravity
*dload
eall,grav,10,0,-1
******  heating on surfaces under and top
*Sfilm
***** h = 100 W/m2K = 100.e-6 W/mm2K and Tair = 100.
under, F, 100., 100.e-6
top, F, 100., 100.e-6
****************************
*el print,position=centroidal,elset=emes,freq=100
mises,s11,s22,s33
*node print,nset=nmes,freq=100
coord,u1,u2,u3
************************************************
*OUTPUT,FIELD,freq=1,var=preselect
*output,history,freq=1,var=preselect
************************************
*monitor,node=10001,dof=11
*END STEP
*************