1. Introduction#

L’opérateur op0070 est constitué de plusieurs parties:

  1. Lecture des données;

  2. Initialisation des données;

  3. Construction et résolution d’un système linéaire;

  4. Mise-à-jour des champs;

  5. Archivage des résultats, post-traitements.

Le tout étant encapsulé dans trois niveaux de boucle: pas de temps, boucles de point fixe (pour le contact) et itérations de Newton avec une gestion de type événementielle.

Ce document propose la description de ces différentes zones de l’opérateur, en insistant particulièrement sur le côté structurant des SD (§ 2 ), des événements (§ 3.3 ) et sur l’algorithme général, boucle par boucle (§ 4 ).

À noter que pour aider au débogage, il est possible de sauvegarder les champs de déplacements dans un fichier MED, à chaque itération de Newton. Pour cela, il faut surcharger la routine dbgcha en changeant le booléen dbg=.false. en dbg=.true. Les champs de déplacements seront alors sauvegardés dans un fichier au format MED, sur l’unité logique 80 (ne pas oublier d’ajouter un fichier sur cette unité en sortie dans astk).