11 #include "libmesh/utility.h" 19 params.
addClassDescription(
"Capped Mohr-Coulomb plasticity stress calculator for the Cosserat " 20 "situation where the host medium (ie, the limit where all Cosserat " 21 "effects are zero) is isotropic. Note that the return-map flow rule " 22 "uses an isotropic elasticity tensor built with the 'host' properties " 23 "defined by the user.");
25 "host_youngs_modulus>0",
26 "Young's modulus for the isotropic host medium");
28 "host_poissons_ratio>=0 & host_poissons_ratio<0.5",
29 "Poisson's ratio for the isotropic host medium");
36 _host_young(getParam<
Real>(
"host_youngs_modulus")),
37 _host_poisson(getParam<
Real>(
"host_poissons_ratio")),
38 _host_E0011(_host_young * _host_poisson / (1.0 + _host_poisson) / (1.0 - 2.0 * _host_poisson)),
39 _host_E0000(_host_E0011 + _host_young / (1.0 + _host_poisson))
45 const std::vector<Real> & ,
47 const std::vector<Real> & ,
48 const std::vector<Real> & ,
51 std::vector<Real> eigvals;
66 for (
unsigned b = 0;
b <
a; ++
b)
74 const std::vector<Real> & stress_params,
76 const std::vector<Real> & ,
82 stress =
RankTwoTensor(stress_params[0], stress_params[1], stress_params[2], 0.0, 0.0, 0.0);
86 stress += 0.5 * (stress_trial - stress_trial.
transpose());
92 const std::vector<Real> & trial_stress_params,
94 const std::vector<Real> & stress_params,
98 bool compute_full_tangent_operator,
99 const std::vector<std::vector<Real>> & dvar_dtrial,
109 compute_full_tangent_operator,
113 if (!compute_full_tangent_operator)
const Real _host_E0011
E0011 = Lame lambda modulus of the host material.
CappedMohrCoulombStressUpdate implements rate-independent nonassociative Mohr-Coulomb plus tensile pl...
void setEffectiveElasticity(const RankFourTensor &Eijkl) override
Sets _Eij and _En and _Cij.
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.
Real _En
normalising factor
Struct designed to hold info about a single yield function and its derivatives, as well as the flow d...
std::vector< std::vector< Real > > _Cij
_Cij[i, j] * _Eij[j, k] = 1 iff j == k
registerMooseObject("SolidMechanicsApp", CappedMohrCoulombCosseratStressUpdate)
CappedMohrCoulombCosseratStressUpdate(const InputParameters ¶meters)
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.
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.
static InputParameters validParams()
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...
const Real _host_poisson
Poisson's of the host material.
CappedMohrCoulombCosseratStressUpdate implements rate-independent nonassociative Mohr-Coulomb plus te...
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.
RankTwoTensorTempl< Real > transpose() const
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static InputParameters validParams()
void symmetricEigenvaluesEigenvectors(std::vector< Real > &eigvals, RankTwoTensorTempl< Real > &eigvecs) const
RankTwoTensor _eigvecs
Eigenvectors of the trial stress as a RankTwoTensor, in order to rotate the returned stress back to s...
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
const Real _host_E0000
E0000 = Lame lambda + 2 * shear modulus of the host material.
static constexpr unsigned _tensor_dimensionality
Internal dimensionality of tensors (currently this is 3 throughout tensor_mechanics) ...
static const std::string k
Real _poissons_ratio
Poisson's ratio.
const unsigned _num_sp
Number of stress parameters.