*HEADING
Round Notch Bar:CAX4; low cycle fatigue for ductile material
using direct cyclic approach
**preprint,model=yes,history=yes
*RESTART,WRITE,FREQ=999
*NODE
1,3,0.
33,5,2
1601,0,0.
1617,0,5
1633,5,5
10000,5,0.
2617,0,37.5
2633,5,37.5
*NGEN,LINE=C
1,33,1,10000
*NGEN
1601,1617,1
1617,1633,1
2617,2633,1
*NSET,NSET=DOWN,GENERATE
1,33,1
*NSET,NSET=UP,GENERATE
1601,1633,1
*NFILL,BIAS=1.0
DOWN,UP,16,100
*NSET,NSET=DOWN1,GENERATE
1617,1633
*NSET,NSET=UP1,GENERATE
2617,2633
*NFILL,BIAS=0.75
DOWN1,UP1,10,100
*NSET,NSET=TOP,GENERATE
2617,2633,1
*NSET,NSET=BOTTOM,GENERATE
1,1601,100
*NSET,NSET=LEFT,GENERATE
1601,1607,1
1617,2617,100
*NSET,NSET=ONE
1601
*nset,nset=nall,generate
1,1633 
*ELEMENT,TYPE=CAX4
1,101,1,2,102
1000,1617,1618,1718,1717
*ELGEN,ELSET=ALLE
1,16,100,1,32,1,16
1000,10,100,1,16,1,10
*ELSET,ELSET=TOPE,GENERATE
1009,1159,10
*elset,elset=junk
498,414
*BOUNDARY
BOTTOM,2
1601,1,1
*MATERIAL,NAME=BIDERMAN
*ELASTIC
**1.86E5, .33
2.0E5,0.33
*PLASTIC, HARDENING=KINEMATIC
400,0
800,2.4
*damage initiation,criterion=hysteresisenergy
10,-1.44
**26800,-1.44
*damage evolution, type=hysteresisenergy
3.5768E-5,1.15
*SOLID SECTION,MATERIAL=BIDERMAN,ELSET=ALLE
1.5,
*AMPLITUDE,NAME=AMP,definition=periodic
1,0.0785,0,0.125
0,0.0625
*STEP,inc=100
*Static
0.1,1
*boundary
top,2,2,0.125
*output,field,freq=100
*node output
u
*element output
s
e
pener
*output, history
*node output,nset=top
u2
*endstep
*STEP,inc=500
*DIRECT CYCLIC,fatigue
.5,80,,,30,30,5,100
3,6,50,1.0
*CONTROLS,TYPE=DIRECT CYCLIC
,100,5E-3,100,5E-3
*boundary,AMP=AMP
top,2,2,1
*output,field,freq=160
*node output
u
*element output
sdeg
status
cycleini,
s
e
pener
*output, history,freq=5
*node output,nset=top
u2
*END STEP