www.mooseframework.org
Public Member Functions | Public Attributes | Protected Member Functions | Protected Attributes | List of all members
TensorMechanicsPlasticWeakPlaneShear Class Reference

Rate-independent associative weak-plane tensile failure with hardening/softening. More...

#include <TensorMechanicsPlasticWeakPlaneShear.h>

Inheritance diagram for TensorMechanicsPlasticWeakPlaneShear:
[legend]

Public Member Functions

 TensorMechanicsPlasticWeakPlaneShear (const InputParameters &parameters)
 
virtual void activeConstraints (const std::vector< Real > &f, const RankTwoTensor &stress, Real intnl, const RankFourTensor &Eijkl, std::vector< bool > &act, RankTwoTensor &returned_stress) const override
 The active yield surfaces, given a vector of yield functions. More...
 
virtual std::string modelName () const override
 
void initialize ()
 
void execute ()
 
void finalize ()
 
virtual unsigned int numberSurfaces () const
 The number of yield surfaces for this plasticity model. More...
 
virtual void yieldFunctionV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &f) const
 Calculates the yield functions. More...
 
virtual void dyieldFunction_dstressV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &df_dstress) const
 The derivative of yield functions with respect to stress. More...
 
virtual void dyieldFunction_dintnlV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &df_dintnl) const
 The derivative of yield functions with respect to the internal parameter. More...
 
virtual void flowPotentialV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &r) const
 The flow potentials. More...
 
virtual void dflowPotential_dstressV (const RankTwoTensor &stress, Real intnl, std::vector< RankFourTensor > &dr_dstress) const
 The derivative of the flow potential with respect to stress. More...
 
virtual void dflowPotential_dintnlV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &dr_dintnl) const
 The derivative of the flow potential with respect to the internal parameter. More...
 
virtual void hardPotentialV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &h) const
 The hardening potential. More...
 
virtual void dhardPotential_dstressV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &dh_dstress) const
 The derivative of the hardening potential with respect to stress. More...
 
virtual void dhardPotential_dintnlV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &dh_dintnl) const
 The derivative of the hardening potential with respect to the internal parameter. More...
 
virtual bool useCustomReturnMap () const
 Returns false. You will want to override this in your derived class if you write a custom returnMap function. More...
 
virtual bool useCustomCTO () const
 Returns false. You will want to override this in your derived class if you write a custom consistent tangent operator function. More...
 
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
 Performs a custom return-map. More...
 
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
 Calculates a custom consistent tangent operator. More...
 
bool KuhnTuckerSingleSurface (Real yf, Real dpm, Real dpm_tol) const
 Returns true if the Kuhn-Tucker conditions for the single surface are satisfied. More...
 

Public Attributes

const Real _f_tol
 Tolerance on yield function. More...
 
const Real _ic_tol
 Tolerance on internal constraint. More...
 

Protected Member Functions

Real yieldFunction (const RankTwoTensor &stress, Real intnl) const override
 The following functions are what you should override when building single-plasticity models. More...
 
RankTwoTensor dyieldFunction_dstress (const RankTwoTensor &stress, Real intnl) const override
 The derivative of yield function with respect to stress. More...
 
Real dyieldFunction_dintnl (const RankTwoTensor &stress, Real intnl) const override
 The derivative of yield function with respect to the internal parameter. More...
 
RankTwoTensor flowPotential (const RankTwoTensor &stress, Real intnl) const override
 The flow potential. More...
 
RankFourTensor dflowPotential_dstress (const RankTwoTensor &stress, Real intnl) const override
 The derivative of the flow potential with respect to stress. More...
 
RankTwoTensor dflowPotential_dintnl (const RankTwoTensor &stress, Real intnl) const override
 The derivative of the flow potential with respect to the internal parameter. More...
 
RankTwoTensor df_dsig (const RankTwoTensor &stress, Real _tan_phi_or_psi) const
 Function that's used in dyieldFunction_dstress and flowPotential. More...
 
virtual Real smooth (const RankTwoTensor &stress) const
 returns the 'a' parameter - see doco for _tip_scheme More...
 
virtual Real dsmooth (const RankTwoTensor &stress) const
 returns the da/dstress(2,2) - see doco for _tip_scheme More...
 
virtual Real d2smooth (const RankTwoTensor &stress) const
 returns the d^2a/dstress(2,2)^2 - see doco for _tip_scheme More...
 
virtual Real cohesion (const Real internal_param) const
 cohesion as a function of internal parameter More...
 
virtual Real dcohesion (const Real internal_param) const
 d(cohesion)/d(internal_param) More...
 
virtual Real tan_phi (const Real internal_param) const
 tan_phi as a function of internal parameter More...
 
virtual Real dtan_phi (const Real internal_param) const
 d(tan_phi)/d(internal_param); More...
 
virtual Real tan_psi (const Real internal_param) const
 tan_psi as a function of internal parameter More...
 
virtual Real dtan_psi (const Real internal_param) const
 d(tan_psi)/d(internal_param); More...
 
virtual Real hardPotential (const RankTwoTensor &stress, Real intnl) const
 The hardening potential. More...
 
virtual RankTwoTensor dhardPotential_dstress (const RankTwoTensor &stress, Real intnl) const
 The derivative of the hardening potential with respect to stress. More...
 
virtual Real dhardPotential_dintnl (const RankTwoTensor &stress, Real intnl) const
 The derivative of the hardening potential with respect to the internal parameter. More...
 

Protected Attributes

const TensorMechanicsHardeningModel_cohesion
 Hardening model for cohesion. More...
 
const TensorMechanicsHardeningModel_tan_phi
 Hardening model for tan(phi) More...
 
const TensorMechanicsHardeningModel_tan_psi
 Hardening model for tan(psi) More...
 
MooseEnum _tip_scheme
 The yield function is modified to f = sqrt(s_xz^2 + s_yz^2 + a) + s_zz*_tan_phi - _cohesion where "a" depends on the tip_scheme. More...
 
Real _small_smoother2
 smoothing parameter for the cone's tip - see doco for _tip_scheme More...
 
Real _cap_start
 smoothing parameter dictating when the 'cap' will start - see doco for _tip_scheme More...
 
Real _cap_rate
 dictates how quickly the 'cap' degenerates to a hemisphere - see doco for _tip_scheme More...
 

Detailed Description

Rate-independent associative weak-plane tensile failure with hardening/softening.

The cone's tip is smoothed.

Definition at line 25 of file TensorMechanicsPlasticWeakPlaneShear.h.

Constructor & Destructor Documentation

◆ TensorMechanicsPlasticWeakPlaneShear()

TensorMechanicsPlasticWeakPlaneShear::TensorMechanicsPlasticWeakPlaneShear ( const InputParameters &  parameters)

Definition at line 61 of file TensorMechanicsPlasticWeakPlaneShear.C.

63  : TensorMechanicsPlasticModel(parameters),
64  _cohesion(getUserObject<TensorMechanicsHardeningModel>("cohesion")),
65  _tan_phi(getUserObject<TensorMechanicsHardeningModel>("tan_friction_angle")),
66  _tan_psi(getUserObject<TensorMechanicsHardeningModel>("tan_dilation_angle")),
67  _tip_scheme(getParam<MooseEnum>("tip_scheme")),
68  _small_smoother2(Utility::pow<2>(getParam<Real>("smoother"))),
69  _cap_start(getParam<Real>("cap_start")),
70  _cap_rate(getParam<Real>("cap_rate"))
71 {
72  // With arbitary UserObjects, it is impossible to check everything, and
73  // I think this is the best I can do
74  if (tan_phi(0) < 0 || tan_psi(0) < 0)
75  mooseError("Weak-Plane-Shear friction and dilation angles must lie in [0, Pi/2]");
76  if (tan_phi(0) < tan_psi(0))
77  mooseError(
78  "Weak-Plane-Shear friction angle must not be less than Weak-Plane-Shear dilation angle");
79  if (cohesion(0) < 0)
80  mooseError("Weak-Plane-Shear cohesion must not be negative");
81 }
MooseEnum _tip_scheme
The yield function is modified to f = sqrt(s_xz^2 + s_yz^2 + a) + s_zz*_tan_phi - _cohesion where "a"...
TensorMechanicsPlasticModel(const InputParameters &parameters)
virtual Real tan_psi(const Real internal_param) const
tan_psi as a function of internal parameter
const TensorMechanicsHardeningModel & _tan_phi
Hardening model for tan(phi)
Real _cap_start
smoothing parameter dictating when the &#39;cap&#39; will start - see doco for _tip_scheme ...
const TensorMechanicsHardeningModel & _cohesion
Hardening model for cohesion.
Real _small_smoother2
smoothing parameter for the cone&#39;s tip - see doco for _tip_scheme
Real _cap_rate
dictates how quickly the &#39;cap&#39; degenerates to a hemisphere - see doco for _tip_scheme ...
virtual Real tan_phi(const Real internal_param) const
tan_phi as a function of internal parameter
virtual Real cohesion(const Real internal_param) const
cohesion as a function of internal parameter
const TensorMechanicsHardeningModel & _tan_psi
Hardening model for tan(psi)

Member Function Documentation

◆ activeConstraints()

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

The active yield surfaces, given a vector of yield functions.

This is used by FiniteStrainMultiPlasticity to determine the initial set of active constraints at the trial (stress, intnl) configuration. It is up to you (the coder) to determine how accurate you want the returned_stress to be. Currently it is only used by FiniteStrainMultiPlasticity to estimate a good starting value for the Newton-Rahson procedure, so currently it may not need to be super perfect.

Parameters
fvalues of the yield functions
stressstress tensor
intnlinternal parameter
Eijklelasticity tensor (stress = Eijkl*strain)
[out]actact[i] = true if the i_th yield function is active
[out]returned_stressApproximate value of the returned stress

Reimplemented from TensorMechanicsPlasticModel.

Definition at line 277 of file TensorMechanicsPlasticWeakPlaneShear.C.

283 {
284  act.assign(1, false);
285  returned_stress = stress;
286 
287  if (f[0] <= _f_tol)
288  return;
289 
290  // in the following i will derive returned_stress for the case smoother=0
291 
292  Real tanpsi = tan_psi(intnl);
293  Real tanphi = tan_phi(intnl);
294 
295  // norm is the normal to the yield surface
296  // with f having psi (dilation angle) instead of phi:
297  // norm(0) = df/dsig(2,0) = df/dsig(0,2)
298  // norm(1) = df/dsig(2,1) = df/dsig(1,2)
299  // norm(2) = df/dsig(2,2)
300  std::vector<Real> norm(3, 0.0);
301  const Real tau = std::sqrt(Utility::pow<2>((stress(0, 2) + stress(2, 0)) / 2.0) +
302  Utility::pow<2>((stress(1, 2) + stress(2, 1)) / 2.0));
303  if (tau > 0.0)
304  {
305  norm[0] = 0.25 * (stress(0, 2) + stress(2, 0)) / tau;
306  norm[1] = 0.25 * (stress(1, 2) + stress(2, 1)) / tau;
307  }
308  else
309  {
310  returned_stress(2, 2) = cohesion(intnl) / tanphi;
311  act[0] = true;
312  return;
313  }
314  norm[2] = tanpsi;
315 
316  // to get the flow directions, we have to multiply norm by Eijkl.
317  // I assume that E(0,2,0,2) = E(1,2,1,2), and E(2,2,0,2) = 0 = E(0,2,1,2), etc
318  // with the usual symmetry. This makes finding the returned_stress
319  // much easier.
320  // returned_stress = stress - alpha*n
321  // where alpha is chosen so that f = 0
322  Real alpha = f[0] / (Eijkl(0, 2, 0, 2) + Eijkl(2, 2, 2, 2) * tanpsi * tanphi);
323 
324  if (1 - alpha * Eijkl(0, 2, 0, 2) / tau >= 0)
325  {
326  // returning to the "surface" of the cone
327  returned_stress(2, 2) = stress(2, 2) - alpha * Eijkl(2, 2, 2, 2) * norm[2];
328  returned_stress(0, 2) = returned_stress(2, 0) =
329  stress(0, 2) - alpha * 2.0 * Eijkl(0, 2, 0, 2) * norm[0];
330  returned_stress(1, 2) = returned_stress(2, 1) =
331  stress(1, 2) - alpha * 2.0 * Eijkl(1, 2, 1, 2) * norm[1];
332  }
333  else
334  {
335  // returning to the "tip" of the cone
336  returned_stress(2, 2) = cohesion(intnl) / tanphi;
337  returned_stress(0, 2) = returned_stress(2, 0) = returned_stress(1, 2) = returned_stress(2, 1) =
338  0.0;
339  }
340  returned_stress(0, 0) =
341  stress(0, 0) - Eijkl(0, 0, 2, 2) * (stress(2, 2) - returned_stress(2, 2)) / Eijkl(2, 2, 2, 2);
342  returned_stress(1, 1) =
343  stress(1, 1) - Eijkl(1, 1, 2, 2) * (stress(2, 2) - returned_stress(2, 2)) / Eijkl(2, 2, 2, 2);
344 
345  act[0] = true;
346 }
virtual Real tan_psi(const Real internal_param) const
tan_psi as a function of internal parameter
virtual Real tan_phi(const Real internal_param) const
tan_phi as a function of internal parameter
const Real _f_tol
Tolerance on yield function.
virtual Real cohesion(const Real internal_param) const
cohesion as a function of internal parameter

◆ cohesion()

Real TensorMechanicsPlasticWeakPlaneShear::cohesion ( const Real  internal_param) const
protectedvirtual

cohesion as a function of internal parameter

Definition at line 184 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by activeConstraints(), TensorMechanicsPlasticWeakPlaneShear(), and yieldFunction().

185 {
186  return _cohesion.value(internal_param);
187 }
const TensorMechanicsHardeningModel & _cohesion
Hardening model for cohesion.
virtual Real value(Real intnl) const

◆ consistentTangentOperator()

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

Calculates a custom consistent tangent operator.

You may choose to over-ride this in your derived TensorMechanicsPlasticXXXX class.

(Note, if you over-ride returnMap, you will probably want to override consistentTangentOpertor too, otherwise it will default to E_ijkl.)

Parameters
stress_oldtrial stress before returning
intnl_oldinternal parameter before returning
stresscurrent returned stress state
intnlinternal parameter
E_ijklelasticity tensor
cumulative_pmthe cumulative plastic multipliers
Returns
the consistent tangent operator: E_ijkl if not over-ridden

Reimplemented in TensorMechanicsPlasticTensileMulti, TensorMechanicsPlasticDruckerPragerHyperbolic, TensorMechanicsPlasticMeanCapTC, and TensorMechanicsPlasticJ2.

Definition at line 253 of file TensorMechanicsPlasticModel.C.

Referenced by TensorMechanicsPlasticJ2::consistentTangentOperator(), TensorMechanicsPlasticDruckerPragerHyperbolic::consistentTangentOperator(), TensorMechanicsPlasticMeanCapTC::consistentTangentOperator(), and TensorMechanicsPlasticTensileMulti::consistentTangentOperator().

260 {
261  return E_ijkl;
262 }

◆ d2smooth()

Real TensorMechanicsPlasticWeakPlaneShear::d2smooth ( const RankTwoTensor &  stress) const
protectedvirtual

returns the d^2a/dstress(2,2)^2 - see doco for _tip_scheme

Definition at line 255 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by dflowPotential_dstress().

256 {
257  Real d2smoother2 = 0;
258  if (_tip_scheme == "cap")
259  {
260  Real x = stress(2, 2) - _cap_start;
261  Real p = 0;
262  Real dp_dx = 0;
263  Real d2p_dx2 = 0;
264  if (x > 0)
265  {
266  const Real exp_cap_rate_x = std::exp(-_cap_rate * x);
267  p = x * (1.0 - exp_cap_rate_x);
268  dp_dx = (1.0 - exp_cap_rate_x) + x * _cap_rate * exp_cap_rate_x;
269  d2p_dx2 = 2.0 * _cap_rate * exp_cap_rate_x - x * Utility::pow<2>(_cap_rate) * exp_cap_rate_x;
270  }
271  d2smoother2 += 2.0 * Utility::pow<2>(dp_dx) + 2.0 * p * d2p_dx2;
272  }
273  return d2smoother2;
274 }
MooseEnum _tip_scheme
The yield function is modified to f = sqrt(s_xz^2 + s_yz^2 + a) + s_zz*_tan_phi - _cohesion where "a"...
Real _cap_start
smoothing parameter dictating when the &#39;cap&#39; will start - see doco for _tip_scheme ...
Real _cap_rate
dictates how quickly the &#39;cap&#39; degenerates to a hemisphere - see doco for _tip_scheme ...

◆ dcohesion()

Real TensorMechanicsPlasticWeakPlaneShear::dcohesion ( const Real  internal_param) const
protectedvirtual

d(cohesion)/d(internal_param)

Definition at line 190 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by dyieldFunction_dintnl().

191 {
192  return _cohesion.derivative(internal_param);
193 }
const TensorMechanicsHardeningModel & _cohesion
Hardening model for cohesion.
virtual Real derivative(Real intnl) const

◆ df_dsig()

RankTwoTensor TensorMechanicsPlasticWeakPlaneShear::df_dsig ( const RankTwoTensor &  stress,
Real  _tan_phi_or_psi 
) const
protected

Function that's used in dyieldFunction_dstress and flowPotential.

Definition at line 93 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by dyieldFunction_dstress(), and flowPotential().

95 {
96  RankTwoTensor deriv; // the constructor zeroes this
97 
98  Real tau = std::sqrt(Utility::pow<2>((stress(0, 2) + stress(2, 0)) / 2.0) +
99  Utility::pow<2>((stress(1, 2) + stress(2, 1)) / 2.0) + smooth(stress));
100  // note that i explicitly symmeterise in preparation for Cosserat
101  if (tau == 0.0)
102  {
103  // the derivative is not defined here, but i want to set it nonzero
104  // because otherwise the return direction might be too crazy
105  deriv(0, 2) = deriv(2, 0) = 0.5;
106  deriv(1, 2) = deriv(2, 1) = 0.5;
107  }
108  else
109  {
110  deriv(0, 2) = deriv(2, 0) = 0.25 * (stress(0, 2) + stress(2, 0)) / tau;
111  deriv(1, 2) = deriv(2, 1) = 0.25 * (stress(1, 2) + stress(2, 1)) / tau;
112  deriv(2, 2) = 0.5 * dsmooth(stress) / tau;
113  }
114  deriv(2, 2) += _tan_phi_or_psi;
115  return deriv;
116 }
virtual Real smooth(const RankTwoTensor &stress) const
returns the &#39;a&#39; parameter - see doco for _tip_scheme
virtual Real dsmooth(const RankTwoTensor &stress) const
returns the da/dstress(2,2) - see doco for _tip_scheme

◆ dflowPotential_dintnl()

RankTwoTensor TensorMechanicsPlasticWeakPlaneShear::dflowPotential_dintnl ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of the flow potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
Returns
dr_dintnl(i, j) = dr(i, j)/dintnl

Reimplemented from TensorMechanicsPlasticModel.

Definition at line 175 of file TensorMechanicsPlasticWeakPlaneShear.C.

177 {
178  RankTwoTensor dr_dintnl;
179  dr_dintnl(2, 2) = dtan_psi(intnl);
180  return dr_dintnl;
181 }
virtual Real dtan_psi(const Real internal_param) const
d(tan_psi)/d(internal_param);

◆ dflowPotential_dintnlV()

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

The derivative of the flow potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
[out]dr_dintnldr_dintnl[alpha](i, j) = dr[alpha](i, j)/dintnl

Reimplemented in TensorMechanicsPlasticMohrCoulombMulti, and TensorMechanicsPlasticTensileMulti.

Definition at line 138 of file TensorMechanicsPlasticModel.C.

141 {
142  return dr_dintnl.assign(1, dflowPotential_dintnl(stress, intnl));
143 }
virtual RankTwoTensor dflowPotential_dintnl(const RankTwoTensor &stress, Real intnl) const
The derivative of the flow potential with respect to the internal parameter.

◆ dflowPotential_dstress()

RankFourTensor TensorMechanicsPlasticWeakPlaneShear::dflowPotential_dstress ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of the flow potential with respect to stress.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
Returns
dr_dstress(i, j, k, l) = dr(i, j)/dstress(k, l)

Reimplemented from TensorMechanicsPlasticModel.

Definition at line 139 of file TensorMechanicsPlasticWeakPlaneShear.C.

141 {
142  RankFourTensor dr_dstress;
143  Real tau = std::sqrt(Utility::pow<2>((stress(0, 2) + stress(2, 0)) / 2.0) +
144  Utility::pow<2>((stress(1, 2) + stress(2, 1)) / 2.0) + smooth(stress));
145  if (tau == 0.0)
146  return dr_dstress;
147 
148  // note that i explicitly symmeterise
149  RankTwoTensor dtau;
150  dtau(0, 2) = dtau(2, 0) = 0.25 * (stress(0, 2) + stress(2, 0)) / tau;
151  dtau(1, 2) = dtau(2, 1) = 0.25 * (stress(1, 2) + stress(2, 1)) / tau;
152  dtau(2, 2) = 0.5 * dsmooth(stress) / tau;
153 
154  for (unsigned i = 0; i < 3; ++i)
155  for (unsigned j = 0; j < 3; ++j)
156  for (unsigned k = 0; k < 3; ++k)
157  for (unsigned l = 0; l < 3; ++l)
158  dr_dstress(i, j, k, l) = -dtau(i, j) * dtau(k, l) / tau;
159 
160  // note that i explicitly symmeterise
161  dr_dstress(0, 2, 0, 2) += 0.25 / tau;
162  dr_dstress(0, 2, 2, 0) += 0.25 / tau;
163  dr_dstress(2, 0, 0, 2) += 0.25 / tau;
164  dr_dstress(2, 0, 2, 0) += 0.25 / tau;
165  dr_dstress(1, 2, 1, 2) += 0.25 / tau;
166  dr_dstress(1, 2, 2, 1) += 0.25 / tau;
167  dr_dstress(2, 1, 1, 2) += 0.25 / tau;
168  dr_dstress(2, 1, 2, 1) += 0.25 / tau;
169  dr_dstress(2, 2, 2, 2) += 0.5 * d2smooth(stress) / tau;
170 
171  return dr_dstress;
172 }
virtual Real d2smooth(const RankTwoTensor &stress) const
returns the d^2a/dstress(2,2)^2 - see doco for _tip_scheme
virtual Real smooth(const RankTwoTensor &stress) const
returns the &#39;a&#39; parameter - see doco for _tip_scheme
virtual Real dsmooth(const RankTwoTensor &stress) const
returns the da/dstress(2,2) - see doco for _tip_scheme

◆ dflowPotential_dstressV()

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

The derivative of the flow potential with respect to stress.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
[out]dr_dstressdr_dstress[alpha](i, j, k, l) = dr[alpha](i, j)/dstress(k, l)

Reimplemented in TensorMechanicsPlasticMohrCoulombMulti, and TensorMechanicsPlasticTensileMulti.

Definition at line 124 of file TensorMechanicsPlasticModel.C.

127 {
128  return dr_dstress.assign(1, dflowPotential_dstress(stress, intnl));
129 }
virtual RankFourTensor dflowPotential_dstress(const RankTwoTensor &stress, Real intnl) const
The derivative of the flow potential with respect to stress.

◆ dhardPotential_dintnl()

Real TensorMechanicsPlasticModel::dhardPotential_dintnl ( const RankTwoTensor &  stress,
Real  intnl 
) const
protectedvirtualinherited

The derivative of the hardening potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
Returns
the derivative

Reimplemented in TensorMechanicsPlasticMeanCapTC.

Definition at line 173 of file TensorMechanicsPlasticModel.C.

Referenced by TensorMechanicsPlasticModel::dhardPotential_dintnlV().

175 {
176  return 0.0;
177 }

◆ dhardPotential_dintnlV()

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

The derivative of the hardening potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
[out]dh_dintnldh_dintnl[alpha] = dh[alpha]/dintnl

Definition at line 179 of file TensorMechanicsPlasticModel.C.

182 {
183  dh_dintnl.resize(numberSurfaces(), dhardPotential_dintnl(stress, intnl));
184 }
virtual unsigned int numberSurfaces() const
The number of yield surfaces for this plasticity model.
virtual Real dhardPotential_dintnl(const RankTwoTensor &stress, Real intnl) const
The derivative of the hardening potential with respect to the internal parameter. ...

◆ dhardPotential_dstress()

RankTwoTensor TensorMechanicsPlasticModel::dhardPotential_dstress ( const RankTwoTensor &  stress,
Real  intnl 
) const
protectedvirtualinherited

The derivative of the hardening potential with respect to stress.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
Returns
dh_dstress(i, j) = dh/dstress(i, j)

Reimplemented in TensorMechanicsPlasticMeanCapTC.

Definition at line 159 of file TensorMechanicsPlasticModel.C.

Referenced by TensorMechanicsPlasticModel::dhardPotential_dstressV().

161 {
162  return RankTwoTensor();
163 }

◆ dhardPotential_dstressV()

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

The derivative of the hardening potential with respect to stress.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
[out]dh_dstressdh_dstress[alpha](i, j) = dh[alpha]/dstress(i, j)

Definition at line 165 of file TensorMechanicsPlasticModel.C.

168 {
169  dh_dstress.assign(numberSurfaces(), dhardPotential_dstress(stress, intnl));
170 }
virtual unsigned int numberSurfaces() const
The number of yield surfaces for this plasticity model.
virtual RankTwoTensor dhardPotential_dstress(const RankTwoTensor &stress, Real intnl) const
The derivative of the hardening potential with respect to stress.

◆ dsmooth()

Real TensorMechanicsPlasticWeakPlaneShear::dsmooth ( const RankTwoTensor &  stress) const
protectedvirtual

returns the da/dstress(2,2) - see doco for _tip_scheme

Definition at line 235 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by df_dsig(), and dflowPotential_dstress().

236 {
237  Real dsmoother2 = 0;
238  if (_tip_scheme == "cap")
239  {
240  Real x = stress(2, 2) - _cap_start;
241  Real p = 0;
242  Real dp_dx = 0;
243  if (x > 0)
244  {
245  const Real exp_cap_rate_x = std::exp(-_cap_rate * x);
246  p = x * (1 - exp_cap_rate_x);
247  dp_dx = (1 - exp_cap_rate_x) + x * _cap_rate * exp_cap_rate_x;
248  }
249  dsmoother2 += 2 * p * dp_dx;
250  }
251  return dsmoother2;
252 }
MooseEnum _tip_scheme
The yield function is modified to f = sqrt(s_xz^2 + s_yz^2 + a) + s_zz*_tan_phi - _cohesion where "a"...
Real _cap_start
smoothing parameter dictating when the &#39;cap&#39; will start - see doco for _tip_scheme ...
Real _cap_rate
dictates how quickly the &#39;cap&#39; degenerates to a hemisphere - see doco for _tip_scheme ...

◆ dtan_phi()

Real TensorMechanicsPlasticWeakPlaneShear::dtan_phi ( const Real  internal_param) const
protectedvirtual

d(tan_phi)/d(internal_param);

Definition at line 202 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by dyieldFunction_dintnl().

203 {
204  return _tan_phi.derivative(internal_param);
205 }
const TensorMechanicsHardeningModel & _tan_phi
Hardening model for tan(phi)
virtual Real derivative(Real intnl) const

◆ dtan_psi()

Real TensorMechanicsPlasticWeakPlaneShear::dtan_psi ( const Real  internal_param) const
protectedvirtual

d(tan_psi)/d(internal_param);

Definition at line 214 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by dflowPotential_dintnl().

215 {
216  return _tan_psi.derivative(internal_param);
217 }
virtual Real derivative(Real intnl) const
const TensorMechanicsHardeningModel & _tan_psi
Hardening model for tan(psi)

◆ dyieldFunction_dintnl()

Real TensorMechanicsPlasticWeakPlaneShear::dyieldFunction_dintnl ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of yield function with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
Returns
the derivative

Reimplemented from TensorMechanicsPlasticModel.

Definition at line 126 of file TensorMechanicsPlasticWeakPlaneShear.C.

128 {
129  return stress(2, 2) * dtan_phi(intnl) - dcohesion(intnl);
130 }
virtual Real dcohesion(const Real internal_param) const
d(cohesion)/d(internal_param)
virtual Real dtan_phi(const Real internal_param) const
d(tan_phi)/d(internal_param);

◆ dyieldFunction_dintnlV()

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

The derivative of yield functions with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
[out]df_dintnldf_dintnl[alpha] = df[alpha]/dintnl

Reimplemented in TensorMechanicsPlasticMohrCoulombMulti, and TensorMechanicsPlasticTensileMulti.

Definition at line 97 of file TensorMechanicsPlasticModel.C.

100 {
101  return df_dintnl.assign(1, dyieldFunction_dintnl(stress, intnl));
102 }
virtual Real dyieldFunction_dintnl(const RankTwoTensor &stress, Real intnl) const
The derivative of yield function with respect to the internal parameter.

◆ dyieldFunction_dstress()

RankTwoTensor TensorMechanicsPlasticWeakPlaneShear::dyieldFunction_dstress ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of yield function with respect to stress.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
Returns
df_dstress(i, j) = dyieldFunction/dstress(i, j)

Reimplemented from TensorMechanicsPlasticModel.

Definition at line 119 of file TensorMechanicsPlasticWeakPlaneShear.C.

121 {
122  return df_dsig(stress, tan_phi(intnl));
123 }
virtual Real tan_phi(const Real internal_param) const
tan_phi as a function of internal parameter
RankTwoTensor df_dsig(const RankTwoTensor &stress, Real _tan_phi_or_psi) const
Function that&#39;s used in dyieldFunction_dstress and flowPotential.

◆ dyieldFunction_dstressV()

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

The derivative of yield functions with respect to stress.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
[out]df_dstressdf_dstress[alpha](i, j) = dyieldFunction[alpha]/dstress(i, j)

Reimplemented in TensorMechanicsPlasticMohrCoulombMulti, and TensorMechanicsPlasticTensileMulti.

Definition at line 83 of file TensorMechanicsPlasticModel.C.

86 {
87  df_dstress.assign(1, dyieldFunction_dstress(stress, intnl));
88 }
virtual RankTwoTensor dyieldFunction_dstress(const RankTwoTensor &stress, Real intnl) const
The derivative of yield function with respect to stress.

◆ execute()

void TensorMechanicsPlasticModel::execute ( )
inherited

Definition at line 46 of file TensorMechanicsPlasticModel.C.

47 {
48 }

◆ finalize()

void TensorMechanicsPlasticModel::finalize ( )
inherited

Definition at line 51 of file TensorMechanicsPlasticModel.C.

52 {
53 }

◆ flowPotential()

RankTwoTensor TensorMechanicsPlasticWeakPlaneShear::flowPotential ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The flow potential.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
Returns
the flow potential

Reimplemented from TensorMechanicsPlasticModel.

Definition at line 133 of file TensorMechanicsPlasticWeakPlaneShear.C.

134 {
135  return df_dsig(stress, tan_psi(intnl));
136 }
virtual Real tan_psi(const Real internal_param) const
tan_psi as a function of internal parameter
RankTwoTensor df_dsig(const RankTwoTensor &stress, Real _tan_phi_or_psi) const
Function that&#39;s used in dyieldFunction_dstress and flowPotential.

◆ flowPotentialV()

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

The flow potentials.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
[out]rr[alpha] is the flow potential for the "alpha" yield function

Reimplemented in TensorMechanicsPlasticMohrCoulombMulti, and TensorMechanicsPlasticTensileMulti.

Definition at line 110 of file TensorMechanicsPlasticModel.C.

113 {
114  return r.assign(1, flowPotential(stress, intnl));
115 }
virtual RankTwoTensor flowPotential(const RankTwoTensor &stress, Real intnl) const
The flow potential.

◆ hardPotential()

Real TensorMechanicsPlasticModel::hardPotential ( const RankTwoTensor &  stress,
Real  intnl 
) const
protectedvirtualinherited

The hardening potential.

Parameters
stressthe stress at which to calculate the hardening potential
intnlinternal parameter
Returns
the hardening potential

Reimplemented in TensorMechanicsPlasticMeanCapTC.

Definition at line 146 of file TensorMechanicsPlasticModel.C.

Referenced by TensorMechanicsPlasticModel::hardPotentialV().

147 {
148  return -1.0;
149 }

◆ hardPotentialV()

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

The hardening potential.

Parameters
stressthe stress at which to calculate the hardening potential
intnlinternal parameter
[out]hh[alpha] is the hardening potential for the "alpha" yield function

Definition at line 151 of file TensorMechanicsPlasticModel.C.

154 {
155  h.assign(numberSurfaces(), hardPotential(stress, intnl));
156 }
virtual Real hardPotential(const RankTwoTensor &stress, Real intnl) const
The hardening potential.
virtual unsigned int numberSurfaces() const
The number of yield surfaces for this plasticity model.

◆ initialize()

void TensorMechanicsPlasticModel::initialize ( )
inherited

Definition at line 41 of file TensorMechanicsPlasticModel.C.

42 {
43 }

◆ KuhnTuckerSingleSurface()

bool TensorMechanicsPlasticModel::KuhnTuckerSingleSurface ( Real  yf,
Real  dpm,
Real  dpm_tol 
) const
inherited

Returns true if the Kuhn-Tucker conditions for the single surface are satisfied.

Parameters
yfYield function value
dpmplastic multiplier
dpm_toltolerance on plastic multiplier: viz dpm>-dpm_tol means "dpm is non-negative"

Definition at line 247 of file TensorMechanicsPlasticModel.C.

Referenced by TensorMechanicsPlasticMohrCoulombMulti::KuhnTuckerOK(), TensorMechanicsPlasticTensileMulti::KuhnTuckerOK(), and TensorMechanicsPlasticModel::returnMap().

248 {
249  return (dpm == 0 && yf <= _f_tol) || (dpm > -dpm_tol && yf <= _f_tol && yf >= -_f_tol);
250 }
const Real _f_tol
Tolerance on yield function.

◆ modelName()

std::string TensorMechanicsPlasticWeakPlaneShear::modelName ( ) const
overridevirtual

Implements TensorMechanicsPlasticModel.

Definition at line 349 of file TensorMechanicsPlasticWeakPlaneShear.C.

350 {
351  return "WeakPlaneShear";
352 }

◆ numberSurfaces()

unsigned TensorMechanicsPlasticModel::numberSurfaces ( ) const
virtualinherited

◆ returnMap()

bool TensorMechanicsPlasticModel::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
virtualinherited

Performs a custom return-map.

You may choose to over-ride this in your derived TensorMechanicsPlasticXXXX class, and you may implement the return-map algorithm in any way that suits you. Eg, using a Newton-Raphson approach, or a radial-return, etc. This may also be used as a quick way of ascertaining whether (trial_stress, intnl_old) is in fact admissible.

For over-riding this function, please note the following.

(1) Denoting the return value of the function by "successful_return", the only possible output values should be: (A) trial_stress_inadmissible=false, successful_return=true. That is, (trial_stress, intnl_old) is in fact admissible (in the elastic domain). (B) trial_stress_inadmissible=true, successful_return=false. That is (trial_stress, intnl_old) is inadmissible (outside the yield surface), and you didn't return to the yield surface. (C) trial_stress_inadmissible=true, successful_return=true. That is (trial_stress, intnl_old) is inadmissible (outside the yield surface), but you did return to the yield surface. The default implementation only handles case (A) and (B): it does not attempt to do a return-map algorithm.

(2) you must correctly signal "successful_return" using the return value of this function. Don't assume the calling function will do Kuhn-Tucker checking and so forth!

(3) In cases (A) and (B) you needn't set returned_stress, returned_intnl, delta_dp, or dpm. This is for computational efficiency.

(4) In cases (A) and (B), you MUST place the yield function values at (trial_stress, intnl_old) into yf so the calling function can use this information optimally. You will have already calculated these yield function values, which can be quite expensive, and it's not very optimal for the calling function to have to re-calculate them.

(5) In case (C), you need to set: returned_stress (the returned value of stress) returned_intnl (the returned value of the internal variable) delta_dp (the change in plastic strain) dpm (the plastic multipliers needed to bring about the return) yf (yield function values at the returned configuration)

(Note, if you over-ride returnMap, you will probably want to override consistentTangentOpertor too, otherwise it will default to E_ijkl.)

Parameters
trial_stressThe trial stress
intnl_oldValue of the internal parameter
E_ijklElasticity tensor
ep_plastic_toleranceTolerance defined by the user for the plastic strain
[out]returned_stressIn case (C): lies on the yield surface after returning and produces the correct plastic strain (normality condition). Otherwise: not defined
[out]returned_intnlIn case (C): the value of the internal parameter after returning. Otherwise: not defined
[out]dpmIn case (C): the plastic multipliers needed to bring about the return. Otherwise: not defined
[out]delta_dpIn case (C): The change in plastic strain induced by the return process. Otherwise: not defined
[out]yfIn case (C): the yield function at (returned_stress, returned_intnl). Otherwise: the yield function at (trial_stress, intnl_old)
[out]trial_stress_inadmissibleShould be set to false if the trial_stress is admissible, and true if the trial_stress is inadmissible. This can be used by the calling prorgram
Returns
true if a successful return (or a return-map not needed), false if the trial_stress is inadmissible but the return process failed

Reimplemented in TensorMechanicsPlasticTensileMulti, TensorMechanicsPlasticMohrCoulombMulti, TensorMechanicsPlasticDruckerPragerHyperbolic, TensorMechanicsPlasticMeanCapTC, and TensorMechanicsPlasticJ2.

Definition at line 220 of file TensorMechanicsPlasticModel.C.

Referenced by TensorMechanicsPlasticJ2::returnMap(), TensorMechanicsPlasticDruckerPragerHyperbolic::returnMap(), TensorMechanicsPlasticMeanCapTC::returnMap(), TensorMechanicsPlasticMohrCoulombMulti::returnMap(), and TensorMechanicsPlasticTensileMulti::returnMap().

230 {
231  trial_stress_inadmissible = false;
232  yieldFunctionV(trial_stress, intnl_old, yf);
233 
234  for (unsigned sf = 0; sf < numberSurfaces(); ++sf)
235  if (yf[sf] > _f_tol)
236  trial_stress_inadmissible = true;
237 
238  // example of checking Kuhn-Tucker
239  std::vector<Real> dpm(numberSurfaces(), 0);
240  for (unsigned sf = 0; sf < numberSurfaces(); ++sf)
241  if (!KuhnTuckerSingleSurface(yf[sf], dpm[sf], 0))
242  return false;
243  return true;
244 }
bool KuhnTuckerSingleSurface(Real yf, Real dpm, Real dpm_tol) const
Returns true if the Kuhn-Tucker conditions for the single surface are satisfied.
virtual unsigned int numberSurfaces() const
The number of yield surfaces for this plasticity model.
virtual void yieldFunctionV(const RankTwoTensor &stress, Real intnl, std::vector< Real > &f) const
Calculates the yield functions.
const Real _f_tol
Tolerance on yield function.

◆ smooth()

Real TensorMechanicsPlasticWeakPlaneShear::smooth ( const RankTwoTensor &  stress) const
protectedvirtual

returns the 'a' parameter - see doco for _tip_scheme

Definition at line 220 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by df_dsig(), dflowPotential_dstress(), and yieldFunction().

221 {
222  Real smoother2 = _small_smoother2;
223  if (_tip_scheme == "cap")
224  {
225  Real x = stress(2, 2) - _cap_start;
226  Real p = 0;
227  if (x > 0)
228  p = x * (1 - std::exp(-_cap_rate * x));
229  smoother2 += Utility::pow<2>(p);
230  }
231  return smoother2;
232 }
MooseEnum _tip_scheme
The yield function is modified to f = sqrt(s_xz^2 + s_yz^2 + a) + s_zz*_tan_phi - _cohesion where "a"...
Real _cap_start
smoothing parameter dictating when the &#39;cap&#39; will start - see doco for _tip_scheme ...
Real _small_smoother2
smoothing parameter for the cone&#39;s tip - see doco for _tip_scheme
Real _cap_rate
dictates how quickly the &#39;cap&#39; degenerates to a hemisphere - see doco for _tip_scheme ...

◆ tan_phi()

Real TensorMechanicsPlasticWeakPlaneShear::tan_phi ( const Real  internal_param) const
protectedvirtual

tan_phi as a function of internal parameter

Definition at line 196 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by activeConstraints(), dyieldFunction_dstress(), TensorMechanicsPlasticWeakPlaneShear(), and yieldFunction().

197 {
198  return _tan_phi.value(internal_param);
199 }
const TensorMechanicsHardeningModel & _tan_phi
Hardening model for tan(phi)
virtual Real value(Real intnl) const

◆ tan_psi()

Real TensorMechanicsPlasticWeakPlaneShear::tan_psi ( const Real  internal_param) const
protectedvirtual

tan_psi as a function of internal parameter

Definition at line 208 of file TensorMechanicsPlasticWeakPlaneShear.C.

Referenced by activeConstraints(), flowPotential(), and TensorMechanicsPlasticWeakPlaneShear().

209 {
210  return _tan_psi.value(internal_param);
211 }
virtual Real value(Real intnl) const
const TensorMechanicsHardeningModel & _tan_psi
Hardening model for tan(psi)

◆ useCustomCTO()

bool TensorMechanicsPlasticModel::useCustomCTO ( ) const
virtualinherited

Returns false. You will want to override this in your derived class if you write a custom consistent tangent operator function.

Reimplemented in TensorMechanicsPlasticTensileMulti, TensorMechanicsPlasticMeanCapTC, TensorMechanicsPlasticDruckerPragerHyperbolic, and TensorMechanicsPlasticJ2.

Definition at line 214 of file TensorMechanicsPlasticModel.C.

215 {
216  return false;
217 }

◆ useCustomReturnMap()

bool TensorMechanicsPlasticModel::useCustomReturnMap ( ) const
virtualinherited

Returns false. You will want to override this in your derived class if you write a custom returnMap function.

Reimplemented in TensorMechanicsPlasticMohrCoulombMulti, TensorMechanicsPlasticTensileMulti, TensorMechanicsPlasticMeanCapTC, TensorMechanicsPlasticDruckerPragerHyperbolic, and TensorMechanicsPlasticJ2.

Definition at line 208 of file TensorMechanicsPlasticModel.C.

209 {
210  return false;
211 }

◆ yieldFunction()

Real TensorMechanicsPlasticWeakPlaneShear::yieldFunction ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The following functions are what you should override when building single-plasticity models.

The yield function

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
Returns
the yield function

Reimplemented from TensorMechanicsPlasticModel.

Definition at line 84 of file TensorMechanicsPlasticWeakPlaneShear.C.

85 {
86  // note that i explicitly symmeterise in preparation for Cosserat
87  return std::sqrt(Utility::pow<2>((stress(0, 2) + stress(2, 0)) / 2.0) +
88  Utility::pow<2>((stress(1, 2) + stress(2, 1)) / 2.0) + smooth(stress)) +
89  stress(2, 2) * tan_phi(intnl) - cohesion(intnl);
90 }
virtual Real tan_phi(const Real internal_param) const
tan_phi as a function of internal parameter
virtual Real smooth(const RankTwoTensor &stress) const
returns the &#39;a&#39; parameter - see doco for _tip_scheme
virtual Real cohesion(const Real internal_param) const
cohesion as a function of internal parameter

◆ yieldFunctionV()

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

Calculates the yield functions.

Note that for single-surface plasticity you don't want to override this - override the private yieldFunction below

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
[out]fthe yield functions

Reimplemented in TensorMechanicsPlasticMohrCoulombMulti, and TensorMechanicsPlasticTensileMulti.

Definition at line 68 of file TensorMechanicsPlasticModel.C.

Referenced by TensorMechanicsPlasticModel::returnMap().

71 {
72  f.assign(1, yieldFunction(stress, intnl));
73 }
virtual Real yieldFunction(const RankTwoTensor &stress, Real intnl) const
The following functions are what you should override when building single-plasticity models...

Member Data Documentation

◆ _cap_rate

Real TensorMechanicsPlasticWeakPlaneShear::_cap_rate
protected

dictates how quickly the 'cap' degenerates to a hemisphere - see doco for _tip_scheme

Definition at line 78 of file TensorMechanicsPlasticWeakPlaneShear.h.

Referenced by d2smooth(), dsmooth(), and smooth().

◆ _cap_start

Real TensorMechanicsPlasticWeakPlaneShear::_cap_start
protected

smoothing parameter dictating when the 'cap' will start - see doco for _tip_scheme

Definition at line 75 of file TensorMechanicsPlasticWeakPlaneShear.h.

Referenced by d2smooth(), dsmooth(), and smooth().

◆ _cohesion

const TensorMechanicsHardeningModel& TensorMechanicsPlasticWeakPlaneShear::_cohesion
protected

Hardening model for cohesion.

Definition at line 41 of file TensorMechanicsPlasticWeakPlaneShear.h.

Referenced by cohesion(), and dcohesion().

◆ _f_tol

const Real TensorMechanicsPlasticModel::_f_tol
inherited

◆ _ic_tol

const Real TensorMechanicsPlasticModel::_ic_tol
inherited

Tolerance on internal constraint.

Definition at line 177 of file TensorMechanicsPlasticModel.h.

◆ _small_smoother2

Real TensorMechanicsPlasticWeakPlaneShear::_small_smoother2
protected

smoothing parameter for the cone's tip - see doco for _tip_scheme

Definition at line 72 of file TensorMechanicsPlasticWeakPlaneShear.h.

Referenced by smooth().

◆ _tan_phi

const TensorMechanicsHardeningModel& TensorMechanicsPlasticWeakPlaneShear::_tan_phi
protected

Hardening model for tan(phi)

Definition at line 44 of file TensorMechanicsPlasticWeakPlaneShear.h.

Referenced by dtan_phi(), and tan_phi().

◆ _tan_psi

const TensorMechanicsHardeningModel& TensorMechanicsPlasticWeakPlaneShear::_tan_psi
protected

Hardening model for tan(psi)

Definition at line 47 of file TensorMechanicsPlasticWeakPlaneShear.h.

Referenced by dtan_psi(), and tan_psi().

◆ _tip_scheme

MooseEnum TensorMechanicsPlasticWeakPlaneShear::_tip_scheme
protected

The yield function is modified to f = sqrt(s_xz^2 + s_yz^2 + a) + s_zz*_tan_phi - _cohesion where "a" depends on the tip_scheme.

Currently _tip_scheme is 'hyperbolic', where a = _small_smoother2 'cap' where a = _small_smoother2 + (p(stress(2,2) - _cap_start))^2 with the function p(x)=x(1-exp(-_cap_rate*x)) for x>0, and p=0 otherwise

Definition at line 69 of file TensorMechanicsPlasticWeakPlaneShear.h.

Referenced by d2smooth(), dsmooth(), and smooth().


The documentation for this class was generated from the following files: