www.mooseframework.org
CappedDruckerPragerStressUpdate.h
Go to the documentation of this file.
1 /****************************************************************/
2 /* MOOSE - Multiphysics Object Oriented Simulation Environment */
3 /* */
4 /* All contents are licensed under LGPL V2.1 */
5 /* See LICENSE for full restrictions */
6 /****************************************************************/
7 #ifndef CAPPEDDRUCKERPRAGERSTRESSUPDATE_H
8 #define CAPPEDDRUCKERPRAGERSTRESSUPDATE_H
9 
13 
15 
16 template <>
18 
59 {
60 public:
61  CappedDruckerPragerStressUpdate(const InputParameters & parameters);
62 
66  bool requiresIsotropicTensor() override { return true; }
67 
68 protected:
71 
74 
77 
79  const Real _small_smoother2;
80 
82  const bool _perfect_guess;
83 
97  enum class StressReturnType
98  {
101  no_tension
103 
109  const bool _small_dilation;
110 
113 
114  virtual void yieldFunctionValues(Real p,
115  Real q,
116  const std::vector<Real> & intnl,
117  std::vector<Real> & yf) const override;
118 
119  virtual void computeAllQ(Real p,
120  Real q,
121  const std::vector<Real> & intnl,
122  std::vector<yieldAndFlow> & all_q) const override;
123 
124  virtual void preReturnMap(Real p_trial,
125  Real q_trial,
126  const RankTwoTensor & stress_trial,
127  const std::vector<Real> & intnl_old,
128  const std::vector<Real> & yf,
129  const RankFourTensor & Eijkl) override;
130 
131  virtual void initializeVars(Real p_trial,
132  Real q_trial,
133  const std::vector<Real> & intnl_old,
134  Real & p,
135  Real & q,
136  Real & gaE,
137  std::vector<Real> & intnl) const override;
138 
139  virtual void setIntnlValues(Real p_trial,
140  Real q_trial,
141  Real p,
142  Real q,
143  const std::vector<Real> & intnl_old,
144  std::vector<Real> & intnl) const override;
145 
146  virtual void setIntnlDerivatives(Real p_trial,
147  Real q_trial,
148  Real p,
149  Real q,
150  const std::vector<Real> & intnl,
151  std::vector<std::vector<Real>> & dintnl) const override;
152 
153  virtual void computePQ(const RankTwoTensor & stress, Real & p, Real & q) const override;
154 
155  virtual void initializeReturnProcess() override;
156 
157  virtual void finalizeReturnProcess(const RankTwoTensor & rotation_increment) override;
158 
159  virtual void setEppEqq(const RankFourTensor & Eijkl, Real & Epp, Real & Eqq) const override;
160 
161  virtual void setStressAfterReturn(const RankTwoTensor & stress_trial,
162  Real p_ok,
163  Real q_ok,
164  Real gaE,
165  const std::vector<Real> & intnl,
166  const yieldAndFlow & smoothed_q,
167  const RankFourTensor & Eijkl,
168  RankTwoTensor & stress) const override;
169 
170  virtual void consistentTangentOperator(const RankTwoTensor & stress_trial,
171  Real p_trial,
172  Real q_trial,
173  const RankTwoTensor & stress,
174  Real p,
175  Real q,
176  Real gaE,
177  const yieldAndFlow & smoothed_q,
178  const RankFourTensor & Eijkl,
179  bool compute_full_tangent_operator,
180  RankFourTensor & cto) const override;
181 
182  virtual RankTwoTensor dpdstress(const RankTwoTensor & stress) const override;
183 
184  virtual RankFourTensor d2pdstress2(const RankTwoTensor & stress) const override;
185 
186  virtual RankTwoTensor dqdstress(const RankTwoTensor & stress) const override;
187 
188  virtual RankFourTensor d2qdstress2(const RankTwoTensor & stress) const override;
189 };
190 
191 #endif // CAPPEDDRUCKERPRAGERSTRESSUPDATE_H
const Real _small_smoother2
The cone vertex is smoothed by this amount.
virtual void setIntnlDerivatives(Real p_trial, Real q_trial, Real p, Real q, const std::vector< Real > &intnl, std::vector< std::vector< Real >> &dintnl) const override
Sets the derivatives of internal parameters, based on the trial values of p and q, their current values, and the old values of the internal parameters.
virtual RankTwoTensor dpdstress(const RankTwoTensor &stress) const override
d(p)/d(stress) Derived classes must override this
virtual void computePQ(const RankTwoTensor &stress, Real &p, Real &q) const override
Computes p and q, given stress.
virtual void yieldFunctionValues(Real p, Real q, const std::vector< Real > &intnl, std::vector< Real > &yf) const override
Computes the values of the yield functions, given p, q and intnl parameters.
virtual void initializeVars(Real p_trial, Real q_trial, const std::vector< Real > &intnl_old, Real &p, Real &q, Real &gaE, std::vector< Real > &intnl) const override
Sets (p, q, gaE, intnl) at "good guesses" of the solution to the Return-Map algorithm.
Struct designed to hold info about a single yield function and its derivatives, as well as the flow d...
TwoParameterPlasticityStressUpdate performs the return-map algorithm and associated stress updates fo...
virtual void finalizeReturnProcess(const RankTwoTensor &rotation_increment) override
Derived classes may use this to perform calculations after the return-map process has completed succe...
const bool _small_dilation
If true, and if the trial stress exceeds the tensile strength, then the user gaurantees that the retu...
InputParameters validParams< CappedDruckerPragerStressUpdate >()
const TensorMechanicsHardeningModel & _tstrength
Hardening model for tensile strength.
CappedDruckerPragerStressUpdate(const InputParameters &parameters)
const TensorMechanicsHardeningModel & _cstrength
Hardening model for compressive strength.
virtual void setIntnlValues(Real p_trial, Real q_trial, Real p, Real q, const std::vector< Real > &intnl_old, std::vector< Real > &intnl) const override
Sets the internal parameters based on the trial values of p and q, their current values, and the old values of the internal parameters.
Rate-independent non-associative Drucker Prager with hardening/softening.
virtual RankFourTensor d2pdstress2(const RankTwoTensor &stress) const override
d2(p)/d(stress)/d(stress) Derived classes must override this
virtual RankTwoTensor dqdstress(const RankTwoTensor &stress) const override
d(q)/d(stress) Derived classes must override this
virtual RankFourTensor d2qdstress2(const RankTwoTensor &stress) const override
d2(q)/d(stress)/d(stress) Derived classes must override this
const TensorMechanicsPlasticDruckerPrager & _dp
Hardening model for cohesion, friction and dilation angles.
bool requiresIsotropicTensor() override
Does the model require the elasticity tensor to be isotropic?
virtual void consistentTangentOperator(const RankTwoTensor &stress_trial, Real p_trial, Real q_trial, const RankTwoTensor &stress, Real p, Real q, Real gaE, const yieldAndFlow &smoothed_q, const RankFourTensor &Eijkl, bool compute_full_tangent_operator, RankFourTensor &cto) const override
Calculates the consistent tangent operator.
virtual void initializeReturnProcess() override
Derived classes may use this to perform calculations before any return-map process is performed...
virtual void setStressAfterReturn(const RankTwoTensor &stress_trial, Real p_ok, Real q_ok, Real gaE, const std::vector< Real > &intnl, const yieldAndFlow &smoothed_q, const RankFourTensor &Eijkl, RankTwoTensor &stress) const override
Sets stress from the admissible parameters.
virtual void computeAllQ(Real p, Real q, const std::vector< Real > &intnl, std::vector< yieldAndFlow > &all_q) const override
Completely fills all_q with correct values.
StressReturnType
This allows some simplification in the return-map process.
CappedDruckerPragerStressUpdate performs the return-map algorithm and associated stress updates for p...
enum CappedDruckerPragerStressUpdate::StressReturnType _stress_return_type
virtual void setEppEqq(const RankFourTensor &Eijkl, Real &Epp, Real &Eqq) const override
Set Epp and Eqq based on the elasticity tensor Derived classes must override this.
const bool _perfect_guess
Initialize the NR proceedure from a guess coming from perfect plasticity.
virtual void preReturnMap(Real p_trial, Real q_trial, const RankTwoTensor &stress_trial, const std::vector< Real > &intnl_old, const std::vector< Real > &yf, const RankFourTensor &Eijkl) override
Derived classes may employ this function to record stuff or do other computations prior to the return...