SolidMechanicsPlasticModel.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 "GeneralUserObject.h"
#include "RankTwoTensor.h"

class SolidMechanicsPlasticModel : public GeneralUserObject
{
public:
  static InputParameters validParams();

  SolidMechanicsPlasticModel(const InputParameters & parameters);

  void initialize();
  void execute();
  void finalize();

  virtual unsigned int numberSurfaces() const;

  virtual void
  yieldFunctionV(const RankTwoTensor & stress, Real intnl, std::vector<Real> & f) const;

  virtual void dyieldFunction_dstressV(const RankTwoTensor & stress,
                                       Real intnl,
                                       std::vector<RankTwoTensor> & df_dstress) const;

  virtual void dyieldFunction_dintnlV(const RankTwoTensor & stress,
                                      Real intnl,
                                      std::vector<Real> & df_dintnl) const;

  virtual void
  flowPotentialV(const RankTwoTensor & stress, Real intnl, std::vector<RankTwoTensor> & r) const;

  virtual void dflowPotential_dstressV(const RankTwoTensor & stress,
                                       Real intnl,
                                       std::vector<RankFourTensor> & dr_dstress) const;

  virtual void dflowPotential_dintnlV(const RankTwoTensor & stress,
                                      Real intnl,
                                      std::vector<RankTwoTensor> & dr_dintnl) const;

  virtual void
  hardPotentialV(const RankTwoTensor & stress, Real intnl, std::vector<Real> & h) const;

  virtual void dhardPotential_dstressV(const RankTwoTensor & stress,
                                       Real intnl,
                                       std::vector<RankTwoTensor> & dh_dstress) const;

  virtual void dhardPotential_dintnlV(const RankTwoTensor & stress,
                                      Real intnl,
                                      std::vector<Real> & dh_dintnl) const;

  virtual void activeConstraints(const std::vector<Real> & f,
                                 const RankTwoTensor & stress,
                                 Real intnl,
                                 const RankFourTensor & Eijkl,
                                 std::vector<bool> & act,
                                 RankTwoTensor & returned_stress) const;

  virtual std::string modelName() const = 0;

  virtual bool useCustomReturnMap() const;

  virtual bool useCustomCTO() const;

  const Real _f_tol;

  const Real _ic_tol;

  virtual bool returnMap(const RankTwoTensor & trial_stress,
                         Real intnl_old,
                         const RankFourTensor & E_ijkl,
                         Real ep_plastic_tolerance,
                         RankTwoTensor & returned_stress,
                         Real & returned_intnl,
                         std::vector<Real> & dpm,
                         RankTwoTensor & delta_dp,
                         std::vector<Real> & yf,
                         bool & trial_stress_inadmissible) const;

  virtual RankFourTensor consistentTangentOperator(const RankTwoTensor & trial_stress,
                                                   Real intnl_old,
                                                   const RankTwoTensor & stress,
                                                   Real intnl,
                                                   const RankFourTensor & E_ijkl,
                                                   const std::vector<Real> & cumulative_pm) const;

  bool KuhnTuckerSingleSurface(Real yf, Real dpm, Real dpm_tol) const;

protected:

  virtual Real yieldFunction(const RankTwoTensor & stress, Real intnl) const;

  virtual RankTwoTensor dyieldFunction_dstress(const RankTwoTensor & stress, Real intnl) const;

  virtual Real dyieldFunction_dintnl(const RankTwoTensor & stress, Real intnl) const;

  virtual RankTwoTensor flowPotential(const RankTwoTensor & stress, Real intnl) const;

  virtual RankFourTensor dflowPotential_dstress(const RankTwoTensor & stress, Real intnl) const;

  virtual RankTwoTensor dflowPotential_dintnl(const RankTwoTensor & stress, Real intnl) const;

  virtual Real hardPotential(const RankTwoTensor & stress, Real intnl) const;

  virtual RankTwoTensor dhardPotential_dstress(const RankTwoTensor & stress, Real intnl) const;

  virtual Real dhardPotential_dintnl(const RankTwoTensor & stress, Real intnl) const;
};