This work is dedicated to the control of the accuracy of computational simulations of sound propagation and scattering. Assuming time-harmonic behaviour, the mathematical models are given as boundary value problems for the Helmholtz equation <i>Delta u+k2u=0 </i> in <i>Oméga</i>. A distinction is made between interior, exterior and coupled problems and this work focuses mainly on interior uncoupled problems for which the Helmholtz equation becomes singular at eigenfrequencies. As in other application fields, error control is an important issue in acoustic computations. It is clear that the numerical parameters (mesh size h and degree of approximation p) must be adapted to the physical parameter k. The well known ‘rule of the thumb’ for the h version with linear elements is to resolve the wavelength <i>lambda=2 pi k-1</i> by six elements characterising the approximability of the finite element mesh. If the numerical model is stable, the quality of the numerical solution is entirely controlled by the approximability of the finite element mesh. The situation is quite different in the presence of singularities. In that case, <i>stability</i> (or the lack thereof) is equally (sometimes more) important. In our application, the solutions are ‘rough’, i.e., highly oscillatory if the wavenumber is large. This is a singularity inherent to the differential operator rather than to the domain or the boundary conditions. This effect is called the <i>k-singularity</i>. Similarly, the discrete operator (“stiffness” matrix) becomes singular at eigenvalues of the discretised interior problem (or nearly singular at damped eigenvalues in solid-fluid interaction). This type of singularities is called the <i>lambda-singularities</i>. Both singularities are of global character. Without adaptive correction, their destabilizing effect generally leads to large error of the finite element results, even if the finite element mesh satisfies the ‘rule of the thumb’. The k- and lambda-singularities are first extensively demonstrated by numerical examples. Then, two <i>a posteriori</i> error estimators are developed and the numerical tests show that, due to these specific phenomena of dynamo-acoustic computations, <i>error control cannot, in general, be accomplished by just ‘transplanting’ methods that worked well in static computations</i>. However, for low wavenumbers, it is necessary to also control the influence of the geometric (reentrants corners) or physical (discontinuities of the boundary conditions) singularities. An <i>h</i>-adaptive version with refinements has been implemented. These tools have been applied to two industrial examples : the GLT, a bi-mode bus from Bombardier Eurorail, and the Vertigo, a sport car from Gillet Automobiles. As a conclusion, it is recommanded to replace the rule of the thumb by a criterion based on the control of the influence of the specific singularities of the Helmholtz operator. As this aim cannot be achieved by the <i>a posteriori</i> error estimators, it is suggested to minimize the influence of the singularities by modifying the formulation of the finite element method or by formulating a “meshless” method.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

L'objectif de ce mémoire de recherche est de déterminer les propriétés acoustiques des milieux poroélastiques à partir d'une démarche multi-échelle et multi-physique. Il traite d'échantillons réels de mousses, à cellules ouvertes ou partiellement fermées, dont les propriétés microstructurales sont caractérisées par des techniques d'imagerie. Cette information est utilisée afin d'identifier une cellule périodique idéalisée tridimensionnelle, qui soit représentative du comportement acoustique du milieu poreux réel. Les paramètres gouvernant les propriétés acoustiques du milieu sont obtenus en appliquant la méthode d'homogénéisation des structures périodiques. Dans une première étape, la structure des mousses est supposée indéformable. Il a été montré que pour le cas d'une distribution étroite de tailles caractéristiques de la géométrie locale, le comportement macroscopique d'une mousse à cellule ouverte peut être calculé à partir des propriétés géométriques locales de manière directe. Dans le cas d'une distribution étendue, le comportement acoustique du milieu est gouverné par des tailles critiques qui sont déterminées à partir de la porosité et de la perméabilité statique pour une mousse à cellules ouvertes ; pour une mousse à cellules partiellement fermées il est nécessaire d'identifier en plus une dimension connue de la géométrie locale. Nos résultats sont comparés avec succès à des données expérimentales obtenues par des mesures au tube d'impédance. Dans une seconde étape, les propriétés élastiques effectives du milieu poreux sont déterminées. Une modélisation par éléments finis de la cellule représentative a été mise en œuvre. Les paramètres élastiques calculés sont finalement comparés avec les données de la littérature, ainsi qu'à des essais mécaniques / This work aims at determining the acoustical properties of poro-elastic media through a multi-scale method. Some imaging techniques (tomography and micrographs) allow to estimate some quantitative microstructure properties of foams containing open or partially closed cells. These properties are used in order to clarify the features of a representative three-dimensional unit cell of a periodic structure, which mimics the behaviour of the real foam. All parameters controlling the acoustical properties of the porous foam are obtained by using the homogenization of periodic structures. In a first step, the structure of the foam is assumed to be rigid. It was shown that, in the case of a narrow distribution of the characteristic size of the local geometry, a direct computation of the macroscopic behaviour from the local geometrical properties is consistent with the measured acoustical properties. For a wide distribution of pore size, the acoustical behaviour is controlled by critical sizes that are obtained from porosity and static permeability for an open-cell foam, while for partially closed cells, the identification of a complementary characteristic dimension within the pores becomes necessary (e.g. closure rate of membranes). Our results compare well with data obtained from an impedance tube set-up. In a second step, effective elastic properties are computed through a modelling of the foam structure by finite elements. The computed elastic parameters are finally compared with data coming from the literature and with results of mechanical tests

Dans de nombreux domaines industriels et scientifiques, en particulier dans les domaines du génie civil et de génie mécanique, des matériaux à l’échelle de la microstructure, un très hétérogène par rapport à la nature du comportement mécanique. Cette fonctionnalité peut faire la prédiction du comportement de la structure soumise à différents types de chargement, nécessaires pour la conception durable, assez difficile. Le contrôle du comportement des ouvrages de génie civil est très complexe en raison de la diversité de la charge à laquelle ils sont soumis. La construction est maintenant réglementée partout dans le monde: les normes sont plus strictes et pris en compte, jusqu’à un état limite, en raison de différentes charges, par exemple des charges sévères tels que l’impact ou tremblement de terre. Modèles de comportement des matériaux et des structures doivent inclure l’élaboration de ces critères de conception et deviennent plus complexe. Ces modèles sont souvent basées sur des approches phénoménologiques, sont capables de reproduire la réponse du matériau au niveau ultime. Réponses de contrainte-déformation des matériaux sous sollicitations cycliques, dont de nombreuses recherches ont été exécutées dans les années précédentes afin de caractériser et le modèle, sont définies par différents types de propriétés de plasticité cycliques tels que l’écrouissageue, l’effet rochet et de de relaxation. En utilisant les modèles de comportement existants, ces réponses mentionnées peuvent être simulés d’une manière raisonnable. Cependant, il peut y avoir échec dans certains simulation des réponses structurelles et la déformation locale et globale. Insuffisance de ces études peut être résolu par le développement de solides modèles de comportement à l’aide d’expériences et de la connaissance des principes de fonctionnement des différents mécanismes de comportement inélastique ensemble. Dans ce travail, nous présentons un modèle phénoménologique constitutive qui est capable de coupler deux principaux mécanismes de comportement inélastique, plasticité et endommagement. Le modèle vise les applications de chargement cycliques. Ainsi, dans une partie de plasticité ou de dommages, les effets de durcissement isotropes et linéaires cinématiques à la fois sont pris en compte. Le principal avantage de ce modèle est l’utilisation de la plasticité indépendante contre les critères de l’endommagement pour décrire les mécanismes inélastiques. Un autre avantage concerne la mise en oeuvre numérique d’un tel modèle fourni en hybride-stress variationnel, obtenu avec une précision très améliorée et calcul efficace du stress et des variables internes dans chaque élément. Plusieurs exemples sont présentés afin de confirmer l’exactitude et l’efficacité de la formulation proposée en application à un chargement cyclique. / In many industrial and scientific domains, especially in civil engineering and mechanical engineering fields, materials that can be used on the microstructure scale, are highly heterogeneous by comparison to the nature of mechanical behavior. This feature can make the prediction of the behavior of the structure subjected to various loading types, necessary for sustainable design, difficult enough. The construction of civil engineering structures is regulated all over the world: the standards are more stringent and taken into account, up to a limit state, due to different loadings, for example severe loadings such as impact or earthquake. Behavior models of materials and structures must include the development of these design criteria and thereby become more complex, highly nonlinear. These models are often based on phenomenological approaches, are capable of reproducing the material response to the ultimate level. Stress-strain responses of materials under cyclic loading, for which many researches have been executed in the previous years in order to characterize and model, are defined by different kind of cyclic plasticity properties such as cyclic hardening, ratcheting and relaxation. By using the existing constitutive models, these mentioned responses can be simulated in a reasonable way. However, there may be failure in some simulation for the structural responses and local and global deformation. Inadequacy of these studies can be solved by developing strong constitutive models with the help of the experiments and the knowledge of the principles of working of different inelastic behavior mechanisms together. This dissertation develops a phenomenological constitutive model which is capable of coupling two basic inelastic behavior mechanisms, plasticity and damage by studying the cyclic inelastic features. In either plasticity or damage part, both isotropic and linear kinematic hardening effects are taken into account. The main advantage of the model is the use of independent plasticity versus damage criteria for describing the inelastic mechanisms. Another advantage concerns the numerical implementation of such model provided in hybrid-stress variational framework, resulting with much enhanced accuracy and efficient computation of stress and internal variables in each element. The model is assessed by simulating hysteresis loop shape, cyclic hardening, cyclic relaxation, and finally a series of ratcheting responses under uniaxial loading responses. Overall, this dissertation demonstrates a methodical and systematic development of a constitutive model for simulating a broad set of cycle responses. Several illustrative examples are presented in order to confirm the accuracy and efficiency of the proposed formulation in application to cyclic loading.

La thèse de doctorat propose la description physiquement cohérente de la limite d'élasticité du matériau sur la base des résultats des expériences associées à des simulations FEM. L'approche théorique est basée sur l'hypothèse de de la limite d'élasticité du matériau de Burzynski. L'hypothèse de l'effort de matériau est dérivée de la définition de l'énergie de fin de course variable de changement de volume et de distorsion. En cas particuliers, les critères comprennent asymétrie de gamme élastique et, sous une forme élargie, il est applicable aux matériaux initialement anisotropes. La thèse a démontré la nécessité de fournir une vue intégrée sur des expériences et des simulations numériques afin d'obtenir la description de l'apparition de la plasticité. Pour atteindre cet objectif, une méthodologie spécifique a été développé- la procédure qui combine des techniques expérimentales avec les simulations numériques pris en charge par la description précise du comportement des matériaux. Comme une partie de l'approche de modélisation numérique, les analyses numériques été effectuées afin de corriger les inexactitudes des résultats expérimentaux présentés. En outre, la Burzynski condition isotrope a été mis en oeuvre dans la méthode des éléments finis via UMAT (user subroutine for ABAQUS). La versatilité du rendement description est illustrée par l'application de matériaux à structure amorphe et différents systèmes de cristaux cubiques: la structure FCC (face centered cubic) et BCC (body centered cubic). A titre d'exemples: - un matériau de type FCC structure - OFHC Cu, Oxygen Free High Conductivity Copper, un matériau de type BCC structure - E335, acier non allié de haute résistance et l'exemple des matériaux amorphes - polycarbonate (PC) et biopolymères (PLA/PBAT). Pour les matériaux les analyses précises ont été effectuées afin d'obtenir des propriétés microstructurales et mécaniques via large gamme de tests de résistance dans diverses conditions de chargement. Le nombre de données exemplaires trouvés dans la littérature des matériaux sensibles à la pression et/ou initialement anisotrope ont été analysés. Il a été prouvé que les critères proposés donnent approximation précise de l'état du début de la plasticité par rapport à d'autres conditions d'usage courant. Les critères résultant semblent être une méthode universelle et largement applicable d'obtenir la description de l'apparition de la plasticité / This thesis proposes physically consistent description of the elastic limit of the material on the basis of the results associated with FEM simulations experiments. The theoretical approach is based on the assumption that the limit of elasticity of material Burzynski. The hypothesis effort material is derived from the definition of the end energy of variable stroke volume change or distortion. In special cases, the criteria include asymmetry elastic range and in an expanded form, it is applicable to initially anisotropic materials. The thesis has demonstrated the need to provide an integrated view on experiments and numerical simulations to obtain the description of the onset of plasticity. To achieve this goal, a specific methodology was developed - the procedure which combines experimental techniques with numerical simulations supported by a precise description of the behavior of materials. As part of the approach to numerical modeling, numerical analyzes were performed to correct the inaccuracies of the experimental results presented. In addition, Burzy?ski isotropically been implemented in the finite element method via UMAT (user subroutine for ABAQUS). The versatility of description performance is illustrated by the application of materials with amorphous structure and cubic crystals of different systems: structure FCC (face centered cubic) and BCC (body centered cubic). Examples : - type material FCC structure - OFHC Cu, Oxygen Free High Conductivity Copper material BCC structure type - E335 , not alloyed high strength and example of the amorphous steel materials - polycarbonate (PC) and biopolymers ( PLA / PBAT ) . Materials for specific analyzes were performed to obtain microstructural and mechanical properties via wide range of stress tests in various loading conditions. The number of copies of data found in the literature of the pressure sensitive and / or anisotropic materials was initially analyzed. It has been proven that the proposed criteria provide accurate approximation of the state of early plasticity in relation to other terms in common use. The resulting criteria appear to be universal and widely applicable method of obtaining a description of the appearance of the plasticity

Ce travail de thèse contribue à étudier le comportement élasto-plastique au rayage du verre métallique massif Zr55Cu30Al10Ni5 par la méthode des éléments finis (MEF) et expérimentaux. Le critère de plasticité de type Drucker-Prager est utilisé, une méthode de remaillage est proposée afin d'éviter le problème de convergence qui vient de la grand contrainte cisaillement. Pour cela, le comportement au rayage est étudié par le changement de l'angle d'attaque, du coefficient de frottement local entre le matériau et l'indenteur. Ces influences sont valorisés par l'évolution des forces rayures, par la morphologie de l'échantillon, et ils montrent le mécanisme de changement des modes de rayage : du labourage à l'usinage avec la formation des copeaux. Les essais sont réalisés à l'aide de l'équipement du LARMAUR : Un nano-triboindenteur de type Hysitron Ti-950. La comparaison des résultats simulés et expérimentaux nous permettre d'évaluer la loi comportement utilisé et d'estimer la valeur de coefficient de frottement entre le matériau et l'indenteur. / This thesis contributes to modeling of the elasto-plastic behavior in scratch test of a Zr55Cu30Al10Ni5 Bulk Metallic Glasses (BMG) by the finite element method and by experimental testing. The Drucker-Prager criterion is used and dedicated re-meshing method is proposed in order to solve the numerical problems classically encountered when modeling such a test and this with a controlled element population. The influence of friction coefficient as well as the angle of attack are investigated to understand the occurrence of two deformation mechanisms : plowing and cutting. The test is realized by using the LARMAUR's equipment : nano triboindenter Hysitron Ti-950. The result of experimental and simulation are confronted the finite element simulations.

Une technique numérique basée sur le couplage d’une approximation par éléments finis et d’une méthode intégrale a été développée pour le calcul du signal induit dans les antennes I.R.M. Ce signal est issu du mouvement de précession libre de l’aimantation transversale du corps à explorer pendant la relaxation. Dans notre modélisation, l’aimantation transversale représente le champ magnétique source. Celui-ci induit dans l’antenne un courant d’une durée très brève (quelques millisecondes) ; il représente le signal contenant toutes les informations de l’échantillon. Notre modélisation des antennes d’I.R.M de type circuit imprimé a été validée par comparaison avec des mesures expérimentales ainsi qu’avec une méthode analytique. Nous l’avons développée en tenant compte de leurs géométries et de leurs caractéristiques électromagnétiques afin d’avoir un meilleur rapport Signal/Bruit. Nous avons pris en considération des principaux facteurs tels que la distance entre l’antenne et l’échantillon à explorer ainsi que les caractéristiques électromagnétiques de l’antenne. / A numerical technique, based on the combination of a finite element method and a boundary integral method, has been developed to compute the induced signal in MRI antennas. This signal rises from a free movement of precession of the transverse magnetization of the sample to explore. In our modeling, the transverse magnetization represents the magnetic source field. Its flux embraces the antenna to give rise to a sinusoidal current which is very quickly attenuated in time (a few ms); it represents the signal containing all the information of the sample. We here want to find the geometrical and electromagnetic characteristics of the antennas which permit to have a signal to noise ratio as great as possible. In our computation, we have taken into account leading factors such as the distance between the probe and the organ to be explored and also the geometrical and electromagnetic characteristics of the probe. Our modeling of printed circuits MRI antenna has been validated by comparing with experimental measurements and also with an anlytical method. We have developped it by taking into account their geometries and their electromagnetical characteristics in order to have a better signal/noise ratio. We have considered principal factors such as the distance between the antenna and the organ to explore and also the electromagnetic characteristics of the antenna.

Cette thèse vise à développer une méthode d'homogénéisation numérique multi-échelle pour la prédiction des propriétés élastoviscoplastiques d'un béton à haute performance (BHP), où les hétérogénéités microstructurales sont analysées par nanoindentation. La procédure d'homogénéisation a été effectuée en deux niveaux en fonction de la microstructure du BHP : le niveau correspondant à la matrice composée de produits d'hydratation, de clinker et de cendre volante, et le niveau correspondant au mortier qui contient la porosité et les agrégats. Le comportement élasto-viscoplastique des phases microstructurales individuelles de la matrice est identifié à partir des données de nanoindentation. Une méthode inverse est réalisée en utilisant la simulation par éléments finis des tests de nanoindentaion combinée à une procédure d'optimisation numérique. Les résultats micromécaniques sont ensuite utilisés comme paramètres d'entrée pour l'homogénéisation élastoviscoplastique numérique à l'échelle microscopique. Le niveau correspondant au mortier est analysé par homogénéisation numérique en utilisant la simulation par éléments finis pour prédire les propriétés globales du béton. Les résultats sont comparés aux résultats expérimentaux et analytiques macroscopiques issus de la littérature montrant un bon accord. La procéduredéveloppée dans ce travail donne des résultats prometteurs pour la prédiction des propriétés élasto-viscoplastiques du BHP et peut être étendue à d'autres lois de comportement, y compris l’endommagement. La dernière partie de la thèse est consacrée à l'étude expérimentale de deux bétons préparés avec deux adjuvants différents. / This thesis aims at developing a numerical multiscale homogenization method for prediction of elasto-viscoplastic properties of a high performance concrete (HPC) where the microstructural heterogeneities are analyzed with the nanoindentation. The homogenization procedure was separated into two-levels according to the microstructure of the HPC: the matrix level composedof hydration products, clinker and fly-ash and the mortar level which accounts for large air porosity and aggregates. The elastoviscoplastic behavior of individual microstructural phases of the matrix are identified from nanoindentation data. An inverse method is applied by using the finite element simulation of the nanoindentaion tests combined with a numerical optimizationprocedure. The micromechanical results are then used as input parameters for numerical elasto-viscoplastic homogenization at microscale. The mortar level is analyzed with numerical homogenization by using the finite element simulation to predict the overall elasto-viscoplastic high performance concrete properties. The results are compared with macroscopic experimental and analytical results from the literature showing a good agreement. The procedure developed in this work gives promising results for the prediction of HPC properties and can be extended to other constitutive laws including damage. The last part of the thesis is dedicated to the experimental study of two concretes prepared with two different adjuvants.

Ces dernières années, l'étude de la rupture localisée des structures massives sous chargement thermomécanique est devenue un enjeu important en Génie Civil du fait de l'augmentation du nombre de constructions endommagées ou totalement effondrées après un feu. Deux questions centrales ont émergé : la modélisation mathématique des phénomènes mis en jeu lors d'un feu d'une part et la simulation numérique de ces problèmes d'autre part. Concernant la modélisation mathématique, la principale difficulté est la mise en place de modèles thermomécaniques capables de modéliser le couplage existant entre les effets thermiques et mécaniques , en particulier dans une zone de rupture localisée. Comment le chargement mécanique affecte la distribution de température dans le matériau et inversement, comment le chargement thermique influence la réponse mécanique ? Sont des questions qui doivent être abordées. Ces questions sont d'autant plus difficiles à aborder que l'on considère une zone de rupture où la mécanique des milieux continus classiques ne peut pas être appliquée du fait de la présence de discontinuités du champ de déplacement et, comme cela est démontré dans ce travail, du flux thermique. Pour ce qui concerne la simulation numérique, la complexité du problème multi-physique posé en termes de système d'équations aux dérivées partielles impose le développement de méthodes de résolution approchées adaptées, efficaces et robustes, la solution analytique n'étant en général pas disponible. Cette thèse contribue sur les deux aspects précédents. En particulier, des modèles thermomécaniques pour le béton et l'acier (les deux principaux matériaux utilisés en Génie Civil) capables de contrôler simultanément les phases d'écrouissage accompagnées de plasticité et/ou d'endommagement diffus ainsi que la phase adoucissante due au développement de macro-fissures sont proposés. Le problème thermomécanique est ensuite résolu par une méthode dite «adiabatic operator split» qui consiste à séparer le problème multiphysique en une partie mécanique et une partie thermique. Chaque partie est résolue séparément en utilisant une fois de plus une méthode «d'operator split» dans le cadre des méthodes à discontinuités fortes. Dans ces dernières, une discontinuité du champ de déplacement ou du flux thermique est introduite et gérée au niveau élémentaire du code de calcul Éléments Finis. Une procédure de condensation statique élémentaire permet de prendre en compte ces discontinuités sans modification de l'architecture globale du code de calcul Éléments Finis fournissant les champs de déplacement et de température. Dans cette thèse est également abordée la question de l'évaluation de la réponse jusqu'à rupture de structures en béton armé de type poteaux/poutres soumises à un feu. L'originalité de la formulation proposée est de tenir compte de la dégradation des propriétés mécaniques du matériau due au chargement thermique pour la détermination de la résistance limite et résiduelle des structures mais également de prendre en compte deux types de rupture caractéristiques des structures poteaux/poutres à savoir les ruptures en flexion et les ruptures en cisaillement. Les travaux présentés dans cette thèse pourront être étendus pour décrire la rupture de structures en béton armé dans des cas bi ou tridimensionnels en tenant compte en particulier du comportement de l'interface acier/béton et/ou d'autres types de rupture comme la rupture par fatigue ou le flambage. Une extension possible est également la prise en compte des effets dynamiques mis en jeu lorsque la structure est sollicitée mécaniquement par un tremblement de terre ou un impact en plus de la sollicitation thermique. / During the last decades, the localized failure of massive structures under thermo-mechanical loads becomes the main interest in civil engineering due to a number of construction damaged and collapsed due to fire accident. Two central questions were carried out concerning the theoretical aspect and the solution aspect of the problem. In the theoretical aspect, the central problem is to introduce a thermo-mechanical model capable of modeling the interaction between these two physical effects, especially in localized failure. Particularly, we have to find the answer to the question: how mechanical loading affect the temperature of the material and inversely, how thermal loading result in the mechanical response of the structure. This question becomes more difficult when considering the localized failure zone, where the classical continuum mechanics theory can not be applied due to the discontinuity in the displacement field and, as will be proved in this thesis, in the heat flow. In terms of solution aspect, as this multi-physical problem is mathematical represented by a differential system, it can not be solved by an ‘exact’ analytical solution and therefore, numerical approximation solution should be carried out. This thesis contributes in both two aspects. Particularly, thermomechanical models for both steel and concrete (the two most important materials in civil engineering), which capable of controling the hardening behavior due to plasticity and/or damage and also the softening behavior due to the localized failure, are carried out and discussed. Then, the thermomechanical problems are solved by ‘adiabatic’ operator split procedure, which ‘separates’ the multi-physical process into the ‘mechanical’ part and the ‘thermal’ part. Each part is solved individually by another operator split procedure in the frame-work of embbed-discontinuity finite element method. In which, the ‘local’ discontinuities of the displacement field and the heat flow is solved in the element level, for each element where localized failure is detected. Then, these discontinuities are brought into the ‘static condensation’ form of the overall equilibrium equation, which is used to solved the displacement field and the temperature field of the structure at the global level. The thesis also contributes to determine the ultimate response of a reinforced concrete frame submitted to fire loading. In which, we take into account not only the degradation of material properties due to temperature but also the thermal effect in identifying the total response of the structure. Moreover, in the proposed method, the shear failure is also considered along with the bending failure in forming the overal failure of the reinforced structure. The thesis can also be extended and completed to solve the behavior of reinforced concrete in 2D or 3D case considering the behavior bond interface or to take into account other type of failures in material such as fatigue or buckling. The proposed models can also be improved to determine the dynamic response of the structure when subjected to earthquake and/or impact.

Une capsule est une goutte liquide entourée d'une fine membrane élastique. La modélisation du comportement mécanique d'une capsule est un problème complexe de couplage entre la mécanique des fluides interne et externe et celle de la paroi de la capsule, qui peut subir des grandes déformations. On a développé une nouvelle méthode numérique pour modéliser le comportement d'une capsule dans un écoulement de Stokes infini. Le problème fluide est résolu à l'aide de la méthode des intégrales de frontière, tandis que la problème solide l'est par des éléments finis de membrane. La stratégie de couplage développée est non-classique, en particulier parce que les déplacements (plutôt que les efforts) sont imposés aux éléments finis de membrane. Cette nouvelle méthode est conditionnellement stable, précise et demeure stable en présence de compression dans le plan de la membrane, ce qui est une situation qui pose problème avec d'autres méthodes. On applique cette méthode à l'étude de capsules de forme initiale ellipsoïdale, dans un écoulement de cisaillement simple infini. Ces capsules ont des applications biomédicales intéressantes, et constituent un premier modèle mécanique de globule rouge. On met en évidence deux régimes distincts, selon l'intensité du cisaillement dans le fluide : le tumbling (quasi-solide) et le swinging (quasi-fluide). On montre également que les propriétés de la paroi (loi de comportement, module de dilatation surfacique) ont une influence cruciale sur le comportement de la capsule. La nouvelle méthode numérique est prometteuse et peut être étendue pour prendre en compte un rigidité de flexion de la paroi.