*heading
Implement isotropic hyperelastic behavior through VUANISOHYPER_INV
*****************************************
**
** Element types: c3d8r, cpe4r
** Material     : Incompressible polynomial, N=2
** Types of test: Plane-strain tension with orientation
** Results      : s11=6.741e2, s22=1.1685e2, s33=4.058e2, s12=-3.370e2
**
*****************************************
*node, nset=alln
** for c3d8r - abaqus isotropic hyper elasticity material
  1,0,0,0
  2,1,0,0
  3,1,1,0
  4,0,1,0
  5,0,0,1
  6,1,0,1
  7,1,1,1
  8,0,1,1
** for c3d8r - vuaniso hyper elasticity material
 11,0,2,0
 12,1,2,0
 13,1,3,0
 14,0,3,0
 15,0,2,1
 16,1,2,1
 17,1,3,1
 18,0,3,1
** for cpe4r - vuaniso hyper elasticity material
 21,0,4
 22,1,4
 23,1,5
 24,0,5
*nset,nset=z0
 1,2,3,4,11,12,13,14
*nset,nset=z1
 5,6,7,8,15,16,17,18
*nset,nset=x0
 1,4,5,8,11,14,15,18,21,24,25,28
*nset,nset=x1
 2,3,6,7,12,13,16,17,22,23,26,27
*nset,nset=y0
 1,2,5,6,11,12,15,16,21,22,25,26
*nset,nset=y1
 3,4,7,8,13,14,17,18,23,24,27,28
*element,type=c3d8r,elset=e_abq
 1,1,2,3,4,5,6,7,8
*element,type=c3d8r,elset=e_vuaniso3d
 2,11,12,13,14,15,16,17,18
*element,type=cpe4r,elset=e_vuaniso2d
 3,21,22,23,24
*elset, elset=eall
 1,2,3
*elset, elset=QA_TEST
 eall
*orientation,name=ori-1,local directions=2
 1,0.5,0,-1,1,0
 3,0
 1,0,0
 1,1,1
*parameter
 density=1.0e-3
*material,name=VUANISO_INVISO_1
*density
 <density>
*anisotropic hyperelastic,user,formulation=invariant, local directions=2,type=incompressible,properties=7
 100.,50.,10.,20.,30.,1.0e-2,0.0
*material,name=polyabq
*density
 <density>
*hyperelastic, n=2
 100.,50.,10.,20.,30.,1.0e-2,0.0
*solid section, elset=e_abq, material=polyabq,orientation=ori-1
 1.
*solid section, elset=e_vuaniso3d, material=VUANISO_INVISO_1,orientation=ori-1
 1.
*solid section, elset=e_vuaniso2d, material=VUANISO_INVISO_1,orientation=ori-1
 1.
*boundary
 z0,3
 z1,3
 y0,2
 x0,1
*amplitude,name=bc,definition=smooth step
 0.0,0.0,1.0,1.0
*step,nlgeom=yes
*dynamic,explicit
 , 1.
*boundary,amp=bc
 x1,1,1,0.5
*output,field,number interval=1
*element output,elset=QA_TEST
 S,E
*node output
 U,RF
*output,history
*element output, elset=eall
 S11,s22,s12,s33,le11
*end step