XFEM Failure Analysis Report

Introduction

The FAR for the XFEM module describes the system functional and non-functional requirements that describe the expected interactions that the software shall provide.

Dependencies

The XFEM module is developed using MOOSE and can itself be based on various MOOSE modules, as such the FAR for the XFEM module is dependent upon the following documents.

Framework Failure Analysis Report
sqa/tensor_mechanics_far.md

Failure Analysis Requirements

The following is a complete list for all the requirements related to failure analysis for the XFEM module.

xfem: Mesh Cut 2D Fracture
23.9.1
The XFEM module shall represent a propagating crack in a 2D mechanics model with XFEM, where the crack is defined using a topologically defined surface cutting mesh by MeshCut2DFunctionUserObject
1. and crack growth and direction is computed by a fracture integral with constant fracture toughness.
2. and crack growth and direction is computed by a fracture integral with fracture toughness sampled from a material at the crack tip.
3. and crack growth is a maximum stress criterion and the growth direction is in the original crack direction.
4. and crack growth and direction is computed by a fracture integral with constant fracture toughness or maximum stress criterion.
5. and produce an error when fracture toughness is specified by more than one input source.
6. and crack growth is defined by a function. The average crack tip stress changes with crack orientation and properly accounts volume averaging near surfaces.
Specification(s): 2D/k_growth, 2D/kvpp_growth, 2D/stress_growth, 2D/k_stress_growth, 2D/k_error, 2D/stress_check
Design: MeshCut2DFunctionUserObject
Issue(s): #23572 #27659 #32418
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunExceptionExodiffCSVDiff

xfem: Nucleation Uo
23.12.4The system shall error if the nucleation length is smaller than the element being cut.
Specification(s): errorSmallCut
Design: MeshCut2DRankTwoTensorNucleation
Issue(s): #27054
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException

xfem: Solid Mechanics Basic
23.17.5
The XFEM module shall represent a propagating crack in a 3D mechanics model with XFEM, where the crack is defined using a topological surface cutting mesh by CrackMeshCut3DUserObject and allowed to propagate with
1. growth speeds determined by the fatigue cracking Paris law, or
2. growth speeds determined by the stress corrosion cracking constant timestep, or
3. growth speeds determined by the stress corrosion cracking maximum growth length, or
4. growth directions determined by computing the maximum hoop stress law using interaction integrals, or
5. growth directions determined by an parsed function, or
6. error when the number of points on the cutter mesh do not match the number of points in the DomainIntegralAction.
Specification(s): mesh_cutter_3d/fatigue, mesh_cutter_3d/scc_parseVec, mesh_cutter_3d/scc_reporter, mesh_cutter_3d/mhs, mesh_cutter_3d/function, mesh_cutter_3d/num_points_error
Design: CrackMeshCut3DUserObject ParsedVectorReporter StressCorrosionCrackingExponential
Issue(s): #6320 #23572 #32418
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunExceptionExodiffJSONDiffCSVDiff


[k_growth]
  type = Exodiff
  input = 'edge_crack_2d_propagation.i kcrit_based_meshCut_uo.i'
  exodiff = 'kcrit_based_meshCut_uo_xfemcutter_XFEMCutMeshOutput.e-s0004'
  abs_zero = 1e-8
  map = false
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' and crack growth and direction is computed by a fracture integral with constant fracture toughness.'
[]


[kvpp_growth]
  type = Exodiff
  input = 'edge_crack_2d_propagation.i kvpp_based_meshCut_uo.i'
  exodiff = 'kvpp_based_meshCut_uo_xfemcutter_XFEMCutMeshOutput.e-s0004'
  abs_zero = 1e-8
  map = false
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' and crack growth and direction is computed by a fracture integral with fracture toughness sampled from a material at the crack tip.'
[]


[stress_growth]
  type = Exodiff
  input = 'edge_crack_2d_propagation.i stress_based_meshCut_uo.i'
  exodiff = 'stress_based_meshCut_uo_xfemcutter_XFEMCutMeshOutput.e-s0004'
  abs_zero = 1e-8
  map = false
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' and crack growth is a maximum stress criterion and the growth direction is in the original crack direction.'
[]


[k_stress_growth]
  type = Exodiff
  input = 'edge_crack_2d_propagation.i kcrit_stress_based_meshCut_uo.i'
  exodiff = 'kcrit_stress_based_meshCut_uo_xfemcutter_XFEMCutMeshOutput.e-s0004'
  abs_zero = 1e-8
  map = false
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' and crack growth and direction is computed by a fracture integral with constant fracture toughness or maximum stress criterion.'
[]


[k_error]
  type = 'RunException'
  input = 'edge_crack_2d_propagation.i kvpp_based_meshCut_uo.i'
  cli_args = "UserObjects/cut_mesh2/k_critical=100"
  expect_err = 'Fracture toughness cannot be specified by both k_critical and k_critical_vectorpostprocessor.'
  detail = ' and produce an error when fracture toughness is specified by more than one input source.'
[]


[stress_check]
  type = CSVDiff
  input = crack_front_stress_function_growth.i
  csvdiff = 'stress_check_CrackFrontNonlocalStress_0001.csv stress_check_CrackFrontNonlocalStress_0002.csv stress_check_CrackFrontNonlocalStress_0003.csv stress_check_CrackFrontNonlocalStress_0004.csv'
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' and crack growth is defined by a function.  The average crack tip stress changes with crack orientation and properly accounts volume averaging near surfaces.'
[]


[errorSmallCut]
  issues = '#27054'
  type = RunException
  input = nucleate_bulkCrack.i
  cli_args = 'UserObjects/nucleate/nucleation_length=.001'
  expect_err = 'Trying to nucleate crack smaller than element length'
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  requirement = 'The system shall error if the nucleation length is smaller than the element being cut.'
[]


[fatigue]
  type = Exodiff
  input = 'edge_crack_3d_domain.i edge_crack_3d_fatigue.i'
  exodiff = 'edge_crack_3d_fatigue_out.e edge_crack_3d_fatigue_out.e-s002'
  abs_zero = 1e-8
  map = false
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' growth speeds determined by the fatigue cracking Paris law, or'
[]


[scc_parseVec]
  type = CSVDiff
  cli_args = 'Executioner/end_time=11'
  input = 'edge_crack_3d_domain.i edge_crack_3d_scc.i'
  csvdiff = 'edge_crack_3d_scc_out_II_KI_1_0012.csv'
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' growth speeds determined by the stress corrosion cracking constant timestep, or'
[]


[scc_reporter]
  type = JSONDiff
  cli_args = 'Executioner/end_time=2'
  input = 'edge_crack_3d_domain.i edge_crack_3d_scc_crit.i'
  jsondiff = 'edge_crack_3d_scc_crit_out.json'
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' growth speeds determined by the stress corrosion cracking maximum growth length, or'
[]


[mhs]
  type = Exodiff
  input = 'edge_crack_3d_domain.i edge_crack_3d_mhs.i'
  exodiff = 'edge_crack_3d_mhs_out.e-s003 edge_crack_3d_mhs_out_XFEMCutMeshOutput.e-s0004'
  abs_zero = 1e-4 #crack face is normal to z so this gets rid of exodiffs caused by z-direction variables
  map = false
  partial = true
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' growth directions determined by computing the maximum hoop stress law using interaction integrals, or'
[]


[function]
  type = Exodiff
  input = 'edge_crack_3d_domain.i edge_crack_3d_propagation.i'
  exodiff = 'edge_crack_3d_propagation_out.e edge_crack_3d_propagation_out.e-s002'
  abs_zero = 1e-8
  map = false
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' growth directions determined by an parsed function, or'
[]


[num_points_error]
  type = RunException
  input = 'edge_crack_3d_domain.i edge_crack_3d_mhs.i'
  cli_args = 'DomainIntegral/number_points_from_provider=3'
  expect_err = 'This must match the number of points provided by the XFEM mesh cutter object.'
  # XFEM requires --enable-unique-ids in libmesh
  capabilities = 'unique_id'
  detail = ' error when the number of points on the cutter mesh do not match the number of points in the DomainIntegralAction.'
[]

Introduction
Dependencies
Failure Analysis Requirements