From a8d8276fbaad37766b62ba0e133925aa26e216ad Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 13 Jul 2011 10:13:37 +0000
Subject: [PATCH] yet another continuum (== non-shell) bugfix
[[Imported from SVN: r7560]]
---
dune/gfe/cosseratenergystiffness.hh | 10 ++++++++--
1 file changed, 8 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index de1a5e91..0257b577 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -205,14 +205,20 @@ energy(const Entity& element,
Dune::FieldMatrix<double,dim,dim> R;
value.q.matrix(R);
- /// \f$ \hat{F} = (\nabla m | \overline{R}_3) \f$
+ /* Compute F, the deformation gradient.
+ In the continuum case this is
+ \f$ \hat{F} = \nabla m \f$
+ In the case of a shell it is
+ \f$ \hat{F} = (\nabla m | \overline{R}_3) \f$
+ */
Dune::FieldMatrix<double,dim,dim> F;
for (int i=0; i<dim; i++)
for (int j=0; j<gridDim; j++)
F[i][j] = derivative[i][j];
for (int i=0; i<dim; i++)
- F[i][dim-1] = R[i][dim-1];
+ for (int j=gridDim; j<dim; j++)
+ F[i][j] = R[i][j];
// U = R^T F
Dune::FieldMatrix<double,dim,dim> U;
--
GitLab