Principes de l' approche modulaire contrast

Après avoir mis en place le cadre et les notations, ce chapitre présente les problèmes posés par les modélisations tridimensionnelles à comportement linéaire : complexité des géométries et des modèles, grand volume de données, encombrement mémoire important et coûts de calcul prohibitifs. Les méthodes de décomposition de domaines, qui essaient de pallier ces problèmes, sont présentées. Elles tentent de découpler le traitement des domaines, permettent une réduction des encombrements et sont fortement orientées vers l’utilisation de calculateurs à architecture parallèle (qui peuvent permettre des gains en temps importants et l’utilisation d’un espace mémoire de très grande taille).

Les principes de l’approche utilisée sont ensuite détaillés. Elle est basée sur une décomposition particulière des assemblages qui donne aux interfaces un rôle majeur. La formulation utilisée pour mettre en place une stratégie de résolution itérative est issue de la méthode LATIN et présente un fort parallélisme. Ce dernier est principalement utilisé, dans ce travail, pour introduire une certaine modularité qui facilite la description et la résolution de gros problèmes. Cette approche a tout d’abord fait l’objet d’une étude de faisabilité pour des problèmes 2D axisymétriques, pour lesquels elle permettait des calculs simplifiés de comportement d’assemblages.

Méthodes de décomposition de domaine

Les méthodes de décomposition de domaine sont développées pour le traitement de problèmes de grande taille (grand nombre de degrés de liberté) pour des comportements linéaires ou localement non-linéaires. Elles sont maintenant utilisées pour une grande variété de problèmes. Avec les possibilités nouvelles apportées par les calculateurs à architecture parallèle, elles ont pris un essor important. En effet, le partionnement du domaine engendré est bien orienté vers le développement d’algorithmes adapté à ces ordinateurs. De plus, il conduit à une réduction de l’encombrement total des matrices de rigidité. Deux grandes classes de méthodes existent : les approches avec recouvrement des domaines entre eux et celles sans recouvrement. Les dernières sont les plus utilisées de nos jours et sont celles auxquelles nous allons nous intéresser. Parmi ces méthodes, trois types d’approches existent : les méthodes primales (ou dites du complément de Schur), les méthodes duales et les méthodes mixtes.

Résolution

Le problème condensé étant de petite taille, il peut être résolu directement [ESCAIG et al, 1994]. Dans ce cas, l’approche est qualifiée de « méthode des sous-structures » et est bien connue en mécanique. Néanmoins, la construction de la matrice du complément de Schur, nécessite l’obtention explicite des matrices [KE]-1 et [KE’]-1. Cette opération n’est pas classique pour un code de calcul et peut s’avérer très coûteuse. De plus, la matrice [S] est très souvent « pleine », contrairement aux matrices utilisées par la méthode des éléments finis. Son traitement est donc particulier. Une idée est donc d’utiliser, pour résoudre le problème condensé, des stratégies itératives qui ne nécessitent pas le calcul explicite de [S], mais uniquement celui de ses produits par des vecteurs. C’est le cas pour les méthodes de type gradient conjugué pour lesquelles les calculs de ces produits sont même complètement parallélisables et dont l’efficacité peut être améliorée par un préconditionnement [De ROECK & LE TALLEC, 1990]. Certains cherchent, néanmoins, à utiliser la méthode directe pour profiter de sa robustesse et de son adaptativité à tous types de problèmes. De plus, elle permet, dans certains cas, de créer des bibliothèques de sous-structures ayant déjà été calculées et donc d’introduire une certaine modularité dans l’utilisation. Pour cela, des méthodes efficaces sont élaborées pour la construction de la matrice du complément de Schur. C’est le cas, par exemple, de la méthode multi-frontale [DUFF, 1986].

Les méthodes mixtes

Dans cette dernière approche, les solutions approchées, obtenues au cours des itérations, ne vérifient a priori ni la continuité en déplacement ni l’équilibre en effort sur les interfaces. C’est le cas de l’approche que nous utilisons pour l’étude d’assemblages de structures tridimensionnelles. D’autres approches mixtes, qui s’identifient à la nôtre dans le cas de liaisons parfaites ou de contact sans frottement entre les sous-domaines, sont proposées. L’une utilise une extension de l’algorithme de Schwarz aux cas de décomposition sans recouvrement [LIONS P.L., 1990] et deux autres utilisent une formulation par lagrangien augmenté [GLOWINSKI & LE TALLEC, 1990 ; LADEVÈZE J., 1985] .

Schéma itératif de résolution

L’approche proposée est basée sur les principes de la méthode LATIN (LArge Time INcrement method) [LADEVEZE, 1996]. Cette méthode, introduite par P. Ladevèze en 1984, est une stratégie générale pour traiter des problèmes non-linéaires d’évolution. Elle n’est pas basée sur les schémas de résolution habituels. Ses nombreuses applications (en élastoplasticité, en viscoplasticité sous chargement cyclique, en dynamique, en grandes transformations, …) ont montré son efficacité et les réductions en coût de calcul qu’elle apporte par rapport aux techniques classiques. La présente étude se place dans une situation dégénérée de la méthode LATIN : le problème d’évolution en temps n’est pas considéré et on ne cherche que la solution du problème à l’instant final. Notre approche n’utilise donc qu’une partie des concepts de la méthode LATIN.

Conclusion

Cette présentation de l’approche proposée pour l’étude d’assemblages tridimensionnels montre la simplicité de mise en œuvre et les premiers avantages engendrés. Les choix effectués pour la discrétisation conduisent à la possibilité d’utiliser un grand nombre d’outils existants pour les méthodes éléments finis classiques. Ainsi, il est naturel et simple de chercher à insérer cette approche dans un code de calcul industriel pour profiter de procédures performantes et optimisées. Ceci fait l’objet du chapitre suivant où l’on présente l’implantation de l’approche CONTRAST dans le code de calcul CASTEM 2000 (développé par le Département des Etudes Mécaniques et Thermiques du Commissariat Français à l’Energie Atomique).

Le rapport de stage ou le pfe est un document d’analyse, de synthèse et d’évaluation de votre apprentissage, c’est pour cela chatpfe.com propose le téléchargement des modèles complet de projet de fin d’étude, rapport de stage, mémoire, pfe, thèse, pour connaître la méthodologie à avoir et savoir comment construire les parties d’un projet de fin d’étude.

Table des matières

INTRODUCTION
I – PRINCIPES DE L’ APPROCHE MODULAIRE CONTRAST
I.1 Problèmes de grande taille
I.1.1 Hypothèses
I.1.2 Notations
I.1.3 Problème
I.1.4 Résolution
I.1.5 Difficultés
I.2 Méthodes de décomposition de domaine
I.2.1 La méthode du complément de Schur
I.2.2 La méthode duale
I.2.3 Les méthodes mixtes
I.3 Principe de l’approche proposée
I.3.1 Décomposition de la structure
I.3.2 Problème mécanique sur les deux entités créées
I.3.3 Schéma itératif de résolution
I.3.4 Discrétisation
I.3.5 Choix des paramètres de direction de recherche
I.3.6 Contrôle du calcul
I.4 Conclusion
II – MISE EN ŒUVRE NUMÉRIQUE ET INFORMATIQUE – PREMIÈRES ILLUSTRATIONS
II.1 Mise en œuvre numérique de la méthode
II.1.1 Implantation dans CASTEM 2000
II.1.2 Organisation des données
II.1.3 Déroulement d’un calcul – Algorithme
II.1.4 Particularités relatives aux interfaces
II.2 Comportement du logiciel – Exemple simple
II.2.1 Comportement de l’algorithme
II.2.2 Influence des paramètres
II.3 Conditions aux limites
II.3.1 Prise en compte des conditions aux limites
II.3.2 Exemple
II.4 Traitement d’un problème à grand nombre de degrés de
liberté
III – TRAITEMENT DES NON LINÉARITÉS DE CONTACT
III.1 Problèmes de contact unilatéral avec et sans frottement
III.1.1 Données et hypothèses du problème
III.1.2 Contact unilatéral sans frottement
III.1.3 Lois de frottement
III.1.4 Problèmes statiques ou monotones quasi-statiques
III.1.5 Approche incrémentale du problème de frottement de Coulomb
III.1.6 Résolution du problème de contact
III.1.7 Analogie avec les problèmes de plasticité
III.1.8 Eléments de contact
III.1.9 Quelques problèmes numériques avec les éléments finis
III.2 Interface de contact unilatéral avec ou sans frottement
III.2.1 Problèmes et hypothèses
III.2.2 Résolution
III.2.3 Interprétation en terme de Lagrangien augmenté
III.3 Validation du traitement des non-linéarités de contact
III.3.1 Sensibilité de l’algorithme au passage adhérenceglissement
III.3.2 Exemple test – Lopin en compression biaxiale
III.4 Non-linéarité de type « joint élastomère unilatéral »
III.4.1 Joint élastomère à comportement non-linéaire unilatéral
III.4.2 Interface « joint élastomère non-linéaire »
III.4.3 Exemple de comportement
III.4.4 Exemple simple
III.5 Bilan sur les modèles d’interface
IV – EXEMPLES INDUSTRIELS D’ ASSEMBLAGES DE STRUCTURES TRIDIMENSIONNELLES – FONCTIONNEMENT ET PERFORMANCES DU LOGICIEL CONTRAST
IV.1 Bride avec vis radiales (comportement)
IV.1.1 Modélisation
IV.1.2 Transmission des efforts par adhérence
IV.1.3 Transmission des efforts par obstacle (cisaillement des vis)
IV.1.4 Comportement itératif particulier
IV.2 Bride avec vis radiales (étude du chargement)
IV.2.1 Modélisation
IV.2.2 Chargement axial
IV.3 Bride conique
IV.3.1 Description du problème.
IV.3.2 Chargement radial.
IV.4 Comparaison avec un code de calcul – Assemblage boulonné
IV.5 Comparaison avec un code de calcul – Bride rivetée
CONCLUSION