ADMultiplePowerLawCreepStressUpdate.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 "RadialReturnCreepStressUpdateBase.h"

class ADMultiplePowerLawCreepStressUpdate : public ADRadialReturnCreepStressUpdateBase
{
public:
  static InputParameters validParams();

  ADMultiplePowerLawCreepStressUpdate(const InputParameters & parameters);

  virtual Real
  computeStrainEnergyRateDensity(const ADMaterialProperty<RankTwoTensor> & stress,
                                 const ADMaterialProperty<RankTwoTensor> & strain_rate) override;

  virtual bool substeppingCapabilityEnabled() override;

protected:
  virtual void computeStressInitialize(const ADReal & effective_trial_stress,
                                       const ADRankFourTensor & elasticity_tensor) override;
  virtual ADReal computeResidual(const ADReal & effective_trial_stress,
                                 const ADReal & scalar) override;
  virtual ADReal computeDerivative(const ADReal & effective_trial_stress,
                                   const ADReal & scalar) override;

  std::size_t stressIndex(const ADReal & effective_trial_stress);

  const ADVariableValue * const _temperature;

  const std::vector<Real> _coefficient;

  const std::vector<Real> _n_exponent;

  const std::vector<Real> _m_exponent;

  const std::vector<Real> _stress_thresholds;

  const std::vector<Real> _activation_energy;

  const Real _gas_constant;

  const Real _start_time;

  ADReal _exponential;

  Real _exp_time;

  unsigned int _stress_index;

  const unsigned int _number_of_models;
};