From c10a79bccf4a8f6a7d7b4aa4bfd20c265ab3f591 Mon Sep 17 00:00:00 2001
From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de>
Date: Mon, 11 Apr 2022 23:54:45 +0200
Subject: [PATCH] Replace traceSquared(sym(.)) by traceSquared(.) as it is the
same
---
dune/gfe/cosseratenergystiffness.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index af54e5f1..9d402a24 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -131,7 +131,7 @@ public:
const std::vector<Rotation<field_type,dim> >& localOrientationConfiguration) const override;
/** \brief The energy \f$ W_{mp}(\overline{U}) \f$, as written in
- * the first equation of (4.4) in Neff's paper
+ * the first equation of (4.4) in Neff's paper from 2006: A geometrically exact planar Cosserat shell model with microstructure: Existence of minimizers for zero Cosserat couple modulus
*/
RT quadraticMembraneEnergy(const Dune::GFE::CosseratStrain<field_type,3,gridDim>& U) const
{
@@ -141,7 +141,7 @@ public:
return mu_ * Dune::GFE::sym(UMinus1).frobenius_norm2()
+ mu_c_ * Dune::GFE::skew(UMinus1).frobenius_norm2()
- + (mu_*lambda_)/(2*mu_ + lambda_) * Dune::GFE::traceSquared(Dune::GFE::sym(UMinus1));
+ + (mu_*lambda_)/(2*mu_ + lambda_) * Dune::GFE::traceSquared(UMinus1); // Dune::GFE::traceSquared(UMinus1) = Dune::GFE::traceSquared(Dune::GFE::sym(UMinus1))
}
/** \brief The energy \f$ W_{mp}(\overline{U}) \f$, as written in
--
GitLab