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 "PorosityFromStrain.h"
      11             : 
      12             : registerMooseObject("TensorMechanicsApp", PorosityFromStrain);
      13             : registerMooseObject("TensorMechanicsApp", ADPorosityFromStrain);
      14             : 
      15             : template <bool is_ad>
      16             : InputParameters
      17         200 : PorosityFromStrainTempl<is_ad>::validParams()
      18             : {
      19         200 :   InputParameters params = Material::validParams();
      20             : 
      21         200 :   params.addClassDescription("Porosity calculation from the inelastic strain.");
      22             : 
      23         400 :   params.addParam<MaterialPropertyName>(
      24             :       "porosity_name", "porosity", "Name of porosity material property");
      25         400 :   params.addRequiredRangeCheckedParam<Real>(
      26             :       "initial_porosity", "initial_porosity>=0.0 & initial_porosity<1.0", "Initial porosity");
      27         400 :   params.addRequiredParam<MaterialPropertyName>("inelastic_strain",
      28             :                                                 "Name of the combined inelastic strain");
      29         400 :   MooseEnum negative_behavior("ZERO INITIAL_CONDITION EXCEPTION", "INITIAL_CONDITION");
      30         400 :   params.addParam<MooseEnum>(
      31             :       "negative_behavior", negative_behavior, "Enum how to handle negative porosities");
      32             : 
      33         200 :   return params;
      34         200 : }
      35             : 
      36             : template <bool is_ad>
      37         150 : PorosityFromStrainTempl<is_ad>::PorosityFromStrainTempl(const InputParameters & parameters)
      38             :   : Material(parameters),
      39         150 :     _porosity(declareGenericProperty<Real, is_ad>(getParam<MaterialPropertyName>("porosity_name"))),
      40         300 :     _porosity_old(getMaterialPropertyOld<Real>("porosity_name")),
      41         300 :     _initial_porosity(getParam<Real>("initial_porosity")),
      42         150 :     _inelastic_strain(getGenericMaterialProperty<RankTwoTensor, is_ad>("inelastic_strain")),
      43         300 :     _inelastic_strain_old(getMaterialPropertyOld<RankTwoTensor>("inelastic_strain")),
      44         300 :     _negative_behavior(this->template getParam<MooseEnum>("negative_behavior")
      45         150 :                            .template getEnum<NegativeBehavior>())
      46             : {
      47         150 : }
      48             : 
      49             : template <bool is_ad>
      50             : void
      51         595 : PorosityFromStrainTempl<is_ad>::initQpStatefulProperties()
      52             : {
      53         595 :   _porosity[_qp] = _initial_porosity;
      54         595 : }
      55             : 
      56             : template <bool is_ad>
      57             : void
      58       41140 : PorosityFromStrainTempl<is_ad>::computeQpProperties()
      59             : {
      60       82228 :   _porosity[_qp] =
      61       82228 :       (1.0 - _porosity_old[_qp]) * (_inelastic_strain[_qp] - _inelastic_strain_old[_qp]).trace() +
      62       41140 :       _porosity_old[_qp];
      63             : 
      64       41140 :   if (_porosity[_qp] < 0.0)
      65             :   {
      66          25 :     if (_negative_behavior == NegativeBehavior::ZERO)
      67          12 :       _porosity[_qp] = 0.0;
      68          25 :     if (_negative_behavior == NegativeBehavior::INITIAL_CONDITION)
      69          12 :       _porosity[_qp] = _initial_porosity;
      70          25 :     if (_negative_behavior == NegativeBehavior::EXCEPTION)
      71           1 :       mooseException("In ", _name, ": porosity is negative.");
      72             :   }
      73             : 
      74       41139 :   if (std::isnan(_porosity[_qp]))
      75           1 :     mooseException("In ", _name, ": porosity is nan. Cutting timestep.");
      76       41138 : }