** **************************************************************
**  Abaqus  6.14-1
********
*HEADING
Bending of cantilever beam: verification of elastic distorsion
beam: section 10x20 mm2, length 200 mm, load : 250 kPa
*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=c3d8
****
*elgen,elset=eall
********************************
*elset,elset=etop,generate
*elset,elset=e100mm,generate
*Surface, type=ELEMENT, name=top
etop, S5
*nset,nset=nmes,generate
*elset,elset=emes
*nset,nset=nfix,generate
*******
*SOLID SECTION, ELSET=eall, MATERIAL=copper
****
*MATERIAL,NAME=copper
*elastic
110.e3,0.32
*density
***** en kg/mm3
8.96e-6
*plastic,hardening=johnson cook
*** pure copper
90.,292.,0.31,0.
***************************************************
*************************************
*STEP,nlgeom=yes,amp=ramp
250 kPa with johnson cook for pure copper
*static
****dt0, step time, min, max
0.1,1.,1.e-5
**************************
*boundary
**** fixed nodes
n0,encastre
*******
*dsload
top,p,0.25
******
****************************
*el print,position=centroidal,elset=emes,freq=100
mises,s11,s22,s33,peeq
*node print,nset=nmes,freq=10
coord,u1,u2,u3
************************************************
*OUTPUT,FIELD,freq=1,var=preselect
*output,history,freq=1,var=preselect
*element output,elset=emes
mises,peeq
************************************
*monitor,node=10001,dof=2
*END STEP
*************
************************
*STEP,nlgeom=yes,amp=ramp
remove pressure
*static
****dt0, step time, min, max
0.1,1.,1.e-5,1.
**************************
*dsload
top,p,0.
*END STEP
*************