ADComputeMultipleInelasticStress.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 "ADComputeFiniteStrainElasticStress.h"
#include "ADRankTwoTensorForward.h"
#include "ADRankFourTensorForward.h"
#include "StressUpdateBase.h"
#include "DamageBase.h"

class ADComputeMultipleInelasticStress : public ADComputeFiniteStrainElasticStress
{
public:
  static InputParameters validParams();

  ADComputeMultipleInelasticStress(const InputParameters & parameters);

  virtual void initialSetup() override;

protected:
  virtual void initQpStatefulProperties() override;

  virtual void computeQpStress() override;

  virtual void computeQpStressIntermediateConfiguration();

  virtual void finiteStrainRotation();

  virtual void updateQpState(ADRankTwoTensor & elastic_strain_increment,
                             ADRankTwoTensor & combined_inelastic_strain_increment);

  virtual void updateQpStateSingleModel(unsigned model_number,
                                        ADRankTwoTensor & elastic_strain_increment,
                                        ADRankTwoTensor & combined_inelastic_strain_increment);

  virtual void computeAdmissibleState(unsigned model_number,
                                      ADRankTwoTensor & elastic_strain_increment,
                                      ADRankTwoTensor & inelastic_strain_increment);

  const unsigned int _max_iterations;
  const Real _relative_tolerance;
  const Real _absolute_tolerance;
  const bool _internal_solve_full_iteration_history;

  const bool _perform_finite_strain_rotations;

  ADMaterialProperty<RankTwoTensor> & _inelastic_strain;

  const MaterialProperty<RankTwoTensor> & _inelastic_strain_old;

  const unsigned _num_models;

  const std::vector<Real> _inelastic_weights;

  const bool _cycle_models;

  MaterialProperty<Real> & _material_timestep_limit;

  std::vector<ADStressUpdateBase *> _models;

  bool _is_elasticity_tensor_guaranteed_isotropic;

  DamageBaseTempl<true> * _damage_model;

  RankTwoTensor _undamaged_stress_old;
};