11 #include "libmesh/utility.h" 21 "A SolidMechanicsHardening UserObject that defines hardening of the tensile strength");
22 params.
addParam<
bool>(
"perfect_guess",
24 "Provide a guess to the Newton-Raphson procedure " 25 "that is the result from perfect plasticity. With " 26 "severe hardening/softening this may be " 29 "Associative, smoothed, tensile (Rankine) plasticity with hardening/softening");
36 _perfect_guess(getParam<bool>(
"perfect_guess")),
43 std::vector<Real> & stress_params)
const 49 std::vector<RankTwoTensor>
53 std::vector<RankTwoTensor> dsp;
58 std::vector<RankFourTensor>
61 std::vector<RankFourTensor> d2;
69 const std::vector<Real> & ,
70 const std::vector<Real> & ,
73 std::vector<Real> eigvals;
79 const std::vector<Real> & stress_params,
81 const std::vector<Real> & ,
87 stress =
RankTwoTensor(stress_params[0], stress_params[1], stress_params[2], 0.0, 0.0, 0.0);
94 const std::vector<Real> & intnl,
95 std::vector<Real> & yf)
const 102 yf[0] = stress_params[2] - ts;
103 yf[1] = stress_params[1] - ts;
104 yf[2] = stress_params[0] - ts;
109 const std::vector<Real> & intnl,
110 std::vector<yieldAndFlow> & all_q)
const 116 all_q[0].f = stress_params[2] - ts;
117 all_q[1].f = stress_params[1] - ts;
118 all_q[2].f = stress_params[0] - ts;
121 for (
unsigned yf = 0; yf <
_num_yf; ++yf)
123 all_q[yf].df[
a] = 0.0;
124 all_q[0].df[2] = 1.0;
125 all_q[1].df[1] = 1.0;
126 all_q[2].df[0] = 1.0;
129 for (
unsigned yf = 0; yf <
_num_yf; ++yf)
130 all_q[yf].df_di[0] = -dts;
134 for (
unsigned yf = 0; yf <
_num_yf; ++yf)
136 all_q[yf].dg[
a] = all_q[yf].df[
a];
139 for (
unsigned yf = 0; yf <
_num_yf; ++yf)
141 all_q[yf].d2g_di[
a][0] = 0.0;
144 for (
unsigned yf = 0; yf <
_num_yf; ++yf)
147 all_q[yf].d2g[
a][
b] = 0.0;
163 for (
unsigned b = 0;
b <
a; ++
b)
170 const std::vector<Real> & intnl_old,
171 std::vector<Real> & stress_params,
173 std::vector<Real> & intnl)
const 177 for (
unsigned i = 0; i <
_num_sp; ++i)
178 stress_params[i] = trial_stress_params[i];
184 stress_params[2] = ts;
185 stress_params[1] = std::min(stress_params[1], ts);
186 stress_params[0] = std::min(stress_params[0], ts);
187 gaE = trial_stress_params[2] - stress_params[2];
194 const std::vector<Real> & current_stress_params,
195 const std::vector<Real> & intnl_old,
196 std::vector<Real> & intnl)
const 198 intnl[0] = intnl_old[0] + (trial_stress_params[2] - current_stress_params[2]) /
_Eij[2][2];
203 const std::vector<Real> & ,
204 const std::vector<Real> & ,
205 std::vector<std::vector<Real>> & dintnl)
const 209 dintnl[0][2] = -1.0 /
_Eij[2][2];
226 const std::vector<Real> & trial_stress_params,
228 const std::vector<Real> & stress_params,
232 bool compute_full_tangent_operator,
233 const std::vector<std::vector<Real>> & dvar_dtrial,
237 if (!compute_full_tangent_operator)
253 if (trial_stress_params[
a] == trial_stress_params[
j])
255 drot_dstress(i,
j,
k, l) +=
258 (trial_stress_params[
j] - trial_stress_params[
a]);
269 dstress_dtrial(i,
j,
k, l) +=
270 drot_dstress(i,
a,
k, l) * stress_params[
a] * eT(
a,
j) +
271 _eigvecs(i,
a) * stress_params[
a] * drot_dstress(
j,
a,
k, l);
280 dstress_dtrial(i,
j,
k, l) +=
TensileStressUpdate(const InputParameters ¶meters)
registerMooseObject("SolidMechanicsApp", TensileStressUpdate)
virtual void preReturnMapV(const std::vector< Real > &trial_stress_params, 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...
void symmetricEigenvalues(std::vector< Real > &eigvals) const
virtual void consistentTangentOperatorV(const RankTwoTensor &stress_trial, const std::vector< Real > &trial_stress_params, const RankTwoTensor &stress, const std::vector< Real > &stress_params, Real gaE, const yieldAndFlow &smoothed_q, const RankFourTensor &Eijkl, bool compute_full_tangent_operator, const std::vector< std::vector< Real >> &dvar_dtrial, RankFourTensor &cto) override
Calculates the consistent tangent operator.
virtual Real dtensile_strength(const Real internal_param) const
d(tensile strength)/d(internal_param) as a function of residual value, rate, and internal_param ...
void dsymmetricEigenvalues(std::vector< Real > &eigvals, std::vector< RankTwoTensorTempl< Real >> &deigvals) const
Real _En
normalising factor
Struct designed to hold info about a single yield function and its derivatives, as well as the flow d...
TensileStressUpdate implements rate-independent associative tensile failure ("Rankine" plasticity) wi...
std::vector< std::vector< Real > > _Cij
_Cij[i, j] * _Eij[j, k] = 1 iff j == k
void setEffectiveElasticity(const RankFourTensor &Eijkl) override
Sets _Eij and _En and _Cij.
void initializeVarsV(const std::vector< Real > &trial_stress_params, const std::vector< Real > &intnl_old, std::vector< Real > &stress_params, Real &gaE, std::vector< Real > &intnl) const override
Sets (stress_params, intnl) at "good guesses" of the solution to the Return-Map algorithm.
virtual Real value(Real intnl) const
const SolidMechanicsHardeningModel & _strength
Hardening model for tensile strength.
RankTwoTensor _eigvecs
Eigenvectors of the trial stress as a RankTwoTensor, in order to rotate the returned stress back to s...
void setIntnlDerivativesV(const std::vector< Real > &trial_stress_params, const std::vector< Real > ¤t_stress_params, 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 stress_params, their current values, and the current values of the internal parameters.
void computeAllQV(const std::vector< Real > &stress_params, const std::vector< Real > &intnl, std::vector< yieldAndFlow > &all_q) const override
Completely fills all_q with correct values.
Real elasticity_tensor(unsigned int i, unsigned int j, unsigned int k, unsigned int l)
std::vector< std::vector< Real > > _Eij
E[i, j] in the system of equations to be solved.
const unsigned _num_yf
Number of yield functions.
void d2symmetricEigenvalues(std::vector< RankFourTensorTempl< Real >> &deriv) const
void yieldFunctionValuesV(const std::vector< Real > &stress_params, const std::vector< Real > &intnl, std::vector< Real > &yf) const override
Computes the values of the yield functions, given stress_params and intnl parameters.
virtual Real derivative(Real intnl) const
RankTwoTensorTempl< Real > transpose() const
Hardening Model base class.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::vector< RankFourTensor > d2stress_param_dstress(const RankTwoTensor &stress) const override
d2(stress_param[i])/d(stress)/d(stress) at given stress
void symmetricEigenvaluesEigenvectors(std::vector< Real > &eigvals, RankTwoTensorTempl< Real > &eigvecs) const
const bool _perfect_guess
Whether to provide an estimate of the returned stress, based on perfect plasticity.
MultiParameterPlasticityStressUpdate performs the return-map algorithm and associated stress updates ...
static InputParameters validParams()
std::vector< RankTwoTensor > dstress_param_dstress(const RankTwoTensor &stress) const override
d(stress_param[i])/d(stress) at given stress
void setIntnlValuesV(const std::vector< Real > &trial_stress_params, const std::vector< Real > ¤t_stress_params, const std::vector< Real > &intnl_old, std::vector< Real > &intnl) const override
Sets the internal parameters based on the trial values of stress_params, their current values...
static InputParameters validParams()
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
void computeStressParams(const RankTwoTensor &stress, std::vector< Real > &stress_params) const override
Computes stress_params, given stress.
virtual Real tensile_strength(const Real internal_param) const
tensile strength as a function of residual value, rate, and internal_param
static constexpr unsigned _tensor_dimensionality
Internal dimensionality of tensors (currently this is 3 throughout tensor_mechanics) ...
static const std::string k
virtual void setStressAfterReturnV(const RankTwoTensor &stress_trial, const std::vector< Real > &stress_params, Real gaE, const std::vector< Real > &intnl, const yieldAndFlow &smoothed_q, const RankFourTensor &Eijkl, RankTwoTensor &stress) const override
Sets stress from the admissible parameters.
const unsigned _num_sp
Number of stress parameters.