12 #include "libmesh/utility.h" 37 mooseError(
"CappedWeakInclinedPlaneStressUpdate: normal_vector must not have zero length");
80 const std::vector<Real> & ,
81 const std::vector<Real> & yf,
107 p = rotated_stress(2, 2);
108 q =
std::sqrt(Utility::pow<2>(rotated_stress(0, 2)) + Utility::pow<2>(rotated_stress(1, 2)));
125 const std::vector<Real> & ,
138 stress(2, 0) = stress(2, 1) = stress(0, 2) = stress(1, 2) = 0.0;
159 bool compute_full_tangent_operator,
171 compute_full_tangent_operator,
RealTensorValue _rot_n_to_z
Rotation matrix that rotates _n to "z".
RealVectorValue _n_input
User-input value of the normal vector to the weak plane.
virtual RankTwoTensor dqdstress(const RankTwoTensor &stress) const override
d(q)/d(stress) Derived classes must override this
auto norm() const -> decltype(std::norm(Real()))
virtual void initQpStatefulProperties() override
virtual RankFourTensor d2qdstress2(const RankTwoTensor &stress) const override
d2(q)/d(stress)/d(stress) Derived classes must override this
RankFourTensor _rotated_Eijkl
Elasticity tensor rotated to the frame where _n points along "z".
RankTwoTensor _rotated_trial
Trial stress rotated to the frame where _n points along "z".
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
Calculates the consistent tangent operator.
CappedWeakInclinedPlaneStressUpdate(const InputParameters ¶meters)
Struct designed to hold info about a single yield function and its derivatives, as well as the flow d...
static constexpr std::size_t dim
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 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.
ADRealEigenVector< T, D, asd > sqrt(const ADRealEigenVector< T, D, asd > &)
RealTensorValue _rot_z_to_n
Rotation matrix that rotates "z" to _n.
void rotate(const RankTwoTensorTempl< Real > &R)
TensorValue< Real > RealTensorValue
Real _in_trial02
trial value of stress(0, 2)
virtual void initializeReturnProcess() override
Derived classes may use this to perform calculations before any return-map process is performed...
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...
CappedWeakInclinedPlaneStressUpdate performs the return-map algorithm and associated stress updates f...
virtual void finalizeReturnProcess(const RankTwoTensor &rotation_increment) override
Derived classes may use this to perform calculations after the return-map process has completed succe...
virtual RankTwoTensor dqdstress(const RankTwoTensor &stress) const override
d(q)/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 initQpStatefulProperties() override
Real _Epp
elasticity tensor in p direction
registerMooseObject("SolidMechanicsApp", CappedWeakInclinedPlaneStressUpdate)
void rotate(const TypeTensor< T > &R)
const bool _perform_finite_strain_rotations
Whether to perform finite-strain rotations.
GenericRealTensorValue< is_ad > rotVecToZ(GenericRealVectorValue< is_ad > vec)
static InputParameters validParams()
virtual void finalizeReturnProcess(const RankTwoTensor &rotation_increment) override
Derived classes may use this to perform calculations after the return-map process has completed succe...
Real _in_q_trial
trial value of q
enum CappedWeakPlaneStressUpdate::StressReturnType _stress_return_type
IntRange< T > make_range(T beg, T end)
void mooseError(Args &&... args) const
static InputParameters validParams()
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
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 MaterialProperty< RealVectorValue > & _n_old
Old value of the normal.
virtual RankFourTensor d2qdstress2(const RankTwoTensor &stress) const override
d2(q)/d(stress)/d(stress) Derived classes must override this
CappedWeakPlaneStressUpdate performs the return-map algorithm and associated stress updates for plast...
virtual void initializeReturnProcess() override
Derived classes may use this to perform calculations before any return-map process is performed...
MaterialProperty< RealVectorValue > & _n
Current value of the normal.
virtual RankTwoTensor dpdstress(const RankTwoTensor &stress) const override
d(p)/d(stress) Derived classes must override this
Real _in_trial12
trial value of stress(1, 2)