MeshCut2DRankTwoTensorNucleation

Nucleate a crack in MeshCut2D UO based on a scalar extracted from a RankTwoTensor

Overview

The XFEM nucleation UserObject MeshCut2DRankTwoTensorNucleation defines new cracks that are added to the cutter mesh used by MeshCut2DFractureUserObject when a user-defined nucleation criterion is met. The nucleation criterion is based on the "scalar_type" extracted from a RankTwoTensor specified by "tensor" (used to store stresses and strains in SolidMechanics), such as a principal stress or a component of stress.

If the scalar exceeds a threshold specified by "nucleation_threshold", a line segment with a length specified by "nucleation_length" will be inserted into the MeshCut2DFractureUserObject cutter mesh. Cracks are only allowed to initiate from elements on boundaries specified by "initiate_on_boundary". Once the nucleation criterion is reached, a line segment of the specified length is inserted into the cutter mesh, centered on the element centroid it nucleates from. The direction of the nucleated crack is normal to the direction returned by the RankTwoTensor scalar. For example, MaxInPlanePrincipal returns the direction of the maximum in-plane principal component and the crack direction is normal to this. The nucleation length should be at least the length of the element it nucleates in so that the nucleated crack will completely cut the element. A crack will only be nucleated if it is at least a distance specified by "nucleation_radius" away from existing or nucleated cracks. If multiple cracks nucleate in the same xfem update and are within the specified nucleation_radius, the crack nucleated from the element with the lowest id will be retained and no other cracks within the nucleation_radius will be nucleated.

MeshCut2DRankTwoTensorNucleation copies several features available in the XFEMRankTwoTensorMarkerUserObject. These include the nucleation threshold being provided as a coupled variable and the computation of the maximum value of the scalar quantity used for nucleation. By providing the nucleation threshold as a coupled variable, it can be specified as either a constant or variable value. Coupled variable input is useful for introducing randomness in the strength by using an AuxVariable that has been initialized with a random initial condition. For example, the Volume Weighted Weibull initial condition is well-suited for this purpose. The determination of crack nucleation is based on the average value of the scalar quantity over all the quadrature points in an element.

Example Input File Syntax

[UserObjects]
  [nucleate]
    type = MeshCut2DRankTwoTensorNucleation
    tensor = stress
    scalar_type = MaxPrincipal
    nucleation_threshold = 180
    initiate_on_boundary = 'left'
    nucleation_length = .2
    nucleation_radius = .21
  []
[]

Input Parameters

initiate_on_boundaryCracks can only initiate on elements adjacent to specified boundaries.
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:Cracks can only initiate on elements adjacent to specified boundaries.
nucleation_lengthNucleated crack length.
C++ Type:double
Controllable:No
Description:Nucleated crack length.
nucleation_thresholdThreshold for the scalar quantity of the RankTwoTensor to nucleate new cracks
C++ Type:std::vector<VariableName>
Controllable:No
Description:Threshold for the scalar quantity of the RankTwoTensor to nucleate new cracks
tensorThe material tensor name.
C++ Type:std::string
Controllable:No
Description:The material tensor name.

Required Parameters

blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
execute_onXFEM_MARKThe list of flag(s) indicating when this object should be executed, the available options include XFEM_MARK, FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
Default:XFEM_MARK
C++ Type:ExecFlagEnum
Options:XFEM_MARK, FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM
Controllable:No
Description:The list of flag(s) indicating when this object should be executed, the available options include XFEM_MARK, FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
nucleation_radius0Cracks will only nucleate if they are outside the nucleation_radius of an existing crack.
Default:0
C++ Type:double
Controllable:No
Description:Cracks will only nucleate if they are outside the nucleation_radius of an existing crack.
point10 0 0Start point for axis used to calculate some cylindrical material tensor quantities
Default:0 0 0
C++ Type:libMesh::Point
Controllable:No
Description:Start point for axis used to calculate some cylindrical material tensor quantities
point20 1 0End point for axis used to calculate some cylindrical material tensor quantities
Default:0 1 0
C++ Type:libMesh::Point
Controllable:No
Description:End point for axis used to calculate some cylindrical material tensor quantities
prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
scalar_typeScalar quantity to be computed from tensor and used as a failure criterion
C++ Type:MooseEnum
Options:VonMisesStress, EffectiveStrain, Hydrostatic, L2norm, MaxPrincipal, MidPrincipal, MinPrincipal, VolumetricStrain, FirstInvariant, SecondInvariant, ThirdInvariant, AxialStress, HoopStress, RadialStress, TriaxialityStress, Direction, MaxShear, StressIntensity
Controllable:No
Description:Scalar quantity to be computed from tensor and used as a failure criterion
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Optional Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.
execution_order_group-1Nucleation UO needs to be completely executed before GeometricCutUserObject.
Default:-1
C++ Type:int
Controllable:No
Description:Nucleation UO needs to be completely executed before GeometricCutUserObject.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Controllable:No
Description:Determines whether this object is calculated using an implicit or explicit form
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Controllable:No
Description:The seed for the master random number generator
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

Input Files

(modules/xfem/test/tests/nucleation_uo/nucleate_edge_crack_2d.i)
(modules/xfem/test/tests/nucleation_uo/nucleate_AllEdgeCracks.i)
(modules/xfem/test/tests/nucleation_uo/nucleate_2edge_cracks_2d.i)

(modules/xfem/test/tests/nucleation_uo/nucleate_AllEdgeCracks.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[XFEM]
  geometric_cut_userobjects = 'cut_mesh2'
  qrule = volfrac
  output_cut_plane = true
[]

[Mesh]
[gen]
  type = GeneratedMeshGenerator
  dim = 2
  nx = 10
  ny = 20
  xmin = 0
  xmax = 1.0
  ymin = 0.0
  ymax = 2.0
  elem_type = QUAD4
[]
[top_left]
  type = BoundingBoxNodeSetGenerator
  new_boundary = pull_top_left
  bottom_left = '-0.01 1.99 0'
  top_right = '0.11 2.01 0'
  input = gen
[]
[]

[DomainIntegral]
  integrals = 'InteractionIntegralKI InteractionIntegralKII'
  displacements = 'disp_x disp_y'
  crack_front_points_provider = cut_mesh2
  2d=true
  number_points_from_provider = 0
  crack_direction_method = CurvedCrackFront
  radius_inner = '0.15'
  radius_outer = '0.45'
  poissons_ratio = 0.3
  youngs_modulus = 207000
  block = 0
  incremental = true
  used_by_xfem_to_grow_crack = true
[]

[UserObjects]
  [nucleate]
    type = MeshCut2DRankTwoTensorNucleation
    tensor = stress
    scalar_type = MaxPrincipal
    nucleation_threshold = 180
    initiate_on_boundary = 'left'
    nucleation_length = .2
    nucleation_radius = .21
  []
  [cut_mesh2]
    type = MeshCut2DFractureUserObject
    mesh_file = make_edge_crack_in.e
    k_critical=500000 #Large so that cracks will not grow
    growth_increment = 0.11
    nucleate_uo = nucleate
  []
[]

[Modules/TensorMechanics/Master]
  [all]
    strain = FINITE
    planar_formulation = plane_strain
    add_variables = true
    generate_output = 'stress_xx stress_yy vonmises_stress max_principal_stress'
  []
[]

[Functions]
  [bc_pull_top]
    type = ParsedFunction
    expression = 0.0005*t
  []
[]

[BCs]
  [top_edges]
      type = FunctionDirichletBC
      boundary = 'pull_top_left'
      variable = disp_y
      function = bc_pull_top
  []
  [bottom_x]
    type = DirichletBC
    boundary = bottom
    variable = disp_x
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    boundary = bottom
    variable = disp_y
    value = 0.0
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 207000
    poissons_ratio = 0.3
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
[]

[Executioner]
  type = Transient

  solve_type = 'PJFNK'
  petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type -pc_hypre_boomeramg_max_iter'
  petsc_options_value = '201                hypre    boomeramg      8'

  line_search = 'none'

  [./Predictor]
    type = SimplePredictor
    scale = 1.0
  [../]

  l_max_its = 100
  l_tol = 1e-2

  nl_max_its = 15
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-9

  start_time = 0.0
  dt = 1.0
  end_time = 5
  max_xfem_update = 100
[]

[Outputs]
  csv=true
  exodus=true
  execute_on = TIMESTEP_END
[]

(modules/xfem/test/tests/nucleation_uo/nucleate_edge_crack_2d.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
  volumetric_locking_correction = true
[]

[XFEM]
  geometric_cut_userobjects = 'cut_mesh2'
  qrule = volfrac
  output_cut_plane = true
[]

[Mesh]
[gen]
  type = GeneratedMeshGenerator
  dim = 2
  nx = 20
  ny = 10
  xmin = 0
  xmax = 2
  ymin = 0.0
  ymax = 1.0
  elem_type = QUAD4
[]
[top_left]
  type = BoundingBoxNodeSetGenerator
  new_boundary = pull_top_left
  bottom_left = '-0.01 0.99 0'
  top_right = '0.11 1.01 0'
  input = gen
[]
[top_right]
  type = BoundingBoxNodeSetGenerator
  new_boundary = pull_top_right
  bottom_left = '1.89 0.99 0'
  top_right = '2.01 1.01 0'
  input = top_left
[]
[top_middle_ss]
  type = SideSetsFromBoundingBoxGenerator
  input = top_right
  bottom_left = '0.79 0.89 0'
  top_right = '1.21 1.01 0'
  block_id = '0'
  boundary_new = top_middle_ss
  boundaries_old = top
[]
[nucleate]
  type = ParsedSubdomainMeshGenerator
  input = top_middle_ss
  combinatorial_geometry = 'y > 0.39 & y < 0.51'
  block_id = 10
[]
[]

[DomainIntegral]
  integrals = 'InteractionIntegralKI InteractionIntegralKII'
  displacements = 'disp_x disp_y'
  crack_front_points_provider = cut_mesh2
  2d=true
  number_points_from_provider = 0
  crack_direction_method = CurvedCrackFront
  radius_inner = '0.15'
  radius_outer = '0.45'
  poissons_ratio = 0.3
  youngs_modulus = 207000
  block = 0
  incremental = true
  used_by_xfem_to_grow_crack = true
[]

[UserObjects]
  [nucleate]
    type = MeshCut2DRankTwoTensorNucleation
    tensor = stress
    scalar_type = MaxPrincipal
    nucleation_threshold = nucleation_threshold
    initiate_on_boundary = 'left right'
    nucleation_length = .2
  []
  [cut_mesh2]
    type = MeshCut2DFractureUserObject
    mesh_file = make_edge_crack_in.e
    k_critical=230
    growth_increment = 0.11
    nucleate_uo = nucleate
  []
[]
[AuxVariables]
  [nucleation_threshold]
    order = CONSTANT
    family = MONOMIAL
  []
[]
[ICs]
   [nucleation_bulk]
     type = ConstantIC
     value = 10000
     variable = nucleation_threshold
     block = 0
   []
   [nucleation_weak]
    type = FunctionIC
    function = nucleation_x
    variable = nucleation_threshold
    block = 10
  []
[]

[Functions]
  [nucleation_x]
    type = ParsedFunction
    expression = '300+x*50'
  []
[]

[Modules/TensorMechanics/Master]
  [./all]
    strain = FINITE
    planar_formulation = plane_strain
    add_variables = true
    generate_output = 'stress_xx stress_yy vonmises_stress max_principal_stress'
  [../]
[]

[Functions]
  [bc_pull_top]
    type = ParsedFunction
    expression = 0.0005*t
  []
[]

[BCs]
  [top_edges]
      type = FunctionDirichletBC
      boundary = 'pull_top_left pull_top_right'
      variable = disp_y
      function = bc_pull_top
  []
  [top_middle]
    type = NeumannBC
    boundary = top_middle_ss
    variable = disp_y
    value = -2000
  []
  [bottom_x]
    type = DirichletBC
    boundary = bottom
    variable = disp_x
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    boundary = bottom
    variable = disp_y
    value = 0.0
  []
[]

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 207000
    poissons_ratio = 0.3
  [../]
  [./stress]
    type = ComputeFiniteStrainElasticStress
  [../]
[]

[Executioner]
  type = Transient

  solve_type = 'PJFNK'
  petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type -pc_hypre_boomeramg_max_iter'
  petsc_options_value = '201                hypre    boomeramg      8'

  line_search = 'none'

  [./Predictor]
    type = SimplePredictor
    scale = 1.0
  [../]

# controls for linear iterations
  l_max_its = 100
  l_tol = 1e-2

# controls for nonlinear iterations
  nl_max_its = 15
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-9

# time control
  start_time = 0.0
  dt = 1.0
  end_time = 5
  max_xfem_update = 100
[]

[Outputs]
  csv=true
  execute_on = TIMESTEP_END
  # [xfemcutter]
  #   type=XFEMCutMeshOutput
  #   xfem_cutter_uo=cut_mesh2
  # []
  # console = false
  [./console]
    type = Console
    output_linear = false
    output_nonlinear = false
  [../]
[]

(modules/xfem/test/tests/nucleation_uo/nucleate_AllEdgeCracks.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[XFEM]
  geometric_cut_userobjects = 'cut_mesh2'
  qrule = volfrac
  output_cut_plane = true
[]

[Mesh]
[gen]
  type = GeneratedMeshGenerator
  dim = 2
  nx = 10
  ny = 20
  xmin = 0
  xmax = 1.0
  ymin = 0.0
  ymax = 2.0
  elem_type = QUAD4
[]
[top_left]
  type = BoundingBoxNodeSetGenerator
  new_boundary = pull_top_left
  bottom_left = '-0.01 1.99 0'
  top_right = '0.11 2.01 0'
  input = gen
[]
[]

[DomainIntegral]
  integrals = 'InteractionIntegralKI InteractionIntegralKII'
  displacements = 'disp_x disp_y'
  crack_front_points_provider = cut_mesh2
  2d=true
  number_points_from_provider = 0
  crack_direction_method = CurvedCrackFront
  radius_inner = '0.15'
  radius_outer = '0.45'
  poissons_ratio = 0.3
  youngs_modulus = 207000
  block = 0
  incremental = true
  used_by_xfem_to_grow_crack = true
[]

[UserObjects]
  [nucleate]
    type = MeshCut2DRankTwoTensorNucleation
    tensor = stress
    scalar_type = MaxPrincipal
    nucleation_threshold = 180
    initiate_on_boundary = 'left'
    nucleation_length = .2
    nucleation_radius = .21
  []
  [cut_mesh2]
    type = MeshCut2DFractureUserObject
    mesh_file = make_edge_crack_in.e
    k_critical=500000 #Large so that cracks will not grow
    growth_increment = 0.11
    nucleate_uo = nucleate
  []
[]

[Modules/TensorMechanics/Master]
  [all]
    strain = FINITE
    planar_formulation = plane_strain
    add_variables = true
    generate_output = 'stress_xx stress_yy vonmises_stress max_principal_stress'
  []
[]

[Functions]
  [bc_pull_top]
    type = ParsedFunction
    expression = 0.0005*t
  []
[]

[BCs]
  [top_edges]
      type = FunctionDirichletBC
      boundary = 'pull_top_left'
      variable = disp_y
      function = bc_pull_top
  []
  [bottom_x]
    type = DirichletBC
    boundary = bottom
    variable = disp_x
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    boundary = bottom
    variable = disp_y
    value = 0.0
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 207000
    poissons_ratio = 0.3
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
[]

[Executioner]
  type = Transient

  solve_type = 'PJFNK'
  petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type -pc_hypre_boomeramg_max_iter'
  petsc_options_value = '201                hypre    boomeramg      8'

  line_search = 'none'

  [./Predictor]
    type = SimplePredictor
    scale = 1.0
  [../]

  l_max_its = 100
  l_tol = 1e-2

  nl_max_its = 15
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-9

  start_time = 0.0
  dt = 1.0
  end_time = 5
  max_xfem_update = 100
[]

[Outputs]
  csv=true
  exodus=true
  execute_on = TIMESTEP_END
[]

(modules/xfem/test/tests/nucleation_uo/nucleate_2edge_cracks_2d.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
  volumetric_locking_correction = true
[]

[XFEM]
  geometric_cut_userobjects = 'cut_mesh2'
  qrule = volfrac
  output_cut_plane = true
[]

[Mesh]
[gen]
  type = GeneratedMeshGenerator
  dim = 2
  nx = 30
  ny = 15
  xmin = -2
  xmax = -.2
  ymin = 0.0
  ymax = 1.0
  elem_type = QUAD4
[]
[dispBlock_top]
  type = BoundingBoxNodeSetGenerator
  new_boundary = pull_top_y
  bottom_left = '-2.1 0.99 0'
  top_right = '-1.9 1.01 0'
  input = gen
[]
[dispBlock_bot]
  type = BoundingBoxNodeSetGenerator
  new_boundary = pull_bot_y
  bottom_left = '-2.1 -.01 0'
  top_right = '-1.9 0.01 0'
  input = dispBlock_top
[]
[]

[DomainIntegral]
  integrals = 'Jintegral InteractionIntegralKI InteractionIntegralKII'
  displacements = 'disp_x disp_y'
  crack_front_points_provider = cut_mesh2
  2d=true
  number_points_from_provider = 1
  crack_direction_method = CurvedCrackFront
  radius_inner = '0.15'
  radius_outer = '0.45'
  poissons_ratio = 0.3
  youngs_modulus = 207000
  block = 0
  incremental = true
  used_by_xfem_to_grow_crack = true
[]

[UserObjects]
  #fixme, nucleate has to be before cut_mesh2 in the input file or cut_mesh2 can't finde the nucleate_uo
  [nucleate]
    type = MeshCut2DRankTwoTensorNucleation
    tensor = stress
    scalar_type = MaxPrincipal
    nucleation_threshold = nucleation_threshold
    initiate_on_boundary = 'left bottom'
    average = true
    nucleation_length = .1
  []
  [cut_mesh2]
    type = MeshCut2DFractureUserObject
    mesh_file = make_edge_crack_in.e
    k_critical=80
    growth_increment = 0.1
    nucleate_uo = nucleate
  []
[]
[AuxVariables]
  [nucleation_threshold]
    order = CONSTANT
    family = MONOMIAL
  []
[]
[ICs]
   [nucleation]
     type = FunctionIC
     function =  nucleation_x_y
     variable = nucleation_threshold
   []
  # [nucleation]
  #   type = VolumeWeightedWeibull
  #   variable = nucleation_threshold
  #   reference_volume = 0.01 #This is the volume of an element for a 100x100 mesh
  #   weibull_modulus = 2
  #   median = 5000.0
  # []
[]

[Functions]
  [nucleation_y]
    type = ParsedFunction
    expression = 'if(y>0.7,10000,if(y<0.5,10000,4000*(1-y)^2-10000))'
  []
  [nucleation_x]
    type = ParsedFunction
    expression = 'if(x>-0.9,10000,if(x<-1.1,10000,1000*(x)^2-10000))'
  []
  [nucleation_x_y]
    type = LinearCombinationFunction
    functions = 'nucleation_x nucleation_y'
    w = '1 1'
  []
[]

[Modules/TensorMechanics/Master]
  [./all]
    strain = FINITE
    planar_formulation = plane_strain
    add_variables = true
    generate_output = 'stress_xx stress_yy vonmises_stress max_principal_stress'
  [../]
[]

[Functions]
  [bc_pull_top]
    type = ParsedFunction
    expression = 'if(t<6,0.0008*t,0.0008*5+0.0004*(t-6))'
  []
  [bc_pull_bot]
    type = ParsedFunction
    expression = 0.0004*t
  []
[]

[BCs]
  [top_left]
      type = FunctionDirichletBC
      boundary = pull_top_y
      variable = disp_y
      function = bc_pull_top
  []
  [bot_left]
    type = FunctionDirichletBC
    boundary = pull_bot_y
    variable = disp_y
    function = bc_pull_bot
[]
  [bottom_x]
    type = DirichletBC
    boundary = right
    variable = disp_x
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    boundary = right
    variable = disp_y
    value = 0.0
  []
[]

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 207000
    poissons_ratio = 0.3
    block = 0
  [../]
  [./stress]
    type = ComputeFiniteStrainElasticStress
    block = 0
  [../]
[]

[Executioner]
  type = Transient

  solve_type = 'PJFNK'
  petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type -pc_hypre_boomeramg_max_iter'
  petsc_options_value = '201                hypre    boomeramg      8'

  line_search = 'none'

  [./Predictor]
    type = SimplePredictor
    scale = 1.0
  [../]

# controls for linear iterations
  l_max_its = 100
  l_tol = 1e-2

# controls for nonlinear iterations
  nl_max_its = 15
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-9

# time control
  start_time = 0.0
  dt = 1.0
  end_time = 55
  max_xfem_update = 2
[]

[Outputs]
  # csv=true
  exodus = true
  execute_on = TIMESTEP_END
  # [xfemcutter]
  #   type=XFEMCutMeshOutput
  #   xfem_cutter_uo=cut_mesh2
  # []
  # console = false
  [./console]
    type = Console
    output_linear = false
    output_nonlinear = false
  [../]
[]

Overview
Example Input File Syntax
Input Parameters
Input Files