Remerciez-le!

Remerciez @Admin pour avoir partagé cet document gratuitement, de la manière la plus simple, en partageant sur les réseaux sociaux.

République Algérienne Démocratique et Populaire Ministère de l’Enseignement Sup

République Algérienne Démocratique et Populaire Ministère de l’Enseignement Supérieur et de la Recherche Scientifique Université de Batna Faculté des Sciences de l’Ingénieur Thèse Préparée au Département de Génie Civil Présentée par LAHBARI Noureddine Pour l’obtention du titre de Docteur d’Etat Spécialité Mécanique des Structures Etude de la Stabilité des Plaques et Coques Métalliques par la Méthode des Eléments finis Soutenue le 15 décembre 2007 Devant le jury composé de : Dr HOUARI Hacene Professeur U. Constantine Président Dr ZEDIRA Hamma Maître de Conférence C.U. Khenchela Rapporteur Dr KARECH Toufik Maître de Conférence U. Batna Examinateur Dr CHEBILI Rachid Maître de Conférence U. Biskra Examinateur Dr MAZOUZ Hamoudi Maître de Conférence U. Batna Examinateur TABLE DES MATIERES Résumé Abstract Introduction 01 Chapitre 1 Etude bibliographique 1-1 Bibliographie des éléments finis de membrane 04 1-2 Bibliographie des éléments finis de plaque 05 1-3 Bibliographie des éléments finis de coque 11 Chapitre 2 Formalisme conceptuel de la stabilité 2-1 Stabilité des plaques 15 2-1-1 Introduction 15 2-1-2 Système de forces dans une plaque rectangulaire 15 2-1-3 Equations non linéaires de stabilité 18 2-1-3-1 Relations cinématiques 18 2-1-3-2 Relations constitutives 23 2-1-3-3 Equations d’équilibre 24 2-1-4 Equations de la stabilité linéaire 28 2-1-5 Application de la méthode énergétique à l’analyse du flambement de plaque 30 2-1-6 Equations de l’élément fini des plaques minces 32 2-1-7 Analyse de la valeur propre par élément fini 37 2-2 Stabilité des coques 38 2-2-1 Introduction 38 2-2-2 Géométrie d’une coque mince 39 2-2-3 Système de forces internes dans une coque mince 41 2-2-4 Comportement structural des coques 44 2-2-5 Aperçu sur le flambement de coque 45 2-2-6 Théories de flambement des coques 47 2-2-7 Analyse de la stabilité des coques cylindriques 48 2-2-7-1 Théorie non linéaire de la stabilité 50 2-2-7-1-1 Relations cinématiques 50 2-2-7-1-2 Relations constitutives 53 2-2-7-1-3 Equations d’équilibre 54 2-2-7-2 Théorie de la stabilité linéaire 58 2-2-7-3 Analyse de la stabilité des coques par élément fini 61 Chapitre 3 Méthodologies Numériques 3-1 Définition des paramètres à l’étude 66 3-1-1 Dimensions géométriques 66 3-1-2 Types de chargement 68 3-2 Outils numériques 70 3-2-1 Techniques de calculs éléments finis 70 3-2-2 Code ABAQUS 71 Chapitre 4 Expérimentations numériques 4-1 Plaques 83 4-1-1 Compression uniaxiale 83 4-1-1-1 Deux côtés simplement appuyés 83 4-1-1-2 Un côté simplement appuyé l’autre libre 102 4-1-1-3 Deux côtés encastrés 113 4-1-1-4 Un côté encastré l’autre libre 120 4-1-2 Compression biaxiale 129 4-1-3 Cisaillement 136 4-1-3-1 Tous les côtes simplement appuyés 136 4-1-3-2 Deux côtés encastrés et deux côtés simplement appuyés 149 4-1-3-3 Tous les côtes encastrés 157 4-2 Coques 166 4-2-1 Compression axiale uniforme 166 4-2-1-1 Les deux extrémités simplement appuyées 166 4-2-1-2 Une extrémité encastrée, l'autre simplement appuyée 169 4-2-1-3 Une extrémité encastrée, l'autre libre 171 4-2-2 Compression axiale avec pression extérieure 173 4-2-2-1 Les deux extrémités simplement appuyées 173 4-2-2-2 Une extrémité encastrée, l'autre simplement appuyée 175 4-2-2-3 Une extrémité encastrée, l'autre libre 177 Conclusions, recommandations et perspectives 180 Bibliographie 182 Chapitre 1 : Etude bibliographique 3 Chapitre 1 Etude bibliographique 1-1 Bibliographie des éléments finis de membrane 1-2 Bibliographie des éléments finis de plaque 1-3 Bibliographie des éléments finis de coque Chapitre 1 : Etude bibliographique 4 Etude bibliographique 1.1. Bibliographie des éléments finis des membranes Le premier élément de membrane quadrilatéral isoparamétrique (Q4) a été développé par Taig et Kerr [61].Cet élément de type déplacement est compatible,mais il est trop rigide surtout dans le cas de flexion dominante à cause de la mauvaise présentation de la déformation de cisaillement dans le plan de l’élément . Depuis 1965,beaucoup d’efforts ont été faits pour améliorer l’élément de Taig en utilisant des principes variationnels variés ou en utilisant des champs de déplacements incompatibles.Une contribution importante a été faite par Wilson et al [74] qui ont introduit deux modes incompatibles supérieurs pour améliorer l’élément Q4, mais cet élément ne passe les patch-tests de déformations constantes que dans le cas de maillages réguliers. Taylor et al [165] ont proposé une modification de l’élément de Wilson afin de passer les patch-tests quand le maillage est distordu. L’avantage d’éléments incompatibles mais convergents a alors été prouvé. Le premier élément de type mixte hybride a été proposé par Pian [138].Bergan et Hansen [38] ont présenté l’approche appelée « Individual Element Test » (IET) où les conditions de patch- tests décrites par Irons et Razzaque [102] ont été utilisés pour une formulation directe de la matrice de rigidité. Bergan et Nygard [39] ont étendu l’approche IET qui a conduit à la formulation libre (FF) (Free Formulation ). Park et Stanley [136] ont décrit et évalué la formulation naturelle (ANS) (Assumed Natural Eements ). Tang et al [164] ont formulé des éléments « quasi conformes ». Pian et al [139] ont présenté une nouvelle formulation d’élément de type mixte hybride. Simo et al [155]ont proposé l’élément EAS (Enhanced Assumed Strain Element ) basé sur un principe variationnel généralisé. Chen et Cheung [56] ont présenté un autre principe variationnel mixte basé sur la théorie de Pian et Sumihara [139] et proposé l’élément de type hybride en contraintes où les approximations des contraintes et des déformations sont choisies telles qu’une condition d’orthogonalité soit satisfaite. Par cette approche, les éléments classiques de type hybride en contraintes sont simplifiés. Zhao et al [177] ont donné une autre formulation généralisée de l’élément Q4 et l’équivalence avec un élément de type hybride en contraintes a été prouvée sous certaines conditions. Très récemment, Zhong et Zeng [180] ont présenté un autre élément quadrilatéral appelé « Rational Element ». Dans la formulation, les champs de déplacements u, v sont décrits par une combinaison de huit solutions de base pour le problème d’élasticité plane dans un élément. L’approximation des déplacements est incompatible, mais avec des modifications cet Chapitre 1 : Etude bibliographique 5 élément peut passer les patch-tests de contraintes constantes et une précision en contraintes et en déplacements supérieure à celle de l’élément Q4 de base peut être observée. Maus [124] a présenté aussi trois types d’éléments quadrilatéraux où les mêmes champs de déplacements sont utilisés pour obtenir les vecteurs des forces et les vecteurs des déplacements correspondant aux variables nodales. La matrice de rigidité élémentaire est établie directement par la multiplication des vecteurs des forces et des déplacements. Une procédure assurant la symétrie de la matrice de rigidité est décrite dans l’article. Ces éléments ont une formulation un peu compliquée. Dans l’article historique sur les éléments finis, Turner et al [170] ont décrit un élément rectangulaire incompatible où un champ de déplacements quadratiques est obtenu à partir d’un champ de contraintes linéaires à 5 paramètres. Les équations d’équilibre et l’équation de compatibilité pour le problème d’élasticité en contraintes planes sont satisfaites dans un élément qui peut ainsi être vu comme un élément hybride en contraintes à déplacements incompatibles. Les éléments RQ4Z et RQ4B présentés par Zheng et Batoz [178] sont basés aussi sur un champ de contraintes à 5 paramètres mais ils sont généralisés au cas de l’élément quadrilatéral et des modifications sont introduites pour satisfaire les patch-tests. Les deux éléments RQ4Z et RQ4B ne sont pas indépendants du choix des axes des coordonnées. Pour surmonter cette imperfection, deux nouveaux éléments quadrilatéraux à 5 nœuds RQ5Z et RQ5B ont été formulés et testés. Ces deux éléments sont basés sur un champ de déplacements quadratiques complets et satisfont au critère d’invariance par rapport au système de coordonnées. De meilleurs résultats que ceux obtenus avec les éléments RQ4Z et RQ4B sont toujours observés. Pour généraliser cette famille d’éléments, deux nouveaux éléments ont été établis de la même manière : un élément pentagonal à 5 nœuds et 10 degrés de liberté RP5 et un élément hexagonal à 6 nœuds et 12 degrés de liberté RH6. 1.2. Bibliographie des éléments finis de plaque Un très grand nombre d’éléments finis de plaque en flexion a été développé. Les éléments finis de plaque peuvent se classer en trois types : éléments de type Kirchhoff, éléments de type Reissner/Mindlin, éléments de type Kirchhoff discret. A cause de l’exigence de continuité C ,les éléments de plaque de type Kirchhoff sont relativement difficiles à formuler. Le premier élément de plaque triangulaire de type Kirchhoff à un déplacement w et deux rotations aux nœuds sommets a été proposé par Zienkiewicz [182]. Mais cet élément ne peut pas passer les patch-tests parceque la continuité de la rotation normale n’est pas satisfaite. Des modifications par Bergan [41] et Specht [157] ont été proposées pour surmonter cet Chapitre 1 : Etude bibliographique 6 inconvénient .Afin d’éviter la difficulté de la continuité C , beaucoup de chercheurs ont formulé des éléments de plaque basés sur la théorie de Reissner/Mindlin où la continuité C est seulement requise pour l’approximation des déplacements et rotations [64,135].Les éléments de plaque isoparamétriques basés sur la théorie de Reissner/Mindlin donne une mauvaise solution quand l’intégration de Gauss complète est employée. Cela est du à la mauvaise représentation des déformations de cisaillement transversal (ce qui conduit au blocage en cisaillement transversal).Des traitements pour éviter le blocage en cisaillement transversal ont été proposés : ▪ Introduction des modes uploads/Ingenierie_Lourd/ gha-lahbari-noureddine 1 .pdf