20template <
int DisplacementDim>
21std::optional<std::tuple<typename MechanicsBase<DisplacementDim>::KelvinVector,
22 std::unique_ptr<
typename MechanicsBase<
23 DisplacementDim>::MaterialStateVariables>,
32 auto const& eps_m = std::get<MPL::SymmetricTensor<DisplacementDim>>(
34 auto const& eps_m_prev = std::get<MPL::SymmetricTensor<DisplacementDim>>(
36 auto const& sigma_prev = std::get<MPL::SymmetricTensor<DisplacementDim>>(
37 variable_array_prev.
stress);
42 KelvinVector sigma = sigma_prev + C * (eps_m - eps_m_prev);
44 return {std::make_tuple(
51template <
int DisplacementDim>
57 return elasticTangentStiffness<DisplacementDim>(_mp.lambda(t, x),
std::optional< std::tuple< typename MechanicsBase< DisplacementDim >::KelvinVector, std::unique_ptr< typename MechanicsBase< DisplacementDim >::MaterialStateVariables >, typename MechanicsBase< DisplacementDim >::KelvinMatrix > > integrateStress(MaterialPropertyLib::VariableArray const &variable_array_prev, MaterialPropertyLib::VariableArray const &variable_array, double const t, ParameterLib::SpatialPosition const &x, double const, typename MechanicsBase< DisplacementDim >::MaterialStateVariables const &material_state_variables) const override
MathLib::KelvinVector::KelvinVectorType< DisplacementDim > KelvinVector
KelvinMatrix getElasticTensor(double const t, ParameterLib::SpatialPosition const &x, double const T) const
MathLib::KelvinVector::KelvinMatrixType< DisplacementDim > KelvinMatrix
std::variant< std::monostate, Eigen::Matrix< double, 4, 1 >, Eigen::Matrix< double, 6, 1 > > mechanical_strain
std::variant< std::monostate, Eigen::Matrix< double, 4, 1 >, Eigen::Matrix< double, 6, 1 > > stress
MathLib::KelvinVector::KelvinMatrixType< DisplacementDim > KelvinMatrix