StressDivergenceTensors mostly copies from StressDivergence. More...

#include <StressDivergenceTensors.h>

Inheritance diagram for StressDivergenceTensors:

Public Types
typedef std::vector< int >	JvarMap

Public Member Functions
	StressDivergenceTensors (const InputParameters &parameters)

virtual void	computeJacobian () override

virtual void	computeOffDiagJacobian (unsigned int jvar) override

unsigned int	mapJvarToCvar (unsigned int jvar)

int	mapJvarToCvar (unsigned int jvar, const JvarMap &jvar_map)

bool	mapJvarToCvar (unsigned int jvar, unsigned int &cvar)

const JvarMap &	getJvarMap ()

const JvarMap &	getParameterJvarMap (std::string parameter_name)

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
	Base name of the material system that this kernel applies to. More...

bool	_use_finite_deform_jacobian

const MaterialProperty< RankTwoTensor > &	_stress
	The stress tensor that the divergence operator operates on. More...

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
	An integer corresponding to the direction this kernel acts in. More...

unsigned int	_ndisp
	Coupled displacement variables. More...

std::vector< unsigned int >	_disp_var
	Displacement variables IDs. More...

std::vector< std::vector< const MaterialProperty< RankTwoTensor > * > >	_deigenstrain_dargs
	eigen strain derivatives wrt coupled variables More...

const bool	_out_of_plane_strain_coupled

const VariableValue *const	_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...

const unsigned int	_n_args

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 24 of file StressDivergenceTensors.h.

Constructor & Destructor Documentation

◆ StressDivergenceTensors()

StressDivergenceTensors::StressDivergenceTensors ( const InputParameters & parameters )

Definition at line 65 of file StressDivergenceTensors.C.

   : JvarMapKernelInterface<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),
     _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"))
 {
   // get coupled displacements
   for (unsigned int i = 0; i < _ndisp; ++i)
     _disp_var[i] = coupled("displacements", i);
 
   // fetch eigenstrain derivatives
   const auto nvar = _coupled_moose_vars.size();
   _deigenstrain_dargs.resize(nvar);
   for (std::size_t i = 0; i < nvar; ++i)
     for (auto eigenstrain_name : getParam<std::vector<MaterialPropertyName>>("eigenstrain_names"))
       _deigenstrain_dargs[i].push_back(&getMaterialPropertyDerivative<RankTwoTensor>(
           eigenstrain_name, _coupled_moose_vars[i]->name()));
 
   // 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 443 of file StressDivergenceTensors.C.

Referenced by computeJacobian(), and computeOffDiagJacobian().

 {
   // 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;
     }
   }
 }

◆ computeAverageGradientTest()

void StressDivergenceTensors::computeAverageGradientTest ( )

protectedvirtual

Reimplemented in StressDivergenceRZTensors.

Definition at line 427 of file StressDivergenceTensors.C.

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

 {
   // 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;
   }
 }

◆ computeFiniteDeformJacobian()

void StressDivergenceTensors::computeFiniteDeformJacobian ( )

protectedvirtual

Definition at line 369 of file StressDivergenceTensors.C.

Referenced by computeJacobian(), and computeOffDiagJacobian().

 {
   usingTensorIndices(i_, j_, k_, l_);
   const auto I = RankTwoTensor::Identity();
   const RankFourTensor I2 = I.times<i_, k_, j_, l_>(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.times<i_, k_, j_, l_>(rot_times_stress) + I.times<j_, k_, i_, l_>(rot_times_stress);
   const RankFourTensor rot_rank_four =
       (*_rotation_increment)[_qp].times<i_, k_, j_, l_>((*_rotation_increment)[_qp]);
   const RankFourTensor drot_dUhatinv = Fhat.times<i_, k_, j_, l_>(I);
 
   const RankTwoTensor A = I - Fhatinv;
 
   // Ctilde = Chat^-1 - I
   const RankTwoTensor Ctilde = A * A.transpose() - A - A.transpose();
   const RankFourTensor dCtilde_dFhatinv =
       -I.times<i_, k_, j_, l_>(A) - I.times<j_, k_, i_, l_>(A) + I2 + I.times<j_, k_, i_, l_>(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 * I2 - 1.0 / 8.0 * (I.times<i_, k_, j_, l_>(Ctilde) + Ctilde.times<i_, k_, j_, l_>(I));
   RankFourTensor drot_dFhatinv = drot_dUhatinv * dUhatinv_dCtilde * dCtilde_dFhatinv;
 
   drot_dFhatinv -= Fhat.times<i_, k_, j_, l_>((*_rotation_increment)[_qp].transpose());
   _finite_deform_Jacobian_mult[_qp] = dstress_drot * drot_dFhatinv;
 
   const RankFourTensor dstrain_increment_dCtilde =
       -0.5 * I2 + 0.25 * (I.times<i_, k_, j_, l_>(Ctilde) + Ctilde.times<i_, k_, j_, l_>(I));
   _finite_deform_Jacobian_mult[_qp] +=
       rot_rank_four * _Jacobian_mult[_qp] * dstrain_increment_dCtilde * dCtilde_dFhatinv;
   _finite_deform_Jacobian_mult[_qp] += Fhat.times<j_, k_, i_, l_>(_stress[_qp]);
 
   const RankFourTensor dFhat_dFhatinv = -Fhat.times<i_, k_, j_, l_>(Fhat.transpose());
   const RankTwoTensor dJ_dFhatinv = dFhat_dFhatinv.innerProductTranspose(Fhat.ddet());
 
   // Component from Jacobian derivative
   _finite_deform_Jacobian_mult[_qp] += _stress[_qp].times<i_, j_, k_, l_>(dJ_dFhatinv);
 
   // Derivative of Fhatinv w.r.t. undisplaced coordinates
   const RankTwoTensor Finv = (*_deformation_gradient)[_qp].inverse();
   const RankFourTensor dFhatinv_dGradu = -Fhatinv.times<i_, k_, j_, l_>(Finv.transpose());
   _finite_deform_Jacobian_mult[_qp] = _finite_deform_Jacobian_mult[_qp] * dFhatinv_dGradu;
 }

◆ computeJacobian()

void StressDivergenceTensors::computeJacobian ( )

overridevirtual

Definition at line 178 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 ( unsigned int jvar )

overridevirtual

Reimplemented from JvarMapKernelInterface< ALEKernel >.

Definition at line 200 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 222 of file StressDivergenceTensors.C.

Referenced by DynamicStressDivergenceTensors::computeQpJacobian().

 {
   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;
     switch (_component)
     {
       case 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);
         break;
 
       case 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);
         break;
 
       case 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);
         break;
     }
 
     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;
 }

◆ computeQpOffDiagJacobian()

Real StressDivergenceTensors::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Reimplemented in StressDivergenceRSphericalTensors, StressDivergenceRZTensors, CosseratStressDivergenceTensors, and DynamicStressDivergenceTensors.

Definition at line 290 of file StressDivergenceTensors.C.

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

 {
   // 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];
 
   // bail out if jvar is not coupled
   if (getJvarMap()[jvar] < 0)
     return 0.0;
 
   // off-diagonal Jacobian with respect to any other coupled variable
   const unsigned int cvar = mapJvarToCvar(jvar);
   RankTwoTensor total_deigenstrain;
   for (const auto deigenstrain_darg : _deigenstrain_dargs[cvar])
     total_deigenstrain += (*deigenstrain_darg)[_qp];
 
   return -((_Jacobian_mult[_qp] * total_deigenstrain) *
            _grad_test[_i][_qp])(_component)*_phi[_j][_qp];
 }

◆ computeQpResidual()

Real StressDivergenceTensors::computeQpResidual ( )

overrideprotectedvirtual

Reimplemented in StressDivergenceRSphericalTensors, StressDivergenceRZTensors, and DynamicStressDivergenceTensors.

Definition at line 160 of file StressDivergenceTensors.C.

Referenced by computeResidual().

 {
   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;
 }

◆ computeResidual()

void StressDivergenceTensors::computeResidual ( )

overrideprotectedvirtual

Definition at line 137 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 124 of file StressDivergenceTensors.C.

 {
   // check if any of the eigenstrains provide derivatives wrt variables that are not coupled
   for (auto eigenstrain_name : getParam<std::vector<MaterialPropertyName>>("eigenstrain_names"))
     validateNonlinearCoupling<RankTwoTensor>(eigenstrain_name);
 
   // make sure the coordinate system is cartesioan
   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 24 of file StressDivergenceTensors.C.

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

 {
   InputParameters params = JvarMapKernelInterface<ALEKernel>::validParams();
   params.addClassDescription("Stress divergence kernel for the Cartesian coordinate system");
   params.addRequiredRangeCheckedParam<unsigned int>("component",
                                                     "component < 3",
                                                     "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");
 
   // maybe this should be deprecated in favor of args
   params.addCoupledVar("temperature",
                        "The name of the temperature variable used in the "
                        "ComputeThermalExpansionEigenstrain.  (Not required for "
                        "simulations without temperature coupling.)");
 
   params.addParam<std::vector<MaterialPropertyName>>(
       "eigenstrain_names",
       {},
       "List of eigenstrains used in the strain calculation. Used for computing their derivatives "
       "for off-diagonal Jacobian terms.");
   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;
 }

Member Data Documentation

◆ _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 85 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 82 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

Base name of the material system that this kernel applies to.

Definition at line 47 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors().

◆ _component

const unsigned int StressDivergenceTensors::_component

protected

An integer corresponding to the direction this kernel acts in.

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 56 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors().

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_deformation_gradient_old

protected

Definition at line 57 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors().

◆ _deigenstrain_dargs

std::vector<std::vector<const MaterialProperty<RankTwoTensor> *> > StressDivergenceTensors::_deigenstrain_dargs

protected

eigen strain derivatives wrt coupled variables

Definition at line 71 of file StressDivergenceTensors.h.

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

◆ _disp_var

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

protected

Displacement variables IDs.

Definition at line 68 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 55 of file StressDivergenceTensors.h.

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

◆ _Jacobian_mult

const MaterialProperty<RankFourTensor>& StressDivergenceTensors::_Jacobian_mult

protected

Definition at line 53 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 76 of file StressDivergenceTensors.h.

Referenced by computeQpOffDiagJacobian(), and StressDivergenceTensors().

◆ _out_of_plane_strain

const VariableValue* const StressDivergenceTensors::_out_of_plane_strain

protected

Definition at line 74 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 73 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 75 of file StressDivergenceTensors.h.

Referenced by computeQpOffDiagJacobian().

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_rotation_increment

protected

Definition at line 58 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors().

◆ _stress

const MaterialProperty<RankTwoTensor>& StressDivergenceTensors::_stress

protected

The stress tensor that the divergence operator operates on.

Definition at line 52 of file StressDivergenceTensors.h.

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

◆ _use_displaced_mesh

const bool StressDivergenceTensors::_use_displaced_mesh

protected

Whether this object is acting on the displaced mesh.

Definition at line 79 of file StressDivergenceTensors.h.

Referenced by computeQpJacobian(), and computeQpResidual().

◆ _use_finite_deform_jacobian

bool StressDivergenceTensors::_use_finite_deform_jacobian

protected

Definition at line 49 of file StressDivergenceTensors.h.

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

◆ _volumetric_locking_correction

bool StressDivergenceTensors::_volumetric_locking_correction

protected

Flag for volumetric locking correction.

Definition at line 88 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:

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

Public Types

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

◆ _avg_grad_phi

◆ _avg_grad_test

◆ _base_name

◆ _component

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _deigenstrain_dargs

◆ _disp_var

◆ _finite_deform_Jacobian_mult

◆ _Jacobian_mult

◆ _ndisp

◆ _out_of_plane_direction

◆ _out_of_plane_strain

◆ _out_of_plane_strain_coupled

◆ _out_of_plane_strain_var

◆ _rotation_increment

◆ _stress

◆ _use_displaced_mesh

◆ _use_finite_deform_jacobian

◆ _volumetric_locking_correction