Anomalie #379
loi ORTHOELA3D dans LOI_CONTRAINTES_PLANES et plissement membrane
0%
Description
Gérard,
je vais tester la loi ORTHOELA3D dans une LOI_CONTRAINTES_PLANES avec trois buts :
1) avoir un coef de poisson supérieur à 0.5 dans le plan de la membrane et peu de variation d'épaisseur
2) appliquer le critère pli
3) pouvoir la pondérer avec des fonctions nD (ce qui a un lien avec le ticket sur les invariants log https://herezh.irdl.fr/issues/378 )
j'ouvre ce ticket avec pour l'instant pas de demande encore concrète. Il faut que je teste ce qui est déjà possible.
J'ai pour l'instant juste remarqué deux erreurs dans la doc théorique (theorie_herezh++.pdf version 7.037).
Section 10 page 103 (Loi type Hooke 3D initialement orthotrope, puis entrainée) :
=> formule (328) : double signe "-" entre le terme nu21*nu13*nu32 et le terme nu12*nu21
=> formules (329) n'a pas l'air bonne. Devrait être :
1 - 2*nu12*nu23*nu13*E3/E1 - nu12^2 * E2/E1 - nu13^2 * E3/E1 - nu23^2 * E3/E2 > 0
ou bien aussi (celle qui est dans la doc Abaqus) :
1 - 2*nu21*nu32*nu13 - nu12*nu21 - nu13*nu31 - nu23*nu32 > 0
NB : je préfère la première version que celle d'Abaqus parce qu'elle fait apparaitre uniquement les coefs qui sont effectivement mis dans le .info (i.e pas nu21, nu31, nu32). Et de toute façon, la deuxième version serait redondante avec l'équation (328)
Fichiers
Mis à jour par Julien Troufflard il y a un jour
remarque non bloquante :
il y a une faille dans la lecture des parametres "avec_parametres_de_reglage_" à partir de la ligne 566 du fichier Loi_ortho3D_entrainee.cc.
Si on met un seul paramètre et que ce paramètre est juste un mot-clé, il n'est pas pris en compte.
J'ai voulu mettre "seule_deviatorique". Ce qui se passe c'est que :
1) lecture d'un mot-clé ligne 584 (variable "nom" contient "seule_deviatorique"
2) passage dans le if ligne 595 car il y a une fin de ligne eof (variable "nom" prend la valeur "fin_parametres_reglage_" car je n'ai qu'un seul paramètre)
3) ne passe pas dans le if ligne 628
4) retour au while ligne 572
=> fin de boucle
c'est non bloquant car :
on peut activer "seule_deviatorique" si on met un autre paramètre (par exemple "sortie_post_ 0") ou bien un entier bidon après seule_deviatorique (par exemple : "seule_deviatorique 1")
Mis à jour par Julien Troufflard il y a un jour
- Fichier ticket_379.tar ticket_379.tar ajouté
je joins un répertoire de test.
quelques remarques ci-dessous sur les lois proposées dans le .info.
J'ai testé 3 lois (toutes encapsulées dans une LOI_CONTRAINTES_PLANES) :
résultat :MAT_ortho3D : une loi ORTHOELA3D
- tourne très bien (très vite)
- fournit les coefs de poisson prévus nu12 (plan de la membrane) et nu13 (dans l'épaisseur). Dans le fichier gnu, j'ai tracé en déformation log 11-22 et 11-33. Les courbes dévissent car loi en déf Almansi mais on voit bien en petite déf que la loi fonctionne correctement
- REMARQUE : si on met nu12=0.8, le calcul ne passe pas du tout (dès le premier incrément, donc en petite déf). Pourtant, la stabilité de la loi est respectée (si je mets E1=E2=200, nu12=0.8, nu13=nu23=0.2 alors il faut avoir E3 strictement inférieur à 500, d'où le E3=499)
résultat : convergence très lente, j'imagine qu'il doit y avoir beaucoup d'itérations de contrainte plane. A éviter donc, ça ne doit pas être très bon d'associer une partie volumique isotrope avec une loi orthotrope.MAT_ortho3D_additive : une loi ISOHYPERBULK3 ajoutée à une loi ORTHOELA3D seule déviatorique
résultat :MAT_hypo_ortho3D : une loi HYPO_ORTHO3D
- tourne très bien
- donne les bons coef de Poisson sur toute la plage de déformation log imposée
- fonctionne même avec nu12=0.8 au lieu de 0.79
Mis à jour par Julien Troufflard il y a environ 23 heures
test avec critère pli :
- ok en traction sur un élément
- ok aussi sur un gonflage de coussin carré
c'est ok au sens où la convergence est atteinte (avec un critère d'arrêt trs honorable de 5e-4).
mais avec cependant pas mal de messages de non convergence contrainte plane double en cours de calcul
extrait 1 :
LoiContraintesPlanesDouble::Calcul_et_Limitation2_h_b(.. **erreur: la deformation eps'_I 0.56333369) est superieur a 0.5, on ne peut pas s'en servir pour mettre a jour la variation de volume
extrait 2 :
**>> non convergence sur l'algo de la resolution de (sigma.V_2=0, sigma.V_3=0), mail: 1, ele= 30, pti=2 , nb_incr_total=6, nb_iter_total=36, val obtenues: eps_P(2,2)= 0.28840284, eps_P(3,3)= -0.064309088, eps_P(1,2)= -0.38861528 precedentes: 0, 0, 0 dernier_residu: 2249.8129, 2202.4613, -60.413879 LoiContraintesPlanesDouble::Calcul_dsigma_deps (...: valeur final de h/h_0(t+dt)= 1 et de b/b_0(t+dt)= 1