Wave propagation in two-dimensional periodic lattices

Plane wave propagation in infinite two-dimensional periodic lattices is investigated using Floquet-Bloch principles. Frequency bandgaps and spatial filtering phenomena are examined in four representative planar lattice topologies: hexagonal honeycomb, Kagomé lattice, triangular honeycomb, and the square honeycomb. These topologies exhibit dramatic differences in their long-wavelength deformation properties. Long-wavelength asymptotes to the dispersion curves based on homogenization theory are in good agreement with the numerical results for each of the four lattices. The slenderness ratio of the constituent beams of the lattice (or relative density) has a significant influence on the band structure. The techniques developed in this work can be used to design lattices with a desired band structure. The observed spatial filtering effects due to anisotropy at high frequencies (short wavelengths) of wave propagation are consistent with the lattice symmetries.     

Exact wave-based Bloch analysis procedure for investigating wave propagation in two-dimensional periodic lattices

This paper presents an exact, wave-based approach for determining Bloch waves in two-dimensional periodic lattices. This is in contrast to existing methods which employ approximate approaches (e.g., finite difference, Ritz, finite element, or plane wave expansion methods) to compute Bloch waves in general two-dimensional lattices. The analysis combines the recently introduced wave-based vibration analysis technique with specialized Bloch boundary conditions developed herein. Timoshenko beams with axial extension are used in modeling the lattice members. The Bloch boundary conditions incorporate a propagation constant capturing Bloch wave propagation in a single direction, but applied to all wave directions propagating in the lattice members. This results in a unique and properly posed Bloch analysis. Results are generated for the simple problem of a periodic bi-material beam, and then for the more complex examples of square, diamond, and hexagonal honeycomb lattices. The bi-material beam clearly introduces the concepts, but also allows the Bloch wave mode to be explored using insight from the technique. The square, diamond, and hexagonal honeycomb lattices illustrate application of the developed technique to two-dimensional periodic lattices, and allow comparison to a finite element approach. Differences are noted in the predicted dispersion curves, and therefore band gaps, which are attributed to the exact procedure more-faithfully modeling the finite nature of lattice connection points. The exact method also differs from approximate methods in that the same number of solution degrees of freedom is needed to resolve low frequency, and arbitrarily high frequency, dispersion branches. These advantageous features may make the method attractive to researchers studying dispersion characteristics, band gap behavior, and energy propagation in two-dimensional periodic lattices.     

Lamb wave propagation in magnetoelectroelastic plates

Based on the three-dimensional linear elastic equations and magnetoelectroelastic constitutive relations, propagation of symmetric and antisymmetric Lamb waves in an infinite magnetoelectroelastic plate is investigated. The coupled differential equations of motion are solved, and the phase velocity equations of symmetric and antisymmetric modes are obtained for both electrically and magnetically open and shorted cases. The dispersive characteristic of wave propogation is explored. The mechanical, electric and magnetic responses of the lowest symmetric and antisymmetric Lamb wave modes are discussed in detailed. Obtained results are valuable for the analysis and design of broadband magnetoelectric transducer using composite materials.     

Elastic guided wave propagation in a periodic array of multi-layered piezoelectric plates with finite cross-sections

In this study, we present a model study of guided wave dispersion and resonance behavior of an array of piezoelectric plates with arbitrary cross-sections. The objective of this work is to analyze the influence of the geometry of an element of a 1D-array ultrasound transducer on generating multi-resonance frequency so as to increase the frequency bandwidth of the transducers. A semi-analytical finite-element (SAFE) method is used to model guided wave propagation in multi-layered 1D-array ultrasound transducers. Each element of the array is composed of LiNbO₃ piezoelectric material with rectangular or subdiced cross-section. Four-node bilinear finite-elements have been used to discretize the cross-section of the transducer. Dispersion curves showing the dependence of phase and group velocities on the frequency, and mode shapes of propagating modes were obtained for different geometry consurations. A parametric analysis was carried out to determine the effect of the aspect ratio, subdicing, inversion layer and matching layers on the vibrational behavior of 1D-array ultrasound transducers. It was found that the geometry with subdiced cross-section causes more vibration modes compared with the rectangular section. Modal analysis showed that the additional modes correspond to lateral modes of the piezoelectric subdiced section. In addition, some modes have strong normal displacements, which may influence the bandwidth and the pressure field in front of the transducer. In addition, the dispersion curves reveal strong coupling between waveguide modes due to the anisotropy of the piezoelectric crystal. The effect of the matching layers was to cluster extensional and flexural modes within a certain frequency range. Finally, inversion layer is found to have a minor effect on the dispersion curves. This analysis may provide a means to analyze and understand the dynamic response of 1D-array ultrasound transducers.     

Harmonic wave propagation in materials with periodic beam-structure

The characteristics of harmonic waves propagating in periodic beam structures are investigated. For this purpose, a very effective method for the calculation of dispersion relations of elastic waves in these materials is developed by applying dynamic theory of crystal lattices to discrete models of periodic beam structures. This method is applicable to general three-dimensional periodic beam structures. Results presented show that the solutions converge to the exact solution as the number of atoms in the discrete models increases. The dispersion relations of plane harmonic waves in the micropolar continuum model, developed in a previous study, are also calculated. They are compared with the exact solutions to examine the applicability of the continuum models to dynamic problems.     

High-directional wave propagation in periodic loss modulated materials

Amplification/attenuation of light waves in artificial materials can become sensitive to the propagation direction by spatially modulating the gain/loss response of the medium on the wavelength scale. We give a numerical proof of the high anisotropy of the gain/loss in two dimensional periodic structures with square and rhombic lattice symmetry by solving the full set of Maxwell's equations using the finite difference time domain method. Anisotropy of amplification/attenuation leads to the narrowing of the angular spectrum of propagating radiation with wavevectors close to the edges of the first Brillouin Zone. The effect provides a novel and useful method to filter out high spatial harmonics from noisy beams.     

Plane wave propagation and frequency band gaps in periodic plates and cylindrical shells with and without heavy fluid loading

Free plane wave propagation in infinitely long periodic elastic structures with and without heavy fluid loading is considered. The structures comprise continuous elements of two different types connected in an alternating sequence. In the absence of fluid loading, an exact solution which describes wave motion in each unboundedly extended element is obtained analytically as a superposition of all propagating and evanescent waves, continuity conditions at the interfaces between elements are formulated and standard Floquet theory is applied to set up a characteristic determinant. An efficient algorithm to compute Bloch parameters (propagation constants) as a function of the excitation frequency is suggested and the location of band gaps is studied as a function of non-dimensional parameters of the structure's composition. In the case of heavy fluid loading, an infinitely large number of propagating or evanescent waves exist in each unboundedly extended elasto-acoustic element of a periodic structure. Wave motion in each element is then presented in the form of a modal decomposition with a finite number of terms retained in these expansions and the accuracy of such an approximation is assessed. A generalized algorithm is used to compute Bloch parameters for a periodic structure with heavy fluid loading as a function of the excitation frequency and, similarly to the previous case, the location of band gaps is studied.     

Theory of electromagnetic wave propagation in a longitudinally periodic cylinder

Thinking of photonic crystals, we investigate the theory of electromagnetic wave propagation in a perfectly conducting semi-cylinder endowed with a periodic permittivity along its axis while its circular base is illuminated by an harmonic Bessel beam, symmetric around the cylinder axis. We prove that the Floquet-Bessel expansions of electromagnetic fields whose Floquet parts are solutions of a Mathieu equation. is a suitable tool to handle this kind of problems.     

Shear-horizontal acoustic wave propagation in piezoelectric bounded plates with metal gratings

In this paper, shear-horizontal (SH) acoustic wave propagation in metal gratings deposited on piezoelectric bounded plates is investigated. The spectral characteristics of the electromechanical coupling coefficient are studied first, which are very important for acoustic wave device designs. And, an effective mathematic method based on even- and odd base functions is also presented for overcoming the large frequency thickness product problem. Then, the characteristics of the grating modes are studied, and the nature and characteristics of the stop bands are investigated fully. The results show that the width and attenuation of the stop bands are dominated by the electromechanical coupling coefficient at the frequency centers of the stop bands.     

基于高级统计能量分析的周期加筋板振动特性研究

统计能量分析(statistical energy analysis, SEA)是复杂耦合系统中、高频动力学特性计算的有力工具. 本文以波传播理论和SEA的基本原理为基础, 研究周期加筋板中弯曲波传播特性. 分析了周期结构的频率带隙特性和加强筋对板上弯曲波的滤波特性对SEA计算结果的影响规律, 发现经典SEA由于忽视了加筋板中物理上不相邻子系统间存在的能量隧穿效应, 而导致响应预测结果产生最高近 40 dB的误差. 为了解决这一问题, 本文应用高级统计能量分析(advanced statistical energy analysis, ASEA)方法, 考虑能量在不相邻子系统间的传递、转移和转化的物理过程, 从而大幅提高子系统响应的预测精度, 将误差在大部分频段降低至小于5 dB. 设计了模拟简支边界条件的加筋板振动测试实验装置, 实验测试结果与有限元结果符合较好, 对理论模型进行了验证.     

Free vibration analysis of planar curved beams by wave propagation

In this paper, a systematic approach for the free vibration analysis of a planar circular curved beam system is presented. The system considered includes multiple point discontinuities such as elastic supports, attached masses, and curvature changes. Neglecting transverse shear and rotary inertia, harmonic wave solutions are found for both extensional and inextensional curved beam models. Dispersion equations are obtained and cut-off frequencies are determined. Wave reflection and transmission matrices are formulated, accounting for general support conditions. These matrices are combined, with the aid of field transfer matrices, to provide a concise and efficient method for the free vibration problem of multi-span planar circular curved beams with general boundary conditions and supports. The solutions are exact since the effects of attenuating wave components are included in the formulation. Several examples are presented and compared with other methods.     

Micro-crack detection of nonlinear Lamb wave propagation in three-dimensional plates with mixed-frequency excitation

We propose a nonlinear ultrasonic technique by using the mixed-frequency signals excited Lamb waves to conduct micro-crack detection in thin plate structures. Simulation models of three-dimensional(3D) aluminum plates and composite laminates are established by ABAQUS software, where the aluminum plate contains buried crack and composite laminates comprises cohesive element whose thickness is zero to simulate delamination damage. The interactions between the S_0 mode Lamb wave and the buried micro-cracks of various dimensions are simulated by using the finite element method.Fourier frequency spectrum analysis is applied to the received time domain signal and fundamental frequency amplitudes,and sum and difference frequencies are extracted and simulated. Simulation results indicate that nonlinear Lamb waves have different sensitivities to various crack sizes. There is a positive correlation among crack length, height, and sum and difference frequency amplitudes for an aluminum plate, with both amplitudes decreasing as crack thickness increased, i.e.,nonlinear effect weakens as the micro-crack becomes thicker. The amplitudes of sum and difference frequency are positively correlated with the length and width of the zero-thickness cohesive element in the composite laminates. Furthermore,amplitude ratio change is investigated and it can be used as an effective tool to detect inner defects in thin 3D plates.     

Reverberation-ray matrix analysis for wave propagation in multiferroic plates with imperfect interfacial bonding

The dispersion behavior of waves in multiferroic plates with imperfect interfacial bonding has been investigated via the method of reverberation-ray matrix, which is directly established from the three-dimensional equations of magneto-electro-elasticity in the form of state space formalism. A generalized spring-layer model is employed to characterize the interfacial imperfection. By introducing a dual system of local coordinates for each single layer, the numerical instability usually encountered in the state space method can be avoided. Based on the proposed method, a typical sandwich plate made of piezoelectric and piezomagnetic phases is considered in numerical examples to calculate the dispersion curves and mode shapes. It is demonstrated that the results obtained by the present method is unconditionally stable as compared to the traditional state space method. The influence of different interfacial bonding conditions on the dispersion characteristics and corresponding mode shapes is investigated.     

The effect of disorder on the wave propagation in one-dimensional periodic optical systems

The influence of disorder on the transmission through periodic waveguides is studied. Using a canonical form of the transfer matrix, we investigate the dependence of the Lyapunov exponent γ on the frequency ν and magnitude of the disorder σ. It is shown that in the bulk of the bands γ?～?σ², while near the band edges it has order γ?～?σ^2/3. This dependence is illustrated by numerical simulations.     

Note on wave propagation through a periodic rectangular wave-guide containing ferrite slabs

Wave propagation through a rectangular wave-guide loaded with periodic ferrite slabs has been considered. Numerical computations for the ratio of amplitudes and the phase shift has been made for the first three harmonics only.I take this opportunity of thanking Dr. S. C.Dasgupta and Dr. S. S.Baral for suggesting this investigation and for their guidance in preparing this note. I also wish to express my gratitude to the C.S.I.R. for financial support as a Senior Research Fellow.     

Analysis of wave propagation in a thin composite cylinder with periodic axial and ring stiffeners using periodic structure theory

Wave propagation characteristics of a thin composite cylinder stiffened by periodically spaced ring frames and axial stringers are investigated by an analytical method using periodic structure theory. It is used for calculating propagation constants in axial and circumferential directions of the cylindrical shell subject to a given circumferential mode or axial half-wave number. The propagation constants corresponding to several different circumferential modes and/or half-wave numbers are combined to determine the vibrational energy ratios between adjacent basic structural elements of the two-dimensional periodic structure. Vibration analyses to validate the theoretical development have been carried out on sufficiently detailed finite element model of the same dimension and configuration as the stiffened cylinder and very good agreement is obtained between the analytical and the dense finite element results. The effects of shell material properties and the length of each periodic element on the wave propagation characteristics are also examined based on the current analytical approach.