*heading
 clay plasticity, uniaxial compression test, cax4r
 test 9-digit labels
*preprint, echo=yes,hist=yes,model=yes
*restart,write
*node,nset=soil
*nset,nset=bottom
*nset,nset=top
*nset,nset=top1
*nset,nset=top2
*nset,nset=top3
*nset,nset=soil
*nset,nset=soil1
*nset,nset=soil2
*nset,nset=soil3
*nset,nset=lhs
*element,type=cax4r,elset=soil
*element,type=cax4r,elset=soil1
*element,type=cax4r,elset=soil2
*element,type=cax4r,elset=soil3
*elset,elset=eout
soil,soil1,soil2,soil3
*solid section,material=sample,elset=soil
*solid section,material=sample1,elset=soil1
*solid section,material=sample1,elset=soil2
*solid section,material=sample2,elset=soil3
*surface,name=side
900000001,s3
*surface,name=side1
900000100,s3
*surface,name=side2
900000200,s3
*surface,name=side3
900000300,s3
*material,name=sample
*density
 1.0e-1
*elastic
 30000.,.3
*clay plasticity,hardening=tab
*clay hardening
*material,name=sample1
*density
 1.0e-1
*elastic
 30000.,.3
*clay plasticity,hardening=tab
*clay hardening,dependencies=1
*material,name=sample2
*density
 1.0e-1
*elastic
 30000.,.3
*clay plasticity,hardening=tab
*clay hardening,dependencies=1
*initial conditions,type=ratio
soil,1.08,0.,1.08,1.
soil1,1.08,0.,1.08,1.
soil2,1.08,0.,1.08,1.
soil3,1.08,0.,1.08,1.
*initial conditions,type=stress,geostatic
soil,-5.,0.,-5.,10.,1.
soil1,-5.,0.,-5.,10.,1.
soil2,-5.,0.,-5.,10.,1.
soil3,-5.,0.,-5.,10.,1.
*initial conditions,type=field,variable=1
soil2,1
*initial conditions,type=field,variable=1
soil3,1
*step, nlgeom=no
  geostatic initial stress state
*dynamic,explicit
 ,1.0e-3
*boundary
top,2
top1,2
top2,2
top3,2
bottom,2
lhs,1
*field,variable=1
soil1,1
soil2,1
soil3,1
*dload
900000001,p2,1.
900000100,p2,1.
900000200,p2,1.
900000300,p2,1.
*output,field
*element output
s,e,pe
*output,history
*element output, elset=eout
pe,peeq,press,mises
*end step
*amplitude, name=amp-1
 0.0,0.0,20.0,1.0
*step, nlgeom=no
  uniaxial compression
*dynamic,explicit
 ,20.
*boundary,amplitude=amp-1
top,2,2,-.5
*end step
*step, nlgeom=no
  uniaxial compression 1
*dynamic,explicit
 ,20.
*field,variable=1
soil1,1
*boundary,amplitude=amp-1
top1,2,2,-.5
*end step
*step, nlgeom=no
  uniaxial compression 2
*dynamic,explicit
 ,20.
*boundary,amplitude=amp-1
top2,2,2,-.5
*end step
*step, nlgeom=no
  uniaxial compression 3
*dynamic,explicit
 ,20.
*boundary,amplitude=amp-1
top3,2,2,-.5
*end step