16 #ifndef dealii_elasticity_standard_tensors_h 17 #define dealii_elasticity_standard_tensors_h 20 #include <deal.II/base/numbers.h> 22 #include <deal.II/base/tensor.h> 24 DEAL_II_NAMESPACE_OPEN
214 template <
typename Number>
229 template <
typename Number>
256 template <
typename Number>
281 template <
typename Number>
299 template <
typename Number>
306 ExcMessage(
"Deformation gradient has a negative determinant."));
317 std::pow(det_F, -2.0 / dim);
325 template <
typename Number>
332 ExcMessage(
"Deformation gradient has a negative determinant."));
343 std::pow(det_F, -2.0 / dim);
351 template <
typename Number>
362 template <
typename Number>
371 for (
unsigned int A = 0; A < dim; ++A)
372 for (
unsigned int B = A; B < dim; ++B)
373 for (
unsigned int C = 0; C < dim; ++C)
374 for (
unsigned int D = C; D < dim; ++D)
375 dC_inv_dC[A][B][C][D] -=
376 0.5 * (C_inv[A][C] * C_inv[B][D] + C_inv[A][D] * C_inv[B][C]);
383 DEAL_II_NAMESPACE_CLOSE
Number determinant(const SymmetricTensor< 2, dim, Number > &)
SymmetricTensor< 4, dim, Number > outer_product(const SymmetricTensor< 2, dim, Number > &t1, const SymmetricTensor< 2, dim, Number > &t2)
SymmetricTensor< 2, dim, Number > invert(const SymmetricTensor< 2, dim, Number > &)
static::ExceptionBase & ExcMessage(std::string arg1)
static const SymmetricTensor< 4, dim > S
#define Assert(cond, exc)
static SymmetricTensor< 2, dim, Number > ddet_F_dC(const Tensor< 2, dim, Number > &F)
static const SymmetricTensor< 2, dim > I
static const SymmetricTensor< 4, dim > dev_P
SymmetricTensor< 2, dim, Number > symmetrize(const Tensor< 2, dim, Number > &t)
SymmetricTensor< rank_, dim, Number > transpose(const SymmetricTensor< rank_, dim, Number > &t)
static const SymmetricTensor< 4, dim > IxI
bool value_is_greater_than(const Number1 &value_1, const Number2 &value_2)
static SymmetricTensor< 4, dim, Number > Dev_P(const Tensor< 2, dim, Number > &F)
static SymmetricTensor< 4, dim, Number > dC_inv_dC(const Tensor< 2, dim, Number > &F)
static SymmetricTensor< 4, dim, Number > Dev_P_T(const Tensor< 2, dim, Number > &F)