ADCompute1DIncrementalStrain.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 "ADCompute1DIncrementalStrain.h"
 
#include "libmesh/quadrature.h"
 
defineADLegacyParams(ADCompute1DIncrementalStrain);
 
template <ComputeStage compute_stage>
InputParameters
ADCompute1DIncrementalStrain<compute_stage>::validParams()
{
  InputParameters params = ADComputeIncrementalSmallStrain<compute_stage>::validParams();
  params.addClassDescription("Compute strain increment for small strains in 1D problems.");
  return params;
}
 
template <ComputeStage compute_stage>
ADCompute1DIncrementalStrain<compute_stage>::ADCompute1DIncrementalStrain(
    const InputParameters & parameters)
  : ADComputeIncrementalSmallStrain<compute_stage>(parameters)
{
}
 
template <ComputeStage compute_stage>
void
ADCompute1DIncrementalStrain<compute_stage>::computeTotalStrainIncrement(
    ADRankTwoTensor & total_strain_increment)
{
  // Deformation gradient calculation for 1D problems
  ADRankTwoTensor A(
      (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]); // Deformation gradient
  RankTwoTensor Fbar((*_grad_disp_old[0])[_qp],
                     (*_grad_disp_old[1])[_qp],
                     (*_grad_disp_old[2])[_qp]); // Old Deformation gradient
 
  // Compute the displacement gradient dUy/dy and dUz/dz value for 1D problems
  A(1, 1) = computeGradDispYY();
  A(2, 2) = computeGradDispZZ();
 
  Fbar(1, 1) = computeGradDispYYOld();
  Fbar(2, 2) = computeGradDispZZOld();
 
  A -= Fbar; // very nearly A = gradU - gradUold, adapted to cylindrical coords
 
  total_strain_increment = 0.5 * (A + A.transpose());
}
 
// explicit instantiation is required for AD base classes
adBaseClass(ADCompute1DIncrementalStrain);