Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://www.mooseframework.org
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "ADComputeIncrementalStrainBase.h"
      11             : #include "MooseMesh.h"
      12             : #include "RankTwoTensor.h"
      13             : #include "RankFourTensor.h"
      14             : #include "SymmetricRankTwoTensor.h"
      15             : #include "SymmetricRankFourTensor.h"
      16             : 
      17             : template <typename R2>
      18             : InputParameters
      19           0 : ADComputeIncrementalStrainBaseTempl<R2>::validParams()
      20             : {
      21        1908 :   InputParameters params = ADComputeStrainBase::validParams();
      22           0 :   return params;
      23             : }
      24             : 
      25             : template <typename R2>
      26        1431 : ADComputeIncrementalStrainBaseTempl<R2>::ADComputeIncrementalStrainBaseTempl(
      27             :     const InputParameters & parameters)
      28             :   : ADComputeStrainBaseTempl<R2>(parameters),
      29        2862 :     _grad_disp_old(3),
      30        1431 :     _strain_rate(this->template declareADProperty<R2>(_base_name + "strain_rate")),
      31        1431 :     _strain_increment(this->template declareADProperty<R2>(_base_name + "strain_increment")),
      32        1431 :     _rotation_increment(
      33        1431 :         this->template declareADProperty<RankTwoTensor>(_base_name + "rotation_increment")),
      34        1431 :     _mechanical_strain_old(
      35        1431 :         this->template getMaterialPropertyOld<R2>(_base_name + "mechanical_strain")),
      36        2862 :     _total_strain_old(this->template getMaterialPropertyOld<R2>(_base_name + "total_strain")),
      37        2862 :     _eigenstrains_old(_eigenstrain_names.size())
      38             : {
      39        1530 :   for (unsigned int i = 0; i < _eigenstrains_old.size(); ++i)
      40          99 :     _eigenstrains_old[i] = &this->template getMaterialPropertyOld<R2>(_eigenstrain_names[i]);
      41        1431 : }
      42             : 
      43             : template <typename R2>
      44             : void
      45        1269 : ADComputeIncrementalStrainBaseTempl<R2>::initialSetup()
      46             : {
      47        1269 :   ADComputeStrainBaseTempl<R2>::initialSetup();
      48        5076 :   for (unsigned int i = 0; i < 3; ++i)
      49             :   {
      50        3807 :     if (this->_fe_problem.isTransient() && i < _ndisp)
      51        3087 :       _grad_disp_old[i] = &this->coupledGradientOld("displacements", i);
      52             :     else
      53         720 :       _grad_disp_old[i] = &_grad_zero;
      54             :   }
      55        1269 : }
      56             : 
      57             : template <typename R2>
      58             : void
      59       78576 : ADComputeIncrementalStrainBaseTempl<R2>::initQpStatefulProperties()
      60             : {
      61       78576 :   _mechanical_strain[_qp].zero();
      62       78576 :   _total_strain[_qp].zero();
      63       78576 : }
      64             : 
      65             : template <typename R2>
      66             : void
      67     8830441 : ADComputeIncrementalStrainBaseTempl<R2>::subtractEigenstrainIncrementFromStrain(ADR2 & strain)
      68             : {
      69     8974249 :   for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
      70             :   {
      71      143808 :     strain -= (*_eigenstrains[i])[_qp];
      72      143808 :     strain += (*_eigenstrains_old[i])[_qp];
      73             :   }
      74     8830441 : }
      75             : 
      76             : template class ADComputeIncrementalStrainBaseTempl<RankTwoTensor>;
      77             : template class ADComputeIncrementalStrainBaseTempl<SymmetricRankTwoTensor>;