IsotropicPlasticityStressUpdate.h

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

#pragma once

#include "RadialReturnStressUpdate.h"

template <bool is_ad>
class IsotropicPlasticityStressUpdateTempl : public RadialReturnStressUpdateTempl<is_ad>
{
public:
  static InputParameters validParams();

  IsotropicPlasticityStressUpdateTempl(const InputParameters & parameters);

  using Material::_qp;
  using RadialReturnStressUpdateTempl<is_ad>::_base_name;
  using RadialReturnStressUpdateTempl<is_ad>::_three_shear_modulus;

  virtual void
  computeStressInitialize(const GenericReal<is_ad> & effective_trial_stress,
                          const GenericRankFourTensor<is_ad> & elasticity_tensor) override;
  virtual GenericReal<is_ad> computeResidual(const GenericReal<is_ad> & effective_trial_stress,
                                             const GenericReal<is_ad> & scalar) override;
  virtual GenericReal<is_ad> computeDerivative(const GenericReal<is_ad> & effective_trial_stress,
                                               const GenericReal<is_ad> & scalar) override;

  virtual void computeYieldStress(const GenericRankFourTensor<is_ad> & elasticity_tensor);

  GenericReal<is_ad> yieldCondition() const { return _yield_condition; }

  virtual void
  computeStressFinalize(const GenericRankTwoTensor<is_ad> & plastic_strain_increment) override;

protected:
  virtual void initQpStatefulProperties() override;
  virtual void propagateQpStatefulProperties() override;

  virtual void iterationFinalize(const GenericReal<is_ad> & scalar) override;

  virtual GenericReal<is_ad> computeHardeningValue(const GenericReal<is_ad> & scalar);
  virtual GenericReal<is_ad> computeHardeningDerivative(const GenericReal<is_ad> & scalar);

  const std::string _plastic_prepend;

  const Function * _yield_stress_function;
  GenericReal<is_ad> _yield_stress;
  const Real _hardening_constant;
  const Function * const _hardening_function;

  GenericReal<is_ad> _yield_condition;
  GenericReal<is_ad> _hardening_slope;

  GenericMaterialProperty<RankTwoTensor, is_ad> & _plastic_strain;

  const MaterialProperty<RankTwoTensor> & _plastic_strain_old;

  GenericMaterialProperty<Real, is_ad> & _hardening_variable;
  const MaterialProperty<Real> & _hardening_variable_old;
  const GenericVariableValue<is_ad> & _temperature;
};

typedef IsotropicPlasticityStressUpdateTempl<false> IsotropicPlasticityStressUpdate;
typedef IsotropicPlasticityStressUpdateTempl<true> ADIsotropicPlasticityStressUpdate;