CrackMeshCut3DUserObject

Creates a UserObject for a mesh cutter in 3D problems

Overview

This class: (1) reads in a mesh describing the crack surface, (2) uses the mesh to do initial cutting of 3D elements, and (3) grows the mesh incrementally based on prescribed growth functions. The code is interfaced with domain integral methods to allow nonplanar crack growth based on empirical propagation direction and speed laws.

Example Input Syntax

This example shows the Mesh block in Listing 1 needed for creating the cutter mesh along with the CrackMeshCut3DUserObject block in Listing 2. The mesh block in Listing 1 defines two separate meshes. The cutter mesh is created in the read_in_cutter_mesh block and must have "save_with_name" set in order to specify this mesh in CrackMeshCut3DUserObject shown in Listing 2 using "mesh_generator_name". The mesh used by the FEM simulation is specifed in the FEM_mesh block in this example and "final_generator"=FEM_mesh must be set because only the FEM_mesh will be used for the finite-element solution and the mesh created by read_in_cutter_mesh will be ignored by the solution.

Listing 1: Setting up the mesh block contain simulation and cutter meshes.

[Mesh<<<{"href": "../../syntax/Mesh/index.html"}>>>]
  [read_in_cutter_mesh]
    type = FileMeshGenerator<<<{"description": "Read a mesh from a file.", "href": "../meshgenerators/FileMeshGenerator.html"}>>>
    file<<<{"description": "The filename to read."}>>> = mesh_edge_crack.xda
    save_with_name<<<{"description": "Keep the mesh from this mesh generator in memory with the name specified"}>>> = mesh_cutter
  []
  [FEM_mesh]
    type = GeneratedMeshGenerator<<<{"description": "Create a line, square, or cube mesh with uniformly spaced or biased elements.", "href": "../meshgenerators/GeneratedMeshGenerator.html"}>>>
    dim<<<{"description": "The dimension of the mesh to be generated"}>>> = 3
    nx<<<{"description": "Number of elements in the X direction"}>>> = 5
    ny<<<{"description": "Number of elements in the Y direction"}>>> = 5
    nz<<<{"description": "Number of elements in the Z direction"}>>> = 2
    xmin<<<{"description": "Lower X Coordinate of the generated mesh"}>>> = 0.0
    xmax<<<{"description": "Upper X Coordinate of the generated mesh"}>>> = 1.0
    ymin<<<{"description": "Lower Y Coordinate of the generated mesh"}>>> = 0.0
    ymax<<<{"description": "Upper Y Coordinate of the generated mesh"}>>> = 1.0
    zmin<<<{"description": "Lower Z Coordinate of the generated mesh"}>>> = 0.0
    zmax<<<{"description": "Upper Z Coordinate of the generated mesh"}>>> = 0.2
    elem_type<<<{"description": "The type of element from libMesh to generate (default: linear element for requested dimension)"}>>> = HEX8
  []
  final_generator = FEM_mesh
[]
(modules/xfem/test/tests/solid_mechanics_basic/edge_crack_3d_mhs.i)

Listing 2: CrackMeshCut3DUserObject that uses the cutter mesh created in Listing 2.

[UserObjects<<<{"href": "../../syntax/UserObjects/index.html"}>>>]
  [./cut_mesh]
    type = CrackMeshCut3DUserObject<<<{"description": "Creates a UserObject for a mesh cutter in 3D problems", "href": "CrackMeshCut3DUserObject.html"}>>>
    mesh_generator_name<<<{"description": "Mesh generator for the XFEM geometric cutter."}>>> = mesh_cutter
    growth_dir_method<<<{"description": "choose from FUNCTION, MAX_HOOP_STRESS"}>>> = MAX_HOOP_STRESS
    size_control<<<{"description": "Criterion for refining elements while growing the crack"}>>> = 1
    n_step_growth<<<{"description": "Number of steps for crack growth"}>>> = 1
    growth_rate<<<{"description": "Function defining crack growth rate"}>>> = growth_func_v
    crack_front_nodes<<<{"description": "Set of nodes to define crack front"}>>> = '7 6 5 4'
    execution_order_group<<<{"description": "Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group."}>>> = -2
  [../]
[]
(modules/xfem/test/tests/solid_mechanics_basic/edge_crack_3d_mhs.i)

Input 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

  • crack_front_nodesSet of nodes to define crack front

    C++ Type:std::vector<unsigned long>

    Controllable:No

    Description:Set of nodes to define crack front

  • growth_dir_methodFUNCTIONchoose from FUNCTION, MAX_HOOP_STRESS

    Default:FUNCTION

    C++ Type:MooseEnum

    Options:MAX_HOOP_STRESS, FUNCTION

    Controllable:No

    Description:choose from FUNCTION, MAX_HOOP_STRESS

  • growth_direction_xFunction defining x-component of crack growth direction

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:Function defining x-component of crack growth direction

  • growth_direction_yFunction defining y-component of crack growth direction

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:Function defining y-component of crack growth direction

  • growth_direction_zFunction defining z-component of crack growth direction

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:Function defining z-component of crack growth direction

  • growth_rateFunction defining crack growth rate

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:Function defining crack growth rate

  • growth_rate_methodFUNCTIONchoose from FUNCTION, FATIGUE

    Default:FUNCTION

    C++ Type:MooseEnum

    Options:FATIGUE, FUNCTION

    Controllable:No

    Description:choose from FUNCTION, FATIGUE

  • heal_alwaysFalseHeal previous cuts at every time step

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Heal previous cuts at every time step

  • mesh_generator_nameMesh generator for the XFEM geometric cutter.

    C++ Type:std::string

    Controllable:No

    Description:Mesh generator for the XFEM geometric cutter.

  • n_step_growth0Number of steps for crack growth

    Default:0

    C++ Type:unsigned int

    Controllable:No

    Description:Number of steps for crack growth

  • size_control0Criterion for refining elements while growing the crack

    Default:0

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Criterion for refining elements while growing the crack

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).

  • execute_onXFEM_MARKThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.

    Default:XFEM_MARK

    C++ Type:ExecFlagEnum

    Options:XFEM_MARK, FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, LINEAR_CONVERGENCE, NONLINEAR, NONLINEAR_CONVERGENCE, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, MULTIAPP_FIXED_POINT_CONVERGENCE, FINAL, CUSTOM, TRANSFER

    Controllable:No

    Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.

  • execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.

    Default:0

    C++ Type:int

    Controllable:No

    Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.

  • 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

Execution Scheduling Parameters

  • 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.

  • 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

  • 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

    Unit:(no unit assumed)

    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.

  • 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.

Material Property Retrieval Parameters

Input Files