Peridynamic Displacement L2 Difference Postprocessor

Description

Postprocessor NodalDisplacementDifferenceL2NormPD is used to calculate the L2 norm of difference between displacements and their analytic solutions using the discrete summation formulation.

$var element = document.getElementById("moose-equation-cd2fffc6-2897-429c-934e-ceafb1b5ea67");katex.render(" \\text{L2 difference} = \\int_{\\Omega} {\\left| \\mathbf{u}_{prediction} - \\mathbf{u}_{analytical} \\right|} d{\\Omega} = \\sum V_{\\mathbf{X}} \\cdot {\\left| \\mathbf{u}_{prediction}^{\\mathbf{X}} - \\mathbf{u}_{analytical}^{\\mathbf{X}} \\right|}", element, {displayMode:true,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ where $var element = document.getElementById("moose-equation-67285bd8-d98d-48a2-a16c-61f9fbfdc063");katex.render("V_{\\mathbf{X}}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ is the area or volume of material point $var element = document.getElementById("moose-equation-2f25e768-5a9b-4929-9102-405c2a2b0697");katex.render("\\mathbf{X}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ , and $var element = document.getElementById("moose-equation-e5dd27f2-0ab4-41c5-aef0-901ac3420f4c");katex.render("\\mathbf{u}_{prediction}^{\\mathbf{X}}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ and $var element = document.getElementById("moose-equation-aea0e72f-8729-4409-b081-24e91de61b33");katex.render("\\mathbf{u}_{analytical}^{\\mathbf{X}}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ are the prediction and analytical solution at material point $var element = document.getElementById("moose-equation-6f8512a9-f733-40aa-84a5-c41729a73940");katex.render("\\mathbf{X}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ .

Input Parameters

analytic_functionsThe known analytic functions for displacements
C++ Type:std::vector
Options:
Description:The known analytic functions for displacements
displacementsNonlinear variable name for the displacements
C++ Type:std::vector
Options:
Description:Nonlinear variable name for the displacements

Required Parameters

blockThe list of block ids (SubdomainID) that this object will be applied
C++ Type:std::vector
Options:
Description:The list of block ids (SubdomainID) that this object will be applied
boundaryThe list of boundary IDs from the mesh where this boundary condition applies
C++ Type:std::vector
Options:
Description:The list of boundary IDs from the mesh where this boundary condition applies
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE INITIAL LINEAR NONLINEAR TIMESTEP_END TIMESTEP_BEGIN FINAL CUSTOM TRANSFER
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
unique_node_executeFalseWhen false (default), block restricted objects will have the execute method called multiple times on a single node if the node lies on a interface between two subdomains.
Default:False
C++ Type:bool
Options:
Description:When false (default), block restricted objects will have the execute method called multiple times on a single node if the node lies on a interface between two subdomains.

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
Options:
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
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
outputsVector of output names were you would like to restrict the output of variables(s) associated with this object
C++ Type:std::vector
Options:
Description:Vector of output names were you would like to restrict the output of variables(s) associated with this object
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Options:
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
Options:
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/peridynamics/test/tests/simple_tests/2D_regularD_constH_BPD.i

modules/peridynamics/test/tests/simple_tests/2D_regularD_constH_BPD.i

# Test for bond-based peridynamic formulation
# for regular grid from generated mesh with const bond constants

# Square plate with Dirichlet boundary conditions applied
# at the left, top and bottom edges

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[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]
  [../]
[]

[BCs]
  [./left_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1003
    value = 0.0
  [../]
  [./top_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1002
    value = 0.0
  [../]
  [./bottom_y]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = 1000
    function = '-0.001*t'
  [../]
[]

[Modules/Peridynamics/Mechanics/Master]
  [./all]
    formulation = BOND
  [../]
[]

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 2e5
    poissons_ratio = 0.33
  [../]

  [./force_density]
    type = ComputeSmallStrainConstantHorizonMaterialBPD
  [../]
[]

[Functions]
  [./disp_x_anal]
    type = PiecewiseLinear
    axis = x
    x = '0 1'
    y = '0 -0.00033'
  [../]
  [./disp_y_anal]
    type = PiecewiseLinear
    axis = y
    x = '0 1'
    y = '-0.001 0'
  [../]
[]

[Postprocessors]
  [./anal_disp_L2]
    type = NodalFunctionsL2NormPD
    functions = 'disp_x_anal disp_y_anal'
  [../]
  [./disp_diff_L2]
    type = NodalDisplacementDifferenceL2NormPD
    analytic_functions = 'disp_x_anal disp_y_anal'
  [../]
[]

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

[Executioner]
  type = Transient
  solve_type = PJFNK
  line_search = none
  start_time = 0
  end_time = 1
[]

[Outputs]
  file_base = 2D_regularD_constH_BPD
  exodus = true
[]

Description
Input Parameters
Input Files

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools