Le présent travail vise à apporter certaines pistes de solution concernant certaines controverses sur l'estimation de la conductivité thermique des nanotubes de carbone par simulation de dynamique moléculaire à l'équilibre avec conditions aux limites périodiques et la formule de Green-Kubo. Entre autre, différents auteurs obtiennent des résultats pouvant parfois varier de plusieurs ordres de grandeur pour un même type de nanotube. H n'y a toutefois que très peu d'études jusqu'à ce jour tentant d'expliquer ces contradictions. Dans la première partie du projet, on détermine les paramètres numériques pouvant influencer la conductivité thermique calculée avec une méthode de dynamique moléculaire à l'équilibre. On effectue ensuite une analyse de sensibilité pour plusieurs de ces paramètres afin de déterminer de quelle manière ils influencent la conductivité thermique calculée (chapitres 3 et 4). Finalement, on présente une étude sur le phénomène de fréquence de coupure lors du calcul de la conductivité thermique (chapitre 5).

Les systèmes complexes interviennent dans de nombreux secteurs d'activité du domaine industriel tels que la mécanique, l'électromécanique, la thermique, l'aérospatiale, la robotique, la biotechnologie, etc., L'asservissement de ces systèmes fait une partie importante des préoccupations actuelles des grandes unités de recherches en automatique et automatisation des procédés industriels à l'échelle internationale. Dans cette étude, une plate forme allant de la modélisation à la base de la commande automatique en temps réel a été développée et mise en ouvre pour ces systèmes. Une nouvelle approche de modélisation numérique des systèmes complexes en vue de leur commande en tempe réel a été proposée et formulée. En se basant sur l'exploitation des Fonctions Opératoires (FOs), le but de cette approche est d'adapter la structure et la paramétrisation du Modèle Interne (MI) ou de connaissance au processus de commande en temps réel. Ces fonctions ont été conçues pour qu'elles deviennent des outils efficaces pour l'estimation d'état et la prédiction des paramètres de la commande. Autrement dit, elles contiennent les informations nécessaires sur le système (caractère non linéaire, changement de phase, interaction interne, etc.) et sur le processus de commande (paramètres de la commande, perturbation, etc.). Un algorithme original d'identification et d'identification en ligne des modèles paramétriques multivariables dans des espaces multidimensionnels a été conçu et incorporé dans la plate-forme. Les contributions majeures de ce nouvel algorithme sont l'identification paramétrique, l'optimisation non linéaire, l'optimisation sous contraintes et l'incertitude paramétrique. Ses performances principales sont la configuration du réseau d'acquisition des données, la minimisation du temps de calcul et la grande précision. Cet algorithme peut être efficacement utilisé pour l'identification et l'ajustement des modèles multivariables avec un grand nombre de paramètres, les modèles de grandes incertitudes etc. Une nouvelle approche appelée Familles des Solutions Dynamiques (FSDs) a été développée comme base de la commande automatique en tempe réel des systèmes complexes. En se basant sur l'utilisation du MI identifié, nous avons conçu une nouvelle technique pour le développement des familles de trajectoires (cibles), des FOs (modèles) et des paramètres de la commande. Une banque mble de données cibles a été créée autour de l'état réel du System asservi pour une estimation d'état en ligne et une future prédiction de la commande automatique en temps réel. Ces familles sont caractérisées par leur stabilité exponentielle (asymptotique), convergence spatiotemporelle et la haute précision. Les FSDs constituent une base complète pour la commande à paramètres répartis et font une approche prometteuse dans le futur de l'automatique et l'automatisation des systèmes industriels de haute performance. / The complex Systems are implied in several modem industrial sectors such as mechanical engineering, mechatronic, aerospace, robotic, thermal engineering, chemical and nuclear engineering, biotechnology, etc. The real time control of these Systems is making the current interests of the advanced international research teams in the automatic and automation of industrial processes. In this thesis, a complete network of the numerical modeling, identification and online identification for a basis of real time automatic control of complex Systems was developed, implemented and tested. A new approach of numerical modeling for the real time control of complex Systems was proposed and formulated. Basing on the use of Operating functions (OFs) and the measures, the aim of this approach is to adapt the structure and the parameterization of the mathematical model or Internal Model (IM) to the real time control process. The OFs were conceived in order to become effective tools of the state estimation and the control prediction. In other words, these functions contain necessary information about the controlled System (nonlinear character, internal interaction, etc.) and the control process (control parameters, disturbances, etc). An original algorithm of identification and online identification of the multivariable models in multidimensional spaces was conceived and incorporated in the network. The major contributions of this algorithm are the online identification, nonlinear optimization, optimization under constraints and parametrical uncertainty. Its principal performances are the fast convergence, the minimal Computing time and the high precision. This algorithm is effectively used for the identification and adjustment of multivariable models, models with a great number of parameters and the models with a high uncertainty. A new basis called the Dynamic Families of Solutions (DFSs) was developed for the real time automatic control of the complex Systems. Using the identified IM, we have conceived a new technique for the development of complete families of tracking trajectories (targets or references) with the corresponding OFs (models) and the necessary control parameters (inputs). These families are characterized by the exponential (asymptotic) stability, fast convergence and high precision. The DFSs were developed as a complete basis for the real time automatic control of complex Systems and they constitute a promising approach in the automatic and automation of the industrial complex Systems such as the distributed parameter control.

La théorie des milieux poreux réactifs fermés est explorée pour modéliser le comportement au jeune âge des bétons pour des couplages chimio-poro-visco-élastiques. Le modèle est construit autour de l'hypothèse que la dépression capillaire est le principal mécanisme responsable du retrait d'autodessiccation (Hua et al. 1995). Ce retrait correspond à la réponse viscoélastique de la matrice solide vieillissante du milieu poreux sous l'action d'une contrainte de compression induite par la condensation capillaire. L'apparition d'une déformation de gonflement dans les premières heures qui suivent la prise du matériau ne peut être expliquée par un phénomène de dessiccation, car l'humidité relative décroît au fur et à mesure de l'évolution de la réaction d'hydratation. Les déformations endogènes (déformations libres du matériau en conditions isothermes, sans chargement extérieur et sans échange hydrique avec l'extérieur) sont donc définies comme la somme d'un gonflement chimique macroscopique dû à la réaction d'hydratation et d'un retrait d'autodessiccation. L'approche chimio-poro-visco-élastique développée dans cette étude pour analyser le comportement au jeune âge des bétons a plusieurs objectifs. Premièrement, une modélisation basée sur une approche de la mécanique des milieux poreux permet d'avoir une base rigoureuse à la compréhension des mécanismes qui ont une influence sur les déformations différées observées à l'échelle macroscopique. Le deuxième objectif consiste en la mise en place d'un outil qui fournit les bases nécessaires à une meilleure compréhension de l'influence de la composition du matériau sur le comportement au jeune âge des bétons. Ce rapport est divisé en quatre parties dont les deux premières font état des connaissances publiées dans la littérature scientifique qui permettront d'arriver à une modélisation basée sur une approche multi-échelles du problème. Les conditions de retrait restreint sont développées dans la troisième partie afin de montrer comment un essai R2D2 peut être analysé en vue d'une modélisation du comportement au jeune âge. La quatrième décrit les bases théoriques du modèle ainsi que l'identification des paramètres de ce dernier. Le rapport finit par une application de l'approche pour étudier l'influence du ciment et du rapport E/C sur le comportement au jeune âge des bétons.

L’objectif de ce projet de recherche est de développer un modèle d’évaporation de l’eau d’une rivière en milieu boréal. Un site expérimental sur la rivière Natashquan dans la région de la Côte-Nord au Québec a été sélectionné pour y installer une station hydrométéorologique faisant un suivi des variables environnementales locales. Cette rivière est marquante dans le paysage de la Minganie par la superficie de 15 930 km2 de son bassin versant et son débit moyen annuel de 349 m3/s. Ce puissant cours d’eau est une des dernières rivières naturelles de la Côte-Nord n’étant pas harnachée par un barrage. Une campagne constituée de 5 sorties de mesures intensives sur le terrain a été effectuée afin d’observer l’évaporation de l’eau de la rivière par la méthode des mini-lysimètres flottants. Ces derniers s’inspirent du processus derrière le bac évaporatoire de classe A pour obtenir une estimation de l’évaporation in-situ en ayant un petit bac évaporatoire portatif à la surface de l’eau de la rivière. Un modèle d’évaporation par transfert de masse a été développé avec les données récoltées et 8 fonctions de vent de la littérature ont été testées. Puis, un modèle d’évaporation par ensemble de réseaux de neurones a été retenu, car il offrait une meilleure performance que les approches classiques. Ce réseau de neurones ayant comme variables d’entrée la tension de vapeur à saturation, le rayonnement incident d’ondes longues, la vitesse du vent et la tension de vapeur de l’air, permet d’explorer une non linéarité entre les variables qui n’est pas accessible à l’approche de transfert de masse. Ce modèle met en évidence le manque de recherche sur les processus gouvernant l’évaporation fluviale dans le cadre de rivière de grande envergure et particulièrement l’évaporation des rivières du milieu boréal. / The objective of this study was to develop a river evaporation model for the boreal biome. A hydrometeorological station was installed on the bank of the Natashquan river to collect environmental data. The Natashquan river is major hydrographic component of the region with its 15 930 km2 watershed and its 349 m3/s mean annual flow. This powerful watercourse is one of the only remaining significant river that has not yet been dammed in the Côte-Nord region. Five intensive field campaigns were carried out targeting different conditions during the summers of 2018 and 2019. River water evaporation was observed using the floating minipan method. This method, inspired by the process behind the class A evaporation pan, allows small portable evaporation pans to be deployed on the water surface and to monitor in-situ evaporation in more holistic conditions. Observations were used to calibrate the wind function of a mass-transfer evaporation model and compare it with different wind functions taken from the literature. A stacked neural network was selected as the river model based on a higher performance compared to the mass transfer approach. The input variables of the network are the saturation vapour pressure, the incoming longwave radiation, the wind speed and the vapour pressure of the above air. The model allows a nonlinear combination of the selected environmental variables that is not accessible to the mass transfer equations. Therefore, the present study highlights the evident lack of research concerning fluvial evaporation in the context of rivers of considerable size, particularly in the boreal biome.

La modélisation de la cryolite, utilisée dans la fabrication de l’aluminium, implique plusieurs défis, notament la présence de discontinuités dans la solution et l’inclusion de la difference de densité entre les phases solide et liquide. Pour surmonter ces défis, plusieurs éléments novateurs ont été développés dans cette thèse. En premier lieu, le problème du changement de phase, communément appelé problème de Stefan, a été résolu en deux dimensions en utilisant la méthode des éléments finis étendue. Une formulation utilisant un multiplicateur de Lagrange stable spécialement développée et une interpolation enrichie a été utilisée pour imposer la température de fusion à l’interface. La vitesse de l’interface est déterminée par le saut dans le flux de chaleur à travers l’interface et a été calculée en utilisant la solution du multiplicateur de Lagrange. En second lieu, les effets convectifs ont été inclus par la résolution des équations de Stokes dans la phase liquide en utilisant la méthode des éléments finis étendue aussi. Troisièmement, le changement de densité entre les phases solide et liquide, généralement négligé dans la littérature, a été pris en compte par l’ajout d’une condition aux limites de vitesse non nulle à l’interface solide-liquide pour respecter la conservation de la masse dans le système. Des problèmes analytiques et numériques ont été résolus pour valider les divers composants du modèle et le système d’équations couplés. Les solutions aux problèmes numériques ont été comparées aux solutions obtenues avec l’algorithme de déplacement de maillage de Comsol. Ces comparaisons démontrent que le modèle par éléments finis étendue reproduit correctement le problème de changement phase avec densités variables. / The modelling of the cryolite bath, used in the smelting of aluminum, offers multiple challenges, particularly the presence of discontinuities in the solution and a difference in density between the solid and liquid phases. To over come these challenges, several novel elements were developed in this thesis. First of all, the phase change problem, commonly named the Stefan problem, was solved in two dimensions using the extended finite element method. A specially designed Lagrange multiplier formulation, using an enriched Lagrange multiplier solution, was implemented to impose the melting temperature on the interface. The interface velocity is determined by the jump in the heat flux across the interface and was calculated using the Lagrange multiplier values. Secondly, convection was included by solving the Stokes equations in the liquid phase using the extended finite element method as well. Thirdly, the density change between solid and liquid phases, usually neglected in the literature, was taken into account by the addition of a non-zero velocity boundary condition at the solid-liquid interface to maintain mass conservation in the system. Benchmark analytical and numerical problems were solved to validated the various components of the model and the coupled system of equations. The solutions to the numerical problems were compared to the solutions obtained using Comsol’s moving mesh algorithm. Theses comparisons show that the extended finite element method correctly solves the phase change problem with non-constant densities.

De nos jours, la plate-forme de Gough-Stewart est utilisée couramment en industrie. L'inconvénient avec ce type de manipulateur, comme c'est le cas pour tous les manipulateurs parallèles, est que l'analyse de ses lieux de singularité est fastidieuse et n'est pas intuitive. Une expression scalaire décrivant ces lieux a déjà été développée, mais elle comporte 20 coefficients complexes, d'où l'intérêt de la récrire en conservant les termes sous forme vectorielle. Après quelques manipulations effectuées sur la matrice jacobienne, on obtient une expression plus compacte qui peut ensuite être utilisée dans des analyses par intervalles. Ces analyses permettent d'étudier les lieux de singularité pour un espace fini et confirment l'intérêt d'une expression vectorielle. De plus, les facteurs qui influencent la précision des résultats sont énoncés. Dans ce mémoire, l'architecture dite générale et celle de type MSSM sont étudiées afin de démontrer à quel point l'expression des lieux de singularité peut être simplifiée.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2007-2008. / L’étude de l’interférence entre un plasmon polariton de surface (SPP) avec le faisceau plan lumineux incident sur une nano-fente métallique de dimension sous-longueur d’onde est devenue récemment un sujet fondamental de recherche dans le domaine de la plasmonique afin de mieux comprendre le mécanisme du rehaussement de la transmission optique de la nano-fente. Par la méthode de différences-finies dans le domaine de temps (FDTD), nous avons étudié le mécanisme de couplage du SPP dans la fente nanométrique. L’objectif de ce projet consiste dans un premier temps à la conception et l’implémentation d’un algorithme FDTD flexible et robuste, capable de simuler l’interaction de sources électromagnétiques avec des nanostructures de métal. L’algorithme développé sur la plate-forme Matlab permet de modéliser dans un espace bi-dimensionel des structures diélectriques dispersives. Dans un deuxième temps, nous avons employé ce simulateur FDTD pour étudier les mécanismes de couplage entre un SPP et une fente nanométrique percée dans une mince couche d’argent. Notre analyse démontre que le SPP incident à la fente est diffusé par le rebord de la fente et que les charges électriques induites sur les rebords ré-irradient l’énergie électromagnétique à l’intérieur de la fente. De plus, ces charges électriques génèrent des nouveaux SPPs sur les parois de la fente, qui contribuent à la formation de modes Fabry-Pérot le long de l’axe central de la fente. Nous démontrons aussi la formation de modes Fabry-Pérot créés par la multi-réflexion des ondes diffusées entre les parois de la fente. La combinaison de ces deux modes Fabry-Pérot produit une distribution de champ asymétrique dans la fente. Nous démontrons que la phase du SPP, relativement au faisceau incident normal, détermine les conditions d’interférence constructive et destructive correspondant respectivement au rehaussement et à la diminution de la transmission à travers la nano-fente. Finalement, nous avons confirmé la théorie d’interférence entre les champs induits par le SPP incident et le faisceau incident normal par l’addition de leurs amplitudes instantanées. / The study of the interference of the surface plasmon polariton (SPP) with the incident plane wave on a subwavelength metallic slit has become recently a fundamental subject of research in the domain of plasmonics. One of the objectives is to better understand the phenomenon of enhanced optical transmission through the nano-slit. Using the numerical method of finite-difference time-domain (FDTD), we have investigated the coupling mechanisms of the SPP inside the nano-slit. The objective of this project first consists in the conception and implementation of a flexible and robust FDTD algorithm capable of simulating the interaction between electromagnetic sources and metallic nanostructures. The algorithm that we have developed on the Matlab platform is able to model two-dimensional dispersive dielectric structures. The second and main objective is to use this FDTD simulator to investigate the coupling mechanisms between the SPP and a nano-slit pierced into a thin silver film. Our analysis demonstrates that the incident SPP is scattered by the inner edge of the slit and that the scattered waves induce oscillating electric charges on the slit edges which reradiate electromagnetic energy inside the slit. Moreover, these electric charges generate new SPPs on the slit walls which contribute to the formation of Fabry-Pérot modes along the central slit axis. We also show the formation of Fabry-Pérot modes set up by the multi-reflections of scattered waves between the slit walls. The combination of these two Fabry-Pérot modes produces an asymmetric field distribution inside the slit. We demonstrate that the phase of the SPP relatively to the normal incident beam determine the conditions of constructive and destructive interference which correspond, respectively, to the enhancement and suppression of the optical transmission through the nano-slit. Finally, we have confirmed the theory of interference between the induced fields by the SPP and incident normal beam via the superposition of their instantaneous fields.

Ce travail de thèse s’est concentré sur la modélisation de l’effet d’un couvert de neige saisonnier sur le comportement thermique de grands ouvrages de génie civil construits dans les régions nordiques, ainsi qu’à son application sur des ouvrages construits avec des matériaux à grande porométrie. À l’heure actuelle, la présence de neige est souvent négligée dans les simulations numériques de régimes thérmiques, sauf pour les études portant sur les avalanches ou le régime nival qui utilisent des modèles conceptuels ou des modèles très détaillés, basés sur la physique. L’utilisation des modèles très détaillés dans l’industrie n’est pas avantageuse pour modéliser de grands ouvrages de génie civil en raison du grand nombre d'équations différentielles partielles et de paramètres parfois difficiles à estimer. De plus, l’introduction de la neige comme un milieu dans le domaine de calcul impose de grands défis pour les simulations en continu sur plusieurs années, car, comme la neige est saisonnière, le domaine qui représente la neige doit être éliminé après sa fonte. À cet égard, l’objectif principal de cette recherche est de développer un nouvel outil pour modéliser, en continu sur plusieurs années, l’effet thermique d’un couvert de neige sur des ouvrages de génie civil, qui puisse être aisément appliqué dans des cas pratiques. Pour cela, les transferts de chaleur et de masse dans la neige ont été étudiés afin d'identifier les modes de transfert qui ont une influence significative sur la température du sol. Ensuite, un outil numérique dans le logiciel FlexPDE a été établi pour modéliser le transfert de chaleur par conduction et convection, entre le sol, la neige et l'atmosphère en continu sur toute l'année. Cet outil considère également l’effet de la neige (frontière fermée ou ouverte) sur la convection dans des matériaux granulaires grossiers, et l’effet de la pluie et du changement de phase de l’eau sur le bilan d’énergie du couvert de neige. Cet outil emploie la méthode du bilan d’énergie de surface qui est considéré comme étant une condition aux limites de Neumann pour la température. L’outil a également servi à réaliser une analyse thermique d’un barrage en remblai et à démontrer la présence des cellules de convection dans l’enrochement et l’influence de la convection d’air sur l’extraction de chaleur de la fondation. Un autre objectif de cette recherche est d'établir un modèle simple et précis de la conductivité thermique pour tous les types de neige, y compris la neige artificiellement manipulée et compactée. La plupart des modèles existants de la conductivité thermique de la neige sont développés grâce à des techniques de régression qui ont l'inconvénient de ne pas respecter les limites physiques de la neige. Pour intégrer ces limites physiques dans un modèle simple, le concept de conductivité thermique relative est utilisé dans cette étude. Ce modèle est examiné avec des données publiées et réévaluées et avec les résultats de tests effectués au laboratoire de l'Université Laval. Le modèle proposé permet d’estimer la conductivité thermique de tous les types de neige, et ce, avec une grande fiabilité. / This work focused on the modeling of the effect of a seasonal snow cover on the thermal behavior of large engineering structures built in the northern regions, and its application for structures built with materials of large porometry. Nowadays, the presence of snow is often neglected in thermal numerical simulations, except for investigating the phenomenon of avalanches or production of water after snowmelt in the mountainous regions, while using conceptual models or very detailed models based on physics. In the industry, modeling a large structure by detailed models have two principals drawback. First, because of large numbers of partial differential equations and parameters, which are sometimes difficult to estimate. The introduction of snow as a medium in the computational domain also imposes great challenges for the continuous simulations for consecutive years, because seasonal snow is present only a few months a year and the domain representing the snow must be removed after melting. In this regard, the main objective of this research is to develop a new tool for the modeling of snow thermal effect on geotechnical structure continuously for several years that can also be applied in simple practical cases. For this, transfers of the heat and mass in snow are first studied to identify the most important transfer modes that significantly affect soil temperature. Then, a numerical tool using the FlexPDE software has been established to model the heat transfer by conduction and convection between soil, snow and atmosphere continuously over the entire year. This tool also considers the effect of snow (close or open boundary) on the air convection in the coarse granular materials, the rain effect and the water phase change in the energy balance of the snow cover. This tool uses the surface energy balance as the Neumann boundary condition for temperature. The tool is also served for thermal analysis of an embankment dam and demonstrates the presence of convection cells and the influence of the air convection on the heat extraction of foundation. Another purposes of this research is to establish a simple and accurate model of thermal conductivity for all types of snow, including artificially manipulated and compacted snow. Most of the existing model of the snow’s thermal conductivity are developed through regression techniques which have the drawback of not respecting the physical limits of snow. To integrate these physical limits in a simple model, the relative thermal conductivity concept is used in this study. This model is verified with published data and further validated with the results of the tests performed in the laboratory of Laval University. The proposed model estimates the thermal conductivity of all types of snow with great reliability.

Cette thèse présente le développement d’une méthode d’ingénierie pour la modélisation de la consolidation au dégel non linéaire à grandes déformations. Les travaux qui y sont présentés s’inscrivent à la suite des théories et de modèles de consolidation existants. Les bases de la théorie de consolidation au dégel à une dimension sont donc reprises. En premier lieu, un modèle numérique de consolidation au dégel non linéaire à grandes déformations est formulé en combinant la théorie de consolidation à grandes déformations de Gibson et des équations de transferts de chaleur. Ces deux composantes sont couplées au sein d’un domaine de modélisation défini en coordonnées lagrangiennes qui s’adaptent aux déformations de sol. Il en résulte donc l’introduction d’une seconde frontière mobile à la surface pour la modélisation de la consolidation au dégel en plus de la frontière mobile au front de dégel. Le modèle utilise des relations non linéaires σʹ_v – e – k_v pour définir les propriétés des sols dégelés. Une étude de cas du pipeline expérimentale d’Inuvik est réalisée avec le modèle numérique. Cette étude de cas permet de démontrer l’utilisation du modèle pour un problème pratique et de valider le modèle. Les résultats obtenus sont comparés aux résultats de la théorie de consolidation au dégel linéaire à petites déformations et aux mesures de terrain. Le modèle non linéaire à grandes déformations offre une augmentation importante de la précision lors du calcul des tassements de fonte, du taux de tassement et des pressions interstitielles excédentaires. En second lieu, un modèle conceptuel pour la consolidation au dégel des sols dégelés à grains fins est proposé. Ce modèle permet de définir les caractéristiques des relations σʹ_v – e – k_v utilisées comme intrants au modèle numérique de consolidation au dégel. Le concept de la contrainte résiduelle y est généralisé aux sols riches en glace en spécifiant qu’il s’agit de la contrainte effective au sein des éléments de sol plutôt que de la contrainte effective globale du sol. Ensuite, des relations empiriques sont formulées pour déterminer les caractéristiques des relations σʹ_v – e – k_v à partir de l’indice des vides initial du sol dégelé et de propriétés d’indice du sol. / This thesis presents the development of an engineering method for the modelling of large strain nonlinear thaw consolidation. The work presented herein follows existing thaw consolidation theories and models. The foundations of one-dimensional thaw consolidation theory are therefore used. First, a numerical model for large strain nonlinear thaw consolidation is formulated by combining the Gibson large strain consolidation theory to heat transfer equations. The two components are coupled in a modelling domain formulated in Lagrangian coordinates that adapts to the soil deformation. This results in the introduction of a second moving boundary at the soil surface to model thaw consolidation in addition to the moving boundary at the thaw front. The model uses nonlinear σʹ_v – e – k_v relationships to define the properties of thawed soils. A case study of the Inuvik experimental pipeline with the numerical model is presented. The case study demonstrates the use of the model for a practical problem and it is used to validate the model. The modelling results are compared to the results obtained with the small strain linear thaw consolidation theory and with the field data. The results obtained with the large strain nonlinear model for thaw settlement, the rate of thaw settlement and the excess pore pressures compare favourably with the field data. Second, a conceptual model for thaw consolidation of thawed fine-grained soils is proposed. The model is used to define the characteristics of the nonlinear σʹ_v – e – k_v relationships used as input for the numerical modelling of thaw consolidation. The concept of the residual stress is generalized to ice rich soils by specifying that it is the effective stress within the soil element rather than the effective stress of the bulk soil. Then, empirical relationships are formulated to determine the characteristics of the σʹ_v – e – k_v relationships as a function of the initial thawed void ratio and soil index properties.

La géométrie et la flexibilité d'un canyon sont les paramètres qui affectent grandement la valeur des périodes naturelles dans les barrages en terre. Le canyon entourant des barrages peut être considéré comme un domaine illimité. Pour prendre en compte ces deux effets, le canyon a été modélisé par SBFEM et le barrage en terre, à géométrie limitée, par FEM. La technique hybride SBFEM-FEM pour l'analyse tridimensionnelle dynamique de l'interaction sol-barrage a été validée avec les résultats disponibles dans la littérature. Comme la matrice de rigidité dynamique du domaine non borné est complexe et dépendante de la fréquence, la méthode classique de superposition de modes n'est pas simple pour le système d'interaction sol-structure. Ainsi, pour obtenir la fréquence propre fondamentale, le barrage a été excité en direction amont-aval. Les périodes naturelles du barrage de terre pour des canyons de formes géométriques et de coefficient de impédance différents ont été obtenues. Ils se sont avérés avoir des effets significatifs sur la période naturelle. Les résultats ont été comparés aux données enregistrées réelles. Il a été constaté que les graphiques proposés dans cette étude peuvent être utilisés par des concepteurs de barrages pour l'estimation des périodes naturelles des barrages en terre dans des canyons de formes et de propriétés matérielles différentes. Plusieurs fonctions d'amplification correspondant à différentes conditions de canyon ont été obtenues en appliquant un déplacement uniforme à la limite du canyon. Une étude approfondie a été réalisée pour examiner les effets de la géométrie et de la flexibilité du canyon sur la réponse en régime permanent du barrage. Ces deux effets ont influencé de manière importante les fonctions d'amplification. Alors que la flexibilité du canyon affecte de manière significative la valeur de la fonction d'amplification maximale, cette valeur ne change pas pour les barrages en terre dans lesquels les canyons ont des formes différentes et la même longueur. De plus, la réponse latérale du barrage de terre dans le domaine temporel a été calculée pour analyser les effets susmentionnés lors d'un tremblement de terre réel. Les fonctions d'amplification proposées ont été utilisées pour comparer les spectres de réponse enregistrés du barrage d'El Infiernillo lors des tremblements de terre de 1966 avec la fonction d'amplification calculée. Un accord raisonnable a été observé entre eux. La méthode linéaire équivalente (EQL) a été implémentée dans le FEM. La technique FEMSBFEM a été étendue pour prendre en compte l'effet du comportement non linéaire des barrages en terre. Il a été observé que le comportement non linéaire affecte grandement la fréquence naturelle, la fonction d'amplification et l'accélération de crête maximale du barrage de terre situé dans les canyons. Les effets de la géométrie et de la flexibilité du canyon sur le comportement non linéaire ont été examinés, et on a vu qu'en augmentant la flexibilité du canyon, l'effet de la non-linéarité était diminué. Le barrage d'El Infiernillo a été modélisé par FEM-SBFEM non linéaire 3D, et une comparaison de la fonction d'amplification de crête obtenue par la méthode proposée avec les données enregistrées montre la précision du FEM-SBFEM non linéaire. / The canyon surrounding a dam can be assumed as an unbounded domain, and the geometry and flexibility of a canyon are parameters that greatly affect the values of natural periods in earth dams. In this thesis, in order to take into account these two effects, canyons are modeled by SBFEM, and earth dams, which have limited geometries, are modeled by FEM. The hybrid FEM-SBFEM technique used for the dynamic three-dimensional analysis of soil-earth dam interactions is validated with results available in the literature. Because the dynamic-stiffness matrix of the unbounded domain is complex and frequency-dependent, the classical mode-superposition method is not straightforward for a soil-structure interaction system, and thus, to obtain their fundamental natural frequencies, the modeled dams were excited in the upstream-downstream direction. The natural periods of earth dams in canyons with different geometries shapes and impedance ratios are obtained, and are found to have significant effects on the dams' natural periods. The results are compared with actual recorded data, and it is found that the graphs put forward in this study may be used by practical engineers for the estimation of natural periods of earth dams in canyons with different shapes and material properties. Several amplification functions corresponding to different canyon conditions are obtained by applying a uniform displacement at the canyons' boundaries. A comprehensive study is performed to examine the effects of canyon geometry and flexibility on the steady-state responses of the dams, and it is found that these two effects significantly influence the amplification functions. While the flexibility of the canyon does affect the maximum amplification function value, this value does not change for earth dams in canyons that have different shapes but the same length. In addition, the lateral responses of earth dams in the time domain are computed in order to analyze the aforementioned effects under an actual earthquake. The proposed amplification functions are used to compare the recorded response spectra of the El Infiernillo dam under the two 1966 earthquakes with the calculated amplification function, and a reasonable agreement is observed between them. The equivalent linear method (EQL) is implemented into the FEM, and the FEM-SBFEM technique is extended in order to take into consideration the effect of earth dams' nonlinear behavior. It is observed that such nonlinear behavior greatly affects the natural frequency, the amplification function, and peak crest acceleration of earth dams located in canyons. The effects of canyon geometry and flexibility on the nonlinear behavior are examined, and it is found that by increasing canyon flexibility, the effect of nonlinearity is decreased. The El Infiernillo dam is modeled by the 3D nonlinear FEM-SBFEM, and comparison of the crest amplification function obtained by the proposed method with the recorded data shows the accuracy of the nonlinear FEM-SBFEM.