Les phénomènes de contact entre les corps, déformables ou non, sont omniprésents dans la vie courante. Leurs modélisations requièrent des outils mathématiques faisant appel à des systèmes d'équations aux dérivées partielles incluant des conditions aux limites non triviales pour décrire le contact. Si les aspects physiques de la mécanique du contact sont connus depuis longtemps, la théorie mathématique qui lui est dédiée reste relativement récente laissant ainsi place à de nombreux problèmes à investiguer. Ce travail porte sur la modélisation, l'analyse et la simulation numérique de tels problèmes. Il se situe à mi-chemin entre la mécanique du contact et les aspects mathématiques inhérents au type de problème qui en découle. L'objectif est ici d'étudier certaines catégories de problèmes faisant intervenir des conditions originales de contact (avec et sans frottement) à la fois d'un point de vue mathématique et numérique, afin d'apporter une contribution à la théorie mathématique, puis de mettre en avant quelques méthodes numériques adaptées à leur résolution dans un cadre spécifique. / Contact phenomena between bodies, whether they are deformable or not, abound in everyday life. Their modellings require mathematical tools using systems of partial differential equations and involving complex boundary conditions, in order to describe the contact. While the physical aspects of such phenomena have been known for a long time, the mathematical theory remains relatively recent which leaves room for numerous problems. This work focuses on the modelling, the analysis and the numerical simulations of such problems. It is located halfway between contact mechanics and the mathematical aspects inherent to the mechanical questions involved. Our aim is to study several groups of problems that include original contact conditions (with or without friction), both from a mathematical and numerical point of view, in order to contribute to the theory, and also to highlight several numerical methods used to solve specific contact problems.

L'objectif de ce travail fut le développement de nouvelles méthodes pour aborder certainsproblèmes inverses en élasticité, en tirant parti de la présence d'un petit paramètre dans ces problèmespour construire des approximation asymptotiques d'ordre élevé.La première partie est consacrée à l'identification de la taille et la position d'une inhomogénéité$BTrue$ enfouie dans un domaine élastique tridimensionnel. Nous nous concentrons sur l'étude defonctions-co^uts $Jbb(Br)$ quantifiant l'écart entre $BTrue$ et une hétérogénéité ``test'' $Br$. Unetelle fonction-co^ut peut en effet être minimisée par rapport à tout ou partie des caractéristiques del'inclusion ``test'' $Br$ (position, taille, propriétés mécaniques ...) pour établir la meilleurecorrespondance possible entre $Br$ et $BTrue$. A cet effet, nous produisons un développement asymptotique de $Jbb$en la taille $incsize$ de $Br$, qui en constitue une approximation polynomiale plus aisée à minimiser. Cedéveloppement, établi jusqu'à l'ordre $O(incsize^6)$, est justifié par une estimation du résidu. Uneméthode d'identification adaptée est ensuite présentée et illustrée par des exemples numériques portant surdes obstacles de formes simples dans l'espace libre $Rbb^3$.L'objet de la seconde partie est de caractériser une inclusion microstructurée de longueur $ltot$, modéliséeen une dimension, composée de couches de deux matériaux alternés périodiquement, en supposant que les plusbasses de ses fréquences propres de transmission (TEs) sont connues. Ces fréquences sont les valeurs propres d'unproblème dit de transmission intérieur (ITP). Afin de disposer d'un modèle propiceà l'inversion, tout en prenant en compte les effets de la microstructure, nous nous reposons sur des approximationsde l'ITP exact obtenues par homogénéisation. A partir du modèle homogénéisé d'ordre 0, nous établissonstout d'abord une méthode simple pour déterminer les paramètres macroscopiques ($ltot$ et contrastes matériaux)d'une telle inclusion. Pour avoir accès à la période de la microstructure, nous nous intéressons ensuite àdes modèles homogénéisés d'ordre élevé, pour lesquels nous soulignons le besoin de conditions aux limitesadaptées. / The purpose of this work was to develop new methods to address inverse problems in elasticity,taking advantage of the presence of a small parameter in the considered problems by means of higher-order asymptoticexpansions.The first part is dedicated to the localization and size identification of a buried inhomogeneity $BTrue$ in a 3Delastic domain. In this goal, we focused on the study of functionals $Jbb(Br)$ quantifying the misfit between $BTrue$and a trial homogeneity $Br$. Such functionals are to be minimized w.r.t. some or all the characteristics of the trialinclusion $Br$ (location, size, mechanical properties ...) to find the best agreement with $BTrue$. To this end, weproduced an expansion of $Jbb$ with respect to the size $incsize$ of $Br$, providing a polynomial approximationeasier to minimize. This expansion, established up to $O(incsize^6)$ in a volume integral equations framework, isjustified by an estimate of the residual. A suited identification procedure is then given and supported by numericalillustrations for simple obstacles in full-space $Rbb^3$.The main purpose of this second part is to characterize a microstructured two-phases layered1D inclusion of length $ltot$, supposing we already know its low-frequency transmission eigenvalues (TEs). Thoseare computed as the eigenvalues of the so-called interior transmission problem (ITP). To provide a convenient invertiblemodel, while accounting for the microstructure effects, we then relied on homogenized approximations of the exact ITPfor the periodic inclusion. Focusing on the leading-order homogenized ITP, we first provide a straightforward method torecover the macroscopic parameters ($ltot$ and material contrast) of such inclusion. To access to the period of themicrostructure, higher-order homogenization is finally addressed, with emphasis on the need for suitable boundaryconditions.

Many problems in continuum mechanics, especially in the theory of elastic materials, lead to nonlinear partial differential equations. The nonconvexity of their underlying energy potential is a challenge for mathematical analysis, since convexity plays an important role in the classical theories of existence and regularity. In the last years one main point of interest was to develop techniques to circumvent these difficulties. One approach was to use different notions of convexity like quasi-- or polyconvexity, but most of the work was done only for static (time independent) equations. In this thesis we want to make some contributions concerning existence, regularity and numerical approximation of nonconvex dynamical problems.

Le thème central de cette thèse est l'étude des résonances pour le problème intérieur en élastodynamique (géométries elliptique et sphéroïdale), et pour le problème de diffusion en acoustique (géométrie elliptique). On s'intéresse en particulier à la levée de dégénérescence des résonances liée à la brisure de symétrie de l'objet lors de la transition du disque circulaire vers le disque elliptique (2D), et de la sphère vers le sphéroïde (3D). Ce phénomène est étudié et interprété d'un point de vue théorique en prenant en compte les symétries de l'objet à l'aide de la théorie des groupes. Cette approche est complétée par une modélisation numérique et une partie expérimentale. En 2D, nous étudions le problème intérieur pour un disque elliptique élastique (étude des modes résonants) et le problème de la diffusion acoustique par des cylindres elliptiques élastiques. Ils sont traités à partir du formalisme modal combiné à la théorie des groupes dans le contexte vectoriel de l'élastodynamique. La levée de dégénérescence est observée théoriquement mais aussi expérimentalement en diffusion. La méthode simplifie considérablement le traitement numérique des problèmes étudiés, fournit une classification des résonances selon les 4 représentations irréductibles du groupe de symétrie C2v (associé à la géométrie elliptique) et donne une interprétation physique de la levée de dégénérescence en termes de brisure de symétrie. Une partie expérimentale en spectroscopie ultrasonore vient compléter l'étude théorique du problème de diffusion. Une série d'expériences en cuve est menée dans le cas de cylindres elliptiques de différentes excentricités en aluminium immergés dans l'eau, dans la bande de fréquence 0 ≤ kr ≤ 50, où kr est le nombre d'onde réduit dans le fluide. Les résultats expérimentaux présentent un très bon accord avec les résultats théoriques, la levée de dégénérescence est observée expérimentalement sur des fonctions de forme et mise en évidence sur des diagrammes angulaires. Le problème intérieur en 3D est traité expérimentalement à partir de la génération et la détection optiques d'ondes élastiques. Une série d'expérimentations sur des objets tridimensionnels (sphère, sphéroïdes oblates et prolates de différentes excentricités) en aluminium est réalisée. Ils sont mis en vibration par impacts laser et les mesures de vitesse et de fréquence s'effectuent par vibrométrie laser. On réalise ainsi une comparaison qualitative entre la théorie 2D et l'expérience 3D. Les mesures sont menées à la fois dans les domaines temporel et fréquentiel pour mettre en évidence la levée de dégénérescence d'une part, et l'onde de Rayleigh qui se propage sur la surface des objets d'autre part. Nous identifions deux trajets pour cette onde en géométrie sphéroïdale, l'un circulaire et l'autre elliptique.Enfin, dans le cadre des problèmes intérieurs 2D et 3D, on donne une interprétation en termes de rayons à travers la dualité entre le spectre des résonances et le spectre des longueurs des orbites périodiques (OPs), avec la mise en évidence du phénomène de conversion de mode et l'identification de l'onde de Rayleigh. Un phénomène, nouveau à notre connaissance, vient s'ajouter au phénomène de bifurcation de certaines orbites. Au cours de la déformation vers le disque elliptique, les orbites avec conversion de mode du disque circulaire se séparent en deux orbites dont les longueurs sont associées aux trajets minimal et maximal qu'elles parcourent. Cette observation s'interprète comme une conséquence du théorème de Fermat. Dans le cas du sphéroïde, on retrouve les orbites du disque circulaire dans le plan équatorial et celles du disque elliptique dans le plan méridien. Nous mettons également en évidence les pics associés aux deux trajets parcourus par l'onde de Rayleigh sur le spectre des OPs. / Resonances for the interior problem in elastodynamics and the acoustic scattering problem are considered in elliptical and spheroidal geometries. Interest is focused on the splitting up of resonances which occurs when the symmetry is broken in the transition from the circular disc to the elliptical one (2D), and from the sphere to the spheroid (3D). From the theoretical point of view, this physical phenomenon is studied and interpreted taking into account the symmetries of the object with the help of group theory. This approach is completed by a numerical modeling and an experimental part. As far as the two dimensional problems are concerned, the interior problem for an elliptical elastic disc (study of resonant modes) and the acoustic scattering problem for infinite elliptical elastic cylinders are studied combining modal formalism and group theory in the vectorial context of elastodynamics. The splitting up of resonances due to the circular symmetry breaking is not only theoretically observed but also experimentally for the scattering problem. The method significantly simplifies the numerical treatment of the problems studied, provides a full classification of resonances over the 4 irreducible representations of the symmetry group C2v (associated with the elliptical geometry) and gives a physical interpretation of the splitting up in terms of symmetry breaking of the symmetry group O(2) (invariance under rotation). An experimental part based on ultrasonic spectroscopy complements the theoretical study of the scattering problem. A series of tank experiments is carried out in the case of aluminum elliptical cylinders immersed in water, in the frequency range 0 ≤ kr ≤ 50, where kr is the reduced wave number in the fluid. The experimental results provide a very good agreement with the theoretical ones, the splitting up is observed on experimental form functions and the split resonant modes are identified on angular diagrams. The interior problem in 3D is studied by means of an experimental approach based on the optical generation and detection of elastic waves. A series of experiments is performed on three-dimensional objects in aluminium. These objects (sphere, prolate and oblate spheroids of various eccentricity) are excited by laser impacts, and the velocity and frequency measurements are carried out by laser vibrometry. Theory and experiments are qualitatively compared. The measurements are performed in both the frequency and time domains to highlight the splitting up of resonances on one hand, and the Rayleigh wave propagating on the surface of the objects on the other hand. We emphasize two paths for this surface wave in spheroidal geometry: a circular one in the equatorial plane and an elliptical one in the meridian plane. Finally, in the context of the interior problems in 2D and 3D, a physical interpretation of resonances in terms of geometrical paths is provided. Mode conversion is highlighted and the Rayleigh wave is identified on the periodic orbits lengths spectrum.In addition to the bifurcations of some periodic orbits, a phenomenon, new to our knowledge, appears. The orbits with mode conversion of the circular disc split in two orbits when the transition to the elliptic disc occurs. The lengths of these orbits are associated with the minimal and maximal travel paths. This observation is interpreted from Fermat's theorem.For the spheroid, orbits of the circular disc and those of the elliptical disc are recovered in the equatorial and meridian planes respectively. We also emphasize the peaks associated with the travel paths of Rayleigh wave in spheroidal geometry appearing on the periodic orbits spectrum.

Le chapitre 1 porte sur une équation des ondes monodimensionnelle soumise à une condition aux limites unilatérale. Sous des hypothèses de régularité appropriées sur les données initiales, une nouvelle preuve d’existence et d’unicité est proposée. La méthode de redistribution de masse qui repose sur une redistribution de la masse d’un corps de telle sorte qu’il n’y ait pas d’inertie au niveau du nœud de contact est introduite et sa convergence est prouvée. Une approximation de ce problème d’évolution combinant la méthode des éléments finis ainsi que la méthode de redistribution de masse est analysée dans le chapitre 2. Puis deux problèmes ainsi que leurs solutions analytiques respectives (l’une étant nouvelle) sont présentés et des discrétisations possibles en utilisant différentes méthodes d’intégration en temps sont décrites. Enfin, des simulations numériques de ces problèmes sont reportées. Dans le chapitre 3, la masse des nœuds de contact est redistribuée sur les autres nœuds et sa convergence ainsi qu’une estimation de l’erreur en temps sont établies. Ensuite, une solution analytique déjà introduite dans le chapitre 3 est comparée aux approximations obtenues pour plusieurs redistributions de masse possibles mettant ainsi en évidence que plus une redistribution de masse d’un corps se fait à proximité des nœuds de contact meilleures sont les solutions approchées obtenues. Les problèmes de contact élastodynamique en dimension d’espace deux et trois sont étudiés dans le chapitre 4. Comme pour les problèmes de contact monodimensionnels, une solution approchée combinant les éléments finis et la redistribution de masse est exposée. Quelques simulations numériques utilisant des méthodes d’intégration en temps mettent en évidence les propriétés de convergence de la méthode de redistribution de masse. / The chapter 1 focuses on a one–dimensional wave equation being subjected to a unilateral boundary condition. Under appropriate regularity assumptions on the initial data, a new proof of existence and uniqueness results is proposed. The mass redistribution method based on a redistribution of the body mass such that there is no inertia at the contact node is introduced and its convergence is proved. An approximation of this evolutionary problem combining the finite element method as well as the mass redistribution method is analyzed in chapter 2. Then two benchmark problems (one being new) with their analytical solutions are presented and some possible discretizations using different time–integration schemes are described. Finally, numerical experiments for these benchmark problems are reported. In chapter 3, the mass of the contact nodes is redistributed on the other nodes and its convergence as well as an error estimate in time are established. Then an analytical solution already introduced in chapter 3 is compared to approximate ones obtained for different choices of mass redistribution highlighting that more a mass redistribution of the body is done near the contact nodes better the approximate solutions are obtained. The two and three–dimensional elastodynamic contact problems are studied in chapter 4. As for the one–dimensional contact problems, an approximated solution combining the finite element and mass redistribution methods is exhibited. Some numerical experiments using time–integration methods highlighted the convergence properties of the mass redistribution method.

Dans cette thèse nous proposons de développer un système d’imagerie ultrasonore innovante de micro- défauts basé sur l’utilisation conjointe de techniques d’acoustique non linéaire et du concept de "transducteur à cavité chaotique". Ce transducteur correspond à la combinaison d’une céramique piézoélectrique collée sur une cavité de forme chaotique et du principe de retournement temporel. La faisabilité et les performances de ce nouveau système sont explorées par des simulations numériques. Des paramètres optimaux d’utilisation pour une implémentation expérimentale sont proposés. Une grande partie des travaux menés dans le cadre de cette thèse se concentre sur le développement d’outils numériques permettant l’amélioration de telles techniques d’imagerie. Un schéma d’éléments finis de type Galerkin Discontinu (GD) est étendu à l’élastodynamique non linéaire. Un type de zone absorbante parfaitement adaptée, appelée "Nearly Perfectly Matched Layer" (NPML) a aussi été développé. Dans le cas de matériaux orthotropes, comme des problèmes de stabilité apparaissent, un mélange de NPML et de zone atténuante, dont on contrôle la proportion respective, est introduit afin de stabiliser les NPML. Une validation expérimentale du concept de "transducteur à cavité chaotique" pour la focalisation dans un milieu solide, réverbérant ou non, en utilisant une seule source est réalisée. Les méthodes de retournement temporel et de filtre inverse sont présentées et comparées. La démonstration expérimentale qu’un "transducteur à cavité chaotique" peut être utilisé conjointement avec les méthodes d’inversion d’impulsion afin de réaliser une image de non linéarités localisées est présentée / In this thesis we propose the development of an innovative micro-damage imaging system based on a combination of Nonlinear Elastic Wave Spectroscopy techniques and “chaotic cavity transducer” concept. It consists of a combination of a PZT ceramic glued to a cavity of chaotic shape with the time reversal principle. The feasibility and capabilities of these new ideas is explored by numerical simulations, and optimal operational parameters for experimental implementation are suggested based on the modelling support. A large part of the research work conducted in this thesis is concentrated on the development of numerical simulation tools to help the improvement of such nonlinear imaging methods. A nodal Discontinuous Galerkin Finite Element Method (DG-FEM) scheme is extended to nonlinear elasto-dynamic including source terms. A Perfectly Matched Layer absorbing boundary condition well adapted to the DG-FEM scheme, called Nearly Perfectly Matched Layer (NPML), is also developed. In the case of orthotropic material as stability problems appear, a mixture of NPML and sponge layer, with a controllable ratio of these two kinds of absorbing layers, is introduced. The experimental validation of “chaotic cavity transducer” to focalize in reverberant and non-reverberant solid media with only one source is made. Classical time reversal, inverse filter and 1 Bit time reversal process are discussed and compared. The experimental demonstration of the use of a “chaotic cavity transducer”, in combination with the pulse inversion and 1-bit methods, to obtain an image of localized nonlinearity is made. This opens the possibility for high resolution imaging of nonlinear defects

Ce travail porte sur l'étude théorique et numérique des cristaux phononiques non linéaires. Les non linéarités étudiées sont celles dues aux constantes élastiques d'ordre deux (quadratiques) et trois (cubiques) des matériaux constituant les cristaux. Les effets non linéaires sont étudiés grâce à des méthodes d'éléments finis en simulant la propagation d'une onde élastique à travers les cristaux.Un premier projet de recherche a porté sur l'étude d'une structure osseuse, et plus spécifiquement sur la dispersion des ondes élastiques dans une structure constituée d'une alternance de couches de collagène et d'hydroxy apatite. Les simulations montrent qu'il existe un lien étroit entre l'hydratation des os et leur capacité à dissiper l'énergie.La seconde étude réalisée concerne un résonateur élastique. Une structure constituée d'inclusions d'acier dans de la silice présente un comportement de commutateur lorsque les non linéarités cubiques de l'acier sont prises en compte. Cet effet fortement non linéaire apparaît lorsque l'amplitude de l'onde incidente dépasse un certain seuil. Un modèle analytique complet est fourni.La dernière étude réalisée montre la conception de matériaux composites possédant de fortes non linéarités cubiques mais de faibles non linéarités quadratiques. La dérivation des lois de mélange des paramètres élastiques d'un matériau non linéaire dans un matériau linéaire est effectuée à l'ordre trois. Les équations montrent une forte amplification des paramètres non linéaires du matériau résultant pour certaines concentrations. Les simulations permettent de conclure que le résonateur mentionné ci-dessus peut effectivement être réalisé. / This work is dedicated to the theoretical and numerical study of nonlinear phononic crystals. The studied nonlinearities are those due to the second (quadratic) and third (cubic) order elastic constants of the materials that constitute the crystals. Nonlinear effects are studied by the means of finite element methods, used to simulate the propagation of an elastic wave through the crystals.A first research project concerns the study of a bone structure, namely the dispersion of elastic waves in a structure composed of collagen and hydroxy apatite alternate constituent layers. Simulations showed that it exists a strong link between bones hydration and their ability to dissipate the energy.The second study relates to an elastic resonator. A structure composed of steel inclusions in a silica matrix shows a switch behavior when the cubic nonlinearities of steel are taken into account. This strong nonlinear effect appears when the amplitude of the incident wave reaches a threshold. A full analytical model is provided.The last study demonstrates the design of composite materials with both strong cubic nonlinearities and weak quadratic nonlinearities. The derivation of the mixing laws of the elastic parameters of a nonlinear material inside a linear one is performed up to order three. Equations show a strong amplification of the nonlinear parameters of the material for some concentrations. Numerical simulations allow to conclude that the above mentioned resonator can be produced.

Cette thèse porte sur l'étude théorique et numérique des méthodes rapides pour résoudre les équations de l'élastodynamique 3D en domaine fréquentiel, et se place dans le cadre d'une collaboration avec la société Shell en vue d'optimiser la convergence des problèmes d'inversion sismique. La méthode repose sur l'utilisation des éléments finis de frontière (BEM) pour la discrétisation et sur les techniques de matrices hiérarchiques (H-matrices) pour l'accélération de la résolution du système linéaire. Dans le cadre de cette thèse on a développé un solveur direct pour les BEMs en utilisant une factorisation LU et un stockage hiérarchique. Si le concept des H-matrices est simple à comprendre, sa mise en oeuvre requiert des développements algorithmiques importants tels que la gestion de la multiplication de matrices représentées par des structures différentes (compressées ou non) qui ne comprend pas mois de 27 sous-cas. Un autre point délicat est l'utilisation des méthodes d'approximations par matrices compressées (de rang faible) dans le cadre des problèmes vectoriels. Une étude algorithmique a donc été faite pour mettre en oeuvre la méthode des H-matrices. Nous avons par ailleurs estimé théoriquement le rang faible attendu pour les noyaux oscillants, ce qui constitue une nouveauté, et montré que la méthode est utilisable en élastodynamique. En outre on a étudié l'influence des divers paramètres de la méthode en acoustique et en élastodynamique 3D, à fin de calibrer leur valeurs numériques optimales. Dans le cadre de la collaboration avec Shell, un cas test spécifique a été étudié. Il s'agit d'un problème de propagation d'une onde sismique dans un demi-espace élastique soumis à une force ponctuelle en surface. Enfin le solveur direct développé a été intégré au code COFFEE développé a POEMS (environ 25000 lignes en Fortran 90) / This thesis focuses on the theoretical and numerical study of fast methods to solve the equations of 3D elastodynamics in frequency-domain. We use the Boundary Element Method (BEM) as discretization technique, in association with the hierarchical matrices (H-matrices) technique for the fast solution of the resulting linear system. The BEM is based on a boundary integral formulation which requires the discretization of the only domain boundaries. Thus, this method is well suited to treat seismic wave propagation problems. A major drawback of classical BEM is that it results in dense matrices, which leads to high memory requirement (O (N 2 ), if N is the number of degrees of freedom) and computational costs.Therefore, the simulation of realistic problems is limited by the number of degrees of freedom. Several fast BEMs have been developed to improve the computational efficiency. We propose a fast H-matrix based direct BEM solver.

The first part of this Thesis deals with the analysis of piecewise-smooth mechanical systems and the definition of special stability criteria in presence of non-conservative follower forces. To illustrate the peculiar stability properties of this kind of dynamical system, a reference 2 d.o.f. structure has been considered, composed of a rigid bar, with one and constrained to slide, without friction, along a curved profile, whereas the other and is subject to a follower force. In particular, the curved constraint is assumed to be composed of two circular profiles, with different and opposite curvatures, defining two separated subsystems. Due to this jump in the curvature, located at the junction point between the curved profiles, the entire mechanical structure can be modelled by discontinuous equations of motion, the differential equations valid in each subsystem can be combined, leading to the definition of a piecewise-smooth dynamical system. When a follower force acts on the structure, an unexpected and counterintuitive behaviour may occur: although the two subsystems are stable when analysed separately, the composed structure is unstable and exhibits flutter-like exponentially-growing oscillations. This special form of instability, previously known only from a mathematical point of view, has been analysed in depth from an engineering perspective, thus finding a mechanical interpretation based on the concept of non-conservative follower load. Moreover, the goal of this work is also the definition of some stability criteria that may help the design of these mechanical piecewise-smooth systems, since classical theorems cannot be used for the investigation of equilibrium configurations located at the discontinuity. In the literature, this unusual behaviour has been explained, from a mathematical perspective, through the existence of a discontinuous invariant cone in the phase space. For this reason, starting from the mechanical system described above, the existence of invariant cones in 2 d.o.f. mechanical systems is investigated through Poincaré maps. A complete theoretical analysis on piecewise-smooth dynamical systems is presented and special mathematical properties have been discovered, valid for generic 2~d.o.f. piecewise-smooth mechanical systems, which are useful for the characterisation of the stability of the equilibrium configurations. Numerical tools are implemented for the analysis of a 2~d.o.f. piecewise-smooth mechanical system, valid for piecewise-linear cases and extendible to the nonlinear ones. A numerical code has been developed, with the aim of predicting the stability of a piecewise-linear dynamical system a priori, varying the mechanical parameters. Moreover, “design maps” are produced for a given subset of the parameters space, so that a system with a desired stable or unstable behaviour can easily be designed. The aforementioned results can find applications in soft actuation or energy harvesting. In particular, in systems devoted to exploiting the flutter-like instability, the range of design parameters can be extended by using piecewise-smooth instead of smooth structures, since unstable flutter-like behaviour is possible also when each subsystem is actually stable. The second part of this Thesis deals with the numerical analysis of an elastic cloak for transient flexural waves in Kirchhoff-Love plates and the design of special metamaterials for this goal. In the literature, relevant applications of transformation elastodynamics have revealed that flexural waves in thin elastic plates can be diverted and channelled, with the aim of shielding a given region of the ambient space. However, the theoretical transformations which define the elastic properties of this “invisibility cloak” lead to the presence of a strong compressive prestress, which may be unfeasible for real applications. Moreover, this theoretical cloak must present, at the same time, high bending stiffness and a null twisting rigidity. In this Thesis, an orthotropic meta-structural plate is proposed as an approximated elastic cloak and the presence of the prestress has been neglected in order to be closer to a realistic design. With the aim of estimating the performance of this approximated cloak, a Finite Element code is implemented, based on a sub-parametric technique. The tool allows the investigation of the sensitivity of specific stiffness parameters that may be difficult to match in a real cloak design. Moreover, the Finite Element code is extended to investigate a meta-plate interacting with a Winkler foundation, to analyse how the substrate modulus transforms in the cloak region. This second topic of the Thesis may find applications in the realization of approximated invisibility cloaks, which can be employed to reduce the destructive effects of earthquakes on civil structures or to shield mechanical components from unwanted vibrations.