*Heading
 Elastic-plastic analysis with FRAME2D 
 PERTURBATION/LOAD CASE     
 BOX section, user defined generalized plastic forces
*node,nset=boun
1,0.,0.
3,10.,0.
*node,nset=load
2,5.,0.
*element, type=frame2d, elset=all
1,1,2
2,2,3
*frame section,section=box,elset=all
1.,1.,.2,.2,.2,.2

3.e+6,1.15e+6
*plastic m1
3300.,.0
4200.,.0015
4200.,.002
*plastic axial
32000.,.0
40000.,.002
40000.,.003
*boundary
1,1,2
3,2
*restart,write,freq=999
**********************************************
*step
first step: concentrated loads in the middle node
*static
0.1,1.0
*cload
2,1,-20000.
2,2,-1000.
*node print, nset=boun
rf, 
*node file, nset=boun
rf, 
*node print, nset=load
u, 
*node file,nset=load
u, 
*el print, elset=all
sf,sep
*el file, elset=all
sf,sep
salpha,
*energy print
*energy file
*end step
**********************************************
*step
second step: bending moments at the boundaries.
*static
0.1,1.0
*cload
1,6,-1000.
3,6,1000.
*node print, nset=boun
rf, 
*node file, nset=boun
rf, 
*node print, nset=load
u, 
*node file,nset=load
u, 
*el print, elset=all
sf,sep
*el file, elset=all
sf,sep
salpha,
*energy print
*energy file
*end step
***********************************************
*STEP,PERTURBATION
*STATIC
*LOAD CASE,NAME="CLOAD IN MIDDLE NODE"
*cload
2,1,-20000.
2,2,-1000.
*END LOAD CASE
*LOAD CASE,NAME="BENDING MOMENT AT BOUNDARIES"
*cload
1,6,-1000.
3,6,1000.
*END LOAD CASE
*END STEP