固-固有限周期声子晶体中SH波全反射隧穿的机理

利用边界条件推导出SH波在一维固-固有限周期声子晶体中的转移矩阵,并用倏逝波的概念和多波束干涉理论,从理论上解释了一维固-固有限周期声子晶体全反射隧穿效应产生的物理机理,得出了全反射隧穿峰带波长满足的解析公式。利用解析公式和转移矩阵研究了一维固-固有限周期声子晶体中SH波的全反射隧穿效应,结果发现两种方法得出的结论是吻合的。     

固-固有限周期声子晶体中SH波全反射隧穿的机理

利用边界条件推导出SH波在一维固-固有限周期声子晶体中的转移矩阵,并用倏逝波的概念和多波束干涉理论,从理论上解释了一维固-固有限周期声子晶体全反射隧穿效应产生的物理机理,得出了全反射隧穿峰带波长满足的解析公式。利用解析公式和转移矩阵研究了一维固-固有限周期声子晶体中SH波的全反射隧穿效应,结果发现两种方法得出的结论是吻合的。     

半无限周期声子晶体的全反射隧穿能带

利用色散法研究了在入射角大于全反射角时弹性波在一维半无限周期声子晶体中的传输特性。结果发现当弹性波大于全反射角入射一维半无限周期声子晶体时出现了两级全反射隧穿导带;并得出了全反射隧穿导带随入射角和周期厚度的响应曲线,分析可知:导带的频率中心随入射角的增加而增加,随周期厚度的增加而降低;导带的频率宽度随入射角的增加而减小,随周期厚度的增加而降低。     

新型声子晶体全反射贯穿效应的滤波特性

利用转移矩阵法计算了弹性波在大于全反射角入射一维声子晶体的透射率,结果发现：在入射角岛为0．33rad、0．4rad、0．9rad、1．3rad附近都出现了明显的透射峰带,即在透射波中出现了全反射贯穿效应;其中在0．4rad附近的这个透射峰带的频率范围小、峰值细密,可以利用它来实现多通道滤波。同时研究了贯穿效应随介质厚度和周期数的变化特征,发现贯穿效应具有优良的滤波特性。利用这些特性可以设计出结构简单、控制方便的新型一维声子晶体的滤波器。     

固-液圆柱掺杂声子晶体中弹性波的缺陷模

利用一维固-液圆柱形掺杂声子晶体中弹性波横向受限的条件,推导出弹性波在一维固-液圆柱形掺杂声子晶体中各种模式满足的关系式,研究了各种模式弹性波的特性。利用转移矩阵计算出弹性波的缺陷模随模式量子数以及圆柱半径的变化规律,得到了一些一维固-液结构圆柱形掺杂声子晶体缺陷模的新特征。研究结果表明,缺陷模的频率和半高宽均与模式量子数和圆柱半径有关。     

镜像对称准周期固/液声子晶体的传输特性

以准周期Fibonacci序列为基础,构建了准周期镜像对称结构的一维固液介质声子晶体。利用转移矩阵法,研究了弹性波在所构建的声子晶体中的透射特性。结果表明:弹性波正入射时,转移矩阵与变介质的特性阻抗有关,由于声子晶体结构的局域失谐,在透射谱的中心频率附近对称出现谐振模;对于Q7序列的声子晶体,带隙内谐振模带宽小于10Hz,能实现超窄带声滤波;弹性波斜入射时,带隙向高频区移动;入射角为5°时,谐振模也向高频移动;入射角为10°时,由于谐振的增强,声子晶体结构对频率的选择性提高,使带隙内的谐振模消失。     

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

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

利用叠加理论研究固-液声子晶体纵波的能带结构

采用相干叠加原理推导出一种研究一维固-液有限周期声子晶体中纵波的透射率公式,建立了一种研究一维固-液有限周期声子晶体能带的新方法——波叠加理论。将波叠加理论和转移矩阵法进行了比较研究,结果表明波叠加理论和转移矩阵法得出的结果是一致的。波叠加理论不仅具有转移矩阵法的优点,而且又克服了转移矩阵法的不足。因此波叠加理论是一种研究一维固-液有限周期声子晶体中纵波能带的更为有效的方法。     

一维声子晶体带隙的非局部辛分析

周期性弹性复合结构(声子晶体)中传播的弹性波存在特殊的色散关系:弹性波只能在某段频率范围内无损耗的传播,该频率范围称为通带.一维声子晶体的色散问题可以看作分层介质中弹性波的传播问题,利用二维弹性理论予以分析.为了研究非局部效应对声子晶体带隙特性的影响,将Eringen的二维非局部弹性理论引入到Hamilton体系下,利用精细积分与扩展的Wittrick Williams算法可获取任意频率范围内的本征解.通过对不同算例的数值计算,分析和对比了非局部理论方法与传统局部理论方法的差别.并进一步指出了该套算法的适用性和优势所在.     

一维固液介质准周期结构声子晶体的传输特性

针对一维固液介质的Fibonacci序列准周期声子晶体,利用传输矩阵法研究了弹性纵波在准周期结构中的透射特性。计算结果表明：准周期结构具有比周期结构更宽的禁带,并伴随带边谐振;序列越高,带边谐振越强,并在较低序列的禁带中心出现谐振模;弹性波以0.1rad斜入射时,透射谱向高频方向移动;当固液介质波阻抗接近时,出现等间距的谐振模,同时禁带消失。     

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.     

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

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

Antiplane elastic wave propagation in pre-stressed periodic structures; tuning,band gap switching and invariance

The effect of nonlinear elastic pre-stress on antiplane elastic wave propagation in a two-dimensional periodic structure is investigated. The medium consists of cylindrical annuli embedded on a periodic square lattice in a uniform host material. An identical inhomogeneous deformation is imposed in each annulus and the theory of small-on-large is used to find the incremental wave equation governing subsequent small-amplitude antiplane waves. The plane-wave-expansion method is employed in order to determine the permissable eigenfrequencies. It is found that pre-stress significantly affects the band gap structure for Mooney–Rivlin and Fung type materials, allowing stop bands to be switched on and off. However, it is also shown that for a specific class of materials, their phononic properties remain invariant under nonlinear deformation, permitting some rather interesting behaviour and leading to the possibility of phononic cloaks.     

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.     

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.     

Absorbing boundary conditions for elastic waves in transversely isotropic media

In the numerical simulation of elastic wave propagation in the solid, it is essential to introduce absorbing boundary conditions to limit the large or unbounded domain of computation. In this paper, the absorbing boundaries for transversely isotropic media are composed of simple first-order partial differential operators, and each of the operators can perfectly absorb a plane wave outgoing at a certain angle. To test the absorbing ability, the reflection coefficient formulas for the quasi-P and quasi-S wave on the absorbing boundary are derived based on the potential functions theory of the elastic wave. Numerical examples show that the absorbing effect is good. The boundary conditions given here have a practical meaning.Supported by National Natural Science Foundation of China.