ADRadialReturnStressUpdate.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "ADStressUpdateBase.h"
#include "ADSingleVariableReturnMappingSolution.h"
 
#define usingRadialReturnStressUpdateMembers                                                       \
  usingStressUpdateBaseMembers;                                                                    \
  usingSingleVariableReturnMappingSolutionMembers;                                                 \
  using ADRadialReturnStressUpdate<compute_stage>::propagateQpStatefulPropertiesRadialReturn;      \
  using ADRadialReturnStressUpdate<compute_stage>::_three_shear_modulus;                           \
  using ADRadialReturnStressUpdate<compute_stage>::_effective_inelastic_strain;                    \
  using ADRadialReturnStressUpdate<compute_stage>::_effective_inelastic_strain_old;                \
  using ADRadialReturnStressUpdate<compute_stage>::_max_inelastic_increment
 
// Forward declarations
template <ComputeStage>
class ADRadialReturnStressUpdate;
 
declareADValidParams(ADRadialReturnStressUpdate);
 
template <ComputeStage compute_stage>
class ADRadialReturnStressUpdate : public ADStressUpdateBase<compute_stage>,
                                   public ADSingleVariableReturnMappingSolution<compute_stage>
{
public:
  static InputParameters validParams();
 
  ADRadialReturnStressUpdate(const InputParameters & parameters);
 
  virtual void updateState(ADRankTwoTensor & strain_increment,
                           ADRankTwoTensor & inelastic_strain_increment,
                           const ADRankTwoTensor & rotation_increment,
                           ADRankTwoTensor & stress_new,
                           const RankTwoTensor & stress_old,
                           const ADRankFourTensor & elasticity_tensor,
                           const RankTwoTensor & elastic_strain_old) override;
 
  virtual Real computeReferenceResidual(const ADReal & effective_trial_stress,
                                        const ADReal & scalar_effective_inelastic_strain) override;
 
  virtual ADReal minimumPermissibleValue(const ADReal & /*effective_trial_stress*/) const override
  {
    return 0.0;
  }
 
  virtual ADReal maximumPermissibleValue(const ADReal & effective_trial_stress) const override;
 
  virtual Real computeTimeStepLimit() override;
 
  bool requiresIsotropicTensor() override { return true; }
 
protected:
  virtual void initQpStatefulProperties() override;
 
  void propagateQpStatefulPropertiesRadialReturn();
 
  virtual void computeStressInitialize(const ADReal & /*effective_trial_stress*/,
                                       const ADRankFourTensor & /*elasticity_tensor*/)
  {
  }
 
  virtual void computeStressFinalize(const ADRankTwoTensor & /*inelasticStrainIncrement*/) {}
 
  void outputIterationSummary(std::stringstream * iter_output,
                              const unsigned int total_it) override;
 
  ADReal _three_shear_modulus;
 
  ADMaterialProperty(Real) & _effective_inelastic_strain;
  const MaterialProperty<Real> & _effective_inelastic_strain_old;
  Real _max_inelastic_increment;
  const bool _apply_strain;
 
  usingStressUpdateBaseMembers;
  usingSingleVariableReturnMappingSolutionMembers;
};