** **************************************************************
**  Abaqus  6.14-1
********
*HEADING
Bending of cantilever beam: verification of elastic distorsion
beam: section 10x20 mm2, length 200 mm, load : 100 kPa
*SYSTEM
***** units: mm, N, kg, s, MPa
***** quadratique elements with plane of symmetry
***** x = 0 est plan de symmetrie: node set = nx
*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=c3d20
**** quadratic element with 20 nodes
*elgen,elset=eall
********************************
*elset,elset=etop,generate
*Surface, type=ELEMENT, name=top
*****, internal
etop,S5
*nset,nset=nmes,generate
*nset,nset=nx,generate
*elset,elset=emes
*******
*SOLID SECTION, ELSET=eall, MATERIAL=alu
****
*MATERIAL,NAME=alu
*density
***** en kg/mm3
2.7e-6
*elastic
75.e3,0.3
***************************************************
*STEP,name=loading 0.1 MPa,nlgeom=no
Bending of cantilever beam: loading 0.2 MPa
******
*static
****dt0, step time, min, max
0.01,1.
**************************
*boundary
**** fixed nodes
n0,1,3,0.
**** symmetry plane
nx,xsymm
*dsload
***** 100 kPa = 0.1 MPa
top,p,0.2
******
*dload
****** garvity alon -ey
eall,grav,10.,0.,-1.,0.
****** centrifugal load
***Eall,centrif,1.e4,0.,5.,0.,1.,0.,0.
*****
****************************
*el print,position=nodes,elset=emes,freq=100
mises,s11,s22,s33
*node print,nset=nmes,freq=10
coord,u1,u2,u3
************************************************
*OUTPUT,FIELD,freq=1,var=preselect
*output,history,freq=1,var=preselect
************************************
*monitor,node=10001,dof=2
*END STEP
*************
***************************************************
*STEP,name=unloading 0 MPa,nlgeom=no
Bending of cantilever beam: unloading
******
*static
****dt0, step time, min, max
0.01,1.
**************************
*boundary
**** fixed nodes
n0,1,3,0.
**** symmetry plane
nx,xsymm
*dsload
*****
top,p,0.
******
*dload
****** garvity alon -ey
eall,grav,10.,0.,-1.,0.
****** centrifugal load
***Eall,centrif,1.e4,0.,5.,0.,1.,0.,0.
****************************
*END STEP
*************