ADViscoplasticityStressUpdateBase.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 usingViscoplasticityStressUpdateBaseMembers                                                \
  usingStressUpdateBaseMembers;                                                                    \
  usingSingleVariableReturnMappingSolutionMembers;                                                 \
  using ADViscoplasticityStressUpdateBase<compute_stage>::updateIntermediatePorosity;              \
  using ADViscoplasticityStressUpdateBase<compute_stage>::computeStressInitialize;                 \
  using ADViscoplasticityStressUpdateBase<compute_stage>::computeStressFinalize;                   \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_intermediate_porosity;                  \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_derivative;                             \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_effective_inelastic_strain;             \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_effective_inelastic_strain_old;         \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_inelastic_strain;                       \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_inelastic_strain_old;                   \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_verbose;                                \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_porosity_old;                           \
  using ADViscoplasticityStressUpdateBase<compute_stage>::_max_inelastic_increment
 
template <ComputeStage>
class ADViscoplasticityStressUpdateBase;
 
declareADValidParams(ADViscoplasticityStressUpdateBase);
 
template <ComputeStage compute_stage>
class ADViscoplasticityStressUpdateBase
  : public ADStressUpdateBase<compute_stage>,
    public ADSingleVariableReturnMappingSolution<compute_stage>
{
public:
  static InputParameters validParams();
 
  ADViscoplasticityStressUpdateBase(const InputParameters & parameters);
 
  virtual Real computeReferenceResidual(const ADReal & effective_trial_stress,
                                        const ADReal & scalar_effective_inelastic_strain) override;
 
  virtual Real computeTimeStepLimit() override;
 
  bool requiresIsotropicTensor() override { return true; }
 
protected:
  virtual void initQpStatefulProperties() override;
 
  virtual void propagateQpStatefulProperties() override;
 
  virtual void computeStressInitialize(const ADReal & /*effective_trial_stress*/,
                                       const ADRankFourTensor & /*elasticity_tensor*/)
  {
  }
 
  virtual void computeStressFinalize(const ADRankTwoTensor & /*plastic_strain_increment*/) {}
 
  virtual ADReal computeDerivative(const ADReal & /*effective_trial_stress*/,
                                   const ADReal & /*scalar*/) override
  {
    return _derivative;
  }
 
  void updateIntermediatePorosity(const ADRankTwoTensor & elastic_strain_increment);
 
  void outputIterationSummary(std::stringstream * iter_output,
                              const unsigned int total_it) override;
 
  const std::string _total_strain_base_name;
 
  const ADMaterialProperty(RankTwoTensor) & _strain_increment;
 
  ADMaterialProperty(Real) & _effective_inelastic_strain;
  const MaterialProperty<Real> & _effective_inelastic_strain_old;
 
  ADMaterialProperty(RankTwoTensor) & _inelastic_strain;
  const MaterialProperty<RankTwoTensor> & _inelastic_strain_old;
 
  Real _max_inelastic_increment;
 
  ADReal _intermediate_porosity;
 
  const MaterialProperty<Real> & _porosity_old;
 
  const bool _verbose;
 
  ADReal _derivative;
 
  usingStressUpdateBaseMembers;
  usingSingleVariableReturnMappingSolutionMembers;
};