*HEADING
 Test import capability with PRF
 Base standard analysis
*RESTART,WRITE,num=20
*PART,NAME=PART_3D
  *NODE,NSET=ALL
   1, 0., 0.
   2, 1., 0.
   3, 1., 1.
   4, 0., 1.
   5, 0., 0. ,1.0
   6, 1., 0. ,1.0
   7, 1., 1. ,1.0
   8, 0., 1. ,1.0
  *NSET,NSET=NFIX1
   1,4,5,8
  *NSET,NSET=NFIX2
   1,2,5,6
  *NSET,NSET=NFIX3
   1,2,3,4
  *NSET,NSET=NLOAD1
   2,3,6,7
  **
  *ELEMENT,TYPE=C3D8R,ELSET=SOLID
   1,1,2,3,4,5,6,7,8
  *ELSET,ELSET=EALL
   SOLID
*END PART
*PART,NAME=PART_PS
  *NODE,NSET=ALL
   1, 0., 0.
   2, 1., 0.
   3, 1., 1.
   4, 0., 1.
  *NSET,NSET=NFIX1
   1,4
  *NSET,NSET=NFIX2
   1,2
  *NSET,NSET=NLOAD1
   2,3
  **
  *ELEMENT,TYPE=CPE4R,ELSET=SOLID
   1,1,2,3,4
  *ELSET,ELSET=EALL
   SOLID
*END PART
*PART,NAME=PART_AXI
  *NODE,NSET=ALL
   1, 0., 0.
   2, 1., 0.
   3, 1., 1.
   4, 0., 1.
  *NSET,NSET=NFIX1
   1,4
  *NSET,NSET=NFIX2
   1,2
  *NSET,NSET=NLOAD1
   2,3
  **
  *ELEMENT,TYPE=CAX4R,ELSET=SOLID
   1,1,2,3,4
  *ELSET,ELSET=EALL
   SOLID
*END PART
*ASSEMBLY,NAME=MODEL
  *INSTANCE,NAME=INST_3D,PART=PART_3D
    *SOLID SECTION, ELSET=SOLID, MATERIAL=VEP_NH_PRF
     1.
  *END INSTANCE
  *INSTANCE,NAME=INST_PS,PART=PART_PS
     2.0,0.0,0.0
    *SOLID SECTION, ELSET=SOLID, MATERIAL=VEP_AB_PRF
     1.
  *END INSTANCE
  *INSTANCE,NAME=INST_AXI,PART=PART_AXI
     4.0,0.0,0.0
    *SOLID SECTION, ELSET=SOLID, MATERIAL=VEP_NH_PRF
     1.
  *END INSTANCE
  **
  *NSET,NSET=NFIX1
   INST_3D.NFIX1,INST_PS.NFIX1,INST_AXI.NFIX1
  *NSET,NSET=NFIX2
   INST_3D.NFIX2,INST_PS.NFIX2,INST_AXI.NFIX2
  *NSET,NSET=NFIX3
   INST_3D.NFIX3
  *NSET,NSET=NLOAD1
   INST_3D.NLOAD1,INST_PS.NLOAD1,INST_AXI.NLOAD1
  *ELSET,ELSET=EALL
   INST_3D.EALL,INST_PS.EALL,INST_AXI.EALL
*END ASSEMBLY
**---------------------------------------------
**
** MATERIAL DEFINITIONS
**
*MATERIAL,NAME=VEP_NH_PRF
*DENSITY
 1.e-6
*HYPERELASTIC,NEO HOOK,MODULI=INSTANT
 10,0.025
*Mullins effect
 ** R       M             beta
 1.9057,  3.77307E-02,  0.59222
*Plastic
 0.19668078,  0
 2.9532671, 0.012512393
 9.5576679, 0.044332625
 92.370193, 0.44332625
*viscoelastic,nonlinear,law=strain,networkid=1,sratio=0.3
 0.1,1.0,0.0
*viscoelastic,nonlinear,law=power law,networkid=2,sratio=0.25
 1.0,1.0,0.0,0.0,1.0
*viscoelastic,nonlinear,law=hyperb,networkid=3,sratio=0.1
 0.03,0.15,2.0
*MATERIAL,NAME=VEP_AB_PRF
*DENSITY
 1.e-6
*HYPERELASTIC,ARRUDA-BOYCE,MODULI=INSTANT
 20,7.0,0.025
*Mullins effect
 ** R       M             beta
 1.9057,  3.77307E-02,  0.59222
*Plastic
 0.19668078,  0
 2.9532671, 0.012512393
 9.5576679, 0.044332625
 92.370193, 0.44332625
*viscoelastic,nonlinear,law=strain,networkid=1,sratio=0.3
 0.1,1.0,0.0
*viscoelastic,nonlinear,law=strain,networkid=2,sratio=0.25
 0.03,2.0,-0.1
*viscoelastic,nonlinear,law=powerlaw,networkid=3,sratio=0.1
 0.03,2.0,-0.1
*MATERIAL,NAME=VEP_OGDEN_PRF
*DENSITY
 1.e-6
*HYPERELASTIC,MOONEY-RIVLIN,MODULI=INSTANT
 20,5.0,0.025
*Mullins effect
 ** R       M             beta
 1.9057,  3.77307E-02,  0.59222
*Plastic
 0.19668078,  0
 2.9532671, 0.012512393
 9.5576679, 0.044332625
 92.370193, 0.44332625
*viscoelastic,nonlinear,law=strain,networkid=1,sratio=0.3
 0.1,1.0,0.0
*viscoelastic,nonlinear,law=strain,networkid=2,sratio=0.25
 0.03,2.0,-0.1
*viscoelastic,nonlinear,law=powerlaw,networkid=3,sratio=0.1
 0.03,2.0,-0.1
*amplitude,name=nominal_strain
0.0, 0.0, 0.5, 1.0, 1.0, 0.05
*STEP,nlgeom=yes
*DYNAMIC,EXPLICIT
 ,1.
*BOUNDARY
 NFIX1,1,1
 NFIX2,2,2
 NFIX3,3,3
*boundary, amplitude=nominal_strain
 NLOAD1,1,1,0.4
*OUTPUT,FIELD,NUMBER=4
*NODE OUTPUT
 U,RF
*ELEMENT OUTPUT
 S,PEEQ,CEEQ,CE,PE,E
*OUTPUT,HISTORY,TIME INTERVAL=0.01
*ELEMENT OUTPUT,ELSET=EALL
 PEEQ,MISES,PRESS,CEEQ
*END STEP