Dans le domaine de la sismologie, la modélisation numérique de la propagation des ondes est un outil essentiel pour, par exemple, la quantification du risque sismique ou l'imagerie du sous-sol par inversion de données sismiques. Ces applications géophysiques nécessitent le développement de méthodes numériques sophistiquées de modélisation des ondes dans des milieux de propagation 3D hétérogènes. Parmi les approches possibles,les méthodes directes fondées sur la représentation du milieu sur des grilles numériques permettent la prise en compte de toute la complexité du champ d'onde dans des milieux fortement hétérogènes. Dans ce contexte, la représentation optimale du milieu de propagation et du champ d'onde propagé sur des grilles numériques irrégulières est actuellement. un axe de recherche fortement développé car elle permet d'atteindre le meilleur compromis entre la précision et la rapidité d'exécution des simulations dans des grilles numériques de taille considérable. Dans cette thèse, je propose une nouvelle méthode de modélisation de la propagation des ondes acoustiques par différences finies dans le domaine des fréquences . Le problème numérique associé est la résolution d'un système matriciel. L'approche développée propose une discrétisation spatiale optimisée de l'équation d'onde rendant possible des modélisations en 3D. Premièrement, je compare des modélisations acoustiques en milieu hétérogène calculées avec des schémas aux différences finies en quinconce et. des schémas optimisés, discrétisés suivant deux systèmes d'axe tournés de 45 degrés (schémas tournés). Pour des modélisations 2D, la seconde approche montre des performances supérieures en terme de temps calcul et de stockage mémoire. J 'en conclus que la discrétisation suivant des grilles tournées doit être utilisée pour des simulations 2D. Dans le cas 3D, les schémas utilisant les deux types de discrétisation contiennent 19 coefficients dans le cas des grilles en quinconce et 27 coefficients dans le cas des schémas tournés. Je suggère que les performances, en termes de précision numérique et de stockage mémoire, des deux types de discrétisation doivent être examinées dans le futur. Deuxièmement, je propose une approche mixte multi-échelles combinant une méthode directe et itérative de résolution du système matriciel qui permet d'effectuer des modélisations réalistes dans des milieux 3D. Dans un premier temps, le champ d'onde exact est calculée via la factorisation LU de la matrice sur une grille de résolution grossière. Deuxièmement, ce champ d'onde approché est interpolé sur une grille numérique de résolution fine pour être injecté comme estimé initial dans l'algorithme de résolution itérative. J 'ai exploré deux approches multi-échelles fondées respectivement sur une" nested iteration ", appellée méthode "Direct-Iterative-Space Solver" (DISS), et une décomposition multigrilles en ondelettes, appellée méthode" Direct-lterative- Wavelet Solver" (DIWS) . Les performances respectives, en termes de temps-CPU et de stockage mémoire, de ces deux méthodes de modélisation des ondes en milieu 2D et 3D sont comparées. L'algorithme de la méthode DISS permet la construction rapide du système matriciel. La convergence de l'algorithme itératif dépend fortement de sa capacité à éliminer des artefacts numériques associés à des phénomènes déphasage. L'importance de ces artefacts dépend fortement de la précision de la solution initiale et de la discrétisation numérique du champ modélisé dans l'algorithme itératif. Un nombre important d'itérations ou l'utilisation de cycles en V et W sont nécessaires pour supprimer ces phénomènes. Dans la méthode DIWS, une représentation multirésolution du système matriciel est développée par projection sur une base d'ondelettes. Malgré le coût induit pour transformer le système dans le domaine spectral, la représentation multirésolution a fourni un outil numérique de préconditonnement permettant de stabiliser et d'accélérer significativement l'algorithme itératif. L'accélération de la convergence est attribuée à la représentation du système sur plusieurs grilles numériques de résolution différente dont les interactions sont prises en compte au cours des itérations. Par ailleurs, la formulation en ondelettes ouvre des perspectives d'optimisation sous forme d'adaptation spatial du maillage en fonction des propriétés locales du milieu et du champ propagé. Plusieurs exemples de simulation d'ondes dans des milieux 2D et 3D de complexité variable sont présentés pour illustrer les performances respectives des méthodes DISS et DIWS. La taille des applications 3D présentées a été limitée par l'utilisation de programmes séquentiels. Finalement, des stratégies de parallélisation en mémoire distribuéedes codes sont proposées.

Le travail réalisé pendant cette thèse a permis d'apporter de nouvelles contraintes sur le développement et la distribution de la déformation dans le manteau supérieur et plus particulièrement au niveau des grandes limites de plaques décrochantes. Grâce à l'apport de l'expérience USArray et d'une dizaine d'années d'enregistrements sismologiques supplémentaires, nous avons pu étudier, de manière précise, les variations d'anisotropie dans le voisinage de la Faille de San Andreas. Nous avons confirmé et étendu l'observation de deux couches anisotropes sous cette limite de plaque. On y observe une première couche localisée dans la lithosphère marquant la déformation induite à la limite de plaque, et une autre, asthénosphérique, cohérente avec l'anisotropie observée loin de la faille et d'origine plus discutée. Nous avons montré que la zone de déformation associée aux failles de San Andreas, Calaveras et d'Hayward a, vraisemblablement, une largeur d'au moins 40 kilomètres en base de lithosphère, sous chacune de ces failles. Nous avons ensuite procédé à la modélisation thermomécanique (ADELI) de la migration d'une limite de plaques décrochante couplée à une modélisation du développement de fabriques cristallographiques par une approche viscoplastique auto-cohérente (VPSC). Ceci nous a permis d'y observer le développement de la déformation et les conséquences des possibles interactions entre la déformation décrochante en surface et le cisaillement en base de lithosphère dû au déplacement horizontal des plaques. Les propriétés élastiques déduites des fabriques cristallographiques modélisées montrent que de telles interactions existent et provoquent, sous la limite de plaques, une rotation des orientations cristallographiques avec la profondeur. Le signal associé à ces rotations progressives n'est toutefois pas cohérent avec la présence de deux couches d'anisotropie comme proposée sous la faille de San Andreas. Nous pensons par conséquent qu'il existe, sous la Californie, une zone de découplage entre la lithosphère et l'asthénosphère, permettant d'individualiser une déformation lithosphérique d'une déformation asthénosphérique. Nous estimons, en outre, que l'anisotropie observée dans l'asthénosphère sous la Californie ne peut être expliquée seulement par le cisaillement induit par le déplacement de la lithosphère Nord Amérique. En effet, les propriétés anisotropes obtenues par modélisation à partir d'une plaque se déplaçant dans une direction et une vitesse proche de celle de la plaque Amérique du Nord montrent qu'on ne peut espérer guère plus que quelques dixièmes de seconde de délai au bout de 10 Ma de déplacement. Les déphasages mesurés en Californie étant de l'ordre de 1,5 s, il est donc nécessaire d'invoquer la présence d'écoulements mantelliques actifs sous cette région / This work provides new constraints on the development and on the distribution of the deformation in the upper mantle and particularly beneath transform plate boundaries. USArray experiment and the remarkable increase of the dataset in California for the past ten years allowed us to scrutinize the lateral variations of the anisotropy in the vicinity of the San Andreas Fault zone. We have confirmed and increased the detection of two layers of anisotropy beneath this plate boundary. The first layer, located in the lithosphere, is related to the deformation induced at the fault, and the other one, located in the asthenosphere, is coherent with the anisotropy observed far from it, its origin is however less clear. We show that the deformation zone associated both to the San Andreas, Calaveras and Hayward Faults, is likely 40 km wide at 70 km depth. We then performed numerical thermomechanical modeling (ADELI) of the displacement of a transform plate boundary associated with the computation of the development of crystallographic fabrics using a viscoplastic self-consistent approach (VPSC). We analyzed the distribution of the deformation in the model ant looked after the possible interactions at depth between deformation caused at surface by the strike-slip dynamic of the fault and the shearing at the base of the lithosphere caused by the horizontal displacement of the plates. Elastic properties derived from the crystallographic fabrics modeled, show that such interactions exist and induce, beneath the fault zone, a progressive rotation of the crystallographic fabrics with depth. Seismological signature of these smooth rotations is however not relevant with the presence of two anisotropic layers as proposed beneath California. We thus consider that a decoupling zone exists between the lithosphere and the asthenosphere beneath the California to account for the sharp separation between a lithospheric and an asthenospheric deformation. We furthermore estimate that anisotropy observed far form the San Andreas Fault in California cannot be explained only by the drag of the asthenosphere by the North America lithosphere as proposed in our article. Indeed, we can only expect few tenths of second of splitting delay from the anisotropic properties derived from the numerical modeling of a plate moving in the same direction and in the same velocity than the North American lithosphere only for 10 Ma of displacement. As delays observed in California rather reach 1.5 s, anisotropy in this region thus requires the existence of an active asthenospheric flow to be explained.

A period of cement hydration is one of critical in the life span of concrete structures. One of the reasons of collapse of concrete structural elements during and after construction is the error in determining the concrete characteristics at early age. Recently, non-destructive test emerged as a popular way to evaluate the properties of cement-based materials. This test offers continuous measurements of concrete properties as well as ability to monitor any changes in the current state of structural materials. In the existing research, some of these methods fail to capture well properties of the materials in the plastics stage. A new piezoelectric pulse testing device (Piezoelectric Ring Actuators Technique), (P-RAT ) was initially developed at the University of Sherbrook as a non-destructive test (NDT ) for soil. This technique is considered a completely new, versatile, advanced and accurate. The development of the new technique (P-RAT ) was done on two main bases: the first was the development of piezoelectric ring actuators set-up and the second is the development of the interpretation method. The setup is composed of two main units; emitter and receiver, and is capable of measuring shear and compression wave velocities in specimens. With this technique, many problems of pulse tests, which make interpretation of results difficult and ambiguous, were solved in soil. The P-RAT overcomes wave reflections at boundaries (end-caps and sides), sample disturbance, weak shear coupling between soil and device (interaction) as well as the fixation problems, low resonant frequency and limited input voltage of the existing device. The previous method is exploited forward to measure the hydration properties of cement-based material. To apply this test method, it is necessary to determine how the evolution of shear wave velocity can be related and sensitive to the hydration of cement-based materials. Validation of the P-RAT with four conventional test setups that can be used to monitor early setting and hydration of cement-based materials is carried out. These tests include penetration resistance to monitor initial and final setting respectively, calorimetric to monitor heat of hydration, electrical conductivity to monitor change in continuity of the pore structure and compressive strength at 24 hours. The phase one of this investigation included trial tests to investigate the possibility of employing the original setup used for soil (P-RAT ) to determine setting and hardening properties of cement-based material. Based on the results of the preliminary test, two modifications were conducted to the previous test device to fit with cement based material and to obtain adequate resonant frequency for cement-based materials. These modifications are the design of the container and changing the dimensions of the rings. The resultant version of P-RAT after the modification was referred to be as P-RAT2 . Calibration of the P-RAT2 with water specimen was undertaken using the compression wave velocity and resulted in 99.33% accuracy. One paste mixture was tested three times to determine the experimental error of the P-RAT2 . The repeatability carried out on the P-ART2 proved the ability of this setup to capture accurate results of the shear wave velocity. This relative error is limited to 9 %. A number of series of validation was performed on cement paste and mortar mixtures proportioned with various water cement ratios (w/cm ) as well as chemical admixtures. The w/cm ratio ranged between (0.35 and 0.50). The investigated chemical admixtures comprise of high-range water-reducing agent, viscosity-modifying agent, set-accelerating agent, and set-retarding agent. The presented validations examine the ability of a P-RAT2 to monitor the hydration of the cement-based materials. The hydration is characterized by setting time, heat of hydration, electrical conductivity, and compressive strength at 24 hours. The results obtained using the P-RAT2 was correlated to those obtained using the conventional tests and strength measurement. The results enable to validate the ability of P-RAT2 to accurately detect variations in the hydration of cement-based materials. In addition, the initial and final time of setting can be determined from the derivation of velocity vs. time curve. The results show that conductivity, resistivity, has a bilinear relationship to shear wave velocity. The compressive strength at 24 hours was correlated to both the shear wave velocity and shear modulus obtained using the P-RAT2 . Furthermore, analytical model was derived to estimate the w/cm in mortar mixture by measuring the shear wave velocity (V[subscript s] ) and the corresponding time (t )||Résumé : La période d'hydratation du ciment est l'une des périodes clé du cycle de vie des structures en béton. L'une des raisons de l'effondrement d'éléments structuraux en béton pendant et après la construction peut être attribuée à une détermination des caractéristiques au jeune âge erronée. Depuis quelques années, l’auscultation des structures est devenue une méthode très populaire pour évaluer les propriétés des matériaux cimentaires. Cette méthode permet d'obtenir les propriétés du béton en continue et possibilité un suivi de changements dans I'état des matériaux structuraux. Dans I'état actuel de la recherche dans ce domaine, certaines méthodes ne sont pas adéquates pour bien mesurer les propriétés des matériaux à I'état plastique. Un nouvel appareil d'essai à impulsions piezoélectriques (Piezoelectric Ring Actuators technique), (P-RAT) a initialement été développé à l’Université de Sherbrooke comme technique d'auscultation des sols. Cette technique est considérée complètement nouvelle, polyvalente, évoluée et précise. Le développement de cette nouvelle technique (P-RAT) a été effectué en deux volets : la première sole est le développement du dispositif de vérin de commande annulaire piezoélectrique et le deuxième est le développement d'une méthode d'interprétation. Le dispositif d'essai est composé de deux unités principales, un émetteur et un récepteur. Et permet de mesurer la vitesse de propagation des ondes de cisaillement et de compression. Grâce à cette technique, plusieurs des problèmes associés aux dispositifs d'essais par impulsion des ondes, qui rendent les résultats ambigus et difficiles à interpréter, ont été résolu pour les sols. Le dispositif P-RAT surmonte les problèmes de réflexion des ondes aux limites (embouts et côtés), la perturbation de l’échantillon, le couplage de cisaillement faible entre le sol et l'appareil (interaction) ainsi que les problèmes de fixation, la faible résonnance des fréquences et le voltage d'entrée limité du dispositif. La méthode décrite a été utilisée pour mesurer les propriétés d'hydratation des matériaux cimentaires. Pour pouvoir appliquer cette méthode, il faut déterminer comment l’évolution de la propagation des ondes de cisaillement peut être reliée à l'hydratation des matériaux cimentaires et être sensible à ces dernières. La validation de la méthode P-RAT est réalisée, à l'aide de quatre configurations conventionnelles que l’on peut utiliser pour faire le suivi de la prise et de l'hydratation des matériaux cimentaires. Ces essais consistent à la résistance à la pénétration afin de pouvoir déterminer la prise initiale et finale, la calorimétrie pour suivre l’evolution de la chaleur d'hydratation, la conductivité électrique pour effectuer le suivi de la structure des pores et la résistance à la compression à 24 heures. La phase 1 de l'étude comprend des essais pour évaluer la possibilité d'utiliser la configuration originale utilisée pour les sols (P-RAT) pour déterminer les propriétés de prise et de durcissement des matériaux cimentaires. Selon les résultats des essais préliminaires, deux modifications ont été effectuées à l'appareil original pour permettre son utilisation avec des matériaux cimentaires et pour obtenir une fréquence de résonnance raisonnable sur les matériaux cimentaires. Les modifications effectuées sont la conception du contenant et un changement de la dimension des anneaux. La version modifiée du P-RAT est designée P-RAT2. La calibration du P-RAT2 à l'aide d'échantillon liquide sous propagation d'ondes de compression a été menée, avec des résultats d'une précision de 99,33 %. Un mélange a été testé trois fois pour déterminer l'erreur expérimentale du P-RAT2. La répétitivité des essais sur le PART2 a démontré la capacité de cet appareil à produire des résultats de cisaillement de propagation des ondes de cisaillement très précis. L'erreur relative se limite à 9 %. Une série d'essais de validation a été menée sur des mélanges de pâte de ciment et de mortier de rapport eau/ciment variés (e/c) ainsi qu'avec des adjuvants. Le rapport e/c variait entre 0,35 et 0,50. Les adjuvants utilisés étaient des supers plastiflants (réducteur d'eau), des agents de viscosité, des agents accélérateurs de prise et des agents retardateurs de prise. Les validations présentées ont pour but de valider la capacité du P-RAT2 à suivre l'hydratation des matériaux cimentaires. L'hydratation est caractérisée par le temps de prise, la chaleur d'hydratation, la conductivité électrique et la résistance à la compression à 24 heures. Les résultats obtenus à l'aide du P-RAT2 ont été comparés à ceux obtenus à l'aide d'essais de mesure des caractéristiques physiques et de résistance traditionnels. Ces résultats permettent de valider la capacité du P-RAT2 à détecter avec précision les variations dans l'hydratation des matériaux cimentaires. De plus, le dispositif P-RAT2 peut avoir une correction avec mesure obteniez avec les appareils traditionnels. II est aussi possible de déterminer le temps de prise initial et final à l'aide d'une courbe de propagation vs le temps. Les résultats montrent que la conductivité et la résistivité ont une relation bilinéaire à la propagation des ondes de cisaillement. La résistance à la compression à 24 heures a été comparee à la fois à la propagation des ondes de cisaillement et au module de cisaillement obtenus avec le P-RAT2. De plus, un modèle analytique a été établi pour estimer le rapport e/c dans le mélange de mortier en mesurant la propagation des ondes de cisaillement (V) correspondant au temps (t).

Le but de cette thèse est d'implémenter une méthode déterministe pour modéliser le canal radio des réseaux wifi. Pour cela la méthode MR-FDPF (Multi Resolution Frequency Domain ParFlow) a été adaptée, Le paramétrage et la calibration du modèle sont présentés, donnant une erreur inférieure a 4dB. Nous montrons comment prendre en compte les autres étages d'un bâtiment en développant la méthode 2.5D. De plus pour bien prendre en compte les diagrammes d'antennes des points d'accès un nouveau procédé est présenté. Enfin l'implémentation de la méthode MR-FDPF en 3D est développée ainsi que les optimisations nécessaires Enfin trois extensions sont proposées: la première concerne une extension de la méthode en milieu urbain, la deuxième présente une méthode de prédiction des évanouissements, et la troisième explique comment le modèle a été intégré dans un logiciel de planification de réseau WiFi.

Les budgets destinés à l'entretien et à la réparation du patrimoine bâti étant très élevés, il est primordial de mettre en place des techniques fiables de suivi et de diagnostic des ouvrages pour réduire ces dépenses. Il existe de nombreuses techniques de contrôle non destructif (CND) employant différentes technologies à des coûts plus ou moins abordables. Elles permettent la détermination des propriétés du matériau renseignant sur l'état de l'ouvrage, et la détection d'inclusions ou de défauts dans le béton. Ces techniques sont connues pour leur rapidité de mise en œuvre et leur grand rendement sur site, mais aussi par leur caractère non intrusif qui permet de limiter le nombre de prélèvements. Cela les rend pertinentes pour la surveillance et l'inspection du patrimoine bâti, puisqu'elles peuvent fournir des informations liées à l'état d'un ouvrage sans l'endommager. Dans un ouvrage en béton armé, la couche d'enrobage des armatures est la plus sujette à de multiples agressions. C'est dans cette couche qu'ont lieu la majorité des réactions chimiques qui causent la dégradation de l'ouvrage, puisqu'avec la présence d'eau cet enrobage est considéré comme la zone de pénétration, de transfert et de fixation des agents agressifs présents dans l'environnement de l'ouvrage, favorisant par exemple la corrosion des armatures. La teneur en eau du béton est considérée comme un facteur très important dans le diagnostic des structures en béton armé, puisqu'elle conditionne le développement de la plupart des pathologies liées à la dégradation des structures affectant ainsi la durabilité du béton. La localisation avec précision des aciers de renforcement est également un élément primordial pour le diagnostic, mais également pour vérifier la conformité des dispositions constructives par rapport aux plans d'exécution. Dans ces travaux de recherche, nous utilisons la technologie radar pour la caractérisation géométrique et physique des éléments en béton armé, cela par le développement d'un algorithme automatique, rapide et sans prérequis. Les objectifs de cette étude sont : de localiser en surface et en profondeur les aciers du premier lit d'armatures, et de déterminer la vitesse de propagation des ondes électromagnétiques dans le béton d'enrobage pour estimer sa teneur en eau. / Since budgets for the maintenance and repair of built heritage structures are very high, it is essential to implement reliable techniques for monitoring and diagnosis of structures to reduce these costs. There are many non-destructive testing (NDT) methods employing different technologies which are more or less expensive. The NDT methods allow the determination of the properties of the material giving also information about the condition of the structure, and the detection of inclusions or defects in the concrete. These techniques are known for their fast implementation and their high performance in-situ, but also by their non-intrusive nature, which makes it possible to limit the number of samples for destructive assessment. This makes them relevant for the monitoring and inspection of built heritage, since they can provide information related to the condition of a structure without any damage. In a reinforced concrete structure, the concrete cover is the most prone to multiple aggressions. The majority of the chemical reactions which cause the degradation of the structure occur in this layer, since with the presence of water this coating is considered as the penetration, transfer and fixing zone of aggressive agents present in the environment of the structure, thus enhancing for example the corrosion of the reinforcements. The water content of concrete is considered as a very important factor in the diagnosis of reinforced concrete structures, since it conditions the development of most pathologies related to structural degradation, thus affecting the durability of concrete. The precise location of the reinforcement in the construction is also an essential element for the diagnosis, but also to check the conformity of the constructive dispositions with the execution plans. In this research, we use GPR technology for the geometrical and physical characterization of reinforced concrete elements, through the development of an automatic algorithm, fast and without pre-requisites. The objectives of this study are to locate the rebars of the first reinforcing layer in the plane and depth and to determine the propagation velocity of the electromagnetic waves in the coating layer to estimate the water content of concrete.

De nombreux travaux concernant les performances des réseaux radio ambiants reposent sur une modélisation assez simpliste de la couche physique, qu'il s'agisse de réseaux d'accès sans fil (RLAN), de réseaux ad hoc ou de réseaux de capteurs. Si une telle modélisation permet d'établir des bornes et des lois simples conduisant au développement d'approches analytiques, le passage direct à un environnement réel en montre les limites. Certains protocoles révèlent en effet des failles importantes et inattendues, alors que certains problèmes s'atténuent naturellement. <br />Le premier axe de notre travail est d'améliorer les outils de simulation et de modélisation de la couche physique afin de réduire l'écart entre ces modèles et l'environnement réel. <br />Le deuxième axe consiste à évaluer l'impact d'une telle modélisation sur les performances des réseaux. Pour cela, nous étudions deux problématiques. La première concerne la planification des réseaux sans fil et la deuxième concerne la connexité des réseaux de capteurs.

L'imagerie des paramètres physiques du sous-sol à partir d'enregistrements sismiques de surface constitue un problème inverse non linéaire. L'inversion des formes d'onde élastique est une méthode d'imagerie quantitative multiparamètres de diffractants, nécessitant au préalable la connaissance précise d'un macromodèle de vitesse. <br /><br />Le problème direct associé, la propagation des ondes élastiques, est résolu dans le domaine fréquentiel, permettant la prise en compte efficace d'acquisitions multisources et multirécepteurs, par une méthode numérique de différences finies modélisant la propagation de tous les types d'onde (ondes de volume, de surface, diffractées, réfractées ...). Le stencil de différences finies choisi permet simule précisément la surface libre et la propagation des ondes de surface. <br /><br />Cette inversion linéarisée s'appuie sur une méthode de gradient, qui minimise une fonction coût contenant les différences entre données observées et calculées. Deux paramètres sismiques sont reconstruits à partir de sismogrammes verticaux et horizontaux. L'inversion est effectuée des basses fréquences vers les hautes fréquences, introduisant des longueurs d'onde de plus en plus courtes dans les images des paramètres. Ces dernières sont sensibles au choix des approximations physiques effectuées pour calculer le gradient de la fonction coût (approximation de Born ou de Rytov), au dispositif d'acquisition, au préconditionnement des données et au choix des paramètres inversés. <br /><br />Cette méthode est ensuite appliquée à des milieux contenant une surface libre. La surface libre est une interface très contrastée qui donne lieu à des ondes de surface très énergétiques dans les sismogrammes. Les milieux de subsurface sont complexes, les ondes qui s'y propagent subissent des réflexions ou diffractions multiples. Lorsque le milieu contient deux anomalies dont le contraste en vitesse vaut $20 \%$, l'inversion les localise correctement dans un macromodèle connu. Les images sont améliorées et l'amplitude des anomalies est très bien reconstruite lorsque le nombre de fréquences inversées augmente et les données sont sélectionnées des faibles déports vers les grands déports. Avec une anomalie fortement contrastée comme une cavité vide, l'inversion retrouve correctement la position, la forme et la taille de l'objet mais son amplitude est sous-estimée.<br /><br />L'application à des données réelles verticales de subsurface acquises dans un milieu complexe contenant une cavité maçonnée montre que le milieu hétérogène issu de l'inversion ajuste mieux les données et permet de bien reproduire les ondes inversées. Néanmoins, la cavité n'est pas imagée.

En sciences et dans l'industrie pour limiter le nombre de maquettes expérimentales dans les projets R&D afin de mieux comprendre et de bien interpréter les phénomènes ultrasonores complexes observés sur site, des simulations de contrôles ultrasonores sont effectuées. Ces simulations sont d'autant plus réalistes que la description des structures à contrôler est précise, en particulier au niveau des constantes d'élasticité et d'atténuation intrinsèque. Les objectifs de cette étude sont d'améliorer d'une part les connaissances sur l'influence des caractéristiques métallurgiques des matériaux anisotropes et hétérogènes sur la propagation ultrasonore et d'autre part les performances des codes de calcul (logiciel ATHENA d'EDF) qui nécessitent de disposer des données d'entrée pertinentes, notamment en ce qui concerne les constantes d'élasticité et l'atténuation ultrasonore. Cette étude est dédiée à la caractérisation des matériaux à gros grains, comme les aciers austéno-ferritiques moulés par centrifugation et les soudures en acier inoxydable austénitique ou en alliages à base nickel. Un système expérimental unique permettant de mesurer les constantes d'élasticité et l'atténuation en incidence oblique à été mis au point. Le point fort de ce dispositif est qu'il permet de travailler au-delà de l'angle critique longitudinal et donc de mesurer les propriétés d'atténuation des ondes transversales. Les constantes d'élasticité sont déduites des vitesses ultrasonores à partir d'un processus d'optimisation basé sur la résolution de problème inverse. Nous avons montré les potentialités d'algorithmes d'optimisation globaux tels que les algorithmes génétiques moins susceptibles de converger vers des minima locaux de la fonction à minimiser. Les résultats obtenus à partir des mesures expérimentales sont en accord avec la littérature. Des résultats de l'atténuation des ondes longitudinales et transversales par décomposition du faisceau en spectre d'ondes planes sont présentés.

La modélisation numérique de la propagation des ondes sismiques fait partie des études principales sur le calcul du mouvement sismique basées sur de différents schèmes numériques. La prise en compte du comportement nonlinéaire du sol est consideré désormais très important afin de pouvoir calculer la réponse du milieu cohéremment aux observations sous les sollications sismiques très fortes. En plus, le paramètre de pression interstitielle, qui pourrait emmener le sol aux phénomènes de liquéfaction, devient très important pour les sols saturés. Dans cette étude, dans un premier temps, la propagation des ondes sismiques est modelisée sur une composante (1C) dans les milieux linéaires et nonlinéaires en utilisant la méthode numérique des éléments spéctraux. Les rhéologies viscoélastique et nonlinéaire sont implementées par le méthode de technique des variables de mémoire et le modèle élastoplastique d’Iwan, respectivement. Ensuite, le modèle 1D - trois composantes (3C) est développé et une comparaison préalable sur l’effet de la considération des approches 1C et 3C est faite. L’effet de pression interstitielle est implementé dans le code 1D-3C et le site américain Wildlife Refuge Liquefaction Array (WRLA), qui a été frappé par le séisme de Superstition Hills en 1987 est étudié. Le changement de la réponse du sol sous les différents hypothèses de rhéologie du sol et de mouvement d’entrée est étudié. Le mouvement calculé est noté d’être amplifié pour les basses fréquences et atténué pour les hautes fréquences en raison de l’excès de pression interstitielle dans les sols liquéfiables. Par ailleurs, le sol devient plus nonlinéaire sous le chargement triaxial dans l’approche 3C et plus dilatant dû à la nonlinéarite élevée. En dépit de la similitude entre les accélérations et les vitesses en surface des approches 1C et 3C, une importante différence dans le déplacement en surface entre les deux approches est notée. Les analyses sont répétées pour deux sites japonais Kushiro Port et Onahama Port, qui ont été influencés par le séisme de Kushiro-Oki en 1993 et le séisme de la côte Pacifique de Tohoku en 2011, respectivement. Il a été montré que les changements apportés par la nonlinéarite ne sont pas identiques dans toute la gamme de fréquence concernée et l’influence du comportement des sols non-cohésives sur la propagation des ondes sismiques dépend fortement des propriétés du modèle et des conditions de chargement. Dernièrement, le code SEM est avancé en 2D en considerant les mêmes modèles implementés en 1D-3C pour la nonlinéarite du sol et les effets de pression interstitielle. Le code SEM 2D est mis en application dans un modèle de bassin sédimentaire dont la géometrie est assymmétrique et le profile du sol est composé des couches possédant différentes propriétés nonlinéaires. Le modèle est étudié par les analyses totale et effective pour les propagation des ondes P-SV et SH. La differentiation du mouvement calculé en surface est très importante sous les chargements avec les signaux d’entrée synthétique et réel. L’analyse effective résulte en plus de déformations dans les couches superficielles par rapport à l’analyse totale.De plus, la durée de propagation des ondes augmente à l’intérieur du bassin et les reflections aux frontières de bassin-rocher entraînent plus de nonlinéarite dans les coins du bassin. Cette thèse révèle la possibilité de la modélisation du comportement nonlinéaire du sol en prenant en compte l’effet de pression interstitielle dans les études de la propagation des ondes sismiques en couplant les modèles différents avec la méthode des éléments spéctreaux. Ces analyses contribuent à l’identification et la compréhension des phénomènes majeures qui se déroulent dans les couches superficielles en respectant les conditions locales et les mouvements d’entrées, ce quirend ce travail très important pour les études spécifiques de sites / Numerical modeling of seismic wave propagation has been a major topic on ground motion studies using a number of different numerical integration schemes. The consideration of soil nonlinearity holds an important place in order to achieve simulations consistent with real observations for strong seismic shaking. Additionally, in the presence of strong ground motion in saturated soils, pore pressure becomes an important parameter to take into account for related phenomena such as flow liquefaction and cyclic mobility. In this study, first, one component (1C) - seismic wave propagation is modeled in linear and nonlinear media in 1D based on the spectral element numerical method. Viscoelastic and nonlinear soil rheologies are implemented by use of the memory variables technique and Iwan’s elastoplastic model, respectively. Then, the same study is extended to a 1D - three component (3C) model and a preliminary comparison on the effect of using 1C and 3C approaches is made. Then, the influence of excess pore pressure development is included in the 1D-3C model and the developped numerical model is applied to realistic case on the site of Wildlife Refuge Liquefaction Array (USA) which is affected by the 1987 Superstition Hills event. The ground motion modification for different assumptions of the soil rheology in the media and different input motions is studied. The calculated motion is found to be amplified on low frequency and damped in high frequency range due to excess pore pressure development. Furthermore, the soil is found to be more nonlinear under triaxial loading in 3C approach and more dilative due to higher nonlinearity. Despite the similitude in surface acceleration and velocity results, significant differences in surface displacement results of 1C and 3C approaches are remarked. Similar analyses are performed on two Japanese sites Kushiro Port and Onahama Port, which are influenced by the 1993 Kushiro-Oki and the 2011 off the Pacific coast of Tohoku earthquakes, respectively. It has been shown that the nonlinearity-related changes are not homogeneous all over the concerned frequency band and the influence of cohesionless soil behavior on wave propagation is highly dependent on model properties and loadingconditions. Lastly, the 2D SEM code is developped by taking into account soil nonlinearity and pore pressure effects similary to 1D-3C SEM code. The developped 2D SEM code is applied to a 2D sedimentary basin site where the basin geometry is asymmetrical and soil profile consists of layers with different nonlinearity properties. Total and effective stress analyses are performed on the 2D basin for P-SV and SH zave propagation models. The calculated surface motion is shown to differ significantly under synthetic and realistic input motion loading conditions and the resultant deformation in superficial layers is found to be very high in effective stress analysis compared to total stress analysis. Also, wave propagation takes longer time inside basin media and the reflections on bedrock-basin boundaries lead the soil in basin edges to higher nonlinearity. This study shows the possibility of modeling nonlinear soil behavior including pore pressure effects in seismic wave propagation studies by coupling different models with spectral element method. These analyses help identifying and understanding dominant phenomena occurring in superficial layers, depending on local conditions and input motions. This is of great importance for site-specific studies

Décrit comme l'un des algorithmes les plus prometteurs du 20ème siècle, le formalisme multipolaire appliqué à la méthode des éléments de frontière, permet de nos jours de traiter de larges problèmes encore inconcevables il y a quelques années. La motivation de ce travail de thèse est d'évaluer la capacité, ainsi que les avantages concernant les ressources numériques, de ce formalisme pour apporter une solution de référence aux problèmes de propagation sonore tri-dimensionnels en environnement urbain, dans l'objectif d'améliorer les algorithmes plus rapides déjà existants. Nous présentons la théorie nécessaire à l'obtention de l'équation intégrale de frontière pour la résolution de problèmes non bornés. Nous discutons également de l'équation intégrale de frontière conventionnelle et hyper-singulière pour traiter les artefacts numériques liés aux fréquences fictives, lorsque l'on résout des problèmes extérieurs. Nous présentons par la suite un bref aperçu historique et technique du formalisme multipolaire rapide et des outils mathématiques requis pour représenter la solution élémentaire de l'équation de Helmholtz. Nous décrivons les principales étapes, d'un point de vue numérique, du calcul multipolaire. Un problème de propagation sonore dans un quartier, composé de 5 bâtiments, nous a permis de mettre en évidence des problèmes d'instabilités dans le calcul par récursion des matrices de translations, se traduisant par des discontinuités sur le champs de pression de surface et une non convergence du solveur. Ceci nous a conduits à considérer le travail très récent de Gumerov et Duraiswamy en lien avec un processus récursif stable pour le calcul des coefficients des matrices de rotation. Cette version améliorée a ensuite été testée avec succès sur un cas de multi diffraction jusqu'à une taille dimensionnelle de problème de 207 longueur d'ondes. Nous effectuons finalement une comparaison entre un algorithme d'élément de frontière, Micado3D, un algorithme multipolaire et un algorithme basé sur le tir de rayons, Icare, pour le calcul de niveaux de pression moyennés dans une cour ouverte et fermée. L'algorithme multipolaire permet de valider les résultats obtenus par tir de rayons dans la cour ouverte jusqu'à 300 Hz (i.e. 100 longueur d'ondes), tandis que concernant la cour fermée, zone très sensible par l'absence de contribution directes ou réfléchies, des études complémentaires sur le préconditionnement de la matrice semblent requises afin de s'assurer de la pertinence des résultats obtenus à l'aide de solveurs itératifs / Described as one of the best ten algorithms of the 20th century, the fast multipole formalism applied to the boundary element method allows to handle large problems which were inconceivable only a few years ago. Thus, the motivation of the present work is to assess the ability, as well as the benefits in term of computational resources provided by the application of this formalism to the boundary element method, for solving sound propagation problems and providing reference solutions, in three dimensional dense urban environments, in the aim of assessing or improving fast engineering tools. We first introduce the mathematical background required for the derivation of the boundary integral equation, for solving sound propagation problems in unbounded domains. We discuss the conventional and hyper-singular boundary integral equation to overcome the numerical artifact of fictitious eigen-frequencies, when solving exterior problems. We then make a brief historical and technical overview of the fast multipole principle and introduce the mathematical tools required to expand the elementary solution of the Helmholtz equation and describe the main steps, from a numerical viewpoint, of fast multipole calculations. A sound propagation problem in a city block made of 5 buildings allows us to highlight instabilities in the recursive computation of translation matrices, resulting in discontinuities of the surface pressure and a no convergence of the iterative solver. This observation leads us to consider the very recent work of Gumerov & Duraiswamy, related to a ``stable'' recursive computation of rotation matrices coefficients in the RCR decomposition. This new improved algorithm has been subsequently assessed successfully on a multi scattering problem up to a dimensionless domain size equal to 207 wavelengths. We finally performed comparisons between a BEM algorithm, extit{Micado3D}, the FMBEM algorithm and a ray tracing algorithm, Icare, for the calculation of averaged pressure levels in an opened and closed court yards. The fast multipole algorithm allowed to validate the results computed with Icare in the opened court yard up to 300 Hz corresponding, (i.e. 100 wavelengths), while in the closed court yard, a very sensitive area without direct or reflective fields, further investigations related to the preconditioning seem required to ensure reliable solutions provided by iterative solver based algorithms