ADStressDivergenceTensors< compute_stage > Class Template Reference

ADStressDivergenceTensors is the automatic differentiation version of StressDivergenceTensors. More...

#include <ADStressDivergenceTensors.h>

Inheritance diagram for ADStressDivergenceTensors< compute_stage >:

Public Member Functions
	ADStressDivergenceTensors (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
void	initialSetup () override
ADReal	computeQpResidual () override
void	precalculateResidual () override
const	ADMaterialProperty (RankTwoTensor) &_stress

Protected Attributes
const std::string	_base_name
const unsigned int	_component
const unsigned int	_ndisp
	Coupled displacement variables. More...
std::vector< unsigned int >	_disp_var
std::vector< ADReal >	_avg_grad_test
	Gradient of test function averaged over the element. Used in volumetric locking correction calculation. More...
const bool	_volumetric_locking_correction
	Flag for volumetric locking correction. More...
	usingKernelMembers

Detailed Description

template<ComputeStage compute_stage>
class ADStressDivergenceTensors< compute_stage >

ADStressDivergenceTensors is the automatic differentiation version of StressDivergenceTensors.

Definition at line 26 of file ADStressDivergenceTensors.h.

Constructor & Destructor Documentation

ADStressDivergenceTensors()

template<ComputeStage compute_stage>

ADStressDivergenceTensors< compute_stage >::ADStressDivergenceTensors

(

const InputParameters &

parameters

)

Definition at line 40 of file ADStressDivergenceTensors.C.

  : ADKernel<compute_stage>(parameters),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
    _stress(getADMaterialProperty<RankTwoTensor>(_base_name + "stress")),
    _component(getParam<unsigned int>("component")),
    _ndisp(coupledComponents("displacements")),
    _disp_var(_ndisp),
    _avg_grad_test(),
    _volumetric_locking_correction(getParam<bool>("volumetric_locking_correction"))
{
  for (unsigned int i = 0; i < _ndisp; ++i)
    // the next line should be _disp_var[i] = coupled("displacements", i);
    // but the Coupleable:: is required to avoid triggering an internal Intel compiler bug
    _disp_var[i] = Coupleable::coupled("displacements", i);
 
  // 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.");
}

Member Function Documentation

ADMaterialProperty()

template<ComputeStage compute_stage>

const ADStressDivergenceTensors< compute_stage >::ADMaterialProperty ( RankTwoTensor  ) &

protected

computeQpResidual()

template<ComputeStage compute_stage>

ADReal ADStressDivergenceTensors< compute_stage >::computeQpResidual ( )

overrideprotected

Definition at line 72 of file ADStressDivergenceTensors.C.

{
  ADReal residual = _stress[_qp].row(_component) * _grad_test[_i][_qp];
 
  // volumetric locking correction
  if (_volumetric_locking_correction)
    residual += (_avg_grad_test[_i] - _grad_test[_i][_qp](_component)) / 3.0 * _stress[_qp].trace();
 
  return residual;
}

initialSetup()

template<ComputeStage compute_stage>

void ADStressDivergenceTensors< compute_stage >::initialSetup ( )

overrideprotected

Definition at line 63 of file ADStressDivergenceTensors.C.

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

precalculateResidual()

template<ComputeStage compute_stage>

void ADStressDivergenceTensors< compute_stage >::precalculateResidual ( )

overrideprotected

Definition at line 85 of file ADStressDivergenceTensors.C.

{
  if (!_volumetric_locking_correction)
    return;
 
  ADReal ad_current_elem_volume = 0.0;
  for (unsigned int qp = 0; qp < _qrule->n_points(); qp++)
    ad_current_elem_volume += _ad_JxW[qp] * _ad_coord[qp];
 
  // 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] = 0.0;
    for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
      _avg_grad_test[_i] += _grad_test[_i][_qp](_component) * _ad_JxW[_qp] * _ad_coord[_qp];
 
    _avg_grad_test[_i] /= ad_current_elem_volume;
  }
}

validParams()

template<ComputeStage compute_stage>

InputParameters ADStressDivergenceTensors< compute_stage >::validParams ( )

static

Definition at line 20 of file ADStressDivergenceTensors.C.

{
  InputParameters params = ADKernel<compute_stage>::validParams();
  params.addClassDescription("Stress divergence kernel with automatic differentiation 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.addParam<std::string>("base_name", "Material property base name");
  params.set<bool>("use_displaced_mesh") = false;
  params.addParam<bool>("volumetric_locking_correction",
                        false,
                        "Set to false to turn off volumetric locking correction");
  return params;
}

Referenced by ADDynamicStressDivergenceTensors< compute_stage >::validParams(), ADStressDivergenceRZTensors< compute_stage >::validParams(), and ADStressDivergenceRSphericalTensors< compute_stage >::validParams().

Member Data Documentation

_avg_grad_test

template<ComputeStage compute_stage>

std::vector<ADReal> ADStressDivergenceTensors< compute_stage >::_avg_grad_test

protected

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

Definition at line 61 of file ADStressDivergenceTensors.h.

_base_name

template<ComputeStage compute_stage>

const std::string ADStressDivergenceTensors< compute_stage >::_base_name

protected

Definition at line 51 of file ADStressDivergenceTensors.h.

_component

template<ComputeStage compute_stage>

const unsigned int ADStressDivergenceTensors< compute_stage >::_component

protected

Definition at line 54 of file ADStressDivergenceTensors.h.

Referenced by ADStressDivergenceRSphericalTensors< compute_stage >::ADStressDivergenceRSphericalTensors().

_disp_var

template<ComputeStage compute_stage>

std::vector<unsigned int> ADStressDivergenceTensors< compute_stage >::_disp_var

protected

Definition at line 58 of file ADStressDivergenceTensors.h.

Referenced by ADStressDivergenceTensors< compute_stage >::ADStressDivergenceTensors().

_ndisp

template<ComputeStage compute_stage>

const unsigned int ADStressDivergenceTensors< compute_stage >::_ndisp

protected

Coupled displacement variables.

Definition at line 57 of file ADStressDivergenceTensors.h.

Referenced by ADStressDivergenceTensors< compute_stage >::ADStressDivergenceTensors().

_volumetric_locking_correction

template<ComputeStage compute_stage>

const bool ADStressDivergenceTensors< compute_stage >::_volumetric_locking_correction

protected

Flag for volumetric locking correction.

Definition at line 64 of file ADStressDivergenceTensors.h.

Referenced by ADStressDivergenceTensors< compute_stage >::ADStressDivergenceTensors().

usingKernelMembers

template<ComputeStage compute_stage>

ADStressDivergenceTensors< compute_stage >::usingKernelMembers

protected

Definition at line 66 of file ADStressDivergenceTensors.h.

---

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

• tensor_mechanics/include/kernels/ADStressDivergenceTensors.h
• tensor_mechanics/src/kernels/ADStressDivergenceTensors.C