There are 3 problems, for now i am working on the first one that is to say the deformation of a 3D rectangular beam clamped on one side and with a pressure applied to the bottom face.
While the deformed configuration given by Getfem is relatively close to the reference(s) solution(s) provided by the benchmark, a visible difference between them still remains and i don't understand where it comes from.
The material is governed by a transversely isotropic constitutive law with an incompressibility constraint, often used in cardiac modelling, where the strain energy function is a function of the components of the Greenâ€“Lagrange strain tensor E.
I tried 2 differents implementations of this law:
-the first use the symmetry of the Green-Lagrange strain tensor to simplify the strain energy function
-The second does not (ergo it is necessary to write the 9 components of the second piola Kirchhoff stress tensor and the 81 components of the fourth order tensor)
Please find enclosed
-the comparison in the first case: Results.png
-the comparison in the second case: Results_nosym.png (slightly better results but 15 times as slow as the first version)
-the python program used to compute the derivative and second derivative of the strain energy function in the first case.
-the implementation of the laws in getfem_nonlinear_elasticity.cc and getfem_nonlinear_elasticity.h
-The program Guccione.cc and Guccione.param used to produce these very pictures
in both pictures, the reference solution is in grey.
The computation uses Q2/Q1 elements (displacement/lagrange multiplier), since there is no restrictions regarding these aspects.
I have tried with a quasi-incompressibility condition instead of the Lagrangian multiplier: same result (which was to be expected).
I have also tried with other meshes (more or less refined) used by other teams but in vain.
Could someone have a look and provide some advices regarding this case/tell me what i am doing wrong?
Thanks a lot.
Yours sincerely,
David.