Le système Géocube est un réseau de capteurs GPS conçu et développé par le Laboratoire d'OptoÉléctronique de Métrologie et d'Instrumentation (LOEMI) de l'Institut National de l'Information Géographique et Forestière (IGN) et maintenu par le même laboratoire et l'entreprise Ophelia- Sensors qui s'occupe de son industrialisation. Il a comme objectif de mesurer les déformations du sol avec une précision millimétrique. Ce réseau de capteurs a la particularité d'être à la fois très peu énergivore, d’un faible coût de revient, simple d’installation et d’utilisation. Il est donc bien adapté à l’usage dans un environnement difficile, comme les volcans. Ce système a déjà été testé avec succès lors d’une précédente thèse sur le glacier d’Argentière et sur un glissement de terrain proche de Super-Sauze en France. La première partie de cette thèse porte sur l’optimisation du système de calcul du Géocube pour l'adapter à des réseaux de tailles plus importantes horizontalement et verticalement en vue de son utilisation dans un contexte volcanique. Cela passe, d’abord, par l’intégration d’une stratégie pour l’estimation du biais troposphérique dans le filtre de Kalman qui constitue le coeur du logiciel de calcul du Géocube. Cette amélioration est ensuite validée en utilisant les données de quelques réseaux GNSS permanents nationaux et internationaux. La deuxième partie consiste à étudier l’apport d’un réseau dense de Géocubes à l’étude du volcanisme à travers une expérience conduite sur le flanc sud-est de l’Etna, où cinq Géocubes ont été déployés entre le 12 Juillet 2016 et le 10 Juillet 2017. Les résultats obtenus et les enseignements tirés de cette expérimentation sont discutés et analysés. Enfin, nous validons les résultats obtenus avec les Géocubes en appliquant une technique PSI (Persistent Scatterer InSAR) sur des interférogrammes RADAR calculés à partir des données des satellites Sentinel-1A/B et qui couvrent la période de déploiement des Géocubes sur l’Etna. Ces deux méthodes (GPS et RADAR) se sont avérées complémentaires puisque le RADAR apporte la densité spatiale des mesures et le système Géocube la précision et la continuité temporelle. / The Geocube system is a network of wireless GPS sensors designed and developed by the Laboratory of Opto-Electronics, Metrology and Instrumentation (LOEMI) of the National Institute of Geographical and Forest Information (IGN) and maintained by the same laboratory and Ophelia-Sensors, the company responsible for its industrialization. Its purpose is to measure ground deformations with millimetre accuracy. This sensor network has the particularity of being very low in energy consumption, low cost, easy to install and easy to use. It is suited for use in harsh environments, such as volcanoes. This system has already been successfully tested in a previous works on the Argentière glacier and a Super-Sauze landslide in France. The first part of this thesis deals with the optimization of the Geocube system for larger networks, horizontally and vertically, in order to use it in a volcanic context. First, a new strategy to estimate the tropospheric bias has been implemented into the Kalman filter (the heart of the Geocube processing software) in real time and in post-processing. This improvement is then validated using data from some national and international permanent GNSS networks. The second part consists in studying the contribution of a dense Geocubes network to the study of volcanism through an experiment conducted on the southeastern flank of Etna, where five Geocubes were deployed between July, 12th 2016 and July, 10th 2017. The results obtained from this experiment are discussed and analysed. Finally, the results obtained with Geocubes are validated by applying a PSI (Persistent Scatterer InSAR) technique on RADAR interferograms calculated from Sentinel-1A/B satellite data covering the period of deployment of the Geocubes on Etna. These two methods (GPS and RADAR) turned out to be complementary since RADAR provides the spatial density of measurements and the Geocube system provides accuracy and temporal continuity.

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.

Le Port de Montréal, la Garde Côtière Canadienne (GCC) et le Service Hydrographique du Canada (SHC) ont mis en place un projet commun pour sécuriser le passage des grands navires sous le pont de Québec. Entre autre, les mesures GNSS, archivées à la seconde de juillet 2012 à juillet 2013, sont utilisées pour étudier les déformations du pont de Québec. Les résultats ont permis de quantifier les déplacements 3D (et leur précision) de la travée centrale du pont en fonction des facteurs comme les conditions météorologiques et les passages des trains et des automobiles. Par exemple, dans les cas extrêmes, nous avons détecté que les trains de marchandise affaissent la travée suspendue du pont de 17 cm et qu’un vent de 100 km/h a poussé transversalement le centre de la travée suspendue de 17 cm. Entre l’hiver et l’été, les déformations de l’acier causées par les changements de température sont compatibles avec les prédictions. Ces mesures de déformation GNSS sont conformes, en général, aux prédictions faites lors de la construction du pont en 1907. / The Port of Montreal, the Canadian Coast Guard (CCG) and the Canadian Hydrographic Service (CHS) have established a joint research project to secure the passage of large ships under the Quebec Bridge. The GNSS measurements (along with radar measurements), archived every second from July 2012 to July 2013, are used to study the deformation of the Quebec Bridge. The results were used to quantify the magnitudes of the 3D movement (and its accuracy) of the central span of the bridge as function of external factors such as meteorological conditions and train and car crossings. In extreme cases, freight trains make the central span of the bridge to subside by 17 cm and a wind of 100 km/h pushed transversally the suspended span by 17 cm. Between summer and winter seasons, due to temperature variations, the steel deformations are compatible with the predictions. These GNSS deformation measurements are generally in good agreements with the predictions made during the construction of the bridge in 1907.

L'objectif du projet de recherche est de connaître le potentiel du piézocomposite MFC™ en tant qu'élément actif pour compenser les déformations induites thermiquement dans les structures composites. Un modèle numérique 3Dest donc développé à l'aide du logiciel de modélisation par la méthode des éléments finis ABAQUS afin de prédire la forme résultante des structures activées par l'actuateur MFC™. Afin d'utiliser le modèle, les propriétés élastiques et thermiques de l'actuateur sont estimées dans les plans 1-3 et 2-3 puisque les propriétés élastiques et thermiques sont données par le manufacturier seulement pour le plan 1-2. Par la suite, les coefficients d'expansion thermique et les coefficients de déformation piézoélectrique de l'actuateur sont mesurés. Une fois les propriétés de l'actuateur connues, sa capacité à déformer une structure simple est étudiée. La déflexion d'une plaque d'aluminium et d'une plaque composite [O4] induite par l'actuateur MFC™ lorsqu'une différence de potentiel est appliquée à ses bornes sont prédites par le modèle numérique. L'essai est ensuite reproduit expérimentalement. Finalement, la possibilité d'utiliser l'actuateur MFC™ afin de compenser les déformations thermiques dans les structures composites est étudiée. Pour ce faire, une plaque composite [0/90/90/90]T sur laquelle deux actuateurs MFC™ sont collés est modélisée dans le logiciel ABAQUS. Cette plaque est par la suite fabriquée et les résultats obtenus pour le modèle numérique et l'essai expérimental sont comparés. L'actuateur MFC™ peut compenser la déformation d'une plaque [0/90/90/90]T induite par une changement de température de 10°C. Le modèle éléments finis développé prédit précisément le comportement d'une structure simple activée par l'actuateur, mais est moins précis dans la prédiction de la forme résultante d'un composite [0/90/90/90]T soumis à un chargement thermique et activé par l'actuateur MFC™.

Simuler numériquement de façon précise les matériaux hyperélastiques en grandes déformations par la méthode des éléments finis est encore un problème difficile. Même avec l’aide d’un maillage très raffiné, la formulation lagrangienne totale peut mener à des problèmes de convergence en raison de la dégénérescence des éléments du maillage. L’estimation d’erreur et le remaillage adaptatif sur la géométrie initiale sont des outils utiles qui peuvent améliorer la précision des solutions (avec moins de degrés de liberté), mais ces outils ne sont malheureusement pas suffisants pour atteindre de très hauts niveaux de déformation. La formulation lagrangienne actualisée où la géométrie du domaine de calcul est périodiquement mise à jour est donc préférée, et ceci même si des étapes de remaillage sont encore nécessaires afin de contrôler la qualité des éléments et d’éviter les éléments trop déformés, voire renversés. Suite à une étape de remaillage, le transfert de données (réinterpolation des variables) de l’ancien maillage vers le nouveau maillage est nécessaire, ce qui est un problème délicat. Si les transferts ne sont pas effectués adéquatement, la précision peut être sérieusement affectée. Dans cette thèse, nous présentons une formulation lagrangienne actualisée où l’erreur sur la solution éléments finis est estimée et combinée au remaillage adaptatif afin de raffiner le maillage dans les régions où l’erreur estimée est grande, et au contraire, enlever des éléments là où l’erreur est considérée petite, le tout en contrôlant la qualité des éléments du maillage. En utilisant cette approche, de très hauts niveaux de déformation peuvent être atteints tout en préservant la précision de la solution. Une attention particulière est portée aux méthodes de transfert de données et une méthode de projection cubique très précise est introduite. La méthode de continuation de Moore-Penrose, qui est très efficace, est aussi utilisée pour piloter automatiquement l’algorithme complet qui inclut l’augmentation de la charge, l’estimation d’erreur, le remaillage adaptatif et le transfert de données. Une nouvelle approche pour l’implémentation de la méthode de continuation de Moore-Penrose, facilitant la détection des points de bifurcation, sera aussi présentée de même que plusieurs exemples. / Accurate simulations of large deformation hyperelastic materials by the finite element method is still a challenging problem. In a total Lagrangian formulation, even when using a very fine initial mesh, the simulation can break down due to severe mesh distortion. Error estimation and adaptive remeshing on the initial geometry are helpful and can provide more accurate solutions (with a smaller number of degrees of freedom) but are not sufficient to attain very large deformations. The updated Lagrangian formulation where the geometry is periodically updated is then preferred. Remeshing may still be necessary to control the quality of the elements and to avoid too severe mesh distortion. It then requires frequent data transfer (reinterpolation) from the old mesh to the new one and this is a very delicate issue. If these transfers are not done appropriately, accuracy can be severely affected. In this thesis, we present an updated Lagrangian formulation where the error on the finite element solution is estimated and adaptive remeshing is performed in order to concentrate the elements of the mesh where the error is large, to coarsen the mesh where the error is small and at the same time to control mesh distortion. In this way, we can reach high level of deformations while preserving the accuracy of the solution. Special attention is given to data transfer methods and a very accurate cubic Lagrange projection method is introduced. As large deformation problems frequently have highly nonlinear solutions, the Moore-Penrose continuation method is used to automatically pilot the complete algorithm including load increase, error estimation, adaptive remeshing and data transfer. A new approach for the implementation of the Moore-Penrose continuation method, facilitating the detection of bifurcation points, will also be presented as well as a number of examples.

Cette thèse contribue au domaine de la modélisation en centrifugeuse géotechnique des barrages en enrochement. Étant donné que les essais œdométriques sont généralement utilisés pour estimer les déformations verticales de la partie centrale des barrages en enrochement, une étude comparative a été réalisée entre les déformations verticales des modèles en centrifugeuses d’assemblages de particules de roche avec des pentes non supportées (2D) et les déformations verticales observées au cours des essais œdométriques (1D). Les deux types d’essais ont été effectués sur des assemblages de particules de roche avec la même granulométrie, tout en tenant compte de l’effet de la compacité. La comparaison a montré que pour une densité relative donnée, les modèles en centrifugeuses ont subi des déformations verticales plus élevées que celles des spécimens œdométriques et que le rapport entre les déformations verticales des sommets des modèles et celles obtenues lors des essais œdométriques diminue quand la densité relative augmente. De plus, une autre série d’essais en centrifugation a été réalisée sur des modèles de barrages en enrochement avec des noyaux centraux minces pour simuler le comportement typique de ce type de barrages pendant la construction et le remplissage du réservoir, tout en explorant les effets de la granulométrie et d’échelle sur la réponse de ces modèles. Ces modèles ont différentes granulométries alors qu’ils ont la même densité relative. Tout au long des essais, la rigidité de réponse a été plus élevée quand le coefficient d’uniformité des granulométries des modèles augmente pour la même taille maximale des particules et quand la taille maximale des particules diminue pour des granulométries parallèles. La taille des particules dans les grandes structures en enrochement telles que les barrages en enrochement empêche les essais de laboratoire sur ces matériaux et par conséquent, la prédiction du comportement sur le terrain nécessite le développement d’une méthode d’analyse afin d’extrapoler les résultats des essais de laboratoire tout en prenant les effets d’échelle et de la granulométrie en compte. Une méthode d’analyse est adoptée dans ce travail afin d’extrapoler les résultats des essais en centrifugeuses menés sur les modèles, déjà mentionnés, des barrages en enrochement avec des noyaux centraux minces, ce qui a permis la prédiction des déformations verticales de différentes zones du barrage de la Romaine 2. Cette méthode d’analyse est basée sur les effets d’échelle sur l’aptitude à la fragmentation des matériaux granulaires et il tient compte des résultats expérimentaux passés relatifs à cette question. Les valeurs prédites des déformations verticales de différentes zones du barrage de la Romaine 2 étaient dans le même ordre des valeurs des mesures in situ pendant la construction et le premier remplissage. / This thesis contributes to the domain of geotechnical centrifuge modeling of rockfill dams. Given that odometer testing is usually used to estimate the vertical deformations of the central part of rockfill dams, a comparative study was performed between the vertical deformations of centrifuge models of rock particle assemblages with unsupported slopes (2D) and the vertical deformations observed during oedometer tests (1D). Both types of testing were conducted on rock particle assemblages with the same gradation while considering the effect of compactness. The comparison showed that for a given relative density, the centrifuge models experienced higher vertical deformations than those of the oedometer specimens and the ratio between the vertical deformation of the top of the models and that obtained during the oedometer tests decreases as relative density increases. Furthermore, another series of centrifuge tests was conducted on models of rockfill dams with thin central cores to simulate the typical behavior of such type of dams during construction and reservoir impoundment while exploring the effects of gradation and size scale on the response of these models. These models have different gradations while they have the same relative density. Throughout the tests, the response stiffness was higher as uniformity coefficient of the models’ gradations increases for the same maximum particle size and as maximum particle size decreases for parallel gradations. Particle size in large rockfill structures such as rockfill dams prevent laboratory testing of these materials and consequently, the prediction of field behaviour requires the development of an analysis method in order to extrapolate the laboratory testing results while taking the effects of size scale and gradation in consideration. An analysis method is adopted in this work in order to extrapolate the results of the centrifuge tests conducted on the already mentioned models of rockfill dams with central thin cores, which allowed the prediction of vertical deformations of different zones of the Romaine 2 dam. This analysis method is based on size scale effects on the crushability of granular materials and it considers past experimental findings relating to this issue. The predicted values of vertical deformations of different zones of the Romaine 2 dam were in the same order of the field measurement values during construction and reservoir impoundment.

Cette thèse contribue au domaine des déformations post-construction des barrages en enrochement. D’abord, en étudiant le comportement in situ en compression de l’enrochement, il a été révélé que la compressibilité à court terme et la compressibilité en fonction du temps de l’enrochement sont reliés. De plus, la résistance à la compression simple (RCS) de la roche-mère et l’indice des vides de l’enrochement sont deux paramètres clés qui déterminent la compressibilité in situ de l’enrochement. Ainsi, ces paramètres, en plus du comportement des matériaux de remplissage pendant la construction (quantifiées par le module de déformation), ont été utilisés pour présenter des lignes de tendance pour prédire les déformations post-construction de différents types de barrages en enrochement. En considérant une base de données de 19 barrages en enrochement avec masque amont en béton et 22 barrages en enrochement zoné avec noyau en terre et en évaluant des tassements pendant et après la construction des lignes de tendance établissant un rapport entre les déformations post-construction et le comportement observé lors de la construction ont été développées. Dans les lignes directrices présentées, les tassements peuvent être obtenus séparément pour chacune des trois phases du cycle de vie des barrages : avant, pendant et après la mise en eau. De plus, le point zéro du temps pour les déformations en fonction du temps a été redéfini pour chaque phase. Lors du calcul des modules de déformation, les géométries de remblai ont été prises en compte en considérant un facteur de réduction pour les contraintes verticales à l'intérieur des remblais. Cette étude fournit une bonne base pour caractériser et estimer le comportement de l’enrochement pendant et après la construction des barrages en enrochement. Les tendances présentées montrent qu’en plus de surmonter les difficultés reliées aux méthodes actuelles (y compris les essais au laboratoire, les modélisations numériques, et des méthodes empiriques), les approches proposées sont précises et très pratiques. Cette recherche permet également une meilleure compréhension des mécanismes de déformation de l’enrochement. Outre l'utilisation de ces méthodes comme outils prédictifs des tassements, les relations présentées peuvent être utilisées afin de détecter un comportement atypique d'un barrage en enrochement au cours de sa durée de vie. / This thesis contributes to the area of post-construction deformations of rockfill dams. First, by studying the in-situ compression behavior of rockfill, it was revealed that the rockfill short-term compressibility and time-dependent compressibility are connected. Additionally, unconfined compressive strength (UCS) of parent rock and the rockfill void ratio are two key parameters that determine the in-situ (short-term and time-dependent) rockfill compressibilities. Thus, these parameters, besides the fill materials behavior during construction (quantified by defining deformation modulus) were used to present guidelines for predicting post-construction deformations of different types of rockfill dams. In other words, by gathering a comprehensive data base of 19 concrete face rockfill dams and 22 earth core rockfill dams and evaluating settlements during and after construction, trend lines were developed which relate the post-construction settlements to the observed behavior during construction. In the presented guidelines, settlements are obtainable separately for each of the three dams’ life cycle phases: before, during and after impoundment. Additionally, the time datum for time-dependent deformations was redefined for each phase. When calculating the material moduli, the embankment geometries were taken into account by introducing a reduction factor for vertical stresses inside embankments. This study provides a good framework for characterizing and estimating the behavior of rockfill during and after construction of embankments. The trend lines presented based on behavior of numerous rockfill dams show that, in addition to addressing some limitations of previous methods (including laboratory testing, numerical modeling, and empirical methods), the approach proposed for predicting post-construction settlements is accurate and highly practical. This research also permits a better understanding of rockfill deformation mechanisms. Apart from using the presented method for predictive purposes, the presented graphs can be used to distinguish unexpected settlement behavior of a rockfill dam during its post-construction lifespan.

Hydro-Québec reçoit plusieurs demandes de la part de compagnies minières pour pouvoir utiliser la crête de leurs ouvrages en remblai comme route d'accès aux ressources naturelles. La problématique liée à ces demandes d'autorisation de passage est le fait que lors de la conception d'un barrage en remblai, la circulation des charges lourdes en crête n'est pas considérée. Il est donc pertinent de développer des outils décisionnels pour les gestionnaires des barrages de manière à ce qu'ils soient en mesure d'octroyer ou d'appliquer certaines restrictions sur les passages demandés. Ces critères décisionnels doivent être basés, entres autres, sur la déformation permanente subie par le noyau du barrage. Concrètement, des instruments de mesure de contraintes et de déformations ont été installés dans l'extrémité supérieure du noyau d'argile du barrage en remblai de St-Timothée. Les contraintes horizontales et verticales mesurées lors du passage de camions en crête permettent de reproduire le cheminement de contraintes en laboratoire lors des essais de déformations permanentes. L'échantillon d'argile recompacté en laboratoire est soumis à des contraintes déviatorique cycliques allant jusqu'à des valeurs très élevées par rapport à ce qui a été enregistré sur le terrain. La contrainte axiale cyclique à laquelle l'argile accumule une certaine déformation permanente suite à sa post-compaction est de 5 kPa et celle à laquelle l'instabilité du matériau est observée correspond à 55 kPa. La charge en crête de l'ouvrage qui engendre ces incréments de contraintes axiales au sein du noyau est déterminée avec les données enregistrées avec les instruments mis en place dans le barrage. Le nombre de passages nécessaire pour engendrer des déformations permanentes de 50 mm et 100 mm au sein du noyau peut être estimé en fonction de différentes charges circulant en crête du barrage. Ceci est effectué en considérant une variation de la profondeur dans le noyau sur laquelle la charge en crête contribue à sa déformation permanente.

Ce mémoire présente une étude sur le contrôle actif, à l'aide de piézocomposites, des déformations induites thermiquement dans les structures composites. L'actionneur utilisé est le MFC™, un piézocomposite récemment développé, innovant et performant. Le cas d'une plaque stratifiée non-symmétrique comportant deux actionneurs est traité. Un modèle par éléments finis est développé afin de prédire le comportement de cette structure et d'étudier les capacités du MFC™. La théorie est tout d'abord validée dans le cas où les actionneurs sont passifs, puis pour la compensation des déformations induites thermiquement dans la structure. Enfin, ce modèle numérique est modifié afin d'intégrer un processus de contrôle actif. Une démarche expérimentale est développée afin d'étudier physiquement le comportement de la structure. La bonne corrélation entre les résultats théoriques et expérimentaux permet d'une part de valider le modèle numérique et d'autre part de mettre en évidence les capacités de l'actionneur MFC™.

La détermination des déformations et des contraintes qui s’exercent au passage d’une charge est une étape centrale des méthodes mécanistiques-empiriques de dimensionnement des chaussées. Ces réponses mécaniques sont utilisées afin d’identifier les modes d’endommagement de la chaussée et d’estimer sa durée de vie. Cette étude se concentre sur la mesure expérimentale et l’analyse des déformations à plusieurs niveaux dans les couches bitumineuses. Quatre structures de chaussées ont été instrumentées à l’aide de capteurs à fibre optique installés à 5 mm de la base du revêtement, à environ 15 mm de la surface du revêtement et à 5 mm de part et d’autre de l’interface entre les couches bitumineuses. Les mesures de déformations ont été effectuées lors d’essais en vraie grandeur en faisant varier les conditions de températures et les paramètres de charges tels que le type de pneu, le type d’essieu, et l’intensité de la charge et la vitesse. Le champ des déformations a été calculé à ces mêmes profondeurs à l’aide du logiciel ViscoRoute© 2.0 afin de les comparer aux mesures expérimentales. L’analyse de la déformation à la base du revêtement a permis de valider les capteurs à fibre optique développés au cours du projet et d’évaluer l’influence des paramètres de la charge. Cette validation repose sur une analyse de la répétabilité et de la dispersion des mesures ainsi que sur une étude comparative effectuée avec des capteurs résistifs usuellement utilisés dans les chaussées. Au sommet de la couche de roulement, les mesures dans la direction longitudinale et transversale confirment que le passage de la charge génère de la contraction sous le pneu et de l’extension à l’extérieur de celui-ci. Dans la direction verticale, la charge impose une forte extension à l’extérieur du pneu tandis que sous celui-ci, le signe et l’amplitude des déformations dépendent de la température et des propriétés viscoélastiques de l’enrobé bitumineux. Dans ces trois directions, l’amplitude et la forme des signaux sont fonction de la structure de chaussée et des caractéristiques de la charge. Les capteurs placés au niveau de l’interface montrent que de l’extension et de la contraction se développent respectivement au bas de la couche de surface et au sommet de la couche de base au passage de la charge. Ces niveaux de déformations peuvent être modélisés en considérant l’interface comme une couche viscoélastique de faible épaisseur. Ce comportement de l’interface influence la distribution des contraintes et des déformations et les cas évalués montrent qu’il doit être pris en compte pour l’identification des modes de sollicitations critiques des couches bitumineuses. / Calculating the stresses and strains induced by traffic loading is one of the main steps of the mechanistic-empirical pavement design method. In this method, structural responses are used to identify the mode of failure and estimated pavement life. The present study focuses on experimental strain measurements and analysis at several depths in the asphalt concrete layers. Fiber optic sensors, retrofitted in the asphalt-bound layers of four pavement structures, were used to measure strains at 5 mm from the bottom, at around 15 mm below the surface and at 5 mm below and above the interface between two bituminous layers. Strain measurements were taken during full-scale pavement testing under several climate conditions and traffic load parameters such as tire type, axle type, load level and speed. A strain field was computed at theses specific depths using the software ViscoRoute© 2.0 in order to compare the results with the experimental strains. Analysis of the strain signal at the bottom of the asphalt concrete layer allowed to validate the fiber optic sensors developed during the project and to characterize the influence of load parameters. This validation is based first on the analysis of the repeatability and the dispersion of the measurements and then on a comparative study with resistive sensors typically used in pavements. Near the surface course, longitudinal and transversal strain measurements confirm that a moving load induced compressive strain under the tire and tensile strain outside the tire edge. High levels of vertical tensile strain are induced outside the tire edge whereas the sign and the magnitude of strains directly under the tire depend on the temperature and the viscoelastic properties of the asphalt layer. In the three directions, the magnitude and form of the signals depend on the pavement structure and parameters of the load. Strain sensors on both sides of the interface show that the moving load induced tensile strain and compressive strain close to the bottom of the surface course and near the top of the base course, respectively. Numerical models, which consider the viscoelasticity of the tack coat at the interface between the asphalt layers, are capable of predicting these strain responses. This viscous-type interface behavior has an impact on strain and stress distribution through the layer and simulations show that it must be considered to identify the failure criteria.