https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | List of all members
GeneralizedPlaneStrainOffDiagOSPD Class Reference

Kernel class for coupled off diagonal Jacobian entries of ordinary state-based peridynamic generalized plane strain model. More...

#include <GeneralizedPlaneStrainOffDiagOSPD.h>

Inheritance diagram for GeneralizedPlaneStrainOffDiagOSPD:
[legend]

Public Types

typedef DerivativeMaterialPropertyNameInterface::SymbolName SymbolName
 

Public Member Functions

 GeneralizedPlaneStrainOffDiagOSPD (const InputParameters &parameters)
 
virtual void computeOffDiagJacobian (unsigned int jvar) override
 
virtual void computeLocalOffDiagJacobian (unsigned int, unsigned int)
 Function to compute local contribution to the off-diagonal Jacobian at the current nodes.
 
virtual void computePDNonlocalOffDiagJacobian (unsigned int, unsigned int)
 Function to compute nonlocal contribution to the off-diagonal Jacobian at the current nodes.
 
virtual void initialSetup () override
 
virtual void prepare () override
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialPropertyByName (const std::string &name)
 
void validateDerivativeMaterialPropertyBase (const std::string &base)
 
const MaterialPropertyName derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual void computeLocalResidual () override
 
virtual void computeOffDiagJacobianScalar (unsigned int jvar_num) override
 
void computeDispFullOffDiagJacobianScalar (unsigned int component, unsigned int jvar_num)
 Function to compute the full off diagonal Jacobian for coupling between displacements and scalar variable.
 
void computeDispPartialOffDiagJacobianScalar (unsigned int component, unsigned int jvar_num)
 Function to compute partial off diagonal Jacobian for coupling between displacements and scalar variable.
 
void computeTempOffDiagJacobianScalar (unsigned int jvar_num)
 Function to compute off disgonal Jacobian for coupling between temperature and scalar variable.
 

Protected Attributes

const MaterialProperty< RankFourTensor > & _Cijkl
 Material point based material property.
 
const unsigned int _scalar_out_of_plane_strain_var_num
 The variable number of the scalar out-of-plane strain variable.
 
DenseMatrix< Number > _ken
 
DenseMatrix< Number > _kne
 
std::vector< MooseVariable * > _disp_var
 displacement variables
 
unsigned int _ndisp
 number of displacement components
 
const std::vector< RealGradient > * _orientation
 Vector of bond in current configuration.
 
std::vector< dof_id_type > _ivardofs
 Current variable dof numbers for nodes i and j.
 
std::vector< Real > _weights
 weights used for the current element to obtain the nodal stress
 
RealGradient _current_vec
 Vector of bond in current configuration.
 
RealGradient _current_unit_vec
 Unit vector of bond in current configuration.
 
const MaterialProperty< Real > & _bond_local_dfdE
 Bond based material properties.
 
const MaterialProperty< Real > & _bond_nonlocal_dfdE
 
const MaterialProperty< Real > & _alpha
 
const bool _temp_coupled
 Temperature variable.
 
MooseVariable * _temp_var
 
const bool _out_of_plane_strain_coupled
 Parameters for out-of-plane strain in weak plane stress formulation.
 
MooseVariable * _out_of_plane_strain_var
 

Private Member Functions

bool haveMaterialProperty (const std::string &prop_name)
 
std::vector< VariableName > buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)
 
void validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)
 
bool isNotObjectVariable (const VariableName &name)
 

Private Attributes

FEProblemBase & _dmi_fe_problem
 

Detailed Description

Kernel class for coupled off diagonal Jacobian entries of ordinary state-based peridynamic generalized plane strain model.

Definition at line 18 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Constructor & Destructor Documentation

◆ GeneralizedPlaneStrainOffDiagOSPD()

GeneralizedPlaneStrainOffDiagOSPD::GeneralizedPlaneStrainOffDiagOSPD ( const InputParameters &  parameters)

Definition at line 31 of file GeneralizedPlaneStrainOffDiagOSPD.C.

33 : MechanicsBasePD(parameters),
34 _bond_local_dfdE(getMaterialProperty<Real>("bond_local_dfdE")),
35 _bond_nonlocal_dfdE(getMaterialProperty<Real>("bond_nonlocal_dfdE")),
36 _alpha(getMaterialProperty<Real>("thermal_expansion_coeff")),
37 _Cijkl(getMaterialProperty<RankFourTensor>("elasticity_tensor")),
38 _scalar_out_of_plane_strain_var_num(coupledScalar("scalar_out_of_plane_strain"))
39{
40 // Consistency check
41 if (_disp_var.size() != 2)
42 mooseError("GeneralizedPlaneStrain only works for two dimensional case!");
43}
mooseError
void mooseError(Args &&... args)
GeneralizedPlaneStrainOffDiagOSPD::_bond_nonlocal_dfdE
const MaterialProperty< Real > & _bond_nonlocal_dfdE
Definition GeneralizedPlaneStrainOffDiagOSPD.h:53
GeneralizedPlaneStrainOffDiagOSPD::_bond_local_dfdE
const MaterialProperty< Real > & _bond_local_dfdE
Bond based material properties.
Definition GeneralizedPlaneStrainOffDiagOSPD.h:52
GeneralizedPlaneStrainOffDiagOSPD::_Cijkl
const MaterialProperty< RankFourTensor > & _Cijkl
Material point based material property.
Definition GeneralizedPlaneStrainOffDiagOSPD.h:58
GeneralizedPlaneStrainOffDiagOSPD::_scalar_out_of_plane_strain_var_num
const unsigned int _scalar_out_of_plane_strain_var_num
The variable number of the scalar out-of-plane strain variable.
Definition GeneralizedPlaneStrainOffDiagOSPD.h:61
MechanicsBasePD
Base kernel class for peridynamic solid mechanics models.
Definition MechanicsBasePD.h:19
MechanicsBasePD::_disp_var
std::vector< MooseVariable * > _disp_var
displacement variables
Definition MechanicsBasePD.h:51

Member Function Documentation

◆ computeDispFullOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar ( unsigned int  component,
unsigned int  jvar_num 
)
protected

Function to compute the full off diagonal Jacobian for coupling between displacements and scalar variable.

Parameters
componentThe index of displacement component
jvar_numThe coupled scalar variable number

Definition at line 68 of file GeneralizedPlaneStrainOffDiagOSPD.C.

70{
71 prepareMatrixTag(_assembly, _var.number(), jvar_num, _ken);
72 prepareMatrixTag(_assembly, jvar_num, _var.number(), _kne);
73 MooseVariableScalar & jvar = _sys.getScalarVariable(_tid, jvar_num);
74
75 // LOCAL jacobian contribution
76 // fill in the column corresponding to the scalar variable from bond ij
77 for (unsigned int i = 0; i < _nnodes; ++i)
78 for (unsigned int j = 0; j < jvar.order(); ++j)
79 _ken(i, j) +=
80 (i == j ? -1 : 1) * _current_unit_vec(component) * _bond_local_dfdE[0] * _bond_status;
81
82 // NONLOCAL jacobian contribution
83 std::vector<RankTwoTensor> shape(_nnodes), dgrad(_nnodes);
84
85 // for off-diagonal low components
86 if (_bond_status > 0.5)
87 {
88 for (unsigned int nd = 0; nd < _nnodes; nd++)
89 {
90 if (_dim == 2)
91 shape[nd](2, 2) = dgrad[nd](2, 2) = 1.0;
92 // calculation of jacobian contribution to current node's neighbors
93 std::vector<dof_id_type> ivardofs(_nnodes);
94 ivardofs[nd] = _ivardofs[nd];
95 std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
96 std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
97
98 Real vol_nb, weight_nb;
99 RealGradient origin_vec_nb, current_vec_nb;
100
101 for (unsigned int nb = 0; nb < neighbors.size(); nb++)
102 if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[nb])) > 0.5)
103 {
104 ivardofs[1 - nd] =
105 _pdmesh.nodePtr(neighbors[nb])->dof_number(_sys.number(), _var.number(), 0);
106 vol_nb = _pdmesh.getNodeVolume(neighbors[nb]);
107
108 // obtain bond nb's origin length and current orientation
109 origin_vec_nb = _pdmesh.getNodeCoord(neighbors[nb]) -
110 _pdmesh.getNodeCoord(_current_elem->node_id(nd));
111
112 for (unsigned int k = 0; k < _dim; k++)
113 current_vec_nb(k) = origin_vec_nb(k) +
114 _disp_var[k]->getNodalValue(*_pdmesh.nodePtr(neighbors[nb])) -
115 _disp_var[k]->getNodalValue(*_current_elem->node_ptr(nd));
116
117 weight_nb = _horizon_radius[nd] / origin_vec_nb.norm();
118 // prepare shape tensor and deformation gradient tensor for current node
119 for (unsigned int k = 0; k < _dim; k++)
120 for (unsigned int l = 0; l < _dim; l++)
121 {
122 shape[nd](k, l) += weight_nb * origin_vec_nb(k) * origin_vec_nb(l) * vol_nb;
123 dgrad[nd](k, l) += weight_nb * current_vec_nb(k) * origin_vec_nb(l) * vol_nb;
124 }
125
126 // cache the nonlocal jacobian contribution
127 _local_ke.resize(_ken.m(), _ken.n());
128 _local_ke.zero();
129 _local_ke(0, 0) = (nd == 0 ? -1 : 1) * current_vec_nb(component) / current_vec_nb.norm() *
130 _bond_nonlocal_dfdE[nd] / origin_vec_nb.norm() * vol_nb * _bond_status;
131 _local_ke(1, 0) = (nd == 0 ? 1 : -1) * current_vec_nb(component) / current_vec_nb.norm() *
132 _bond_nonlocal_dfdE[nd] / origin_vec_nb.norm() * vol_nb * _bond_status;
133
134 addJacobian(_assembly, _local_ke, ivardofs, jvar.dofIndices(), _var.scalingFactor());
135 }
136
137 // finalize the shape tensor and deformation gradient tensor for node_i
138 if (shape[nd].det() == 0.)
139 mooseError("Singular shape tensor is detected in GeneralizedPlaneStrainOffDiagOSPD! Use "
140 "SingularShapeTensorEliminatorUserObjectPD to avoid singular shape tensor");
141
142 shape[nd] = shape[nd].inverse();
143 dgrad[nd] = dgrad[nd] * shape[nd];
144 }
145 }
146
147 // off-diagonal jacobian entries on the row
148 Real dEidUi =
149 -_node_vol[1] * _horizon_radius[0] / _origin_vec.norm() *
150 (_Cijkl[0](2, 2, 0, 0) * (_origin_vec(0) * shape[0](0, 0) + _origin_vec(1) * shape[0](1, 0)) *
151 dgrad[0](component, 0) +
152 _Cijkl[0](2, 2, 1, 1) * (_origin_vec(0) * shape[0](0, 1) + _origin_vec(1) * shape[0](1, 1)) *
153 dgrad[0](component, 1));
154 Real dEjdUj =
155 _node_vol[0] * _horizon_radius[1] / _origin_vec.norm() *
156 (_Cijkl[0](2, 2, 0, 0) * (_origin_vec(0) * shape[1](0, 0) + _origin_vec(1) * shape[1](1, 0)) *
157 dgrad[1](component, 0) +
158 _Cijkl[0](2, 2, 1, 1) * (_origin_vec(0) * shape[1](0, 1) + _origin_vec(1) * shape[1](1, 1)) *
159 dgrad[1](component, 1));
160
161 // fill in the row corresponding to the scalar variable
162 _kne(0, 0) += (dEidUi * _node_vol[0] - dEjdUj * _node_vol[1]) * _bond_status; // node i
163 _kne(0, 1) += (dEjdUj * _node_vol[1] - dEidUi * _node_vol[0]) * _bond_status; // node j
164
165 accumulateTaggedLocalMatrix(_assembly, _var.number(), jvar_num, _ken);
166 accumulateTaggedLocalMatrix(_assembly, jvar_num, _var.number(), _kne);
167}
MechanicsBasePD::_current_unit_vec
RealGradient _current_unit_vec
Unit vector of bond in current configuration.
Definition MechanicsBasePD.h:79
MechanicsBasePD::_ivardofs
std::vector< dof_id_type > _ivardofs
Current variable dof numbers for nodes i and j.
Definition MechanicsBasePD.h:70
MooseVariableScalar::dofIndices
virtual const std::vector< dof_id_type > & dofIndices() const
MooseVariableScalar::order
libMesh::Order order() const
libMesh::DenseMatrix::n
unsigned int n() const
libMesh::DenseMatrix::m
unsigned int m() const
RankTwoScalarTools::component
T component(const RankTwoTensorTempl< T > &r2tensor, unsigned int i, unsigned int j)
Definition RankTwoScalarTools.h:67
libMesh::RealGradient
RealVectorValue RealGradient
Definition ElasticityTensorTools.h:18
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

Referenced by computeOffDiagJacobianScalar().

◆ computeDispPartialOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeDispPartialOffDiagJacobianScalar ( unsigned int  component,
unsigned int  jvar_num 
)
protected

Function to compute partial off diagonal Jacobian for coupling between displacements and scalar variable.

Parameters
componentThe index of displacement component
jvar_numThe coupled scalar variable number

Definition at line 170 of file GeneralizedPlaneStrainOffDiagOSPD.C.

172{
173 prepareMatrixTag(_assembly, _var.number(), jvar_num, _ken);
174 prepareMatrixTag(_assembly, jvar_num, _var.number(), _kne);
175 MooseVariableScalar & jvar = _sys.getScalarVariable(_tid, jvar_num);
176
177 // fill in the column corresponding to the scalar variable from bond ij
178 for (unsigned int i = 0; i < _nnodes; ++i)
179 for (unsigned int j = 0; j < jvar.order(); ++j)
180 {
181 _ken(i, j) +=
182 (i == j ? -1 : 1) * _current_unit_vec(component) * _bond_local_dfdE[0] * _bond_status;
183 _kne(j, i) +=
184 (i == j ? -1 : 1) * _current_unit_vec(component) * _bond_local_dfdE[0] * _bond_status;
185 }
186
187 accumulateTaggedLocalMatrix(_assembly, _var.number(), jvar_num, _ken);
188 accumulateTaggedLocalMatrix(_assembly, jvar_num, _var.number(), _kne);
189}

Referenced by computeOffDiagJacobianScalar().

◆ computeLocalOffDiagJacobian()

virtual void MechanicsBasePD::computeLocalOffDiagJacobian ( unsigned int  ,
unsigned int   
)
inlinevirtualinherited

Function to compute local contribution to the off-diagonal Jacobian at the current nodes.

Parameters
coupled_componentThe coupled variable number

Reimplemented in ForceStabilizedSmallStrainMechanicsNOSPD, HorizonStabilizedFormIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormISmallStrainMechanicsNOSPD, WeakPlaneStressNOSPD, MechanicsBPD, HorizonStabilizedFormIISmallStrainMechanicsNOSPD, and MechanicsOSPD.

Definition at line 31 of file MechanicsBasePD.h.

33 {
34 }

Referenced by MechanicsBasePD::computeOffDiagJacobian().

◆ computeLocalResidual()

virtual void GeneralizedPlaneStrainOffDiagOSPD::computeLocalResidual ( )
inlineoverrideprotectedvirtual

Definition at line 26 of file GeneralizedPlaneStrainOffDiagOSPD.h.

26{};

◆ computeOffDiagJacobian()

void MechanicsBasePD::computeOffDiagJacobian ( unsigned int  jvar)
overridevirtualinherited

Definition at line 73 of file MechanicsBasePD.C.

74{
75 prepare();
76
77 if (jvar_num == _var.number())
78 computeJacobian();
79 else
80 {
81 unsigned int coupled_component = 0;
82 bool active = false;
83
84 for (unsigned int i = 0; i < _dim; ++i)
85 if (jvar_num == _disp_var[i]->number())
86 {
87 coupled_component = i;
88 active = true;
89 }
90
91 if (_temp_coupled && jvar_num == _temp_var->number())
92 active = true;
93
94 if (_out_of_plane_strain_coupled && jvar_num == _out_of_plane_strain_var->number())
95 active = true;
96
97 if (active)
98 {
99 prepareMatrixTag(_assembly, _var.number(), jvar_num);
100 computeLocalOffDiagJacobian(jvar_num, coupled_component);
101 accumulateTaggedLocalMatrix();
102
103 if (_use_full_jacobian)
104 computePDNonlocalOffDiagJacobian(jvar_num, coupled_component);
105 }
106 }
107}
MechanicsBasePD::_out_of_plane_strain_coupled
const bool _out_of_plane_strain_coupled
Parameters for out-of-plane strain in weak plane stress formulation.
Definition MechanicsBasePD.h:62
MechanicsBasePD::prepare
virtual void prepare() override
Definition MechanicsBasePD.C:51
MechanicsBasePD::_temp_var
MooseVariable * _temp_var
Definition MechanicsBasePD.h:55
MechanicsBasePD::computePDNonlocalOffDiagJacobian
virtual void computePDNonlocalOffDiagJacobian(unsigned int, unsigned int)
Function to compute nonlocal contribution to the off-diagonal Jacobian at the current nodes.
Definition MechanicsBasePD.h:41
MechanicsBasePD::_out_of_plane_strain_var
MooseVariable * _out_of_plane_strain_var
Definition MechanicsBasePD.h:63
MechanicsBasePD::computeLocalOffDiagJacobian
virtual void computeLocalOffDiagJacobian(unsigned int, unsigned int)
Function to compute local contribution to the off-diagonal Jacobian at the current nodes.
Definition MechanicsBasePD.h:31
MechanicsBasePD::_temp_coupled
const bool _temp_coupled
Temperature variable.
Definition MechanicsBasePD.h:54
MooseVariableFE::number
unsigned int number() const

◆ computeOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeOffDiagJacobianScalar ( unsigned int  jvar_num)
overrideprotectedvirtual

Definition at line 46 of file GeneralizedPlaneStrainOffDiagOSPD.C.

47{
48 if (jvar_num == _scalar_out_of_plane_strain_var_num)
49 {
50 prepare();
51
52 if (_var.number() == _disp_var[0]->number())
53 if (_use_full_jacobian)
54 computeDispFullOffDiagJacobianScalar(0, jvar_num);
55 else
56 computeDispPartialOffDiagJacobianScalar(0, jvar_num);
57 else if (_var.number() == _disp_var[1]->number())
58 if (_use_full_jacobian)
59 computeDispFullOffDiagJacobianScalar(1, jvar_num);
60 else
61 computeDispPartialOffDiagJacobianScalar(1, jvar_num);
62 else if (_temp_coupled ? _var.number() == _temp_var->number() : 0)
63 computeTempOffDiagJacobianScalar(jvar_num);
64 }
65}
GeneralizedPlaneStrainOffDiagOSPD::computeTempOffDiagJacobianScalar
void computeTempOffDiagJacobianScalar(unsigned int jvar_num)
Function to compute off disgonal Jacobian for coupling between temperature and scalar variable.
Definition GeneralizedPlaneStrainOffDiagOSPD.C:192
GeneralizedPlaneStrainOffDiagOSPD::computeDispPartialOffDiagJacobianScalar
void computeDispPartialOffDiagJacobianScalar(unsigned int component, unsigned int jvar_num)
Function to compute partial off diagonal Jacobian for coupling between displacements and scalar varia...
Definition GeneralizedPlaneStrainOffDiagOSPD.C:170
GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar
void computeDispFullOffDiagJacobianScalar(unsigned int component, unsigned int jvar_num)
Function to compute the full off diagonal Jacobian for coupling between displacements and scalar vari...
Definition GeneralizedPlaneStrainOffDiagOSPD.C:68

◆ computePDNonlocalOffDiagJacobian()

virtual void MechanicsBasePD::computePDNonlocalOffDiagJacobian ( unsigned int  ,
unsigned int   
)
inlinevirtualinherited

Function to compute nonlocal contribution to the off-diagonal Jacobian at the current nodes.

Parameters
jvar_numThe number of the first coupled variable
coupled_componentThe component number of the second coupled variable

Reimplemented in ForceStabilizedSmallStrainMechanicsNOSPD, HorizonStabilizedFormIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormIISmallStrainMechanicsNOSPD, HorizonStabilizedFormISmallStrainMechanicsNOSPD, MechanicsOSPD, and WeakPlaneStressNOSPD.

Definition at line 41 of file MechanicsBasePD.h.

43 {
44 }

Referenced by MechanicsBasePD::computeOffDiagJacobian().

◆ computeTempOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeTempOffDiagJacobianScalar ( unsigned int  jvar_num)
protected

Function to compute off disgonal Jacobian for coupling between temperature and scalar variable.

Parameters
jvar_numThe coupled scalar variable number

Definition at line 192 of file GeneralizedPlaneStrainOffDiagOSPD.C.

193{
194 // off-diagonal jacobian entries on the row
195 prepareMatrixTag(_assembly, jvar_num, _var.number(), _kne);
196
197 // calculate number of active neighbors for node i and j
198 std::vector<unsigned int> active_neighbors(_nnodes, 0);
199 for (unsigned int nd = 0; nd < _nnodes; nd++)
200 {
201 std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
202 for (unsigned int nb = 0; nb < bonds.size(); ++nb)
203 if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[nb])) > 0.5)
204 active_neighbors[nd]++;
205
206 if (active_neighbors[nd] == 0) // avoid dividing by zero
207 active_neighbors[nd] = 1;
208 }
209
210 // one-way coupling between the strain_zz and temperature. fill in the row corresponding to the
211 // scalar variable strain_zz
212 _kne(0, 0) += -_alpha[0] *
213 (_Cijkl[0](2, 2, 0, 0) + _Cijkl[0](2, 2, 1, 1) + _Cijkl[0](2, 2, 2, 2)) *
214 _node_vol[0] / active_neighbors[0] * _bond_status; // node i
215 _kne(0, 1) += -_alpha[0] *
216 (_Cijkl[0](2, 2, 0, 0) + _Cijkl[0](2, 2, 1, 1) + _Cijkl[0](2, 2, 2, 2)) *
217 _node_vol[1] / active_neighbors[1] * _bond_status; // node j
218
219 accumulateTaggedLocalMatrix(_assembly, jvar_num, _var.number(), _kne);
220}

Referenced by computeOffDiagJacobianScalar().

◆ initialSetup()

void MechanicsBasePD::initialSetup ( )
overridevirtualinherited

Definition at line 45 of file MechanicsBasePD.C.

46{
47 _orientation = &_assembly.getFE(FEType(), 1)->get_dxyzdxi();
48}
MechanicsBasePD::_orientation
const std::vector< RealGradient > * _orientation
Vector of bond in current configuration.
Definition MechanicsBasePD.h:67

◆ prepare()

void MechanicsBasePD::prepare ( )
overridevirtualinherited

Definition at line 51 of file MechanicsBasePD.C.

52{
53 PeridynamicsKernelBase::prepare();
54
55 _ivardofs.resize(_nnodes);
56 _weights.resize(_nnodes);
57 for (unsigned int nd = 0; nd < _nnodes; ++nd)
58 {
59 _ivardofs[nd] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
60 _weights[nd] = _pdmesh.getNeighborWeight(
61 _current_elem->node_id(nd),
62 _pdmesh.getNeighborIndex(_current_elem->node_id(nd), _current_elem->node_id(1 - nd)));
63 }
64
65 for (unsigned int i = 0; i < _dim; ++i)
66 _current_vec(i) = _origin_vec(i) + _disp_var[i]->getNodalValue(*_current_elem->node_ptr(1)) -
67 _disp_var[i]->getNodalValue(*_current_elem->node_ptr(0));
68
69 _current_unit_vec = _current_vec / _current_vec.norm();
70}
MechanicsBasePD::_current_vec
RealGradient _current_vec
Vector of bond in current configuration.
Definition MechanicsBasePD.h:76
MechanicsBasePD::_weights
std::vector< Real > _weights
weights used for the current element to obtain the nodal stress
Definition MechanicsBasePD.h:73
PeridynamicsKernelBase::prepare
virtual void prepare()
Function to precalculate data which will be used in the derived classes.
Definition PeridynamicsKernelBase.C:43

Referenced by MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), and computeOffDiagJacobianScalar().

◆ validParams()

InputParameters GeneralizedPlaneStrainOffDiagOSPD::validParams ( )
static

Definition at line 17 of file GeneralizedPlaneStrainOffDiagOSPD.C.

18{
19 InputParameters params = MechanicsBasePD::validParams();
20 params.addClassDescription(
21 "Class for calculating the off-diagonal Jacobian corresponding to "
22 "coupling between displacements (or temperature) and the scalar out-of-plane strain for "
23 "the generalized plane strain using the OSPD formulation");
24
25 params.addCoupledVar("scalar_out_of_plane_strain",
26 "Scalar variable for strain in the out-of-plane direction");
27
28 return params;
29}
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)
InputParameters::addCoupledVar
void addCoupledVar(const std::string &name, const std::string &doc_string)
MechanicsBasePD::validParams
static InputParameters validParams()
Definition MechanicsBasePD.C:13

Member Data Documentation

◆ _alpha

const MaterialProperty<Real>& GeneralizedPlaneStrainOffDiagOSPD::_alpha
protected

Definition at line 54 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeTempOffDiagJacobianScalar().

◆ _bond_local_dfdE

const MaterialProperty<Real>& GeneralizedPlaneStrainOffDiagOSPD::_bond_local_dfdE
protected

Bond based material properties.

Definition at line 52 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), and computeDispPartialOffDiagJacobianScalar().

◆ _bond_nonlocal_dfdE

const MaterialProperty<Real>& GeneralizedPlaneStrainOffDiagOSPD::_bond_nonlocal_dfdE
protected

Definition at line 53 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar().

◆ _Cijkl

const MaterialProperty<RankFourTensor>& GeneralizedPlaneStrainOffDiagOSPD::_Cijkl
protected

Material point based material property.

Definition at line 58 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), and computeTempOffDiagJacobianScalar().

◆ _current_unit_vec

RealGradient MechanicsBasePD::_current_unit_vec
protectedinherited

Unit vector of bond in current configuration.

Definition at line 79 of file MechanicsBasePD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), computeDispPartialOffDiagJacobianScalar(), MechanicsBPD::computeLocalJacobian(), MechanicsOSPD::computeLocalJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), MechanicsBPD::computeLocalResidual(), MechanicsOSPD::computeLocalResidual(), MechanicsOSPD::computeNonlocalJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

◆ _current_vec

RealGradient MechanicsBasePD::_current_vec
protectedinherited

Vector of bond in current configuration.

Definition at line 76 of file MechanicsBasePD.h.

Referenced by MechanicsBPD::computeLocalJacobian(), MechanicsOSPD::computeLocalJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), and MechanicsBasePD::prepare().

◆ _disp_var

std::vector<MooseVariable *> MechanicsBasePD::_disp_var
protectedinherited

displacement variables

Definition at line 51 of file MechanicsBasePD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalResidual(), MechanicsOSPD::computeNonlocalJacobian(), MechanicsOSPD::computeNonlocalResidual(), MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), computeOffDiagJacobianScalar(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::GeneralizedPlaneStrainOffDiagNOSPD(), GeneralizedPlaneStrainOffDiagOSPD(), MechanicsBasePD::MechanicsBasePD(), and MechanicsBasePD::prepare().

◆ _ivardofs

std::vector<dof_id_type> MechanicsBasePD::_ivardofs
protectedinherited

Current variable dof numbers for nodes i and j.

Definition at line 70 of file MechanicsBasePD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeNonlocalJacobian(), MechanicsOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalResidual(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeNonlocalResidual(), MechanicsOSPD::computeNonlocalResidual(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

◆ _ken

DenseMatrix<Number> GeneralizedPlaneStrainOffDiagOSPD::_ken
protected

Definition at line 63 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), and computeDispPartialOffDiagJacobianScalar().

◆ _kne

DenseMatrix<Number> GeneralizedPlaneStrainOffDiagOSPD::_kne
protected

Definition at line 64 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), computeDispPartialOffDiagJacobianScalar(), and computeTempOffDiagJacobianScalar().

◆ _ndisp

unsigned int MechanicsBasePD::_ndisp
protectedinherited

number of displacement components

Definition at line 59 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::MechanicsBasePD().

◆ _orientation

const std::vector<RealGradient>* MechanicsBasePD::_orientation
protectedinherited

Vector of bond in current configuration.

Definition at line 67 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::initialSetup().

◆ _out_of_plane_strain_coupled

const bool MechanicsBasePD::_out_of_plane_strain_coupled
protectedinherited

Parameters for out-of-plane strain in weak plane stress formulation.

Definition at line 62 of file MechanicsBasePD.h.

Referenced by HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

◆ _out_of_plane_strain_var

MooseVariable* MechanicsBasePD::_out_of_plane_strain_var
protectedinherited

Definition at line 63 of file MechanicsBasePD.h.

Referenced by HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

◆ _scalar_out_of_plane_strain_var_num

const unsigned int GeneralizedPlaneStrainOffDiagOSPD::_scalar_out_of_plane_strain_var_num
protected

The variable number of the scalar out-of-plane strain variable.

Definition at line 61 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeOffDiagJacobianScalar().

◆ _temp_coupled

const bool MechanicsBasePD::_temp_coupled
protectedinherited

Temperature variable.

Definition at line 54 of file MechanicsBasePD.h.

Referenced by ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), computeOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBaseNOSPD::MechanicsBaseNOSPD().

◆ _temp_var

MooseVariable* MechanicsBasePD::_temp_var
protectedinherited

Definition at line 55 of file MechanicsBasePD.h.

Referenced by ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), computeOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBaseNOSPD::MechanicsBaseNOSPD().

◆ _weights

std::vector<Real> MechanicsBasePD::_weights
protectedinherited

weights used for the current element to obtain the nodal stress

Definition at line 73 of file MechanicsBasePD.h.

Referenced by WeakPlaneStressNOSPD::computeLocalJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalResidual(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

The documentation for this class was generated from the following files:
Generated on Fri Aug 21 2026 13:42:50 for https://mooseframework.inl.gov by doxygen 1.9.8