ADStressDivergenceShell< compute_stage > Class Template Reference

ADStressDivergenceShell computes the stress divergence term for shell elements. More...

#include <ADStressDivergenceShell.h>

Inheritance diagram for ADStressDivergenceShell< compute_stage >:

Public Member Functions
	ADStressDivergenceShell (const InputParameters &parameters)

Protected Member Functions
virtual ADReal	computeQpResidual () override

Protected Attributes
const unsigned int	_component
const bool	_large_strain
std::vector< const ADMaterialProperty(RankTwoTensor) * >	_stress
std::vector< const MaterialProperty< RankTwoTensor > * >	_stress_old
std::vector< const ADMaterialProperty(DenseMatrix< Real >) * >	_B_mat
std::vector< const ADMaterialProperty(DenseMatrix< Real >) * >	_B_nl
std::vector< const ADMaterialProperty(Real) * >	_J_map
std::unique_ptr< QGauss >	_t_qrule
	Quadrature rule in the out of plane direction. More...
std::vector< Real >	_t_weights
	Quadrature weights in the out of plane direction in isoparametric coordinate system. More...
std::vector< Real >	_q_weights
	Qrule weights in isoparametric coordinate system. More...
unsigned int	_qp_z
	qp index in out of plane direction More...
	usingKernelMembers

Detailed Description

template<ComputeStage compute_stage>
class ADStressDivergenceShell< compute_stage >

ADStressDivergenceShell computes the stress divergence term for shell elements.

Definition at line 19 of file ADStressDivergenceShell.h.

Constructor & Destructor Documentation

ADStressDivergenceShell()

template<ComputeStage compute_stage>

ADStressDivergenceShell< compute_stage >::ADStressDivergenceShell

(

const InputParameters &

parameters

)

Definition at line 45 of file ADStressDivergenceShell.C.

  : ADKernel<compute_stage>(parameters),
    _component(getParam<unsigned int>("component")),
    _large_strain(getParam<bool>("large_strain"))
{
  _t_qrule = libmesh_make_unique<QGauss>(
      1, Utility::string_to_enum<Order>(getParam<std::string>("through_thickness_order")));
  _t_weights = _t_qrule->get_weights();
 
  _stress.resize(_t_weights.size());
  _stress_old.resize(_t_weights.size());
  _B_mat.resize(_t_weights.size());
  if (_large_strain)
    _B_nl.resize(_t_weights.size());
  _J_map.resize(_t_weights.size());
 
  for (unsigned int i = 0; i < _t_weights.size(); ++i)
  {
    _stress[i] = &getADMaterialProperty<RankTwoTensor>("stress_t_points_" + std::to_string(i));
    _stress_old[i] =
        &getMaterialPropertyOldByName<RankTwoTensor>("stress_t_points_" + std::to_string(i));
    _B_mat[i] = &getADMaterialProperty<DenseMatrix<Real>>("B_t_points_" + std::to_string(i));
    if (_large_strain)
      _B_nl[i] = &getADMaterialProperty<DenseMatrix<Real>>("B_nl_t_points_" + std::to_string(i));
 
    _J_map[i] = &getADMaterialProperty<Real>("J_mapping_t_points_" + std::to_string(i));
  }
}

Member Function Documentation

computeQpResidual()

template<ComputeStage compute_stage>

ADReal ADStressDivergenceShell< compute_stage >::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 76 of file ADStressDivergenceShell.C.

{
  _q_weights = _qrule->get_weights();
  ADReal residual = 0.0;
  ADReal residual1 = 0.0;
  for (_qp_z = 0; _qp_z < _t_weights.size(); ++_qp_z)
  {
    residual1 = (*_stress[_qp_z])[_qp](0, 0) * (*_B_mat[_qp_z])[_qp](0, _i + _component * 4) +
                (*_stress[_qp_z])[_qp](1, 1) * (*_B_mat[_qp_z])[_qp](1, _i + _component * 4) +
                2.0 * (*_stress[_qp_z])[_qp](0, 1) * (*_B_mat[_qp_z])[_qp](2, _i + _component * 4) +
                2.0 * (*_stress[_qp_z])[_qp](0, 2) * (*_B_mat[_qp_z])[_qp](3, _i + _component * 4) +
                2.0 * (*_stress[_qp_z])[_qp](1, 2) * (*_B_mat[_qp_z])[_qp](4, _i + _component * 4);
 
    if (_large_strain)
      residual1 +=
          ((*_stress_old[_qp_z])[_qp](0, 0) * (*_B_nl[_qp_z])[_qp](0, _i + _component * 4) +
           (*_stress_old[_qp_z])[_qp](1, 1) * (*_B_nl[_qp_z])[_qp](1, _i + _component * 4) +
           2.0 * (*_stress_old[_qp_z])[_qp](0, 1) * (*_B_nl[_qp_z])[_qp](2, _i + _component * 4) +
           2.0 * (*_stress_old[_qp_z])[_qp](0, 2) * (*_B_nl[_qp_z])[_qp](3, _i + _component * 4) +
           2.0 * (*_stress_old[_qp_z])[_qp](1, 2) * (*_B_nl[_qp_z])[_qp](4, _i + _component * 4));
 
    residual += residual1 * (*_J_map[_qp_z])[_qp] * _q_weights[_qp] * _t_weights[_qp_z] /
                (_ad_JxW[_qp] * _ad_coord[_qp]);
  }
  return residual;
}

Member Data Documentation

_B_mat

template<ComputeStage compute_stage>

std::vector<const ADMaterialProperty(DenseMatrix<Real>) *> ADStressDivergenceShell< compute_stage >::_B_mat

protected

Definition at line 54 of file ADStressDivergenceShell.h.

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

_B_nl

template<ComputeStage compute_stage>

std::vector<const ADMaterialProperty(DenseMatrix<Real>) *> ADStressDivergenceShell< compute_stage >::_B_nl

protected

Definition at line 55 of file ADStressDivergenceShell.h.

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

_component

template<ComputeStage compute_stage>

const unsigned int ADStressDivergenceShell< compute_stage >::_component

protected

Definition at line 49 of file ADStressDivergenceShell.h.

_J_map

template<ComputeStage compute_stage>

std::vector<const ADMaterialProperty(Real) *> ADStressDivergenceShell< compute_stage >::_J_map

protected

Definition at line 56 of file ADStressDivergenceShell.h.

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

_large_strain

template<ComputeStage compute_stage>

const bool ADStressDivergenceShell< compute_stage >::_large_strain

protected

Definition at line 50 of file ADStressDivergenceShell.h.

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

_q_weights

template<ComputeStage compute_stage>

std::vector<Real> ADStressDivergenceShell< compute_stage >::_q_weights

protected

Qrule weights in isoparametric coordinate system.

Definition at line 65 of file ADStressDivergenceShell.h.

_qp_z

template<ComputeStage compute_stage>

unsigned int ADStressDivergenceShell< compute_stage >::_qp_z

protected

qp index in out of plane direction

Definition at line 68 of file ADStressDivergenceShell.h.

_stress

template<ComputeStage compute_stage>

std::vector<const ADMaterialProperty(RankTwoTensor) *> ADStressDivergenceShell< compute_stage >::_stress

protected

Definition at line 52 of file ADStressDivergenceShell.h.

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

_stress_old

template<ComputeStage compute_stage>

std::vector<const MaterialProperty<RankTwoTensor> *> ADStressDivergenceShell< compute_stage >::_stress_old

protected

Definition at line 53 of file ADStressDivergenceShell.h.

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

_t_qrule

template<ComputeStage compute_stage>

std::unique_ptr<QGauss> ADStressDivergenceShell< compute_stage >::_t_qrule

protected

Quadrature rule in the out of plane direction.

Definition at line 59 of file ADStressDivergenceShell.h.

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

_t_weights

template<ComputeStage compute_stage>

std::vector<Real> ADStressDivergenceShell< compute_stage >::_t_weights

protected

Quadrature weights in the out of plane direction in isoparametric coordinate system.

Definition at line 62 of file ADStressDivergenceShell.h.

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

usingKernelMembers

template<ComputeStage compute_stage>

ADStressDivergenceShell< compute_stage >::usingKernelMembers

protected

Definition at line 70 of file ADStressDivergenceShell.h.

---

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

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