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 "libmesh/quadrature.h"
 
registerADMooseObject("TensorMechanicsApp", ADComputeSmallStrain);
 
defineADLegacyParams(ADComputeSmallStrain);
 
template <ComputeStage compute_stage>
InputParameters
ADComputeSmallStrain<compute_stage>::validParams()
{
  InputParameters params = ADComputeStrainBase<compute_stage>::validParams();
  params.addClassDescription("Compute a small strain.");
  return params;
}
 
template <ComputeStage compute_stage>
ADComputeSmallStrain<compute_stage>::ADComputeSmallStrain(const InputParameters & parameters)
  : ADComputeStrainBase<compute_stage>(parameters)
{
}
 
template <ComputeStage compute_stage>
void
ADComputeSmallStrain<compute_stage>::computeProperties()
{
  ADReal volumetric_strain = 0.0;
 
  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
  {
    // strain = (grad_disp + grad_disp^T)/2
    ADRankTwoTensor grad_tensor(
        (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]);
 
    _total_strain[_qp] = (grad_tensor + grad_tensor.transpose()) / 2.0;
 
    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](0, 0) += correction;
      _total_strain[_qp](1, 1) += correction;
      _total_strain[_qp](2, 2) += 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];
  }
 
  copyDualNumbersToValues();
}
 
// explicit instantiation is required for AD base classes
adBaseClass(ADComputeSmallStrain);