diff --git a/dune/fufem/mechanics/isotropictensor.hh b/dune/fufem/mechanics/isotropictensor.hh index ba25117834e78c983871320e73e27cabcb70b864..5feeefd3d411bd8ec59943294737ede6f344fa2b 100644 --- a/dune/fufem/mechanics/isotropictensor.hh +++ b/dune/fufem/mechanics/isotropictensor.hh @@ -17,6 +17,15 @@ template <> class IsotropicTensor<3> : public ElasticityTensor<3> { public: + /* Same representation as in symmetrictensor.hh: + + [s11] [ * * * ] [e11] + [s22] [ * * * ] [e22] + [s33] = [ * * * ] * [e33] + [s12] [ * ] [e12] + [s13] [ * ] [e13] + [s23] [ * ] [e23] + */ IsotropicTensor(double E, double nu) { ElasticityTensor<3>::operator=(0.0); @@ -33,9 +42,9 @@ public: (*this)[2][1] = nu; (*this)[2][2] = 1 - nu; - (*this)[3][3] = 0.5 - nu; - (*this)[4][4] = 0.5 - nu; - (*this)[5][5] = 0.5 - nu; + (*this)[3][3] = 1 - 2*nu; + (*this)[4][4] = 1 - 2*nu; + (*this)[5][5] = 1 - 2*nu; (*this) *= E/(1 + nu)/(1 - 2*nu); } @@ -55,7 +64,7 @@ public: (*this)[1][0] = nu; (*this)[1][1] = 1 - nu; - (*this)[2][2] = 0.5 - nu; + (*this)[2][2] = 1 - 2*nu; (*this) *= (E/(1 + nu)/(1 - 2*nu)); }