CappedMohrCoulombStressUpdate.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 "MultiParameterPlasticityStressUpdate.h"
#include "SolidMechanicsHardeningModel.h"

class CappedMohrCoulombStressUpdate : public MultiParameterPlasticityStressUpdate
{
public:
  static InputParameters validParams();

  CappedMohrCoulombStressUpdate(const InputParameters & parameters);

  bool requiresIsotropicTensor() override { return true; }

  bool isIsotropic() override { return true; };

protected:
  const SolidMechanicsHardeningModel & _tensile_strength;

  const SolidMechanicsHardeningModel & _compressive_strength;

  const SolidMechanicsHardeningModel & _cohesion;

  const SolidMechanicsHardeningModel & _phi;

  const SolidMechanicsHardeningModel & _psi;

  const bool _perfect_guess;

  Real _poissons_ratio;

  const Real _shifter;

  RankTwoTensor _eigvecs;

  void computeStressParams(const RankTwoTensor & stress,
                           std::vector<Real> & stress_params) const override;

  void dstressparam_dstress(const RankTwoTensor & stress,
                            std::vector<RankTwoTensor> & dsp) const override;

  void d2stressparam_dstress(const RankTwoTensor & stress,
                             std::vector<RankFourTensor> & d2sp) const override;

  virtual void setStressAfterReturnV(const RankTwoTensor & stress_trial,
                                     const std::vector<Real> & stress_params,
                                     Real gaE,
                                     const std::vector<Real> & intnl,
                                     const yieldAndFlow & smoothed_q,
                                     const RankFourTensor & Eijkl,
                                     RankTwoTensor & stress) const override;

  virtual void preReturnMapV(const std::vector<Real> & trial_stress_params,
                             const RankTwoTensor & stress_trial,
                             const std::vector<Real> & intnl_old,
                             const std::vector<Real> & yf,
                             const RankFourTensor & Eijkl) override;

  void setEffectiveElasticity(const RankFourTensor & Eijkl) override;

  void yieldFunctionValuesV(const std::vector<Real> & stress_params,
                            const std::vector<Real> & intnl,
                            std::vector<Real> & yf) const override;

  void computeAllQV(const std::vector<Real> & stress_params,
                    const std::vector<Real> & intnl,
                    std::vector<yieldAndFlow> & all_q) const override;

  void initializeVarsV(const std::vector<Real> & trial_stress_params,
                       const std::vector<Real> & intnl_old,
                       std::vector<Real> & stress_params,
                       Real & gaE,
                       std::vector<Real> & intnl) const override;

  void setIntnlValuesV(const std::vector<Real> & trial_stress_params,
                       const std::vector<Real> & current_stress_params,
                       const std::vector<Real> & intnl_old,
                       std::vector<Real> & intnl) const override;

  void setIntnlDerivativesV(const std::vector<Real> & trial_stress_params,
                            const std::vector<Real> & current_stress_params,
                            const std::vector<Real> & intnl,
                            std::vector<std::vector<Real>> & dintnl) const override;

  virtual void consistentTangentOperatorV(const RankTwoTensor & stress_trial,
                                          const std::vector<Real> & trial_stress_params,
                                          const RankTwoTensor & stress,
                                          const std::vector<Real> & stress_params,
                                          Real gaE,
                                          const yieldAndFlow & smoothed_q,
                                          const RankFourTensor & Eijkl,
                                          bool compute_full_tangent_operator,
                                          const std::vector<std::vector<Real>> & dvar_dtrial,
                                          RankFourTensor & cto) override;

private:
  mutable std::vector<RankTwoTensor> _dsp_trial_scratch;

  mutable std::vector<Real> _eigvals_scratch;

  mutable std::vector<Real> _dga_shear_scratch;

  mutable std::vector<Real> _dshear_correction_scratch;
};