[Getfem-commits] r5355 - /trunk/getfem/src/getfem_nonlinear_elasticity.cc

[logari81]

Author: logari81
Date: Tue Jul 26 23:37:45 2016
New Revision: 5355

URL: http://svn.gna.org/viewcvs/getfem?rev=5355&view=rev
Log:
fix neo-Hookean constitutive law (cf. Eqs. 3.116, 3.118 and 3.120 in Nonlinear 
Finite Element Methods, Wriggers, 2008)

Modified:
    trunk/getfem/src/getfem_nonlinear_elasticity.cc

Modified: trunk/getfem/src/getfem_nonlinear_elasticity.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_nonlinear_elasticity.cc?rev=5355&r1=5354&r2=5355&view=diff
==============================================================================
--- trunk/getfem/src/getfem_nonlinear_elasticity.cc     (original)
+++ trunk/getfem/src/getfem_nonlinear_elasticity.cc     Tue Jul 26 23:37:45 2016
@@ -643,7 +643,7 @@
                             (lambda/2 * log(ci.i3()) - mu) / ci.i3()), result);
     else
        gmm::add(gmm::scaled(ci.grad_i3(),
-                            lambda/4 - lambda/(4*ci.i3()) - mu / ci.i3()), 
result);
+                            lambda/2 - lambda/(2*ci.i3()) - mu / ci.i3()), 
result);
     
     if (det_trans <= scalar_type(0))
       gmm::add(gmm::scaled(C, 1e200), result);
@@ -672,9 +672,9 @@
       coeff = (lambda + 2 * mu - lambda * logi3) / gmm::sqr(ci.i3());
     } else {
       gmm::copy(gmm::scaled(ci.sym_grad_grad_i3().as_vector(),
-                            lambda / 2  - (lambda / 2 + 2 * mu) / ci.i3()),
+                            lambda  - (lambda + 2 * mu) / ci.i3()),
                 result.as_vector());
-      coeff = (lambda / 2 + 2 * mu) / gmm::sqr(ci.i3());
+      coeff = (lambda + 2 * mu) / gmm::sqr(ci.i3());
     }
 
     const base_matrix &di = ci.grad_i3();