Ce travail de thèse porte essentiellement sur la détermination des propriétés mécanique et thermique des biocomposites HDPE-particules de bois et PET-fibres de chanvre, en utilisant des approches expérimentales, numériques et analytiques. La recherche des propriétés effectives de ces biocomposites prend en compte différents paramètres tels que la morphologie et l’orientation des fibres, d'une part, la porosité, et l’interphase, d'autre part. En effet, une étude basée sur la technique d'homogénéisation numérique a été réalisée en vue de vérifier l'influence des fibres courtes et logues sur les propriétés apparentes des différents matériaux étudiés. Aussi, des calculs numériques ont permis d'évaluer les propriétés élastiques ainsi que la conductivité thermique des biocomposites en fonction de la fraction volumique des fibres et des particules. L'objectif ultime de ce travail consiste à une modélisation de la mise en forme des biocomposites par thermoformage. En effet, une étude a été axée sur la thermoformabilité d’un coffrage en biocomposite HDPE-particules de bois, en utilisant des données expérimentales pour mieux décrire le comportement mécanique et thermique en vue d'une modélisation numérique des différentes étapes du thermoformage, en fonction de la teneur en particules de bois. / This thesis work concerns the determination of mechanical and thermal properties of HDPE-wood particles and PET-hemp fibers biocomposites using experimental, numerical and analytical approaches. The search for the effective properties of these biocomposites involves different parameters such as the fiber morphology and orientation on the one hand, and the porosity and the interphase on the other. A study based on the numerical homogenization technique has been carried out in order to verify the influence of short and long fibers on the apparent properties of the different materials. Numerical calculations have also allowed the estimation of the elastic properties as well as the thermal conductivity of the biocomposites in relation to the fiber and particle volume fraction. The ultimate objective of this work consists in the modelling of the thermoforming procedure on the biocomposites. A study on the thermoformability of the HDPE-wood particles biocomposite into a formwork in relation to the wood particle content has been carried out by using experimental data of the mechanical and thermal behavior for the numerical simulation of the thermoforming procedure.

Le design des structures en treillis est souvent contrôlé par des charges dynamiques comme le vent, les séismes ou des chocs dus à des ruptures soudaines comme les bris de conducteurs pour un pylône servant au transport d'énergie. Afin d'être en mesure de bien évaluer la réponse de ces structures dans le but d'optimiser leur conception, il faut être en mesure de bien définir les paramètres dynamiques tels que la rigidité et l'amortissement en fonction de la sollicitation. Il semble y avoir un manque de données dans la littérature pour définir ces paramètres. Afin de combler une partie du manque de données sur le sujet et de tester différentes techniques de modélisation en fonction de la sollicitation, une section d'un pylône électrique de 8 m de haut a été érigée au laboratoire et testée en régime d'oscillation libre résultant d'un lâché à plusieurs niveaux d'effort allant de 5 à 100% de la capacité théorique du pylône. Plusieurs modélisations numériques ont aussi été réalisées et comparées aux résultats expérimentaux. Steel lattice structure design is often controlled by dynamic loads such as wind and earthquake or shocks resulting from broken component such as guy cables or conductors. In order to properly evaluate the response of these structures under this type of loads, it is necessary to model accurately the structure. However, there is a lack of data in litterature defining those parameter. As well, the variation of stiffness and damping with the load level is not very well researched. In order to provide additional data and guidance on the modelling techniques to be used for lattice structures, and specifically to evaluate the relation of stiffness and damping with the load level, a 8 m long section of a transmission line tower was erected, pulled at different level of solicitation and brought to free vibration after the load was suddenly released. Numerical modelling were also conduct and compared to experimental results. Clear guidance are provided for practicing engineer interested in the modelling of such structures.

L’amélioration des sols par la technique de vibrosubstitution est habituellement utilisée pour les sols cohésifs, mais peut toutefois être utilisée dans les sols pulvérulents avec potentiel de liquéfaction. Cette technique permet la création d’un réseau de colonnes ballastées (CB) qui réduit le tassement et augmente la capacité portante d’un sol. L'intérêt est certes bien affiché. En effet, Geopac Inc. (entreprise spécialisée) se propose d'évaluer sa méthode d'amélioration des sols par CB en vue de réduire (voire éliminer) le potentiel de liquéfaction des sols grenus. Geopac a fait recours à une collaboration avec l'UdeS pour analyser expérimentalement et numériquement l'influence de la mise en place des CB sur les propriétés d'un sol pulvérulent. Un programme expérimental a été élaboré dans le but d’étudier l’effet de la mise en place des CB sur le sol encaissant par cette technique. À cet effet, un modèle réduit unique a été conçu, instrumenté et calibré à l’UdeS pour permettre dans un premier temps, de simuler l’incorporation de pierre nette (ballast) dans un sol lâche non saturé tout en maintenant une contrainte verticale et en mesurant la contrainte latérale générée durant l’incorporation. Le modèle réduit permet l’étude de l’expansion de cavité telle qu’obtenue durant la construction des CB par vibrosubstitution. Les colonnes sont construites à l’intérieur d’une cellule de PVC de 600 mm de haut et de 382 mm de diamètre. Un total de trois colonnes ballastées ont été construites sous une contrainte de 60, 80 et 100 kPa respectivement représentant un taux d’incorporation de près de 4 %. L’incorporation démontre une augmentation de la pression latérale caractérisée par le coefficient de pression latérale K jusqu’à 1.37 et une augmentation de la densité de 44 % à 55 % en moyenne. Une modélisation numérique a ensuite permis de confirmer les résultats expérimentaux et de vérifier le comportement du montage expérimental. La modélisation numérique a été réalisée avec le logiciel PLAXIS 2D 2012 en axisymétrie de révolution et considérant une loi comportementale de type Hardening Soil Model pour le sol encaissant et de type Mohr-Coulomb pour la colonne de pierre. La modélisation numérique a permis de déterminer que la forme de la contrainte verticale sur le sol est parabolique, c’est-à-dire maximale au centre et nulle à l’interface sol-paroi. L’expansion du sol a été simulée par l’imposition d’une déformation volumique anisotropique en cinq couches dans le sens radial et circonférentiel équivalent aux dimensions de la colonne obtenue expérimentalement. Les résultats numériques ont démontré une erreur de < 1 % pour la CB #1, de pour la CB #2 et de < 30 % pour la CB #3. L’erreur rencontrée sur la dernière colonne effectuée sous une contrainte verticale de 60 kPa peut-être associés à une mesure erronée lors de l'expérimentation. L’obtention des résultats expérimentaux par modélisation numérique a confirmé la fonctionnalité du montage expérimental et de la procédure de construction des CB en laboratoire. Cela ouvre la porte à des études subséquentes de l’effet de la mise en place des CB sur le sol encaissant en incorporant la vibration et la saturation du sable reproduisant la technique de vibrosubstitution.

L'objectif principal de cette thèse est la modélisation de matériaux supraconducteurs, et plus particulièrement, de ceux dits à haute température critique. Ce travail est divisé en deux parties distinctes : dans un premier temps, nous nous intéressons à la reconstruction de la distribution dynamique du courant dans la section d'un ruban supraconducteur à l'aide d'un modèle décrivant son comportement global. Ensuite, nous nous focalisons sur la modélisation électromagnétique et thermique par éléments finis des supraconducteurs en régime sur-critique, c'est à dire au-delà du courant critique Ic. Une introduction à la supraconductivité est faite dans le premier chapitre de ce travail, présentant les propriétés physiques principales que possède un matériau supraconducteur. Ensuite, nous décrirons les différents modèles existants permettant de décrire le comportement d'un échantillon supraconducteur. En se basant sur le modèle de Bean, nous proposons une nouvelle méthode permettant d'évaluer la distribution du courant dans la section d'un ruban monofilamentaire à partir de la mesure du profil magnétique à sa surface. Les résultats obtenus sont comparés avec ceux données par une modélisation certes plus précise, mais aussi plus complexe : l'utilisation de la méthode des éléments finis. Cette dernière sera largement utilisée dans la seconde partie de ce travail, qui traite de la modélisation de matériaux en régime sur-critique. Ceci nécessite d'introduire une formulation différente de l'habituelle loi de puissance Ec(J/Jc)^n pour décrire le comportement électrique du matériau dans une plage de courant plus importante, incluant les régimes sur-critiques. Pour cela, nous proposons une méthode originale, permettant d'approcher des mesures effectuées sur des échantillons d'YBCO et autorisant la description de la transition de l'état supraconducteur vers l'état normal. La dépendance en température a aussi été introduite afin d'étudier les phénomènes thermiques qui se produisent. Il est ainsi possible de résoudre de manière couplée les équations électromagnétiques et thermiques du problème : la résolution d'un pas de temps électromagnétique conduit à la connaissance des pertes locale. Ces pertes sont injectées dans le modèle thermique, et permet de calculer l'augmentation de température correspondante. Ce résultat est alors utilisé pour modifier les paramètres électriques en vue de la résolution du prochain pas de temps. Finalement, cette méthode a été utilisée afin de simuler le comportement d'un limiteur de courant supraconducteur. Le comportement global du système peut être reproduit par le modèle numérique implémenté, autorisant ainsi l'étude des variables locales, telles que la densité de courant ou le profil de température. Les résultats obtenus peuvent être utilisés afin d'optimiser les performances du dispositif selon des critères spécifiques. C'est ainsi que nous avons proposé une modification de la géométrie permettant d'éviter un possible emballement thermique qui pourrait conduire à la destruction du système. Sur la base de cette proposition, des modifications ont été apportées au design du dispositif, et sont actuellement testées.

La rupture des pentes et l’éclatement des roches, qui représentent deux types de risques naturels fréquents dans le monde, peuvent engendrer des conséquences économiques importantes et des pertes en vie humaine. Malgré que les phénomènes soient étudiés depuis de longues années, il reste encore des questions ouvertes et sans réponse et il est donc encore nécessaire poursuivre les recherches sur cette thématique. Le présent travail de thèse est consacré la modélisation numérique de la stabilité des grandes pentes et de l’éclatement des massifs rocheux en utilisant des méthodes basées sur l’intelligence artificielle en proposant des modifications et des améliorations de telles méthodes. En se basant sur des observations de déplacements de terrain, le glissement de terrain, qui est phénomène commun de la rupture de pentes, est étudié par le processus de Gauss afin de prédire son apparition temporelle. Ensuite, la question d’évaluation de la stabilité des pentes est abordée en utilisant la stratégie de machine à vecteurs de pertinence (RVM) avec des hyper-paramètres adaptatifs. Une approche itérative est proposée afin de déterminer les valeurs optimales des hyper-paramètres. Afin d’améliorer la prédiction, l’évaluation complète de la stabilité des pentes est réalisée en proposant un modèle basé sur la théorie de flou (CM) associé à un processus analytique d’hiérarchisation pondérée (WAHP). Ce modèle est utilisé à l’évaluation de la stabilité de la pente de rive gauche de la centrale hydroélectrique de Jinping 1, dans la région Sud-Ouest de Chine. Enfin, dans la dernière partie, la problématique de l’éclatement des massifs rocheux est abordée en utilisant des modèles basés sur la théorie du flou, en se basant sur une synthèse de 164 cas réels. Des comparaisons entre les résultats numériques et des données de terrain sont présentées pour de différents cas étudiés dans cette thèse. / Slope failures and rock burst, which are two typical types of geological hazards, create tremendous economic damages and cause massive losses to the health of domestic humans or animals every year throughout the world. The severe situation implies that they are still in need to be further studied despite the fact that they have been discussed for a long time. The present work is devoted to presenting the analysis of slope failures and rock burst using some computational intelligent models with modifications and improvements. Then landslide, a common type of slope failure, is analyzed for time occurrence prediction using the Gaussian Process by means of field-observed displacement series. After that, the problem of slope stability evaluation is discussed using the strategy of relevance vector machine (RVM) with adaptive hyper-parameter. An iteration approach is presented to find optimal hyper-parameter values in this chapter. Afterwards, the comprehensive evaluation of slope stability is carried out with the cloud model (CM) and weighted analytical hierarchy process (WAHP) closely related to the left abutment slope of Jinping 1 Hydropower Station, southwest of China. Finally, prediction of rock burst classification is engaged using the cloud models synthesized with the attribution weights on the basis of 164 rock burst cases. In each modeling of the associated problems, comparisons are given on the performance of each strategy as well as some evaluations.

Les phénomènes d’instabilité générés par les séismes dans les pentes représentent un risque naturel majeur à l’échelle mondiale. Ce risque a, de plus, tendance à croître du fait d’une urbanisation croissante dans des zones à fort aléa, comme en témoignent de nombreux évènements catastrophique à travers le passé. Pour garantir la sécurité des personnes et des biens face à ce type d’aléa, il faut améliorer les prédictions des déplacements générés par les séismes dans les pentes.Les premiers travaux scientifiques visant à corréler les caractéristiques des instabilités dans les pentes aux paramètres sismiques sont de nature empirique et remontent aux années 1980. Ces méthodes ne permettent pas d’expliquer tous les mouvements de terrain observés à travers le monde, notamment lorsque les effets de site modifient sensiblement la distribution des mouvements dans les pentes et génèrent des déplacements plus grands ou plus petits que ceux prédits par les lois empiriques qui négligent les effets de site.La première partie de cette thèse présente une nouvelle base de données mondiale de 277 glissements de terrain d’origine sismique ou non. L’analyse statistique conduite sur les données de cette base a montré que la forme géométrique moyenne des glissements de terrain reste stable lorsque le volume des instabilités augmente.La deuxième partie de la thèse est dédiée à l’analyse de la stabilité du glissement de terrain de Diezma (Espagne) au moyen de trois méthodes : la Méthode de NEWMARK (1965) basée sur le principe de l’équilibre limite en 2D qui suppose des conditions elasto-parfaitement-plastiques ; l’analyse par différences finies dans le domaine temporel avec le code FLAC en 2D et en 3D qui suppose des conditions visco-elasto-parfaitement-plastiques ; et l’analyse par recombinaison modale dans le domaine fréquentiel avec le code CESAR-LCPC en 2D et en 3D qui suppose des conditions visco-élastiques.Situé dans une région sismiquement active, le glissement de terrain de Diezma est susceptible d’être affecté par des séismes. Une analyse des déplacements induits dans ce versant par onze signaux sismiques différents a été menée dans l’objectif d’établir des corrélations entre les paramètres macro-sismiques des scenarios sismiques appliqués et les déformations calculées par les trois méthodes.Les résultats montrent des différences marquées à la fois qualitatives et quantitatives. Les conclusions principales sont : (a) Les déplacements obtenus par les trois méthodes ont des ordres de grandeur différents et peuvent donc conduire à une surestimation ou à une sous-estimation des déplacements ; (b) Les deux méthodes numériques montrent que les déplacements sont fortement conditionnés par la géométrie des modèles ; seule la méthode par différences finies fait apparaître une réponse en termes de déplacements qui dépend du scénario sismique considéré; (c) Les résultats des simulations2D et 3D ne sont pas comparables ; des analyses complémentaires doivent encore être menées pour guider l’utilisateur dans le choix de la méthode la plus appropriée; (d) La Méthode de NEWMARK (1965), dont l’utilisation reste très répandue de nos jours, est tout à fait appropriée à l’étude de cas simples mais elle peut se révéler inexacte lorsque la structure géologique / topographique du versant conduit à un fort effet de site car ce dernier n’est pas pris en compte par cette méthode / Seismically induced slope deformation is a worldwide common phenomenon that poses an increasing and considerable threat to fast expanding urbanization, and a great number of catastrophic events throughout the past attest thereof. For this reason, displacement predictions allowing for proper slope surveillance became a major concern.Early attempts to relate slope failures to seismic parameters are of empirical nature and date back to the 1980s. Although having proven stable, these relations are frequently disturbed by site effects causing outliers in terms of smaller or greater displacements than expected.The first part of this thesis presents a newly build chronological database of 277 globally distributed seismically and non-seismically induced landslides. A comprehensive statistical analysis was conducted on the data of this database with the main result that – statistically seen – the average geometrical shapes of landslides differ only proportionally.The second part of the thesis is dedicated to a comparative slope stability analysis of the Diezma Landslide (Spain) by means of three methods: the limit-equilibrium based NEWMARK-Method (1965) in 2D under elasto-perfectly-plastic conditions; finite-difference analysis in the time domain with the code FLAC in 2D and 3D under visco-elasto-perfectly-plastic conditions; and modal recombination analysis in the frequency domain with the finite-element code CESAR-LCPC in 2D and 3D under visco-elastic conditions.Located in a seismically active region, the Diezma Landslide is likely to be affected by earthquakes and was therefore considered as a suitable model case. A broad analysis of expected displacements was conducted using eleven strong-motion seismic scenarios. Moreover, the thesis searches for potential relations between macro-seismic parameters of the applied earthquake scenarios and the predicted deformations obtained from the three methods.It appeared that – due to the functionality of the methods – results differ quantitatively as well as qualitatively, and so does their suitability. Major findings are: (a) Results from the three methods are of different orders of magnitude and, thus, can easily lead to over- or under-estimations of displacements; (b) Both numerical methods reveal a strong influence of the model-geometry on the predicted displacements, whereas a scenario-dependent slope behavior manifested itself only within the finite difference analysis; (c) The switch from 2D to 3D does not necessarily result in a similar performance in each dimension and results must be critically judged before further use; (d) The NEWMARK-Method (1965) has proven itself once more to be appropriate for first slope assessments but not for sophisticated evaluations of ground response to seismic shaking

Les matériaux supraconducteurs présentent des propriétés physiques et géométriques particulières qui exigent des approches de modélisation spatio-temporelle fines, où les méthodes classiques trouvent rapidement leurs limites en termes de convergence, de précision et de temps de calcul. Ce dernier peut être très conséquent, ce qui est incompatible avec les problèmes de dimensionnement et d’optimisation. Dans ce contexte, ce travail a pour objectif de développer des approches de modélisation multiphysique rapides pour le dimensionnement et l’optimisation des systèmes à base de supraconducteurs. Un intérêt particulier est porté pour les méthodes intégrales. Les verrous scientifiques à lever, qui constituent également l’originalité du travail, résident dans l’intégration des lois de comportement E(J) des supraconducteurs dans les schémas numériques de ce type de méthodes. Dans ce travail on développe un modèle numérique afin d’étudier la distribution de la densité de courant et d’estimer les pertes AC dans les supraconducteurs à haute température (HTS). Le modèle développé est basé sur une formulation intégro-différentielle en termes de potentiel vecteur électrique dans les deux domaines fréquentiel et temporel. Une campagne de test est menée afin de valider et de bien cerner les possibilités offertes et les limites de cette approche pour la modélisation des supraconducteurs. / Superconducting materials have particular physical and geometric properties that require spatial-temporal modeling approaches fines, where conventional methods quickly reach their limits in terms of convergence, precision and computational time. The latter can be very consistent, which is incompatible with the design and optimization problems. In this context, this work aims to develop rapid multiphysics modeling approaches for the design and optimization of superconductor-based systems. Particular attention is paid to the integral methods. Scientific obstacles to overcome, which also constitute the originality of the work lies in the integration of behavior laws E (J) of superconducting digital patterns of such methods. A numerical model is developed for a rapid computation of eddy currents in multifilamentary high temperature superconductive (HTS) for the evaluation of AC losses. The developed model is based on an integro-differential formulation in terms of the electric vector potential in the frequency and temporal domains. A test campaign is conducted to validate and clearly identify the possibilities and limitations of this approach for modeling superconductors.

Ce travail porte sur la mise en oeuvre en 3D de la méthode des éléments naturels contrainte CNEM en vue de son utilisation pour la simulation du cisaillage à grande vitesse. La CNEM est une approche à mi-chemin des approches sans maillage et des éléments finis. La construction de son interpolation utilise le diagramme de Voronoï contraint (dual du maillage de Delaunay contraint) associé à un nuage de noeuds répartis sur le domaine étudié muni d'une description de sa frontière. La mise en oeuvre de la CNEM comporte trois aspects principaux : i) la construction du diagramme de Voronoï contraint, ii) le calcul des fonctions de forme éléments naturels Sibson, iii) la discrétisation d'une formulation variationnelle générique par utilisation de l'intégration nodale stabilisée conforme, SCNI, introduite par Chen et Al en 2001. Une partie importante de ce travail concerne les deux derniers points. Pour le calcul des fonctions de formes Sibson 3D cinq algorithmes sont présentés, dont deux développés au cours de la thèse, et sont comparés en terme de performance. Par ailleurs, une discrétisation est proposée pour être applicable au cas des domaines fortement non convexes. La mise en oeuvre proposée est validée sur des exemples en élasticité linéaire 3D en petites perturbations (vis à vis de solutions analytiques et de résultats éléments finis) puis en grandes transformations (test de la barre de Taylor). L'application de la CNEM au cisaillage grande vitesse est finalement abordée. Les développements effectués ont été intégrés à la plateforme logicielle Nessy. Cette plateforme a pour objectif la capitalisation du savoir faire du LMSP en simulation numérique.

La géophysique est une science fondamentalement liée aux observations, et son évolution rapide dans de nombreux domaines au cours des deux dernières décennies est due en grande partie à des avancées métrologiques importantes qui améliorent la qualité et la quantité des observations recueillies. Par exemple, dans le domaine de la sismologie, le réseau global (Global Seismographic Network - GSN) contient plus de 120 sismomètres, le premier instrument permanent de fond de mer (ocean-bottom seismometer - H2O) a été installé, et dans un futur proche des projets ambitieux tels que le réseau ultra dense américain de sismomètres ‘USArray' permettront d'obtenir des quantités très importantes de données à l'échelle d'un continent ou d'une région, enregistrées par des réseaux haute résolution de stations large bande. Dans le domaine de la géodésie satellitaire, que l'on peut aussi considérer en un sens comme la sismologie ultra longue période, les progrès métrologiques de positionnement GPS et l'interférométrie radar ont permis de mesurer les déformations et d'en déduire les contraintes à l'échelle régionale pendant et après de grands séismes (par exemple en Californie du Sud, Landers en 1992 ou Hector Mine en 1999), et d'étudier l'évolution de ces régions liée aux phénomènes post-sismiques. Bien évidemment, les progrès réalisés en instrumentation et acquisition de données doivent être accompagnés d'une meilleure compréhension théorique des phénomènes mis en jeu, ainsi que d'une meilleure modélisation de ceux-ci, au moyen de techniques de calcul numérique tridimensionnelles. Si nous comparons à l'évolution d'autres domaines tels que la mécanique des fluides, la météorologie ou encore l'astrophysique, dans lesquels les progrès en calcul scientifique ont suivi de près les progrès en instrumentation (cas des souffleries numériques, étude numérique de phénomènes climatiques tels que les ouragans etc.), le domaine de la sismologie n'a pas encore bénéficié pleinement des progrès rapides de l'analyse numérique, du calcul scientifique et de la technologie des ordinateurs modernes (calcul parallèle et vectoriel). Le présent mémoire résume les recherches que j'ai effectuées jusqu'à présent pour faire le lien entre différents domaines de la géophysique et les progrès du calcul scientifique. Nous avons d'une part développé des outils sophistiqués de calcul, nous les avons validés de manière rigoureuse, et nous les avons ensuite appliqués à des cas concrets intéressant la sismologie locale, régionale et globale. Dans la conclusion de ce document, nous exposons également quelques thèmes de recherche que nous souhaitons développer dans le futur.

La modélisation numérique des écoulements de glace est indispensable pour prédire l'évolution des calottes polaires suite au réchauffement climatique. De récentes études ont souligné l'importance des modèles d'écoulement dits d'ordre supérieur voir même de Stokes au lieu de la traditionnelle approximation de couche mince dont les hypothèses ne sont pas valables dans certaines zones critiques mais à l'étendue limitée. Cependant, ces modèles d'ordre supérieur sont difficiles à utiliser à l'échelle d'un continent en raison de leurs temps de calculs prohibitifs. Ce travail de thèse propose une nouvelle technique qui permet de réduire les temps de calculs tout en maximisant la précision des modèles. Plusieurs modèles d'écoulement de glace de complexité variables ont été mis en place dans ISSM (Ice Sheet System Model), un code élément fini massivement parallèle développé par le Jet Propulsion Laboratory. L'analyse et la comparaison des différents modèles, à la fois sur des cas théoriques et réels, montrent que l'utilisation des modéles les plus complets est principalement nécessaire au voisinage de la zone d'échouage, transition entre les parties flottantes et posées de la glace, mais aussi que des modèles plus simples peuvent être utilisés sur la majeure partie des glaciers. Coupler différents modèles présente donc un avantage significatif en terme de temps de calcul mais aussi d'amélioration de la physique utilisées dans les modèles. Plusieurs méthodes de couplage de modèles existent et sont présentées dans ce manuscrit. Une nouvelle technique, dite de tuilage, particulièrement adaptée au couplage de modèles d'écoulement de glace est décrite ici : son principe repose sur la superposition et le raccordement de plusieurs modèles mécaniques. Une analyse mathématique est effectuée afin de définir les conditions d'utilisation de cette méthode de tuilage. Le traitement du couplage entre un modèle de Stokes et des modèles simplifiés, pour lesquels le calcul des vitesses horizontales et verticales est découplé, est ensuite présenté. Cette technique a été mise en place dans ISSM afin de pouvoir créer des modèles hybrides combinant plusieurs modèles d'écoulement de complexité variable. Après avoir été validée sur des cas synthétiques, cette technique est utilisée sur des glaciers réels comme Pine Island Glacier, dans l'Antarctique de l'Ouest, afin d'illustrer sa pertinence. Les modèles hybrides ont le potentiel d'améliorer la précision des résultats en combinant différents modèles mécaniques, utilisés chacun dans les zones où leurs approximations sont valides, tout en réduisant les temps de calcul et en étant compatibles avec les ressources informatiques actuelles.