diff --git a/dune/elasticity/materials/mooneyrivlindensity.hh b/dune/elasticity/materials/mooneyrivlindensity.hh index 0e7a76658c6be1d2db1bb9b2ac19eef0d6a7389e..6871bc92b168156a50038d23d92057ba7d7f9339 100644 --- a/dune/elasticity/materials/mooneyrivlindensity.hh +++ b/dune/elasticity/materials/mooneyrivlindensity.hh @@ -95,8 +95,7 @@ public: { using std::log; return mooneyrivlin_a*normFSquared + mooneyrivlin_b*normFinvSquared*detF + mooneyrivlin_c*detF*detF - ((dim-1)*mooneyrivlin_a + mooneyrivlin_b + 2*mooneyrivlin_c)*log(detF); - } - else { + } else { strainEnergy = mooneyrivlin_10 * trCTildeMinus3 + mooneyrivlin_01 * c2TildeMinus3 + mooneyrivlin_20 * trCTildeMinus3 * trCTildeMinus3 + @@ -110,7 +109,7 @@ public: using std::log; field_type logDetF = log(detF); return strainEnergy + 0.5 * mooneyrivlin_k* logDetF * logDetF; - } else if (mooneyrivlin_energy == "square") { + } else { //mooneyrivlin_energy is "square" field_type detFMinus1 = detF - 1; return strainEnergy + mooneyrivlin_k* detFMinus1 * detFMinus1; }