Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       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 "SolidMechanicsHardeningGaussian.h"
      11             : 
      12             : registerMooseObject("SolidMechanicsApp", SolidMechanicsHardeningGaussian);
      13             : registerMooseObjectRenamed("SolidMechanicsApp",
      14             :                            TensorMechanicsHardeningGaussian,
      15             :                            "01/01/2025 00:00",
      16             :                            SolidMechanicsHardeningGaussian);
      17             : 
      18             : InputParameters
      19          32 : SolidMechanicsHardeningGaussian::validParams()
      20             : {
      21          32 :   InputParameters params = SolidMechanicsHardeningModel::validParams();
      22          64 :   params.addRequiredParam<Real>(
      23             :       "value_0", "The value of the parameter for all internal_parameter <= internal_0");
      24          64 :   params.addParam<Real>("value_residual",
      25             :                         "The value of the parameter for internal_parameter = "
      26             :                         "infinity.  Default = value_0, ie perfect plasticity");
      27          64 :   params.addParam<Real>(
      28          64 :       "internal_0", 0, "The value of the internal_parameter when hardening begins");
      29          64 :   params.addParam<Real>("rate",
      30          64 :                         0,
      31             :                         "Let p = internal_parameter.  Then value = value_0 for "
      32             :                         "p<internal_0, and value = value_residual + (value_0 - "
      33             :                         "value_residual)*exp(-0.5*rate*(p - internal_0)^2)");
      34          32 :   params.addClassDescription("Hardening is Gaussian");
      35          32 :   return params;
      36           0 : }
      37             : 
      38          16 : SolidMechanicsHardeningGaussian::SolidMechanicsHardeningGaussian(const InputParameters & parameters)
      39             :   : SolidMechanicsHardeningModel(parameters),
      40          16 :     _val_0(getParam<Real>("value_0")),
      41          48 :     _val_res(parameters.isParamValid("value_residual") ? getParam<Real>("value_residual") : _val_0),
      42          32 :     _intnl_0(getParam<Real>("internal_0")),
      43          48 :     _rate(getParam<Real>("rate"))
      44             : {
      45          16 : }
      46             : 
      47             : Real
      48       32396 : SolidMechanicsHardeningGaussian::value(Real intnl) const
      49             : {
      50       32396 :   Real x = intnl - _intnl_0;
      51       32396 :   if (x <= 0)
      52        1460 :     return _val_0;
      53             :   else
      54       30936 :     return _val_res + (_val_0 - _val_res) * std::exp(-0.5 * _rate * x * x);
      55             : }
      56             : 
      57             : Real
      58       12072 : SolidMechanicsHardeningGaussian::derivative(Real intnl) const
      59             : {
      60       12072 :   Real x = intnl - _intnl_0;
      61       12072 :   if (x <= 0)
      62             :     return 0;
      63             :   else
      64       11760 :     return -_rate * x * (_val_0 - _val_res) * std::exp(-0.5 * _rate * x * x);
      65             : }
      66             : 
      67             : std::string
      68           0 : SolidMechanicsHardeningGaussian::modelName() const
      69             : {
      70           0 :   return "Gaussian";
      71             : }