FiniteStrainPlasticMaterial.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

// Original class author: A.M. Jokisaari,  O. Heinonen

#include "ComputeStressBase.h"

class FiniteStrainPlasticMaterial : public ComputeStressBase
{
public:
  static InputParameters validParams();

  FiniteStrainPlasticMaterial(const InputParameters & parameters);

protected:
  virtual void computeQpStress();
  virtual void initQpStatefulProperties();
  std::vector<Real> _yield_stress_vector;
  MaterialProperty<RankTwoTensor> & _plastic_strain;
  const MaterialProperty<RankTwoTensor> & _plastic_strain_old;
  MaterialProperty<Real> & _eqv_plastic_strain;
  const MaterialProperty<Real> & _eqv_plastic_strain_old;
  const MaterialProperty<RankTwoTensor> & _stress_old;
  const MaterialProperty<RankTwoTensor> & _strain_increment;
  const MaterialProperty<RankTwoTensor> & _rotation_increment;
  const std::string _elasticity_tensor_name;
  const MaterialProperty<RankFourTensor> & _elasticity_tensor;
  Real _rtol;
  Real _ftol;
  Real _eptol;

  // outer and mixed product of the delta function tensors
  RankFourTensor _deltaOuter, _deltaMixed;

  virtual void returnMap(const RankTwoTensor & sig_old,
                         const Real eqvpstrain_old,
                         const RankTwoTensor & plastic_strain_old,
                         const RankTwoTensor & delta_d,
                         const RankFourTensor & E_ijkl,
                         RankTwoTensor & sig,
                         Real & eqvpstrain,
                         RankTwoTensor & plastic_strain);

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

  virtual RankTwoTensor dyieldFunction_dstress(const RankTwoTensor & stress);

  virtual Real dyieldFunction_dinternal(const Real equivalent_plastic_strain);

  virtual RankTwoTensor flowPotential(const RankTwoTensor & stress);

  virtual Real internalPotential();

  Real getSigEqv(const RankTwoTensor & stress);

  virtual void getJac(const RankTwoTensor & sig,
                      const RankFourTensor & E_ijkl,
                      Real flow_incr,
                      RankFourTensor & dresid_dsig);

  Real getYieldStress(const Real equivalent_plastic_strain);

  Real getdYieldStressdPlasticStrain(const Real equivalent_plastic_strain);

  usingTensorIndices(i_, j_, k_, l_);
};