StressDivergenceTensors mostly copies from StressDivergence. More...

#include <StressDivergenceTensors.h>

Inheritance diagram for StressDivergenceTensors:

Public Member Functions
	StressDivergenceTensors (const InputParameters &parameters)

virtual void	computeJacobian () override

virtual void	computeOffDiagJacobian (MooseVariableFEBase &jvar) override

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	initialSetup () override

virtual void	computeResidual () override

virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

virtual void	computeFiniteDeformJacobian ()

virtual void	computeAverageGradientTest ()

virtual void	computeAverageGradientPhi ()

Protected Attributes
const std::string	_base_name

bool	_use_finite_deform_jacobian

const MaterialProperty< RankTwoTensor > &	_stress

const MaterialProperty< RankFourTensor > &	_Jacobian_mult

std::vector< RankFourTensor >	_finite_deform_Jacobian_mult

const MaterialProperty< RankTwoTensor > *	_deformation_gradient

const MaterialProperty< RankTwoTensor > *	_deformation_gradient_old

const MaterialProperty< RankTwoTensor > *	_rotation_increment

const unsigned int	_component

unsigned int	_ndisp
	Coupled displacement variables. More...

std::vector< unsigned int >	_disp_var

const bool	_temp_coupled

const unsigned int	_temp_var

const MaterialProperty< RankTwoTensor > *const	_deigenstrain_dT
	d(strain)/d(temperature), if computed by ComputeThermalExpansionEigenstrain More...

const bool	_out_of_plane_strain_coupled

const VariableValue *	_out_of_plane_strain

const unsigned int	_out_of_plane_strain_var

const unsigned int	_out_of_plane_direction

const bool	_use_displaced_mesh
	Whether this object is acting on the displaced mesh. More...

std::vector< std::vector< Real > >	_avg_grad_test
	Gradient of test function averaged over the element. Used in volumetric locking correction calculation. More...

std::vector< std::vector< Real > >	_avg_grad_phi
	Gradient of phi function averaged over the element. Used in volumetric locking correction calculation. More...

bool	_volumetric_locking_correction
	Flag for volumetric locking correction. More...

Assembly &	_assembly_undisplaced
	undisplaced problem More...

MooseVariable &	_var_undisplaced
	Reference to this Kernel's undisplaced MooseVariable object. More...


const VariablePhiGradient &	_grad_phi_undisplaced
	Shape and test functions on the undisplaced mesh. More...

const VariableTestGradient &	_grad_test_undisplaced

Detailed Description

StressDivergenceTensors mostly copies from StressDivergence.

There are small changes to use RankFourTensor and RankTwoTensors instead of SymmElasticityTensors and SymmTensors. This is done to allow for more mathematical transparancy.

Definition at line 27 of file StressDivergenceTensors.h.

Constructor & Destructor Documentation

◆ StressDivergenceTensors()

StressDivergenceTensors::StressDivergenceTensors ( const InputParameters & parameters )

Definition at line 62 of file StressDivergenceTensors.C.

   : ALEKernel(parameters),
     _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
     _use_finite_deform_jacobian(getParam<bool>("use_finite_deform_jacobian")),
     _stress(getMaterialPropertyByName<RankTwoTensor>(_base_name + "stress")),
     _Jacobian_mult(getMaterialPropertyByName<RankFourTensor>(_base_name + "Jacobian_mult")),
     _component(getParam<unsigned int>("component")),
     _ndisp(coupledComponents("displacements")),
     _disp_var(_ndisp),
     _temp_coupled(isCoupled("temperature")),
     _temp_var(_temp_coupled ? coupled("temperature") : 0),
     _deigenstrain_dT(_temp_coupled ? &getMaterialPropertyDerivative<RankTwoTensor>(
                                          getParam<std::string>("thermal_eigenstrain_name"),
                                          getVar("temperature", 0)->name())
                                    : nullptr),
     _out_of_plane_strain_coupled(isCoupled("out_of_plane_strain")),
     _out_of_plane_strain(_out_of_plane_strain_coupled ? &coupledValue("out_of_plane_strain")
                                                       : nullptr),
     _out_of_plane_strain_var(_out_of_plane_strain_coupled ? coupled("out_of_plane_strain") : 0),
     _out_of_plane_direction(getParam<MooseEnum>("out_of_plane_direction")),
     _use_displaced_mesh(getParam<bool>("use_displaced_mesh")),
     _avg_grad_test(_test.size(), std::vector<Real>(3, 0.0)),
     _avg_grad_phi(_phi.size(), std::vector<Real>(3, 0.0)),
     _volumetric_locking_correction(getParam<bool>("volumetric_locking_correction"))
 {
   for (unsigned int i = 0; i < _ndisp; ++i)
     _disp_var[i] = coupled("displacements", i);
  
   // Checking for consistency between mesh size and length of the provided displacements vector
   if (_out_of_plane_direction != 2 && _ndisp != 3)
     mooseError("For 2D simulations where the out-of-plane direction is x or y coordinate "
                "directions the number of supplied displacements must be three.");
   else if (_out_of_plane_direction == 2 && _ndisp != _mesh.dimension())
     mooseError("The number of displacement variables supplied must match the mesh dimension");
  
   if (_use_finite_deform_jacobian)
   {
     _deformation_gradient =
         &getMaterialProperty<RankTwoTensor>(_base_name + "deformation_gradient");
     _deformation_gradient_old =
         &getMaterialPropertyOld<RankTwoTensor>(_base_name + "deformation_gradient");
     _rotation_increment = &getMaterialProperty<RankTwoTensor>(_base_name + "rotation_increment");
   }
  
   // Error if volumetric locking correction is turned on for 1D problems
   if (_ndisp == 1 && _volumetric_locking_correction)
     mooseError("Volumetric locking correction should be set to false for 1-D problems.");
  
   // Generate warning when volumetric locking correction is used with second order elements
   if (_mesh.hasSecondOrderElements() && _volumetric_locking_correction)
     mooseWarning("Volumteric locking correction is not required for second order elements. Using "
                  "volumetric locking with second order elements could cause zigzag patterns in "
                  "stresses and strains.");
 }

Member Function Documentation

◆ computeAverageGradientPhi()

void StressDivergenceTensors::computeAverageGradientPhi ( )

protectedvirtual

Reimplemented in StressDivergenceRZTensors.

Definition at line 423 of file StressDivergenceTensors.C.

 {
   // Calculate volume average derivatives for phi
   _avg_grad_phi.resize(_phi.size());
   for (_i = 0; _i < _phi.size(); ++_i)
   {
     _avg_grad_phi[_i].resize(3);
     for (unsigned int component = 0; component < _mesh.dimension(); ++component)
     {
       _avg_grad_phi[_i][component] = 0.0;
       for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
         _avg_grad_phi[_i][component] += _grad_phi[_i][_qp](component) * _JxW[_qp] * _coord[_qp];
  
       _avg_grad_phi[_i][component] /= _current_elem_volume;
     }
   }
 }

Referenced by computeJacobian(), and computeOffDiagJacobian().

◆ computeAverageGradientTest()

void StressDivergenceTensors::computeAverageGradientTest ( )

protectedvirtual

Reimplemented in StressDivergenceRZTensors.

Definition at line 407 of file StressDivergenceTensors.C.

 {
   // Calculate volume averaged value of shape function derivative
   _avg_grad_test.resize(_test.size());
   for (_i = 0; _i < _test.size(); ++_i)
   {
     _avg_grad_test[_i].resize(3);
     _avg_grad_test[_i][_component] = 0.0;
     for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
       _avg_grad_test[_i][_component] += _grad_test[_i][_qp](_component) * _JxW[_qp] * _coord[_qp];
  
     _avg_grad_test[_i][_component] /= _current_elem_volume;
   }
 }

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

◆ computeFiniteDeformJacobian()

void StressDivergenceTensors::computeFiniteDeformJacobian ( )

protectedvirtual

Definition at line 350 of file StressDivergenceTensors.C.

 {
   const RankTwoTensor I(RankTwoTensor::initIdentity);
   const RankFourTensor II_ijkl = I.mixedProductIkJl(I);
  
   // Bring back to unrotated config
   const RankTwoTensor unrotated_stress =
       (*_rotation_increment)[_qp].transpose() * _stress[_qp] * (*_rotation_increment)[_qp];
  
   // Incremental deformation gradient Fhat
   const RankTwoTensor Fhat =
       (*_deformation_gradient)[_qp] * (*_deformation_gradient_old)[_qp].inverse();
   const RankTwoTensor Fhatinv = Fhat.inverse();
  
   const RankTwoTensor rot_times_stress = (*_rotation_increment)[_qp] * unrotated_stress;
   const RankFourTensor dstress_drot =
       I.mixedProductIkJl(rot_times_stress) + I.mixedProductJkIl(rot_times_stress);
   const RankFourTensor rot_rank_four =
       (*_rotation_increment)[_qp].mixedProductIkJl((*_rotation_increment)[_qp]);
   const RankFourTensor drot_dUhatinv = Fhat.mixedProductIkJl(I);
  
   const RankTwoTensor A = I - Fhatinv;
  
   // Ctilde = Chat^-1 - I
   const RankTwoTensor Ctilde = A * A.transpose() - A - A.transpose();
   const RankFourTensor dCtilde_dFhatinv =
       -I.mixedProductIkJl(A) - I.mixedProductJkIl(A) + II_ijkl + I.mixedProductJkIl(I);
  
   // Second order approximation of Uhat - consistent with strain increment definition
   // const RankTwoTensor Uhat = I - 0.5 * Ctilde - 3.0/8.0 * Ctilde * Ctilde;
  
   RankFourTensor dUhatinv_dCtilde =
       0.5 * II_ijkl - 1.0 / 8.0 * (I.mixedProductIkJl(Ctilde) + Ctilde.mixedProductIkJl(I));
   RankFourTensor drot_dFhatinv = drot_dUhatinv * dUhatinv_dCtilde * dCtilde_dFhatinv;
  
   drot_dFhatinv -= Fhat.mixedProductIkJl((*_rotation_increment)[_qp].transpose());
   _finite_deform_Jacobian_mult[_qp] = dstress_drot * drot_dFhatinv;
  
   const RankFourTensor dstrain_increment_dCtilde =
       -0.5 * II_ijkl + 0.25 * (I.mixedProductIkJl(Ctilde) + Ctilde.mixedProductIkJl(I));
   _finite_deform_Jacobian_mult[_qp] +=
       rot_rank_four * _Jacobian_mult[_qp] * dstrain_increment_dCtilde * dCtilde_dFhatinv;
   _finite_deform_Jacobian_mult[_qp] += Fhat.mixedProductJkIl(_stress[_qp]);
  
   const RankFourTensor dFhat_dFhatinv = -Fhat.mixedProductIkJl(Fhat.transpose());
   const RankTwoTensor dJ_dFhatinv = dFhat_dFhatinv.innerProductTranspose(Fhat.ddet());
  
   // Component from Jacobian derivative
   _finite_deform_Jacobian_mult[_qp] += _stress[_qp].outerProduct(dJ_dFhatinv);
  
   // Derivative of Fhatinv w.r.t. undisplaced coordinates
   const RankTwoTensor Finv = (*_deformation_gradient)[_qp].inverse();
   const RankFourTensor dFhatinv_dGradu = -Fhatinv.mixedProductIkJl(Finv.transpose());
   _finite_deform_Jacobian_mult[_qp] = _finite_deform_Jacobian_mult[_qp] * dFhatinv_dGradu;
 }

Referenced by computeJacobian(), and computeOffDiagJacobian().

◆ computeJacobian()

void StressDivergenceTensors::computeJacobian ( )

overridevirtual

Reimplemented from ALEKernel.

Definition at line 168 of file StressDivergenceTensors.C.

 {
   if (_volumetric_locking_correction)
   {
     computeAverageGradientTest();
     computeAverageGradientPhi();
   }
  
   if (_use_finite_deform_jacobian)
   {
     _finite_deform_Jacobian_mult.resize(_qrule->n_points());
  
     for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
       computeFiniteDeformJacobian();
  
     ALEKernel::computeJacobian();
   }
   else
     Kernel::computeJacobian();
 }

◆ computeOffDiagJacobian()

void StressDivergenceTensors::computeOffDiagJacobian ( MooseVariableFEBase & jvar )

overridevirtual

Reimplemented from ALEKernel.

Definition at line 190 of file StressDivergenceTensors.C.

 {
   if (_volumetric_locking_correction)
   {
     computeAverageGradientPhi();
     computeAverageGradientTest();
   }
  
   if (_use_finite_deform_jacobian)
   {
     _finite_deform_Jacobian_mult.resize(_qrule->n_points());
  
     for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
       computeFiniteDeformJacobian();
  
     ALEKernel::computeOffDiagJacobian(jvar);
   }
   else
     Kernel::computeOffDiagJacobian(jvar);
 }

◆ computeQpJacobian()

Real StressDivergenceTensors::computeQpJacobian ( )

overrideprotectedvirtual

Reimplemented in StressDivergenceRSphericalTensors, StressDivergenceRZTensors, and DynamicStressDivergenceTensors.

Definition at line 212 of file StressDivergenceTensors.C.

 {
   if (_use_finite_deform_jacobian)
     return ElasticityTensorTools::elasticJacobian(_finite_deform_Jacobian_mult[_qp],
                                                   _component,
                                                   _component,
                                                   _grad_test[_i][_qp],
                                                   _grad_phi_undisplaced[_j][_qp]);
  
   Real sum_C3x3 = _Jacobian_mult[_qp].sum3x3();
   RealGradient sum_C3x1 = _Jacobian_mult[_qp].sum3x1();
  
   Real jacobian = 0.0;
   // B^T_i * C * B_j
   jacobian += ElasticityTensorTools::elasticJacobian(
       _Jacobian_mult[_qp], _component, _component, _grad_test[_i][_qp], _grad_phi[_j][_qp]);
  
   if (_volumetric_locking_correction)
   {
     // jacobian = Bbar^T_i * C * Bbar_j where Bbar = B + Bvol
     // jacobian = B^T_i * C * B_j + Bvol^T_i * C * Bvol_j +  Bvol^T_i * C * B_j + B^T_i * C * Bvol_j
  
     // Bvol^T_i * C * Bvol_j
     jacobian += sum_C3x3 * (_avg_grad_test[_i][_component] - _grad_test[_i][_qp](_component)) *
                 (_avg_grad_phi[_j][_component] - _grad_phi[_j][_qp](_component)) / 9.0;
  
     // B^T_i * C * Bvol_j
     jacobian += sum_C3x1(_component) * _grad_test[_i][_qp](_component) *
                 (_avg_grad_phi[_j][_component] - _grad_phi[_j][_qp](_component)) / 3.0;
  
     // Bvol^T_i * C * B_j
     RankTwoTensor phi;
     if (_component == 0)
     {
       phi(0, 0) = _grad_phi[_j][_qp](0);
       phi(0, 1) = phi(1, 0) = _grad_phi[_j][_qp](1);
       phi(0, 2) = phi(2, 0) = _grad_phi[_j][_qp](2);
     }
     else if (_component == 1)
     {
       phi(1, 1) = _grad_phi[_j][_qp](1);
       phi(0, 1) = phi(1, 0) = _grad_phi[_j][_qp](0);
       phi(1, 2) = phi(2, 1) = _grad_phi[_j][_qp](2);
     }
     else if (_component == 2)
     {
       phi(2, 2) = _grad_phi[_j][_qp](2);
       phi(0, 2) = phi(2, 0) = _grad_phi[_j][_qp](0);
       phi(1, 2) = phi(2, 1) = _grad_phi[_j][_qp](1);
     }
  
     jacobian += (_Jacobian_mult[_qp] * phi).trace() *
                 (_avg_grad_test[_i][_component] - _grad_test[_i][_qp](_component)) / 3.0;
   }
  
   if (_ndisp != 3 && _out_of_plane_strain_coupled && _use_displaced_mesh)
   {
     const Real out_of_plane_thickness = std::exp((*_out_of_plane_strain)[_qp]);
     jacobian *= out_of_plane_thickness;
   }
  
   return jacobian;
 }

Referenced by DynamicStressDivergenceTensors::computeQpJacobian().

◆ computeQpOffDiagJacobian()

Real StressDivergenceTensors::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Reimplemented in StressDivergenceRSphericalTensors, StressDivergenceRZTensors, CosseratStressDivergenceTensors, and DynamicStressDivergenceTensors.

Definition at line 277 of file StressDivergenceTensors.C.

 {
   // off-diagonal Jacobian with respect to a coupled displacement component
   for (unsigned int coupled_component = 0; coupled_component < _ndisp; ++coupled_component)
     if (jvar == _disp_var[coupled_component])
     {
       if (_out_of_plane_direction != 2)
       {
         if (coupled_component == _out_of_plane_direction)
           continue;
       }
  
       if (_use_finite_deform_jacobian)
         return ElasticityTensorTools::elasticJacobian(_finite_deform_Jacobian_mult[_qp],
                                                       _component,
                                                       coupled_component,
                                                       _grad_test[_i][_qp],
                                                       _grad_phi_undisplaced[_j][_qp]);
  
       const Real sum_C3x3 = _Jacobian_mult[_qp].sum3x3();
       const RealGradient sum_C3x1 = _Jacobian_mult[_qp].sum3x1();
       Real jacobian = 0.0;
  
       // B^T_i * C * B_j
       jacobian += ElasticityTensorTools::elasticJacobian(_Jacobian_mult[_qp],
                                                          _component,
                                                          coupled_component,
                                                          _grad_test[_i][_qp],
                                                          _grad_phi[_j][_qp]);
  
       if (_volumetric_locking_correction)
       {
         // jacobian = Bbar^T_i * C * Bbar_j where Bbar = B + Bvol
         // jacobian = B^T_i * C * B_j + Bvol^T_i * C * Bvol_j +  Bvol^T_i * C * B_j + B^T_i * C *
         // Bvol_j
  
         // Bvol^T_i * C * Bvol_j
         jacobian += sum_C3x3 * (_avg_grad_test[_i][_component] - _grad_test[_i][_qp](_component)) *
                     (_avg_grad_phi[_j][coupled_component] - _grad_phi[_j][_qp](coupled_component)) /
                     9.0;
  
         // B^T_i * C * Bvol_j
         jacobian += sum_C3x1(_component) * _grad_test[_i][_qp](_component) *
                     (_avg_grad_phi[_j][coupled_component] - _grad_phi[_j][_qp](coupled_component)) /
                     3.0;
  
         // Bvol^T_i * C * B_i
         RankTwoTensor phi;
         for (unsigned int i = 0; i < 3; ++i)
           phi(coupled_component, i) = _grad_phi[_j][_qp](i);
  
         jacobian += (_Jacobian_mult[_qp] * phi).trace() *
                     (_avg_grad_test[_i][_component] - _grad_test[_i][_qp](_component)) / 3.0;
       }
  
       return jacobian;
     }
  
   // off-diagonal Jacobian with respect to a coupled out_of_plane_strain variable
   if (_out_of_plane_strain_coupled && jvar == _out_of_plane_strain_var)
     return _Jacobian_mult[_qp](
                _component, _component, _out_of_plane_direction, _out_of_plane_direction) *
            _grad_test[_i][_qp](_component) * _phi[_j][_qp];
  
   // off-diagonal Jacobian with respect to a coupled temperature variable
   if (_temp_coupled && jvar == _temp_var)
     return -((_Jacobian_mult[_qp] * (*_deigenstrain_dT)[_qp]) *
              _grad_test[_i][_qp])(_component)*_phi[_j][_qp];
  
   return 0.0;
 }

Referenced by CosseratStressDivergenceTensors::computeQpOffDiagJacobian(), and DynamicStressDivergenceTensors::computeQpOffDiagJacobian().

◆ computeQpResidual()

Real StressDivergenceTensors::computeQpResidual ( )

overrideprotectedvirtual

Reimplemented in StressDivergenceRSphericalTensors, StressDivergenceRZTensors, and DynamicStressDivergenceTensors.

Definition at line 149 of file StressDivergenceTensors.C.

 {
  
   Real residual = _stress[_qp].row(_component) * _grad_test[_i][_qp];
   // volumetric locking correction
   if (_volumetric_locking_correction)
     residual += _stress[_qp].trace() / 3.0 *
                 (_avg_grad_test[_i][_component] - _grad_test[_i][_qp](_component));
  
   if (_ndisp != 3 && _out_of_plane_strain_coupled && _use_displaced_mesh)
   {
     const Real out_of_plane_thickness = std::exp((*_out_of_plane_strain)[_qp]);
     residual *= out_of_plane_thickness;
   }
  
   return residual;
 }

Referenced by computeResidual().

◆ computeResidual()

void StressDivergenceTensors::computeResidual ( )

overrideprotectedvirtual

Definition at line 126 of file StressDivergenceTensors.C.

 {
   prepareVectorTag(_assembly, _var.number());
  
   if (_volumetric_locking_correction)
     computeAverageGradientTest();
  
   precalculateResidual();
   for (_i = 0; _i < _test.size(); ++_i)
     for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
       _local_re(_i) += _JxW[_qp] * _coord[_qp] * computeQpResidual();
  
   accumulateTaggedLocalResidual();
  
   if (_has_save_in)
   {
     Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
     for (const auto & var : _save_in)
       var->sys().solution().add_vector(_local_re, var->dofIndices());
   }
 }

◆ initialSetup()

void StressDivergenceTensors::initialSetup ( )

overrideprotectedvirtual

Reimplemented in StressDivergenceRSphericalTensors, and StressDivergenceRZTensors.

Definition at line 118 of file StressDivergenceTensors.C.

 {
   if (getBlockCoordSystem() != Moose::COORD_XYZ)
     mooseError(
         "The coordinate system in the Problem block must be set to XYZ for cartesian geometries.");
 }

◆ validParams()

InputParameters StressDivergenceTensors::validParams ( )

static

Definition at line 26 of file StressDivergenceTensors.C.

 {
   InputParameters params = ALEKernel::validParams();
   params.addClassDescription("Stress divergence kernel for the Cartesian coordinate system");
   params.addRequiredParam<unsigned int>("component",
                                         "An integer corresponding to the direction "
                                         "the variable this kernel acts in. (0 for x, "
                                         "1 for y, 2 for z)");
   params.addRequiredCoupledVar("displacements",
                                "The string of displacements suitable for the problem statement");
   params.addCoupledVar("temperature",
                        "The name of the temperature variable used in the "
                        "ComputeThermalExpansionEigenstrain.  (Not required for "
                        "simulations without temperature coupling.)");
   params.addParam<std::string>(
       "thermal_eigenstrain_name",
       "thermal_eigenstrain",
       "The eigenstrain_name used in the ComputeThermalExpansionEigenstrain.");
   params.addCoupledVar("out_of_plane_strain",
                        "The name of the out_of_plane_strain variable used in the "
                        "WeakPlaneStress kernel.");
   MooseEnum out_of_plane_direction("x y z", "z");
   params.addParam<MooseEnum>(
       "out_of_plane_direction",
       out_of_plane_direction,
       "The direction of the out_of_plane_strain variable used in the WeakPlaneStress kernel.");
   params.addParam<std::string>("base_name", "Material property base name");
   params.set<bool>("use_displaced_mesh") = false;
   params.addParam<bool>(
       "use_finite_deform_jacobian", false, "Jacobian for corotational finite strain");
   params.addParam<bool>("volumetric_locking_correction",
                         false,
                         "Set to false to turn off volumetric locking correction");
   return params;
 }

Referenced by DynamicStressDivergenceTensors::validParams(), CosseratStressDivergenceTensors::validParams(), StressDivergenceRZTensors::validParams(), and StressDivergenceRSphericalTensors::validParams().

Member Data Documentation

◆ _assembly_undisplaced

Assembly& ALEKernel::_assembly_undisplaced

protectedinherited

undisplaced problem

Definition at line 34 of file ALEKernel.h.

◆ _avg_grad_phi

std::vector<std::vector<Real> > StressDivergenceTensors::_avg_grad_phi

protected

Gradient of phi function averaged over the element. Used in volumetric locking correction calculation.

Definition at line 86 of file StressDivergenceTensors.h.

Referenced by StressDivergenceRZTensors::calculateJacobian(), StressDivergenceRZTensors::computeAverageGradientPhi(), computeAverageGradientPhi(), computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _avg_grad_test

std::vector<std::vector<Real> > StressDivergenceTensors::_avg_grad_test

protected

Gradient of test function averaged over the element. Used in volumetric locking correction calculation.

Definition at line 83 of file StressDivergenceTensors.h.

Referenced by StressDivergenceRZTensors::calculateJacobian(), StressDivergenceRZTensors::computeAverageGradientTest(), computeAverageGradientTest(), computeQpJacobian(), computeQpOffDiagJacobian(), DynamicStressDivergenceTensors::computeQpResidual(), StressDivergenceRZTensors::computeQpResidual(), and computeQpResidual().

◆ _base_name

const std::string StressDivergenceTensors::_base_name

protected

Definition at line 50 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors().

◆ _component

const unsigned int StressDivergenceTensors::_component

protected

Definition at line 62 of file StressDivergenceTensors.h.

Referenced by StressDivergenceRZTensors::computeAverageGradientTest(), computeAverageGradientTest(), StressDivergenceRZTensors::computeQpJacobian(), computeQpJacobian(), StressDivergenceRSphericalTensors::computeQpJacobian(), CosseratStressDivergenceTensors::computeQpOffDiagJacobian(), StressDivergenceRZTensors::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), StressDivergenceRSphericalTensors::computeQpOffDiagJacobian(), DynamicStressDivergenceTensors::computeQpResidual(), StressDivergenceRZTensors::computeQpResidual(), computeQpResidual(), and StressDivergenceRSphericalTensors::StressDivergenceRSphericalTensors().

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_deformation_gradient

protected

Definition at line 57 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors().

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_deformation_gradient_old

protected

Definition at line 58 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors().

◆ _deigenstrain_dT

const MaterialProperty<RankTwoTensor>* const StressDivergenceTensors::_deigenstrain_dT

protected

d(strain)/d(temperature), if computed by ComputeThermalExpansionEigenstrain

Definition at line 72 of file StressDivergenceTensors.h.

◆ _disp_var

std::vector<unsigned int> StressDivergenceTensors::_disp_var

protected

Definition at line 66 of file StressDivergenceTensors.h.

Referenced by DynamicStressDivergenceTensors::computeQpOffDiagJacobian(), StressDivergenceRZTensors::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), StressDivergenceRSphericalTensors::computeQpOffDiagJacobian(), and StressDivergenceTensors().

◆ _finite_deform_Jacobian_mult

std::vector<RankFourTensor> StressDivergenceTensors::_finite_deform_Jacobian_mult

protected

Definition at line 56 of file StressDivergenceTensors.h.

Referenced by computeFiniteDeformJacobian(), computeJacobian(), computeOffDiagJacobian(), computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _grad_phi_undisplaced

const VariablePhiGradient& ALEKernel::_grad_phi_undisplaced

protectedinherited

Shape and test functions on the undisplaced mesh.

Definition at line 40 of file ALEKernel.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _grad_test_undisplaced

const VariableTestGradient& ALEKernel::_grad_test_undisplaced

protectedinherited

Definition at line 41 of file ALEKernel.h.

◆ _Jacobian_mult

const MaterialProperty<RankFourTensor>& StressDivergenceTensors::_Jacobian_mult

protected

Definition at line 54 of file StressDivergenceTensors.h.

Referenced by StressDivergenceRSphericalTensors::calculateJacobian(), StressDivergenceRZTensors::calculateJacobian(), computeFiniteDeformJacobian(), computeQpJacobian(), CosseratStressDivergenceTensors::computeQpOffDiagJacobian(), StressDivergenceRZTensors::computeQpOffDiagJacobian(), and computeQpOffDiagJacobian().

◆ _ndisp

unsigned int StressDivergenceTensors::_ndisp

protected

Coupled displacement variables.

Definition at line 65 of file StressDivergenceTensors.h.

Referenced by computeQpJacobian(), DynamicStressDivergenceTensors::computeQpOffDiagJacobian(), StressDivergenceRZTensors::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), StressDivergenceRSphericalTensors::computeQpOffDiagJacobian(), computeQpResidual(), and StressDivergenceTensors().

◆ _out_of_plane_direction

const unsigned int StressDivergenceTensors::_out_of_plane_direction

protected

Definition at line 77 of file StressDivergenceTensors.h.

Referenced by computeQpOffDiagJacobian(), and StressDivergenceTensors().

◆ _out_of_plane_strain

const VariableValue* StressDivergenceTensors::_out_of_plane_strain

protected

Definition at line 75 of file StressDivergenceTensors.h.

Referenced by computeQpJacobian(), and computeQpResidual().

◆ _out_of_plane_strain_coupled

const bool StressDivergenceTensors::_out_of_plane_strain_coupled

protected

Definition at line 74 of file StressDivergenceTensors.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

◆ _out_of_plane_strain_var

const unsigned int StressDivergenceTensors::_out_of_plane_strain_var

protected

Definition at line 76 of file StressDivergenceTensors.h.

Referenced by computeQpOffDiagJacobian().

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_rotation_increment

protected

Definition at line 59 of file StressDivergenceTensors.h.

Referenced by computeFiniteDeformJacobian(), and StressDivergenceTensors().

◆ _stress

const MaterialProperty<RankTwoTensor>& StressDivergenceTensors::_stress

protected

Definition at line 53 of file StressDivergenceTensors.h.

Referenced by computeFiniteDeformJacobian(), DynamicStressDivergenceTensors::computeQpResidual(), StressDivergenceRZTensors::computeQpResidual(), computeQpResidual(), and StressDivergenceRSphericalTensors::computeQpResidual().

◆ _temp_coupled

const bool StressDivergenceTensors::_temp_coupled

protected

Definition at line 68 of file StressDivergenceTensors.h.

Referenced by DynamicStressDivergenceTensors::computeQpOffDiagJacobian(), StressDivergenceRZTensors::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), and StressDivergenceRSphericalTensors::computeQpOffDiagJacobian().

◆ _temp_var

const unsigned int StressDivergenceTensors::_temp_var

protected

Definition at line 69 of file StressDivergenceTensors.h.

Referenced by DynamicStressDivergenceTensors::computeQpOffDiagJacobian(), StressDivergenceRZTensors::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), and StressDivergenceRSphericalTensors::computeQpOffDiagJacobian().

◆ _use_displaced_mesh

const bool StressDivergenceTensors::_use_displaced_mesh

protected

Whether this object is acting on the displaced mesh.

Definition at line 80 of file StressDivergenceTensors.h.

Referenced by computeQpJacobian(), and computeQpResidual().

◆ _use_finite_deform_jacobian

bool StressDivergenceTensors::_use_finite_deform_jacobian

protected

Definition at line 51 of file StressDivergenceTensors.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), computeQpJacobian(), computeQpOffDiagJacobian(), and StressDivergenceTensors().

◆ _var_undisplaced

MooseVariable& ALEKernel::_var_undisplaced

protectedinherited

Reference to this Kernel's undisplaced MooseVariable object.

Definition at line 37 of file ALEKernel.h.

◆ _volumetric_locking_correction

bool StressDivergenceTensors::_volumetric_locking_correction

protected

Flag for volumetric locking correction.

Definition at line 89 of file StressDivergenceTensors.h.

Referenced by StressDivergenceRZTensors::calculateJacobian(), computeJacobian(), computeOffDiagJacobian(), computeQpJacobian(), computeQpOffDiagJacobian(), DynamicStressDivergenceTensors::computeQpResidual(), StressDivergenceRZTensors::computeQpResidual(), computeQpResidual(), computeResidual(), and StressDivergenceTensors().

The documentation for this class was generated from the following files:

tensor_mechanics/include/kernels/StressDivergenceTensors.h
tensor_mechanics/src/kernels/StressDivergenceTensors.C

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ StressDivergenceTensors()

Member Function Documentation

◆ computeAverageGradientPhi()

◆ computeAverageGradientTest()

◆ computeFiniteDeformJacobian()

◆ computeJacobian()

◆ computeOffDiagJacobian()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ computeResidual()

◆ initialSetup()

◆ validParams()

Member Data Documentation

◆ _assembly_undisplaced

◆ _avg_grad_phi

◆ _avg_grad_test

◆ _base_name

◆ _component

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _deigenstrain_dT

◆ _disp_var

◆ _finite_deform_Jacobian_mult

◆ _grad_phi_undisplaced

◆ _grad_test_undisplaced

◆ _Jacobian_mult

◆ _ndisp

◆ _out_of_plane_direction

◆ _out_of_plane_strain

◆ _out_of_plane_strain_coupled

◆ _out_of_plane_strain_var

◆ _rotation_increment

◆ _stress

◆ _temp_coupled

◆ _temp_var

◆ _use_displaced_mesh

◆ _use_finite_deform_jacobian

◆ _var_undisplaced

◆ _volumetric_locking_correction