SolidMechanicsHardeningGaussian.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 "SolidMechanicsHardeningGaussian.h"

registerMooseObject("SolidMechanicsApp", SolidMechanicsHardeningGaussian);
registerMooseObjectRenamed("SolidMechanicsApp",
                           TensorMechanicsHardeningGaussian,
                           "01/01/2025 00:00",
                           SolidMechanicsHardeningGaussian);

InputParameters
SolidMechanicsHardeningGaussian::validParams()
{
  InputParameters params = SolidMechanicsHardeningModel::validParams();
  params.addRequiredParam<Real>(
      "value_0", "The value of the parameter for all internal_parameter <= internal_0");
  params.addParam<Real>("value_residual",
                        "The value of the parameter for internal_parameter = "
                        "infinity.  Default = value_0, ie perfect plasticity");
  params.addParam<Real>(
      "internal_0", 0, "The value of the internal_parameter when hardening begins");
  params.addParam<Real>("rate",
                        0,
                        "Let p = internal_parameter.  Then value = value_0 for "
                        "p<internal_0, and value = value_residual + (value_0 - "
                        "value_residual)*exp(-0.5*rate*(p - internal_0)^2)");
  params.addClassDescription("Hardening is Gaussian");
  return params;
}

SolidMechanicsHardeningGaussian::SolidMechanicsHardeningGaussian(const InputParameters & parameters)
  : SolidMechanicsHardeningModel(parameters),
    _val_0(getParam<Real>("value_0")),
    _val_res(parameters.isParamValid("value_residual") ? getParam<Real>("value_residual") : _val_0),
    _intnl_0(getParam<Real>("internal_0")),
    _rate(getParam<Real>("rate"))
{
}

Real
SolidMechanicsHardeningGaussian::value(Real intnl) const
{
  Real x = intnl - _intnl_0;
  if (x <= 0)
    return _val_0;
  else
    return _val_res + (_val_0 - _val_res) * std::exp(-0.5 * _rate * x * x);
}

Real
SolidMechanicsHardeningGaussian::derivative(Real intnl) const
{
  Real x = intnl - _intnl_0;
  if (x <= 0)
    return 0;
  else
    return -_rate * x * (_val_0 - _val_res) * std::exp(-0.5 * _rate * x * x);
}

std::string
SolidMechanicsHardeningGaussian::modelName() const
{
  return "Gaussian";
}