This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivativesA program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively.First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined.The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Ce travail de thèse est consacré à l'étude des milieux granulaires humides. On étudie prin- cipalement le cas de faibles teneurs en eau. L'influence du phénomène de capillarité sur la rétention d'eau et sur le comportement mécanique du milieu est analysée. A chaque état d'équilibre, la distribution d'eau dans le milieu granulaire est obtenue en appliquant une succion homogène dans un volume élémentaire représentatif. La méthode multi-échelles utilisée se nourrit d'expérimentations/modélisation à l'échelle des interactions entre les grains, et produit des simulations en éléments discrets du comportement ma- croscopique qui sont comparées aux résultats d'expérimentation. A l'échelle des interactions capillaires, une étude expérimentale du pont liquide pilotée en succion permet la validation d'un modèle basé sur l'approximation toroïdale du profil du pont liquide. Dans le cas d'un essai à succion constante, on note que la distance de rupture des ponts liquides n'obéit plus au critère de Lian valide dans le cas d'un essai de traction à volume constant. A l'échelle macroscopique, les courbes de rétention d'eau simulées sont proches des courbes de rétention d'eau expérimentales réalisées sur des milieux modèles composés de billes de verre. Enfin, le comportement mécanique macroscopique est simulé. On note une fragilisation du matériau en fonction de la succion dans le domaines des faibles teneurs en eau qui s'explique par une diminution de la densité des ponts liquides. La prise en compte d'une rugosité des grains permet une meilleure description de la transition entre un état humide et un état sec. L'analyse des contraintes montre la pertinence du tenseur des contraintes associé aux interactions de traction pour analyser l'évolution de la contrainte à la rupture (compression simple) en fonction de la succion.

Cette thèse à caractère biographique est une étude de la carrière et de l'œuvre scientifiques de Gaspard-Gustave de Coriolis (1792-1843), polytechnicien, ingénieur des ponts et chaussées, dont le nom est universellement connu (force de Coriolis), mais dont le parcours et l'œuvre multiforme sont peu connus. Le fait que Coriolis fasse l'intégralité de sa carrière comme enseignant, puis directeur des études à l'École polytechnique, et que par ailleurs il participe à l'évolution de l'École et du Corps des ponts et chaussées à partir de 1830, nous donne un éclairage intéressant sur ces institutions (notamment grâce à la correspondance privée de Coriolis de 1838 à sa mort). Concernant son œuvre scientifique, sont rappelés ses apports de mathématicien. Son œuvre se caractérise principalement, toutefois, par une approche avant tout mathématique et théorique de la mécanique appliquée aux machines - il établit les fondements d'une théorie du travail. Cette approche et les résultats importants qui peuvent être attribués à ce savant, comme la définition physique du travail, les forces d'entraînement ou les forces centrifuges composées (forces de Coriolis), témoignent du lien tissé par Coriolis entre la mécanique rationnelle des géomètres et la mécanique appliquée à l'industrie naissante des machines.

Les analyses structurale et mécanique des pentes rocheuses constituent des éléments clés pour l'évaluation de leur stabilité. L'utilisation complémentaire de la photogrammétrie et des modèles numériques qui couplent les réseaux discrets de discontinuités (DFN selon son sigle en anglais) avec la méthode des éléments discrets (DEM selon son sigle en anglais), présente une méthodologie qui peut être utilisée pour évaluer le comportement mécanique des configurations tridimensionnelles de terrain pour lesquelles l'existence de discontinuités non persistantes peut être supposée. La stabilité des masses rocheuses est généralement supposée être contrôlée par la résistance au cisaillement le long des plans de discontinuité. Si les discontinuités sont non persistantes, avec leur continuité interrompue par la présence de ponts rocheux (portions de roche intacte reliant la masse rocheuse au massif), leur résistance apparente augmente considérablement. Dans ce cas, la contribution des ponts rocheux localisés entre ces discontinuités doit être prise en compte dans l'analyse de stabilité. La déstabilisation progressive des massifs rocheux dans lesquels des discontinuités non persistantes sont présentes, peut être étudiée par des simulations numériques réalisées à l'aide de l'approche DEM. La roche intacte est représentée comme un assemblage de particules (ou éléments discrets) liées entre elles par des contacts dont les lois de comportement spécifiques peuvent être calibrées pour représenter correctement le comportement de la roche. L'intérêt de la méthode est qu'elle permet de simuler l'initiation de la rupture et sa propagation à l'intérieur de la matrice rocheuse du fait de la rupture des contacts cohésifs entre les particules. De plus, les discontinuités préexistantes peuvent être prises en compte explicitement dans le modèle en utilisant une loi de contact ad hoc qui assure un comportement mécanique représentatif des plans de discontinuité. Des analyses de stabilité ont été effectuées et ont mis en évidence le rôle des ponts rocheux dans la génération de nouvelles surfaces de rupture qui peuvent se développer à travers des mécanismes de rupture mixte en traction et en cisaillement. On peut considérer la formulation de Jennings comme l'une des premières méthodes d'analyse de la stabilité des pentes rocheuses qui évaluent la résistance au glissement comme une combinaison pondérée des résistances mécaniques des ponts rocheux et des plans de discontinuité. Sa validité a été discutée et systématiquement comparée aux résultats obtenus à partir de simulations numériques. Il a pu être montré que la formulation de Jennings perd sa validité dès que la rupture des ponts rocheux intervient majoritairement par des mécanismes de traction. Une formulation complémentaire a alors été proposée. En ce qui concerne l'étude de la stabilité des massifs rocheux sur site, il a été montré que l'association entre les données issues de la photogrammétrie en haute résolution et l'approche DFN-DEM peut être utilisée pour identifier des scénarios de rupture. L'analyse en retour de cas réels a montré que les surfaces de rupture peuvent être simulées comme le résultat de mécanismes combinant la fracturation des ponts rocheux et le glissement le long des discontinuités préexistantes. La rupture d'un dièdre qui a eu lieu dans une mine de charbon australienne, a été utilisée pour valider cette méthodologie. Des simulations numériques ont été réalisées pour déterminer les scénarios pour lesquels les surfaces de rupture simulées et celles repérées sur le terrain, peuvent être utilisés pour calibrer les paramètres de résistance du modèle numérique. Le travail présenté ici répond à un besoin plus général visant à améliorer la gestion des risques naturels et miniers liés aux masses rocheuses instables. La méthodologie proposée constitue une alternative robuste dédiée à renforcer la fiabilité des analyses de stabilité pour les pentes rocheuses fracturées à structure complexe. / Structural and mechanical analyses of rock mass are key components for rock slope stability assessment. The complementary use of photogrammetric techniques and numerical models coupling discrete fracture networks (DFN) with the discrete element method (DEM) provides a methodology that can be applied to assess the mechanical behaviour of realistic three-dimensional (3D) configurations for which fracture persistence cannot be assumed. The stability of the rock mass is generally assumed to be controlled by the shear strength along discontinuity planes present within the slope. If the discontinuities are non–persistent with their continuity being interrupted by the presence of intact rock bridges, their apparent strength increases considerably. In this case, the contribution of the rock bridges located in-between these discontinuities have to be accounted for in the stability analysis. The progressive failure of rock slope involving non–persistent discontinuities can be numerically investigated based upon simulations performed using a DEM approach. The intact material is represented as an assembly of bonded particles interacting through dedicated contact laws that can be calibrated to properly represent the behaviour of the rock material. The advantage of the method is that it enables to simulate fracture initiation and propagation inside the rock matrix as a result of inter-particle bond breakage. In addition, pre–existing discontinuities can be explicitly included in the model by using a modified contact logic that ensures an explicit and constitutive mechanical behaviour of the discontinuity planes. Stability analyses were carried out with emphasis on the contribution of rock bridges failure through a mixed shear-tensile failure process, leading to the generation of new failure surfaces. Jennings' formulation being considered to be one of the first rock slope stability analysis that evaluates the resistance to sliding as a weighted combination of both, intact rock bridges and discontinuity planes strengths, its validity was discussed and systematically compared to results obtained from numerical simulations. We demonstrate that the validity of Jennings' formulation is limited as soon as tensile failure becomes predominant and an alternative formulation is proposed to assess the resulting equivalent strength. Regarding field slope stability, we show that the combination of high resolution photogrammetric data and DFN-DEM modelling can be used to identify valid model scenarios of unstable wedges and blocks daylighting at the surface of both natural and engineered rock slopes. Back analysis of a real case study confirmed that failure surfaces can be simulated as a result of both fracture propagation across rock bridges and sliding along pre-existing discontinuities. An identified wedge failure that occurred in an Australian coal mine was used to validate the methodology. Numerical simulations were undertaken to determine in what scenarios the measured and predicted failure surfaces can be used to calibrate strength parameters in the model. The work presented here is part of a more global need to improve natural and mining hazards management related to unstable rock masses. We believe that the proposed methodology can strengthen the basis for a more comprehensive stability analysis of complex fractured rock slopes.

L’étude de la voirie antique du site archéologique de Saint-Romain-en-Gal (Rhône), menée depuis vingt-cinq ans, a permis de mettre en lumière le rôle essentiel de la rue du Portique dans le développement de ce quartier de Vienne antique, situé sur la rive droite du Rhône. Traversant la plaine d’ouest en est, cette rue qui s’interrompait brutalement sur la berge du Rhône en face du cœur ancien de la capitale des Allobroges, a soulevé la question du franchissement du fleuve. Le nombre des ponts attribués à l’époque romaine a connu des évolutions au fil du temps. Une recherche historiographique importante a donc été menée afin d’établir la connaissance effective des indices prouvant la présence d’un ou de plusieurs ouvrages d’art antiques à la hauteur de Vienna, duplex urbs. Mais pour comprendre leur logique d’implantation, il a fallu également suivre la genèse d’installation des points de franchissement plus récents qui ont été soumis aux mêmes contraintes topographiques et urbaines que les ouvrages antiques. Il ressort de cette étude que deux ponts peuvent être attribués à l’époque romaine. Le troisième, plus monumental, le pont de pierre, considéré dans la mémoire collective comme le Pont Romain de la ville, semble toutefois remonter au Haut Moyen Âge. / The study of the ancient road to the archaeological site of Saint-Romain-en-Gal (Rhône), conducted for twenty-five years, helped to highlight the essential role of the Portico street in the development of this area of Antique Vienna, on the right bank of the Rhône. Across the plain from west to east, the street would strop abruptly on the bank of the Rhône opposite the ancient heart of the capital of the Allobroges, raised the issue of crossing the river. The number of bridges attributed to the Roman era has seen changes over time. An important historiographical research has been conducted to determine actual knowledge clues proving the presence of one or more works of ancient art at the height of Vienna, duplex urbs. But to understand their logic implementation, it took also follow the genesis of installing newer crossing points that were subject to the same topographical and urban and ancient books constraints. It appears from this study that two bridges can be attributed to the Roman era. The third, more monumental, stone bridge, seen in the collective memory as the roman bridge in the city, however, seems to go back to the High Middle Age.

Cette thèse est dévolue à l'étude de certaines propriétés asymptotiques du processus autorégressif d'ordre p. Ce dernier qualifie communément une suite aléatoire $(Y_{n})$ définie sur $\dN$ ou $\dZ$ et entièrement décrite par une combinaison linéaire de ses $p$ valeurs passées, perturbée par un bruit blanc $(\veps_{n})$. Tout au long de ce mémoire, nous traitons deux problématiques majeures de l'étude de tels processus : l'\textit{autocorrélation résiduelle} et la \textit{stationnarité}. Nous proposons en guise d'introduction un survol nécessaire des propriétés usuelles du processus autorégressif. Les deux chapitres suivants sont consacrés aux conséquences inférentielles induites par la présence d'une autorégression significative dans la perturbation $(\veps_{n})$ pour $p=1$ tout d'abord, puis pour une valeur quelconque de $p$, dans un cadre de stabilité. Ces résultats nous permettent d'apposer un regard nouveau et plus rigoureux sur certaines procédures statistiques bien connues sous la dénomination de \textit{test de Durbin-Watson} et de \textit{H-test}. Dans ce contexte de bruit autocorrélé, nous complétons cette étude par un ensemble de principes de déviations modérées liées à nos estimateurs. Nous abordons ensuite un équivalent en temps continu du processus autorégressif. Ce dernier est décrit par une équation différentielle stochastique et sa solution est plus connue sous le nom de \textit{processus d'Ornstein-Uhlenbeck}. Lorsque le processus d'Ornstein-Uhlenbeck est lui-même engendré par une diffusion similaire, cela nous permet de traiter la problématique de l'autocorrélation résiduelle dans le processus à temps continu. Nous inférons dès lors quelques propriétés statistiques de tels modèles, gardant pour objectif le parallèle avec le cas discret étudié dans les chapitres précédents. Enfin, le dernier chapitre est entièrement dévolu à la problématique de la stationnarité. Nous nous plaçons dans le cadre très général où le processus autorégressif possède une tendance polynomiale d'ordre $r$ tout en étant engendré par une marche aléatoire intégrée d'ordre $d$. Les résultats de convergence que nous obtenons dans un contexte d'instabilité généralisent le \textit{test de Leybourne et McCabe} et certains aspects du \textit{test KPSS}. De nombreux graphes obtenus en simulations viennent conforter les résultats que nous établissons tout au long de notre étude.

Ces travaux de thèse consistent en la conception et la réalisation d’une chaîne de tractionintégrée pour véhicule hybride plug-in. L’étude s’oriente vers une solution de convertisseur mutualisé,dans l’objectif de partager la traction et les modes chargeurs de batteries, la structure en NPC à 3niveaux est retenue. Le chargeur monophasé se base une topologie de redresseur à MLI monophaséavec trois bras entrelacés, avec l’utilisation des enroulements du moteur pour le filtrage. En chargeurtriphasé nous adaptons la topologie pour réaliser un montage en double boost triphasé. Pour chaqueconfiguration, les passifs sont dimensionnés pour répondre aux contraintes en courant BF et HF. Lecontrôle adopté se base sur les correcteurs résonants. Enfin, un prototype de 5 kW a été réalisé pourvalider les différents modes de l’application.Dans une seconde partie, nous proposons une solution de chargeur isolé sans étage continu auprimaire à double ponts actifs (DAB). La topologie est modélisée au premier harmonique et unecommande assurant l’absorption sinusoïdale est étudiée. Une configuration isolée triphasée permetl’accès aux puissances plus élevées ainsi que la réduction des ondulations de courant BF en sortie. / This thesis consists on the design and realization of a plug-in hybrid vehicle integrated tractiondrive supply. The work turns to a solution of a mutualized converter, in the objective to imagine asolution which shared drive and battery chargers modes, the three-level NPC topology has beenretained. The single phase charger is based on an interleaved PWM rectifier, and motor windings areused as smoothing inductors. A double-boost PFC configuration is introduced to ensure the threephasecharger. Passives are sized in each configuration in order to take in account the whole currentconstraints (LF and HF). The PFC behavior is based on the resonant controllers. Then, a 5 kWprototype has been realized to validate the different application modes.In a second part, a single-stage isolated charger based on a Dual-Active-Bridge (DAB) isproposed. The topology is modeled to the fundamental and the PFC control law is studied. A threephaseconfiguration is simulated in order to achieve higher charging powers and to reduce batterycurrent low-frequency ripple.

Les ponts roulants sont des engins de levage situés en haut des bâtiments qu'ils équipent. Ils servent à manutentionner des charges très lourdes et parfois critiques. Pendant un séisme, un pont roulant est exposé à des chocs multiples. Ces impacts peuvent causer des dommages importants dans la structure pouvant conduire à une chute de la charge manutentionnée ou du pont roulant lui-même. Donc, la vérification de la tenue des ponts roulants au séisme est une question primordiale. Actuellement, cette vérification est basée sur des méthodes de calcul statiques. Ces méthodes font l'hypothèse d'un comportement purement linéaire des ponts roulants ce qui les rend très conservatives. Depuis quelques années les niveaux sismiques imposés par les autorités nationales augmentent chaque année, et les constructeurs de ponts roulants se trouvent dans l'incapacité de construire à partir des efforts sur-estimés fournis par les méthodes statiques. L'objectif de la thèse est l'étude de la réponse dynamique d'un pont roulant pendant un séisme en prenant en compte les non-linéarités géométriques et matériau. Afin de modéliser ces phénomènes, une analyse dynamique temporelle avec une approche multi-échelle en temps est adoptée. Pour prendre en compte l'aspect haute fréquence des chocs, un intégrateur temporel variationnel explicite, basé sur la méthode des multiplicateurs de Lagrange et dédié au contact/impact, est développé. Ensuite, un intégrateur hétérogène (différents schémas d'intégration) asynchrone (différents pas de temps), basé sur la méthode de couplage GC, est appliqué au problème du pont roulant. Cette stratégie multi-échelle en temps permet d'adapter le schéma d'intégration et le pas de temps au sous domaine considéré. Par conséquent, l'intégrateur explicite est adopté dans les zones de contact et un schéma implicite de type accélération moyenne, est adopté dans le reste de la structure. Finalement, un démonstrateur de co-simulation entre les logiciels Cast3M et Europlexus est mis en place pour montrer le gain très significatif en temps de calcul dans le cas d'un modèle élément finis tridimensionnel d'un pont roulant industriel. / Bridge cranes are hoisting appliances located overhead in buildings. They are used to handle very heavy and sometimes critical loads. During an earthquake, a bridge crane may be subjected to multiple impacts between wheels and rails. These impacts can cause significant damage to the structure leading to a fall of the handled load or the bridge crane itself. Therefore, the qualification of such equipment, subjected to an earthquake, is very important. Currently, this qualification is based on static methods. These methods assume a purely linear behavior of the bridge cranes, which leads to a very conservative forces. Consequently, the bridge cranes manufacturers are sometimes unable to design the equipement from the over-estimated efforts provided by the static methods. The aim of this work is to study the dynamic response of a bridge crane during an earthquake by taking into account the geometric and material non-linearities. In order to model such phenomena, a time-history dynamic analysis with a multi-scale approach is performed. To take into account the high frequency aspect of the impacts between wheels and rails, a Lagrange explicit contact/impact time integrator is proposed. This work has also led to the development of an explicit–implicit HATI (Heterogeneous Asynchronous Time Integrator) for contact/impact dynamics. This method allows us to adopt an explicit contact/impact time integrator in the contact area and an implicit time integrator with a coarse mesh in the rest of the domain. Finally, a co-simulation demonstrator between Cast3M and Europlexus softwares is set up to show the very significant gain in computation time for a three-dimensional finite element model of an industrial bridge crane.

Abstract: Highway bridges have great values in a country because in case of any natural disaster they may serve as lines to save people’s lives. Being vulnerable under significant seismic loads, different methods can be considered to design resistant highway bridges and rehabilitate the existing ones. In this study, base isolation has been considered as one efficient method in this regards which in some cases reduces significantly the seismic load effects on the structure. By reducing the ductility demand on the structure without a notable increase of strength, the structure is designed to remain elastic under seismic loads. The problem associated with the isolated bridges, especially with elastomeric bearings, can be their excessive displacements under service and seismic loads. This can defy the purpose of using elastomeric bearings for small to medium span typical bridges where expansion joints and clearances may result in significant increase of initial and maintenance cost. Thus, supplementing the structure with dampers with some stiffness can serve as a solution which in turn, however, may increase the structure base shear. The main objective of this thesis is to provide a simplified method for the evaluation of optimal parameters for dampers in isolated bridges. Firstly, performing a parametric study, some directions are given for the use of simple isolation devices such as elastomeric bearings to rehabilitate existing bridges with high importance. Parameters like geometry of the bridge, code provisions and the type of soil on which the structure is constructed have been introduced to a typical two span bridge. It is concluded that the stiffness of the substructure, soil type and special provisions in the code can determine the employment of base isolation for retrofitting of bridges. Secondly, based on the elastic response coefficient of isolated bridges, a simplified design method of dampers for seismically isolated regular highway bridges has been presented in this study. By setting objectives for reduction of displacement and base shear variation, the required stiffness and damping of a hysteretic damper can be determined. By modelling a typical two span bridge, numerical analyses have followed to verify the effectiveness of the method. The method has been used to identify equivalent linear parameters and subsequently, nonlinear parameters of hysteretic damper for various designated scenarios of displacement and base shear requirements. Comparison of the results of the nonlinear numerical model without damper and with damper has shown that the method is sufficiently accurate. Finally, an innovative and simple hysteretic steel damper was designed. Five specimens were fabricated from two steel grades and were tested accompanying a real scale elastomeric isolator in the structural laboratory of the Université de Sherbrooke. The test procedure was to characterize the specimens by cyclic displacement controlled tests and subsequently to test them by real-time dynamic substructuring (RTDS) method. The test results were then used to establish a numerical model of the system which went through nonlinear time history analyses under several earthquakes. The outcome of the experimental and numerical showed an acceptable conformity with the simplified method. / Résumé: Les ponts routiers ont une grande valeur dans un pays parce qu’en cas de catastrophe naturelle, ils peuvent servir comme des lignes pour sauver des vies. Étant vulnérable sous des charges sismiques importantes, on peut considérer différentes méthodes pour concevoir des ponts routiers résistants et également pour réhabiliter des ponts existants. Dans cette étude, l'isolation de la base a été considérée comme une méthode efficace qui peut réduire significativement les effets des charges sismiques sur la structure. En réduisant la demande en ductilité sur la structure sans une augmentation notable de force, la structure est conçue pour rester élastique sous des charges sismiques. Le problème associé aux ponts isolés, particulièrement avec des appuis en élastomère, peut être leurs déplacements excessifs sous les charges de service et de séisme. Ceci peut défier l’objectif d'utiliser des appuis en élastomère pour les ponts typiques de petite portée où les joints de dilatation et les dégagements peuvent aboutir à une augmentation significative des frais d'exploitation et de maintenance. Ainsi, supplémenter la structure avec des amortisseurs d’une certaine rigidité peut servir de solution, ce qui peut cependant augmenter l’effort tranchant transmis à la sous-structure. Cette étude a pour but de fournir une méthode simplifiée afin d’évaluer les paramètres optimaux des amortisseurs dans les ponts isolés. Dans cette thèse, premièrement, basé sur une étude paramétrique, quelques directions sont données pour l'utilisation de dispositifs d'isolation simples, dont les appuis en élastomère, afin de réhabiliter des ponts existant avec une haute importance. Les paramètres comme la géométrie du pont, les clauses des normes et le type de sol sur lequel la structure est construite ont été appliqués sur un pont typique de deux portées. Il est conclu que les paramètres mentionnés peuvent déterminer l'emploi d'isolement de la base des ponts routiers. À la deuxième phase, basé sur le coefficient de réponse élastique des ponts isolés, une méthode de conception simplifiée d’amortisseur pour des ponts routiers réguliers isolés à la base a été présentée dans cette étude. En sélectionnant des objectifs pour la réduction du déplacement et la variation de l’effort tranchant, la rigidité et l'amortissement exigés d'un amortisseur hystérétique peuvent être déterminés. L’étude s’est poursuivie par une modélisation numérique d’un pont à deux portées pour vérifier l'efficacité de la méthode. Pour un modèle numérique d'un pont isolé typique, la méthode a été utilisée pour identifier des paramètres linéaires équivalents pour un certain déplacement et effort tranchant désigné. Par la suite, assumant un amortisseur de type hystérétique, les paramètres non linéaires de l’amortisseur ont été calculés et utilisés. La comparaison des résultats du modèle numérique sans amortisseur et avec l'amortisseur a démontré que la méthode proposée est suffisamment précise. Par la suite, un nouvel amortisseur hystérétique simple en acier a été conçu. Cinq spécimens ont été fabriqués de deux différents grades d’acier et ont été testés en combinaison avec un isolateur à l’échelle réelle dans le laboratoire de structures de l'Université de Sherbrooke. La procédure comprenait la caractérisation des spécimens par des tests cycliques en contrôle de déplacement et par la suite la réalisation d’essais par la méthode de sous-structuration dynamique en temps réel. Les résultats des essais ont été utilisés pour établir un modèle numérique du système qui a subi des analyses temporelles non linéaires sous plusieurs séismes. Le résultat des essais expérimentaux et numériques montrent une conformité acceptable avec la méthode simplifiée.

Le phénomène d'instabilité de la structure de fluide se produit fréquemment dans le système d'écoulement polyphasé des régimes sous-critiques de Reynolds, et le phénomène de vibrations de câbles de sustentation de ponts provoquées par la pluie (sigle RWIV en anglais) est considéré comme un exemple en génie civil pour caractériser l'instabilité causée par les interactions fluide-structure. Le RWIV est supposée être un nouveau type de vibrations; régulièrement accompagné de deux phénomènes significatifs: le ruisseau supérieur vibre circonférentiellement et la fréquence du vortex de Von Karman se déplace vers une valeur beaucoup plus faible par rapport à la fréquence attendue. Les phénomène est observé habituellement sur les haubans du pont à hauban immobilisé. En raison des mécanismes de couplage complexes, le mécanisme de RWIV n'a pas été complètement décrit par les chercheurs précédents. La plupart ont mis l'accent sur les observations sur le terrain et les aspects expérimentaux en soufflerie, mais rarement sur la simulation numérique. Pour élaborer un cadre numérique systématique, nous abordons le modèle non couplé, le modèle faiblement couplé, le modèle couplé, et le modèle multiphasique multi-échelle (MMM). L'objectif est de mettre en place un modèle numérique avec une grande exactitude et précision pour RWIV, et de reconnaître et clarifier le mécanisme de RWIV, diverses enquêtes numériques ont été réalisés dans cette thèse.Pour simuler les effets de la pluie/l’eau pluviale (eng. rainwater) comme un ruisseau artificiel (fixé / solide mobile attaché / oscillé le long de la circonférence de hauban immobilisé) lorsque RWIV se produit, la méthode séparée est mise en œuvre sur la base des équations incompressibles de Navier-Stokes en combinaison avec la simulation monotone intégré des grandes échelles (MILES) pour évaluer les sous-grille termes de pression. Les effets des ruisseaux artificiels de différentes positions le long de la circonférence de hauban immobilisé sur la structure de formation de tourbillons derrière le hauban, la distribution de pression à travers le hauban, et la force aérodynamique fréquence dominante du hauban sont analysés. Cependant, les enquêtes indiquent les positions de ruisseau artificiel le long de la circonférence du hauban affectent très faiblement la fréquence du vortex de Von Karman proche le sillage du hauban.Pour capturer l'évolution dynamique de la morphologie de l’eau pluviale , le modèle semi-couplé simplifie les équations incompressibles de Navier-Stokes avec la théorie de lubrification sur l'hypothèse qu'un mince film d'eau environnante autour du hauban. Les enquêtes indiquent que l’eau pluviale rassemble aux endroits près des points de séparation, et forme deux ruisseaux symétriques le long de la circonférence du hauban. Cependant, la vibration circonférentielle du ruisseau supérieur et les phénomènes de décalage de fréquence accompagnant RWIV ne peuvent être résolus et expliqués clairement.Afin d'améliorer le modèle de semi-couplé à suivre l'évolution de la morphologie de l’eau pluviale, la méthode du Volume-de-Fluide (VOF) combinée avec les équations incompressibles de Navier-Stokes est utilisée dans le modèle couplé. L ‘évolution hautement non linéaire du ruisseau de la pluie le long de la circonférence du hauban immobilisé et les caractéristiques aérodynamiques du hauban de séjour peuvent être obtenus et analysés. Les résultats indiquent que le ruisseau de la pluie est formé près des points de séparation le long de la circonférence du hauban; la zone de pression négative le long de la circonférence du hauban est pré-requise pour la formation de ruisseau supérieur. / A fluid structure instability phenomenon frequently occurs in the subcritical Reynolds regimes multiphase flow system, and rain--wind-induced vibration (RWIV) is taken as an example in civil engineering to characterize the aeroelastic instability caused by fluid-structure interactions. RWIV is hypothesized to be a new type of vibration; regularly accompanied by two significant phenomena: the circumferentially vibrating upper rivulet and the Von Karman vortex shedding frequency shift to a much lower value compared with the convectional evaluation; and customarily observed from the stay cables of cable--stayed bridge. Due to the complicated interactions mechanisms in the liquid-gas-solid system, the mechanism of RWIV has not been thoroughly solved and recognized by the previous researchers. Most have focused on the research topic from the field observation, the analytical dynamic model, and the wind tunnel experiment aspects, but rarely on numerical investigation aspect. To develop a systematic numerical framework, including the separated model, the semi-coupled model, the coupled model, and the multiphase multi-scale model (MMM) distinguished by different ways to simulate the rain effects when RWIV occurs, to establish highly accurate and precise numerical model for RWIV, and to recognize and clarify the mechanism of RWIV, various numerical investigations have been made in this thesis.To simulate the rain effects as an artificial rivulet (fixed/moving solid attaching/oscillating along the circumference of stay cable) when RWIV occurs, the separated method is implemented based on the incompressible Navier-Stokes equations in combination with the monotone integrated large eddy simulation (MILES) to evaluate the sub-grid stress terms. The effects of various artificial rivulet positions along the circumference of stay cable on the vortex shedding structure behind the cable, pressure distribution around the cable, and the aerodynamic force of the cable are analyzed. However, investigations indicate the positions of artificial rivulet along the circumference of cable extremely weakly affect Von Karman vortex shedding frequency near the wake of the cable.To capture the dynamic rainwater morphology evolution, the semi-coupled model simplifies the incompressible Navier-Stokes equations with the lubrication theory on the assumption that a thin water film surrounding around the cable. The investigations indicate the rainwater gathers at the locations near the separation points, and forms two symmetrical rivulets along the circumference of cable. However, both the circumferentially vibrating upper rivulet and the frequency shift phenomena accompanying RWIV cannot be solved and explained detailedly and clearly. To improve the semi-coupled model on tracking the rainwater morphology evolution, volume-of-fluid (VOF) method combined with incompressible Navier-Stokes equations is employed in the coupled model. Both the high-nonlinear rainwater rivulets evolution along the circumference of cable and the aerodynamic characteristics of stay cable can be obtained and analyzed. The results indicate rainwater rivulet are formed near the separation points along the circumference of cable; the negative pressure zone along the circumference of cable provides a prerequisite for the formation of upper rivulet. However, the computational efficiency is reduced due to the smaller droplets scatter in the surrounding air, furthermore, the assumptions, surrounding the stay cable with the constant volume of rainwater, cannot reflect the real physical conditions (i.e., rain infall process) and cannot obtain the real aerodynamic force from physical aspect.