www.mooseframework.org
FiniteStrainPlasticMaterial.h
Go to the documentation of this file.
1 //* This file is part of the MOOSE framework
2 //* https://www.mooseframework.org
3 //*
4 //* All rights reserved, see COPYRIGHT for full restrictions
5 //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
6 //*
7 //* Licensed under LGPL 2.1, please see LICENSE for details
8 //* https://www.gnu.org/licenses/lgpl-2.1.html
9 
10 #pragma once
11 // Original class author: A.M. Jokisaari, O. Heinonen
12 
13 #include "ComputeStressBase.h"
14 
24 {
25 public:
27 
29 
30 protected:
31  virtual void computeQpStress();
32  virtual void initQpStatefulProperties();
33  std::vector<Real> _yield_stress_vector;
43  const std::string _elasticity_tensor_name;
49 
50  // outer and mixed product of the delta function tensors
52 
70  virtual void returnMap(const RankTwoTensor & sig_old,
71  const Real eqvpstrain_old,
72  const RankTwoTensor & plastic_strain_old,
73  const RankTwoTensor & delta_d,
74  const RankFourTensor & E_ijkl,
75  RankTwoTensor & sig,
76  Real & eqvpstrain,
77  RankTwoTensor & plastic_strain);
78 
85  virtual Real yieldFunction(const RankTwoTensor & stress, const Real yield_stress);
86 
90  virtual RankTwoTensor dyieldFunction_dstress(const RankTwoTensor & stress);
91 
95  virtual Real dyieldFunction_dinternal(const Real equivalent_plastic_strain);
96 
102  virtual RankTwoTensor flowPotential(const RankTwoTensor & stress);
103 
107  virtual Real internalPotential();
108 
112  Real getSigEqv(const RankTwoTensor & stress);
113 
122  virtual void getJac(const RankTwoTensor & sig,
123  const RankFourTensor & E_ijkl,
124  Real flow_incr,
125  RankFourTensor & dresid_dsig);
126 
131  Real getYieldStress(const Real equivalent_plastic_strain);
132 
136  Real getdYieldStressdPlasticStrain(const Real equivalent_plastic_strain);
137 
139  usingTensorIndices(i_, j_, k_, l_);
140 };
virtual void computeQpStress()
Compute the stress and store it in the _stress material property for the current quadrature point...
FiniteStrainPlasticMaterial implements rate-independent associative J2 plasticity with isotropic hard...
ComputeStressBase is the base class for stress tensors computed from MOOSE&#39;s strain calculators...
const MaterialProperty< RankTwoTensor > & _strain_increment
MaterialProperty< Real > & _eqv_plastic_strain
virtual Real dyieldFunction_dinternal(const Real equivalent_plastic_strain)
Derivative of yieldFunction with respect to the equivalent plastic strain.
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)
Implements the return map.
virtual RankTwoTensor flowPotential(const RankTwoTensor &stress)
Flow potential, which in this case is just dyieldFunction_dstress because we are doing associative fl...
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Elasticity tensor material property.
const MaterialProperty< Real > & _eqv_plastic_strain_old
virtual Real yieldFunction(const RankTwoTensor &stress, const Real yield_stress)
Calculates the yield function.
Real getYieldStress(const Real equivalent_plastic_strain)
yield stress as a function of equivalent plastic strain.
const std::string _elasticity_tensor_name
Name of the elasticity tensor material property.
FiniteStrainPlasticMaterial(const InputParameters &parameters)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const MaterialProperty< RankTwoTensor > & _plastic_strain_old
Real getSigEqv(const RankTwoTensor &stress)
Equivalent stress.
const InputParameters & parameters() const
virtual RankTwoTensor dyieldFunction_dstress(const RankTwoTensor &stress)
Derivative of yieldFunction with respect to the stress.
virtual Real internalPotential()
The internal potential.
const MaterialProperty< RankTwoTensor > & _rotation_increment
Real getdYieldStressdPlasticStrain(const Real equivalent_plastic_strain)
d(yieldstress)/d(equivalent plastic strain)
MaterialProperty< RankTwoTensor > & _plastic_strain
virtual void getJac(const RankTwoTensor &sig, const RankFourTensor &E_ijkl, Real flow_incr, RankFourTensor &dresid_dsig)
Evaluates the derivative d(resid_ij)/d(sig_kl), where resid_ij = flow_incr*flowPotential_ij - (E^{-1}...
usingTensorIndices(i_, j_, k_, l_)
make i,j,k,l available as tensor indices for mixedProduct
static InputParameters validParams()
const MaterialProperty< RankTwoTensor > & _stress_old
The old stress tensor.