**
**  HGO model with 3 fibers, C3D8, C3D20R, C3D10 elements
**  theoretical result: s11 = 1486.
*Heading
 HGO Model
*Preprint, echo=NO, model=NO, history=NO, contact=NO
** ----------------------------------------------------------------
** 
** PART INSTANCE: Part-1-1
** 
*Node
      1,           1.,           1.,           1.
      2,           1.,           0.,           1.
      3,           1.,           1.,           0.
      4,           1.,           0.,           0.
      5,           0.,           1.,           1.
      6,           0.,           0.,           1.
      7,           0.,           1.,           0.
      8,           0.,           0.,           0.
*Element, type=C3D8
1, 5, 6, 8, 7, 1, 2, 4, 3
*Nset, nset=Part-1-1__PickedSet2, generate
 1,  8,  1
*Elset, elset=Part-1-1__PickedSet2
 1,
** Section: Section-1
*Solid Section, elset=Part-1-1__PickedSet2, material=holzapfel,orientation=ori-1
1.
*Nset, nset=y0, generate
 2,  8,  2
*Nset, nset=x0, generate
 5,  8,  1
*Nset, nset=z0
 3, 4, 7, 8
*Nset, nset=x1, generate
 1,  4,  1
*Nset, nset=z1
 1, 2, 5, 6
*Nset, nset=y1, generate
 1,  7,  2
** 
** MATERIALS
** 
*orientation,name=ori-1,local directions=3
 1.0,0.0,0.0,0.0,1.0,0.0
 3,0.0
 0.577350269, 0.577350269, 0.577350269
-0.577350269, 0.577350269,-0.577350269
 1,0,0
**
*solid section,elset=solid,material=holzapfel,orientation=ori-1
 1.0,
**
*parameter
 density=5.0e-4
 C10=7.64
 D=1.e-8
 K1=996.6
 K2=524.6
 Kappa=0.226
*material,name=holzapfel
*density
 <density>
*anisotropic hyperelastic,holzapfel, local directions=3
 <C10>, <D>, <K1>, <K2>, <Kappa>
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
x0, XSYMM
** Name: BC-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary
y0, YSYMM
** Name: BC-3 Type: Symmetry/Antisymmetry/Encastre
*Boundary
z0, ZSYMM
z1, ZSYMM
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*Step, name=Step-1, nlgeom=YES
*Static
0.1, 1., 1e-05, 0.1
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-4 Type: Displacement/Rotation
*Boundary
x1, 1, 1, 0.1
** 
** OUTPUT REQUESTS
** 
*Restart, write, frequency=0
** 
** FIELD OUTPUT: F-Output-1
** 
*Output, field
*element output
 S,E
*node output
 U,RF
** 
*End Step