Fault slip test case

inputFiles/lagrangianContactMechanics/ContactMechanics_PassingCrack_base.xml

@@ -143,11 +143,11 @@
     <TableFunction
       name="apertureTable"
       coordinates="{ -1.0e-3, 0.0 }"
-      values="{ 0.0, 0.0 }"/>
+      values="{ 1.0e-3, 1.0e-3 }"/>
   </Functions>
 
   <Outputs>
-    <Silo
+    <VTK
       name="siloOutput"
       plotLevel="3"/>
 

inputFiles/lagrangianContactMechanics/ContactMechanics_slippingFault_base.xml

@@ -0,0 +1,128 @@
+<?xml version="1.0" ?>
+
+<Problem>
+<!-- SPHINX_SOLVER -->
+  <Solvers
+    gravityVector="{0.0, 0.0, 0.0}">
+    <LagrangianContact
+      name="lagrangiancontact"
+      solidSolverName="lagsolve"
+      stabilizationName="TPFAstabilization"
+      logLevel="1"
+      discretization="FE1"
+      targetRegions="{ Region, Fracture }"
+      contactRelationName="fractureMaterial"
+      fractureRegionName="Fracture">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-8"
+        logLevel="2"
+        newtonMaxIter="10"
+        maxNumConfigurationAttempts="10"
+        lineSearchAction="Require"
+        lineSearchMaxCuts="2"
+        maxTimeStepCuts="2"/>
+      <LinearSolverParameters
+        solverType="direct"
+        directParallel="0"
+        logLevel="0"/>
+    </LagrangianContact>
+
+    <SolidMechanics_LagrangianFEM
+      name="lagsolve"
+      timeIntegrationOption="QuasiStatic"
+      logLevel="0"
+      discretization="FE1"
+      targetRegions="{ Region, Fracture }">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-6"
+        newtonMaxIter="5"/>
+      <LinearSolverParameters
+        krylovTol="1.0e-10"
+        logLevel="0"/>
+    </SolidMechanics_LagrangianFEM>
+
+    <SurfaceGenerator
+      name="SurfaceGen"
+      logLevel="0"
+      fractureRegion="Fracture"
+      targetRegions="{ Region }"
+      rockToughness="1.0e6"
+      mpiCommOrder="1"/>
+  </Solvers>
+
+  <Events
+    maxTime="1.0">
+    <SoloEvent
+      name="preFracture"
+      target="/Solvers/SurfaceGen"/>
+
+    <PeriodicEvent
+      name="solverApplications"
+      beginTime="0.0"
+      endTime="1.0"
+      forceDt="0.1"
+      target="/Solvers/lagrangiancontact"/>
+
+    <PeriodicEvent
+      name="outputs"
+      beginTime="0.0"
+      endTime="1.0"
+      forceDt="0.1"
+      target="/Outputs/vtkOutput"/>
+  </Events> 
+
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace
+        name="FE1"
+        order="1"/>
+    </FiniteElements>
+
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="TPFAstabilization"/>
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Region"
+      cellBlocks="{ cb1 }"
+      materialList="{ rock }"/>
+
+    <SurfaceElementRegion
+      name="Fracture"
+      faceBlock="faceElementSubRegion"
+      defaultAperture="0.0"
+      materialList="{ fractureMaterial, rock }"/>
+  </ElementRegions>
+
+<!-- SPHINX_MATERIAL -->
+  <Constitutive>
+    <ElasticIsotropic
+      name="rock"
+      defaultDensity="2700"
+      defaultBulkModulus="2e10"
+      defaultShearModulus="9e9"/>
+
+    <Coulomb
+      name="fractureMaterial"
+      cohesion="0.0"
+      frictionCoefficient="0.1"
+      apertureTableName="apertureTable"/>
+  </Constitutive>
+
+  <Functions>
+    <TableFunction
+      name="apertureTable"
+      coordinates="{ -1.0e-3, 0.0 }"
+      values="{ 1.0e-6, 9.15e-6 }"/>   
+
+    <TableFunction
+      name="dispFunction"
+      inputVarNames="{ time }"
+      coordinates="{ 0.0, 1.0 }"
+      values="{ 0.0, 1.0 }"/>
+  </Functions>
+
+</Problem>

inputFiles/lagrangianContactMechanics/ContactMechanics_slippingFault_horizontal_smoke.xml

@@ -0,0 +1,125 @@
+<?xml version="1.0" ?>
+
+<Problem>
+
+  <Included>
+    <File name="./ContactMechanics_slippingFault_base.xml"/>
+  </Included>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ -100, 100 }"
+      yCoords="{ -25, 0, 25 }"
+      zCoords="{ 0, 1 }"
+      nx="{ 40 }"
+      ny="{ 5, 5 }"
+      nz="{ 1 }"
+      cellBlockNames="{ cb1 }"/>
+  </Mesh>
+
+  <Geometry>
+    <BoundedPlane
+      name="fracture"
+      normal="{0.0, 1.0, 0.0}"
+      origin="{0.0, 0.0, 0.0}"
+      lengthVector="{1.0, 0.0, 0.0}"
+      widthVector="{0.0, 0.0, 1.0}"
+      dimensions="{ 180, 10 }"/>
+  </Geometry>
+
+  <FieldSpecifications>
+    <FieldSpecification
+      name="frac"
+      initialCondition="1"
+      setNames="{ fracture }"
+      objectPath="faceManager"
+      fieldName="ruptureState"
+      scale="1"/>
+
+    <FieldSpecification
+      name="separableFace"
+      initialCondition="1"
+      setNames="{ core }"
+      objectPath="faceManager"
+      fieldName="isFaceSeparable"
+      scale="1"/>
+
+    <FieldSpecification
+      name="xconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="0"
+      scale="0.0"
+      setNames="{ yneg }"/>
+
+    <FieldSpecification
+      name="yconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="1"
+      scale="0.0"
+      setNames="{ yneg, ypos, xneg, xpos }"/>
+
+    <FieldSpecification
+      name="zconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="2"
+      scale="0.0"
+      setNames="{ zneg, zpos }"/>
+
+    <FieldSpecification 
+      name="Sigmay"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Region"
+      fieldName="rock_stress"
+      component="1"
+      scale="-4.25e6"/>
+
+    <!-- <Traction
+      name="normalLoad"
+      objectPath="faceManager"
+      tractionType="vector"
+      direction="{0,1,0}"
+      scale="-4.25e6"
+      setNames="{ ypos }"/>  -->
+
+    <!-- <Traction
+      name="shearLoad"
+      objectPath="faceManager"
+      tractionType="normal"
+      scale="2e6"
+      setNames="{ ypos }"
+      functionName="tractionFunction"/>  -->
+
+    <!-- <FieldSpecification
+      name="normalLoad"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="1"
+      scale="-0.00821"
+      setNames="{ ypos }"/> -->
+
+    <FieldSpecification
+      name="shearLoad"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="0"
+      scale="0.01"
+      setNames="{ ypos }"
+      functionName="dispFunction"/>
+  </FieldSpecifications>
+<!-- SPHINX_BC_END -->
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      plotFileRoot="faultSlip_horizontal"/>
+
+    <Restart
+      name="restartOutput"/>      
+    </Outputs>
+</Problem>

inputFiles/lagrangianContactMechanics/ContactMechanics_slippingFault_vertical_smoke.xml

@@ -0,0 +1,100 @@
+<?xml version="1.0" ?>
+
+<Problem>
+
+  <Included>
+    <File name="./ContactMechanics_slippingFault_base.xml"/>
+  </Included>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ -25, 0, 25 }"
+      yCoords="{ -100, 100 }"
+      zCoords="{ 0, 1 }"
+      nx="{ 5, 5 }"
+      ny="{ 40 }"
+      nz="{ 1 }"
+      cellBlockNames="{ cb1 }"/>
+  </Mesh>
+
+  <Geometry>
+    <BoundedPlane
+      name="fracture"
+      normal="{1.0, 0.0, 0.0}"
+      origin="{0.0, 0.0, 0.0}"
+      lengthVector="{0.0, 1.0, 0.0}"
+      widthVector="{0.0, 0.0, 1.0}"
+      dimensions="{ 180, 10 }"/>
+  </Geometry>
+
+  <FieldSpecifications>
+    <FieldSpecification
+      name="frac"
+      initialCondition="1"
+      setNames="{ fracture }"
+      objectPath="faceManager"
+      fieldName="ruptureState"
+      scale="1"/>
+
+    <FieldSpecification
+      name="separableFace"
+      initialCondition="1"
+      setNames="{ core }"
+      objectPath="faceManager"
+      fieldName="isFaceSeparable"
+      scale="1"/>
+
+    <FieldSpecification
+      name="yconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="1"
+      scale="0.0"
+      setNames="{ xneg }"/>
+
+    <FieldSpecification
+      name="xconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="0"
+      scale="0.0"
+      setNames="{ yneg, ypos, xneg, xpos }"/>
+
+    <FieldSpecification
+      name="zconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="2"
+      scale="0.0"
+      setNames="{ zneg, zpos }"/>
+
+    <FieldSpecification 
+      name="Sigmax"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Region"
+      fieldName="rock_stress"
+      component="0"
+      scale="-4.25e6"/>
+
+    <FieldSpecification
+      name="shearLoad"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="1"
+      scale="0.01"
+      setNames="{ xpos }"
+      functionName="dispFunction"/>
+  </FieldSpecifications>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      plotFileRoot="faultSlip_vertical"/>
+
+    <Restart
+      name="restartOutput"/>      
+    </Outputs>
+</Problem>