Herezh++

À partir de la version V 6.719:

  - Mise en place de la possibilité d'utiliser plusieurs matrices de raideurs et/ou plusieurs méthodes différentes de résolution, dans le cas d'un comportement transitoire non dynamique: algorithme: "non_dynamique" . 

Fonctionnement: Si la résolution du système linéaire avec la matrice normale échoue, la seconde matrice est alors utilisée. Si la seconde échoue, la troisième est utilisée, etc. La doc sera mise à jour dès que possible, mais en attendant voici un exemple de syntaxe qui permet de déclarer une matrice principale de type Lapack symétrique, et une matrice secondaire bande symétrique classique.

para_syteme_lineaire #----------------------------------------------

TYPE_MATRICE BANDE_SYMETRIQUE_LAPACK

MATRICE_S_SECONDAIRE_S_  BANDE_SYMETRIQUE   FIN_TYPE_MAT_SECONDAIRES_

TYPE_RESOLUTION_S_SECONDAIRE_S_  CHOLESKY FIN_TYPE_RESOLUTION_SECONDAIRES_

A voir dans les essais futurs !

---------------------------------

Version 6.720

- mise en place d'une loi d'hystérésis sur un comportement purement volumique. Dans le cas d'un matériau, il est donc nécessaire d'adjoindre à cette loi, une partie déviatorique. En fait la nouvelle loi : Hysteresis_Bulk est prévue pour représenter le comportement de matériaux poreux type mousses - matériaux architecturés. La doc sera mise à jour dès que possible, en attendant, il est possible d'avoir une information via la constitution interactive d'un fichier info et en choisissant comme loi Hysteresis_Bulk, un exemple de loi que j'ai utilisé pour la mise au point suit:

#------------------------------------------

# Nom Materiau    |      Type loi         |

#------------------------------------------

  hyste_bulk         HYSTERESIS_BULK

#-------------------------------------------------------------------

# .......  loi_de_comportement d'hysteresis spherique ........      |

# para de prager(>=0)  :         mu         : limite de plasticite |

#-------------------------------------------------------------------

   np= 2       mu= 8000     Qzero= 100 avec_parametres_de_reglage_

     type_de_resolution_  2 

     cas_kutta_  5  

#     erreurAbsolue_  1.e-3  erreurRelative_ 1.e-5

     erreurAbsolue_  1.e-8  erreurRelative_ 1.e-8

     nbMaxiAppel_  10000

     nb_iteration_maxi_  20  nb_dichotomie_maxi_  1  

     tolerance_residu_ 1.e-8 tolerance_residu_rel_ 1.e-8#5

     depassement_Q0_ 1000

     sortie_post_ 1

#     permet_affichage_ 8 

   fin_parametres_reglage_Hysteresis_

---------------------------------

Version 6.722

Pour l'ensemble des lois hyperélastiques 3D, on introduit la possibilité de sortie de variable de travail. La doc est mise à jour.

Par exemple, il est possible en post-traitement, d'avoir accès à l'intensité du potentiel, sous forme de valeur aux pti, ou sous forme

d'isovaleurs.

--------------------------------

version 6.724

Concernant l'hystérésis 3D: modification des différents algorithmes de gestion des coïncidences et inversion, correction de bugs.

Version 6.725

Lorsqu'il y a une erreur au niveau de la loi de comportement, il y avait affichage de l'erreur et maintenant il y a également affichage du numéro de l'élément et du point d'intégration où est calculée la loi.

-------------------------------

version 6.726

- correction d'un bug sur l'énuméré des coque-poutre et plaque qui faisait que les membranes passaient certaine fois sous les procédures des coques ,poutres et plaques, dans le cas des sorties en post-traitement.

- optimisation des méthodes d'initialisation du post-traitement, appelées lors d'une sortie au fil du calcul => gain de temps sur les sorties au fil du calcul

- introduction de nouvelles grandeurs en sortie concernant le contact => possibilité de visualiser les noeuds en contact, les éléments en contact, les forces de contact aux noeuds esclaves, les forces aux noeuds des facettes maîtres, etc. (cf. les différentes possibilités offertes en interactif après une étude avec contact)

- au niveau des mouvements solides initiaux, possibilités d'effectuer avant le calcul, une opération d'homothétie, avec éventuellement des rapports différents suivant x, y et z  

-------------------------------

version 6.728

- Amélioration de la lecture: normalement maintenant on peut avoir un fichier qui s'arrête à la fin d'une ligne de donnée sans ligne vide qui suit. Sortie dans le .reac final de l'ensemble des réactions des ddl bloqués, sous forme de vecteur (et non de composantes) avec en commentaire une référence associée qui globalise tous les noeuds bloqués. 

- Mise en place de la possibilité d'introduire des chargements sous forme de champ de valeurs scalaires ou vectorielles (suivant les types de chargement). Disponible pour les forces ponctuelles, les pressions, les forces linéiques, suiveuses ou pas. 

- Correction de petits bugs.

Mise à jour de la documentation: ajout des infos pour les lois de type Mooney-Rivlin concernant la mise en donnée de la prise en compte d'un raidissement pour de fortes élongations.

-------------------------------

version 6.729

Modification du contact:

- mise à jour de la documentation au niveau des paramètres de contact, 

- mise en place d'un cas 4 (TYPE_PENALISATION_PENETRATION   4) qui permet d'ajuster la pénalisation pour satisfaire une pénétration maxi fixée. 

- petites améliorations du fonctionnement des différents algos au niveau du calcul automatique du facteur de pénalisation en fonction des raideurs matériaux. 

-------------------------------

version 6.730

Correction d'un bug sur le post-traitement de calcul d'erreur et de remonté au contrainte dans le cas d'un maillage quadrangle.

-------------------------------

version 6.732

- mise en place de la possibilité d'utiliser plusieurs matrices de raideurs dans le cas

  d'un calcul dynamique implicite (comme c'était déjà le cas en statique). L'objectif visé

  est d'utiliser un type de matrice principal qui est performant. Si jamais la résolution 

  du système linéaire global échoue, Herezh switch alors sur un type de matrice secondaire

  éventuellement moins performant mais par exemple plus robuste. Le switch s'effectue sans

  recalculer les raideurs locales ce qui fait que l'impact en temps de calcul est minimisé.

-------------------------------

version 6.733

- amélioration de la prise en compte des conditions linéaires lorsque celles-ci concernent

  soit les ddl d'un noeud, soit le contact solide-déformable par pénalisation.

  En particulier cela concerne les conditions linéaires de type "Déplacement ou 

  positionnement dans un plan (3D) ou sur une droite (2D)". Les plans et droites pouvant

  par ailleurs évoluer pendant le calcul.

  exemple de gain de calcul:

  ancienne méthode sur un maillage triangles de 346 noeuds, pour un calcul d'équilibre 

  (loi élastique sur membrane) en 

  relaxation dynamique avec 146 conditions linéaires: 

  -> environ 44 minutes sur ma machine (2 incréments, 22000 itérations)

  nouvelle méthode: -> environ 5 minutes sur ma machine (2 incréments, 22000 itérations)

  a priori avec les mêmes résultats ... mais reste à confirmer. 

-----------------------------

Mise à jour de la documentation d'utilisation d'Herezh. En particulier:

- au niveau de l'ajout de matrices secondaires (cf. section "Résolution des systèmes linéaires", mots clés MATRICE_S_SECONDAIRE_S_ , FIN_TYPE_MAT_SECONDAIRES_ , TYPE_RESOLUTION_S_SECONDAIRE_S_  , FIN_TYPE_RESOLUTION_SECONDAIRES_)

- au niveau des conséquences des conditions linéaires sur les matrices de raideur et de masse (cf. 1) section "Résolution des systèmes linéaires", remarques en fin de section, 2) section "Dynamique", remarques en fin de section, 3) section "Conséquences des CLL sur le stockage matriciel (largeur de bande)" )

-------------------------------

version 6.734

Introduction de la mesure de certains temps de calcul. On se reportera à la documentation pour plus de détail, car c'est un peu long à décrire.

Un point important est qu'à partir de cette version il est nécessaire d'avoir installé la bibliothèque boost qui est disponible sur toutes les plateformes: linux, window et osx

Dans tous les cas il est préférable d'installer "toute la bibliothèque", car dans les futures versions, l'utilisation de nouveaux éléments de la bibliothèque est prévue en particulier pour le multi-streads

-> les versions linux et osx sont pour l'instant installées dans la partie test de manière à supprimer les bugs de jeunesse.

-------------------------------

version 6.735

Correction d'un bug sur le contact en axisymétrique. La surface associée au calcul du facteur de pénalisation était erronée. Correction d'un bug sur la sortie des forces: en axi, seules les directions x et y ont des composantes non nulles

-------------------------------

version 6.736

Introduction du temps de recherche du contact. Mise à jour de la documentation.

-------------------------------

version 6.737

Introduction de la sortie systématique des temps cpu globaux dans un fichier spécifique: <nom du .info>_temps.cpu

Mise à jour de la documentation

-------------------------------

version 6.738-9

Correction de bug sur la sortie en fichier des temps cpu.

-------------------------------

version 6.740

- Modification de la sortie des incréments (sortie au fil du calcul). Pour les grandeurs autres que celles calculées

  aux points d'intégration, il y avait une sortie pour l'incrément 0. Cette sortie n'était pas cohérente avec les autres, elle a donc été supprimée. Ceci étant, lorsque l'on veut avoir l'état des grandeurs pour t=0 la solution est 

d'effectuer un premier incrément de temps avec un chargement nul.

- Introduction de la possibilité de créer interactivement une condition linéaire de type projection sur plan ou droite,

ou de type général. Cette possibilité était jusqu'à maintenant absente dans le menu de création interactif d'un nouveau fichier de commande .info

-------------------------------

version 6.741

- Versions associées à une gestion complète des dépendances entre les sources -> refonte complète des makefiles. Devrait (?) permettre par la suite de générer rapidement les nouvelles versions Linux (au même titre que celles du système osX)

- correction d'un petit bug sur la génération interactive (via la création d'un nouveau .info) de l'entête correspondant à l'algorithme "information", nécessaire par exemple pour générer de nouvelles références.

-------------------------------

version 6.742

Modification du code pour la sortie interactive de référence: une étonnante différence de fonctionnement entre osX et linux ?? 

-------------------------------

version 6.743

- amélioration de l'algorithme de contact pour le type 4: le facteur multiplicatif (qui dépend de la pénétration maxi imposée) est lissé au cours des incréments, sous forme d'une moyenne pondérée sur 2 pas -> meilleure stabilité 

-------------------------------

V 6.744 

- Dans le cas de l'utilisation d'une UMAT via un dialogue entre Herezh et Herezh ou entre Herezh et Abaqus il faut que le nom de la loi de comportement soit cohérente entre les deux programmes. La version 6.744 intègre un message d'erreur dans le cas où les noms ne sont pas identiques.

-------------------------------

V 6.745

-  dans le cas de l'utilisation d'une loi externe à Herezh via une Umat (cf. doc), extension du cas 3D actuel à la prise en compte du cas des calculs en axisymétrique.  L'extension est également valable lorsque Herezh sert pour calculer la loi de comportement en tant qu'Umat extérieur. En fait, dans le cas d'un dialogue Herezh--Herezh, il n'y avait pas de pb, par contre dans le cas d'un dialogue Herezh--Abaqus, il y avait un pb à cause du fait qu'Abaqus ne fournit que 4 composantes des tenseurs (les composantes non nulles). La version est pour l'instant en test. Le programme c (ou c++) permettant de relier Herezh avec un programme externe et l'exemple d'une liaison avec Abaqus est également accessible sur le site. 

- passage à l'utilisation de la librairie boost version 1.58

-------------------------------

version 6.746

-ajout dans les utilitaires Herezh, de la possibilité de supprimer des éléments à 2 noeuds très proches: au lieu de supprimer les noeuds dans les références, les anciens numéros sont remplacés par les numéros des noeuds restants correspondants

- amélioration de la mise à jour des références de noeuds après le soudage de noeuds très proches

- amélioration de la mise à jour des références d'éléments après la suppression d'éléments superposés: au lieu de supprimer les éléments dans les références, les anciens numéros sont remplacés par les numéros des noeuds restants correspondants

- intégration de la température dans les lois de type UMAT, en interne et en externe.

-------------------------------

version 6.747

- suite à la demande #98, modification de la signification (et du fonctionnement interne) de la  variable température passée en paramètre.

- suite à la demande #94, modification de la signification (et du fonctionnement interne) des variables temps

- mise à jour de la documentation au niveau des chargements: on précise que l'on peut utiliser en option une courbe de charge + échelle + temps mini et maxi

-------------------------------

version 6.748

- ajout de la possibilité de générer automatiquement à la lecture des maillages d'une référence sur les noeuds non référencés par des éléments. cf. la documentation chap: 2.1

- correction bugs sur la prise en compte de la viscosité pour les lois Umat (dépendance au temps)

-------------------------------

version 6.749

- correction de bugs et amélioration de la loi de comportement HYSTERESIS_BULK ... en espérant que tout est ok maintenant !

-------------------------------

version 6.750

- introduction d'une nouvelle classe de fonctions 1D quelconques qui peuvent être définies par une expression littérale. L’implantation utilise la bibliothèque muParser. La mise à jour de la documentation présente les différentes possibilités offertes dans l'écriture des expressions littérales.

-------------------------------

version 6.751

- introduction de la résolution d'un problème thermique via des éléments spécifiques thermiques qui ont pour vocation de pouvoir être couplés faiblement et fortement avec des éléments mécaniques. L'ensemble des classes génériques sont implantées et un premier élément "biellette thermique" fonctionne pour un cas simple de condition de température bloquée en statique.

La partie post-traitement fonctionne également, en format maple et en post traitement graphique: grandeurs disponibles pour l'instant: température, gradient de température, flux thermique, vitesse de gradient thermique. 

bon... maintenant il faut étoffer ! La doc sera mise à jour dès que possible.

-------------------------------

version 6.752

- introduction d'un chargement volumique dit "pseudo-massique" correspondant en fait à un chargement volumique relativement au volume initial.

-------------------------------

version  6.753-54 

- modification et amélioration de la sortie des réactions de contact (mais le travail n'est pas terminé)

- introduction du critère plis sur les biellettes. Ce critère vient compléter ce qui existe déjà avec des membranes ce qui permet de coupler des maillages de membranes et de biellettes tout en garantissant que chaque élément satisfera aux critères plis respectifs.

-------------------------------

versions 6.755-58

- Introduction du pilotage des critères. 

- Introduction de variables globales accessibles de partout, avec un cycle de temps d'accès et/modification rapide. 

- Première utilisation des fonctions multidimensionnelles

- qq modifs sur le calcul des masses virtuelles en relaxation dynamique

-------------------------------

versions 6.759

- Introduction de nouvelles grandeurs  globales accessibles par toutes les classes, et qui peuvent donc être utilisées par la loi critère. Il s'agit pour l'ensemble du système, des grandeurs: 

  . concernant les lois de comportement: l'énergie élastique, l'énergie plastique, l'énergie visqueuse, 

  . concernant les énergies "numériques  : l'énergie visqueuse numérique, l'énergie du bulk viscosity, l'énergie d'hourglass.

-------------------------------

versions 6.760

- correction d'un bug relatif à l'utilisation d'une loi additive à l'intérieur d'une loi critère.

-------------------------------

versions 6.761

- correction d'un bug sur le pilotage de la loi critère

-------------------------------

versions 6.762

- correction d'un bug sur le calcul de valeurs propres sur les tenseurs d'ordre 2

-------------------------------

versions 6.763-4

- introduction des variables globales volume matière et volumes avec des plans. Pour l'instant ce n'est pas testé.

- introduction du critère mixte: relaxation dynamique avec amortissement cinétique puis visqueux (type 4). Fonctionnement par défaut via la précision -> testé ok, -> mise à jour de la documentation

- introduction pour le type 4 de la relaxation dynamique, de l'utilisation d'un pilotage avec une fonction dépendant du temps et/ou une fonction nD dépendant des variables globales. Implantation terminée. Pour l'instant ce n'est pas testé et il faudra mettre à jour la doc.

- introduction pour tous les types de chargement classiques de l'utilisation d'une fonction nD. Implantation terminée. Pour l'instant ce n'est pas testé et il faudra mettre à jour la doc.

- correction de bugs: loi additive, loi contrainte plane.

-------------------------------

versions 6.765

- inversion de l'orientation de l'effort produit par un chargement hydrostatique.

-------------------------------

versions 6.766

- correction du pb de fonctionnement des lois non linéaires 1D ISO_ELAS_SE1D et ISO_ELAS_ESPO1D suite à [Herezh++ - Anomalie #117] 

Introduction dans les lois ISO\_ELAS\_SE1D et ISO\_ELAS\_ESPO1D du calcul du module de compressibilité et de cisaillement, utilisées en particulier pour la mise à jour éventuelle des sections. Vérification sur le cas d'école fourni dans [Herezh++ - Anomalie #117]

-------------------------------

versions 6.767

- intégration de la méthode "Calcul_dsigma_deps" pour la loi Iso_elas_expo3D ce qui permet d'utiliser les routines CP et CP2

-------------------------------

versions 6.768

- demande #115 -> correction d'un bug sur la sortie des grandeurs à la suite de l'utilisation de _suite_point_info_

-------------------------------

versions 6.769

- demande #102: correction d'un bug sur la sortie de la déformation d'épaisseur pour les lois contraintes planes de doublement plane, il manquait une sortie. Cependant ce qui était sortie en non nul, était correct.

-------------------------------

versions  6.770 

- amélioration du restart en particulier pour éléments 1D et 2D. Intégration du recalcul de la métrique à l'instant t. Cela impacte le cas des variations d'épaisseur pour les éléments 2D et les variations de section pour les éléments 1D

-------------------------------

versions  6.771

- introduction de l'opérateur tangent d_sig/d_eps pour l'hypo-élasticité 3D. Mise à jour de la documentation théorique associée. Modification du test de validation en conséquence: on obtient maintenant des résultats quasi identiques au cas 2D contraintes planes lorsque l'on utilise la loi 3D. La différence résiduelle est a priori du au fait que le calcul de l'épaisseur est différent suivant que l'on utilise une loi de contrainte plane hypo ou que l'on utilise une loi 3D avec un critère de contrainte plane (cf. doc théorique sur le calcul de l'épaisseur). 

-------------------------------

versions 6.772-74

- loi additive: introduction d'une fonction nD de pondération fonction de grandeurs locales à une des lois membres

- correction  bug de restart lorsque celui-ci commençait à la fin d'un .BI. Bug étrange, mais une solution a été trouvée.

- correction d'un bug d'initialisation de la  déformation thermique sur loi additive + loi des mélanges

- mise à jour de la documentation utilisateur concernant la loi d'hystérésis bulk, et pour la loi additive, l'utilisation d'une fonction nD de pondération utilisant des grandeurs locales aux lois internes. 

-------------------------------

versions 6.775-77

- cor bug sur matrices secondaires avec contact: les matrices peuvent changer de dimensions

- cor bug init def thermique sur toutes les lois combinées: l'ordre d'initialisation était erronée

- amélioration de l'option _suite_point_info_ pour 2 algos: relaxation dynamique et statique. Il sera possible de modifier également tous les autres algos, mais j'attends qu'il y ait des demandes !

-------------------------------

versions 6.778

- cor bug contact axisymétrique, 

- amélioration du restart dans le cas de contact et de l'utilisation du mot clé _suite_point_info_

- correction d'un bug sur les éléments points pour Umat,

- mise à jour de la doc avec introduction du chapitre Multistep

-------------------------------

versions 6.779-81

- cor bug contact 1D

- amélioration des éléments axi, du contact en axi(cor bug), restart en axi

- amélioration du contact axi avec triangles esclaves quadratiques, intro d'une nouvelle répartition de pti pour trai axi quadracomp, chgt des messages warning et erreur,

- intro loi des mélanges avec pilotage via des grandeurs globales, intro d'une nouvelle mise en données et de nouvelles fonctionnalités (voir doc)

------------------------------

versions 6.782-84

- cor bug lors du passage d'incrément sur la loi CP2

- introduction de bornes mini et maxi pour les variations de h/h0 et de b/b0, ceci pour les lois CP et CP2

- prise en compte des conditions limites cinématiques suivant l'axe 3, par neutralisation, en axisymétrique et en dynamique.

- Mise à jour de la doc : pour la loi des mélanges, chapitre "Les grandeurs", chapitre "Pondérations par fonctions".  

------------------------------

versions 6.785

- Cor bug sur plusieurs bornes TEMPS_MIN= ET TEMPS_MAX= à suivre, dans les conditions limites. 

- Introduction des fonctions nD "FONCTION_COURBE1D" et "FONC_SCAL_COMBINEES_ND". 

- Dans les fonctions nD, amélioration de la prise en compte mixte des variables globales et locales.

- Amélioration de l'utilitaire de création d'un .info: validation par retour chariot des paramètres par défaut à définir avant l'arrivé au menu principal. 

- Ajout de la variable globale "temps_courant".

------------------------------

versions 6.786

- validation avec correction de bug du fonctionnement des fonctions combinées nD

- correction d'un bug sur les chargement avec fonction nD sans courbe de charge.

- correction d'un bug sur les sorties .maple : dans le cas où une variable demandée n'existait pas, Herezh++ dans la version fast, sortait une valeur incontrôlée contrairement à la version lente ou c'était ok. 

------------------------------

versions 6.787

- correction de plusieurs bugs sur la sortie des directions principales (ticket #128)

- mise à jour doc a propos du mot clé les_fonctions_nD qui remplace les_Fonctions_nD dans le .info (ticket #129)

----------

test de notification

------------------------------

versions 6.788-91

- fin cor bug sur la sortie des directions principales

- introduction d’une nouvelle mise à jour et calcul du facteur lambda dans l’algorithme de relaxation dynamique

- introduction de la  vérification automatique des noeuds non référencés par des éléments. Le fait d’avoir des noeuds sans raideur venant d’élément n’est pas interdit, par contre lorsque ces noeuds ne sont pas bloqués cela peut engendrer des erreurs difficiles à voir. Et lorsqu’ils sont bloqués et coïncidant avec d’autres noeuds, on peut être étonné des comportements obtenus.

- cor bug pression  hydro: correction d’une erreur d’orientation (ticket #136)

- cor bug sur les conditions linéaires dans le cas d’une géométrie axisymétrique

------------------------------

versions 6.792-93

- introduction du calcul d'intégrales de volumes (mais non validé pour l'instant)

- correction loi Mooney Rivlin 1D (ticket #140)

------------------------------

version 6.794:

- introduction d'une initialisation de température pour le calcul du temps critique de stabilité, utilisé par les algos de type explicite (ticket #141)

------------------------------

version 6.795:

- introduction d'un nouvel utilitaire permettant la création d'un nouveau maillage à partir d'un maillage déjà existant et d'une référence d'éléments -> le nouveau "sous-maillage" contient également les références réduites.

------------------------------

version  6.796: 

- correction d'un bug sur l'utilisation de variable globale par les fonctions nD existantes, 

- amélioration de l'algorithme d'amortissement cinétique: intro d'un sous-cas qui "semble" bien améliorer la convergence globale, lors d'une variation pendant un même pas de temps, du chargement et/ou de la loi de comportement... à valider plus en détail par retour d'expérience. 

------------------------------

version  6.797: 

- prise en compte automatique des référence i_tout (i= N, F, A, E, G, P) du premier maillage de 

manière à récupérer ces références dans le maillage final,

- introduction d'une nouvelle option : "calcul_geometrique", pour l'algorithme "informations": 

   -> un premier calcul: distance initiale entre deux noeuds

------------------------------

version  6.798-6.801 : 

 - corr bug cas init certains contacts déformables rigides 

 - modification du  pilotage de lambda dans algo relax, amélioration création SFE3: prise en compte de la distorsion du maillage pour un choix automatique des SFE1 dans le cas où les SFE3 sont mis en défaut

 - extension du pilotage des différents chargements à l'aide d'une fonction nD pouvant dépendre de la position du point d'application pour les forces ponctuelles, du point d'intégration de la surface chargée pour les chargements de surface, du point d'intégration de la ligne pour les chargements linéiques 

 - validation sur des cas d'école, des intégrales de volume (ou de surface, ou linéique suivant le type d'élément fini support) de fonction nD (utilisateur), idem pour une intégration en volume + en temps, 

 - introduction des grandeurs intégrées, au niveau global,  au temps courant "et" au dernier pas de temps qui a convergé

 - validation sur des cas d'école, de l'utilisation des grandeurs globales issues des intégrales de volume et en temps, dans des fonctions nD permettant ainsi de piloter tout ce qui est pilotable par des fonctions nD, en particulier les chargements et certaines parties des lois de comportement combiné. 

 - sortie dans le .BI de la valeur des intégrales de volume et en temps

==> (beaucoup de changement de code => demandera peut-être des ajustements à l'usage)

------------------------------

version  6.802 :

 - implantation de 3 nouveaux éléments:

   . élément biel quadratique, 

   . élément biel linéaire axisymétrique,

   . élément biel quadratique axisymétrique,

  Il faut maintenant les tester,

 - modification de la loi Maheo_hyper pour qu'elle fonctionne en dynamique, et en relaxation dynamique

------------------------------

version 6.803

 - modification relaxation dynamique: jusqu'à présent le temps fin de calcul en relaxation dynamique n'était pas respecté exactement (d'une manière analogue à l'explicite). C'est maintenant modifié pour garantir des temps de fin strictement satisfaits. 

 - correction d'un bug sur l'utilisation des fonctions nD dépendantes du point pour les différents chargements.

------------------------------

version 6.804

 - correction d'un bug qui apparaissait uniquement à la fin de l'exécution, lorsque l'on faisait une succession d'opérations sur les maillages avec création et suppressions d'éléments et de maillages. La surcharge de l'opérateur de copie des frontières de l'élément était incomplète d'où un pb au moment de la destruction.

------------------------------

version 6.805 - 6

- modification du statut des données imposées c-a-d lorsque on met en place une condition limite (ou une initialisation) sur une grandeur qui n'est pas directement un ddl primaire (= un déplacement pour un problème mécanique) cette grandeur prend automatiquement le statut de donnée en service (auparavant elle avait le statut de donnée hors service). Ensuite, la donnée est "lisible libre" si elle n'est pas active (le temps est inf au temps mini ou sup au temps maxi de la condition), et bloquée sinon. Auparavant il fallait qu'un objet particulier (exemple : une loi de comportement) active la donnée pour la rendre lisible sinon elle restait HS. A priori ce fct n'est pas suffisamment souple d'où la modif.

Normalement cela ne doit rien changer au fct actuel, mais... à voir

------------------------------

version 6.807-8

- amélioration des éléments biel axi: les éléments biel axi linéaires et quadratiques étaient implantés,

mais n'avaient pas été testés donc plusieurs bugs au rendez-vous qui ont été réparés.

 . les éléments ont été testés sur un cas basique de stabilité d'un ballon sphérique soumis à une pression uniforme, dont on a une solution théorique. On retrouve la solution théorique avec néanmoins une petite déviation au niveau des conditions limites... c'est peut-être normal compte tenu de la discrétisation ?

 . les éléments ont également été testés sur une forme isotensoïde proposée par Anne-Sophie Lectez et on obtient quelque chose qui parait correct ...(à compléter néanmoins) 

- amélioration de l'hystérésis bulk: dans le cas de boucles il peut arriver que sur un même pas de temps on ait une inversion + une coïncidence. Ce phénomène peut apparaître  particulièrement dans le cas de petites boucles et au niveau du début du chargement. Ce cas n'était pas pris en compte d'où des erreurs observées par Émilie Vieville dans ses calculs de thèses. Normalement maintenant c'est ok. 

------------------------------

version 6.809

- nouvelle modification de l'avancement temporel pour l'hystérésis Bulk, avec amélioration du filtrage des dépassements de la saturation

- amélioration du contact: le paramètre PENETRATION_CONTACT_MAXI qui contrôle le maximum de pénétration autorisée dans le cas d'un "TYPE_PENALISATION_PENETRATION     4 " est maintenant satisfait de manière plus rigoureuse.

------------------------------

version 6.810 et 11

- correction bug 1D contact + améliorations commentaires contact

- correction bug sortie des déformations logarithmiques 

------------------------------

version 6.812

- amélioration du contrôle de la loi critère plis: introduction de la possibilité de prendre en compte des paramètres différents pour les lois en contraintes planes (CP) et contraintes doublement planes (CP2)

- modification des initialisations de déformations transversales lors du passage CP à CP2 et vis-versa

- modification de la prise en compte de la non-convergence des algo de Newton en CP et CP2: un appel direct à la loi 3D, est fait avec les variations de déformations transversales inchangées. On suppose qu'à convergence, l'algo de Newton ne sera pas mis en défaut (ce qui semble le cas dans les tests effectués).

------------------------------

version 6.813

- amélioration du calcul des épaisseurs et du post-traitement pour le critère plis

- ajout chapitre 12.3 sur la doc théorique, explication des variations d'épaisseur calculées avec le critère plis

------------------------------

version 6.814-16

- mise en place d'un nouveau fonctionnement des fct nD, via des grandeurs sous forme scalaires, mais également vectorielles/tensorielles. En pratique pour l'utilisateur cela se traduit par la possibilité de piloter les chargements via une fonction nD qui peut dépendre:

 . de grandeurs globales : soit celles existantes systématiquement comme le volume, l'énergie..., soit de nouvelles grandeurs qui apparaissent au cours du calcul sous forme d'intégrales sur un domaine particulier (cf. définition d'intégrale de fonction nD)

 . de grandeurs cinématiques: positions, delta de positions, vitesse et accélérations s'il s'agit d'un calcul dynamique

 . de données interpolées, mises en place par la mise en données: Température par exemple

- Dans le cas des chargements volumiques, il est possible d'utiliser toutes les grandeurs existantes aux points d'intégration. Pour les autres chargements, les grandeurs disponibles sont restreintes à celles existantes aux noeuds.

------------------------------

version 6.817-18

- cor d'un bug sur l'utilisation de fonction nD dans le chargement

- mise en place de la possibilité d'utiliser des fonctions nD conjointement avec des courbes de charges pour appliquer des conditions cinématiques: pour des ddl fixés directement par les fonctions ou pour des mouvements solides. Herezh est également à jour pour créer la mise en données correspondantes en interactif (cf. option -n au lancement d'herezh)

------------------------------

version 6.819: (version test)

- prise en compte pour le calcul des temps cpu, d'une sortie inopinée: les temps cpu était erronés, car le comptage n'était pas correctement finalisé.

- modification de la prise en compte du cas type de contact = 1, dans le dimensionnement de la matrice de masse en relaxation dynamique. Dans le cas où les surfaces maîtres ont une cinématique imposée, les conditions linéaires associées au contact type 1, sont réduites aux ddl de chaque noeud indépendamment . En conséquence, la matrice masse est transformée par défaut en une matrice bande tridiagonale. Mais on peut également indiquer un autre type de stockage via la mise en donnée habituelle. 

-> concernant la demande #115, la comparaison sur les 1000 premiers incréments entre les versions 6.812 (qui comportait un bug qui permettait l'utilisation d'une matrice strictement diagonale) et la version 6.819 avec une bande tridiagonale ne présente pas de différence en temps cpu (la partie résolution du système linéaire est négligeable / aux calculs des forces généralisées)

---------------------------------

Version 6.820: (version test)

- correction bug fonction nD (ticket #156)

---------------------------------

Version 6.821-25: (versions tests)

- corrections de bugs sur les fonctions nD (cf. ticket #156)

- optimisation de l’accès en interne aux grandeurs nécessaires pour la sortie de résultats et pour l'alimentation des fonctions nD (possibilités de bug de jeunesse, car un nombre important de modifs)

- introduction de la possibilité dans les fonctions nD, d'utiliser les coordonnées à t=0 et à t à l'aide des noms: X1_t, X2_t, X3_t, X1_t0, X2_t0, X3_t0

- introduction de la possibilité dans les fonctions nD, d'utiliser les normales à t=0 et à t à l'aide des noms: N_surf_1_t, N_surf_2_t, N_surf_3_t, N_surf_1_t0, N_surf_2_t0, N_surf_3_t0

---------------------------------

Version 6.826-27: (versions tests)

Mise en place d'un restart pour des fonctions nD et des intégrales associées d'où plusieurs ajouts dont:

1) Concernant les intégrations

- le stockage (.BI) au niveau des éléments, des intégrations locales

- le stockage (.BI) des intégrations globales 

- la possibilité de geler une intégration après restart

2) au niveau des fonctions nD:

- un mécanisme récursif de sortie d'erreur lors de l'utilisation de fonctions emboitées, avec une information détaillée sur les paramètres d'appel à chaque niveau,

- une amélioration des tests d'erreurs sur le cas de manipulation de grandeurs très grandes

- la possibilité d'introduire un niveau d'impression particulier, sur chaque fonction de manière séparée: en particulier on peut ainsi afficher la valeur des paramètres d'appel et celles de la fonction

---------------------------------

Version 6.828: (versions tests)

- correction bug #166

- ajout dans la loi Maxwell 3D d'un niveau local de sortie d'info

- validation de l'utilisation d'une umat interne (sans utilisation de pipe) en contrainte plane: la suite via les pipes est en cours.

---------------------------------

Version 6.829-33

- mise en place d'une première version opérationnelle d'Umat CP

- corrections  bug sur algo Newmark, 

- correction d'un bug sur la sortie des .her et .lis lorsque l'on demande une sortie sans calcul (via " plus " mais sans le "avec" après la définition de l'algorithme).  Auparavant sortait des positions nulles !!

- au niveau de la définition de nouvelles références via Herezh: 

    .ajout de la possibilité de choisir des numéros de noeuds au lieu de positions pour la définition des plans, des droites, des points dans le cas des créations de références

    . Correction d'un bug sur l'affichage des frontières lors de la sortie d'informations pour l'utilisateur

    . Mise en place des 3 nouvelles méthodes entre_plan_dist, entre_droite_dist, entre_point_dist, par exemple entre_plan_dist propose une condition entre plans qui est définie à partir d'un plan médian et d'une distance +- via la normale du plan.

- modifications relatives au contact:

   . Mise en place d'une nouvelle gestion des zones de contact plus précise, optimisée --> beaucoup de changements et d'essais d'améliorations (?)

   .prise en compte de contacts multiples

      -> premier test sur le contact collant ok

      ... correction des autres bugs:

   . Mise en place d'un contact collant: mise en oeuvre d'une pénalisation autonome (-> paramètres supplémentaires) sur le frottement, validation sur un cas test

   . Mise en place de la possibilité de suppression des gaps de contact avant le début du calcul en déplaçant les noeuds esclaves sur les surfaces présumées de contact.

   . Mise en place de la possibilité de renuméroter à la volée en fonction des contacts mise en jeux. Application à la possibilité d'optimiser en continu les matrices de raideurs.

---------------------------------

Version 6.834

- ajout d'un message d'erreur supplémentaire dans le cas de la définition de zones de contact qui sont incomplètes ou inutilisables

- ajout dans la loi de Maxwell 3D, dans l'opérateur tangent dsigma/deps , et dans le cas de l'utilisation de fonction d'évolution de E et/ou des mu en fonction de la déformation de mises et/ou de Q_Deps, du cas de repère curviligne (cas normal dans Herezh fonctionnant seul) le cas auparavant implanté était uniquement dédié aux repères particuliers orthonormés : cas d'un fonctionnement d'Herezh en Umat par exemple.

- correction d'un bug (en fait il manquait une partie) sur la définition d'un plan de ref à partir d'un noeud et d'une normale dans la partie : définition interactive de référence

---------------------------------

Version 6.835-6.841

- corr bug Suppression_gap_pour_noeud_collant + d'autres bugs de jeunesse vis-à-vis des nouveaux ajouts

- introduction de la prise en compte de repères d'anisotropie au niveau des éléments 

- travaux concernant la mise en place d'une loi élastique orthotrope entraînée. Pour l'instant, pas encore opérationnelle, la mise en données est lue et exploitée correctement (paramètres fixes de la loi + repère d'orthotropie) ... à suivre.  

---------------------------------

Version 6.842 (version de développement)

- mise à jour de la documentation d'utilisation de la loi d'orthotropie élastique entraînée

- modif dans l'insertion d'un repère d'anisotropie,

- intro de la sortie d'info en post-traitement pour la loi d'orthotropie élastique entraînée: ajout d'une grandeur quelconque type repère

... mise au point de la loi d'orthotropie élastique entraînée: en cours 

- mise à jour des sorties de grandeurs tensoriels,

- introduction d'un niveau de commentaire particulier à la loi critère, mot clé : permet_affichage_

- corr bug sur la sortie en .maple des grandeurs globales créées après le démarrage du calcul par les algos (ex: le compteur d'itération)

---------------------------------

Version 6.843 (version de développement)

- mise à jour et amélioration des sorties de commentaires dans loi critère +  plis, dispo pour la version non fast et non actif pour la version fast

- correction de bugs, et amélioration du fct critère plis 

- amélioration de la prise en compte des masses nulles

- doc utilisateur complété sur paramètre algo relaxation dynamique et présentation des mots clés pour les algos

- doc théorique complété sur des détails

---------------------------------

Version 6.844 (version de développement)

- pour la relaxation dynamique, introduction de la variable globale AMOR_CINET_VISQUEUX qui permet de savoir si on est en amortissement cinétique ou visqueux

- correction d'un bug d'écriture temporaire sur le comptage du temps cpu dans les lois de comp

- introduction dans l'écriture des temps cpu, du type de version : rapide ou non

- doc mise à jour (utilisateur et dans Herezh) sur les paramètres spécifiques d'amortissement mixte en relaxation dynamique

---------------------------------

Version 6.847 

- correction de bugs

- améliorations de la loi maxwell 3D

- ajout de la possibilité de moduler la précision d'équilibre global à l'aide d'une fonction nD qui utilise uniquement des variables globales

- amélioration de la loi critère: l'idée est de figer la direction des plis pendant plusieurs itérations (voir sur un incrément ?)

   . introduction d'un paramètre qui permet d'activer le recalcul de la direction: par défaut le recalcul s'effectue à chaque itération: maintenant il peut-être piloté par une fonction nD

- amélioration des messages d'erreur dans la méthode interne de Newton, utilisée par exemple dans les lois complexes

---------------------------------

Version 6.848 

- mise en place d'une dépendance possible des paramètres d'une loi isotrope élastique, à une fonction nD

  . possibilité d'obtenir en sortie les valeurs des paramètres matériau réellement calculés à chaque point d'intégration

- mise en place pour la loi de Maxwell 3D de la possibilité d'obtenir en sortie les valeurs des paramètres matériau (E, nu, mu, mu_p) réellement calculés à chaque point d'intégration

---------------------------------

Version 6.849

- premières validations de la loi d'orthotropie entraînée:

  1) cas d'un cube et d'un hexaèdre dont les cotés sont alignés avec les directions initiales d'orthotropie:

    . vérification de la convergence en implicite via l'opérateur tangent d_sigma/d_ddl,

    . vérification du bon calcul de chacune des contraintes (tractions, cisaillements)

  2) cas d'un hexaèdre avec plusieurs éléments et des directions initiales d'orthotropie, à 45° des côtés de l'hexaèdre:

    . vérification du calcul du repère via les fonctions nD et de son transfert dans la loi de comportement aux points d'intégration,

    . vérification de la bonne convergence du calcul en implicite

  3) cas d'une poutre en flexion maillée en éléments quadratiques avec intégration complète, et un repère d'orthotropie à 45°

    . vérification de  la bonne convergence du calcul en implicite

  NB: pour l'instant l'opérateur tangent d_sigma/d_eps n'est pas encore opérationnel, la mise au point de l'opérateur d_sigma/d_ddl a mis en évidence des modifications a effectuer qui sont en cours.   

---------------------------------

Version 6.850

- mise en place de l'opérateur tangent d_sig/d_eps pour la loi de comportement d'orthotropie entraînée -> possibilité d'utiliser un comportement ce contrainte plane.

 . mise à jour de la Documentation:

  .. partie théorique : Loi de Hooke initialement ortotrope puis matériellement entraînée.

     theorie_herezh++.pdf : version 6.850: chapitre 4: page 35 - 46

  .. manuel d’utilisation: 

     utilisation_herezh++.pdf : version 6.849: chapitre 36.15.1: page 290 - 292

- Tests effectués: en 3D volumique, en 2D contraintes planes

  . sollicitations basiques : traction, cisaillement: comparaison analytique,

  . convergence sur un problème plus complexe: 

     . flexion d’une poutre avec des éléments complets quadratiques, algorithme Newton-Raphson en quasi-statique

     . gonflement d’une membrane via une géométrie de coussin carré (sans la prise en compte de plis): algorithme de Relaxation dynamique

- validation de la version 6.850 vis-à-vis des tests actuellement enregistrés pour Herezh

---------------------------------

Version 6.851

- modif de l'algo utilitaire de renumérotation lorsqu'il y a plusieurs maillages avec éventuellement des CLL: cf. #175 (il y a encore de nouvelles modifs prévues pour automatiser les choix par défaut)

- modif de la sortie des contraintes et déformations locales dans le cas d'une loi orthotrope entraînée (les contraintes et déformations exprimées dans le repère d'orthotropie)

- cor d'un bug sur les sorties en format maple des moyennes de grandeurs quelconques et particulières

---------------------------------

Version 6.852-53 (version de développement)

- complément dans la loi orthotrope entraînée concernant le transport des grandeurs associées à chaque point d'intégration

- ajout de la prise en compte des raideurs de contact dans le calcul de la stabilité pour l'algorithme de relaxation dynamique

- mise en place de plusieurs grandeurs disponibles en post-traitement du contact, accessibles avec l'option debug et au fil du calcul (force de contact, pénétration, glissement, force de glissement, énergies, statut ... voir doc chap: 50.6 Remarques sur le post-traitement )

- modification exploratoire du type 8 de contact... à suivre

---------------------------------

Version 6.854 (version de développement)

- amélioration de l'algorithme de calcul d'erreur et sur la sortie des résultats d'erreur 

- correction d'un bug sur le dimensionnement du conteneur des volumes déplacés en 2D et 1D, calcul qui n'existe pas en 2D et 1D.

---------------------------------

Version 6.855 (version de développement)

- correction erreur sur la sortie des vecteurs propres pour les tenseurs (il y avait inversion des directions)

- introduction d'un addimentionnement au niveau du calcul des inverses de matrices de dimension 1,2, 3

- modification de la loi ortho vis-à-vis de son utilisation en CP2 + plis

---------------------------------

Version 6.856 (version de développement)

- correction erreur sur l'utilisation de force linéique avec des éléments pentaédriques

- correction erreur sur la mise en données d'un calcul 1D

---------------------------------

Version 6.857 (version de développement)

- correction d'un bug un peu gênant sur le constructeur de copie de noeuds: qui induisait certaine fois un crash lors de la fusion de maillages

---------------------------------

Version 6.858-6.864 (versions de développement)

- corr d'un bug de vérification

- amélioration loi CP2

- modif recalcul dir plis

- ajout de nouvelle sortie (vecteur résidu avant résolution)

- amélioration sortie grandeurs particulières loi plis, et de l'affichage d'infos

- modif affichage incrément sur l'algo global pour loi Umat

- ajout de la prise en compte de la température au noeud pour un élément point (utilisé pour la loi Umat)

- correction d'un bug sur le constructeur de copie de l'hystérésis 3D (utilisé pour un critère plis)

---------------------------------

Version 6.865

- introduction dans l'algorithme de relaxation dynamique du paramètre "et_recalcul_masse_a_la_transition_" qui pour un amortissement mixte (type 4) permet de recalculer ou non, la matrice masse juste au moment de la transition: Hugo Le Meitour a montré que le recalcul de la masse améliore la convergence (pour certains cas, sans cette option, pas de convergence)

- introduction dans l'algorithme de relaxation dynamique du paramètre "opt_visqueux_recal_mass=" qui permet d'avoir conjointement deux types de re-calcul de masse, un type pour l'amortissement visqueux et un type pour l'amortissement cinétique

---------------------------------

Version 6.866-69

- correction de plusieurs bugs d'initialisation

  . initialisation des algos utilitaires, lorsqu'il y a des fonctions nD

  . initialisation de tous les algos, dans le cas de l'utilisation d'une loi ortho entraînées avec une en particulier une sortie d'info (la création du CVisu plantait !)

- introduction d'un message explicite d'erreur en cas de taille de matrice de raideur trop importante pour l'ordinateur utilisé, ceci dans le cas d'un stockage matriciel par défaut (stockage bande symétrique)

- introduction de la possibilité de piloter le paramètre inter_nb_entre_relax (cf. amortissement cinétique) via une fonction nD

- introduction de la possibilité de piloter le paramètre de recalcul de la masse via une fonction nD (voir la documentation)