SYMPLECTIC ANALYSIS FOR WAVE PROPAGATION OF HIERARCHICAL HONEYCOMB STRUCTURES

Wave propagation in two-dimensional hierarchical honeycomb structures with twoorder hierarchy is investigated by using the symplectic algorithm. By applying the variational principle to the dual variables, the wave propagation problem is transformed into a two-dimensional symplectic eigenvalue problem. The band gaps and spatial filtering phenomena are examined to find the stop bands and directional stop bands. Special attention is directed to the effects of the relative density and the length ratio on the band gaps and phase constant surfaces. This work provides new opportunities for designing hierarchical honeycomb structures in sound insulation applications.     

Tunable low frequency band gaps and sound transmission loss of a lever-type metamaterial plate

A novel metamaterial plate with subwavelength lever-type resonators is proposed to obtain low frequency broadband band gaps and good sound insulation performance. The band structure is theoretically derived, and the validity of the theoretical method is verified by the finite element method. The formation mechanisms of the band gaps are illustrated by the analysis of the effective dynamic mass density and group velocity. The effect of the lever ratio on the band gaps is analyzed. The results indicate that as the lever ratio increases, the first band gap shifts to lower frequencies, while the bandwidth is widened. Moreover, the sound insulation performance of the proposed metamaterial plate is evaluated via examining the sound transmission loss (STL). Compared with the metamaterial plates without lever accessories, the proposed metamaterial plates with a suitable lever ratio have better sound insulation performance at low frequencies.     

蜂窝材料的弹性波传播特性

通过研究蜂窝材料的弹性波频散关系,分析了其弹性波传播特性. 采用基于小波理论的分析方法,将材料参数和位移均展开为双正交周期小波基函数的形式,利用Bloch定理将波动方程转化为特征值方程,求解该方程得到3种典型结构------正方、三角与六角排列的铝(Al)和聚丙烯(PP)蜂窝材料的频散关系,并进行了比较分析. 结果显示:结构形式的不同显著地影响其波动特性,而制作材料的不同则没有影响;3种结构形式都不存在完全带隙,但正方和三角形结构在一定的传播方向范围内存在方向带隙,而六角形结构则在任何方向都不存在方向带隙;与正方结构相比,三角结构在相同孔隙率下,在更广的传播方向内和更低的频率下,能产生较宽的方向带隙.      

Design of radial sonic crystal for sound attenuation from divergent sound source

Sonic crystals are periodic arrangement of sound hard scatterers, generally in square or triangular lattice configuration. The periodic obstruction to the sound wave by the scatterers leads to an interesting phenomenon of the band gap, which results in a high sound attenuation in the band gap region. In this work, a design of sonic crystal called as the radial sonic crystal is presented, which consists of periodic structures in polar coordinates. Such a structure attenuates divergent sound source. The radial sonic crystal is designed based on the Webster horn equation and using the property of invariance of governing equation from one unit cell to another. The designed radial sonic crystal is tested experimentally and by the finite element simulation. The experimental results are in good agreement with the simulation and show high sound attenuation of 30 dB. The high sound attenuation is due to the presence of the band gap in the radial sonic crystal.     

弯曲波彩虹捕获效应及其在能量俘获中的应用

弹性波在色散关系经过设计的梯度结构中传播时会产生空间分频现象和波场能量增强现象,即不同频率的弹性波会在结构的不同位置停止向前传播并发生能量聚集,这就是弹性波彩虹捕获效应.其相关研究成果可以促进结构健康监测、振动控制以及能量俘获等领域的发展.本文通过所设计的梯度结构梁,系统地研究了弯曲波彩虹捕获效应及其在压电能量俘获中的应用.首先,利用传递矩阵法获得了梯度结构梁元胞能带结构的解析解,进而分析了弯曲波彩虹捕获效应的产生机理:不同频率的弯曲波会在不同元胞附近群速度减小到零,从而停止向前传播并发生反射;入射波和反射波的叠加,以及群速度减小带来的能量聚集,会显著增强反射处的波场能量.其次,通过有限元仿真和实验验证了弯曲波彩虹捕获效应的空间分频现象和波场能量增强现象.最后,通过有限元多物理场耦合仿真和实验,研究了粘贴PVDF压电薄膜的梯度结构梁相对于均匀梁的弯曲波能量俘获效果及其随入射波频率的变化规律.结果表明,在弯曲波彩虹捕获效应发生频带内,粘贴在梯度结构梁上的PVDF压电薄膜的输出电压约为粘贴在均匀梁相应位置处的PVDF压电薄膜的输出电压的2倍.     

夹层圆柱壳中弹性波传播的辛特性分析

论文研究了正交各向异性夹层圆柱壳中轴对称自由简谐波的传播问题.通过对变量合理的组织变换,将结构本构方程化为状态空间形式,采用分段平均假设得到哈密顿矩阵,进而利用哈密顿系统下的辛数学方法,扩展的Wittrick-Williams算法及精细积分方法,得到各种夹层结构波传播问题的频散关系,并将该方法与多项式方法进行对比,验证了该方法在多孔结构波传播问题中的优越性.     

Local resonance phononic band gaps in modified two-dimensional lattice materials

In this paper,modified two-dimensional periodic lattice materials with local resonance phononic bandgaps are designed and investigated.The design concept isto introduce some auxiliary structures into conventional periodic lattice materials.Elastic wave propagation in this kindof modified two-dimensional lattice materials is studied using a combination of Bloch’s theorem with finite elementmethod.The calculated frequency band structures of illustrative modified square lattice materials reveal the existenceof frequency band gaps in the low frequency region due tothe introduction of the auxiliary structures.The mechanismunderlying the occurrence of these frequency band gaps isthoroughly discussed and natural resonances of the auxiliarystructures are validated to be the origin.The effect of geometric parameters of the auxiliary structures on the width ofthe local resonance phononic band gaps is explored.Finally,a conceptual broadband vibration-insulating structure basedon the modified lattice materials is designed and its capability is demonstrated.The present work is anticipated to beuseful in designing structures which can insulate mechanicalvibrations within desired frequency ranges.     

PWAS激励Lamb波的谱有限元建模及参数研究

基于谱有限元法发展了一种由压电晶片主动传感器(PWAS)、胶层和主结构组成的三层模型,来模拟PWAS激励结构中Lamb波的传播。首先在各层使用不同的梁理论,推导PWAS-胶层-主结构三层模型的控制方程和力的边界条件,建立谱有限元模型。通过和传统的有限单元法进行比较,表明了在显著提高计算效率的同时,所发展谱有限元模型在分析结构中Lamb波传播上仍具有较高的精度。分析了激励频率、PWAS长度与厚度、胶层厚度等参数变化对输出电压信号的影响,可以为基于PWAS和Lamb波的主动健康监测技术提供参考。     

Wave propagation in double helical rods

The propagation of waves in helical rods has been studied extensively. However, studying the wave propagation in double helical rods have received less attention although this can be useful in multiple fields of science and engineering. Obtaining an analytical model for a double helical rod is challenging since the curvature and tortuosity are not constant. Thus, resolving the wave behaviour analytically is nearly impossible. In this paper, wave propagation in a double helical rod will be studied using the wave and finite element method which is a technique that can be used to model homogeneous and periodic one and two dimensional structures based on the periodic structure theory. For modelling a double helical rod, the finite element model of a single turn is processed using Bloch waves. The dispersion curves and wavemodes are obtained and the similarities and differences of waves in helical and double helical rods are highlighted.     

Vibration modelling of structural networks using a hybrid finite element/wave and finite element approach

The vibration modelling of waveguide structures is considered. These structures comprise waveguides connected via joints. Traditionally, analytical models of the wave behaviour of such structures can be developed if they are simple (beams or rods connected at point joints, etc.). However, if the waveguides are of complicated constructions (truss-like, layered media, etc.) or the joints are complicated (e.g. of significant physical dimensions), obtaining the wave characteristics might be a formidable task. In this paper, such structures are modelled using a hybrid finite element/wave and finite element (FE/WFE) approach. The waveguides are modelled using the WFE method and thus their wave characteristics are obtained regardless of the complexity of their cross-section. The joints are modelled using standard FE, and the WFE and FE models are coupled to yield the scattering properties of the joints. The propagation and scattering models are assembled to describe the behaviour of the structure using relatively small models, while also providing information for other applications such as structure-borne sound, statistical energy analysis, etc. Numerical examples are presented to illustrate the approach.     

三维非均匀介质中弹性波传播的数值模拟

提出了一种三维非均匀介质中弹性波传播数值模拟的方法,文中称为三维格子法。该算法是二维格子法（一种二维非均匀介质中P－SV波传播的数值模拟算法）向三维非均匀介质情况的推广。在空间离散上该文方法与有限元方法类似,容许根据连续体的形状和介质分界面任意剖面网格,且自然满足自由表面边界条件。不同于常规有限差分法在各个节点上满足动力学微分方程,该算法通过满足各节点周围格子的整体平衡（积分平衡方程）来对问题进行求解,三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当。算例表明,该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射和绕射。     

阶梯压电层合梁的波动动力学特性

采用行波理论系统地研究了压电阶梯梁的自由振动分析以及强迫响应的分析方法. 基于分布 参数理论研究了压电阶梯梁的波传播特性，忽略柔性梁横向剪切和转动惯量的影响，给出了 梁的轴向和横向的简谐波解. 将压电阶梯梁离散化为单元，考虑压电片的刚度和质量的影响， 建立了节点散射模型. 应用位移连续和力平衡条件，推导了节点的波反射和波传递矩阵，在 此基础上，引入波循环矩阵的概念，给出波循环矩阵、波传递系数矩阵的确定方法. 应用波 循环矩阵可以有效地计算结构的固有频率. 另外，应用波传递系数研究了压电陶瓷作动器位 置对其驱动能力的影响. 得出两个主要结论：1)作动器靠近悬臂梁固定端将有较强的驱动 能力，悬臂梁边界反射行波产生弯曲消失波有利于增大压电波的模态传递系数；2)模态传 递系数与固有频率的灵敏度密切相关，波传递系数越大, 对应该处固有频率变化灵敏度越大. 另外，数值算例表明了行波方法比有限元方法具有更高的计算精度.     

Propagation of a stress wave through a virtual functionally graded foam

Stress wave propagation through a Functionally Graded Foam Material (FGFM) is analysed in this paper using the finite element method. A finite element model of the Split Hopkinson Pressure Bar (SHPB) is developed to apply realistic boundary conditions to a uniform density foam and is validated against laboratory SHPB tests. Wave propagation through virtual FGFMs with various gradient functions is then considered. The amplitude of the stress wave is found to be shaped by the gradient functions, i.e., the stress can be amplified or diminished following propagation through the FGFMs. The plastic dissipation energy in the specimens is also shaped by the gradient functions. This property of FGFMs provides significant potential for such materials to be used for cushioning structures.     

动力学平衡方程的辛两步求解算法

基于线性多步方法的构造格式和辛变换,给出了动力学方程的两种辛两步法求解格式,它们分别具有四阶精度和二阶精度,但都只有二阶格式的计算量,因此四阶辛两步法具有较大的应用价值。对两种辛两步法和解析解进行了数值比较,证明了二阶精度辛两步格式在一定条件下就是欧拉中点保辛算法,或δ=0.5和α=0.25的Newmark辛格式。     

截锥型薄壁结构声振耦合动力特性分析

采用大型通用软件ANSYS,建立截锥型薄壁结构的实体有限元动力学模型,通过与相关实验数据的对比验证了模型合理性。据此,利用无限元模拟自由声场边界,建立声场-截锥型薄壁结构的直接耦合有限元动力学模型。通过数值仿真分析研究了声场中截锥壳结构的振动特性,并讨论了声振动对结构动力特性的影响。研究结果表明:数值仿真结果和截锥壳声振实验数据比较一致。在考虑声场影响后发现:结构位移共振频率值大多有所降低,结构位移共振频率数量显著增多;在低频下,结构位移响应峰值在声场的影响下明显增大;在高频下则明显减小。     

Diffraction of waves by semi-infinite breakwater using finite and infinite elements

A finite and infinite element model is derived to predict wave patterns around a semi-infinite breakwater in water of constant depth. Both circular and square meshes of elements are used. The wave theory used is that of Berkhoff. The appropriate boundary conditions for finite and infinite boundaries are described. The singularity in the velocity at the breakwater tip is modelled effectively using the technique of Henshell and Shaw originally developed in elasticity. The results agree well with the analytical solution. In addition the problem of waves incident upon a semi-infinite breakwater and parabolic shoal, where both diffraction and refraction are present, is solved. There is no analytical solution for this case. The combination of finite and infinite elements is found to be an effective and accurate technique for such problems.     

Optimum design of band-gap beam structures

The design of band-gap structures receives increasing attention for many applications in mitigation of undesirable vibration and noise emission levels. A band-gap structure usually consists of a periodic distribution of elastic materials or segments, where the propagation of waves is impeded or significantly suppressed for a range of external excitation frequencies. Maximization of the band-gap is therefore an obvious objective for optimum design. This problem is sometimes formulated by optimizing a parameterized design model which assumes multiple periodicity in the design. However, it is shown in the present paper that such an a priori assumption is not necessary since, in general, just the maximization of the gap between two consecutive natural frequencies leads to significant design periodicity.The aim of this paper is to maximize frequency gaps by shape optimization of transversely vibrating Bernoulli–Euler beams subjected to free, standing wave vibration or forced, time-harmonic wave propagation, and to study the associated creation of periodicity of the optimized beam designs. The beams are assumed to have variable cross-sectional area, given total volume and length, and to be made of a single, linearly elastic material without damping. Numerical results are presented for different combinations of classical boundary conditions, prescribed orders of the upper and lower natural frequencies of maximized natural frequency gaps, and a given minimum constraint value for the beam cross-sectional area.To study the band-gap for travelling waves, a repeated inner segment of the optimized beams is analyzed using Floquet theory and the waveguide finite element (WFE) method. Finally, the frequency response is computed for the optimized beams when these are subjected to an external time-harmonic loading with different excitation frequencies, in order to investigate the attenuation levels in prescribed frequency band-gaps. The results demonstrate that there is almost perfect correlation between the band-gap size/location of the emerging band structure and the size/location of the corresponding natural frequency gap in the finite structure.     

The mutual variational principle of free wave propagation in periodic structures

By taking infinite periodic beams as examples, the mutual variational principle for analyzing the free wave propagation in periodic structures is established and demonstrated through the use of the propagation constant in the present paper, and the corresponding hierarchical finite element formulation is then derived. Thus, it provides the numerical analysis of that problem with a firm theoretical basis of variational principles, with which one may conveniently illustrate the mathematical and physical mechanisms of the wave propagation in periodic structures and the relationship with the natural vibration. The solution is discussed and examples are given. Supported by Doctorate Training Fund of National Education Commission of China     

Numerical modeling of PZT-induced Lamb wave-based crack detection in plate-like structures

This paper presents the numerical modeling and simulations of PZT-induced Lamb wave propagation in plate-like structures by using the spectral finite element method. A novel spectral plate finite element, which can efficiently model the three-dimensional (3D) behavior of Lamb waves, is proposed. In the formulation, linear displacement distributions in the thickness direction are assumed for both the PZT layer and the base plate. A way to avoid the thickness locking is proposed and used in the formulations. Two examples, one for the validation of the proposed two-dimensional (2D) spectral finite element and the other for the demonstration of crack detection in plates, are presented and discussed. The contact between the two faces of crack is considered. Numerical results show that (1) only the anti-symmetric mode is prone to thickness locking thus remedy should be made only on this part, (2) the proposed 2D spectral finite element can adequately model the Lamb wave propagation in plate-like structures and the complex scattering for the crack, and (3) crack location can be well determined by a PZT-induced Lamb wave-based diagnosis algorithm.     

Determination of stress intensity factors by the finite element discretized symplectic method

When rewriting the governing equations in Hamiltonian form, analytical solutions in the form of symplectic series can be obtained by the method of separation of variable satisfying the crack face conditions. In theory, there exists sufficient number of coefficients of the symplectic series to satisfy any outer boundary conditions. In practice, the matrix relating the coefficients to the outer boundary conditions is ill-conditioned unless the boundary is very simple, e.g., circular. In this paper, a new two-level finite element method using the symplectic series as global functions while using the conventional finite element shape functions as local functions is developed. With the available classical finite elements and symplectic series, the main unknowns are no longer the nodal displacements but are the coefficients of the symplectic series. Since the first few coefficients are the stress intensity factors, post-processing is not required. A number of numerical examples as well as convergence studies are given.