Ajouté par Julien Troufflard il y a plus de 5 ans. Mis à jour il y a plus de 5 ans.
Description
Gérard,
la loi POLY_HYPER3D ne rend pas un bon résultat. Je t'ai mis un exemple pour lequel on devrait obtenir la même chose que MOONEY_RIVLIN3D (calcul de traction sur un cube). Mais on obtient une contrainte nulle. Le fichier gnu sert pour comparer le résultat du calcul avec ce qu'on obtient pour MOONEY_RIVLIN3D (traction_princ_MOONEY_RIVLIN.maple).
Fichiers
| pb_POLY_HYPER3D.tar (16,4 ko) pb_POLY_HYPER3D.tar | Julien Troufflard, 23/04/2019 15:05 |
Oui effectivement il y avait un bug pas banal: un ligne de switch qui était doublée... pas de message d'erreur, cela fonctionnait mais ... mal !
Ce sera ok à partir de la version 6.894
Maintenant on trouve a priori les mêmes résultats !
NB: j'en ai profité pour mettre en place des sorties de messages disponibles avec un niveau de commentaires > 4 , exemple de mise en données
MAT MOONEY_RIVLIN_3D
C01= 4 C10= 0.4 K= 3000
permet_affichage_ 5
et exemple de sortie:
MooneyRivlin3D::Calcul_DsigmaHH_tdt(...
Hyper_W_gene_3D::InvariantsEtVar1(...
I_B= 3.03023, II_B= 3.06045, III_B= 1.03023
J_r(1)= 3.0003, J_r(2)= 3.00029, J_r(3)= 1.03023
d_J_r_epsBB_HH(1)= Tenseur3HH 0.15492526872916484 -0.079803942488252469 -0.079803942488252469 0 4.4300087203159849e-18 0 , d_J_r_epsBB_HH(2)= Tenseur3HH 0.15339512699304206 -0.079015747353204352 -0.079015747353204352 0 4.3862551009747805e-18 0 , d_J_r_epsBB_HH(3)= Tenseur3HH 7.9999999999999938 8.2417999999999907 8.2417999999999907 0 -4.5751180621777995e-16 0
MooneyRivlin3D::Potentiel_et_var2(...
W_d= 0.0012969801369379752, W_d_J1= 0.40000000000000002, W_d_J2= 4
W_v= 0.32922415640179725, W_v_J3= 21.357349711498202, W_v_J3J3= 665.09910042808224
sigHH Tenseur3HH 168.99935793837369 173.07883772001728 173.07883772001728 0 -9.6078055356062043e-15 0