ADComputeSmallStrain.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "ADComputeSmallStrain.h"
#include "RankTwoTensor.h"
#include "SymmetricRankTwoTensor.h"
#include "libmesh/quadrature.h"

registerMooseObject("SolidMechanicsApp", ADComputeSmallStrain);
registerMooseObject("SolidMechanicsApp", ADSymmetricSmallStrain);

template <typename R2>
InputParameters
ADComputeSmallStrainTempl<R2>::validParams()
{
  InputParameters params = ADComputeStrainBase::validParams();
  params.addClassDescription("Compute a small strain.");
  return params;
}

template <typename R2>
ADComputeSmallStrainTempl<R2>::ADComputeSmallStrainTempl(const InputParameters & parameters)
  : ADComputeStrainBaseTempl<R2>(parameters)
{
}

template <typename R2>
void
ADComputeSmallStrainTempl<R2>::computeProperties()
{
  ADReal volumetric_strain = 0.0;

  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
  {
    // strain = (grad_disp + grad_disp^R2)/2
    _total_strain[_qp] = Moose::ADType<R2>::type::initializeSymmetric(
        (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]);

    if (_volumetric_locking_correction)
      volumetric_strain += _total_strain[_qp].trace() * _JxW[_qp] * _coord[_qp];
  }

  if (_volumetric_locking_correction)
    volumetric_strain /= _current_elem_volume;

  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
  {
    if (_volumetric_locking_correction)
    {
      ADReal correction = (volumetric_strain - _total_strain[_qp].trace()) / 3.0;
      _total_strain[_qp].addIa(correction);
    }

    if (_global_strain)
      _total_strain[_qp] += (*_global_strain)[_qp];

    _mechanical_strain[_qp] = _total_strain[_qp];

    // Remove the Eigen strain
    for (auto es : _eigenstrains)
      _mechanical_strain[_qp] -= (*es)[_qp];
  }
}

template class ADComputeSmallStrainTempl<RankTwoTensor>;
template class ADComputeSmallStrainTempl<SymmetricRankTwoTensor>;