CappedMohrCoulombStressUpdate.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 "MultiParameterPlasticityStressUpdate.h"
#include "TensorMechanicsHardeningModel.h"
 
class CappedMohrCoulombStressUpdate;
 
template <>
InputParameters validParams<CappedMohrCoulombStressUpdate>();
 
class CappedMohrCoulombStressUpdate : public MultiParameterPlasticityStressUpdate
{
public:
  static InputParameters validParams();
 
  CappedMohrCoulombStressUpdate(const InputParameters & parameters);
 
  bool requiresIsotropicTensor() override { return true; }
 
  bool isIsotropic() override { return true; };
 
protected:
  const TensorMechanicsHardeningModel & _tensile_strength;
 
  const TensorMechanicsHardeningModel & _compressive_strength;
 
  const TensorMechanicsHardeningModel & _cohesion;
 
  const TensorMechanicsHardeningModel & _phi;
 
  const TensorMechanicsHardeningModel & _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;
 
  std::vector<RankTwoTensor> dstress_param_dstress(const RankTwoTensor & stress) const override;
 
  std::vector<RankFourTensor> d2stress_param_dstress(const RankTwoTensor & stress) 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;
};