RadialReturnStressUpdate.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 "StressUpdateBase.h"
#include "SingleVariableReturnMappingSolution.h"
 
// Forward declaration
class RadialReturnStressUpdate;
 
template <>
InputParameters validParams<RadialReturnStressUpdate>();
 
class RadialReturnStressUpdate : public StressUpdateBase, public SingleVariableReturnMappingSolution
{
public:
  static InputParameters validParams();
 
  RadialReturnStressUpdate(const InputParameters & parameters);
 
  virtual void updateState(RankTwoTensor & strain_increment,
                           RankTwoTensor & inelastic_strain_increment,
                           const RankTwoTensor & rotation_increment,
                           RankTwoTensor & stress_new,
                           const RankTwoTensor & stress_old,
                           const RankFourTensor & elasticity_tensor,
                           const RankTwoTensor & elastic_strain_old,
                           bool compute_full_tangent_operator,
                           RankFourTensor & tangent_operator) override;
 
  virtual Real computeReferenceResidual(const Real effective_trial_stress,
                                        const Real scalar_effective_inelastic_strain) override;
 
  virtual Real minimumPermissibleValue(const Real /*effective_trial_stress*/) const override
  {
    return 0.0;
  }
 
  virtual Real maximumPermissibleValue(const Real effective_trial_stress) const override;
 
  virtual Real computeTimeStepLimit() override;
 
  bool requiresIsotropicTensor() override { return true; }
 
  bool isIsotropic() override { return true; };
 
protected:
  virtual void initQpStatefulProperties() override;
 
  void propagateQpStatefulPropertiesRadialReturn();
 
  virtual void computeStressInitialize(const Real /*effective_trial_stress*/,
                                       const RankFourTensor & /*elasticity_tensor*/)
  {
  }
 
  virtual Real computeStressDerivative(const Real /*effective_trial_stress*/, const Real /*scalar*/)
  {
    return 0.0;
  }
 
  virtual void computeStressFinalize(const RankTwoTensor & /*inelasticStrainIncrement*/) {}
 
  void outputIterationSummary(std::stringstream * iter_output,
                              const unsigned int total_it) override;
 
  Real _three_shear_modulus;
 
  MaterialProperty<Real> & _effective_inelastic_strain;
  const MaterialProperty<Real> & _effective_inelastic_strain_old;
  Real _max_inelastic_increment;
 
  const RankTwoTensor _identity_two;
 
  const RankFourTensor _identity_symmetric_four;
 
  const RankFourTensor _deviatoric_projection_four;
};