Commit cdd47c98 authored by DIAZ Jerome's avatar DIAZ Jerome Committed by GILLES Sebastien

#1418 Correct factor for the local stiffnesses.

parent 65a8a9b7
......@@ -208,8 +208,8 @@ namespace MoReFEM
double stiffnessloc = weightI4 * density_interpolate_I4 + weightI6 * density_interpolate_I6;
// calculation of the fiber traction and of its derivative with respect ot the square elongation
double raideur_c = stiffnessloc * (long1 - 1.) / long1;
double draideur_c = 0.5 * stiffnessloc / (long1 * long1 * long1);
double raideur_c = 2. * stiffnessloc * (long1 - 1.) / long1;
double draideur_c = 2. * stiffnessloc / (long1 * long1 * long1);
// update of gradient and of second derivative with respect to Cauchy Greene
dW += raideur_c * tauxtau;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment