Anomalie #117
Erreur loi de comportement elastique non linéaire
Description
Bonjour Gérard
Les lois de comportement 1D élastique non linéaire ISO ELAS ESPO1D et ISO ELAS SE1D donnent des résultats faux.
Cela est vrai également avec la loi ISO ELAS ESPO3D encapsulé dans une loi contrainte plane double.
Ci-joint un test élémentaire avec un élément bielle.
PS. Je pense souvent à toi ces derniers jours. Lorsque je lis les infos il n'est question que de Rio :-)
Il y a 2 semaines on a parlé d'un certain Petitjean pour des raisons plus tragiques.
Fichiers
Mis à jour par Gérard Rio il y a plus de 8 ans
- Statut changé de Nouveau à En cours
- % réalisé changé de 0 à 80
A priori cela fonctionne sur ma version (6.766 osX).
On a E=1.5e-9 N/m^2, F= -2.e4 N, S_0 = 1.52e-4m^2, l_0 = 1m
La loi est iso-élastique 1D.
Le coef de Poisson = 0, donc pas de variation de la section: S = S_0
Avec une mesure de déformation logarithmique:
sigma = F/S = F/S_0 = E * log (l/l_0)
d'où on obtient l = exp(2./(15.2*1.5)) = 1.0917
à comparer avec la valeur donnée par Herezh: 1.0914
ce qui doit être correct à la précision demandée près
Avec une mesure de déformation par défaut (Almansi): eps = 0.5(1.- (l_0/l)^2)
on obtient :
l = sqrt(1./(1.-2.*a)) = 1.1013 avec a = 2./(15.2*1.5)
à comparer avec la valeur donnée par Herezh: 1.10199 ce qui est correct à e-3 près.
Il y a peut-être quelque chose que je n'ai pas bien vu ?
Mis à jour par Gérard Rio il y a plus de 8 ans
désolé, j'ai répondu un peu vite, en fait je n'ai fait que vérifier la loi isoélas...
bon... je regarde maintenant les lois non linéaires !
@suivre ...
Mis à jour par Gérard Rio il y a plus de 8 ans
- Statut changé de En cours à Résolu
- % réalisé changé de 80 à 100
Frank,
effectivement les lois non linéaires 1D ne fonctionnaient plus. En cause le fait que les modules de compressibilité et de cisaillement n'étaient pas calculés (oublie regrettable !). Ces modules sont maintenant utilisés systématiquement pour mettre à jour la section.
A priori c'est réparé et j'ai mis à jour les doc d'utilisation et théorique (explication de la méthode de calcul).
J'ai également intégré un paramètre optionnel nu= pour ISO_ELAS_SE1D ce qui permet de tenir compte éventuellement d'une variation de section.
La dernière version 6.766 fonctionne sur l'exemple fourni. Je te laisse l'intégrer dans les cas test automatiques !!
@bientôt
Mis à jour par Frank Petitjean il y a plus de 8 ans
Merci pour cette correction concernant les lois 1D. C'est maintenant fonctionnel chez moi.
Le problème demeure avec la loi ISO_ELAS_ESPO3D encapsulé dans une loi contrainte plane double.
Pour le vérifier, tu peux reprendre mon cas test et ajouter cette loi :
loi_sangle4 LOI_CONTRAINTES_PLANES_DOUBLE
NEWTON_LOCAL avec_parametres_de_reglage
nb_iteration_maxi_ 50
nb_dichotomie_maxi_ 50
tolerance_residu_ 1e-0
tolerance_residu_rel_ 1e-2
fin_parametres_reglage_Algo_Newton_
ISO_ELAS_ESPO3D
E= 1.5e9 nu= 0.
f_coefficient CPL1D DdlP 0. 1. 1. 1. FdlP
type_de_deformation DEFORMATION_LOGARITHMIQUE
fin_loi_contrainte_plane_double_
J'ai le message :_ ERREUR: methode non encore implante !!!!
_
Mis à jour par Gérard Rio il y a plus de 8 ans
Bonjour Frank,
oui, ce n'est pas vraiment un bug mais effectivement la loi Iso_elas_expo3D n'est pas prévue actuellement pour fonctionner avec la loi de contrainte plane et doublement plane. En fait cela vient du fait qu'une méthode particulière "Calcul_dsigma_deps" n'était pas disponible. Après une recherche rapide dans le code, j'ai l'impression que c'est la seule loi pour laquelle la méthode n'était pas finalisée...
Bon... j'ai fait l'implantation et normalement dans la nouvelle version 6.767, cela fonctionne. Par contre je te conseille de ne pas utiliser la mesure logarithmique qui pour l'instant peut poser des pb... dans certains cas de calcul.
Tiens moi au courant.
@+
NB: le mot clé "avec_parametres_de_reglage" n'est pas correct il faut "avec_parametres_de_reglage_"