ADComputeMultiplePorousInelasticStress.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 "ADComputeMultiplePorousInelasticStress.h"
 
#include "RankTwoTensor.h"
 
registerMooseObjectDeprecated("SolidMechanicsApp",
                              ADComputeMultiplePorousInelasticStress,
                              "01/30/2021 24:00");
 
InputParameters
ADComputeMultiplePorousInelasticStress::validParams()
{
  InputParameters params = ADComputeMultipleInelasticStress::validParams();
  params.addClassDescription(
      "Compute state (stress and internal parameters such as plastic "
      "strains and internal parameters) using an iterative process. A porosity material property "
      "is defined and is calculated from the trace of inelastic strain increment.");
 
  params.addParam<MaterialPropertyName>(
      "porosity_name", "porosity", "Name of porosity material property");
  params.addRequiredRangeCheckedParam<Real>(
      "initial_porosity", "initial_porosity>0.0 & initial_porosity<1.0", "Initial porosity");
  return params;
}
 
ADComputeMultiplePorousInelasticStress::ADComputeMultiplePorousInelasticStress(
    const InputParameters & parameters)
  : ADComputeMultipleInelasticStress(parameters),
    _porosity(declareADProperty<Real>(getParam<MaterialPropertyName>("porosity_name"))),
    _porosity_old(getMaterialPropertyOld<Real>(getParam<MaterialPropertyName>("porosity_name"))),
    _initial_porosity(getParam<Real>("initial_porosity"))
{
}
 
void
ADComputeMultiplePorousInelasticStress::initQpStatefulProperties()
{
  ADComputeStressBase::initQpStatefulProperties();
 
  _porosity[_qp] = _initial_porosity;
}
 
void
ADComputeMultiplePorousInelasticStress::computeQpProperties()
{
  ADComputeStressBase::computeQpProperties();
 
  _porosity[_qp] =
      (1.0 - _porosity_old[_qp]) * (_inelastic_strain[_qp] - _inelastic_strain_old[_qp]).trace() +
      _porosity_old[_qp];
  if (_porosity[_qp] < 0.0)
    _porosity[_qp] = 0.0;
}