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