Projet

Général

Profil

Actions

Anomalie #385

ouvert

V7.039 : modif cal. compressibilité loi aniso

Ajouté par Frank Petitjean il y a 3 mois. Mis à jour il y a environ 2 mois.

Statut:
En cours
Priorité:
Normal
Assigné à:
Version cible:
Début:
06/08/2025
Echéance:
% réalisé:

0%

Temps estimé:
Temps passé:

Description

Gérard,
Les modifications sur les lois aniso apportées à la version V7.039 ne sont pas sans conséquences sur Omher2025 !
La version 2024 était basée sur la 7.024. Pour la dernière édition j'ai pris la dernière version de Herezh (7.053)et j'ai lancé ma batterie de tests (j'en ai plus de 100 !) ce qui a révélé des problèmes.

Pour modéliser le BSO en axi j'ai recours à une loi ORTHOELA2D_C avec des coef particuliers pour représenter le comportement 1D des rubans mécaniques, rubans qui sont intégrés au film (2D) par une loi des mélanges.

ORTHOELA2D_C  demarre_a_prop_non_nulle_puis_strictement_decroissante_
E1= 5000e6 E2= 1e-2 E3= 1e-2 nu12= 0 nu13= 0 nu23= 0 G12= 1e-7
nom_repere_associe_ repere

Seul le coef E1 est significatif. Cette loi donne les résultats attendus jusqu'à la version 7.037. Avec la version V7.039 il y a une raideur dans le sens 2 et le résultat est faux. J'ai dû réaliser de nombreux calculs entre les versions 7.024 et 7.053 pour isoler le problème.

Est-ce que tu comprends pourquoi les modifications apportées au calcul de la compressibilité dans le cas de la loi aniso ne permet plus d'utiliser mes coefs ? Comment représenter une loi quasi 1D avec la loi ORTHOELA2D_C modifiée 7.039 ?

Merci
Frank


Fichiers

maillage.her (3,43 ko) maillage.her Frank Petitjean, 24/08/2025 11:52
modele.CVisu (21,7 ko) modele.CVisu Frank Petitjean, 24/08/2025 11:52
modele.info (6,36 ko) modele.info Frank Petitjean, 24/08/2025 11:52

Mis à jour par Gérard Rio il y a 2 mois

bonjour Frank,
oui j'ai une idée (cf. #379)
Avec les versions < 7.039 on avait le code suivant pour le calcul de la compressibilité:

{if (log_var_vol > ConstMath::petit)
   {module_compressibilite = untiers * sigBH.Trace() / (log_var_vol);}
else // si la variation de volume est trop faible on passe par la moyenne
// des compressibilités dans les 2 directions d'orthotropie (cf. théorie) {double unsurKs1 = 1./E1 - nu12/E1 -nu13/E1;
double unsurKs2 = -nu12/E1+1./E2-nu23/E2;
double unsurKs3 = -nu13/E1-nu23/E2+1./E3;
module_compressibilite = untiers * untiers * (1./unsurKs1+1./unsurKs2+1./unsurKs3);
};

donc en traction log_var_vol est < 0 et la compressibilité était calculée avec la seconde formule (cf. théorie), c-a-d la formule des petites def

à partir de la version 7.039, le test devient:

if (Dabs(log_var_vol) > ConstMath::petit)

du coup c'est la première formule qui est utilisée.

Donc voilà pour l'explication du comportement.

Maintenant sur le fond:
- Je pense que la différence vient des différences de mesure de def : Almansi - log. A priori on doit avoir des résultats semblables en petites def (proche de 5%) par contre en grandes def c'est différent.

- dans le cas de tes coef, on obtient en traction, versions < 7.039 :
module_compressibilite = 1/9 * E1
ce qui est peut-être assez différent de la formule en log
Du coup quand l'épaisseur est mise à jour, on a également une différence de traitement.

Peut-être que ce serait bien que tu mettes un petit exemple de mise en évidence de ton pb, et voir comment on peut paramétrer correctement ...

Mis à jour par Frank Petitjean il y a 2 mois

Bonjour Gérard,
Je comprends bien en quoi cette modification du calcul de la compressibilité change le résultat mais je ne vois quelles modifications apporter aux paramètres pour retrouver le comportement attendu.
Avec la version 7.037 le ballon se déforme correctement tandis qu'il reste ouvert avec les versions supérieures.
Ci-joint la version minimale (moins de 2 min).
Merci
Frank

Mis à jour par Gérard Rio il y a environ 2 mois

  • Statut changé de Nouveau à En cours
  • Version cible mis à décembre 2025

Bonjour Frank,
En utilisant tes fichiers avec la version j'ai réussi à obtenir effectivement un résultat différent entre les deux versions.
La différence apparait uniquement au 4ième incrément, au niveau du sommet. A priori, la partie au somment rentre en compression d'où un comportement différent au niveau du module de compressibilité.
Il y a énormément de relaxation, même quand j'augmente fortement le lambda, c'est un peu étrange.
Pour l'instant j'en suis au niveau des constatations...

Mis à jour par Frank Petitjean il y a environ 2 mois

Merci d'avoir pris le temps de tester. Le profil axi initial est obtenue par la forme à plat des fuseaux. C'est la présence des talons aux extrémités (sur largeur) qui donne cette forme évasée. Le déplacement forcé des noeuds extrémités permet de retrouver la forme classique avec donc un fort plissement au pôle. Le critère plis est fortement sollicité ! Avec le choix des coef de la loi orthotrope j'ai essayé de faire en sorte que les rubans ne travaillent que dans le sens longi.

Frank

Actions

Formats disponibles : Atom PDF

Redmine Appliance - Powered by TurnKey Linux