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 "ADComputeStrainIncrementBasedStress.h"
      11             : 
      12             : registerMooseObject("TensorMechanicsApp", ADComputeStrainIncrementBasedStress);
      13             : 
      14             : InputParameters
      15          12 : ADComputeStrainIncrementBasedStress::validParams()
      16             : {
      17          12 :   InputParameters params = ADComputeStressBase::validParams();
      18          12 :   params.addClassDescription("Compute stress after subtracting inelastic strain increments");
      19          24 :   params.addParam<std::vector<MaterialPropertyName>>(
      20             :       "inelastic_strain_names", {}, "Names of inelastic strain properties");
      21             : 
      22          12 :   return params;
      23           0 : }
      24             : 
      25           9 : ADComputeStrainIncrementBasedStress::ADComputeStrainIncrementBasedStress(
      26           9 :     const InputParameters & parameters)
      27             :   : ADComputeStressBase(parameters),
      28           9 :     _elasticity_tensor_name(_base_name + "elasticity_tensor"),
      29           9 :     _elasticity_tensor(getADMaterialPropertyByName<RankFourTensor>(_elasticity_tensor_name)),
      30          18 :     _stress_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "stress")),
      31           9 :     _mechanical_strain_old(
      32           9 :         getMaterialPropertyOldByName<RankTwoTensor>(_base_name + "mechanical_strain")),
      33          36 :     _inelastic_strain_names(getParam<std::vector<MaterialPropertyName>>("inelastic_strain_names"))
      34             : {
      35           9 :   _num_inelastic_strain_models = _inelastic_strain_names.size();
      36             : 
      37           9 :   if (_num_inelastic_strain_models > 0)
      38             :   {
      39           0 :     _inelastic_strains.resize(_num_inelastic_strain_models);
      40           0 :     _inelastic_strains_old.resize(_num_inelastic_strain_models);
      41             : 
      42           0 :     for (unsigned int i = 0; i < _num_inelastic_strain_models; ++i)
      43             :     {
      44           0 :       _inelastic_strains[i] = &getADMaterialProperty<RankTwoTensor>(_inelastic_strain_names[i]);
      45           0 :       _inelastic_strains_old[i] =
      46           0 :           &getMaterialPropertyOld<RankTwoTensor>(_inelastic_strain_names[i]);
      47             :     }
      48             :   }
      49           9 : }
      50             : 
      51             : void
      52         960 : ADComputeStrainIncrementBasedStress::computeQpStress()
      53             : {
      54             :   ADRankTwoTensor elastic_strain_increment =
      55         960 :       (_mechanical_strain[_qp] - _mechanical_strain_old[_qp]);
      56             : 
      57         960 :   for (unsigned int i = 0; i < _num_inelastic_strain_models; ++i)
      58           0 :     elastic_strain_increment -= (*_inelastic_strains[i])[_qp] - (*_inelastic_strains_old[i])[_qp];
      59             : 
      60         960 :   _stress[_qp] = _stress_old[_qp] + _elasticity_tensor[_qp] * elastic_strain_increment;
      61         960 : }