SolidMechanicsPlasticMohrCoulomb.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 SolidMechanicsPlasticMohrCoulomb : public SolidMechanicsPlasticModel
{
public:
  static InputParameters validParams();
 
  SolidMechanicsPlasticMohrCoulomb(const InputParameters & parameters);
 
  virtual std::string modelName() const override;
 
protected:
  Real yieldFunction(const RankTwoTensor & stress, Real intnl) const override;
 
  RankTwoTensor dyieldFunction_dstress(const RankTwoTensor & stress, Real intnl) const override;
 
  Real dyieldFunction_dintnl(const RankTwoTensor & stress, Real intnl) const override;
 
  RankTwoTensor flowPotential(const RankTwoTensor & stress, Real intnl) const override;
 
  RankFourTensor dflowPotential_dstress(const RankTwoTensor & stress, Real intnl) const override;
 
  RankTwoTensor dflowPotential_dintnl(const RankTwoTensor & stress, Real intnl) const override;
 
  const SolidMechanicsHardeningModel & _cohesion;
 
  const SolidMechanicsHardeningModel & _phi;
 
  const SolidMechanicsHardeningModel & _psi;
 
  MooseEnum _tip_scheme;
 
  Real _small_smoother2;
 
  Real _cap_start;
 
  Real _cap_rate;
 
  Real _tt;
 
  Real _costt;
 
  Real _sintt;
 
  Real _cos3tt;
 
  Real _sin3tt;
 
  Real _cos6tt;
 
  Real _sin6tt;
 
  Real _lode_cutoff;
 
  virtual Real smooth(const RankTwoTensor & stress) const;
 
  virtual Real dsmooth(const RankTwoTensor & stress) const;
 
  virtual Real d2smooth(const RankTwoTensor & stress) const;
 
  virtual Real cohesion(const Real internal_param) const;
 
  virtual Real dcohesion(const Real internal_param) const;
 
  virtual Real phi(const Real internal_param) const;
 
  virtual Real dphi(const Real internal_param) const;
 
  virtual Real psi(const Real internal_param) const;
 
  virtual Real dpsi(const Real internal_param) const;
 
private:
  void abbo(const Real sin3lode, const Real sin_angle, Real & aaa, Real & bbb, Real & ccc) const;
 
  void
  dabbo(const Real sin3lode, const Real sin_angle, Real & daaa, Real & dbbb, Real & dccc) const;
 
  RankTwoTensor df_dsig(const RankTwoTensor & stress, const Real sin_angle) const;
};