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