L'objectif de cette thèse est de développer des techniques de modélisation et de simulation multi-échelle avancées et efficaces pour étudier les phénomènes d'instabilité dans trois structures d'ingénierie courantes: membrane, film/substrat et structures sandwich, en combinant la technique des coefficients de Fourier lentement variables (TSVFC) et la méthode numérique asymptotique (ANM). À cette fin, basée sur les équations de la plaque de Von Karman, la TSVFC été utilisée pour développer un modèle de Fourier à bidimensionnel (2D) qui a également été implémenté dans ABAQUS via sa sous-routine UEL. Ensuite, un 2D modèle de Fourier est construit pour le film/substrat. En outre, en utilisant leurs caractéristiques de déformation, un 1D modèle de Fourier est développé en utilisant à la fois le TSVFC et le CUF. Par la suite, sur la base d'une cinématique Zig-Zag d'ordre supérieur, un 2D modèle de Fourier est déduit pour une plaque sandwich. Les équations directrices pour les modèles ci-dessus sont discrétisées par la méthode des éléments finis, et les systèmes non linéaires résultants sont résolus par le solveur non linéaire efficace et robuste ANM. Ces modèles sont ensuite adoptés pour étudier les instabilités dans ces structures. Les résultats montrent que les modèles établis peuvent simuler avec précision et efficacité divers phénomènes d'instabilité. En outre, on constate que l'instabilité membranaire est sensible aux conditions aux limites et qu'il existe un paramètre sans dimension presque constant près du point de bifurcation pour différents cas de charge et paramètres géométriques, ce qui peut être utile pour prédire rapidement l'apparition des rides / The main aim of this thesis is to develop advanced and efficient multi-scale modeling and simulation techniques to study instability phenomena in three common engineering structures, i.e., membrane, film/substrate and sandwich structures, by combining the Technique of Slowly Variable Fourier Coefficients (TSVFC) and the Asymptotic Numerical Method (ANM). Towards this end, based on the Von Karman plate equations, the TSVFC has been firstly used to develop a two-dimensional (2D) Fourier double-scale model for membrane, which has also been implemented into ABAQUS via its subroutine UEL. Then a 2D Fourier model is constructed for film/substrate. Further, making use of deformation features of the film/substrate, a 1D Fourier model is developed by using both the TSVFC and the Carrera’s Unified Formulation (CUF). Subsequently, based on high-order kinematics belonging to Zig-Zag theory, a 2D Fourier model is deduced for sandwich plate. The governing equations for the above models are discretized by the Finite Element Method, and the resulting nonlinear systems are solved by the efficient and robust nonlinear solver ANM. These models are then adopted to study instabilities in these structures. Results show that the established models could accurately and efficiently simulate various instability phenomena. Besides, it’s found that the membrane instability is very sensitive to boundary conditions, and there exists a dimensionless parameter that is almost constant near bifurcation point for various loading cases and geometric parameters, which may be helpful for fast predicting the occurrence of wrinkles

Malgré le cumul de connaissances sur le sujet du flambage des coques minces, des questions essentielles demeurent. En effet, malgré les avancées et une meilleure compréhension de l’effet des défauts, le dimensionnement continue de recourir aux anciennes règles qui découlent d’une démarche empirique et qui s’avèrent souvent trop conservatives comme c’est le cas par exemple pour la NASA SP8007 (1968). Cette règle est utilisée notamment pour le dimensionnement du lanceur Ariane 5. L’Etage Principal Cryogénique est en effet constitué de coques cylindriques minces, qui sont soumises à une combinaison de chargements et donc sujettes au déclenchement d’instabilités pouvant être catastrophiques. Il est ainsi nécessaire d’avoir une meilleure compréhension du phénomène pour pouvoir améliorer ces méthodes de dimensionnement dans le cas de coques moyennement longues (1 < L/R < 3) et minces (250 < R/t <1500). Nous employons pour cela une approche à la fois numérique et expérimentale. L’outil numérique est utilisé, via une modélisation pertinente, afin de construire une nouvelle règle de dimensionnement et d’étudier l’influence des différents paramètres (géométriques, matériau). L’aspect expérimental a pris une place prépondérante, une large campagne nous permettant de valider les résultats simulations pour différentes configurations, mais également d’avoir une bonne compréhension du phénomène. Ces deux aspects de notre recherche nous ont également permis de mieux déterminer l’interaction entre les différents chargements (pression interne, compression, flexion, cisaillement), plus ou moins comprise selon les cas. / The stability of thin cylindrical shells has been studied for decades, but despite the accumulated knowledge on the subject, important questions remain, in particular regarding the design of such structures. The effect of geometrical imperfections is now well understood, but old design recommandations resulting from purely empirical approaches, are still used. The NASA SP8007 recommendation (1968) is one them, and leads to conservative designs. This rule is used in particular for the design of the Ariane 5 launcher. Indeed, the Principal Cryogenic Stage is made of thin cylindrical shells, subjected throughout the structure’s life cycle to a combination of loadings, and hence prone to possibly disastrous structural instabilities. It is is then necessary to better understand those phenomena in order to improve the design methods for relatively long (1 < L/R < 3) and thin (250 < R/t < 1500) shells. Our work is based on both numerical and experimental approaches. Numerical tools are used, through a relevant choice of the critical imperfections, to build a new design recommendation and study the influence of different parameters (geometrical or material). Experiments took an important place in our research, and the large experimental campaign allowed us to validate the numerical results for different configurations, and to better understand the physics of the problem as well. These two aspects of our study also enabled us to clarify the interaction between different loadings (internal pressure, compression, bending and transverse shear), which is today more or less understood, according to the case.

Lorsqu’on découpe un feuillard à l’aide d’un outil, ou lorsqu’on découpe un tube selon son axe, au fur et à mesure que l’on propage la fissure qui traduit la découpe il arrive que des ondulations de flambage perturbent les deux bords libres générés par la propagation de la fissure. Cette étude vise à analyser les origines de ces ondulations. Nous avons mené une campagne expérimentale, dans laquelle des tubes en acier inox avec différentes géométries (rayon/épaisseur) sont « découpés » selon une génératrice. Une instrumentation adéquate, couplant des mesures ponctuelles, à l’aide de jauges de déformation, et une méthode champ par corrélation d’image, nous a permis de correctement mettre en exergue la phénoménologie, en particulier les cinématiques induites à l’échelle géométrique de la fissure (front de fissure) ainsi qu’à l’échelle du tube, avec les longueurs d’onde de flambage observées à l’aval de la fissure. La modélisation numérique menée en non linéaire géométrique (flambage), matériau (déchirure ductile), et conditions aux limites (contact) est aussi abordée à l’aide du code de calcul Abaqus/Standard. Pour la gestion de la propagation de la fissure, deux modèles de rupture sont proposés. Le premier modèle dit zone cohésive est développé et implanté dans le code Abaqus via la subroutine UEL. Pour la deuxième modélisation, nous avons utilisé le modèle dit « d’endommagement ductile » du code Abaqus. La modélisation via des éléments massifs ou des éléments coques volumiques ainsi que l’utilisation de ces modèles de rupture permettent de corroborer les observations expérimentales. Ces travaux montrent que l’augmentation de la charge inhérente au déplacement de l’outil de « découpe », induit une extension dans la direction circonférentielle et donc une striction dans la direction radiale amenant finalement la rupture. Lors de la rupture, un « sillage plastique » apparait, relativement large, près et parallèle aux bords de la fissure. « Confiné » par les autres parties du tube qui restent élastiques, des contraintes de compression axiale résiduelles apparaissent dans ce sillage plastique, à l’aval de la fissure, leur intensité est suffisante pour produire les ondulations des bords libres qui traduisent un flambage local. Les contraintes résiduelles liées à l’opération de découpe induisent donc le flambage. / When one uses a tool to cut a sheet metal, or a tube according to his axis, as one propagates the crack which translates cutting it arrives that undulations of buckling disturb the two free edges generated by the propagation of the crack. This study aims at analyzing the origins of this behavior. We conducted an experimental campaign, in which stainless steel tubes with various geometries (radius/thickness) are « cut out » according to a generator. An adequate instrumentation, coupling of specific measurements, using strain gauges, and a field method, by digital image correlation, allowed us accurately to put forward phenomenology, in particular the kinematics at the scale of the crack (ahead of crack tip) and at the level of tube, with the wavelengths of buckling observed at the downstream of the crack tip. The numerical modeling taking into account nonlinearities of material (ductile tear), geometry (buckling) and boundary conditions (contact) is also approached using the code Abaqus/Standard. For the management of the crack propagation, two rupture models are proposed. The first model called cohesive zone is developed and implemented in the Abaqus code via the user routine UEL. For the second modeling, we used the model called “ductile damage model” in the Abaqus code. Modeling via solid elements or shell continuum elements as well as the use of these rupture models make it possible to corroborate the experimental observations. These studies show that the increase of the load inherent in the displacement of the tool of « cutting » induced a circumferential extension of the tube that leads to a local necking in the radial direction bringing the rupture finally. During the failure, a “plastic wake” appears, relatively wide, close and parallel to crack lips. Constrained by other parts of the tube which remain elastic, sufficient axial residual compressive stresses produced in this plastic wake produce the undulations which represents a local buckling. The residual stresses related to the operation of cutting thus induce buckling.

De nombreuses structures spatiales telles que les lanceurs sont équipées d’une mousse de protection thermique. Cette couche de matériau extrêmement léger présente d’excellentes propriétés d’isolation thermique mais des caractéristiques mécaniques très faibles. Il est ici proposé d’analyser le comportement au flambage de ce type de construction où une fine et légère structure, coque cylindrique en aluminium avec un ratio R/t de 665, est recouverte par une épaisse couche de mousse de faible densité. Afin d'évaluer l'effet du comportement au flambage de ce type de coque multicouche, des essais expérimentaux et numériques sont réalisés sur des cylindres faiblement pressurisés soumis à du cisaillement et de la flexion. Ce cas de charge représente la configuration « cas sol », d’un lanceur sur le pas de tir, dans l’attente du lancement, et soumis à des charges mécaniques extrêmes générées par le vent. / Many space structures such as launchers are equipped with a thermal foam protection barrier. This layer, of extremely light material generally exhibits excellent properties of heat insulation but very weak mechanical characteristics. This study is devoted to the analysis of the buckling behavior of this hybrid wall construction where the skin of a very thin light-weight structure, an aluminum cylindrical shell with an R/t ratio of about 665, is coated with a thick layer of foam with a low material density. To gauge the effect on the buckling behavior of this kind of multilayered shells, experimental and numerical tests are conducted on slightly pressurized cylindrical shells submitted to shear load. This load configuration represents the case of a rocket on the launching pad, waiting to be launched, and submitted to mechanical loads induced by the wind.

Dans cette thèse, nous nous intéressons à la biomécanique de deux mouvements parmi les plus rapides du règne végétal. La première partie porte sur l’étude théorique du mouvement permettant la capture de proies par les pièges de l’utriculaire, des outres déformables de taille millimétrique refermées par une porte flexible. Un modèle dynamique, élaboré à partir d’ingrédients mécaniques, hydrodynamiques et élastiques, permet de relier la différence de pression entre l’intérieur et l’extérieur du piège à la position de la porte via deux équations différentielles couplées. Le modèle permet de capturer la dynamique de fonctionnement du piège et de prédire l’intégralité des comportements observés dans la nature via l’ajout d’un bruit stochastique. La seconde partie est consacrée au mouvement d’éjection des spores par les sporanges de fougères. Notre étude à la fois expérimentale et théorique permet de révéler le caractère remarquable de l’anneau, une structure spécialisée du sporange, dont la nature poroélastique lui confère un comportement comparable à celui d’une catapulte autonome. Les différentes phases du mouvement : ouverture, déclenchement par cavitation, fermeture rapide et recharge éventuelle sont observées via imagerie ultra-rapide. La courbure de l’anneau au cours du mouvement est mesurée expérimentalement puis comparée aux prédictions théoriques pour chacune des phases du mouvement. Cette étude nous permet d’identifier les différents processus physiques à l’origine du mouvement et de déterminer les paramètres caractéristiques de l’anneau : raideur, perméabilité membranaire et pression osmotique interne, ainsi que la pression négative de cavitation. / In this PhD work, we focus on the biomechanics of two motions among the fastest in plant kingdom. The first part is a theoretical study of the motion leading to a prey capture by the bladderwort’s traps, elastic millimeter-sized bladders closed by a flexible door. A dynamical model, based on mechanical, elastic and hydrodynamic ingredients, links the pressure difference between the trap and its surroundings with the door position by the means of two coupled ordinary differential equations. The model captures the dynamics of the trap and predicts all the range of behaviors found in nature by including stochastic noise in the system. The second part focuses on the fern sporangium motion allowing the spores dispersal. Our experimental and theoretical studies point out the remarkable character of the annulus, a specialized structure of the sporangium, as it behaves as an autonomous catapult due to its poroelastic nature. The different stages of the motion: opening, triggering by cavitation, fast closure and eventual reloading are observed using ultra-fast imaging. The annulus curvature is measured experimentally and then compared to the theoretical predictions for each motion phase. This study enables us to identify the physical processes governing the complete motion and to determine the characteristic parameters of the annulus: stiffness, permeability and internal osmotic pressure, and also the negative cavitation pressure.

Dans l’optique d’aller toujours vers le plus petit, un réacteur plasma micro-ondes en cavité résonnante fonctionnant à la pression atmosphérique a été développé à l’Institut Jean-Lamour (IJL, UMR 7198) dans l’équipe ESPRITS (Expériences et Simulations des Plasmas Réactifs - Interaction plasma-surface et Traitement des Surfaces).Ce réacteur plasma présente la particularité de pouvoir fournir une post-décharge de petite dimension (de l’ordre du millimètre). Un mélange Ar-O2 est utilisé comme gaz composant le plasma. L’aspect filamentaire de ce type de plasma est une des caractéristiques majeures du dispositif. De plus, la post-décharge se trouve notamment composée d’oxygène atomique, espèce active très utilisée par exemple en PECVD ou pour la fonctionnalisation de surface.Notre travail a consisté à mettre au point ce réacteur plasma et à caractériser la post-décharge résultante. Après avoir maîtrisé les principaux paramètres pouvant jouer notamment sur la forme de la post-décharge, nous nous sommes intéressés à son utilisation en tant que source plasma pour des applications de traitements de surfaces.Nous avons ainsi étudié des dépôts de films minces de SiOx à partir d’hexaméthyldisiloxane (HMDSO) en utilisant la post-décharge comme source d’oxygène atomique permettant la décomposition du précurseur organosilicié. L’HMDSO est utilisé ici car il est assez simple à manipuler et il a été souvent étudié, nous permettant d’obtenir sur lui de nombreuses informations. Différentes structures auto-organisées sont apparues au sein de ces dépôts et ont été étudiées afin de comprendre leur mécanisme de formation / In order to always go to the smallest, a plasma reactor microwave in resonant cavity operating at atmospheric pressure was developed at the Institut Jean Lamour (IJL, UMR 7198) in the team ESPRITS (Expériences et Simulations des Plasmas Réactifs - Interaction plasma-surface et Traitement des Surfaces).This plasma reactor has the particularity to provide post-discharge of small size (about one millimeter). An Ar-O2 gas is used as the plasma component. The filamentous appearance of this type of plasma is a major feature of the device. Furthermore, the afterglow is notably composed of atomic oxygen, active species such as widely used PECVD or surface functionalization.Our job was to develop the plasma reactor and to characterize the resulting afterglow. After mastering the key parameters that can play on the particular shape of the afterglow, we were interested in its use as a plasma source for surface treatment applications.We studied deposition of SiOx thin films from hexamethyldisiloxane (HMDSO) using post-discharge source of atomic oxygen for the decomposition of the organosilicon precursor. The HMDSO is used here because it is fairly easy to handle and it has often been studied, allowing us to obtain much information on him. Various self-organized structures appeared in these deposits and have been investigated to understand their formation mechanism

L'optimisation topologique vise à concevoir une structure en recherchant la disposition optimale du matériau dans un espace de conception donné, permettant ainsi de proposer des designs optimaux innovants. Cette thèse est centrée sur l'optimisation topologique pour des problèmes de conception de structures prenant en compte des contraintes de flambage. Dans une large variété de domaines de l'ingénierie, la conception innovante de structures est cruciale. L'allègement des structures lors la phase de conception tient une place prépondérante afin de réduire les coûts de fabrication. Ainsi l'objectif est souvent la minimisation de la masse de la structure à concevoir. En ce qui concerne les contraintes, en plus des contraintes mécaniques classiques (compression, tension), il est nécessaire de prendre en compte des phénomènes dits de flambage, qui se caractérisent par une amplification des déformations de la structure et une potentielle annihilation des capacités de la structure à supporter les efforts appliqués. Dans le but d'adresser un large panel de problèmes d'optimisation topologique, nous considérons les deux types de représentation d'une structure : les structures treillis et les structures continues. Dans le cadre de structures treillis, l'objectif est de minimiser la masse en optimisant le nombre d'éléments de la structure et les dimensions des sections transversales associées à ces éléments. Nous considérons les structures constituées d'éléments poutres et nous introduisons une formulation du problème comme un problème d'optimisation non-linéaire en variables mixtes. Afin de prendre en compte des contraintes de manufacturabilité, nous proposons une fonction coût combinant la masse et la somme des seconds moments d'inertie de chaque poutre. Nous avons développé un algorithme adapté au problème d'optimisation considéré. Les résultats numériques montrent que l'approche proposée mène à des gains de masses significatifs par rapport à des approches existantes. Dans le cas des structures continues, l'optimisation topologique vise à discrétiser le domaine de conception et à déterminer les éléments de ce domaine discrétisé qui doivent être composés de matière, définissant ainsi un problème d'optimisation discret. [...] / Topology optimization aims to design a structure by seeking the optimal material layout within a given design space, thus making it possible to propose innovative optimal designs. This thesis focuses on topology optimization for structural problems taking into account buckling constraints. In a wide variety of engineering fields, innovative structural design is crucial. The lightening of structures during the design phase holds a prominent place in order to reduce manufacturing costs. Thus the goal is often the minimization of the mass of the structure to be designed. Regarding the constraints, in addition to the conventional mechanical constraints (compression, tension), it is necessary to take into account buckling phenomena which are characterized by an amplification of the deformations of the structure and a potential annihilation of the capabilities of the structure to support the applied efforts. In order to adress a wide range of topology optimization problems, we consider the two types of representation of a structure: lattice structures and continuous structures. In the framework of lattice structures, the objective is to minimize the mass by optimizing the number of elements of the structure and the dimensions of the cross sections associated to these elements. We consider structures constituted by a set of frame elements and we introduce a formulation of the problem as a mixed-integer nonlinear problem. In order to obtain a manufacturable structure, we propose a cost function combining the mass and the sum of the second moments of inertia of each frame. We developed an algorithm adapted to the considered optimization problem. The numerical results show that the proposed approach leads to significant mass gains over existing approaches. In the case of continuous structures, topology optimization aims to discretize the design domain and to determine the elements of this discretized domain that must be composed of material, thus defining a discrete optimization problem. [...]

Le sujet de cette thèse est l'hydrodynamique de corps minces (feuilles) et élancés (filamenteux) de fluide visqueux immergés dans un second fluide ayant une viscosité différente. Nous nous concentrons sur deux exemples : la subduction de la lithosphère océanique et le flambage de fils visqueux dans microcanaux divergents, les deux ont un nombre de Reynolds caractéristique Re<<1. Pour le cas de la subduction d'une feuille mince, nous proposons une hybride méthode «boundary integral & thin sheet» (BITS). Après la validation en comparant ses prévisions avec celles de la boundary-element méthode, deux solutions instantanées et dépendant du temps sont effectués pour analyser la subduction avec la méthode BITS. L'analyse à l'échelle de la vitesse d'immersion normalisée en fonction de «la rigidité en flexion» de la feuille est confirmée par nos prédictions numériques. Pour des rapports de viscosité modérée (≈100), la feuille amincit sensiblement quand elle coule, mais pas assez pour conduire à la «rupture de la dalle» que l'on observe dans plusieurs zones de subduction sur Terre. Ensuite, le code BLEU parallèle pour écoulements polyphasiques est utilisé à simuler pliage visqueux tridimensionnel dans des microcanaux divergent. Nous avons réalisé une étude paramétrique comprenant cinq simulations dans lequel le rapport de débit volumétrique, le rapport de viscosité, le nombre de Reynolds, et la forme de la chaîne ont été modifiées par rapport à un modèle de référence. Le fil devient instable à une instabilité de pliage en raison de la contrainte de compression longitudinale. L'axe de pliage initial peut être parallèle ou perpendiculaire à la dimension étroite de la chambre. Dans le premier cas, le pliage transforme lentement au pliage perpendiculaire au moyen d'une torsion, ou peut disparaître totalement. / The hydrodynamics of thin (sheet-like) and slender (filamentary) bodies of viscous fluid immersed in a second fluid with a different viscosity is studied. Here we focuses on two examples: the subduction of oceanic lithosphere and the buckling of viscous threads in diverging microchannels, both have a characteristic Reynolds number Re<<1. A hybrid boundary integral & thin sheet method (BITS) is build for the subduction of 2D immersed sheet. After the validation by comparing with the results of full boundary elements method, both instantaneous and time-dependant soloutions are done to analyze the subduction with the BITS method. The scaling analysis of the normalized sinking speed V/V_Stokes as a function of the sheet's 'flexural stiffness' is confirmed by our numerical predictions. For moderate viscosity ratios (≈100), the sheet thins substantially as it sinks, but not enough to lead to the ‘slab breakoff’ that is observed in several subduction zones on Earth. Next, the parallel code BLUE for multi-phases flows is used to simulate the 3-dimensional viscous folding in diverging microchannels. We performed a parameter study comprising five simulations in which the flow rate ratio, the viscosity ratio, the Reynolds number, and the shape of the channel were varied relative to a reference model. The thread becomes unstable to a folding instability due to the longitudinal compressive stress. The initial folding axis can be either parallel or perpendicular to the narrow dimension of the chamber. In the former case, the folding slowly transforms via twisting to perpendicular folding , or may disappear totally.

Nous proposons dans cette thèse des modèles éléments finis pour décrire les phénomènes de flambage que l'on rencontre souvent en laminage des tôles minces. Partant d'un modèle simplifié qui suppose le mode de flambage harmonique dans le sens du laminage, le code permet de détecter des points de bifurcations, de décrire le comportement en post-flambage ou encore d'analyser l'influence de la traction globale sur les tailles de défauts et les modes de flambage. Le modèle n'étant pas prévu pour tenir compte des chargements complexes en laminage, nous avons proposé un autre modèle plus complet tenant compte de toutes les composantes des contraintes résiduelles et capable de coupler les phénomènes en amont comme en aval de l'emprise (emprise-flambage). Les modèles existants traitent généralement un couplage itératif entre l'emprise et le flambage (Thèse d’Abdelkhalek) ou un couplage direct, mais ce dernier est limité pour représenter les modes de flambage (Thèse de Counhaye). Dans cette thèse, nous proposons un couplage direct entre l'emprise et le flambage utilisant un code de laminage LAM3 pour décrire l'emprise et un modèle coque pour décrire le flambage. Nous avons utilisé la méthode Arlequin pour coupler les deux modèles. Cette méthode de couplage très prometteuse, est l'une des plus flexibles qui traite le couplage par superposition ou collage des modèles possédant des propriétés différentes. L'originalité du modèle développé réside essentiellement dans le déplacement de la zone de couplage à chaque incrément de temps. Pour valider le modèle développé, nous avons effectué des cas tests notamment un cas industriel et des cas académiques de laminage pouvant engendrer des défauts bords longs tout comme des plis longitudinaux que l'on observe souvent à la sortie de l'emprise. Les résultats issus de ce code ont été validés avec des mesures expérimentales et avec des modèles de références. Le modèle prédit bien les contraintes relaxées après flambage et montre bien les défauts de planéité correspondants / We propose in this thesis finite element models to describe buckling phenomena that are often encountered in thin sheet rolling processes. Starting with a simplified model assuming buckling mode as being harmonic in the rolling direction, the code can detect the bifurcation points and describe post-buckling behavior. The model is not intended to reflect complex rolling loads, we proposed another more complete model taking into account all components of the residual stresses and able to couple the phenomena at the upstream of the roll mill with the buckling phenomena at the downstream domain. Existing models generally treat iterative couplings between the zone under the bite and the buckling phenomena (Abdelkhalek's thesis) or direct coupling but it is limited to represent buckling modes (Counhaye's thesis). In this thesis, we propose a direct coupling between the upstream of the roll mill and the downstream domain using a rolling code LAM3 to describe the bite and a shell model to describe buckling phenomena in the downstream domain of the sheet. We used Arlequin method which is one of the most flexible coupling techniques to couple both models. This method leads to a partition of the space, each model being valid in a part of the domain. Both models should be considered valid in intersection of two zones. The key points are the definitions of a moving coupling zone, of a relevant coupling operator and of a simple procedure to build varying meshes.To validate the proposed model, we performed some test cases including an industrial case and academic rolling test cases including edge-wave defects or local folds out of the roll mill. The results have been validated by comparison with experimental measurements and with reference models

Ce travail de recherche répond aux besoins actuels et futurs dans le domaine de l’Aérospatial de faire évoluer les règles de dimensionnement au flambage des réservoirs structuraux de l’Etage Principal Cryogénique (EPC) des lanceurs. Ces réservoirs, composés de coques cylindriques, peuvent être associés à un faible raidissage en termes de masse ajoutée. L’objectif d’alléger le lanceur pour optimiser la charge utile, conduit au choix de coques constitutives de plus en plus minces, le risque de flambage sous diverses sollicitations est donc d’autant plus accru. Le dimensionnement au flambage de l’EPC est basé principalement sur la norme NASA SP8007 qui date de 1968, et qui semble trop conservative, notamment aux basses pressions. Précisons aussi, que l’EPC est équipé d’une couche de protection thermique (PT) qui n’est pas prise en compte dans le design au vue de sa très faible rigidité de membrane. La contribution de cette couche à la capacité de flambage de la coque est cependant un sujet ouvert. / This research work responds to the current and future requirements in Aerospatiale are to improve the design for buckling of the tanks of the Cryogenic Main Stage (EPC) of the launcher. These tanks are composed of cylindrical shells and can be associated with weak stiffening, which are becoming thinner and therefore more susceptible to a risk of buckling. The buckling design of the EPC based mainly on the NASA SP8007 standard, which is ac-cording to many specialists too preservative, especially under low pressure. Moreover, the EPC is equipped with a thermal protection layer (PT), which is extremely light and has an excellent thermal insulation property but very low mechanical properties. The contribution of this layer to the buckling capacity of a lightly pressurized thin cylindrical shell under var-ious solicitations, therefore, appears to be a major subject.