复杂结构冲击声合成及实验验证

以复杂结构受击振动响应的时域计算为目的, 讨论了结构阻尼的计算方法, 给出一种用于冲击声合成的综合数值方法, 并进行了实验验证. 首先, 考虑到阻尼是影响瞬态振动时变特性的重要因素, 详细讨论了两种模态阻尼的计算方法; 其次, 对阻尼板的受击振动和声辐射进行了时域仿真, 并与时域有限差分法的计算结果进行对比, 显示出两种声音合成方法的计算结果具有高度的一致性; 最后, 针对有限长圆柱壳的受击振动, 将合成声与实验录音进行了对比研究. 结果表明, 合成声与实际录音的时域包络、频谱结构以及衰减趋势基本一致, 证明了采用数值方法进行冲击声合成的有效性. 关键词： 声音合成 模态阻尼 冲击声 数值方法     

传统笙的物理模型及声音合成方法

为研究传统笙的物理结构与音色之间的关系,提出了传统笙的物理模型及声音合成方法。在簧舌振动模型的基础上,结合指孔和音窗的传输矩阵计算方法,建立了笙管的等效电路,并将二者结合,得到笙的完整物理模型。为了验证模型的有效性,设计了一套笙的实验系统,实现了对笙管发声过程中特征量的多通道同步测量。通过对比实验和模型仿真结果的时频域特点以及与音色相关的特征量,分析了模型的性能.结果表明,该模型可较好地模拟传统笙管的发声过程,合成笙样本能够较准确的反映真实笙样本的音色特征。      

Structural damping identification based on an iterative regularization method

A time domain structural damping identification method is presented in this paper. The damping of a structure is assumed to be time-invariant Rayleigh damping, time-variant Rayleigh damping and modal damping, respectively, and an iterative regularization method is proposed to identify these three types of damping in turn. A sensitivity approach is adopted in the formulation of the ill-posed inverse problem. A new constraint is imposed on the identified iterative increment as well as on the unknown structural parameters to ensure their physical meaning in the identification process is not lost. The proposed method is verified in numerical studies with a space frame structure and laboratory measurements of accelerations with accurate results.     

Estimation of interfacial forces in time domain for linear systems

A structural path rank ordering process under transient excitations requires a good knowledge of the interfacial path forces, which are difficult to directly measure. Four time domain methods to estimate the interfacial forces are proposed and comparatively evaluated with application to linear time-invariant, proportionally damped discrete systems. First, the transient response is derived by modal analysis and a direct time domain technique to calculate the interfacial forces is outlined. Next, the frequency domain estimation methods, based on the sub-system concept are reviewed, and an inverse Fourier transform scheme is introduced. An indirect method of estimating interfacial force in transient state is then developed through an inverse procedure of modal analysis. The sub-system approach is employed to obtain the interfacial forces based on the forced vibration response of the original system and modal data of the sub-system. Finally, an approximate time domain scheme is suggested that could be used only if the system properties are known or precisely estimated. Although the proposed indirect methods are designed for eventual experimental applications, this article provides numerical feasibility studies via a simple source-path-receiver system (with parallel vibration paths) that has five translational degrees of freedom. The proposed methods are compared under ideal impulse force excitation input and a periodic sawtooth load (without and with Gaussian noise) to observe the starting transients as well as subsequent motions and interfacial forces. Preliminary comparisons with a laboratory experiment are very promising.     

Efficient model order reduction for dynamic systems with local nonlinearities

In the nonlinear structural analysis, the nonlinear effects are commonly localized and the rest of the structure behaves in a linear manner. Considering this fact, this research work proposes a harmonic balance solution in order to determine the nonlinear response of the structures. The solution is simplified by using an exact dynamic reduction along with the modal expansion technique. This novel approach, which is applicable to both discrete and continuous systems, converts the system equations of motion in each harmonic to a small set of nonlinear algebraic equations. The full set of system equations is reduced to a discrete system with a few generalized degrees of freedom (DOFs) confined to the localized nonlinear regions. The resultant reduced order model is shown to be accurate enough for determining the periodic response. To demonstrate the capability of the proposed method, numerical case studies for continuous and discrete systems, including systems with internal resonance, have been studied and the outcomes are validated with benchmark studies. In addition, the method is applied in the identification process of an experimental test setup with unknown frictional support parameters, and the results are presented and discussed.     

Subjective evaluations of sound radiated by impacted plates, using the design of experiments method

This study combines physical and subjective approaches to evaluate the sound quality of impacted plates. A numerical model, based on the Hertz law of contact, has been used to determine the sound pressure level at any point in space resulting from an impact. Sounds synthesized using this model and those from experiments can then be exploited in a physical analysis and/or a subjective analysis of the effects associated with variations in parameters. Here the influence of certain physical parameters on the sound perception of impacted plates is evaluated through a design of experiments method and a subjective test of preference. The subjective test is based on an estimate of preference between two pairs of synthetic sounds by varying several structural parameters at the same time. This differs from other studies that vary only one parameter at a time and fix the other parameters.     

Expansion-chamber muffler for impulse noise of pneumatic frictional clutch and brake in mechanical presses

A compound expansion-chamber muffler, which consists of a sound absorbing chamber and a switch valve, the chamber integrating structural features of impedance muffler and micropunch plate muffler, is proposed to diminish impulse exhaust noise of pneumatic friction clutch and pneumatic friction brake (PFC/B) in mechanical presses. The structure decreases the impulse exhaust noise of PFC/B over 30 dB(A). A one-dimensional flow model is applied to study the aerodynamic characteristics of compound exhaust process of the single acting cylinder and muffler because the exhaust time is a critical factor for application of muffler in PFC/B. The volume of sound absorbing chamber is found to be an important design parameter to minimize the exhaust resistance of pneumatic cylinder. Experiments are also conducted to validate analytical results. Then the effects of diameter of exhaust ducts and volume of muffler on the exhaust time are discussed in detail. The proposed one-dimensional computational method, which considers the coupling of air-flow field and sound field, gives satisfactory results for the preliminary design of an expansion-chamber muffler. This method has been applied to an existing model HKM3-40MN to reduce its impulse noise.     

响应面法的杆式超声电机有限元模型修正

现有杆式超声电机的有限元模型大都采用连续复合材料的简化结构形式,忽略了实际超声电机各部件间的接触特性,导致计算结果与实验结果具有较大偏差。针对该问题,提出了一种采用响应面法的杆式超声电机有限元模型修正方法,来获得高精度的有限元模型。该方法考虑了超声电机螺栓预紧力及各部件接触界面的法向接触刚度和摩擦系数,筛选出显著影响电机工作模态的参数,建立响应面模型替代超声电机的有限元模型实现快速计算结构响应的目的,并以实验模态分析结果为目标对模型进行修正。修正结果表明,修正后的模型模态频率的平均误差由修正前的1.20%降到0.21%,模型精度得到明显改善,表明以响应面的有限元模型修正方法对杆式超声电机以及类似夹心式压电振子的设计具有应用价值。      

A modal test method using sound pressure transducers based on vibro-acoustic reciprocity

A modal test method that uses sound pressure transducers at fixed locations and an impact hammer roving over a test structure is developed in this work. Since sound pressure transducers are used, the current method deals with a coupled structural–acoustic system. Based on the vibro-acoustic reciprocity, the method is equivalent to one, where acoustic excitations at fixed locations are given and the resulting acceleration of the test structure is measured. The current method can eliminate mass loading due to use of accelerometers, which can destroy existence of repeated or close natural frequencies of a symmetric structure. It can also avoid effects of a nodal line of a mode and an inactive area of a local mode, and measure all the out-of-plane modes within a frequency range of interest, including global and local ones. The coupling between the structure and the acoustic field in a structural–acoustic system introduces asymmetry in the model formulation. An equivalent state space formulation is used for a damped structural–acoustic system and the associated eigenvalue problem is derived. The biorthonormality relations between the left and right eigenvectors and the relations between the structural and acoustic components in the left and right eigenvectors are proved. The frequency response functions associated with the current method are derived and their physical meanings are explained. The guidelines for using the current method, including the types of structures that are suitable for the method, the positions of the sound pressure transducers, and the orientation of the test structure relative to the transducers, are provided. Modal tests were carried out on an automotive disk brake using the traditional and current methods, where multiple accelerometers and microphones were used to measure its dynamic responses induced by impacts, respectively. The differences between the measured natural frequencies using the current method and those from the finite element model of the disk brake are less than 3 percent for the first 18 elastic modes, and the modal assurance criterion values of the associated mode shapes are all above 90 percent. The current method was also used to measure the natural frequencies, damping ratios, and mode shapes of a light circuit board.     

Finite element model updating of a rod-type ultrasonicmotor based on response surface method

The conventional finite element model(FEM) of a rod-type ultrasonic motor is usually simplified by means of continuous composite structure. Because the actual contact characteristics between the parts of the ultrasonic motor is ignored, there is bigger error between the calculated values and experimental results. Aiming at solving problem, a new modeling method of a rod-type ultrasonic motor is presented to obtain a high-accuracy FEM. The bolt pretension and the normal contact stiffness and friction coefficient of the contact surface of ultrasonic motor are all considered in this method, and the significant parameters of working mode of the motor are selected by the response surface method, and the goal of calculating the structural response rapidly is realized by building the response surface model to replace the FEM. The result of finite element model updating shows that the average error of modal frequencies of updated model drops to 0.21% from 1.20%. The accuracy of FEM is obviously improved, which indicates that the FEM updating based on response surface method is of great application value on the design for a rod-type ultrasonic motor.     

基于超二次曲面的颗粒材料缓冲性能离散元分析

自然界或工业中普遍是由非球形颗粒组成的复杂体系,与球形颗粒相比,非球形颗粒间的高离散和咬合互锁可使冲击载荷引起的能量有效衰减实现缓冲作用.基于连续函数包络的超二次曲面单元能准确地描述非球形颗粒的几何形态,并可精确地计算单元间的接触碰撞作用.本文采用离散元方法对冲击载荷作用下非球形颗粒物质的缓冲性能进行数值分析,并与圆柱体冲击的理论结果和球体冲击的实验结果进行对比验证.在此基础之上,进一步研究了筒底作用力在不同颗粒层厚度和形状等因素影响下的变化规律.计算结果表明:不同颗粒形状都存在一个临界厚度H_c.当HH_c时,缓冲率随H的增加而增加;当HH_c时,缓冲率的变化不再显著并趋于稳定值.此外,减小颗粒表面尖锐度和增加或减小圆柱形和长方形颗粒的长宽比都会提高颗粒材料的缓冲效果.     

三维绝热简正波-抛物方程理论及应用

复杂海域通常存在环境参数的水平变化,这会导致声波在传播过程中发生水平折射,呈现出三维效应.利用绝热简正波-抛物方程理论进行三维声场建模,在垂直方向上使用标准简正波模型KRAKEN求解本征值和本征函数,水平方向上使用宽角抛物方程模型RAM求解简正波幅度.该模型物理意义清晰,计算效率高,但由于忽略了各号简正波之间的耦合,只适用于环境参数水平变化缓慢的问题.使用该模型分析了内波环境和大陆架楔形波导中的声波水平折射现象,结果表明,声波的水平折射将水平平面分为不同区域,每个区域内的声场结构明显不同.此外,声强在水平平面内的分布与声源频率和简正波号数有关,这种依赖关系是导致声信号频谱变化、波形畸变以及声场时空扰动的主要原因.     

液滴撞击圆柱内表面的数值研究

针对液滴撞击圆柱内表面的过程,利用基于相场的格子Boltzmann方法模拟液滴以不同初速度、从不同初始高度、撞击不同大小的圆柱内表面时液滴的形态变化,分析了液滴自身物性(如密度和黏性等)和圆柱内表面润湿性等因素对撞击现象的具体影响.研究发现:撞击韦伯数、密度比及动力黏性比、圆柱半径等对液滴撞击后沿圆柱内表面的铺展均有一定影响,较高的韦伯数下液滴可能会发生分裂;液滴初始高度对大密度比和动力黏性比的撞击影响较小;液滴反弹现象可能出现在接触角较大时;重力作用会抑制撞击后液滴的振荡.     

考虑任意阻抗壁面条件管腔结构声场特性分析

针对任意阻抗壁面条件一维管腔声学系统建模,对系统动力学特性进行预报。为了满足阻抗边界条件对声压一阶导数连续性要求,管腔声压函数通过在标准傅里叶级数端点位置引入边界光滑辅助多项式进行构建。结合壁面阻抗声学边界和管腔声学Helmholtz控制微分方程得到强形式标准特征值问题,获得相应的声学模态信息。在数值仿真中,通过算例给出各种边界条件下管腔声学模态频率、声压振型、声压和质点振速频率响应曲线,与现有文献中相关结果进行对比,充分验证了本文求解方法的正确性和有效性,证明该方法可对任意阻抗壁面条件管腔系统声学特性进行准确预报。     

Acoustic characteristics analysis of slender cavity system with arbitrary impedance end conditions

A modeling method for the dynamic characteristics analysis of a slender acoustical cavity with impedance end conditions is established. In order to satisfy the continuity requirement at impedance ends for the first order differential of sound pressure, field function is constructed as the standard Fourier series supplemented by boundary smoothed auxiliary polynomials. System characteristic equation is derived by solving the governing differential equation and impedance acoustic boundary of slend...     

Modal analysis by free vibration response only for discrete and continuous systems

This work aims to develop the algorithm for modal analysis by free vibration response only (MAFVRO), in particular for the general or non-proportional viscous damping system model. If the structural displacement or acceleration response due to free vibration can be measured, the system response matrices, including the displacement, velocity and acceleration, can be obtained through numerical differential or integration methods. These response matrices can then be applied to the developed MAFVRO method to determine the structural modal parameters. The numerical differential and integration methods are introduced and adopted to establish the modal parameter prediction program for the non-proportional damping model of MAFVRO. This work also shows the applications of MAFVRO to the multiple degree-of-freedom (mdof) systems and the cantilever beam, respectively. Both the discrete and continuous systems are demonstrated for the feasibility of the MAFVRO algorithm. The developed method uses the free vibration output response only and can obtain the structural modal parameters successfully.     

Hybrid holographic-numerical method for modal analysis of complex structures

This paper presents a hybrid holographic-numerical method for modal analysis of complex structures. A continuous structure is first lumped into a number of discrete elements to form an elastically connected lumped linear system. The matrix of influence coefficients of the lumped linear system are then determined by exerting a static load to the element centers and measuring the corresponding whole-field displacement using digital holographic interferometry. The eigenvalues and eigenvectors of the influence coefficients, which in a physical sense represent the natural frequencies and mode shapes of the structure, are then calculated using the numerical method. A major advantage of the proposed hybrid method is that it is not necessary to know the Young's modulus, the Poisson's ratio of the material and the boundary conditions, as the displacement field measured by the optical method has automatically reflected the real boundary conditions and the material properties, which makes this method particularly useful for studying objects made from anisotropic materials such as composites. Another advantage of the proposed method is that structures of any complex and irregular shape will not increase the complexity of the characterization process. The proposed is also suitable for experimentally validating the modal analysis results from finite element method models.     

一种双正交心音小波的构造方法

为了提高小波分析在心音信号处理中的性能, 在分析小波构造理论的基础上, 构造了一种专门用于心音信号处理的小波基. 首先提出一种构造滤波器长度为偶数的紧支撑双正交小波的一般方法; 然后根据心音信号的特点, 讨论心音小波的构造原则和一种基于心音 小波族的心音信号合成模型, 并且在此基础上构造出心音小 波. 为了突出使用心音小波处理心音信号的先进性和实用性, 对心音小波进行了比较全面的理论和数值仿真分析. 实验结果表明, 相比常用的db, bior系列小波, 运用心音小波对心音信号进行处理, 能够获得更好的去噪效果、 更精确的心音分类信息以及更小的重构误差率, 为心音特征提取和身份识别的深入研究提供了一种新方法, 在表征心音个体特征的细节方面具有积极的意义. 本文根据应用对象设计专用小波的方法也为工程应用中小波基的选择提供了一种新途径. 关键词： 双正交小波 心音小波 构造方法 心音合成模型     

A hybrid method of modal synthesis using vibration tests

The assembly of structures along continuous boundaries presents great difficulty in the context of modal synthesis. In order to solve such problems, a method is proposed in which a hybrid model is defined reflecting the dynamic behavior of a structure loaded along a boundary. It is based on Weinstein's method and corresponds to a generalization of the impedance matrix method. Generalized boundary co-ordinates are defined from branch modes obtained by introducing mass loading along the boundary. Thus, the hybrid model can be derived from testing as a result of two independent modal identifications. The method permits high precision prediction of the influence of strong structural modifications. Thus, in the case of rectangular plates, it has been possible to find the modes of a cantilever plate and of a plate with stiffeners from the free modes.     

An efficient broadband coupled-mode model using the Hamiltonian method for modal solutions

A numerically efficient broadband, range-dependent propagation model is proposed, which incorporates the Hamiltonian method into the coupled-mode model DGMCM. The Hamiltonian method is highly efficient for finding broadband eigenvalues, and DGMCM is an accurate model for range-dependent propagation in the frequency domain. Consequently, the proposed broadband model combining the Hamiltonian method and DGMCM has significant virtue in terms of both efficiency and accuracy. Numerical simulations are also provided. The numerical results indicate that the proposed model has a better performance over the broadband model using the Fourier synthesis and COUPLE, while retaining the same level of accuracy.