SolidMechanicsPlasticDruckerPrager.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 "SolidMechanicsPlasticModel.h"
#include "SolidMechanicsHardeningModel.h"
 
class SolidMechanicsPlasticDruckerPrager : public SolidMechanicsPlasticModel
{
public:
  static InputParameters validParams();
 
  SolidMechanicsPlasticDruckerPrager(const InputParameters & parameters);
 
  virtual std::string modelName() const override;
 
  void bothAB(Real intnl, Real & aaa, Real & bbb) const;
 
  void dbothAB(Real intnl, Real & daaa, Real & dbbb) const;
 
  enum FrictionDilation
  {
    friction = 0,
    dilation = 1
  };
 
  void onlyB(Real intnl, int fd, Real & bbb) const;
 
  void donlyB(Real intnl, int fd, Real & dbbb) const;
 
protected:
  const SolidMechanicsHardeningModel & _mc_cohesion;
 
  const SolidMechanicsHardeningModel & _mc_phi;
 
  const SolidMechanicsHardeningModel & _mc_psi;
 
  virtual Real yieldFunction(const RankTwoTensor & stress, Real intnl) const override;
 
  virtual RankTwoTensor dyieldFunction_dstress(const RankTwoTensor & stress,
                                               Real intnl) const override;
 
  virtual Real dyieldFunction_dintnl(const RankTwoTensor & stress, Real intnl) const override;
 
  virtual RankTwoTensor flowPotential(const RankTwoTensor & stress, Real intnl) const override;
 
  virtual RankFourTensor dflowPotential_dstress(const RankTwoTensor & stress,
                                                Real intnl) const override;
 
  virtual RankTwoTensor dflowPotential_dintnl(const RankTwoTensor & stress,
                                              Real intnl) const override;
 
  const MooseEnum _mc_interpolation_scheme;
 
  const bool _zero_cohesion_hardening;
 
  const bool _zero_phi_hardening;
 
  const bool _zero_psi_hardening;
 
  virtual RankTwoTensor df_dsig(const RankTwoTensor & stress, Real bbb) const;
 
private:
  Real _aaa;
  Real _bbb;
  Real _bbb_flow;
 
  void initializeAandB(Real intnl, Real & aaa, Real & bbb) const;
 
  void initializeB(Real intnl, int fd, Real & bbb) const;
};