*HEADING
**--------------------------------------------------------------------**
**--------------------------------------------------------------------**
**Parameters
*PARAMETER
**--------------------------------------**
  Xl0      =  1.0
  Rate     =  0.5
  Strain   =  0.15
  Period   =  Strain/Rate
  Disp     =  Xl0*exp(Rate*Period)-Xl0
  Velocity =  Disp/Period
**--------------------------------------**
  Density       = 7.850E-09
  Mass_Scaling  = 100.
**--------------------------------------**
**  Period  =  log((Disp+Xl0)/Xl0)/Rate
**--------------------------------------**
**--------------------------------------**
** Material Parameters for VUMAT BCJ Model (vumat_dmg-55p.f)
**----------  Material Parameters
  G         = 17000.
  a         = 6.160e+05
  K         = 50000.
  b         = 5.86
  Tmelt     = 435.
**----------  BCJ Parameters
  C01       = 0.
  C02       = 0.
  C03       = 92.82
  C04       = 47.93
  C03       = 0.
  C05       = 1.000e-05
  C06       = 5.991e-07
  C07       = 14.8
  C08       = 6.911e-07
  C09       = 4.077e+04
  C10       = 102.4
  C11       = 0.
  C12       = 0.
  C13       = 1.929e+07
  C14       = 6.868e+03
  C15       = 1.577e+03
  C16       = 5.931e-01
  C17       = 6.529e-05
  C18       = 1.064e+06
  C19       = 0.
  C20       = 0.
  Ca        = 1.883
  Cb        = 8.272e-03
**----------  Temperature Parameters
  InitTemp  = 293.
  HeatGen   = 0.
**----------  Microstructure Parameters
  VoidGrthn = 0.01
  InitVRad  = 2.0000e-4
  InitVRad  = 0.
  TorsCsta  = 1.
  NuclCsta  = 1.
  NuclCstc  = 1.
  NuclCoef  = 0.
  KIc       = 17.3
  PartSize  = 4.000e-04
  PartVF    = 0.07
  cd1       = 5.5
  cd2       = 5.
**----------  Grain Size Parameters
  Dcs0      = 20.
  Dcs       = 20.
  zz        = 0.0509
  InitVVF   = 0.01
**----------  Other BCJ Parameters
  C21       = 0.
  C22       = 0.
  C23       = 0.
  C24       = 0.
  C25       = 0.
  C26       = 0.
**----------  Temperature Parameters for Nucleation and Coalescence
  Tnucl     = 0.
  Tcoal     = 0.
**----------  Flag, Cacon and Elastic Modulus porosity exponent
  Flag      = 0.
  Cacon     = 18.
  Zeta      = 0.
**
** END PARAMETERS
**---------------------------------------------------------------**
**---------------------------------------------------------------**
**---------------------------------------------------------------**
*NODE
2, <Xl0>,    0.,    0.
3, <Xl0>, <Xl0>,    0.
4,    0., <Xl0>,    0.
5,    0.,    0., <Xl0>
6, <Xl0>,    0., <Xl0>
7, <Xl0>, <Xl0>, <Xl0>
8,    0., <Xl0>, <Xl0>
*NSET,NSET=FACE1
*NSET,NSET=FACE2
*NSET,NSET=FACE3
*NSET,NSET=FACE4
*NSET,NSET=FACE5
*NSET,NSET=FACE6
*ELEMENT,TYPE=C3D8R,ELSET=EVR
**-------------------------------------------------------------**
**-------------------------------------------------------------**
*SOLID SECTION,ELSET=EVR,MATERIAL=METAL
**HOURGLASS STIFFNESS
**49.248
**-------------------------------------------------------------**
*MATERIAL,NAME=METAL
*DENSITY
 <Density>
*USER MATERIAL,CONSTANTS=55
**---------------------------------------------------------------------------------------------------**
**---------------------------------------------------------------------------------------------------**
        <G>,         <a>,         <K>,         <b>,     <Tmelt>,       <C01>,      <C02>,       <C03>,
      <C04>,       <C05>,       <C06>,       <C07>,       <C08>,       <C09>,      <C10>,       <C11>,
      <C12>,       <C13>,       <C14>,       <C15>,       <C16>,       <C17>,      <C18>,       <C19>,
      <C20>,        <Ca>,        <Cb>,  <InitTemp>,   <HeatGen>, <VoidGrthn>, <InitVRad>,  <TorsCsta>,
 <NuclCsta>,  <NuclCstc>,  <NuclCoef>,       <KIc>,  <PartSize>,    <PartVF>,      <cd1>,       <cd2>,
     <Dcs0>,       <Dcs>,        <zz>,   <InitVVF>,       <C21>,       <C22>,      <C23>,       <C24>,
      <C25>,       <C26>,     <Tnucl>,     <Tcoal>,      <Flag>,     <Cacon>,     <Zeta>
**---------------------------------------------------------------------------------------------------**
*DEPVAR
**---------------------------------------------------------------------------------------------------**
**---------------------------------------------------------------------------------------------------**
**--------------------------------------------------------**
**--------   AMPLITUDE FOR CONSTANT STRAIN RATE   --------**
**--------------------------------------------------------**
**AMPLITUDE,NAME=RATE,VALUE=ABSOLUTE,INPUT=disp.inp
**-------------------------------**
**--    BOUNDARY CONDITIONS    --**
**-------------------------------**
*BOUNDARY
1,2,3
FACE1,3,
FACE3,2,
FACE6,1,
**----------------------------------------**
*STEP,NLGEOM=YES
*DYNAMIC,EXPLICIT
**DYNAMIC,EXPLICIT,DIRECT USER CONTROL
,<Period>
**----------------------------------------**
*BULK VISCOSITY
**----------------------------------------**
** Mass Scaling:
*Fixed Mass Scaling, factor=<Mass_Scaling>
**----------------------------------------**
**----------  TENSION  B.C.  -------------**
**----------------------------------------**
**BOUNDARY,TYPE=VELOCITY
**FACE4,1,1,<Velocity>
**---------------------**
**---  User Boundary Conditions for constant strain rate (vdisp.f)
*BOUNDARY,USER,TYPE=VELOCITY
FACE4,1,1,
**---------------------**
*MONITOR,DOF=1,NODE=2
**----------------------------------------**
**----------------------------------------**
*OUTPUT, FIELD, NUMBER INTERVAL=200
*ELEMENT OUTPUT,ELSET=EVR
S, SDV, LE
*NODE OUTPUT
COORD, U, V
**----------------------------------------**
*FILE OUTPUT,TIMEMARKS=YES,NUM=1000
*NODE FILE
U,RF
*EL FILE
SDV,LE,S
*END STEP
**--------------------------------------------------------------------**
**--------------------------------------------------------------------**
**--------------------------------------------------------------------**