WAVELET METHOD FOR CALCULATING THE DEFECT STATES OF TWO-DIMENSIONAL PHONONIC CRYSTALS

Based on the variational theory, a wavelet-based numerical method is developed to calculate the defect states of acoustic waves in two-dimensional phononic crystals with point and line defects. The supercell technique is applied. By expanding the displacement field and the material constants （mass density and elastic stiffness） in periodic wavelets, the explicit formulations of an eigenvalue problem for the plane harmonic bulk waves in such a phononic structure are derived. The point and line defect states in solid-liquid and solid-solid systems are calculated. Comparisons of the present results with those measured experimentally or those from the plane wave expansion method show that the present method can yield accurate results with faster convergence and less computing time.     

Application of plane wave expansion and stiffness matrix methods to study transmission properties and guided mode of phononic plates

Numerical procedure based on plane wave expansion and stiffness matrix method is developed to calculate the transmission factor of a micro two-dimensional phononic plate. Calculations of the dispersion curve have been achieved by introducing particular functions which transform motion equations into an eigenvalue problem. The state vector has been generalized to a phononic material, it leads to a comparatively convenient matrix formulation. The influence of the layer number on the transmission factor is studied. In addition, our interest is focused on the observed gap and how it behaves when phononic structure undergoes a slight change. The result shows that if the central phononic layer is replaced by one or two homogeneous layers, guided modes originate inside the frequency band gaps.     

INFLUENCES OF ANISOTROPY ON BAND GAPS OF 2D PHONONIC CRYSTAL

Band gaps of 2D phononic crystal with orthotropic cylindrical fillers embedded in the isotropic host are studied in this paper. Two kinds of periodic structures, namely, the square lattice and the triangle lattice, are considered. For anisotropic phononic crystal, band gaps not only depend on the periodic lattice but also the angle between the symmetry axis of orthotropic material and that of the periodic structure. Rotating these cylindrical fillers makes the angle changing continuously; as a result, pass bands and forbidden bands of the phononic crystal are changed. The plane wave expansion method is used to reduce the band gap problem to an eigenvalue problem. The numerical example is given for YBCO/Epoxy composites. The location and the width of band gaps are estimated for different rotating angles. The influence of anisotropy on band gaps is discussed based on numerical results.     

含有功能梯度材料的一维声子晶体弹性波带隙研究

将功能梯度材料引入声子晶体之中,构成由均匀材料和功能梯度材料交替复合而成的材料常数宏观上呈连续变化的一维声子晶体结构.应用平面波展开法研究了该结构中存在的弹性波带隙特征,并得到前四阶归一化频率带隙结构随材料长度和功能梯度材料指数因子变化的关系图.通过对其带隙结构与不含功能梯度材料的常规声子晶体带隙结构比较可知:功能梯度材料使原有带隙结构发生变化,含有功能梯度材料的声子晶体较常规声子晶体前两阶带隙归一化起始频率升高.这些结果扩宽了功能梯度材料的使用范围,提升了声子晶体工程可靠性,为含有功能梯度材料的声子晶体在工程实际中的应用提供了理论依据和指导.     

Band gaps of elastic waves in three-dimensional piezoelectric phononic crystals with initial stress

In this paper, the stop band properties of elastic waves in three-dimensional piezoelectric phononic crystals with initial stress are studied taking the mechanical and electrical coupling into account. The band gap characteristics for three kinds of lattice arrangements (i.e. sc, bcc and fcc) are investigated by the plane wave expansion (PWE) method. Regarding the variables of mechanical and electrical fields as the elements of the generalized state vector, the expression of the generalized eigenvalue equation for three-dimensional piezoelectric periodic structures is derived. Numerical calculations are performed for the PZT-2/polymer and ZnO/polymer phononic crystals. It can be observed from the results that the fcc lattice is more favorable to create the stop band than the sc and bcc lattices for the piezoelectric phononic crystals, which has also been proved for the pure elastic periodic structures. Compared with the PZT-2/polymer systems, the band gap of the sc lattice for the ZnO/polymer structures is narrower. However, the widths of the bcc and fcc lattices for the ZnO/polymer phononic crystals are much larger than those for the PZT-2/polymer structures. The lattice arrangements and the piezoelectricity have remarkable influences on the stop band behaviors.     

Propagation of bending waves in phononic crystal thin plates with a point defect

In this paper, the band structures of bending waves in a phononic crystal thin plate with a point defect are studied using an improved plane wave expansion method combined with the supercell technique. In particular, a phononic crystal thin plate composed of an array of circular crystalline Al₂O₃ cylinders embedded in an epoxy matrix with a square lattice is considered in detail. Full band gaps are shown. When a point defect is introduced, the bending waves are highly localized at or near the defect, resulting in defect modes. The frequency and number of the defect modes are strongly dependent on the filling fraction of the system and the size of the point defect. The defect bands appear from the upper edge of the gap and move to the middle of the gap as the defect size is reduced. For a given defect, the frequencies of the defect bands increase as the filling fraction increases.     

Optimum design of phononic crystal perforated plate structures for widest bandgap of fundamental guided wave modes and maximized in-plane stiffness

This paper presents a topology optimization of single material phononic crystal plate (PhP) to be produced by perforation of a uniform background plate. The primary objective of this optimization study is to explore widest exclusive bandgaps of fundamental (first order) symmetric or asymmetric guided wave modes as well as widest complete bandgap of mixed wave modes (symmetric and asymmetric). However, in the case of single material porous phononic crystals the bandgap width essentially depends on the resultant structural integration introduced by achieved unitcell topology. Thinner connections of scattering segments (i.e. lower effective stiffness) generally lead to (i) wider bandgap due to enhanced interfacial reflections, and (ii) lower bandgap frequency range due to lower wave speed. In other words higher relative bandgap width (RBW) is produced by topology with lower effective stiffness. Hence in order to study the bandgap efficiency of PhP unitcell with respect to its structural worthiness, the in-plane stiffness is incorporated in optimization algorithm as an opposing objective to be maximized. Thick and relatively thin Polysilicon PhP unitcells with square symmetry are studied. Non-dominated sorting genetic algorithm NSGA-II is employed for this multi-objective optimization problem and modal band analysis of individual topologies is performed through finite element method. Specialized topology initiation, evaluation and filtering are applied to achieve refined feasible topologies without penalizing the randomness of genetic algorithm (GA) and diversity of search space. Selected Pareto topologies are presented and gradient of RBW and elastic properties in between the two Pareto front extremes are investigated. Chosen intermediate Pareto topology, even not extreme topology with widest bandgap, show superior bandgap efficiency compared with the results reported in other works on widest bandgap topology of asymmetric guided waves, available in the literature. Finally, steady state and transient frequency response of finite thin PhP structures of selected Pareto topologies are studied and validity of obtained bandgaps is confirmed.     

Electro-mechanical coupling wave propagating in a locally resonant piezoelectric/elastic phononic crystal nanobeam with surface effects

The model of a "spring-mass" resonator periodically attached to a piezoelectric/elastic phononic crystal(PC) nanobeam with surface effects is proposed, and the corresponding calculation method of the band structures is formulized and displayed by introducing the Euler beam theory and the surface piezoelectricity theory to the plane wave expansion(PWE) method. In order to reveal the unique wave propagation characteristics of such a model, the band structures of locally resonant(LR) elastic PC Eul...     

Simulation of elastic wave diffraction by multiple strip-like cracks in a layered periodic composite

The problem of numerical simulation of the steady-state harmonic vibrations of a layered phononic crystal (elastic periodic composite) with a set of strip-like cracks parallel to the layer boundaries is solved, and the accompanying wave phenomena are considered. The transfer matrix method (propagator matrix method) is used to describe the incident wave field. It allows one not only to construct the wave fields but also to calculate the pass bands and band gaps and to find the localization factor. The wave field scattered by multiple defects is represented by means of an integral approach as a superposition of the fields scattered by all cracks. An integral representation in the form of a convolution of the Fourier symbols of Green’s matrices for the corresponding layered structures and a Fourier transform of the crack opening displacement vector is constructed for each of the scattered fields. The crack opening displacements are determined by the boundary integral equation method using the Bubnov-Galerkin scheme, where Chebyshev polynomials of the second kind, which take into account the behavior of the solution near the crack edges, are chosen as the projection and basis systems. The system of linear algebraic equations with a diagonal predominance of components arising when the system of integral equations is discretized has a block structure. The characteristics describing qualitatively and quantitatively the wave processes that take place under the diffraction of plane elastic waves by multiple cracks in a phononic crystal are analyzed. The resonant properties of a system of defects and the influence of the relative positions and sizes of defects in a layered phononic crystal on the resonant properties are studied. To obtain clearer results and to explain them, the energy flux vector is calculated and the energy surfaces and streamlines corresponding to them are constructed.     

PHONONIC BAND GAPS IN TWO-DIMENSIONAL HYBRID TRIANGULAR LATTICE

<正>Absolute phononic band gaps can be substantially improved in two-dimensional lattices by using a symmetry reduction approach.In this paper,the propagation of elastic waves in a two-dimensional hybrid triangular lattice structure consisting of stainless steel cylinders in air is investigated theoretically.The band structure is calculated with the plane wave expansion (PWE)method.The hybrid triangular Bravais lattice is formed by two kinds of triangular lattices. Different from ordinary triangular lattices,the band gap opens at low frequency(between the first and the second bands)regime because of lifting the bands degeneracy at high symmetry points of the Brillouin zone.The location and width of the band gaps can be tuned by the position of the additional rods.     

Modeling Phononic Crystals via the Weighted Relaxed Micromorphic Model with Free and Gradient Micro-Inertia

In this paper the relaxed micromorphic continuum model with weighted free and gradient micro-inertia is used to describe the dynamical behavior of a real two-dimensional phononic crystal for a wide range of wavelengths. In particular, a periodic structure with specific micro-structural topology and mechanical properties, capable of opening a phononic band-gap, is chosen with the criterion of showing a low degree of anisotropy (the band-gap is almost independent of the direction of propagation of the traveling wave). A Bloch wave analysis is performed to obtain the dispersion curves and the corresponding vibrational modes of the periodic structure. A linear-elastic, isotropic, relaxed micromorphic model including both a free micro-inertia (related to free vibrations of the microstructures) and a gradient micro-inertia (related to the motions of the microstructure which are coupled to the macro-deformation of the unit cell) is introduced and particularized to the case of plane wave propagation. The parameters of the relaxed model, which are independent of frequency, are then calibrated on the dispersion curves of the phononic crystal showing an excellent agreement in terms of both dispersion curves and vibrational modes. Almost all the homogenized elastic parameters of the relaxed micromorphic model result to be determined. This opens the way to the design of morphologically complex meta-structures which make use of the chosen phononic material as the basic building block and which preserve its ability of “stopping” elastic wave propagation at the scale of the structure.     

ELASTIC WAVE LOCALIZATION IN TWO-DIMENSIONAL PHONONIC CRYSTALS WITH ONE-DIMENSIONAL QUASI-PERIODICITY AND RANDOM DISORDER

The band structures of both in-plane and anti-plane elastic waves propagating in two-dimensional ordered and disordered （in one direction） phononic crystals are studied in this paper. The localization of wave propagation due to random disorder is discussed by introducing the concept of the localization factor that is calculated by the plane-wave-based transfer-matrix method. By treating the quasi-periodicity as the deviation from the periodicity in a special way, two kinds of quasi phononic crystal that has quasi-periodicity （Fibonacci sequence） in one direction and translational symmetry in the other direction are considered and the band structures are characterized by using localization factors. The results show that the localization factor is an effective parameter in characterizing the band gaps of two-dimensional perfect, randomly disordered and quasi-periodic phononic crystals. Band structures of the phononic crystals can be tuned by different random disorder or changing quasi-periodic parameters. The quasi phononic crystals exhibit more band gaps with narrower width than the ordered and randomly disordered systems.     

固-固掺杂结构声子晶体中弹性波的缺陷模

为了研究弹性波在固-固掺杂结构声子晶体中缺陷模的特性,推导出弹性波斜入射固-固掺杂结构声子晶体的转移矩阵。利用它计算了固-固掺杂结构声子晶体中弹性波的透射系数,得到由于掺杂在禁带的中间出现了多条缺陷模透射峰带,缺陷模透射峰带随杂质厚度的变化呈周期性地出现,在同一缺陷模透射峰带上缺陷模的频率随杂质厚度增加近似呈线性减少。缺陷模的这些特征可以作为设计声子晶体滤波器的理论依据。     

On the quasistatic effective elastic moduli for elastic waves in three-dimensional phononic crystals

Effective elastic moduli for 3D solid–solid phononic crystals of arbitrary anisotropy and oblique lattice structure are formulated analytically using the plane-wave expansion (PWE) method and the recently proposed monodromy-matrix (MM) method. The latter approach employs Fourier series in two dimensions with direct numerical integration along the third direction. As a result, the MM method converges much quicker to the exact moduli in comparison with the PWE as the number of Fourier coefficients increases. The MM method yields a more explicit formula than previous results, enabling a closed-form upper bound on the effective Christoffel tensor. The MM approach significantly improves the efficiency and accuracy of evaluating effective wave speeds for high-contrast composites and for configurations of closely spaced inclusions, as demonstrated by three-dimensional examples.     

Optimal 2D auxetic micro-structures with band gap

A systematic investigation is presented that explores band gap properties of periodic micro-structures architected for maximum auxeticity. The design of two-dimensional auxetic cells is addressed using inverse homogenization. A non-convex optimization problem is formulated that is solved through mathematical programming. Different starting guesses are used to explore local minima when distributing material and void or two materials and void. The same numerical tool succeeds in capturing re-entrant, chiral and anti-chiral layouts with negative Poisson’s ratio, retrieving solutions originally found through other approaches as well as generating variations. A Floquet–Bloch approach is then applied to the achieved periodic cells to investigate possible band gaps characterizing the in-plane wave propagation. Directional and full band-pass filters are found in the case of micro-structures whose auxetic behavior comes from the arising of a rotational deformation of the periodic cell. Such kind of topologies could be exploited to design tunable wave guides and tunable phononic crystals, respectively.

     

基于人工神经网络的声子晶体逆向设计

声子晶体是一种人工周期性复合材料, 其带隙特性使其在减振、隔声、滤波和声学功能器件等领域具有潜在的应用价值. 如何准确操纵声波和机械波是声子晶体设计的主要挑战. 现有设计方法是基于对结构几何参数与材料参数的分析调整使其匹配特定的应用特性, 设计效率不高且无法达到最佳性能. 为此, 本文以一维层状声子晶体为例, 提出了一种基于Softmax逻辑回归和多任务学习的人工神经网络声子晶体逆向设计方法, 其中, Softmax逻辑回归实现分层结构各区域材料种类的选择, 通过多任务学习确定各区域材料的分布, 从而, 将声子晶体逆向设计问题转化为对单位胞元拓扑结构多组分材料的分类问题. 首先, 随机生成大量声子晶体拓扑结构样本; 然后, 采用有限元法进行并行计算得到所有样本的带隙分布; 接着, 通过神经网络建立带隙分布和拓扑结构之间的映射关系; 最后, 利用训练好的神经网络设计具有目标带隙特性的声子晶体, 即以目标带隙作为神经网络的输入, 网络将直接输出对应的声子晶体单元胞元拓扑结构. 算例表明本方法可根据应用需求快速高效地得到具有目标带隙的一维声子晶体. 该方法为声子晶体的逆向设计提供了一种新颖思路.      

Wave localization in randomly disordered layered three-component phononic crystals with thermal effects

In this paper, the propagation and localization of elastic waves in randomly disordered layered three-component phononic crystals with thermal effects are studied. The transfer matrix is obtained by applying the continuity conditions between three consecutive sub-cells. Based on the transfer matrix method and Bloch theory, the expressions of the localization factor and dispersion relation are presented. The relation between the localization factors and dispersion curves is discussed. Numerical simulations are performed to investigate the influences of the incident angle on band structures of ordered phononic crystals. For the randomly disordered ones, disorders of structural thickness ratios and Lamé constants are considered. The incident angles, disorder degrees, thickness ratios, Lamé constants and temperature changes have prominent effects on wave localization phenomena in three-component systems. Furthermore, it can be observed that stopbands locate in very low-frequency regions. The localization factor is an effective way to investigate randomly disordered phononic crystals in which the band structure cannot be described.     

WAVE PROPAGATION IN TWO-DIMENSIONAL DISORDERED PIEZOELECTRIC PHONONIC CRYSTALS

The wave propagation is studied in two-dimensional disordered piezoelectric phononic crystals using the finite-difference time-domain （FDTD） method. For different cases of disorder, the transmission coefficients are calculated. The influences of disorders on band gaps are investigated. The results show that the disorder in the piezoelectric phononic crystals has more significant influences on the band gap in the low frequency regions than in the high frequency ones. The relation between the width of band gap and the direction of position disorder is also discussed. When the position disorder is along the direction perpendicular to the wave transmission, the piezoelectric phononic crystals have wider band gaps at low frequency regions than the case of position disorder being along the wave transmission direction. It can also be found that the effect of. size disorder on band gaps is analogous to that of location disorder. When the perturbation coefficient is big, it has more pronounced effects on the pass bands in the piezoelectric phononic crystals with both size and location disorders than in the piezoelectric phononic crystals with single disorder. In higher frequency regions the piezoelectric effect reduces the transmission coefficients. But for larger disorder degree, the effects of the piezoelectricity will be reduced.     

方形固-液声子晶体的全反射隧穿及其多通道滤波

为了研究一维固-液结构方形声子晶体中弹性波的全反射隧穿效应,利用传输矩阵法计算了弹性波在大于全反射角入射一维固-液结构方形声子晶体的透射率。在透射波中发现了全反射隧穿效应,并且全反射隧穿峰具有优良的多通道滤波特性。得出了一维固-液结构声子晶体的全反射隧穿的滤波特性随周期数、模式量子数、边长的变化规律。为设计性能优良的多通道声子晶体滤波器提供了理论依据。     

直角平面内弹性圆夹杂对入射平面SH波的散射

利用复变函数法、多极坐标移动技术及傅立叶级数展开求解二维直角平面内圆形弹性夹杂对稳态入射平面SH波的散射问题。首先写出直角平面内不含夹杂时的入射波场和反射波场;其次建立直角平面内含夹杂时夹杂外的散射波解和夹杂内的驻波解,并利用叠加原理写出问题的总波场,借助夹杂边界处应力和位移的连续条件建立求解散射波解和驻波解中未知系数的无穷代数方程组并求解,通过算例具体讨论了直角平面水平边界点的位移幅度比和夹杂边界处径向应力集中系数随不同无量纲波数、入射角及圆孔位置的变化情况,结果表明了算法的有效实用性。