No description available.

On étudie quelques problèmes singuliers d'homogénéisation en élasticité linéarisée ou électricté non-linéaire où, en plus de la période de distribution des hétérogénéités, apparaissent d'autres paramètres comme la très forte (ou faible) rigidité ou permittivité d'une des phases et la taille relative de celle-ci. Sont visés des milieux fibrés ou stratifiés et des maçonneries planes et minces. / We study some singular problems of homogenization in linearized elasticity or nonlinear electricty where, in addition to the period of distribution of heterogenities, other parameters do appear as the very high (or low) stiffness or permittivity of one phase and its the relative size. Fibered media, stratified media and flat and thin masonries are considered.

We are concerned with forward wave motion simulations in two-dimensional elastic, heterogeneous, semi-infinite media. We use Perfectly Matched Layers (PMLs) to truncate the semi-infinite extent of the physical domain to arrive at a finite computational domain. We use a recently developed hybrid formulation, where the Navier equations for the interior domain are coupled with a mixed formulation for an unsplit-field PML. Here, we implement the hybrid formulation using spectral elements, and report on its performance. The motivation stems from the following considerations: Of concern is the long-time instability that has been reported even in homogeneous and isotropic cases, when the standard complex-stretching function is used in the PML. The onset of the instability is always within the PML zone, and it manifests as error growth in time. It has been suggested that the instability arises when waves impinge at grazing angle on the PML-interior domain interface. Yet, the instability does not always appear. Furthermore, different values of the various PML parameters (mesh density, attenuation strength, order of attenuation function, etc) can either hinder or delay the onset of the instability. It is thus conjectured that the instability is associated with the spectral properties of the discrete operators. In this thesis, we report numerical results based on both Lagrange interpolants, and results based on spectral elements. Spectral elements are explored since they lead to diagonal mass matrices, have improved dispersion error, and, more importantly, have different spectral properties than Lagrangian-based finite elements. Spectral elements are thus used in an attempt to explore whether the reported instability issues could be alleviated. We design numerical experiments involving explosive sources situated at varying depths from the surface, capable of inducing grazing-angle waves. We use the energy decay as the primary metric for reporting the results of comparisons between various spectral element orders and classical Lagrange interpolants. We also report the results of parametric studies. Overall, it is shown that the spectral elements alone are not capable of removing the instability, though, on occasion, they can. Careful parameterization of the PML could also either remove it or alleviate it. The issue remains open. / text

We introduce a mathematical framework for the inverse medium problem arising commonly in geotechnical site characterization and geophysical probing applications, when stress waves are used to probe the material composition of the interrogated medium. Specifically, we attempt to recover the spatial distribution of Lame's parameters ( and μ) of an elastic semi-infinite arbitrarily heterogeneous medium, using surface measurements of the medium's response to prescribed dynamic excitations. The focus is on characterizing near-surface deposits, and to this end, we develop a method that is implemented directly in the time-domain, is driven by the full waveform response collected at receivers on the surface, while the domain of interest is truncated using Perfectly-Matched-Layers (PMLs) to limit the originally semi-infinite extent of the physical domain. There are two key issues associated with the problem at hand: (a) the forward problem, namely the numerical simulation of the wave motion in the domain of interest; and (b) the framework and strategies for tackling the inverse problem. To address the forward problem, it is necessary that the domain of interest be truncated, and the resulting finite domain be forced to mimic the physics of the original problem: to this end, we introduce unsplit-field PMLs, and develop and implement two new formulations, one fully-mixed and one hybrid (mixed coupled with a non-mixed approach) that model wave motion within the, now PML-truncated, domain. To address the inverse problem, we adopt a partial-differential-equation-constrained optimization framework that results in the usual triplet of an initial-and-boundary-value forward problem, a final-and-boundary-value adjoint problem, and a time-independent boundary-value control problem. This triplet of boundary-value-problems is used to guide the optimizer to the target profile of the spatially distributed Lame parameters. Given the multiplicity of solutions, we assist the optimizer, by deploying regularization schemes, continuation schemes (regularization factor and source-frequency content), as well as a physics-driven simple procedure to bias the search directions. We report numerical examples attesting to the quality, stability, and efficiency of the forward wave modeling. We also report moderate success with numerical experiments targeting inversion of both smooth and sharp profiles in two dimensions. / text

As many of the world’s nuclear power plants are near the end of their supposed life span a need arise to assess the components crucial to the safety of these plants. One of these crucial components is the concrete reactor confinement; to assess its condition, non-destructive testing (NDT) is an attractive method. Traditional testing of concrete structures has comprised of drilling out a sample and performing stress tests on it, but because of the radioactive environment inside the containment this method is far from ideal. NDT is of course possible to use at any structure but at reactor containments the benefits from not creating holes in the structure are prominent; NDT is also an attractive option from an esthetical point of view because it leaves the structure intact. The NDT method pertaining to this study is the impact echo method which comprise of applying a force on the structure, usually a hammer blow, and measuring the response with a receiver. The impact will excite waves propagating in the structure which gives rise to Lamb modes. Lamb modes are structural oscillations of the wall and it is the frequency of these modes that are used to determine the thickness of the wall. The elastic properties of the structure can in turn be obtained by measuring the velocities of the waves propagation. It is also possible to use the impact echo method to detect irregularities in the structure such as cracks or delamination. To simulate the dynamics of a system using NDT numerical methods such as finite element modeling (FEM) is often used. The purpose of this study is to assess the possibility to utilize absorbing layers using increasing damping (ALID) in models to reduce the computational time of FEM analyses. ALIDs are used at the edges to simulate an infinite system and are thus supposed to cancel out incoming waves to prevent unwanted reflection from the edges. The models in this study have all pertained to two dimensional plates utilizing infinitesimal strain theory; the decrease in computational time is significant when using ALIDs and for three dimensional models it would be even more so. The ALIDs are specified by length and maximum mass proportional Rayleigh damping (CMmax), in this study three different lengths are tested, 0.5, 1.5 and 4.5 m for CMmax ranging from 103 to 2*105 Ns/m. The damping is increased with increasing distance into the ALID with specified maximum value at the back edge. However, it should be noted that the increase in damping causes difference in impedance between elements and if this difference is too large it will cause reflections of waves at the boundary between the elements. The ALID must thus be defined so that it sufficiently cancels out the wave without causing unwanted reflections due to impedance differences. The conclusion is that the 0.5 m long ALID does not provide good results for any choice of maximum mass proportional Rayleigh damping. Both the 1.5 and 4.5 m long ALIDs are, however, concluded to be applicable; the 1.5 m ALID having 2*104 < CMmax <5*104 Ns/m and the 4.5 m ALID having 5*103 < CMmax < 104 Ns/m are choices that have shown promise in the performed simulations. The hope is that the results obtained in this study will aid in the development of numerical analysis techniques for NDT methods that can be used in the construction of new reactor confinements and/or maintenance of existing reactor confinements and other thick concrete structures. / Många av världens kärnkraftverk närmar sig slutet på sin beräknade livslängd och ett behov uppstår då att kunna utvärdera de komponenter som är väsentliga för säkerheten på dessa verk. Reaktoromslutningen i betong är en av dessa komponenter och oförstörande provning (NDT) är en attraktiv metod för att bedöma dess tillstånd. Traditionellt har utvärdering av betongkonstruktioner bestått av stresstester på borrprover men p.g.a. den radioaktiva miljön på insidan av omslutningen är denna metod ej att föredra. NDT är självklart möjligt att använda på allsköns betongkonstruktioner då det ger både konstruktionsmässiga och estetiska fördelar. NDT metoden som rör denna studie kallas impact echo och går ut på att man med en hammare slår till en punkt på väggen och mäter responsen en bit därifrån. Lasten ger upphov till vågor i form av deformation som propagerar i väggen och dessa ger i sin tur upphov till Lamb moder. Lamb moderna är strukturella oscillationer av väggen och genom att studera dess frekvenser kan väggens tjocklek bestämmas. Elastiska egenskaper i väggen erhålls utifrån de olika vågornas propageringshastigheter. Impact echo metoden kan även användas för att finna strukturella oegentligheter inuti väggen så som sprickor och delaminering. För att utföra numeriska simuleringar av dynamiska system med NDT-metoder är finita elementmetoden (FEM) användbar. Syftet med denna studie är att bedöma vilka möjligheter som finns för att implementera absorberande ränder med ökande dämpning (ALID) i datamodeller för att minska beräkningstiden av FEM-analyser. ALID används vid kanterna för att simulera ett oändligt system, dess uppgift är att dämpa bort inkommande vågor så att dessa ej reflekteras tillbaka och stör mätningarna. Samtliga modeller i denna studie är två-dimensionella med antagen oändligt liten spänning i normalriktningen. Vinsten i beräkningstid av att använda ALID är stor och ökar ytterligare om modellen utökas till tre dimensioner. Ett ALID definieras genom dess längd och maximala massproportionerlig Rayleigh-dämpning (CMmax). I denna rapport har längderna 0.5, 1.5 and 4.5 m använts med CMmax i intervallet från 103 till 2*105 Ns/m. Dämpningen ökar med ökat avstånd in i ALID med det specificerade maxvärdet vid den bakre kanten. Det bör noteras att skillnad i dämpning mellan element leder till skillnad i impedans; reflektioner av vågorna uppstår vid övergång från ett element med lägre impedans till ett med högre impedans. Ett ALID måste således vara definierat så att det dämpar bort tillräckligt av de inkommande vågorna utan att oönskade reflektioner i ALID uppstår. Studien pekar på att ett 0.5 m långt ALID inte åstadkommer önskvärda resultat för något av valen för CMmax som använts i denna rapport. Både det 1.5 och 4.5 m långa ALID har däremot get bra resultat; ett 1.5 m långt ALID bör ha 2*104 < CMmax <5*104 Ns/m och ett 4.5 m långt ALID 5*103 < CMmax < 104 Ns/m. Förhoppningen med studien är att resultaten skall underlätta utvecklingen av NDT-metoder som kan användas vid konstruktion och underhåll av reaktoromslutningar och andra tjocka betongkonstruktioner.

L'identification passive d'un système réside dans l'estimation des paramètres qui décrivent ce système uniquement à l'aide de sollicitations ambiantes. Dans le génie civil, cette discipline est appliquée pour le suivi de l'état de santé des structures, on parle de SHM (Structural Health Monitoring) passif. Le SHM passif est généralement réalisé à l'aide d'une instrumentation déployée en surface. La thèse s'est intéressée aux possibilités offertes par une instrumentation qui serait enfouie. Dans une première partie, on établit les résultats associés à l'identification passive des milieux visco-élastiques. L'originalité de ces travaux réside dans la prise en compte d'un modèle de dissipation réaliste, la viscosité, ainsi que du caractère vectoriel des ondes élastiques. Ces résultats théoriques sont validés expérimentalement et démontrent la portabilité du SHM passif en surface au SHM passif en volume. Dans une deuxième partie, on s'intéresse à deux problèmes attachés à l'enfouissement de capteurs: l'estimation passive de leur position (problème SNL, Sensor Network Location problem) et de leur attitude (problème SNA, Sensor Network Attitude problem). Ces problèmes sont résolus grâce à l'identification passive qui fournit, en plus d'information physique sur le milieu, des informations géométriques sur le réseau. En particulier, on peut estimer des distances et des attitudes relatives entre capteurs. A l'aide de ces informations partielles et bruitées, des algorithmes de résolution des problèmes SNL et SNA ont été développés puis validés expérimentalement. Enfin, on synthétise l'apport de la thèse et on identifie les verrous technologiques à lever afin de justifier la faisabilité de l'enfouissement d'un réseau de capteurs dans le but de faire du SHM passif. / Passive identification of a system relies on the estimation of the parameters which describe that system only by using ambient sources. In civil engineering, we can apply this technique to monitor the state of health of structures. This is called passive SHM (Structural Health Monitoring). Passive SHM is generally realised by using an instrumentation distributed on the surface. This thesis focuses on the possibility given by the use of an embedded instrumentation.In the first part, we establish new results associated to passive identification in visco-elastic media. The originality of this work relies on the consideration of a realistic dissipation model, the viscosity, and the vectorial aspect of elastic waves. Those theoretical results, which are experimentally validated, prove the portability of surface passive SHM to volume passive SHM. In the second part, we focus on two problems related to an embedded sensors network: the passive estimation of sensors position (SNL problem - Sensor Network Location problem) and attitude (SNA problem - Sensor Network Attitude problem). Those problems are solved by using passive identification which gives, besides physical information on the medium, geometrical information on the network. In particular, we can estimate the distances and relative attitudes between sensors. With that partial and noisy information we have developed algorithms solving SNL and SNA problems and we have validated them experimentally.At last, we synthesize the contribution of the thesis and we identify the technological locks to release in order to justify the feasibility of passive SHM using an embedded instrumentation.