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;
};