[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-commits] r5355 - /trunk/getfem/src/getfem_nonlinear_elasticity.c
From: |
logari81 |
Subject: |
[Getfem-commits] r5355 - /trunk/getfem/src/getfem_nonlinear_elasticity.cc |
Date: |
Tue, 26 Jul 2016 21:37:46 -0000 |
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();
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Getfem-commits] r5355 - /trunk/getfem/src/getfem_nonlinear_elasticity.cc,
logari81 <=