L'objectif de ce travail est de proposer une méthodologie simple pour évaluer l'intégrité et la durée de vie d'échangeurs thermiques soudés. Une approche à deux échelles est proposée. Une description macroscopique avec la prise en compte de la structure de l'échangeur est menée pour permettre des calculs thermoélastiques par éléments finis. La réponse mécanique de l'échangeur pour deschargements thermiques, cycliques, simples est évaluée. Notamment, les zones de concentration decontraintes sont repérées. A partir cette étude, une étude micromécanique du comportement dumatériau composant l'échangeur est menée. Le matériau considéré est un acier 316L. Soncomportement élastoplastique est identifié avec un écrouissage isotrope et cinématique. La tenuemécanique pour des chargements en fatigue oligocyclique est évaluée à l'aide d'un dispositif deflexion 4 points alternée et un critère de Manson-Coffin est identifié. Ce critère est utilisé pour évaluerle nombre de cycle admissible par l'échangeur pour une amplitude de chargement donnée. Pour cela,la déformation plastique attendue dans l'échangeur est évaluée à partir d'une équivalence en énergieaux endroits où la contrainte se concentre. Les prédictions du modèle ont été comparées de manièresatisfaisante avec les résultats expérimentaux menés sur un échangeur test, pour la réponsethermoélastique que pour l'évaluation du nombre de cycles à rupture.

Dans cette thèse nous nous intéressons à la modélisation, analyse mathématique et simulation numérique d'une classe de modèles combinant différents phénomènes dissipatifs liés à la plasticité, rupture et dissipation visqueuse.Tout d'abord, nous construisons des modèles d'évolution contenant plasticité, viscoplasticité, écrouissage cinématique linéaire et rupture. En particulier, nous montrons une inégalité thermodynamique de type Clausius-Duhem pour nos modèles. Ensuite, nous montrons l'existence d'évolutions pour deux modèles: celui d'élasto-visco-plasticité avec la rupture approchée via la fonctionnelle Ambrosio-Tortorelli et celui d'élasto-viscoplasticité avec écrouissage cinématique linéaire et rupture approchée basée sur l'utilisation de la fonctionnelle d'Ambrosio-Tortorelli avec un r-Laplacien. Enfin, nous étudions numériquement nos modèles en fonction de différents paramètres mécaniques. Nous proposons aussi une extension de la méthode numérique de backtracking aux matériaux à mémoire. Au final, nous effectuons des comparaisons numériques entre un de nos modèles et l'expérience géophysique de plasticine de Peltzer et Tapponnier qui modélise la propagation des failles dans la crôute terrestre. / In this work, we are interested in modeling, mathematical analysis and numerical simulation of a class of models that combine several mecanisms of dissipation: plasticity, fracture and viscous dissipation. Firslty, we construct evolution models containing plasticity, viscoplasticity, linear kinematic hardening and fracture. In particular, we show for our models a Clausius-Duhem like thermodynamical inequality. Then, we prove an existence result for evolutions for an elasto-visco-plastic model with regularized fracture using the Ambrosio-Tortorelli functional and for an elasto-viscoplastic model with kinematic hardening and fractures regularized with the modified r-Laplacian Ambrosio-Tortorelli functional. Finally, we study from a numerical point of view our models in function of various mecanical parameters. We also propose an extension of the backtracking algorithm for materials with memory. In the end, we test numerically one of our models on a geophysical Peltzer and Tapponnier's experiment of plasticine that models failure propagation in the Earth crust.

Les capacités de résistance de l'argilite Callovo-Oxfordian (COx), qui est une roche hôte potentielle pour le dépôt souterrain profond de déchets radioactifs de haute activité en France, sont étudiées. À une échelle microscopique, des micros pores peuvent être observés dans la matrice. Une première étape d'homogénéisation a été réalisée afin d'évaluer le critère de résistance de la matrice. L'analyse microstructurale de ce matériau à quelques centaines d'échelle, référencée échelle échelle mésoscopique, montre une matrice argileuse et une distribution aléatoire d'inclusions minérales (quartz et calcite).Dans le but de déterminer le domaine de résistance à l'argilite COx, un premier outil numérique a été développé dans le contexte du comportement élastoplastique de la matrice. Plusieurs modèles morphologiques du volume élémentaire représentatif ont été considérés, et soumis à un chargement incrémental dans des conditions périodiques jusqu'à la charge limite. A la suite de ce calcul élastoplastique, un point de la frontière du domaine de résistance est obtenu. Ce dernier est alors obtenu par des calculs élastoplastiques successifs.Une alternative aux simulations élastoplastique directes, des approches cinématiques et statiques du calcul à la rupture sont réalisées. Une méthode du type éléments finis basée sur la construction d'un champ de contrainte (dans l'approche statique) et d'un champ de vitesse (dans l'approche cinématique) est développé dans un outil numérique permettant de calculer une limite inférieure et une limite supérieure de domaine de résistance / The strength capacities of Callovo-Oxfordian (COx) argillite which is a potential host rock for the deep underground repository of high-level radioactive waste in France are investigated. At a micro-scale, micro-pores can be observed in the matrix. A first strength homogenization step has been performed in order to evaluate the matrix strength criteria. The microstructure analysis of this material at some hundreds of micromet scale, referred at meso-scale, shows a clay matrix and a random distribution of mineral inclusions (quartz and calcite).Aiming to the determination of COx argillite strength domain, an FEM numerical tool has been developed in the context of the elastoplastic behavior of the matrix. Several morphological patterns of the representative elementary volume have been considered and subjected to an incremental loading in periodic conditions until collapse occurs. As a result of such elastoplastic calculation, one point of the boundary of the strength domain is obtained. The latter then could be reached by successive elastoplastic calculations.As an alternative to direct elastoplastic simulations, kinematic and static approaches of limit analysis are performed. The stress-based (static approach) and the velocity-based (kinematic approach) finite element method are used to develop a numerical tool able to derive a lower bound and upper bound of strength domain, respectively

Le récent développement des menaces terroristes renforce l'effort de recherche du CEA et d'EDF pour la protection des citoyens et des installations. De nombreux scénarios doivent être envisagés comme, par exemple, la chute d'un avion de ligne sur une structure de génie civil. La dispersion du carburant dans la structure, son embrasement sous forme de boule de feu et les effets thermiques associés sont des éléments essentiels du problème. L'utilisation de modèles numériques est indispensable car des expériences seraient difficiles à mettre en œuvre, coûteuses et dangereuses. Le problème type que l'on cherche à modéliser est donc l'impact d'un réservoir rempli de fluide, sa déchirure et la dispersion de son contenu. C'est un problème complexe qui fait intervenir une structure mince avec un comportement fortement non-linéaire allant jusqu'à rupture, un fluide dont la surface libre peut varier drastiquement et des interactions fluide-structure non permanentes. L'utilisation des méthodes numériques traditionnelles pour résoudre ce problème semble difficile, essentiellement parce qu'elles reposent sur un maillage. Cela complique la gestion des grandes déformations, la modélisation des interfaces variables et l'introduction de discontinuités telles que les fissures. Afin de s'affranchir de ces problèmes, la méthode sans maillage SPH (\og Smoothed Particle Hydrodynamics \fg) a été utilisée pour modéliser le fluide et la structure. Ce travail, inscrit dans la continuité de recherches précédentes, a permis d'étendre un modèle de coque SPH à la modélisation des ruptures. Un algorithme de gestion des interactions fluide-structure a également été adapté à la topologie particulière des coques. Afin de réduire les coûts de calcul importants liés à ce modèle, un couplage avec la méthode des éléments finis a également été élaboré. Il permet de n'utiliser les SPH que dans les zones d'intérêt où la rupture est attendue. Finalement, des essais réalisés par l'ONERA sont étudiés pour valider la méthode. Ces travaux ont permis de doter le logiciel de dynamique rapide Europlexus d'un outil original et efficace pour la simulation des impacts de structures minces en interaction avec un fluide. Un calcul démonstratif montre enfin la pertinence de l'approche et sa mise en œuvre dans un cadre industriel. / The recent development of terrorist threats increases the research effort of the french Atomic Energy Commission (CEA) and the French Electricity company (EDF) for the protection of citizens and facilities. Many scenarios should be considered as, for example, the fall of an airliner on a civil engineering structure. The dispersion of fuel in the structure, the formation of a fireball and associated thermal effects are essential elements of the problem. The use of numerical models is essential because experiences would be difficult to organize, costly and dangerous. The typical problem that we want to model is the impact of a tank filled with fluid, its rupture and the dispersion of its contents. It is a complex problem which involves a thin structure with a highly non-linear behavior up to rupture, a fluid with a free surface that can vary drastically and non permanent fluid-structure interactions. The use of traditional numerical methods to solve this kind of problems is difficult, mainly because they rely on a mesh. This complicates the management of large deformations, the modeling of moving interfaces and the introduction of discontinuities such as cracks. To overcome these problems, the meshfree method SPH (Smoothed Particle Hydrodynamics) was used to model both the fluid and the structure. This work, which is a continuation of previous research, has extended a model of SPH shell to the modeling of ruptures. An algorithm for managing fluid-structure interactions has also been adapted to the particular topology of shells. To reduce the important computational costs associated with this model, a coupling with the finite element method was also developed. It allows the use of SPH in areas of interest where the rupture is expected. Finally, tests performed by the french Aerospace Lab (ONERA) are studied to validate the method. This work helped to provide fast dynamic software Europlexus an original and effective tool for the simulation of the impact of thin structures interacting with fluid. A demonstrative calculation finally shows the relevance of the approach and its use within an industrial framework.

Le développement actuel de l’industrie vise à prévoir l’intégrité des structures dans le temps ou dans le cas de sollicitation extrême. Les risques liés à la propagation des fissures dans le cas de chocs ou d’impacts sont encore difficiles à prévoir. Les codes de calcul dans ce domaine regroupent plusieurs méthodes de simulation au sein d’un même code de calcul. Afin de présenter les différentes méthodes numériques mises en oeuvre, ce mémoire a été découpé en trois parties distinctes. Dans la première partie, nous présentons la bibliographie, puis notre apport aux méthodes de simulation numérique en l’appliquant au cas de la propagation de fissure dynamique et enfin les résultats obtenus à partir des méthodes proposées. Nous comparons ces simulations à des résultats expérimentaux ou à des simulations 2D trouvés dans la littérature. À travers la bibliographie, nous présenterons la théorie de la mécanique de la rupture pour arriver à un critère de propagation de fissure adapté à la dynamique transitoire. Ce critère a déjà été utilisé pour la fissuration dynamique en 2 dimensions. Nous décrirons la méthode des éléments finis étendus utilisée jusqu’ici principalement en quasi-statique. Nous donnerons les avantages mais aussi les limites de mise en oeuvre de cette méthode, notamment à travers le choix des enrichissements et de l’intégration des éléments coupés par la fissure. La méthode des level-sets est ensuite présentée : elle permet de décrire et faire évoluer la fissure indépendamment de la structure. On met en évidence le besoin de robustesse pour faire évoluer la fissure en dynamique explicite. La seconde partie est consacrée au développement et à l’extension de la méthode en 3D. Après avoir rappelé le critère de propagation en 3D fragile et avec plasticité, on cherche à proposer des schémas d’intégration spatiale plus économiques. Une nouvelle stratégie de propagation des level-sets basé sur la géométrie est proposée pour la dynamique explicite 3D. Enfin dans la troisième partie, nous appliquerons les méthodes à des cas de propagation de fissure bidimensionnelle puis tridimensionnelle. Nous simulerons dans un premier temps des cas 2D en mode I puis en mode mixte, afin de vérifier que l’on arrive à résultats proches des cas déjà simulés en 2D. Pour terminer par des simulations de propagation tridimensionnelle de fissure avec arrêt et redémarrage de la fissure. Tous ces développements on été implémentés dans le code de calcul de dynamique explicite EUROPLEXUS, co-propriété du CEA et de la Commission Européenne. / The current development of the industry focus on structural integrity over time or in the case of extremes stresses. Risks related to the cracks propagation in the event of shocks or impacts are still difficult to predict. Computing codes in this area groups several methods of simulation within the same computer software. To present the various numerical methods used, this thesis was divided into three distinct parts. In the first part we present the literature. Then, in second part, our contribution to the numerical simulation methods are presented by applying it to the case of dynamic crack propagation. Finally the results obtained from the proposed methods are described. We compare these simulations with experimental results or 2D simulations found in the literature. Through the first part, we present the theory of fracture mechanics to reach a criterion of crack propagation adapted to the transient dynamics. This criterion has been used for dynamic cracks in two dimensions. We describe the extended finite element method mainly used for quasi-static problems. We give the advantages but also the limits of this method: the choice of enrichment and the integration method are particularly important. The level-sets method is then presented: it allows to describe and develop the crack regardless of the structure. It highlights the need of robustness due to explicit dynamics scheme. The second part is devoted to the development and extension of the method in 3D. After reminding the propagation criterion in 3D, we try to offer more economic patterns of spatial integration. A new strategy of level-sets propagation based on geometrical approach is proposed for the explicit dynamic and applied in 3D. In the third part, we apply the methods to the case of two-dimensional crack propagation and three-dimensional. We initially simulate 2D mode I then mixed mode, to ensure that we arrive at results close to earlier 2D simulations. To finish, we present three-dimensional simulations of crack propagation with stopping and restarting crack. All these developments have been implemented in the computing software EUROPLEXUS , co-owned by the CEA and the European Commission.

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.

Les résines époxy renforcées par des nodules thermoplastiques sont largement utilisées dans les matériaux composites à matrice organique de l’industrie aéronautique. Dans le cas particulier des stratifiés, les résines époxy déterminent la résistance du milieu inter-laminaire au délaminage. Ce processus de décohésion macroscopique apparaît classiquement lors d’un impact dynamique et compromet fortement l’intégrité des structures composites légères. Pourtant, les modèles cohésifs actuels ne sont pas prévus pour la simulation du délaminage sous sollicitation rapide. Notamment, l’influence présumée de la dynamique sur la rupture du milieu inter-laminaire riche en résine époxy n’est pas prise en compte. Par conséquent, cette recherche consiste à caractériser l’initiation et la propagation de fissure dans une résine époxy sous chargement dynamique. Dans cet objectif, des éprouvettes entaillées sont fabriquées à partir de la résine époxy Hexply RM21 spécifique du domaine aéronautique. Un protocole expérimental est développé et validé afin d’estimer les ténacités d’initiation et de propagation d’une éprouvette chargée en flexion 3-points par un vérin hydraulique. Les essais réalisés démontrent que la ténacité dynamique augmente fortement avec la vitesse de propagation de fissure, à l’inverse de la ténacité d’initiation qui diminue avec la vitesse d’impact. Ces évolutions sont cohérentes avec les micrographies réalisées post-mortem qui présentent une augmentation des mécanismes de rupture. Finalement, une loi cohésive est définie et identifiée pour prendre en compte l’influence de la dynamique sur l’initiation et la propagation d’une fissure dans la résine époxy étudiée. / Epoxy resins toughened with embedded thermoplastic particles are classically used in the aeronautical composites manufacturing. As matrix into laminated composites, epoxy resins determine the delamination resistance of the interlaminar field. The delamination process is a large decohesion that often occurs during impact loadings and it critically compromises the integrity of the light weight composite structures. However, most of the cohesive zone models are not intended to simulate delamination under high rate loadings. In particular, these models do not consider the expected effects of loading rate on fracture behaviour of the resin-rich interlaminar field. Therefore, this research aims at characterizing crack initiation and propagation in epoxy resins under dynamic loadings. For that purpose, we manufacture notched specimens made of pure Hexply RM21 epoxy resin which is mainly used in aeronautics. An experimental procedure is developed and validated in order to estimate initiation and propagation fracture toughness on a 3-points bending test associated to a servo-hydraulic testing machine. Experimental results demonstrate that the dynamic fracture toughness highly increases with the crack propagation speed while the initiation toughness decreases with the impact velocity. These evolutions are consistent with post-mortem micrographic observations that show intensifying fracture mechanisms. Finally, a cohesive law is defined and identified to take into account the dynamic dependency of crack initiation and propagation in the considered epoxy resin.

Les dômes en treillis sont utilisés dans la construction de stades et d'aéroports pour accueillir un grand nombre de personnes. Ces structures entrent dans la catégorie d'importance sismique élevée car elles servent d'abris en cas d'événement sismique. Il existe des informations basées sur la recherche concernant la performance sismique des dômes ; cependant, l'application pratique de ces informations pour la conception reste un défi de taille. Essentiellement, il n'y a pas de directives établies pour la conception sismique des dômes. L'utilisation d'alliages d'aluminium peut présenter une solution efficace pour les dômes exposés à des environnements agressifs tels que les installations de stockage de produits chimiques, les dômes couvrant de grands espaces tels que les stades, et pour les dômes où des formes extrudées spécialisées peuvent faciliter les connexions efficaces entre les éléments structurels. La préférence pour ce matériau s'explique par sa durabilité, son faible rapport résistance/poids et son extrudabilité. Cependant, les alliages d'aluminium sont également connus pour être plus susceptibles à la rupture par fatigue oligocyclique. La fatigue oligocyclique peut conduire à la rupture d'un élément structurel sous l'effet de cycles de déformation induits par un tremblement de terre. La vulnérabilité des dômes en aluminium aux mouvements du sol dus aux tremblements de terre peut être affectée par la résistance à la fatigue oligocyclique du matériau. Cette étude de recherche évalue la performance sismique et la vulnérabilité sismique d'un dôme en treillis en alliage d'aluminium, en comparaison avec un dôme en treillis en acier. La vulnérabilité sismique des dômes en aluminium et en acier sous la même charge de gravité est comparée en développant des fonctions de fragilité basées sur des analyses dynamiques incrémentales. Pour les analyses dynamiques incrémentales, une suite de onze mouvements de sol non échelonnés avec une accélération maximale du sol allant de 0,2 g à 0,82 g, a été considérée. Il a été observé que les caractéristiques modales des deux dômes étaient similaires. Il n'y avait pas de déformation plastique dans les deux dômes soumis à la suite sélectionnée de mouvements de sol non échelonnés. Cependant, il y avait des déformations plastiques dans les deux dômes soumis à la suite de mouvements du sol mis à l'échelle à des intensités sismiques plus élevées. Les dômes en aluminium ont subi un déplacement deux fois plus important que le dôme en acier pour toutes les intensités sismiques considérées. Les résultats ont montré que la rupture par fatigue est attendue dans les dômes en aluminium qui sont soumis à des tremblements de terre sévères, avec une accélération spectrale supérieure à 1,5 g. Il a également été observé que le dôme en aluminium démontré une bonne résistance sismique pour les intensités de mouvement du sol représentatives du spectre de conception de l'Ouest du Canada, en particulier Vancouver, BC. / Lattice domes are used to construct stadiums and airports that accommodate many people. These structures fit into the high seismic importance category because they serve as shelters during seismic events. Information based on research regarding the seismic performance of domes exists; however, the practical application of this information for design remains a daunting challenge. Essentially, there are no established guidelines for the seismic design of domes. The use of aluminium alloys may present an effective solution for domes exposed to aggressive environments such as chemical storage facilities, domes covering large spaces such as stadiums, and for domes where specialized extruded shapes may facilitate efficient connections between structural elements. The preference for this material is due to its durability, low strength-to-weight ratio and extrudability. However, aluminium alloys are also known to be more susceptible to failure under low-cycle fatigue. Low cycle fatigue may lead to the rupture of a structural member under earthquake-induced strain cycles. The vulnerability of aluminium domes under earthquake ground motions may be affected by the low cycle fatigue resistance of the material. This research study assesses the seismic performance and seismic vulnerability of an aluminium lattice dome compared to a steel lattice dome. The seismic vulnerability of aluminium and steel domes under the same gravity load are compared by developing fragility functions based on incremental dynamic analyses. For the incremental dynamic analyses, a suite of eleven unscaled ground motions with peak ground acceleration ranging from 0.2 g and 0.82 g was considered. It was observed that the modal characteristics of both domes were similar. There was no plastic deformation in both domes subjected to the selected suite of unscaled ground motions. However, there were plastic deformations in both domes subjected to the suite of ground motions scaled to higher seismic intensities. The aluminium domes experienced twice as much displacement as the steel dome for all seismic intensities considered. The results also showed that fatigue failure is expected in aluminium domes subjected to severe earthquakes, with a spectral acceleration greater than 1.5 g. It was also observed that the aluminium dome showed a good seismic resistance for ground motion intensities representative of the design spectrum of western Canada, particularly Vancouver, BC.

On s'intéresse plus particulièrement à la modélisation du comportement d'un matériau hétérogène méso-fissuré (béton, roche,...), soumis à une sollicitation thermo-hydro-mécanique avec prise en compte du couplage géométrique. Pour conduire cette étude, on s'appuie notamment sur les approches micro-mécaniques du comportement des milieux méso-fissurés non saturés développées depuis quelques années au Laboratoire des Matériaux et des Structures du Génie Civil - Ur Navier - Université Paris Est. Le milieu fissuré non saturé traité ici est constitué d'une matrice solide homogène élastique linéaire et de fissures connectées saturées par deux fluides immiscibles : un liquide et un gaz séparés par une surface capillaire. La fissure est traditionnellement considérée comme une cavité ellipsoïdale (cas 3D) ou elliptique (cas 2D) dont le rapport d'aspect tend vers zéro. Deux morphologies typiques de matériau sont considérés dans ce travail : la situation où les fissures sont toutes orientées dans la même direction et la situation où les fissures possèdent des orientations aléatoires. Dans une première étape, on rappelle brièvement les résultats disponibles concernant la modélisation des fissures non saturées par des cavités ellipsoïdales aplaties. A la fin de cette première partie, on complète les résultats déjà disponibles en étudiant l'influence de l'histoire de chargement sur la réponse de matériau. Dans une deuxième étape, on s'attache à valider une partie des résultats obtenus en utilisant une description des efforts capillaires dans les fissures par une précontrainte homogène en seréférant aux solutions analytiques exactes disponibles dans la littérature permettant de décrire le comportement d'une fissure isolée au sein d'une matrice élastique. Dans une troisième étape, on s'intéresse aux phénomènes de propagation des fissures en condition non saturée. Les lois de propagation sous critique et le phénomène de branchement des fissures sont également prises en compte dans cette approche. La dernière partie de la thèse concerne l'influence de la température sur le comportement des milieux poreux non saturés / The main topic of my work is the development of a micromechanical model for the behaviour of unsaturated mesocracks in media (concrete, rock...) in which the thermo-hydro-mechanical loadingsand thermo-hydro-mechanical couplings are taken into account. For this, we used the micromechanical approach model of behaviour of cracked porous media recently developed at LMSGC. My thesis is focused on the equilibrium configurations of a porous material whose pore space is saturated by a vapour and a liquid phase. The behaviour of an elastic medium containing unsaturated mesocracks is studied in the framework of a micromechanical approach. The cracks are filled by two immiscible fluids, namely a liquid and a gas, separated by a capillary interface. Furthermore, it is assumed that the set of cracks constitutes a connected network ; the capillary pressure is uniform over a representative elementary volume. The cracks are modelled as flat oblate spheroid cavities. Several geometrical configurations of cracks in porous media are considered in the framework of Eshelby-based homogenization methods (parallel cracks, randomly oriented cracks). First, a previously developed model showed that when coupling between the deformation of the cracks and the capillary forces is taken into account, there is no more a one-to-one relationship between the loading parameters and the state-variables. Thus, we describe the loading history prescribed to the material in order to compute its response. Second, we validate these results referring to the exact solutions available in the literature to describe the behaviour of a unsaturated crack within an elastic matrix. Third, the description of crack propagation in unsaturated media is considered in the framework of linear elastic fracture mechanics. The phenomenon of subcritical crack growth due to stress corrosion cracking is taken into account in this approach. Mixed mode fracture in the plane is also examined. Finally, we are interested in the influence of the temperature on the behavior of unsaturated porous media in the framework of the micromechanical approach

L’utilisation des goujons collés dans les structures bois répond au souci de conservation du bâti et de discrétion de l’intervention. A ce jour, plusieurs procédés de caractérisation et de dimensionnement sont disponibles, sans pour autant donner un socle commun à l’évaluation de la résistance des assemblages collés. L’utilisation des goujons collés suscite des interrogations concernant leur tenue au feu. Bien que le matériau bois entourant les tiges de renforcement soit isolant, il est nécessaire de fournir de plus amples estimations concernant la tenue mécanique des polymères pour différentes températures que l’on pourrait trouver au sein d’une liaison au cours d’un incendie. La présente étude permet de coupler méthodes expérimentales et simulations numériques, afin d’appréhender les mécanismes gouvernant la rupture de ces assemblages. La caractérisation expérimentale permet d’estimer les propriétés mécaniques locales des assemblages suivant divers paramètres (les longueurs de scellement, l’orientation du fil du bois, l’essence de bois, ou encore la température d’exposition), ainsi que les propriétés intrinsèques aux matériaux constitutifs. Cette base de données expérimentale est indispensable pour l’ajustement du modèle aux éléments finis en élasticité linéaire. L’approche numérique met en évidence une présence importante de contraintes normales en tête de collage, avant l’apparition de contraintes de cisaillement. Les outils de la Mécanique Linéaire Élastique de la Rupture équivalente permettent d’établir des courbes de résistance liées à chaque mode de ruine (mode I et mode II). Enfin, afin de décrire précisément le processus complet de rupture de ces assemblages, un critère de rupture en mode mixte (mode I et mode II) est utilisé. Une formulation analytique permettant d’estimer la charge au pic est proposée et permettra la réalisation d’abaques de dimensionnement des assemblages par goujons collés, utilisables en bureau d’études. / Glued-in-rods lead to overcome the use of traditional bolted connections, preserve a large part of original timber and offer aesthetic benefits (since the repair is hidden in the cross sections of the members). Despite previous research programs in many countries, some design rules, predicting the axial strength of such connections, are available, but a common criterion is still lacking. However, the durability of this process according to fluctuating temperature is not well known. During fire exposure, connections are not directly in contact with flames, since they are isolated by surrounding wood. The current study combines experiments with finite element computations, in order to lead a better to a better knowledge about their mechanical and fracture behaviors. An important experimental campaign is carried out on such connections and provides the influence of various parameters, such as the anchorage length, the rod-to-grain angle, the specie of wood or the temperature exposure, on their mechanical behaviors. Moreover, the inherent mechanical properties of the rod and the adhesives used are also studied. The finite element modeling reproduces the experimental configuration, and reveals significant tensile stresses in comparison with shear stresses. Within the framework of equivalent linear elastic fracture mechanics, R-curves in mode I and mode II can be estimated for each specimen. Finally, a fracture criterion in mixed mode is used to describe the complete fracture process of glued-in rods. An analytical formulation is then proposed and allows the evaluation of peak load of each specimen. This approach leads to realize design tables, usable by design offices.