Un système SHM par ondes guidées a pour but d'évaluer l'intégrité d'une grande variété de structures fines, telles que les fuselages d'avions, les tuyaux, les réservoirs, etc. Un tel système est basé sur un réseau de capteurs piézoélectriques pour l'excitation et la mesure des ondes guidées. Cette méthode de SHM par ondes guidées est prometteuse pour l'inspection de structures de grande taille, ces ondes se propageant sur de grandes distances avec peu d'atténuation, tout en étant sensibles aux défauts surfaciques et subsurfaciques.Cette thèse présente les travaux menés dans le but de développer un système de SHM par ondes guidées capable de détecter, localiser et dimensionner efficacement les défauts dans des structures aéronautiques assimilables à des plaques, en matériaux composites ou en aluminium. Des simulations et des méthodes d'apprentissage sont utilisées pour déterminer les caractéristiques principales des ondes guidées propagées, notamment les vitesses de phase et de groupe ainsi que la fonction de Green 3D. Celles-ci sont ensuite utilisées pour traiter les signaux des ondes guidées afin de produire des images représentant l'intégrité des structures étudiées. Ce travail comprend également une étude approfondie des algorithmes d'imagerie DAS, MV et Excitelet, les plus prometteurs parmi ceux de la littérature, une évaluation de leurs performances par analyse statistique sur une grande base de données de résultats de simulations d'imagerie par ondes guidées et propose une méthode d'imagerie parcimonieuse.Alors que la détection et la localisation des défauts à partir de l'analyse des images est aisée, le dimensionnement du défaut est un problème plus complexe en raison de sa forte dimensionnalité et de sa non-linéarité. Il est démontré que ce problème peut être résolu par des méthodes d'apprentissage automatique sur une grande base de données de résultats de simulations d'imagerie par ondes guidées. Ces méthodes d'imagerie nécessitent cependant une référence, mesurée sur la structure dans un état supposé sain. Elles sont efficaces dans des conditions opérationnelles stationnaires mais sont sensibles aux variations de l'environnement et notamment aux fluctuations de température.Ce travail présente donc l'étude de la robustesse face aux effets thermiques des méthodes d'imagerie par ondes guidées et propose un modèle de détection de défauts capable d'analyser des résultats d'imagerie détériorés. Plusieurs techniques de compensation des effets thermiques sont étudiées et des améliorations sont proposées. Leur efficacité est validée pour les plaques d'aluminium mais des améliorations supplémentaires sont nécessaires pour les étendre aux plaques de composites. / A guided wave-based structural health monitoring (SHM) system aims at determining the integrity of a wide variety of plate-like structures, including aircraft fuselages, pipes, tanks etc. It relies on a sparse array of piezoelectric transducers for guided waves (GWs) excitation and sensing. With a number of benefits, these waves are standing out among other methods as a promising method for the inspection of large structures. They can propagate on significant distances with small attenuation while being sensitive to surface and subsurface defects.This thesis presents studies conducted with the purpose of developing such a GWs-based SHM system that is capable of efficient defect detection, localization and sizing aeronautical plate-like structures made of aluminum and composite materials. Simulation and data-driven approaches are presented for determining principal characteristics of propagating GWs, namely modal group and phase velocities, 3D Green's functions etc. in structures of interest. They are then used for GWs signals processing in order to compute images representing the integrity of studied structures. This work also provides a comprehensive overview of DAS, MV and Excitelet defect imaging algorithms, determines their performance using statistical analysis of an extensive dataset of simulated guided waves imaging (GWI) results and proposes a method for sparse defect imaging.While defect detection and localization are straightforward from the image analysis, the defect sizing is a more complex problem due to its high dimensionality and non-linearity. It is demonstrated that this problem can be solved by means of machine learning methods, relying on an extensive database of simulated GWI results. Aforementioned defect imaging methods are baseline demanding. They are efficient under stationary operational conditions but vulnerable to environmental variations, especially to the temperature fluctuation.Finally, this work presents studies on the robustness of GWI methods against thermal effects, and a defect detection model capable of analyzing deteriorated GWI results is proposed. Different techniques for thermal effects compensation are reviewed, and improvements are proposed. Their effectiveness is validated for aluminum plates but further improvements are required to translate these techniques to composite plates.

L'objectif de ce travail de thèse est de prendre en compte à la fois l'hétérogénéité, l'anisotropie et des incertitudes dans la simulation 3D de la propagation d'ondes élastiques. Pour ce faire, dans un premier temps, on modélise le champ de propriétés mécaniques, ici le champ de tenseur d'élasticité, par un modèle de champ stochastique 3D des matrices définie-positives. La construction de ce modèle de champ est essentiellement fondée sur celle de Soize [2008]. Notre modèle conserve ainsi les propriétés principales du modèle de Soize comme le paramétrage minimal contrôlant l'amplitude de la fluctuation et la taille caractéristique de la variabilité spatiale, le comportement local a priori arbitrairement anisotrope (anisotropie triclinique) et les propriétés mathématiques fondamentales. De plus, un nouveau paramètre est introduit dans ce modèle pour imposer un niveau d'anisotropie moyen souhaité. Dans un deuxième temps, on effectue des adaptations du code de calcul d'éléments finis spectraux, à savoir le code parallèle SPEC3D, afin d'une part de générer les réalisations du champ stochastique du tenseur d'élasticité et d'autre part de prendre en compte l'anisotropie dans la résolution numérique du problème élastodynamique. Des études paramétriques utilisant SPEC3D sont ensuite réalisées mettant en évidence les influences de l'anisotropie et des paramètres d'hétérogénéité sur la propagation d'ondes sismiques. En particulier, elles démontrent une dépendance directe entre la longueur de corrélation du champ de propriétés et le temps caractéristique d'apparition de la diffusion. Ce régime se manifeste par l'équipartition d'énergie entre les mouvements irrotationnels et rotationnels.

L'objectif de ce travail de thèse est de prendre en compte à la fois l'hétérogénéité, l'anisotropie et des incertitudes dans la simulation 3D de la propagation d'ondes élastiques. Pour ce faire, dans un premier temps, on modélise le champ de propriétés mécaniques, ici le champ de tenseur d'élasticité, par un modèle de champ stochastique 3D des matrices définie-positives. La construction de ce modèle de champ est essentiellement fondée sur celle de Soize [2008]. Notre modèle conserve ainsi les propriétés principales du modèle de Soize comme le paramétrage minimal contrôlant l'amplitude de la fluctuation et la taille caractéristique de la variabilité spatiale, le comportement local a priori arbitrairement anisotrope (anisotropie triclinique) et les propriétés mathématiques fondamentales. De plus, un nouveau paramètre est introduit dans ce modèle pour imposer un niveau d'anisotropie moyen souhaité. Dans un deuxième temps, on effectue des adaptations du code de calcul d'éléments finis spectraux, à savoir le code parallèle SPEC3D, afin d'une part de générer les réalisations du champ stochastique du tenseur d'élasticité et d'autre part de prendre en compte l'anisotropie dans la résolution numérique du problème élastodynamique. Des études paramétriques utilisant SPEC3D sont ensuite réalisées mettant en évidence les influences de l'anisotropie et des paramètres d'hétérogénéité sur la propagation d'ondes sismiques. En particulier, elles démontrent une dépendance directe entre la longueur de corrélation du champ de propriétés et le temps caractéristique d'apparition de la diffusion. Ce régime se manifeste par l'équipartition d'énergie entre les mouvements irrotationnels et rotationnels.

Les travaux de la présente thèse portent sur l’estimation d'erreur a posteriori pour les solutions numériques par éléments finis de l'équation des ondes élastiques dans les milieux hétérogènes. Deux types d’estimation ont été développés. Le premier considère directement l’équation élastodynamique et conduit à un nouvel estimateur d'erreur a posteriori explicite en norme L∞ en temps. Les principales caractéristiques de cet estimateur explicite sont l'utilisation de la méthode de résidus et le développement de reconstructions en temps et en espace selon les différentes régularités exigées par les différents termes contribuant à l’obtention d’une borne supérieure. L’analyse numérique de cet estimateur dans le cas des maillages uniformes montre qu’il assure bien une borne supérieure mais avec une propriété asymptotique qui reste à améliorer. Le deuxième type d’estimateur d’erreur est développé dans le contexte de la propagation des ondes à haute fréquence dans des milieux hétérogènes à l’échelle mésoscopique. Il s’agit d’une nouvelle erreur en résidus basée sur l'équation de transfert radiatif, qui est obtenue par un développement asymptotique multi-échelle de l'équation d'onde en utilisant la transformation de Wigner en espace-temps. Les résidus sont exprimés en termes de densités énergétiques calculés dans l’espace des phases pour les solutions d’onde numériques transitoires par éléments finis. L’analyse numérique de cette erreur appliquée aux milieux homogènes et hétérogènes en 1D a permis de valider notre approche. Les champs d’application visés sont la propagation des ondes sismiques dans les milieux géophysiques ou la propagation des ondes ultrasonores dans les milieux polycristallins. / This thesis work deals with a posteriori error estimates for finite element solutions of the elastic wave equation in heterogeneous media. Two different a posteriori estimation approaches are developed. The first one, in a classical way, considers directly the elastodynamic equation and results in a new explicit error estimator in a non-natural L∞ norm in time. Its key features are the use of the residual method and the development of space and time reconstructions with respect to regularities required by different residual operators contributing to the proposed error bound. Numerical applications of the error bound with different mesh sizes show that it gives rise to a fully computable upper bound. However, its effectivity index and its asymptotic accuracy remain to be improved. The second error estimator is derived for high frequency wave propagation problem in heterogeneous media in the weak coupling regime. It is a new residual-type error based on the radiative transfer equation, which is derived by a multi-scale asymptotic expansion of the wave equation in terms of the spatio-temporal Wigner transforms of wave fields. The residual errors are in terms of angularly resolved energy quantities of numerical solutions of waves by finite element method. Numerical calculations of the defined errors in 1D homogeneous and heterogeneous media allow validating the proposed error estimation approach. The application field of this work is the numerical modelling of the seismic wave propagation in geophysical media or the ultrasonic wave propagation in polycrystalline materials.

La méthode d’imagerie la plus répandue dans l’industrie pétrolière est la RTM (Reverse Time Migration) qui repose sur la simulation de la propagation des ondes dans le sous-sol. Nous nous sommes concentrés sur un propagateur d'ondes élastiques 3D en milieu anisotrope de type TTI (Tilted Transverse Isotropic). Nous avons directement travaillé dans le code de recherche de Total DIVA (Depth Imaging Velocity Analysis), basé sur une discrétisation par la méthode de Galerkin Discontinue et le schéma Leap-Frog, et développé pour le calcul parallèle intensif – HPC (High Performance Computing). Nous avons ciblé plus particulièrement deux contributions possibles qui, si elles supposent des compétences très différentes, ont la même finalité : réduire les coûts de calculs requis pour la simulation. D'une part, les conditions aux limites classiques de type PML (Perfectly Matched Layers) ne sont pas stables dans des milieux TTI. Nous avons proposé de formuler une CLA (Conditions aux Limites Absorbantes) stable dans des milieux anisotropes. La méthode de construction repose sur les propriétés des courbes de lenteur, ce qui donne à notre approche un caractère original. D'autre part, le parallélisme initial, basé sur une décomposition de domaine et des communications par passage de messages à l'aide de la bibliothèque MPI, conduit à un déséquilibrage de charge qui détériore son efficacité parallèle. Nous avons corrigé cela en remplaçant le paradigme parallélisme par l'utilisation de la programmation à base de tâches sur support d'exécution. Cette thèse a été réalisée dans le cadre de l'action de recherche DIP (Depth Imaging Partnership) qui lie la compagnie pétrolière Total et Inria. / The most common method of Seismic Imaging is the RTM (Reverse Time Migration) which depends on wave propagation simulations in the subsurface. We focused on a 3D elastic wave propagator in anisotropic media, more precisely TTI (Tilted Transverse Isotropic). We directly worked in the Total code DIVA (Depth Imaging Velocity Analysis) which is based on a discretization by the Discontinuous Galerkin method and the Leap-Frog scheme, and developed for intensive parallel computing – HPC (High Performance Computing). We choose to especially target two contributions. Although they required very different skills, they share the same goal: to reduce the computational cost of the simulation. On one hand, classical boundary conditions like PML (Perfectly Matched Layers) are unstable in TTI media. We have proposed a formulation of a stable ABC (Absorbing Boundary Condition) in anisotropic media. The technique is based on slowness curve properties, giving to our approach an original side. On the other hand, the initial parallelism, which is based on a domain decomposition and communications by message passing through the MPI library, leads to load-imbalance and so poor parallel efficiency. We have fixed this issue by replacing the paradigm for parallelism by the use of task-based programming through runtime system. This PhD thesis have been done in the framework of the research action DIP (Depth Imaging Partnership) between the Total oil company and Inria.

Les instabilités paramétriques non linéaires (NL) ont été observées sur les ondes magnéto-élastiques dans le cas d’un couplage de trois quasi-phonons sous pompage électromagnétique. La théorie en prédit une dynamique supercritique explosive, mais limitée expérimentalement par le décalage de fréquence dû aux fortes nonlinéarités. La dynamique supercritique des instabilités paramétriques NL est étudiée dans deux matériaux antiferromagnétiques "plan facile" (AFEP): l’hématite α-Fe2O3 et le borate de fer FeBO3. Ces matériaux possèdent une très grande NL acoustique effective en raison du couplage magnéto-élastique élevé. Les mécanismes de limitation de la dynamique explosive ont été analysés à l'aide de l'approximation anharmonique. La compensation du décalage fréquentiel NL par une modulation de phase singulière du pompage a été proposée et théoriquement vérifiée, puis utilisée pour l’observation expérimentale de la dynamique supercritique explosive des excitations de trois quasi-phonons dans les résonateurs magnéto-élastiques. Les études sur FeBO3 ont été réalisées dans la gamme de température 77 K - 293 K où les paramètres magnéto-élastiques du cristal varient de façon significative. Un modèle fortement non linéaire des excitations de trois quasi-phonons dans les AFEPs a été développé. Les simulations numériques sont en accord avec les résultats expérimentaux. Les études théoriques de couplage de trois ondes magnéto-élastiques progressives ont été effectuées sur la base de modèles théoriques prenant en compte la non-linéarité cubique des cristaux AFEP réels. Les simulations numériques prévoient un comportement explosif et une localisation spatiale des triades générées / Recently discovered nonlinear parametric instabilities occur when nonlinear parameter of a system is modulated. These instabilities were reported on magnetoelastic waves as three quasi-phonon coupling under electromagnetic pumping. Theoretical studies predicted supercritical explosive dynamics of these instabilities. Experimentally such singular behavior is limited by strong nonlinear frequency shift.Presented work studies supercritical dynamics of nonlinear parametric instabilities in two easy plane antiferromagnets (AFEP): hematite α-Fe2O3 and iron borate FeBO3. These materials possess unprecedented effective acoustic nonlinearity due to high magneto-elastic coupling. Limiting mechanisms of explosive dynamics were analyzed with the help of anharmonic approximation. Nonlinear frequency shift compensation via singular pumping field phase modulation was suggested and theoretically approbated. This technique was successfully used for experimental observation and investigation of supercritical explosive dynamics of three quasi-phonon excitations in magnetoelastic resonators. Iron borate studies were performed in the temperature range 77 K – 293 K where magnetoelastic parameters of the crystal vary essentially. Strongly nonlinear model of three quasi-phonon excitations in AFEPs was developed. Numerical simulations of the model showed good agreement with experimental results.Theoretical studies of three travelling magnetoelastic waves coupling are performed on the basis of suggested theoretical models that take into account cubic nonlinearity of real AFEP crystals. Numerical simulations of the models suggest explosive behavior and spatial localization of generated triads

La connaissance des structures internes de la Terre, à différentes échelles, présente des enjeux majeurs d'ordres économiques, humains, environnementaux et scientifiques. Diverses méthodes d'imagerie ont été développées en utilisant les informations contenues dans les ondes sismiques. La méthode d'inversion des formes d'ondes construit des images quantitative haute résolution des paramètres physiques du sous-sol, en exploitant le champ d'onde complet, sous la forme d'un problème d'optimisation. Dans ce travail de thèse, je présente l'application de l'inversion des formes d'ondes en domaine fréquentiel, pour imager les paramètres visco-élastiques dans des géometries à deux dimensions à grands offsets. Dans un premier temps les développements méthodologiques et algorithmiques sont présentés. La modélisation de la propagation des ondes P-SV en domaine fréquentiel, le problème direct du processus d'imagerie, est assurée par une méthode d'éléments finis Galerkin discontinus, assurant une grande flexibilité dans le choix des ordres d'interpolation et dans l'utilisation de maillages triangulaires non-structurés. Le problème inverse est résolu sous une forme linéarisée, afin de limiter le nombre de simulations directes, et utilise l'algorithme quasi-Newton L-BFGS permettant de tirer bénéfice de l'estimation "économique" du Hessien. Le processus global d'imagerie est implémenté sous la forme d'un algorithme massivement parallèle destiné aux calculateurs modernes à mémoire distribuée. Dans un deuxième temps, les algorithmes développés sont appliqués à des cas d'étude. Des applications sont menées dans des modèles synthétiques réalistes représentatifs d'environnements terrestres et marins. Ces études montrent les difficultés associées à la reconstruction des paramètres élastiques à partir de données mettant en jeu des phénomènes de propagations complexes (ondes converties, multiples, ondes de surfaces...). Des solutions sont proposées sous forme de processus hiérarchiques multi-échelles, afin de limiter les effets des non-linéarités du problème inverse et ainsi d'améliorer la convergence du processus vers le minimum global. Enfin, la sensibilité de différentes normes et critères de minimisation est analysée, à partir de données bruités issues de modèles synthétiques réalistes, ainsi que sous l'approximation acoustique pour un jeu de données réelles pétrolière. Ces tests montrent certaines limites du formalisme classique basé sur la norme L2 dans l'espace des données, tandis que la norme L1 apparaît comme alternative robuste pour l'inversion de données décimées en domaine fréquentiel.

Une modélisation numérique des ondes poroélastiques, décrites par le modèle de Biot, est proposée dans le domaine temporel. La dissipation visqueuse à l'intérieur des pores est décrite par le modèle de perméabilité dynamique, développé par Johnson-Koplik-Dashen (JKD). Certains coefficients du modèle de Biot-JKD sont proportionnels à la racine carrée de la fréquence : dans le domaine temporel, ces coefficients introduisent des dérivées fractionnaires décalées d'ordre 1/2, qui reviennent à un produit de convolution. Basé sur une représentation diffusive, le produit de convolution est remplacé par un nombre fini de variables de mémoire, dont la relaxation est gouvernée par une équation différentielle ordinaire locale en temps, ce qui mène au modèle de Biot-DA (approximation diffusive). Les propriétés du modèle de Biot-JKD et du modèle de Biot-DA sont analysées : hyperbolicité, décroissance de l'énergie, dispersion. Pour déterminer les coefficients de l'approximation diffusive, différentes méthodes de quadrature sont proposées : quadratures de Gauss, procédures d'optimisation linéaire ou non-linéaire sur la plage de fréquence d'intérêt. On montre que l'optimisation non-linéaire est la meilleure méthode de détermination. Le système est modélisé numériquement en utilisant une méthode de splitting : la partie propagative est discrétisée par un schéma aux différences finies ADER, d'ordre 4 en espace et en temps, et la partie diffusive est intégrée exactement. Une méthode d'interface immergée est implémentée pour discrétiser la géometrie sur une grille cartésienne et pour discrétiser les conditions de saut aux interfaces. Des simulations numériques sont présentées, pour des milieux isotropes et isotropes transverses. Des comparaisons avec des solutions analytiques montrent l'efficacité et la précision de cette approche. Des simulations numériques en milieux complexes sont réalisées : influence de la porosité d'os spongieux, diffusion multiple en milieu aléatoire.

Un système microfluidique à onde élastique de surface a été développé dans le but de réaliser une réaction d'amplification de brins d'ADN par PCR. Nous avons principalement étudié la température et l'uniformité de l'échauffement des gouttes irradiées par des ondes de type Rayleigh. Ces dernières sont générées à la surface d'un substrat piézoélectrique de Niobate de Lithium (LiNbO3). Nous avons pensé un système consommant le moins d'énergie électrique possible pour atteindre les températures désirées et permettant une meilleure uniformité de la température des gouttes. Pour cela, un dispositif à transducteur enterré a été réalisé sous une couche isolante de silice. Les gouttes sont ainsi directement posées sur le transducteur ce qui limite les pertes et améliore la répartition de la chaleur au sein des gouttes. Nous avons ensuite vérifié que les réactifs de la PCR ne sont pas affectés par les ondes de Rayleigh ce qui laisse présager que la PCR peut être réalisée à l'aide d'un système d'échauffement par ondes de Rayleigh. Par ailleurs, le déplacement de goutte sur ce type de substrat de LiNbO3 est important pour des applications de type laboratoire sur puce. Ce substrat ayant des propriétés hydrophiles, des revêtements ont été développés afin de minimiser la force nécessaire à l'actionnement des gouttes. Dans ces travaux, nous proposons un nouveau type de couche à base de copolymère P(VDF-TrFe) dont la fabrication est simplement réalisée par dissolution et étalement par spin-coating. Nous avons montré que ce type de couche n'affecte que très peu la propagation des ondes de Rayleigh et les propriétés hydrophobes sont équivalentes à d'autres revêtements / In this work, a microfluidic system based on surface acoustic wave has been developed in order to achieve the amplification of DNA strands by temperature cycling (PCR). We studied mainly the temperature and the heat uniformity of microdroplets irradiated by Rayleigh waves. These waves are generated at the surface of a lithium niobate substrate. We propose a system allowing better temperature uniformity within microdroplets with an optimal energy consumption. For this, a device with buried transducer has been developed under an insulating layer (Silice). The droplets are then placed directly on the transducer which limits losses and improves the distribution of heat within the microdroplets. We then verified that the PCR reagents are not affected by the Rayleigh waves which suggests that PCR can be performed using a heating system by Rayleigh waves. Moreover, the move of microdroplets on this kind of LiNbO3 substrate is important for lab on chip applications. This substrate having hydrophilic properties, some coatings have been developed to minimize the required force to actuate the microdroplets. In this work, we developed a new hydrophobic layer based on copolymer P(VDF-TrFe) whose production is simply made by dissolving and spreading by spin-coating. We have shown that this kind of layer is compatible with Rayleigh waves and that the hydrophobic properties are equivalent to other coatings

Les cristaux granulaires sont des structures périodiques de particules disposées en réseau cristallin. Les interactions entre ces billes peuvent être modélisées par leurs contacts, qui ont des dimensions et des masses effectives beaucoup plus petites que celles des billes. Ceci induit une propagation d'ondes élastiques dans les structures granulaires avec des vitesses significativement plus lentes que dans le matériau des grains individuels. En outre, en raison de forces de cisaillement non centrales, les rotations de particules peuvent être initiées, conduisant à des modes de phononiques supplémentaires dans ces cristaux. Dans ce manuscrit, on étudie la propagation d’ondes dans les cristaux granulaires monocouche bidimensionnels avec un mouvement des particules hors-plan ou dans le plan. Les propriétés phononiques sont étudiées, y compris les points de Dirac, les modes de fréquence nulle, les modes à vitesse de groupe nulle et leur transformation en modes de propagation lente. En outre, en présence de bords, on peut prévoir également des ondes de bord élastiques à fréquence nulle et extrêmement lentes dans des cristaux granulaires en « nid d'abeille » (graphène granulaire). En outre, les propriétés topologiques des ondes de bord rotationelles-transverses dans un graphène granulaire sont théoriquement démontrées. En induisant une transition topologique, qui transforme l'ordre topologique du graphène granulaire de trivial en non trivial, on peut observer le transport de bord topologique dans le graphène granulaire. Les théories développées pourraient mener potentiellement à des applications sur le contrôle des ondes élastiques par des structures granulaires. / Granular crystals are spatially periodic structures of elastic particles arranged in crystal lattices. The interactions between particles take place via their elastic interconnections, which are of much smaller dimensions and weights than the beads. This induces propagation of elastic waves in granular structures at significantly slower velocities than in the individual grains. In addition, due to the existence of non-central shear forces, rotations of particles can be initiated, leading to extra phononic modes in the crystals. In the manuscript, wave dynamics in two-dimensional monolayer granular crystals with either out-of-plane or in-plane particle motion is studied. The phononic properties are investigated, including Dirac points, zero-frequency modes, zero-group-velocity modes and their transformation into slow propagating phononic modes. Furthermore, in the presence of edges/boundaries, zero-frequency and extremely slow elastic edge waves can be also predicted in mechanical granular honeycomb crystals (granular graphene). In addition, topological properties of rotational edge waves in a granular graphene are theoretically demonstrated. By inducing topological transition, which turns the topological order of granular graphene from trivial to nontrivial, topological edge transport in the granular graphene can be observed. The developed theories could promote the potential applications of designed granular structures with novel elastic wave propagation properties.