elasticHollowCylinder_base.xml

<?xml version="1.0" ?>

<Problem>
  <!-- SPHINX_Gravity -->
  <Solvers
    gravityVector="{ 0.0, 0.0, 0.0 }">
    <!-- SPHINX_GravityEnd -->
    <!-- SPHINX_SolidMechanicsSolver -->
    <SolidMechanicsLagrangianSSLE
      name="lagsolve"
      timeIntegrationOption="QuasiStatic"
      discretization="FE1"
      targetRegions="{ Region1 }">
      <!-- SPHINX_SolidMechanicsSolverEnd -->
      <NonlinearSolverParameters
        newtonTol="1.0e-6"
        newtonMaxIter="8"/>
      <LinearSolverParameters
        directParallel="0"/>
    </SolidMechanicsLagrangianSSLE>
  </Solvers>

  <!-- SPHINX_WellboreMesh -->
  <Mesh>
    <InternalWellbore
      name="mesh1"
      elementTypes="{ C3D8 }"
      radius="{ 0.1, 0.11 }"
      theta="{ 0, 360 }"
      zCoords="{ 0, 1 }"
      nr="{ 4 }"
      nt="{ 16 }"
      nz="{ 1 }"
      trajectory="{ { 0.0, 0.0, 0.0 },
                    { 0.0, 0.0, 1.0 } }"
      cellBlockNames="{ cb1 }"/>
  </Mesh>

  <Geometry>
    <Box
      name="xConstraint"
      xMin="{ -0.01, -1e9, -1e9 }"
      xMax="{ 0.01, 1e9, 1e9 }"/>

    <Box
      name="yConstraint"
      xMin="{ -1e9, -0.01, -1e9 }"
      xMax="{ 1e9, 0.01, 1e9 }"/>
  </Geometry>

  <!-- SPHINX_WellboreMeshEnd -->
  <Events
    maxTime="5.0">
    <!-- This event is applied every cycle, and overrides the
    solver time-step request -->
    <!-- SPHINX_Execution -->
    <PeriodicEvent
      name="solverApplications"
      forceDt="1.0"
      target="/Solvers/lagsolve"/>

    <!-- SPHINX_ExecutionEnd -->
    <!-- This event is applied every 5.0e-5s.  The targetExactTimestep
    flag allows this event to request a dt modification to match an
    integer multiple of the timeFrequency. -->
    <PeriodicEvent
      name="outputs"
      timeFrequency="1.0"
      target="/Outputs/vtkOutput"/>

    <PeriodicEvent
      name="restarts"
      timeFrequency="2.0"
      targetExactTimestep="0"
      target="/Outputs/restartOutput"/>
  </Events>

  <!-- SPHINX_FEDiscretization -->
  <NumericalMethods>
    <FiniteElements>
      <FiniteElementSpace
        name="FE1"
        order="1"/>
    </FiniteElements>
  </NumericalMethods>

  <!-- SPHINX_FEDiscretizationEnd -->

  <Constitutive>
    <!-- SPHINX_Material -->
    <ElasticIsotropic
      name="isotropicMaterial"
      defaultDensity="2700"
      defaultBulkModulus="175e9"
      defaultShearModulus="80.8e9"/>

    <!-- SPHINX_MaterialEnd -->

    <ElasticTransverseIsotropic
      name="transverseIsotropicMaterial"
      defaultDensity="2700"
      defaultPoissonRatioAxialTransverse="0.3"
      defaultPoissonRatioTransverse="0.3"
      defaultYoungModulusAxial="210.0e9"
      defaultYoungModulusTransverse="210.0e9"
      defaultShearModulusAxialTransverse="80.8e9"/>

    <ElasticOrthotropic
      name="orthotropicMaterial"
      defaultDensity="2700"
      defaultC11="283e9"
      defaultC12="121e9"
      defaultC13="121e9"
      defaultC22="283e9"
      defaultC23="121e9"
      defaultC33="283e9"
      defaultC44="80.8e9"
      defaultC55="80.8e9"
      defaultC66="80.8e9"/>

  </Constitutive>

  <!-- SPHINX_BoundaryConditions -->
  <FieldSpecifications>
    <FieldSpecification
      name="xnegConstraint"
      objectPath="nodeManager"
      fieldName="totalDisplacement"
      component="0"
      scale="0.0"
      setNames="{ xConstraint }"/>

    <FieldSpecification
      name="ynegConstraint"
      objectPath="nodeManager"
      fieldName="totalDisplacement"
      component="1"
      scale="0.0"
      setNames="{ yConstraint }"/>

    <FieldSpecification
      name="zconstraint"
      objectPath="nodeManager"
      fieldName="totalDisplacement"
      component="2"
      scale="0.0"
      setNames="{ zneg, zpos }"/>

    <Traction
      name="insituStress"
      objectPath="faceManager"
      tractionType="stress"
      scale="1.0e6"
      inputStress="{ -1, -1.5, -2, 0, 0, 0 }"
      setNames="{ rpos }"/>

    <Traction
      name="innerPressure"
      objectPath="faceManager"
      tractionType="normal"
      scale="-1.0e6"
      functionName="timeFunction"
      setNames="{ rneg }"/>
  </FieldSpecifications>

  <!-- SPHINX_BoundaryConditionsEnd -->
  <!-- SPHINX_TableFunction -->
  <Functions>
    <TableFunction
      name="timeFunction"
      inputVarNames="{ time }"
      coordinates="{ 0.0, 1.0, 10.0 }"
      values="{ 0.0, 0.0, 10.0 }"/>
  </Functions>

  <!-- SPHINX_TableFunctionEnd -->
  <Outputs>
    <VTK 
      name="vtkOutput"    
    />

    <Restart
      name="restartOutput"/>
  </Outputs>
</Problem>