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 "ComputeStrainIncrementBasedStress.h"
      11             : 
      12             : registerMooseObject("TensorMechanicsApp", ComputeStrainIncrementBasedStress);
      13             : 
      14             : InputParameters
      15         132 : ComputeStrainIncrementBasedStress::validParams()
      16             : {
      17         132 :   InputParameters params = ComputeStressBase::validParams();
      18         132 :   params.addClassDescription("Compute stress after subtracting inelastic strain increments");
      19         264 :   params.addParam<std::vector<MaterialPropertyName>>(
      20             :       "inelastic_strain_names", {}, "Names of inelastic strain properties");
      21             : 
      22         132 :   return params;
      23           0 : }
      24             : 
      25          99 : ComputeStrainIncrementBasedStress::ComputeStrainIncrementBasedStress(
      26          99 :     const InputParameters & parameters)
      27             :   : ComputeStressBase(parameters),
      28          99 :     _elasticity_tensor_name(_base_name + "elasticity_tensor"),
      29          99 :     _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_elasticity_tensor_name)),
      30         198 :     _stress_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "stress")),
      31          99 :     _mechanical_strain_old(
      32          99 :         getMaterialPropertyOldByName<RankTwoTensor>(_base_name + "mechanical_strain")),
      33         396 :     _inelastic_strain_names(getParam<std::vector<MaterialPropertyName>>("inelastic_strain_names"))
      34             : {
      35          99 :   _num_inelastic_strain_models = _inelastic_strain_names.size();
      36             : 
      37          99 :   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] = &getMaterialProperty<RankTwoTensor>(_inelastic_strain_names[i]);
      45           0 :       _inelastic_strains_old[i] =
      46           0 :           &getMaterialPropertyOld<RankTwoTensor>(_inelastic_strain_names[i]);
      47             :     }
      48             :   }
      49          99 : }
      50             : 
      51             : void
      52       27800 : ComputeStrainIncrementBasedStress::computeQpStress()
      53             : {
      54       27800 :   RankTwoTensor elastic_strain_increment = (_mechanical_strain[_qp] - _mechanical_strain_old[_qp]);
      55             : 
      56       27800 :   for (unsigned int i = 0; i < _num_inelastic_strain_models; ++i)
      57           0 :     elastic_strain_increment -= (*_inelastic_strains[i])[_qp] - (*_inelastic_strains_old[i])[_qp];
      58             : 
      59       27800 :   _stress[_qp] = _stress_old[_qp] + _elasticity_tensor[_qp] * elastic_strain_increment;
      60             : 
      61       27800 :   computeQpJacobian();
      62       27800 : }
      63             : 
      64             : void
      65       27800 : ComputeStrainIncrementBasedStress::computeQpJacobian()
      66             : {
      67       27800 :   _Jacobian_mult[_qp] = _elasticity_tensor[_qp];
      68       27800 : }