Quelles sont les déformations et les contraintes que subissent une structure gonflée précontrainte lorsqu'elle est soumise à des perturbations de son environnement extérieur? Cette thèse répond à cette interrogation dans le cas d'un lobe formé de deux membranes orthotropes, fixées sur leurs longueurs à deux axes rigides, successivement rectilignes puis légèrement arrondis, et sur leurs largeurs à des arcs de cercle également rigides. L'intérieur est gonflé par un gaz supposé parfait.<br /> A partir d'une modélisation mécanique en élasticité linéaire, nous avons établi un modèle mathématique rigoureux en 3D. Puis en faisant tendre l'épaisseur des membranes vers zéro, nous avons obtenu un modèle asymptotique 2D, bien posé pour certaines forces, dites admissibles, dans des espaces obtenus par complétion. Nous avons alors démontré que la suite des valeurs moyennes dans l'épaisseur des solutions des problèmes 3D converge fortement vers l'unique solution du problème asymptotique. Nous avons de plus mis en évidence des conditions suffisantes d'admissibilité des forces extérieures.<br /> L'analyse numérique du modèle asymptotique a montré que les estimations d'erreur a priori se font dans la norme de l'énergie. Ce qui, comme l'ont confirmé quelques essais numériques, ne va pas sans poser de problèmes quand on s'intéresse aux déplacements, mais est tout à fait satisfaisant du point de vue des contraintes.

L'informatique en nuage révolutionne complètement la façon de gérer les ressources. Grâce à l'élasticité, les ressources peuvent être provisionnées en quelques minutes pour satisfaire un niveau de qualité de service (QdS) formalisé par un accord de niveau de service (SLA) entre les différents acteurs du nuage. Le principal défi des fournisseurs de services est de maintenir la satisfaction de leurs consommateurs tout en minimisant le coût de ces services. Du point de vue SaaS, ce défi peut être résolu d'une manière ad-hoc par l'allocation/-libération des ressources selon un ensemble de règles prédéfinies avec Amazon Auto Scaling par exemple. Cependant, implémenter finement ces règles d'élasticité n'est pas une tâche triviale. D'une part, la difficulté de profiler la performance d'un service compromet la précision de la planification des ressources. D'autre part, plusieurs paramètres doivent être pris en compte, tels que la multiplication des types de ressources, le temps non-négligeable d'initialisation de ressource et le modèle de facturation IaaS. Cette thèse propose une solution complète pour la gestion des contrats de service du nuage. Nous introduisons CSLA (Cloud ServiceLevel Agreement), un langage dédié à la définition de contrat de service en nuage. Il adresse finement les violations SLA via la dégradation fonctionnelle/QdS et des modèles de pénalité avancés. Nous proposons, ensuite, HybridScale un framework de dimensionnement automatique dirigé par les SLA. Il implémente l'élasticité de façon hybride : gestion réactive-proactive, dimensionnement vertical horizontalet multi-couches (application-infrastructure). Notre solution est validée expérimentalement sur Amazon EC2.

Les travaux effectués au cours de cette thèse portent principalement sur la construction de nouvelles bornes du comportement effectif des matériaux biphasés de type matrice-inclusions à comportement linéaire élastique. Dans un premier temps, afin d'encadrer le comportement effectif, nous présentons une nouvelle approche numérique, inspirée des travaux de Huet (J. Mech. Phys. Solids 1990 ; 38:813-41), qui repose sur le calcul des comportements apparents associés à des volumes élémentaires (VE) non-carrés construits à partir d'assemblages de cellules de Voronoï, chaque cellule contenant une inclusion entourée de matrice. De tels VE non-carrés permettent d'éviter l'application directe des CL sur les inclusions à l'origine d'une surestimation artificielle des comportements apparents. En utilisant les théorèmes énergétiques de l'élasticité linéaire et des procédures de moyennisation appropriées portant sur les comportements apparents, un nouvel encadrement du comportement effectif est obtenu. Son application au cas d'un composite biphasé, constitué d'une matrice isotrope et de fibres cylindriques parallèles et identiques distribuées aléatoirement dans le plan transverse, conduit à des bornes plus resserrées que celles obtenues par Huet. En nous appuyant sur cette nouvelle procédure numérique, nous avons ensuite réalisé une étude statistique des comportements apparents à l'aide de simulations de type Monté Carlo. Puis, à partir des tendances issues de cette étude statistique, nous avons proposé et mis en œuvre de nouveaux critères de tailles de VER.

Ce travail de thèse a pour objet l'étude de l'influence de la plasticité sur le délaminage et le flambage de films minces déposés sur substrats. Il repose sur une approche mixte combinant des simulations atomistiques et des calculs analytiques basés sur la théorie des plaques minces de Föppl-von Kármán (FvK). Les simulations ont permis de caractériser, au cours de la formation d'une ride droite, un mécanisme de glissement localisé dans l'interface en pied de cloque entraînant une augmentation de la déflexion maximale de la ride. Ce mécanisme de glissement est également présent lorsque le délaminage piloté par le flambage du film mince est lui aussi observé. En l'intégrant dans le modèle élastique de FvK, la forme de la ride droite ainsi que le processus de délaminage ont ensuite été caractérisés. Le bon accord trouvé entre les simulations atomistiques et le modèle explique notamment le délaminage des cloques sans introduire de dépendance entre l'énergie d'adhésion et la mixité modale. L'initiation du cloquage à partir d'une marche d'interface créée par des dislocations venant du substrat a également été étudiée. Les simulations révèlent qu'avant flambage, le film se décolle à la fois sur le haut et sur le bas de la marche. Un mécanisme de glissement est là aussi identifié. Une déformation critique de flambage qui tient compte de ces phénomènes a été déterminée en modélisant le film mince sur la marche dans le formalisme de FvK. Les résultats des simulations couplés au modèle élastique expliquent, comme il est par ailleurs observé expérimentalement, pourquoi les cloques se forment préférentiellement au-dessus de défauts tels que des marches.

La prédiction de la géométrie d'une pièce mise en forme par déformation plastique grâce à un logiciel de calcul par éléments finis (EF) s'effectue en suivant séquentiellement différentes étapes : la caractérisation thermo-mécanique du matériau, la modélisation de son comportement et son intégration dans un logiciel EF, puis la mise en données et la simulation de l'opération de formage. La phase de modélisation consiste entre autre à identifier quel type de modèle de comportement est le plus approprié pour prédire les réactions du matériau lors de l'opération de formage. Ces modèles sont essentiellement développés dans le cadre de la mécanique des milieux continus (MMC). L'hypothèse forte, si ce n'est centrale, de la MMC consiste à considérer que les variables qui servent à déterminer le comportement du matériau sont continues et dérivables. Cependant, les connaissances les plus élémentaires de métallurgie indiquent que les grandeurs locales dans les matériaux métalliques sont discontinues. La majorité des modèles de comportement mécanique des matériaux métalliques repose sur la définition d'un volume élémentaire représentatif dont la taille est assez grande pour permettre une homogénéisation de la description du comportement en gommant l'influence des hétérogénéités localesL'objet de ces travaux est de montrer que la prise en compte des hétérogénéités locales dans la modélisation du comportement des matériaux métalliques est pertinente et contribue à l'amélioration de la prédiction des simulations d'opérations de mise en forme en élargissant le potentiel prédictif des modèles ainsi construits. Un modèle élasto-plastique prenant en compte les hétérogénéités locales est alors proposé. / To obtain a relevant shape of a formed part during its finite element simulation, several steps are needed: thermo-mechanical caracterization of the material, definition of the most relevant model and integration of this model in the FE software and finally after data converting and computing processes. The modelling step include, among other things, the identification of the most appropriate model to fit the experimental material behaviour. Those models are essentially developped within the framework of continuum mechanics (CM). A strong, if not the main assumption of the CM consists in considering that mechanical description variables are continuous and differentiable. However, the basic knowledge of metallurgy indicates that local data in metallic materials are discontinuous. For metallic materials, the majority of constitutive models are based on the definition of a representative elementary volume (REV). This REV is supposed to be large enough to erase the incidence of local heterogeneities. Then those constitutive models are assumed to be homogeneous.The aim of this work is to show that introducing local heterogeneities in the description of constitutive models is relevant and contribute to improve the simulation accuracy. Those models also provide an enlargement of the simulation predictive potential. Then an elasto-plastic model, based on local heterogeneities description, is proposed.

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 matériaux super-élastiques sous forme de couche minces présentent un domaine d’applications assez large. Grâce à leur aptitude à accommoder de grandes déformations de manière réversible, les films super-élastiques peuvent être utilisés comme sous couches entre le substrat et les revêtements dures. Cela permet d’améliorer significativement de nombreuses propriétés telles que : l’adhésion, le taux d’usure et le coefficient de frottement. Elles peuvent également être utilisées dans la fabrication des dispositifs médicaux tels que les stents. Bien que les alliages β-métastable à base de titane présentent des propriétés remarquables, une grande partie de la littérature publiée sur les films super-élastiques est consacrée aux alliages Ni-Ti. Ainsi, le but de ce travail est de développer de nouvelles couches super-élastiques de type β-métastable, obtenues par pulvérisation cathodique magnétron. Dans ce sens, deux types d’alliages ont été étudiés, les alliages binaires Ti-Nb et l’alliage quaternaire Ti-Nb-Ta-Zr. Après une description du dispositif expérimental utilisé pour la synthèse des couches, nous étudions plus particulièrement l’influence des paramètres de dépôt sur les propriétés des revêtements. La morphologie des revêtements a été caractérisée par microscope électronique à balayage. La diffraction de rayons X a été employée pour l’étude de la microstructure et la texture des films déposés. Enfin, l’effet super-élastique a été caractérisé par deux méthodes, l’indentation sphérique et la compression de micro-piliers usinés par faisceau d’ions focalisés / Superelastic thin films have a fairly wide field of applications. Thanks to their ability to accommodate a large deformation reversibly, superelastic films may be used as intermediate layers between the substrate and the hard layers. This significantly improves many properties such as adhesion, wear and coefficient of friction. They can also be used in the manufacture of medical devices such as stents. Although β-metastable titanium-based alloys have excellent properties, much of the published literature on the superelastic films is devoted to Ni-Ti. Thus, the purpose of this work is to develop new superelastic β-metastable titanium-based films, obtained by magnetron sputtering. After describing the experimental device used for the synthesis of thin films, we study more specifically the influence of deposition parameters on the properties of coatings. The morphology was characterized by the observation of fractured cross-sectional by scanning electron microscope. The X-ray diffraction was used to study the microstructure and texture of the deposited films. Finally, the superelastic effect was characterized by two methods, the spherical indentation and micro-pillars compression machined by focused ion beam

Dans un contexte d’augmentation constante des volumes de boues d’épuration à traiter, l’optimisation du traitement des boues est un enjeu primordial. Les étapes de traitement, et de transport mettent en jeu des écoulements qu’il est nécessaire de comprendre et de prédire afin, par exemple, de pour pouvoir estimer les pertes de charges en conduite ou bien pour dimensionner les installations de pompage. D’un point de vue physique, les boues peuvent être considérées comme une suspension de particules dans un gel suspendant. Ainsi, le comportement rhéologique des boues d’épuration présente des similitudes importantes avec les suspensions colloïdales et les gels polymériques. Ces trois types de matériaux, i.e. les boues d’épuration, les gels colloïdaux et les suspensions polymériques, présentent un comportement rhéologique complexe dépendant du temps et de la sollicitation imposée. Ils présentent un comportement dual, solide aux contraintes faibles, et liquide pour des contraintes élevées. La transition solide-liquide est généralement modélisée par la définition d’un seuil de contrainte ou de déformation, supposé séparer les régimes solide et liquide. Cependant, cette notion de seuil suppose une transition abrupte, et s’oppose aux observations expérimentales qui mettent en évidence une transition continue et progressive. L’étude de la littérature a permis de mettre en évidence une nécessité d’améliorer la compréhension et la modélisation du phénomène de transition solide-liquide. De plus, il est nécessaire d’unifier la description des régimes solide et liquide sous un même modèle, afin de mettre en lien une continuité mathématique avec le caractère continu et progressif du phénomène physique modélisé. Une analyse des résultats disponibles dans la littérature nous a permis de construire un modèle mathématique unique pour décrire le comportement solide et le comportement liquide des matériaux étudiés. Les hypothèses posées à partir de la littérature pour construire ce modèle ont ensuite été validées expérimentalement. Le modèle proposé est basé sur la décomposition de la complaisance du matériau en la somme d’une contribution solide et d’une contribution liquide, dépendant du temps, de la sollicitation appliquée et de l’histoire du matériau. Ce modèle permet une description commune des comportements solides et liquides du matériau, en tenant compte de l’existence d’une élasticité résiduelle y compris pour des contraintes élevées, et d’une dissipation visqueuse faible pour les contraintes faibles, conformément aux résultats expérimentaux. Ces travaux de thèse ont permis de mettre en évidence le fait que le mécanisme de transition solide-liquide était piloté non pas par la contrainte ou par la déformation, mais par la complaisance du matériau. De plus, ils ont permis d’ouvrir la voie à une nouvelle manière d’appréhender la thixotropie et la transition solide-liquide des matériaux pâteux. En effet, le comportement d’un matériau pâteux est piloté par deux paramètres : un module élastique plateau correspondant à un état totalement structuré, et une viscosité infinie correspondant à un état totalement déstructuré. Ces paramètres intrinsèques au matériau sont alors pondérés par des évolutions de la microstructure, menant à une compétition entre les effets élastiques et les effets visqueux. Ainsi, la différence entre un comportement de type loi de puissance et un comportement de type loi de puissance à seuil peut être expliquée simplement par l’apparition d’effets élastiques non négligeables. / In a context of constant increasing volumes of wastewater treatment sludge, optimizing the treatment of sludge appears to be crucial. Each step of treatment and transportation involves flows. It appears necessary to understand and predict these flows in order, for example, to estimate pressure drops in pipes or to size properly pumping facilities. In a physical point of view, sludge can be considered as a suspension of particles in a gel. Thus, its rheological behaviour presents significant similarities to that of colloidal suspensions of polymeric gels. These three types of materials, i.e. wastewater treatment sludge, colloidal suspensions and polymeric gels, present a complex rheological behaviour which depend on both time and the applied solicitation. They exhibit a dual behaviour, solid at low shear stresses, and liquid when the applied shear stress is high. The solid-liquid behaviour is generally modelled by defining a critical shear stress or a critical strain, supposed to be the limit between the solid and liquid regimes. Nevertheless, this concept implies an abrupt transition, unlike experimental observations showing a continuous and progressive transition. The study of the literature permitted to highlight the need to improve the understanding and modelling of the solid-liquid transition. Moreover, it appears necessary to unify the description of the solid and liquid regime in a unique model, in order to link a mathematical continuity with thecontinuous and progressive nature of the physical phenomenon to model. The study of the results available in the literature permited us to build a unique mathematical model to describe both the solid behaviour and the liquid behaviour of the studied materials. The assumptions made from the literature results have thus been experimentally validated. The proposed model is based on the decomposition of the compliance of the material in the sum of a solid contribution and a liquid contribution, depending on time, the applied solicitation and the story of the material. This model permits a unique description of solide and liquid regimes of the material, taking into account the existence of a residual elasticity at high shear stresses, and a viscous dissipation for low shear stresses, in accordance with experimental results. This work permitted to highlight the fact that the solid-liquid transition mecanism is controlled by the compliance of the material, and not the shear stress or the strain. Moreover, it opened the way to a new way of understanding the thixotropy and the solid-liquid transition of pasty materials. Thus, the behaviour of a pasty material is controlled by two parameters : a plateau elastic modulus corresponding to a totally structured state, and an infinite viscosity corresponding to a totally destructured state. These parameters intrinsic to the material are pondered by the evolutions of the microstructure, leading to a competition between elastic and viscous effects. Thus, the difference between the power law behaviour and the Herschel-Bulkley behaviour can be simply explained by the apparition of elastic effects that can’t be neglected.

La caractérisation des propriétés matérielles basée sur le dialogue essai-calcul requiert des précautions quant à la prise en compte des grandeurs physiques mises en jeu, aussi bien mesurées que calculées. Les mesures de champs cinématiques délivrent des informations qui permettent de rendre compte du comportement local du matériau et constituent une entrée riche pour ce type de dialogue axé sur la confrontation modèle - mesure. L’objet de ce travail est de discuter l’identification de propriétés élastiques homogènes ou hétérogènes, basée sur l’exploitation de mesures de champs, qui soulève la problématique de la prise en compte de mesures perturbées et d’hypothèses fortes de modélisation. Le but de cette étude est de discuter deux manières différentes pour aborder le problème. La première consiste à traiter les mesures en amont de l’étape d’identification, afin de filtrer les perturbations de mesures, et est basée sur l’outil d’approximation diffuse. La seconde consiste à développer une méthode d’identification inverse à partir d’essais quasi-statiques cherchant à tenir compte du caractère incertain des mesures, dès sa formulation, en utilisant l’ensemble des informations théoriques et expérimentales disponibles, en particulier en l’absence d’informations sur les conditions aux limites. Cette méthode est basée sur les principes de l’erreur en relation de comportement modifiée, amenant à résoudre des problèmes éléments finis non standards, et a conduit à mettre en œuvre des méthodes de résolution adaptées. L’étude a permis de souligner la robustesse de la méthode face aux incertitudes de mesure, notamment en comparaison avec le recalage par éléments finis. / The identification of material properties based on the confrontation between the experimental data and their numerical counterpart needs to pay special attention to their reliability and then to take some precautions. Especially, when dealing with full-field measurements which provide rich kinematic information enabling new possibilities for the characterization of material constitutive parameters. The aim of this work is to discuss the identification of the elastic parameters, for either homogeneous or heterogeneous materials, leading to the more general problem of the taking into account of noisy experimental data and strong model hypotheses. Two different approaches are addressed in this study. The first one consists in the processing of the experimental data upstream of the identification stage, by filtering them through a diffuse approximation algorithm. The second approach focuses on the development of an inverse identification strategy from static tests. This method takes into account the uncertainty of the measurements from the beginning of the formulation and is based on the modified constitutive relation error principles. It involves the resolution of non standard finite element problems using specific resolution methods. The study pointed out the robustness of the proposed method towards measurement errors, in particular when compared to the finite element model updating method.

Les bétons sont des matériaux composites à la microstructure complexe et constitués de phases dont le contraste des propriétés physiques et mécaniques peut être très grand. Ces matériaux posent des difficultés aux approches macroscopiques lorsqu'il s'agit de maîtriser leurs comportements effectifs comme celui du fluage. Malgré ces difficultés, EDF doit se doter d'outils permettant de modéliser de façon prédictive l'évolution des bétons des ouvrages en service ou de prescrire lecahier des charges des bétons de nouvelles installations. Ayant pour objectif de contribuer à la résolution de ce problème, ce travail de thèse développe des méthodes numériques multi échelle pour le calcul des structures constituées de matériaux fortement hétérogènes élastiques ou viscoélastiques. Plus précisément, ce travail de thèse comporte trois parties. Dans la première partie, nous nous intéressons à un composite constitué d'une matrice élastique renforcée par des inclusionsélastiques dont les formes géométriques peuvent être quelconques et dont la fraction volumique peut être importante. Pour modéliser ce matériau composite, une première approche numérique consistant à combiner la méthode des éléments finis étendus (XFEM) standard et la méthode "level-set" (LS) classique est d'abord utilisée. Nous montrons que cette première approche numérique, qui apparaît naturelle, induit en fait plusieurs artefacts numériques non rapportés dans la littérature, conduisant en particulier à une convergence non optimale par rapport à la finessedu maillage. Par suite, nous élaborons une nouvelle approche numérique ($mu $-XFEM) basée sur la description des interfaces par des courbes de niveaux multiples et sur un enrichissement augmenté permettant de prendre en compte plusieurs interfaces dans un même élément. Nous démontrons au travers des comparaisons et exemples que la convergence est améliorée de manière substantielle par rapport à la première approche numérique. Dans la deuxième partie, nous proposons une nouvelle méthode pour calculer les déformations différées des structures composées de matériaux hétérogènes viscoélastiques linéaires. Contrairement aux approches proposées jusqu'à présent, notre méthode opère directement dans l'espace temporel et permet d'extraire de manière séquentielle le comportement homogénéisé d'un matériau hétérogène viscoélastique linéaire. Concrètement, les composantes du tenseur de relaxation effectif du matériau sont d'abord obtenues à partir d'un volume élémentaire représentatif et échantillonnées au cours du temps. Une technique d'interpolation et un algorithme implicite permettent ensuite d'évaluer numériquement la réponse temporelle du matériau par le biais d'un produit de convolution. Les déformations différées des structures sont enfin calculées par la méthode des éléments finis classique. Différents tests sont effectués pour évaluer la qualité et l'efficacité de la méthode proposée, montrant que cette dernière permet d'avoir un gain en temps de l'ordre de plusieurs centaines par rapport aux approches de type éléments finis multiniveaux. La troisième partie est consacrée à l'étude de la structure de l'enceinte de confinement d'un réacteur nucléaire. Nous prenons en compte les quatre niveaux d'échelles associés à la pâte deciment, au mortier, au béton et à la structure en béton précontraint par des câbles en acier. La méthode numérique d'homogénéisation élaborée dans la seconde partie est appliquée afin de construire les lois de comportement pour chacun des trois premiers niveaux. Les résultats obtenus présentent un intérêt pratique pour résoudre des problèmes posés par EDF / Concretes are composite materials having complex microstructure and consisting of phases whose physical and mechanical properties can exhibit high contrast. Difficulties arise when macroscopic approaches are used to evaluate their effective behaviors such as the creeping one. Despite thesedifficulties, EDF has to be endowed with tools allowing to model in a predictive way the evolution of concrete structures in service or to prescribe the specifications of concrete for new facilities. Aimed to contribute to solving this problem, this thesis develops multi-scale numerical methods for the computation of structures made of highly heterogeneous elastic or viscoelastic materials. More precisely, this thesis comprises three parts. In the first part, we focus on a composite consisting of an elastic matrix reinforced by elastic inclusions which may be of any geometric shapes and whose volume fraction can be significant. To model this composite material, a first numerical approach which combines the standard extended finite element method (XFEM) and the classical level-set method (LSM) is used. We show that this numerical approach, which appears natural, leads in fact to several numerical artefacts having not been reported in the literature, giving rise in particular to non-optimal convergence with respect to the fineness of the mesh. In view of this, we develop a new numerical approach based on the description of the interfaces with multiple level sets and augmented enrichment so as to account for multiple interfaces in a single finite element. We demonstrate through examples and benchmarks that the proposed method significantly improves convergence compared to the first numerical approach. In the second part, we elaborate a new method to compute the creeping of structures formed of linearly viscoelastic heterogeneous materials. Unlike the approaches proposed up to now, our method operates directly in the time space and allows to sequentially extract the homogenized properties of a linearly viscoelastic heterogeneous material. Precisely, the components of the effective relaxation tensor of the material is first obtained from a representative volume element and sampled over time. An interpolation technique and an implicit algorithm are then used to numerically evaluate the time-dependent response of the material through a convolution product. The creeping of structures is finally calculated by the classical finite element method. Various tests are performed to assess the quality and effectiveness of the proposed method, showing that it can give a gain of time of the order of several hundreds compared to approaches such as multi-level finite element. The third part of the thesis is devoted to the study of the structure of containment of a nuclear reactor. We consider four scale levels associated with cement paste, mortar, concrete and a pre-stressed concrete structure with steel cables. The numerical method of homogenization developed in the second part is applied to construct the constitutive equations for each of the first three scale levels. The results thus obtained are useful for solving some practical problems posed by EDF