Correspondence Small Strain Material

Description

Material ComputeSmallStrainNOSPD calculates the bond-associated deformation gradient used in self-stabilized non-ordinary state-based peridynamic mechanics model under small strain assumptions.

Given the bond-associated deformation gradient $var element = document.getElementById("moose-equation-52d1b2fd-1af5-4d56-9eed-98f6345ae0a7");katex.render("\\mathbf{F}_{\\mathbf{\\xi}}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ , the bond-associated total strain tensor can be obtained as

$var element = document.getElementById("moose-equation-fad6a00f-2e9f-4d19-be2c-510682dcda40");katex.render(" \\boldsymbol{\\epsilon}_{\\mathbf{\\xi}, \\text{total}} = \\frac{1}{2} \\left( \\mathbf{F}_{\\mathbf{\\xi}}^{T} \\mathbf{F}_{\\mathbf{\\xi}} - \\mathbf{I} \\right)", element, {displayMode:true,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$

and the elastic strain tensor can be obtained by subtracting the eigen strain tensor from the total strain tensor as

$var element = document.getElementById("moose-equation-806f6a73-3db2-4d8c-94b4-562b5bc3902d");katex.render(" \\boldsymbol{\\epsilon}_{\\mathbf{\\xi}, \\text{elastic}} = \\boldsymbol{\\epsilon}_{\\mathbf{\\xi}, \\text{total}} - \\boldsymbol{\\epsilon}_{\\mathbf{\\xi}, \\text{eigen}}", element, {displayMode:true,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$

The computed bond-associated elastic strain tensor is then used in stress calculator to calculate bond-associated stress tensor.

Note that Material ComputeSmallStrainNOSPD can be used in general 3D and plane strain modeling and simulation.

Input Parameters

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
computeTrueWhen false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
Default:True
C++ Type:bool
Options:
Description:When false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeSubdomainProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
Default:NONE
C++ Type:MooseEnum
Options:NONE ELEMENT SUBDOMAIN
Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeSubdomainProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
eigenstrain_namesList of eigenstrains to be applied in this strain calculation
C++ Type:std::vector
Options:
Description:List of eigenstrains to be applied in this strain calculation
plane_strainFalsePlane strain problem or not, this will affect the mechanical stretch calculation for problem with eigenstrains
Default:False
C++ Type:bool
Options:
Description:Plane strain problem or not, this will affect the mechanical stretch calculation for problem with eigenstrains
temperatureNonlinear variable name for the temperature
C++ Type:std::vector
Options:
Description:Nonlinear variable name for the temperature

Optional Parameters

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.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Options:
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
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

output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)
C++ Type:std::vector
Options:
Description:List of material properties, from this material, to output (outputs must also be defined to an output type)
outputsnone Vector of output names were you would like to restrict the output of variables(s) associated with this object
Default:none
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

Outputs Parameters

Input Files

modules/peridynamics/test/tests/failure_tests/2D_stress_failure_NOSPD.i
modules/peridynamics/test/tests/plane_stress/conventional_planestress_NOSPD.i
modules/peridynamics/test/tests/jacobian_check/3D_mechanics_smallstrain_NOSPD.i
modules/peridynamics/test/tests/simple_tests/2D_small_strain_NOSPD.i

modules/peridynamics/test/tests/failure_tests/2D_stress_failure_NOSPD.i


[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  type = PeridynamicsMesh
  horizon_number = 3
  cracks_start = '0.25 0.5 0'
  cracks_end = '0.75 0.5 0'

  [./gmg]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 8
    ny = 8
  [../]
  [./gpd]
    type = MeshGeneratorPD
    input = gmg
    retain_fe_mesh = false
  [../]
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
[]

[AuxVariables]
  [./damage]
  [../]
  [./intact_bonds_num]
  [../]
  [./critical_stress]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./bond_status]
    type = RankTwoBasedFailureCriteriaNOSPD
    variable = bond_status
    rank_two_tensor = stress
    critical_variable = critical_stress
    failure_criterion = VonMisesStress
  [../]
[]

[UserObjects]
  [./damage]
    type = NodalDamageIndexPD
    variable = damage
  [../]
  [./intact_bonds]
    type = NodalNumIntactBondsPD
    variable = intact_bonds_num
  [../]
[]

[ICs]
  [./critical_stretch]
    type = ConstantIC
    variable = critical_stress
    value = 150
  [../]
[]

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

  [./rbm_x]
    type = RBMPresetOldValuePD
    variable = disp_x
    boundary = 999
  [../]
  [./rbm_y]
    type = RBMPresetOldValuePD
    variable = disp_y
    boundary = 999
  [../]
[]

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

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 2e5
    poissons_ratio = 0.33
  [../]
  [./strain]
    type = ComputeSmallStrainNOSPD
  [../]
  [./stress]
    type = ComputeLinearElasticStress
  [../]
[]

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

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

  [./Quadrature]
    type = GAUSS_LOBATTO
    order = FIRST
  [../]
[]

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

modules/peridynamics/test/tests/plane_stress/conventional_planestress_NOSPD.i

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  type = PeridynamicsMesh
  horizon_number = 3

  [./gmg]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 10
    ny = 10
  [../]
  [./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
  [../]
  [./left_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1003
    value = 0.0
  [../]
  [./right_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1001
    value = 0.001
  [../]
[]

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

[Materials]
  [./elasticity_tensor]
    type = ComputePlaneStressIsotropicElasticityTensor
    youngs_modulus = 2.1e8
    poissons_ratio = 0.3
  [../]
  [./strain]
    type = ComputeSmallStrainNOSPD
  [../]
  [./stress]
    type = ComputeLinearElasticStress
  [../]
[]

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

[Executioner]
  type = Transient
  solve_type = PJFNK
  line_search = none
  start_time = 0
  end_time = 1
  nl_rel_tol = 1e-10

  [./Quadrature]
    type = GAUSS_LOBATTO
    order = FIRST
  [../]
[]

[Outputs]
  file_base = conventional_planestress_NOSPD
  exodus = true
[]

modules/peridynamics/test/tests/jacobian_check/3D_mechanics_smallstrain_NOSPD.i

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
  full_jacobian = true
[]

[Mesh]
  type = PeridynamicsMesh
  horizon_number = 2

  [./gmg]
    type = GeneratedMeshGenerator
    dim = 3
    nx = 3
    ny = 3
    nz = 3
  [../]
  [./gpd]
    type = MeshGeneratorPD
    input = gmg
    retain_fe_mesh = false
  [../]
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]

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

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 2e5
    poissons_ratio = 0.0
  [../]
  [./strain]
    type = ComputeSmallStrainNOSPD
  [../]
  [./stress]
    type = ComputeLinearElasticStress
  [../]
[]

[Preconditioning]
  [./SMP]
    type = SMP
    full = true
    petsc_options_iname = '-ksp_type -pc_type -snes_type'
    petsc_options_value = 'bcgs bjacobi test'
  [../]
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  end_time = 1
  dt = 1
  num_steps = 1

  [./Quadrature]
    type = GAUSS_LOBATTO
    order = FIRST
  [../]
[]

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


[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  type = PeridynamicsMesh
  horizon_number = 3

  [./gmg]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 8
    ny = 8
  [../]
  [./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
  [../]
  [./left_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1003
    value = 0.0
  [../]
  [./right_x]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = 1001
    function = '0.001*t'
  [../]
[]

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

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 2.1e8
    poissons_ratio = 0.3
  [../]
  [./strain]
    type = ComputeSmallStrainNOSPD
  [../]
  [./stress]
    type = ComputeLinearElasticStress
  [../]
[]

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

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

  [./Quadrature]
    type = GAUSS_LOBATTO
    order = FIRST
  [../]
[]

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

Child Objects

modules/peridynamics/include/materials/ComputeForceStabilizedSmallStrainNOSPD.h
modules/peridynamics/include/materials/ComputePlaneSmallStrainNOSPD.h

modules/peridynamics/include/materials/ComputeForceStabilizedSmallStrainNOSPD.h

// This file is part of the MOOSE framework
// https://www.mooseframework.org
//
// All rights reserved, see COPYRIGHT for full restrictions
// https://github.com/idaholab/moose/blob/master/COPYRIGHT
//
// Licensed under LGPL 2.1, please see LICENSE for details
// https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "ComputeSmallStrainNOSPD.h"

class ComputeForceStabilizedSmallStrainNOSPD;

template <>
InputParameters validParams<ComputeForceStabilizedSmallStrainNOSPD>();

/**
 * Material class for conventional correspondence material model with fictitious force stabilization
 * for small strain
 */
class ComputeForceStabilizedSmallStrainNOSPD : public ComputeSmallStrainNOSPD
{
public:
  ComputeForceStabilizedSmallStrainNOSPD(const InputParameters & parameters);

protected:
  virtual void computeQpDeformationGradient() override;

  /// Elasticity tensor
  const MaterialProperty<RankFourTensor> & _Cijkl;

  /// Material property for fictitious force
  MaterialProperty<Real> & _sf_coeff;
};

modules/peridynamics/include/materials/ComputePlaneSmallStrainNOSPD.h

// This file is part of the MOOSE framework
// https://www.mooseframework.org
//
// All rights reserved, see COPYRIGHT for full restrictions
// https://github.com/idaholab/moose/blob/master/COPYRIGHT
//
// Licensed under LGPL 2.1, please see LICENSE for details
// https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "ComputeSmallStrainNOSPD.h"
#include "RankTwoTensor.h"

class ComputePlaneSmallStrainNOSPD;

template <>
InputParameters validParams<ComputePlaneSmallStrainNOSPD>();

/**
 * Material class for 2D correspondence material model for small strain: plane strain, generalized
 * plane strain, weak plane stress
 */
class ComputePlaneSmallStrainNOSPD : public ComputeSmallStrainNOSPD
{
public:
  ComputePlaneSmallStrainNOSPD(const InputParameters & parameters);

protected:
  virtual void computeQpTotalStrain() override;

  /**
   * Function to compute out-of-plane component of strain tensor for generalized plane strain and
   * weak plane stress
   * @return The value of out-of-plane strain
   */
  Real computeOutOfPlaneStrain();

private:
  ///@{ Scalar out-of-plane strain for generalized plane strain
  const bool _scalar_out_of_plane_strain_coupled;
  const VariableValue & _scalar_out_of_plane_strain;
  ///@}

  ///@{ Out-of-plane strain for weak plane stress
  const bool _out_of_plane_strain_coupled;
  const VariableValue & _out_of_plane_strain;
  ///@}
};

Description
Input Parameters
Input Files
Child Objects

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools