Academic literature on the topic 'Méthodes des éléments finis multi-Échelles'

L'objectif de ce travail est d'étendre un code industriel, résolvant les équations de Navier-Stokes stationnaires moyennées (RANS), à l'instationnaire puis à la simulation des grandes échelles (LES). Une approximation éléments finis en espace fondée sur la formulation Galerkin/Moindres Carrés est utilisée. L'intégration temporelle est assurée par un schéma explicite Runge Kutta à quatre pas. Les différentes formulaions spatiales et le schéma en temps sont décrits, analysés sur l'équation d'advection monodimensionnelle, puis validés sur des calculs d'amplification temporelle d'instabilités (couche de mélange et couche limite bidimensionnelles). Le modèle de Smagorinsky et une variante sélective sont implémentés. La validation des méthodes numériques et des différents modèles de sous-maille est effectuée sur des simulations de turbulence homogène isotrope et de couches de mélange temporelles tridimensionnelles à différents nombres de Reynolds. L'application finale est l'étude de l'effet d'un jet supersonique sur une couche de mélange plane confinée subsonique supersonique. Les résultats des différentes simulations des grandes échelles sont comparés à une expérience réalisée au CEAT de Poitiers ainsi qu'avec des simulations moyennées k-e.

En dynamique explicite, l’analyse des structures est souvent très coûteuse à cause du grand nombre de pas de temps nécessaire à la stabilité des modèles. Le but est ici de proposer une méthode originale basée sur les techniques de décomposition de domaine, afin de réduire les temps de calcul. Une méthode de gestion du contact unilatéral entre sous-domaines a tout d’abord été introduite (méthode de la peau artificielle). Celle-ci repose sur la superposition d’une interface de contact (par pénalité) et d’une interface de liaison ; elle renvoie ainsi le traitement du contact inter-domaine à un contact interne au domaine. Par ailleurs, une méthode de raffinement a été intégrée à l’approche multi-domaines afin d’enrichir localement et temporairement le maillage au cours du calcul. La méthode étant rendue particulièrement instable par l’utilisation du schéma explicite, l’accent a été mis sur la manière de rééquilibrer la solution avant la reprise des calculs sur la nouvelle discrétisation<br>In explicit dynamics, structural analysis is often very expensive because of the large number of time steps necessary to the stability of the models. The goal is here to propose an original method, based on the domain decomposition algorithms, in order to reduce the time of calculation. A method of handling unilateral contact between domains was first introduced. This method rests on the superposition of a interface contact (by penalty) and of a link interface; it thus returns the treatment of the inter-domain contact to an internal contact to the domain. In addition, a refinement method was integrated into the multi-domain approach in order to enrich locally and temporarily the mesh during calculation. The method being made particularly unstable by the use of the explicit scheme, the accent was put on the way of rebalancing the solution before the resumption of calculations on the new discretization

Le travail de cette thèse a porté sur la simulation numérique des matériaux multi-échelles. On considère des matériaux hétérogènes dont les propriétés physiques ou mécaniques (conductivité thermique, tenseur d’élasticité, ...) varient à une échelle petite par rapport à la taille du matériau. La thèse s'articule en deux parties qui correspondent à deux aspects différents des problèmes multi-échelles.Dans la première partie, on se place dans le cadre de l'homogénéisation aléatoire et on s’intéresse à une question plus fine que la caractérisation d'un comportement moyen : on cherche à étudier les fluctuations de la réponse. Plus généralement, nous visons à comprendre : (i) quels paramètres de la distribution des coefficients du matériau à l'échelle fine affectent la distribution de la réponse à l'échelle macroscopique, et (ii) s’il est possible d’estimer cette distribution sans utiliser une méthode type Monte-Carlo, très couteuse. Sur le plan théorique, nous avons considéré un matériau faiblement aléatoire (micro-structure périodique avec ajout d’une perturbation aléatoire petite). Nous avons montré qu’en utilisant le correcteur standard issu de la théorie de l’homogénéisation aléatoire, nous sommes capables de calculer un tenseur Q qui gouverne complètement les fluctuations de la réponse. Ce tenseur, défini par une formule explicite, permet d’estimer la fluctuation de la réponse sans résoudre le problème fin pour de nombreuses réalisations. Une stratégie d’approximation numérique de ce tenseur a ensuite été développée et testée numériquement dans des cas plus généraux.Dans la deuxième partie de la thèse, on considère un matériau hétérogène déterministe fixé où les hypothèses classiques d'homogénéisation (périodicité, ...) ne sont pas vérifiées. Les méthodes de résolution standard type Éléments Finis donnent de mauvaises approximations. Pour pallier cette difficulté, la Méthode des Éléments Finis Multi-échelles (MsFEM) a été introduite il y une vingtaine d'année. La méthode MsFEM se décompose en deux étapes : (i) créer un espace d'approximation grossier engendré par les solutions de problèmes locaux bien choisis; (ii) approximer la solution avec une approche de Galerkine peu couteuse sur l'espace construit dans (i). Dans cette deuxième partie, plusieurs taches ont été réalisées. Tout d'abord, une implémentation de plusieurs variantes MsFEM a été effectuée sous forme de template dans le logiciel de calcul Éléments Finis FreeFem++.Par ailleurs, plusieurs variantes des MsFEM pâtissent d'une erreur dite de résonance : lorsque la taille des hétérogénéités est proche de la taille du maillage grossier, la méthode devient très imprécise. Pour pallier ce problème, une méthode MsFEM enrichie a été développée : à la base MsFEM classique on rajoute des solutions de problèmes locaux ayant pour conditions aux limites des polynômes de haut degré. L'utilisation de polynômes nous permet d'obtenir une convergence de l'approche à des coûts de calcul raisonnables<br>This thesis is about the numerical approximation of multi-scale materials. We consider heterogeneous materials whose physical or mechanical (thermal conductivity, elasticity tensor, ...) vary on a small scale compared to the material length. This thesis is composed of two parts describing two different aspects of multi-scale problems.In the first part, we consider the stochastic homogenization framework. The aim here is to go beyond the identification of an effective behavior, and to characterize the fluctuations of the response. Generally speaking we strive to understand: (i) what parameters of the distribution of the material coefficient affect the distribution of the response and (ii) if it is possible to approximate this distribution without resorting to a costly Monte-Carlo method. On the theoretical standpoint, we consider a weakly random material (the micro-structure is periodic and presents some small random defects).We show that we are able to compute a tensor Q that governs completely the fluctuations of the response, thanks to the use of standard corrector functions from the stochastic homogenization theory. This tensor is defined by an explicit formula and allows us to estimate the fluctuation of the response without solving the fine problem for many realizations. A numerical approximation of this tensor has been proposed and numerical experiments have been performed in broader random frameworks to assess the effectiveness of the approach.In the second part, we consider a heterogeneous deterministic material where classical homogenization (periodicity, ...) assumptions are not satisfied. Standard methods such as Finite Elements give bad approximations. In order to solve this issue the Multi-scale Finite Element Method (MsFEM) can be used.This approach proceeds in two steps: (i) design a coarse approximation space spanned by solutions to well-chosen local problems; (ii) approximate the solution by an inexpensive Galerkin approach on the space designed in (i). On this topic, we first implemented the main variants of the MsFEM methods in the Finite Element software FreeFem++ on template form.Second, many MsFEM approaches suffer from resonance error: when the size of the heterogeneities is close to the coarse mesh size the accuracy decreases. In order to circumvent this issue, we designed an enriched MsFEM method: to the classical MsFEM basis, we add solutions to local problems with high degree polynomial boundary conditions. The use of polynomials allows us to obtain a converging approach for a limited computational cost

Le cœur d'un réacteur nucléaire est un milieu très hétérogène encombré de nombreux obstacles solides et les phénomènes thermohydrauliques à l'échelle macroscopique sont directement impactés par les phénomènes locaux. Toutefois les ressources informatiques actuelles ne suffisent pas à effectuer des simulations numériques directes d'un cœur complet avec la précision souhaitée. Cette thèse est consacré au développement de méthodes d'éléments finis multi-échelles (MsFEMs) pour simuler les écoulements incompressibles dans un milieu hétérogène avec un coût de calcul raisonnable. Les équations de Navier-Stokes sont approchées sur un maillage grossier par une méthode de Galerkin stabilisé, dans laquelle les fonctions de base sont solutions de problèmes locaux sur des maillages fins prenant précisément en compte la géométrie locale. Ces problèmes locaux sont définis par les équations de Stokes ou d'Oseen avec des conditions aux limites ou des termes sources appropriés. On propose plusieurs méthodes pour améliorer la précision des MsFEMs, en enrichissant l'espace des fonctions de base locales. Notamment, on propose des MsFEMs d'ordre élevée dans lesquelles ces conditions aux limites et termes sources sont choisis dans des espaces de polynômes dont on peut faire varier le degré. Les simulations numériques montrent que les MsFEMs d'ordre élevés améliorent significativement la précision de la solution. Une chaîne de simulation multi-échelle est construite pour simuler des écoulements dans des milieux hétérogènes de dimension deux et trois<br>The nuclear reactor core is a highly heterogeneous medium crowded with numerous solid obstacles and macroscopic thermohydraulic phenomena are directly affected by localized phenomena. However, modern computing resources are not powerful enough to carry out direct numerical simulations of the full core with the desired accuracy. This thesis is devoted to the development of Multiscale Finite Element Methods (MsFEMs) to simulate incompressible flows in heterogeneous media with reasonable computational costs. Navier-Stokes equations are approximated on the coarse mesh by a stabilized Galerkin method, where basis functions are solutions of local problems on fine meshes by taking precisely local geometries into account. Local problems are defined by Stokes or Oseen equations with appropriate boundary conditions and source terms. We propose several methods to improve the accuracy of MsFEMs, by enriching the approximation space of basis functions. In particular, we propose high-order MsFEMs where boundary conditions and source terms are chosen in spaces of polynomials whose degrees can vary. Numerical simulations show that high-order MsFEMs improve significantly the accuracy of the solution. A multiscale simulation chain is constructed to simulate successfully flows in two- and three-dimensional heterogeneous media

Mes travaux de recherche portent sur l'analyse numérique des intégrateurs géométriques et multi-échelles pour les équations différentielles déterministes ou stochastiques. Les modèles d'équations différentielles issus de la physique ou la chimie possèdent souvent une structure géométrique ou multi-échelles particulière (par exemple, les structures hamiltoniennes, les intégrales premières, les structures multi-échelles en temps ou en espace, les systèmes hautement oscillatoires), mais leur complexité est souvent telle qu'une solution satisfaisante est hors de portée en utilisant seulement des méthodes numériques standards à usage général. L'objectif est donc d'identifier les propriétés géométriques ou multi-échelles pertinentes de ces problèmes, et d'en tirer avantage pour concevoir et analyser de nouveaux intégrateurs efficaces, fiables et précis, reproduisant fidèlement le comportement qualitatif de la solution exacte des modèles considérés.

Les structures maçonnées sont largement connues en génie civil comme constituant d’une partie des bâtiments, mais également en tant que garnissages réfractaires dans des structures utilisées à hautes températures, par exemple en sidérurgie. Malheureusement, les outils actuels ne sont pas suffisamment puissants pour prédire le comportement de ces structures avec l’apparition de fissures et pour tenir compte du comportement non linéaire d’un des deux constituants (le mortier par exemple). Ce travail de thèse contribue à la modélisation multi-échelles des maçonneries classiques et des garnissages réfractaires avec un faible coût numérique grâce à la technique d’homogénéisation périodique. Les techniques de modélisation et de simulation du comportement des maçonneries sont présentées et développées. L’influence des lois d’interface entre briques et mortier, des paramètres géométriques et matériels, ainsi que de la densité des fissures sur le comportement effectif des maçonneries est étudié. Trois approches (une extension analytique de Cecchi et Tralli, une approche numérique et un modèle micromécanique) sont proposées pour la détermination du comportement effectif d’une cellule périodique dans le cas de maçonneries avec mortier viscoélastique microfissuré et briques saines élastiques ou rigides. Les résultats des calculs sur deux exemples de maçonneries (1D et 2D) ont confirmé que l’approche multi-échelle est une solution appropriée avec une grande capacité à exprimer le comportement des maçonneries viscoélastiques microfissurées. Ce travail, limité au cas sans propagation de fissures, peut être étendu aux mortiers à comportement viscoplastique<br>Masonry structures are widely used in civil engineering as part of buildings or in refractory linings of structures working at high temperatures, like in steel industry. Unfortunately, the present tools are not powerful enough to predict the behavior of these structures at their micro-cracked state and/or if one of their constituents behaves nonlinearly (e.g. the mortar). This research contributes to the multi-level modeling of classical masonries and refractory linings with a low numerical cost using basically the periodic homogenization technique. Modeling and simulation techniques of masonry behavior are presented and developped. The influence of interface law between bricks and mortar, of geometrical and material parameters, and of crack density on the effective masonry behavior is studied. Three approaches (analytical extension of Cecchi and Tralli, numerical approach and micromechanical modeling) were proposed to determine the effective behavior of a periodic masonry cell with micro-cracked viscoelastic mortar and safe elastic or rigid bricks. The results obtained on two examples of masonry (1D and 2D) confirmed that the multi-scale approach is a suitable solution with a great ability to model the effective behavior of microcracked viscoelastic masonry. This work, actually limited to the case without crack propagation, could be extended to mortars with viscoplastic behavior

Ce travail de thèse a pour but de modéliser la macroségrégation des produits obtenus par solidification en prenant en compte le transport des grains équiaxes. Le modèle de solidification à deux phases (solide et liquide) est traité par une méthode d'éléments finis, consistant à résoudre les équations de conservation moyennées de l'énergie, de la quantité de mouvement et de la masse, dans lesquelles les évolutions multi-échelles de la masse des phases et des solutés sont modélisées en utilisant une approche de splitting. D'après cette technique, la variation des quantités est considérée comme résultant de la contribution de deux étapes : l'étape de croissance et l'étape de transport. L'implémentation numérique du modèle a été réalisée avec trois opérations principales : tout d'abord implémenter le modèle de croissance des grains, ensuite intégrer des phénomènes de transport résultant de la convection thermo-solutale du liquide et du mouvement du solide, enfin mettre en œuvre le modèle complet en combinant les étapes de croissance et de transport. Lors de ces opérations, une investigation attentive a été consacrée à l'établissement de la résolution par éléments finis pour les équations de transport avec champs discontinus de vitesse à divergence non nulle, afin de surmonter des problèmes numériques en respectant la qualité des solutions physiques. Parallèlement à ces travaux, différents tests de simulation 2D ont été effectués à chaque étape d'implémentation. De bons accords ont été globalement obtenus entre les solutions données par le modèle présent et celles de référence dans la littérature. Finalement, des applications industrielles et des simulations 3D ont été menées, pour lesquelles les résultats numériques reproduisent les configurations caractéristiques des mesures expérimentales : un profil typique de macroségrégation des lingots composé d'une ségrégation négative occupant de la zone inférieure et d'une ségrégation positive en zone supérieure. Ce profil est la signature caractéristique de la sédimentation des cristaux et de la convection thermo-solutale. En outre, un modèle à trois phases étendu à partir du modèle à deux phases précédent – en distinguant la phase liquide interdendritique – s'est avéré capable de décrire la morphologie des grains dendritiques<br>The present work aims at modeling macrosegregation of castings, accounting for the transport of equiaxed grains. A two-phase (solid and liquid) finite element solidification model is presented, consisting in solving a system of volume-averaged conservation equations of energy, momentum, solute, in which the multi-scale evolutions of phase and solute mass are modeled by using a splitting method. According to this approach, the variation of quantities is considered as due to the contribution of two stages: the growth stage and the transport stage. The numerical implementation was realized with three principal steps: first implementing growth processes, then integrating transport phenomena including the thermo-solutal liquid convection and the solid movement, lastly combing the growth and transport stages to achieve a complete growth-transport model. Of these steps, solving the transport equations with discontinuous and non-divergence-free velocity fields by using finite element method required an attentive investigation in order to overcome numerical issues while respecting for physical solutions. Parallel to these works, various two-dimensional simulation tests were carried out in each implementation step. Agreements were globally found between results obtained from the present model and those of reference from the literature. Finally, industrial applications and three-dimensional simulations were performed, which show that computational solutions can predict essential features of experimental measurements. In particular, a typical macrosegregation profile of steel ingots, containing a negative segregation in the lower zone and a positive segregation in the upper zone, which is predominantly characterized by crystals sedimentation and fluid circulation was retrieved. Moreover, a three-phase model considering two different liquid phases, extended from the above-mentioned two-phase model, which enables to describe the morphology of dendritic solid crystals was implemented

Dans cette thèse, des techniques d’homogénéisation multi-échelles sont proposées pour l’analyse des vibrations des matériaux composites viscoélastiques. Dans la première partie, la Méthode Asymptotique à Deux Echelles (MADE) est proposée pour la modélisation des vibrations des longues structures sandwichs viscoélastiques répétitives. Pour ce type de structures les pulsations amorties correspondant aux modes propres de vibration sont regroupées en paquets bien distincts. La MADE décompose le problème initial de grande taille en deux problèmes de petites tailles. Le premier est défini sur quelques cellules de base et le second est une équation différentielle d’amplitude à coefficients complexes. La résolution de ces problèmes permet de déterminer les propriétés amortissantes correspondant aux modes de début et de fin de paquet de la structure tout en évitant la discrétisation de toute la structure. Pour les structures dont les coeurs ont un module d’Young dépendant de la fréquence, le problème non linéaire formulé sur les cellules de bases est résolu par l’approche diamant. Les modèles ADF et à dérivées fractionnaires ont été considérés dans les tests numériques. En utilisant la MADE, on évite la discrétisation de toute la structure, ce qui permet donc de réduire considérablement le temps de calcul ainsi que l’espace mémoire CPU nécessaires. L’approche proposée a été validée en comparant les résultats à ceux de la simulation éléments finis basée sur la discrétisation de toute la structure, et utilisant l’approche diamant. Dans la seconde partie de cette thèse, la méthode des éléments finis multi-échelles (EF2) a été développée pour le calcul des propriétés modales des structures à matériaux hétérogènes viscoélastiques en terme de fréquences amorties et amortissements modaux. Dans le principe de l’approche EF2, le problème de vibration est formulé à deux échelles : l’échelle de la structure globale (échelle macroscopique) et l’échelle d’un VER minutieusement choisi (échelle microscopique). Le problème à résoudre à l’échelle microscopique est un problème non linéaire alors que le problème à résoudre à l’échelle macroscopique est un problème linéaire. La non linéarité à l’échelle microscopique est introduite par la dépendance en fréquence du module d’Young des matériaux des phases viscoélastiques. Le problème non linéaire ainsi généré à l’échelle microscopique est résolu grâce à la MAN et ses outils de différentiation automatique réalisés sous Matlab, Fortran et C++. Un outil numérique, générique, robuste, peu coûteux en temps de calcul et espace mémoire CPU, de résolution des problèmes de vibrations non amorties des structures composites viscoélastique est ainsi mis en place. Le modèle viscoélastique à module constant ainsi que des modèles à modules dépendant de la fréquence notamment le modèle ADF et le modèle à dérivées fractionnaires ont été considérés pour les tests numériques de validation. Les comparaisons avec les résultats ABAQUS ont confirmé l’efficacité du code propos é. Le modèle est ensuite utilisé pour le calcul des propriétés amortissantes des structures sandwichs viscoélastiques à coeur composite. Les capacités de la nouvelle approche à concevoir des structures sandwichs viscoélastiques à coeur composite et à haut pouvoir amortissant ont été testées avec succès à travers l’étude de l’influence des différents paramètres des inclusions sur les propriétés amortissantes d’une structure sandwich viscoélastique à coeur composite<br>In this thesis, multiscale homogenization techniques are proposed for vibration analysis of structures with viscoelastic composite materials. In the first part, the Double Scale Asymptotic Method is proposed for vibration modeling of large repetitive viscoelastic sandwich structures. For this kind of structures, la eigenfrequencies are closely located in well separated packets. The DSAM splits the initial problem of large size into two problems of relatively small sizes. The first problem is posed on few basic cells, and the second one is an amplitude equation with complex coefficients. The resolution of these equations permits to compute the damping properties that correspond to the beginning and the end of every packets of eigenmodes. In case of structure with frequency dependent Young modulus in the core, the diamant approach is used to solve the nonlinear problem posed on basic cells. The ADF and fractional derivative models are considered in numerical tests. By using the DSAM, one avoid the discretization of the whole structure, and the computation time and needed CPU memory are thus reduced. The proposed method is validated by comparing its results with those of the direct finite element method using the diamant approach. In the second part of this thesis, the multiscale finite element method (FE2) is proposed for computation of modal properties (resonant frequency and modal loss factors) of structures with composite materials. In the principle of the (FE2) method, the vibration problem is formulated at two scales: the scale of the whole structure (macroscopic scale) and the scale of a Representative Volume Element (RVE) considered as the microscopic scale. The microscopic problem is a nonlinear one and the macroscopic problem is linear. The nonlinearity at the microscopic scale is introduced by the frequency dependence of the Young modulus of the viscoelastic phases. This nonlinear problem is solved by the Asymptotic Numerical Method and its automatic differentiation tools realizable in Matlab, Fortran or C++. From this approach, numerical tool that is generic, flexible, robust and inexpensive in term of CPU time and memory is proposed for vibration analysis of viscoelastic structures. The constant Young modulus and frequency dependent Young modulus are considered in validation tests. The results of numerical simulation with ABAQUS are used are reference. The model is then used to compute the modal properties of sandwich structure with viscoelastic composite core. To test the capacities of the proposed approach to design sandwich viscoelastic structure with high damping properties, the influence of parameters of the inclusions are studied

Lors de la conception des structures du génie civil et du génie mécanique il est très important de pouvoir prévoir comportement près de leur état ultime. Or, les matériaux utilisés sont souvent fortement hétérogènes, comme par exemple le béton et les métaux composites, et sont en conséquence difficiles à modéliser en cas de sollicitations extrêmes. Dans ce travail nous proposons une nouvelle stratégie multi échelles adaptée pour ce type des problèmes. Selon l'approche on couple un modèle basé sur la méthode des éléments finis (MEF) à l'échelle de la structure (macro) avec un modèle très fin à l'échelle de la microstructure (micro), lui aussi basé sur la MEF. Cela nous permet de modéliser l'échelle micro beaucoup plus fidèlement que ce qui est possible avec les approches phénoménologiques ou avec des méthodes d'homogénéisation analytique. L'inconvénient principal d'une telle approche étant en général le coût de calcul très élévé, nous considérons trois aspects rendant possible une application aux situations réelles. Premièrement, la formulation variationnelle permet de n'utiliser le modèle fin que dans les parties restreintes de la structure, là où les échelles sont vraiment fortement couplées, tout en gardant un modèle homogénéisé ailleurs. Deuxièment, nous adaptons la MEF sur l'échelle micro en employant l'approche dite structurée, avec laquelle on peut réduire le nombre de degrés de liberté d'une manière importante en gardant la même précision. Le troisième aspect porte sur l'implantation de la méthode dans un code de calcul, adaptée pour des machines parallèles, ce qui nous a permis de lancer les analyses d'une taille comparable à celle des situations réelles. Enfin, nous montrons de quelle manière cette approche peut être utilisée dans l'optimisation de forme des inclusions d'un matériau composite<br>In design of civil and mechanical engineering structures it is very important to be able to adequately predict their failure mechanisms. However, materials used are often strongly heterogeneous, like for instance the concrete or metalic composites, and are consequently difficult to model when submitted to extreme loading. In this work we propose a new multi scale strategy adapted for this kind of problems. According to the approach we couple a finite element method (FEM) model at the structural scale (micro) with a very fine model at the scale of the microstructure (micro), also based on the FEM. It allows us to model the micro scale much more accurately than what is possible with phenomenological approaches or analytical homogenisation methods. Since in general the major incovenience of such an approach is a very high computational cost, we consider three aspects which enable its application to realistic situations. First of all, the variational formulation allows for using the fine model only in parts of the structure, where the scales are really strongly coupled, whereas we use a homogenised model elsewhere. In addition, we adapt the FEM on the micro scale by employing the so called structured approach, with which we can significanlty reduce the number of degrees of freedom by retaining the came accuracy. The third aspect concerns the implementation of the method in a finite element code, adapted for parallel machines, with which we could run analyses of a size comparable to realistic situations. Finally, we show in which way this approach can bc used for inclusion shape optimisation of a composite material

Lors de la conception des structures du génie civil et du génie mécanique il est très important de pouvoir prévoir comportement près de leur état ultime. Or, les matériaux utilisés sont souvent fortement hétérogènes, comme par exemple le béton et les métaux composites, et sont en conséquence difficiles à modéliser en cas de sollicitations extrêmes. Dans ce travail nous proposons une nouvelle stratégie multi échelles adaptée pour ce type des problèmes. Selon l'approche on couple un modèle basé sur la méthode des éléments finis (MEF) à l'échelle de la structure (macro) avec un modèle très fin à l'échelle de la microstructure (micro), lui aussi basé sur la MEF. Cela nous permet de modéliser l'échelle micro beaucoup plus fidèlement que ce qui est possible avec les approches phénoménologiques ou avec des méthodes d'homogénéisation analytique. L'inconvénient principal d'une telle approche étant en général le coût de calcul très élévé, nous considérons trois aspects rendant possible une application aux situations réelles. Premièrement, la formulation variationnelle permet de n'utiliser le modèle fin que dans les parties restreintes de la structure, là où les échelles sont vraiment fortement couplées, tout en gardant un modèle homogénéisé ailleurs. Deuxièment, nous adaptons la MEF sur l'échelle micro en employant l'approche dite structurée, avec laquelle on peut réduire le nombre de degrés de liberté d'une manière importante en gardant la même précision. Le troisième aspect porte sur l'implantation de la méthode dans un code de calcul, adaptée pour des machines parallèles, ce qui nous a permis de lancer les analyses d'une taille comparable à celle des situations réelles. Enfin, nous montrons de quelle manière cette approche peut être utilisée dans l'optimisation de forme des inclusions d'un matériau composite.