Herezh++

Je réponds à Julien,

J'ai bien conscience de la brutalité de la chose et en ce moment un calcul tourne avec, comme tu le proposes une loi simple pour appliquer les CLL+plis puis bascule en HH.

Avec le .info pour correspondant à la 1ere proposition

Frank

J'ai lancé un calcul avec un premier équilibre avec une loi HEo + plis et CLL. A l'incrément suivant j'active la loi HH par une loi des mélanges. Le calcul plante à la 1ere itération :

erreur de loi de comportement element= 1 point d'integration= 1

Pourtant à l'incrément convergé précédent (HEo+plis+CLL) les résultats sont correctes et identiques au cas HEo seul.

Mais à quoi bon cette stratégie puisque avec une loi des mélanges les 2 lois sont simultanément calculées donc la loi HH mémorise bien l'histoire du matériaux même si le choix du alpha ne la rend pas active dans la contrainte combinée.

Je ne vois donc pas de "stratégie" pour faire tourner la loi HH dans le cas de BSO !

Je joins la mise en données, seul le .info est différent du fichier précédent.

Frank,

dans ta loi de mélange, il faut mettre l'un mot-clé suivant après chaque loi de comportement du mélange :
demarre_a_prop_non_nulle_puis_strictement_decroissante_
calcule_si_prop_non_nulle_

L'un ou l'autre empêchera de calculer toutes les lois dans le mélange.

ci-joint un exemple de loi de comportement. Choisis l'un ou l'autre (j'ai mis demarre_a_prop_non_nulle_puis_strictement_decroissante_ et ensuite en commentaire calcule_si_prop_non_nulle_).

Dans ton cas, si le mélange est gouverné par le temps, alors il n'y a pas de risque que cette grandeur décroit.
Donc je pense => calcule_si_prop_non_nulle_ suffira

si c'est le critère d'arrêt par exemple qui gouverne le mélange, alors il vaut mieux mettre demarre_a_prop_non_nulle_puis_strictement_decroissante_

Frank,

je suis en train de tester ça (cf pièces jointes) :
temps 1 => deploiement en isoelas sans pli (critere arret 1e-3)
temps 2 => application des CLL en isoelas sans pli (critere arret 1e-3)
temps 3 => application loi HE sans pli (critere arrêt à 1e-4)
temps 4 => applciation loi HH + plis (critere arrêt à 1e-4)
temps 5 => tous les noeuds bloques

à noter qu'au temps 4, la loi HH avec pli a une fonction f_critere qui dépend de la norme de convergence. C'est une faille dans mon calcul car cette norme va osciller. d'où la demande suivante à Gérard :

-------------------------------------
POUR GERARD :
-------------------------------------
ça serait vraiment utile d'avoir une grandeur de nom "norme_de_convergence_mini" (norme de convergence minimum atteinte au cours des itérations). Le but serait de toujours utiliser la meilleure convergence atteinte dans une fonction 1D ou nD. Par exemple, ça permettrait de piloter un avecNiveauSigmaI_mini_pour_plis_ dans la loi pli sans revenir en arrière.

voilà, à voir si mon calcul converge... et si oui quand? car mon ordi n'est pas très puissant.

erreur dans test1.info précédent (la loi HH n'était jamais activée)

ci-joint une version corrigée de test1.info et un fichier test2.info qui est le même calcul mais avec une fonction f_critere échelonnée différemment.

a priori, c'est parti pour calculer toute la nuit. Les itérations en loi HH s'égrenne au rythme de 1 itération toutes les 2 secondes. ça plante pas pour l'instant mais il va falloir être patient.
à suivre...

j'obtiens des erreurs avec les 2 fichiers précédents.

Finalement, c'est le fichier pièce jointe (test4.info) qui a l'air d'aller le plus loin. J'ai mis une fonction f_critere constante et égale à -1e-4. J'ai changé les paramètres d'hystérésis :
nb_iteration_maxi_ => 30
nb_dichotomie_maxi_ => 10
nb_maxInvCoinSurUnPas_ => -20

si tu veux le tester chez toi. J'ai été obligé de l'arrêter, il avait atteint 0.05 en précision.

je n'ai pas de raison particulière à ce test. Je me suis dit qu'il pouvait y avoir pas mal d'inversion avec la loi pli. Mais je n'en sais rien. Difficile de savoir quoi faire quand on ne sait pas d'où vient l'erreur. Je sais juste que l'erreur vient de l'élément 1 qui est un des éléments d'un pole.

Bonjour à tous,
quelques commentaires:
- l'idée de julien d'appliquer les CLL sur une loi simple tout d'abord puis HE avec plis puis loi finale, me paraît effectivement astucieuse, car les CLL sont actuellement appliquées de manière exacte (= des multiplicateurs de Lagrange condensés) du coup c'est tout sauf progressif, donc effectivement avec une loi HH la convergence peut-être difficile.

- les mots clé demarre_a_prop_non_nulle_puis_strictement_decroissante_ et calcule_si_prop_non_nulle_ ont effectivement été mis en place pour éviter de calculer systématiquement toutes les lois.

- "ça serait vraiment utile d'avoir une grandeur de nom "norme_de_convergence_mini" (norme de convergence minimum atteinte au cours des itérations). Le but serait de toujours utiliser la meilleure convergence atteinte dans une fonction 1D ou nD. Par exemple, ça permettrait de piloter un avecNiveauSigmaI_mini_pour_plis_ dans la loi pli sans revenir en arrière."
J'ai en prévision de pouvoir stocker au niveau des fonction nD, des grandeurs que l'utilisateur choisis. Je ne sais pas trop encore dans le détail comment cela fonctionnera mais l'idée serait:
. l'utilisateur définit en plus de la fonction initiale, une seconde qui donne un booleen : si valide la valeur calculée par la fonction est stocké, sinon rien.
. La clé permettant d'accéder ensuite à la valeur stockée est également définit par l'utilisateur via une variable dont la valeur au moment du stockage, sert de clé : si la valeur de la clé est identique à une valeur déjà stockée, l'ancienne valeur est écrasée
. le stockage s'effectuerait via une table de hachage ou ou une map C++
Dans le contexte de cette nouvelle possibilité des fonctions nD, l'utilisateur pourrait, je pense, lui-même définir la norme_de convergence mini. En fait il pourrait définir n'importe quelle grandeur mini...via une fonction ad hoc

Je vais aussi regarder les tests...
@suivre !

Suit à la remarque justifiée de Julien, j'ai introduit dans la loi isoelas2D_C une sortie conditionnelle d'infos à l'écran via le paramètre permet_affichage_ de manière à avoir plus d'infos sur les grandeurs calculées.
Cette option n'est dispo qu'avec la version non fast (pour limiter les tests).

J'ai repris le test4 et sauvegardé le résultat à t=1, puis fait un restart à 1 pour aller à 2s, avec comme paramètre de la loi élastique les mêmes que ceux initiales + permet_affichage_ ->:

#-- loi 1 melange principal (isoelas)
    ISOELAS2D_C calcule_si_prop_non_nulle_
        800.e6  0.45
        permet_affichage_  5

et voici un extrait de ce que j'obtiens:

1) cas sans pb:
...
--- loi de comportement iso elas 2D CP --- mail: 1 ele: 23 npti: 1
module_compressibilite= 2.66667e+09 module cisaillement= 2.75862e+08 eps33= -0.00206603 ener_elas= 3199.54
Ieps 2D = 0.00252514 Isig= 3.67293e+06
...

2) cas avec erreur :
...
--- loi de comportement iso elas 2D CP --- mail: 1 ele: 13 npti: 3
module_compressibilite= 2.66667e+09 module cisaillement= 2.75862e+08 eps33= 10.1419 ener_elas= 7.70862e+10
Ieps 2D = 12.3957 Isig= -1.80301e+10
-- loi de comportement iso elas 2D CP --- mail: 1 ele: 13 npti: 4
module_compressibilite= 2.66667e+09 module cisaillement= 2.75862e+08 eps33= 39.4869 ener_elas= 1.16832e+12
Ieps 2D = -48.2617 Isig= -7.01989e+10

• pb dans le calcul de la nouvelle epaisseur: le coef de compressibilite moyen = 2.66667e+09 est par rapport a delta la trace de sigma = -4.39271e+10 10 fois inferieur !! ce n'est pas normal a priori on continue sans changer l'epaisseur ...
  QuadraMemb::CalEpaisseurMoyenne_et_vol_pti(...)
  ...

L'élément 13 est tout près du bord inférieur.

On obtiens une def mécanique d'épaisseur de 39.5 !! là on est complètement en dehors des clous de la loi isoelas... d'où le pb

Donc l'introduction des CLL de manière exacte introduit un choc numérique vraiment brutal !!
Actuellement, comme le message d'erreur l'indique, quand un comportement anormal est détecté, l'épaisseur de la membrane n'est pas remise à jour. Mais il y a tous les cas où l'erreur n'est pas forcément très visible et qui ne doivent pas être très bons pour la qualité du calcul...
En élasticité, la convergence arrive à ce faire, par contre en HH, cela doit entraîner pas mal de difficultés.

bon... je continue à réfléchir

NB: Julien je mets la version 6.918 en osX sur le serveur si tu veux t'en servir