Peridynamic Generalized Plane Strain Action System

Description

This action sets up a generalized plane strain model, including kernel to provide coupled off-diagonal Jacobian entries, scalar kernel to provide residual and diagonal Jacobian, and user object to compute residual and diagonal Jacobian for scalar variable.

Constructed MooseObjects

Table 1: Correspondence Among Action Functionality and MooseObjects for the GeneralizedPlaneStrainPD Action

Functionality	Replaced Classes	Associated Parameters
Out-of-plane scalar variable equilibrium condition	Generalized Plane Strain ScalarKernel	`generalized_plane_strain_uo`: UserObject name of the GeneralizedPlaneStrainUserObjectBasePD
Residual and diagonal Jacobian calculation for scalar out-of-plane strain variables	Ordinary State-based Generalized Plane Strain UserObject or Non-Ordinary State-based Generalized Plane Strain UserObject	`out_of_plane_stress_variable`: Auxiliary variable name for out-of-plane stress for ordinary state-based generalized plane strain model; or `none` is required for non-ordinary state-based generalized plane strain model\|
Scalar out-of-plan strain coupling with in-plane field variables	Ordinary State-based Generalized Plane Strain Off-diagonal Kernel or Non-Ordinary State-based Generalized Plane Strain Off-diagonal Kernel	`scalar_out_of_plane_strain`: a string of scalar variable for the out-of-plane strain direction; `displacements` : a string of the displacement field (in-plane) variables; `temperature`: a string of the temperature field variable

Example Input Syntax

Subblocks

The subblocks of the GeneralizedPlaneStrain action are what triggers MOOSE objects to be built. If none of the mechanics is subdomain restricted a single subblock will be used

[Modules]
  [Peridynamics]
    [Mechanics]
      [GeneralizedPlaneStrain]
        [all]
          formulation = ORDINARY_STATE
          out_of_plane_stress_variable = stress_zz
        []
      []
    []
  []
[]

If different generalized plane strain models are needed, multiple subblocks with subdomain restrictions can be used.

[Modules]
  [Peridynamics]
    [Mechanics]
      [GeneralizedPlaneStrain]
        [block1]
          formulation = ORDINARY_STATE
          scalar_out_of_plane_strain = scalar_strain_zz1
          out_of_plane_stress_variable = stress_zz1
          block = 1001
        []
        [block2]
          formulation = ORDINARY_STATE
          scalar_out_of_plane_strain = scalar_strain_zz2
          out_of_plane_stress_variable = stress_zz2
          block = 1002
        []
      []
    []
  []
[]

Parameters supplied at the [Modules/Peridynamics/GeneralizedPlaneStrain] level act as defaults for the GeneralizedPlaneStrain action subblocks.

Input Parameters

displacementsNonlinear variable name for the displacements
C++ Type:std::vector<VariableName>
Controllable:No
Description:Nonlinear variable name for the displacements
scalar_out_of_plane_strainScalar variable for strain in the out-of-plane direction
C++ Type:VariableName
Controllable:No
Description:Scalar variable for strain in the out-of-plane direction

Required Parameters

active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified only the blocks named will be visited and made active
blockList of ids of the blocks (subdomains) that the GeneralizedPlaneStrainActionPD will be applied to
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:List of ids of the blocks (subdomains) that the GeneralizedPlaneStrainActionPD will be applied to
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.
eigenstrain_namesList of eigenstrains to be applied in this strain calculation
C++ Type:std::vector<MaterialPropertyName>
Controllable:No
Description:List of eigenstrains to be applied in this strain calculation
factor1Scale factor applied to prescribed out-of-plane pressure
Default:1
C++ Type:double
Controllable:No
Description:Scale factor applied to prescribed out-of-plane pressure
formulationNONORDINARY_STATEPeridynamic formulation options
Default:NONORDINARY_STATE
C++ Type:MooseEnum
Options:ORDINARY_STATE, NONORDINARY_STATE
Controllable:No
Description:Peridynamic formulation options
full_jacobianFalseParameter to set whether to use the nonlocal full Jacobian formulation for the scalar components
Default:False
C++ Type:bool
Controllable:No
Description:Parameter to set whether to use the nonlocal full Jacobian formulation for the scalar components
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified blocks matching these identifiers will be skipped.
out_of_plane_pressure0Function used to prescribe pressure in the out-of-plane direction
Default:0
C++ Type:FunctionName
Controllable:No
Description:Function used to prescribe pressure in the out-of-plane direction
out_of_plane_stress_variableName of out-of-plane stress auxiliary variable
C++ Type:VariableName
Controllable:No
Description:Name of out-of-plane stress auxiliary variable
strainSMALLStrain formulation
Default:SMALL
C++ Type:MooseEnum
Options:SMALL, FINITE
Controllable:No
Description:Strain formulation
temperatureNonlinear variable for the temperature
C++ Type:VariableName
Controllable:No
Description:Nonlinear variable for the temperature

Optional Parameters

Associated Actions

Available Actions

Peridynamics App
GeneralizedPlaneStrainActionPDClass for setting up the Kernel, ScalarKernel, and UserObject for peridynamic generalized plane strain model

(modules/peridynamics/test/tests/generalized_plane_strain/generalized_plane_strain_OSPD.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
  scalar_out_of_plane_strain = scalar_strain_zz
[]

[Mesh]
  type = PeridynamicsMesh
  horizon_number = 3

  [gmg]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 4
    ny = 4
  []
  [gpd]
    type = MeshGeneratorPD
    input = gmg
    retain_fe_mesh = false
  []
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []
  [scalar_strain_zz]
    order = FIRST
    family = SCALAR
  []
[]

[AuxVariables]
  [temp]
    order = FIRST
    family = LAGRANGE
  []

  [stress_zz]
    order = FIRST
    family = LAGRANGE
  []
[]

[Modules/Peridynamics/Mechanics]
  [Master]
    [all]
      formulation = ORDINARY_STATE
    []
  []
  [GeneralizedPlaneStrain]
    [all]
      formulation = ORDINARY_STATE
      out_of_plane_stress_variable = stress_zz
    []
  []
[]

[AuxKernels]
  [tempfuncaux]
    type = FunctionAux
    variable = temp
    function = tempfunc
    use_displaced_mesh = false
  []

  [stress_zz]
    type = NodalRankTwoPD
    variable = stress_zz

    poissons_ratio = 0.3
    youngs_modulus = 1e6
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5

    rank_two_tensor = stress
    output_type = component
    index_i = 2
    index_j = 2
  []
[]

[Postprocessors]
  [react_z]
    type = NodalVariableIntegralPD
    variable = stress_zz
  []
[]

[Functions]
  [tempfunc]
    type = ParsedFunction
    expression = '(1-x)*t'
  []
[]

[BCs]
  [bottom_x]
    type = DirichletBC
    boundary = 1000
    variable = disp_x
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    boundary = 1000
    variable = disp_y
    value = 0.0
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = 0.3
    youngs_modulus = 1e6
  []

  [force_density]
    type = ComputeSmallStrainVariableHorizonMaterialOSPD
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5
  []
[]

[Preconditioning]
  [SMP]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient

  solve_type = PJFNK
  line_search = none

  nl_rel_tol = 1e-15
  nl_abs_tol = 1e-09

  start_time = 0.0
  end_time = 1.0

  use_pre_SMO_residual = true
[]

[Outputs]
  exodus = true
  file_base = generalized_plane_strain_OSPD
[]

(modules/peridynamics/test/tests/generalized_plane_strain/generalized_plane_strain_squares_OSPD.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  type = PeridynamicsMesh
  horizon_number = 3

  [fmg]
    type = FileMeshGenerator
    file = squares.e
  []
  [gpd]
    type = MeshGeneratorPD
    input = fmg
    retain_fe_mesh = false
  []
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []

  [scalar_strain_zz1]
    order = FIRST
    family = SCALAR
  []
  [scalar_strain_zz2]
    order = FIRST
    family = SCALAR
  []
[]

[AuxVariables]
  [temp]
    order = FIRST
    family = LAGRANGE
  []

  [stress_zz1]
    order = FIRST
    family = LAGRANGE
  []
  [stress_zz2]
    order = FIRST
    family = LAGRANGE
  []
[]

[Modules/Peridynamics/Mechanics]
  [Master]
    [block1]
      formulation = ORDINARY_STATE
      block = 1001
    []
    [block2]
      formulation = ORDINARY_STATE
      block = 1002
    []
  []
  [GeneralizedPlaneStrain]
    [block1]
      formulation = ORDINARY_STATE
      scalar_out_of_plane_strain = scalar_strain_zz1
      out_of_plane_stress_variable = stress_zz1
      block = 1001
    []
    [block2]
      formulation = ORDINARY_STATE
      scalar_out_of_plane_strain = scalar_strain_zz2
      out_of_plane_stress_variable = stress_zz2
      block = 1002
    []
  []
[]

[AuxKernels]
  [tempfuncaux]
    type = FunctionAux
    variable = temp
    function = tempfunc
    use_displaced_mesh = false
  []

  [stress_zz1]
    type = NodalRankTwoPD
    variable = stress_zz1
    rank_two_tensor = stress
    scalar_out_of_plane_strain = scalar_strain_zz1

    poissons_ratio = 0.3
    youngs_modulus = 1e6
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5

    output_type = component
    index_i = 2
    index_j = 2
    block = 1001
  []

  [stress_zz2]
    type = NodalRankTwoPD
    variable = stress_zz2
    scalar_out_of_plane_strain = scalar_strain_zz2

    poissons_ratio = 0.3
    youngs_modulus = 1e6
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5

    rank_two_tensor = stress
    output_type = component
    index_i = 2
    index_j = 2
    block = 1002
  []
[]

[Postprocessors]
  [react_z1]
    type = NodalVariableIntegralPD
    variable = stress_zz1
    block = 1001
  []
  [react_z2]
    type = NodalVariableIntegralPD
    variable = stress_zz2
    block = 1002
  []
[]

[Functions]
  [tempfunc]
    type = ParsedFunction
    expression = '(1-x)*t'
  []
[]

[BCs]
  [bottom1_x]
    type = DirichletBC
    boundary = 1001
    variable = disp_x
    value = 0.0
  []
  [bottom1_y]
    type = DirichletBC
    boundary = 1001
    variable = disp_y
    value = 0.0
  []

  [bottom2_x]
    type = DirichletBC
    boundary = 1002
    variable = disp_x
    value = 0.0
  []
  [bottom2_y]
    type = DirichletBC
    boundary = 1002
    variable = disp_y
    value = 0.0
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = 0.3
    youngs_modulus = 1e6
    block = '1001 1002'
  []

  [force_density1]
    type = ComputeSmallStrainVariableHorizonMaterialOSPD
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5
    scalar_out_of_plane_strain = scalar_strain_zz1
    block = 1001
  []
  [force_density2]
    type = ComputeSmallStrainVariableHorizonMaterialOSPD
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5
    scalar_out_of_plane_strain = scalar_strain_zz2
    block = 1002
  []
[]

[Preconditioning]
  [SMP]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient

  solve_type = PJFNK
  line_search = none

  l_tol = 1e-8
  nl_rel_tol = 1e-15
  nl_abs_tol = 1e-09

  start_time = 0.0
  end_time = 1.0

  use_pre_SMO_residual = true
[]

[Outputs]
  exodus = true
  file_base = generalized_plane_strain_squares_OSPD
[]

Description
Constructed MooseObjects
Example Input Syntax
Input Parameters
Associated Actions
Available Actions