CombinedNonlinearHardeningPlasticity.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 "MaterialProperty.h"
#include "RadialReturnBackstressStressUpdateBase.h"
 
template <bool is_ad>
class CombinedNonlinearHardeningPlasticityTempl
  : public RadialReturnBackstressStressUpdateBaseTempl<is_ad>
{
public:
  static InputParameters validParams();
 
  CombinedNonlinearHardeningPlasticityTempl(const InputParameters & parameters);
 
  using Material::_qp;
  using RadialReturnBackstressStressUpdateBaseTempl<is_ad>::_base_name;
  using RadialReturnBackstressStressUpdateBaseTempl<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> computeIsotropicHardeningValue(const GenericReal<is_ad> & scalar);
  virtual GenericReal<is_ad> computeIsotropicHardeningDerivative(const GenericReal<is_ad> & scalar);
  virtual GenericReal<is_ad> computeKinematicHardeningValue(const GenericReal<is_ad> & scalar);
 
  const Function * _yield_stress_function;
  GenericReal<is_ad> _yield_stress;
  const Real _isotropic_hardening_constant;
  const Function * const _isotropic_hardening_function;
 
  GenericReal<is_ad> _yield_condition;
  GenericReal<is_ad> _isotropic_hardening_slope;
  GenericReal<is_ad> _kinematic_hardening_slope;
  GenericRankTwoTensor<is_ad> stress_new;
  GenericRankTwoTensor<is_ad> elastic_strain_old;
  GenericRankTwoTensor<is_ad> strain_increment;
 
  GenericMaterialProperty<RankTwoTensor, is_ad> & _plastic_strain;
 
  const MaterialProperty<RankTwoTensor> & _plastic_strain_old;
 
  const Real _kinematic_hardening_modulus;
  const Real _gamma;
  const Real _q;
  const Real _b;
 
  GenericMaterialProperty<Real, is_ad> & _isotropic_hardening_variable;
  const MaterialProperty<Real> & _isotropic_hardening_variable_old;
  GenericMaterialProperty<Real, is_ad> & _kinematic_hardening_variable;
  const MaterialProperty<Real> & _kinematic_hardening_variable_old;
  const GenericVariableValue<is_ad> & _temperature;
};
 
typedef CombinedNonlinearHardeningPlasticityTempl<false> CombinedNonlinearHardeningPlasticity;
typedef CombinedNonlinearHardeningPlasticityTempl<true> ADCombinedNonlinearHardeningPlasticity;