https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
HorizonStabilizedFormIISmallStrainMechanicsNOSPD Class Reference

Kernel class for Form II of the horizon-associated peridynamic correspondence material model for small strain. More...

#include <HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h>

Inheritance diagram for HorizonStabilizedFormIISmallStrainMechanicsNOSPD:
[legend]

Public Types

typedef DerivativeMaterialPropertyNameInterface::SymbolName SymbolName
 

Public Member Functions

 HorizonStabilizedFormIISmallStrainMechanicsNOSPD (const InputParameters &parameters)
 
virtual void computeOffDiagJacobian (unsigned int jvar) override
 
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 computeNonlocalResidual () override
 
virtual void computeLocalJacobian () override
 
virtual void computeNonlocalJacobian () override
 
virtual void computeLocalOffDiagJacobian (unsigned int, unsigned int coupled_component) override
 Function to compute local contribution to the off-diagonal Jacobian at the current nodes. More...
 
virtual void computePDNonlocalOffDiagJacobian (unsigned int jvar_num, unsigned int coupled_component) override
 Function to compute nonlocal contribution to the off-diagonal Jacobian at the current nodes. More...
 
virtual RankTwoTensor computeDSDU (unsigned int component, unsigned int nd)
 Function to compute derivative of stress with respect to displacements for small strain problems. More...
 

Protected Attributes

const unsigned int _component
 The index of displacement component. More...
 
std::vector< MooseVariable * > _disp_var
 displacement variables More...
 
unsigned int _ndisp
 number of displacement components More...
 
const std::vector< RealGradient > * _orientation
 Vector of bond in current configuration. More...
 
std::vector< dof_id_type > _ivardofs
 Current variable dof numbers for nodes i and j. More...
 
std::vector< Real > _weights
 weights used for the current element to obtain the nodal stress More...
 
RealGradient _current_vec
 Vector of bond in current configuration. More...
 
RealGradient _current_unit_vec
 Unit vector of bond in current configuration. More...
 
const MaterialProperty< Real > & _multi
 Material point based material properties. More...
 
const MaterialProperty< RankTwoTensor > & _stress
 
const MaterialProperty< RankTwoTensor > & _shape2
 
const MaterialProperty< RankTwoTensor > & _dgrad
 
const MaterialProperty< RankTwoTensor > & _ddgraddu
 
const MaterialProperty< RankTwoTensor > & _ddgraddv
 
const MaterialProperty< RankTwoTensor > & _ddgraddw
 
const MaterialProperty< RankFourTensor > & _Jacobian_mult
 
const std::vector< MaterialPropertyName > _eigenstrain_names
 
std::vector< const MaterialProperty< RankTwoTensor > * > _deigenstrain_dT
 
const bool _temp_coupled
 Temperature variable. More...
 
MooseVariable * _temp_var
 
const bool _out_of_plane_strain_coupled
 Parameters for out-of-plane strain in weak plane stress formulation. More...
 
MooseVariable * _out_of_plane_strain_var
 

Detailed Description

Kernel class for Form II of the horizon-associated peridynamic correspondence material model for small strain.

Definition at line 20 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h.

Constructor & Destructor Documentation

◆ HorizonStabilizedFormIISmallStrainMechanicsNOSPD()

HorizonStabilizedFormIISmallStrainMechanicsNOSPD::HorizonStabilizedFormIISmallStrainMechanicsNOSPD ( const InputParameters &  parameters)

Definition at line 30 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

32  : MechanicsBaseNOSPD(parameters), _component(getParam<unsigned int>("component"))
33 {
34 }
HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
const unsigned int _component
The index of displacement component.
Definition: HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h:40
MechanicsBaseNOSPD::MechanicsBaseNOSPD
MechanicsBaseNOSPD(const InputParameters &parameters)
Definition: MechanicsBaseNOSPD.C:27

Member Function Documentation

◆ computeDSDU()

RankTwoTensor MechanicsBaseNOSPD::computeDSDU ( unsigned int  component,
unsigned int  nd 
)
protectedvirtualinherited

Function to compute derivative of stress with respect to displacements for small strain problems.

Parameters
componentThe index of displacement component
ndThe local index of element node (either 1 or 2 for Edge2 element)
Returns
The calculated derivative

Reimplemented in MechanicsFiniteStrainBaseNOSPD.

Definition at line 47 of file MechanicsBaseNOSPD.C.

Referenced by GeneralizedPlaneStrainOffDiagNOSPD::computeDispFullOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagNOSPD::computeDispPartialOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalJacobian(), computeLocalJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), and computeLocalOffDiagJacobian().

48 {
49  // compute the derivative of stress w.r.t the solution components for small strain
50  RankTwoTensor dSdU;
51  if (component == 0)
52  dSdU = _Jacobian_mult[nd] * 0.5 * (_ddgraddu[nd].transpose() + _ddgraddu[nd]);
53  else if (component == 1)
54  dSdU = _Jacobian_mult[nd] * 0.5 * (_ddgraddv[nd].transpose() + _ddgraddv[nd]);
55  else if (component == 2)
56  dSdU = _Jacobian_mult[nd] * 0.5 * (_ddgraddw[nd].transpose() + _ddgraddw[nd]);
57 
58  return dSdU;
59 }
MechanicsBaseNOSPD::_ddgraddu
const MaterialProperty< RankTwoTensor > & _ddgraddu
Definition: MechanicsBaseNOSPD.h:39
MechanicsBaseNOSPD::_Jacobian_mult
const MaterialProperty< RankFourTensor > & _Jacobian_mult
Definition: MechanicsBaseNOSPD.h:42
NS::component
static const std::string component
Definition: NS.h:153
MechanicsBaseNOSPD::_ddgraddw
const MaterialProperty< RankTwoTensor > & _ddgraddw
Definition: MechanicsBaseNOSPD.h:41
MechanicsBaseNOSPD::_ddgraddv
const MaterialProperty< RankTwoTensor > & _ddgraddv
Definition: MechanicsBaseNOSPD.h:40

◆ computeLocalJacobian()

void HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalJacobian ( )
overrideprotectedvirtual

Definition at line 111 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

112 {
113  for (unsigned int nd = 0; nd < _nnodes; ++nd)
114  {
115  std::vector<dof_id_type> ivardofs(_nnodes);
116  ivardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
117  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
118  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
119 
120  dof_id_type nb_index =
121  std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
122  neighbors.begin();
123  std::vector<dof_id_type> dg_neighbors =
124  _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
125 
126  RealGradient origin_vec_nb;
127  Real node_vol_nb;
128 
129  for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
130  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
131  {
132  ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
133  ->dof_number(_sys.number(), _var.number(), 0);
134  origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
135  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
136  node_vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
137 
138  for (unsigned int i = 0; i < _nnodes; ++i)
139  for (unsigned int j = 0; j < _nnodes; ++j)
140  _local_ke(i, j) =
141  (i == 0 ? -1 : 1) * (j == 0 ? 1 : 0) * _multi[nd] * _horizon_radius[nd] /
142  origin_vec_nb.norm() *
143  (computeDSDU(_component, nd) * _shape2[nd].inverse()).row(_component) *
144  origin_vec_nb * node_vol_nb * _bond_status;
145 
146  addJacobian(_assembly, _local_ke, ivardofs, ivardofs, _var.scalingFactor());
147  }
148  }
149  _local_ke.zero();
150 }
MechanicsBaseNOSPD::computeDSDU
virtual RankTwoTensor computeDSDU(unsigned int component, unsigned int nd)
Function to compute derivative of stress with respect to displacements for small strain problems...
Definition: MechanicsBaseNOSPD.C:47
VectorValue< Real >::norm
auto norm() const -> decltype(std::norm(Real()))
HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
const unsigned int _component
The index of displacement component.
Definition: HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h:40
MechanicsBaseNOSPD::_shape2
const MaterialProperty< RankTwoTensor > & _shape2
Definition: MechanicsBaseNOSPD.h:37
MechanicsBaseNOSPD::_multi
const MaterialProperty< Real > & _multi
Material point based material properties.
Definition: MechanicsBaseNOSPD.h:35
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
dof_id_type
uint8_t dof_id_type

◆ computeLocalOffDiagJacobian()

void HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian ( unsigned int  ,
unsigned int   
)
overrideprotectedvirtual

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

Parameters
coupled_componentThe coupled variable number

Reimplemented from MechanicsBasePD.

Definition at line 241 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

243 {
244  if (_temp_coupled && jvar_num == _temp_var->number()) // temperature is coupled
245  {
246  std::vector<RankTwoTensor> dSdT(_nnodes);
247  for (unsigned int nd = 0; nd < _nnodes; ++nd)
248  for (unsigned int es = 0; es < _deigenstrain_dT.size(); ++es)
249  dSdT[nd] = -_Jacobian_mult[nd] * (*_deigenstrain_dT[es])[nd];
250 
251  for (unsigned int nd = 0; nd < _nnodes; ++nd)
252  {
253  std::vector<dof_id_type> ivardofs(_nnodes), jvardofs(_nnodes);
254  ivardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
255  jvardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), jvar_num, 0);
256  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
257  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
258 
259  dof_id_type nb_index =
260  std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
261  neighbors.begin();
262  std::vector<dof_id_type> dg_neighbors =
263  _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
264 
265  RealGradient origin_vec_nb;
266  Real node_vol_nb;
267 
268  for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
269  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
270  {
271  ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
272  ->dof_number(_sys.number(), _var.number(), 0);
273  jvardofs[1] =
274  _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])->dof_number(_sys.number(), jvar_num, 0);
275  origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
276  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
277  node_vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
278 
279  for (unsigned int i = 0; i < _nnodes; ++i)
280  for (unsigned int j = 0; j < _nnodes; ++j)
281  _local_ke(i, j) = (i == 0 ? -1 : 1) * (j == 0 ? 1 : 0) * _multi[nd] *
282  _horizon_radius[nd] / origin_vec_nb.norm() *
283  (dSdT[nd] * _shape2[nd].inverse()).row(_component) * origin_vec_nb *
284  node_vol_nb * _bond_status;
285 
286  addJacobian(_assembly, _local_ke, ivardofs, jvardofs, _var.scalingFactor());
287  }
288  }
289  _local_ke.zero();
290  }
291  else if (_out_of_plane_strain_coupled &&
292  jvar_num == _out_of_plane_strain_var
293  ->number()) // weak plane stress case, out_of_plane_strain is coupled
294  {
295  std::vector<RankTwoTensor> dSdE33(_nnodes);
296  for (unsigned int nd = 0; nd < _nnodes; ++nd)
297  for (unsigned int i = 0; i < 3; ++i)
298  for (unsigned int j = 0; j < 3; ++j)
299  dSdE33[nd](i, j) = _Jacobian_mult[nd](i, j, 2, 2);
300 
301  for (unsigned int nd = 0; nd < _nnodes; ++nd)
302  {
303  std::vector<dof_id_type> ivardofs(_nnodes), jvardofs(_nnodes);
304  ivardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
305  jvardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), jvar_num, 0);
306  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
307  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
308 
309  dof_id_type nb_index =
310  std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
311  neighbors.begin();
312  std::vector<dof_id_type> dg_neighbors =
313  _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
314 
315  RealGradient origin_vec_nb;
316  Real node_vol_nb;
317 
318  for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
319  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
320  {
321  ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
322  ->dof_number(_sys.number(), _var.number(), 0);
323  jvardofs[1] =
324  _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])->dof_number(_sys.number(), jvar_num, 0);
325  origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
326  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
327  node_vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
328 
329  for (unsigned int i = 0; i < _nnodes; ++i)
330  for (unsigned int j = 0; j < _nnodes; ++j)
331  _local_ke(i, j) = (i == 0 ? -1 : 1) * (j == 0 ? 1 : 0) * _multi[nd] *
332  _horizon_radius[nd] / origin_vec_nb.norm() *
333  (dSdE33[nd] * _shape2[nd].inverse()).row(_component) *
334  origin_vec_nb * node_vol_nb * _bond_status;
335 
336  addJacobian(_assembly, _local_ke, ivardofs, jvardofs, _var.scalingFactor());
337  }
338  }
339  _local_ke.zero();
340  }
341  else // off-diagonal Jacobian with respect to other displacement variables
342  {
343  for (unsigned int nd = 0; nd < _nnodes; ++nd)
344  {
345  std::vector<dof_id_type> ivardofs(_nnodes), jvardofs(_nnodes);
346  ivardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
347  jvardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), jvar_num, 0);
348  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
349  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
350 
351  dof_id_type nb_index =
352  std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
353  neighbors.begin();
354  std::vector<dof_id_type> dg_neighbors =
355  _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
356 
357  RealGradient origin_vec_nb;
358  Real node_vol_nb;
359 
360  for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
361  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
362  {
363  ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
364  ->dof_number(_sys.number(), _var.number(), 0);
365  jvardofs[1] =
366  _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])->dof_number(_sys.number(), jvar_num, 0);
367  origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
368  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
369  node_vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
370 
371  for (unsigned int i = 0; i < _nnodes; ++i)
372  for (unsigned int j = 0; j < _nnodes; ++j)
373  _local_ke(i, j) =
374  (i == 0 ? -1 : 1) * (j == 0 ? 1 : 0) * _multi[nd] * _horizon_radius[nd] /
375  origin_vec_nb.norm() *
376  (computeDSDU(coupled_component, nd) * _shape2[nd].inverse()).row(_component) *
377  origin_vec_nb * node_vol_nb * _bond_status;
378 
379  addJacobian(_assembly, _local_ke, ivardofs, jvardofs, _var.scalingFactor());
380  }
381  }
382  _local_ke.zero();
383  }
384 }
MechanicsBasePD::_temp_coupled
const bool _temp_coupled
Temperature variable.
Definition: MechanicsBasePD.h:50
MechanicsBaseNOSPD::computeDSDU
virtual RankTwoTensor computeDSDU(unsigned int component, unsigned int nd)
Function to compute derivative of stress with respect to displacements for small strain problems...
Definition: MechanicsBaseNOSPD.C:47
VectorValue< Real >::norm
auto norm() const -> decltype(std::norm(Real()))
MechanicsBaseNOSPD::_Jacobian_mult
const MaterialProperty< RankFourTensor > & _Jacobian_mult
Definition: MechanicsBaseNOSPD.h:42
MooseVariableFE< Real >::number
unsigned int number() const
HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
const unsigned int _component
The index of displacement component.
Definition: HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h:40
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:58
MechanicsBaseNOSPD::_deigenstrain_dT
std::vector< const MaterialProperty< RankTwoTensor > * > _deigenstrain_dT
Definition: MechanicsBaseNOSPD.h:44
MechanicsBasePD::_out_of_plane_strain_var
MooseVariable * _out_of_plane_strain_var
Definition: MechanicsBasePD.h:59
MechanicsBaseNOSPD::_shape2
const MaterialProperty< RankTwoTensor > & _shape2
Definition: MechanicsBaseNOSPD.h:37
MechanicsBaseNOSPD::_multi
const MaterialProperty< Real > & _multi
Material point based material properties.
Definition: MechanicsBaseNOSPD.h:35
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
MechanicsBasePD::_temp_var
MooseVariable * _temp_var
Definition: MechanicsBasePD.h:51
dof_id_type
uint8_t dof_id_type

◆ computeLocalResidual()

void HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalResidual ( )
overrideprotectedvirtual

Definition at line 37 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

38 {
39  // For small strain assumptions, stress measures, i.e., Cauchy stress (Sigma), the first
40  // Piola-Kirchhoff stress (P), and the second Piola-Kirchhoff stress (S) are approximately the
41  // same. Thus, the nodal force state tensors are calculated using the Cauchy stresses,
42  // i.e., T = Sigma * inv(Shape) * xi * multi.
43  // Cauchy stress is calculated as Sigma = C * E in the ComputeSmallStrainNOSPD material class.
44 
45  for (unsigned int nd = 0; nd < _nnodes; ++nd)
46  for (unsigned int i = 0; i < _nnodes; ++i)
47  _local_re(i) += (i == 0 ? -1 : 1) * _multi[nd] * _horizon_radius[nd] / _origin_vec.norm() *
48  (_stress[nd] * _shape2[nd].inverse()).row(_component) * _origin_vec *
49  _node_vol[1 - nd] * _bond_status;
50 }
HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
const unsigned int _component
The index of displacement component.
Definition: HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h:40
MechanicsBaseNOSPD::_stress
const MaterialProperty< RankTwoTensor > & _stress
Definition: MechanicsBaseNOSPD.h:36
MechanicsBaseNOSPD::_shape2
const MaterialProperty< RankTwoTensor > & _shape2
Definition: MechanicsBaseNOSPD.h:37
MechanicsBaseNOSPD::_multi
const MaterialProperty< Real > & _multi
Material point based material properties.
Definition: MechanicsBaseNOSPD.h:35

◆ computeNonlocalJacobian()

void HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeNonlocalJacobian ( )
overrideprotectedvirtual

Definition at line 153 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

154 {
155  // includes dTi/dUj and dTj/dUi contributions
156  for (unsigned int nd = 0; nd < _nnodes; ++nd)
157  {
158  // calculation of jacobian contribution to current_node's neighbors
159  std::vector<dof_id_type> ivardofs(_nnodes), jvardofs(_nnodes);
160  ivardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
161  jvardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
162  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
163  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
164 
165  dof_id_type nb_index =
166  std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
167  neighbors.begin();
168  std::vector<dof_id_type> dg_neighbors =
169  _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
170 
171  RealGradient origin_vec_nb1;
172  Real node_vol_nb1;
173 
174  for (unsigned int nb1 = 0; nb1 < dg_neighbors.size(); ++nb1)
175  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb1]])) > 0.5)
176  {
177  ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb1]])
178  ->dof_number(_sys.number(), _var.number(), 0);
179  origin_vec_nb1 = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb1]]) -
180  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
181  node_vol_nb1 = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb1]]);
182 
183  Real vol_nb2;
184  RealGradient origin_vec_nb2;
185  RankTwoTensor dFdUk, dPxdUkx;
186 
187  for (unsigned int nb2 = 0; nb2 < dg_neighbors.size(); ++nb2)
188  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb2]])) > 0.5)
189  {
190  jvardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb2]])
191  ->dof_number(_sys.number(), _var.number(), 0);
192  vol_nb2 = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb2]]);
193 
194  origin_vec_nb2 = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb2]]) -
195  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
196 
197  dFdUk.zero();
198  for (unsigned int i = 0; i < _dim; ++i)
199  dFdUk(_component, i) =
200  _horizon_radius[nd] / origin_vec_nb2.norm() * origin_vec_nb2(i) * vol_nb2;
201 
202  dFdUk *= _shape2[nd].inverse();
203  dPxdUkx = _Jacobian_mult[nd] * 0.5 * (dFdUk.transpose() + dFdUk);
204 
205  for (unsigned int i = 0; i < _nnodes; ++i)
206  for (unsigned int j = 0; j < _nnodes; ++j)
207  _local_ke(i, j) = (i == 0 ? -1 : 1) * (j == 0 ? 0 : 1) * _multi[nd] *
208  _horizon_radius[nd] / origin_vec_nb1.norm() *
209  (dPxdUkx * _shape2[nd].inverse()).row(_component) *
210  origin_vec_nb1 * node_vol_nb1 * _bond_status;
211 
212  addJacobian(_assembly, _local_ke, ivardofs, jvardofs, _var.scalingFactor());
213 
214  if (_has_diag_save_in)
215  {
216  unsigned int rows = _nnodes;
217  DenseVector<Real> diag(rows);
218  for (unsigned int i = 0; i < rows; ++i)
219  diag(i) = _local_ke(i, i);
220 
221  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
222  for (unsigned int i = 0; i < _diag_save_in.size(); ++i)
223  {
224  std::vector<dof_id_type> diag_save_in_dofs(2);
225  diag_save_in_dofs[0] = _current_elem->node_ptr(nd)->dof_number(
226  _diag_save_in[i]->sys().number(), _diag_save_in[i]->number(), 0);
227  diag_save_in_dofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb2]])
228  ->dof_number(_diag_save_in[i]->sys().number(),
229  _diag_save_in[i]->number(),
230  0);
231 
232  _diag_save_in[i]->sys().solution().add_vector(diag, diag_save_in_dofs);
233  }
234  }
235  }
236  }
237  }
238 }
VectorValue< Real >::norm
auto norm() const -> decltype(std::norm(Real()))
MechanicsBaseNOSPD::_Jacobian_mult
const MaterialProperty< RankFourTensor > & _Jacobian_mult
Definition: MechanicsBaseNOSPD.h:42
HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
const unsigned int _component
The index of displacement component.
Definition: HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h:40
MechanicsBaseNOSPD::_shape2
const MaterialProperty< RankTwoTensor > & _shape2
Definition: MechanicsBaseNOSPD.h:37
MechanicsBaseNOSPD::_multi
const MaterialProperty< Real > & _multi
Material point based material properties.
Definition: MechanicsBaseNOSPD.h:35
RankTwoTensorTempl< Real >::transpose
RankTwoTensorTempl< Real > transpose() const
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
dof_id_type
uint8_t dof_id_type

◆ computeNonlocalResidual()

void HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeNonlocalResidual ( )
overrideprotectedvirtual

Definition at line 53 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

54 {
55  for (unsigned int nd = 0; nd < _nnodes; ++nd)
56  {
57  // calculation of residual contribution to current_node's neighbors
58  std::vector<dof_id_type> ivardofs(_nnodes);
59  ivardofs[0] = _ivardofs[nd];
60  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
61  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
62 
63  dof_id_type nb_index =
64  std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
65  neighbors.begin();
66  std::vector<dof_id_type> dg_neighbors =
67  _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
68 
69  RealGradient origin_vec_nb;
70  Real node_vol_nb;
71 
72  for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
73  if (neighbors[dg_neighbors[nb]] != _current_elem->node_id(1 - nd) &&
74  _bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
75  {
76  ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
77  ->dof_number(_sys.number(), _var.number(), 0);
78  origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
79  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
80  node_vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
81 
82  for (unsigned int i = 0; i < _nnodes; ++i)
83  _local_re(i) = (i == 0 ? -1 : 1) * _multi[nd] * _horizon_radius[nd] /
84  origin_vec_nb.norm() *
85  (_stress[nd] * _shape2[nd].inverse()).row(_component) * origin_vec_nb *
86  node_vol_nb * _bond_status;
87 
88  // cache the residual contribution
89  addResiduals(_assembly, _local_re, ivardofs, _var.scalingFactor());
90 
91  if (_has_save_in)
92  {
93  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
94  for (unsigned int i = 0; i < _save_in.size(); ++i)
95  {
96  std::vector<dof_id_type> save_in_dofs(_nnodes);
97  save_in_dofs[0] = _current_elem->node_ptr(nd)->dof_number(
98  _save_in[i]->sys().number(), _save_in[i]->number(), 0);
99  save_in_dofs[1] =
100  _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
101  ->dof_number(_save_in[i]->sys().number(), _save_in[i]->number(), 0);
102 
103  _save_in[i]->sys().solution().add_vector(_local_re, save_in_dofs);
104  }
105  }
106  }
107  }
108 }
VectorValue< Real >::norm
auto norm() const -> decltype(std::norm(Real()))
HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
const unsigned int _component
The index of displacement component.
Definition: HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h:40
MechanicsBaseNOSPD::_stress
const MaterialProperty< RankTwoTensor > & _stress
Definition: MechanicsBaseNOSPD.h:36
MechanicsBasePD::_ivardofs
std::vector< dof_id_type > _ivardofs
Current variable dof numbers for nodes i and j.
Definition: MechanicsBasePD.h:66
MechanicsBaseNOSPD::_shape2
const MaterialProperty< RankTwoTensor > & _shape2
Definition: MechanicsBaseNOSPD.h:37
MechanicsBaseNOSPD::_multi
const MaterialProperty< Real > & _multi
Material point based material properties.
Definition: MechanicsBaseNOSPD.h:35
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
dof_id_type
uint8_t dof_id_type

◆ 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::_temp_coupled
const bool _temp_coupled
Temperature variable.
Definition: MechanicsBasePD.h:50
MooseVariableFE< Real >::number
unsigned int number() const
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::_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:58
MechanicsBasePD::_disp_var
std::vector< MooseVariable * > _disp_var
displacement variables
Definition: MechanicsBasePD.h:47
MechanicsBasePD::prepare
virtual void prepare() override
Definition: MechanicsBasePD.C:51
MechanicsBasePD::_out_of_plane_strain_var
MooseVariable * _out_of_plane_strain_var
Definition: MechanicsBasePD.h:59
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:39
MechanicsBasePD::_temp_var
MooseVariable * _temp_var
Definition: MechanicsBasePD.h:51

◆ computePDNonlocalOffDiagJacobian()

void HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian ( unsigned int  ,
unsigned int   
)
overrideprotectedvirtual

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 from MechanicsBasePD.

Definition at line 387 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

389 {
390  if (_temp_coupled && jvar_num == _temp_var->number())
391  {
392  // no nonlocal contribution from temperature
393  }
394  else if (_out_of_plane_strain_coupled && jvar_num == _out_of_plane_strain_var->number())
395  {
396  // no nonlocal contribution from out of plane strain
397  }
398  else
399  {
400  for (unsigned int nd = 0; nd < _nnodes; ++nd)
401  {
402  // calculation of jacobian contribution to current_node's neighbors
403  // NOT including the contribution to nodes i and j, which is considered as local off-diagonal
404  std::vector<dof_id_type> ivardofs(_nnodes), jvardofs(_nnodes);
405  ivardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
406  jvardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), jvar_num, 0);
407  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
408  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
409 
410  dof_id_type nb_index =
411  std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
412  neighbors.begin();
413  std::vector<dof_id_type> dg_neighbors =
414  _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
415 
416  RealGradient origin_vec_nb1;
417  Real node_vol_nb1;
418 
419  for (unsigned int nb1 = 0; nb1 < dg_neighbors.size(); ++nb1)
420  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb1]])) > 0.5)
421  {
422  ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb1]])
423  ->dof_number(_sys.number(), _var.number(), 0);
424  origin_vec_nb1 = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb1]]) -
425  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
426  node_vol_nb1 = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb1]]);
427 
428  Real vol_nb2;
429  RealGradient origin_vec_nb2;
430  RankTwoTensor dFdUk, dPxdUky;
431 
432  for (unsigned int nb2 = 0; nb2 < dg_neighbors.size(); ++nb2)
433  if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb2]])) >
434  0.5)
435  {
436  jvardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb2]])
437  ->dof_number(_sys.number(), jvar_num, 0);
438  vol_nb2 = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb2]]);
439 
440  origin_vec_nb2 = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb2]]) -
441  _pdmesh.getNodeCoord(_current_elem->node_id(nd));
442 
443  dFdUk.zero();
444  for (unsigned int i = 0; i < _dim; ++i)
445  dFdUk(coupled_component, i) =
446  _horizon_radius[nd] / origin_vec_nb2.norm() * origin_vec_nb2(i) * vol_nb2;
447 
448  dFdUk *= _shape2[nd].inverse();
449  dPxdUky = _Jacobian_mult[nd] * 0.5 * (dFdUk.transpose() + dFdUk);
450 
451  for (unsigned int i = 0; i < _nnodes; ++i)
452  for (unsigned int j = 0; j < _nnodes; ++j)
453  _local_ke(i, j) = (i == 0 ? -1 : 1) * (j == 0 ? 0 : 1) * _multi[nd] *
454  _horizon_radius[nd] / origin_vec_nb1.norm() *
455  (dPxdUky * _shape2[nd].inverse()).row(_component) *
456  origin_vec_nb1 * node_vol_nb1 * _bond_status;
457 
458  addJacobian(_assembly, _local_ke, ivardofs, jvardofs, _var.scalingFactor());
459  }
460  }
461  }
462  }
463 }
MechanicsBasePD::_temp_coupled
const bool _temp_coupled
Temperature variable.
Definition: MechanicsBasePD.h:50
VectorValue< Real >::norm
auto norm() const -> decltype(std::norm(Real()))
MechanicsBaseNOSPD::_Jacobian_mult
const MaterialProperty< RankFourTensor > & _Jacobian_mult
Definition: MechanicsBaseNOSPD.h:42
MooseVariableFE< Real >::number
unsigned int number() const
HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
const unsigned int _component
The index of displacement component.
Definition: HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h:40
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:58
MechanicsBasePD::_out_of_plane_strain_var
MooseVariable * _out_of_plane_strain_var
Definition: MechanicsBasePD.h:59
MechanicsBaseNOSPD::_shape2
const MaterialProperty< RankTwoTensor > & _shape2
Definition: MechanicsBaseNOSPD.h:37
MechanicsBaseNOSPD::_multi
const MaterialProperty< Real > & _multi
Material point based material properties.
Definition: MechanicsBaseNOSPD.h:35
RankTwoTensorTempl< Real >::transpose
RankTwoTensorTempl< Real > transpose() const
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
MechanicsBasePD::_temp_var
MooseVariable * _temp_var
Definition: MechanicsBasePD.h:51
dof_id_type
uint8_t dof_id_type

◆ 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:63

◆ prepare()

void MechanicsBasePD::prepare ( )
overridevirtualinherited

Definition at line 51 of file MechanicsBasePD.C.

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

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 }
VectorValue< Real >::norm
auto norm() const -> decltype(std::norm(Real()))
MechanicsBasePD::_current_unit_vec
RealGradient _current_unit_vec
Unit vector of bond in current configuration.
Definition: MechanicsBasePD.h:75
MechanicsBasePD::_weights
std::vector< Real > _weights
weights used for the current element to obtain the nodal stress
Definition: MechanicsBasePD.h:69
MechanicsBasePD::_disp_var
std::vector< MooseVariable * > _disp_var
displacement variables
Definition: MechanicsBasePD.h:47
MechanicsBasePD::_current_vec
RealGradient _current_vec
Vector of bond in current configuration.
Definition: MechanicsBasePD.h:72
MechanicsBasePD::_ivardofs
std::vector< dof_id_type > _ivardofs
Current variable dof numbers for nodes i and j.
Definition: MechanicsBasePD.h:66
PeridynamicsKernelBase::prepare
virtual void prepare()
Function to precalculate data which will be used in the derived classes.
Definition: PeridynamicsKernelBase.C:43

◆ validParams()

InputParameters HorizonStabilizedFormIISmallStrainMechanicsNOSPD::validParams ( )
static

Definition at line 16 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.C.

17 {
18  InputParameters params = MechanicsBaseNOSPD::validParams();
19  params.addClassDescription(
20  "Class for calculating the residual and the Jacobian for Form II "
21  "of the horizon-stabilized peridynamic correspondence model under small strain assumptions");
22 
23  params.addRequiredParam<unsigned int>(
24  "component",
25  "An integer corresponding to the variable this kernel acts on (0 for x, 1 for y, 2 for z)");
26 
27  return params;
28 }
InputParameters::addRequiredParam
void addRequiredParam(const std::string &name, const std::string &doc_string)
MechanicsBaseNOSPD::validParams
static InputParameters validParams()
Definition: MechanicsBaseNOSPD.C:13
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _component

const unsigned int HorizonStabilizedFormIISmallStrainMechanicsNOSPD::_component
protected

The index of displacement component.

Definition at line 40 of file HorizonStabilizedFormIISmallStrainMechanicsNOSPD.h.

Referenced by computeLocalJacobian(), computeLocalOffDiagJacobian(), computeLocalResidual(), computeNonlocalJacobian(), computeNonlocalResidual(), and computePDNonlocalOffDiagJacobian().

◆ _current_unit_vec

RealGradient MechanicsBasePD::_current_unit_vec
protectedinherited

Unit vector of bond in current configuration.

Definition at line 75 of file MechanicsBasePD.h.

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

◆ _current_vec

RealGradient MechanicsBasePD::_current_vec
protectedinherited

Vector of bond in current configuration.

Definition at line 72 of file MechanicsBasePD.h.

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

◆ _ddgraddu

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_ddgraddu
protectedinherited

Definition at line 39 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsFiniteStrainBaseNOSPD::computeDinvFTDU(), MechanicsFiniteStrainBaseNOSPD::computeDJDU(), MechanicsFiniteStrainBaseNOSPD::computeDSDU(), and MechanicsBaseNOSPD::computeDSDU().

◆ _ddgraddv

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_ddgraddv
protectedinherited

Definition at line 40 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsFiniteStrainBaseNOSPD::computeDinvFTDU(), MechanicsFiniteStrainBaseNOSPD::computeDJDU(), MechanicsFiniteStrainBaseNOSPD::computeDSDU(), and MechanicsBaseNOSPD::computeDSDU().

◆ _ddgraddw

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_ddgraddw
protectedinherited

Definition at line 41 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsFiniteStrainBaseNOSPD::computeDinvFTDU(), MechanicsFiniteStrainBaseNOSPD::computeDJDU(), MechanicsFiniteStrainBaseNOSPD::computeDSDU(), and MechanicsBaseNOSPD::computeDSDU().

◆ _deigenstrain_dT

std::vector<const MaterialProperty<RankTwoTensor> *> MechanicsBaseNOSPD::_deigenstrain_dT
protectedinherited

Definition at line 44 of file MechanicsBaseNOSPD.h.

Referenced by WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeTempOffDiagJacobianScalar(), and MechanicsBaseNOSPD::MechanicsBaseNOSPD().

◆ _dgrad

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_dgrad
protectedinherited

Definition at line 38 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsFiniteStrainBaseNOSPD::computeDinvFTDU(), MechanicsFiniteStrainBaseNOSPD::computeDJDU(), MechanicsFiniteStrainBaseNOSPD::computeDSDFhat(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalResidual(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

◆ _disp_var

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

displacement variables

Definition at line 47 of file MechanicsBasePD.h.

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

◆ _eigenstrain_names

const std::vector<MaterialPropertyName> MechanicsBaseNOSPD::_eigenstrain_names
protectedinherited

Definition at line 43 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsBaseNOSPD::MechanicsBaseNOSPD().

◆ _ivardofs

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

Current variable dof numbers for nodes i and j.

Definition at line 66 of file MechanicsBasePD.h.

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

◆ _Jacobian_mult

const MaterialProperty<RankFourTensor>& MechanicsBaseNOSPD::_Jacobian_mult
protectedinherited

Definition at line 42 of file MechanicsBaseNOSPD.h.

Referenced by GeneralizedPlaneStrainOffDiagNOSPD::computeDispFullOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagNOSPD::computeDispPartialOffDiagJacobianScalar(), MechanicsFiniteStrainBaseNOSPD::computeDSDFhat(), MechanicsBaseNOSPD::computeDSDU(), WeakPlaneStressNOSPD::computeLocalJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeNonlocalJacobian(), computeNonlocalJacobian(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), computePDNonlocalOffDiagJacobian(), and GeneralizedPlaneStrainOffDiagNOSPD::computeTempOffDiagJacobianScalar().

◆ _multi

const MaterialProperty<Real>& MechanicsBaseNOSPD::_multi
protectedinherited

Material point based material properties.

Definition at line 35 of file MechanicsBaseNOSPD.h.

Referenced by GeneralizedPlaneStrainOffDiagNOSPD::computeDispFullOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagNOSPD::computeDispPartialOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalJacobian(), computeLocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalResidual(), computeLocalResidual(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalResidual(), ForceStabilizedSmallStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), computeNonlocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), computeNonlocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalResidual(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

◆ _ndisp

unsigned int MechanicsBasePD::_ndisp
protectedinherited

number of displacement components

Definition at line 55 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 63 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 58 of file MechanicsBasePD.h.

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

◆ _out_of_plane_strain_var

MooseVariable* MechanicsBasePD::_out_of_plane_strain_var
protectedinherited

Definition at line 59 of file MechanicsBasePD.h.

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

◆ _shape2

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_shape2
protectedinherited

Definition at line 37 of file MechanicsBaseNOSPD.h.

Referenced by GeneralizedPlaneStrainOffDiagNOSPD::computeDispFullOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagNOSPD::computeDispPartialOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalJacobian(), computeLocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalResidual(), computeLocalResidual(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalResidual(), ForceStabilizedSmallStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), computeNonlocalJacobian(), computeNonlocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalResidual(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and computePDNonlocalOffDiagJacobian().

◆ _stress

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_stress
protectedinherited

Definition at line 36 of file MechanicsBaseNOSPD.h.

Referenced by HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalResidual(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalResidual(), computeLocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalResidual(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalResidual(), computeNonlocalResidual(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

◆ _temp_coupled

const bool MechanicsBasePD::_temp_coupled
protectedinherited

Temperature variable.

Definition at line 50 of file MechanicsBasePD.h.

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

◆ _temp_var

MooseVariable* MechanicsBasePD::_temp_var
protectedinherited

Definition at line 51 of file MechanicsBasePD.h.

Referenced by MechanicsBPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagOSPD::computeOffDiagJacobianScalar(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), 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 69 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 Thu Jul 17 2025 01:48:19 for https://mooseframework.inl.gov by   doxygen 1.8.14