*Heading
 Elastic-plastic analysis with FRAME2D      
 PIPE cross-section with user defined 
 generalized plastic loads
*node,nset=fixed
1,0.,0.
*node,nset=frame
2,0.,10.
3,10.,10.
4,10.,4.
*element, type=frame2d, elset=all
1,1,2
2,2,3
3,3,4
*frame section,section=pipe,elset=all
1.,.174355

3.e+6,1.15e+6
*plastic m1
21000.,0.
28000.,.002
28000.,.003
*plastic axial
50000.,0.
60000.,.002
60000.,.003
*boundary
1,encastre
*restart,write,freq=999
**********************************************
*step
first step: concentrated load at node 4
*static
0.1,1.0
*cload
4,1,-3800.
*node print, nset=fixed
rf, 
*node file, nset=fixed
rf, 
*node print, nset=frame
u, 
*node file,nset=frame
u, 
*el print, elset=all
sf,sep
*el file, elset=all
sf,sep
salpha,
*energy print
*energy file
*end step
***********************************************
*step
second step: concentrated load at node 2
*static
0.1,1.0
*cload
2,2,-30000.
*node print, nset=fixed
rf, 
*node file, nset=fixed
rf, 
*node print, nset=frame
u, 
*node file,nset=frame
u, 
*el print, elset=all
sf,sep
*el file, elset=all
sf,sep
salpha,
*energy print
*energy file
*end step
***********************************************
*step,inc=1000
third step: second load at node 4
*static
0.1,1.0
*cload
4,2,850.
*node print, nset=fixed
rf, 
*node file, nset=fixed
rf, 
*node print, nset=frame
u, 
*node file,nset=frame
u, 
*el print, elset=all
sf,sep
*el file, elset=all
sf,sep
salpha,
*energy print
*energy file
*end step
***********************************************