正交各向异性板中非主轴方向的Lamb波

基于线性三维弹性理论,采用勒让德正交多项式展开法,推导了波沿正交各向异性材料非主轴方向传播时的Lamb波耦合波动方程,并对耦合波动方程进行了数值求解。为验证该方法的适用性和正确性,首先将此方法应用于各向同性材料,并与已知的数据结果进行了比较;然后以单向纤维增强复合材料为例,计算了耦合Lamb波沿不同的非主轴方向传播时的相速度频散曲线,并分别研究了传播方向改变时低阶模态Lamb波和高阶模态Lamb波频散特性的变化。最后,针对潜在用于各向异性复合材料结构健康监测的耦合Lamb波低阶模态,给出了其在不同传播方向时的相速度分布和群速度分布。同时,结合低阶模态Lamb波的位移分布特性和材料的各向异性特点,阐释了S0模态对波的传播方向变化最为敏感的原因。     

连续体结构动力拓扑优化中局部模态处理的新方法

采用固体各向同性材料惩罚模型（solid isotropic material with penalization, SIMP） 进行动力拓朴优化通常在优化过程中会出现虚假的局部振动模态,为消除这种虚假模态产生的不利影响,提出了移频与虚假模态识别相结合的通用方法. 研究中考虑以材料体积为约束、结构基频最大化为目标的优化模型,并采用节点设计变量描述设计域内材料分布. 基于虚假模态的特性,首先在特征值分析中应用移频方法排除特征值接近于零的低阶虚假模态,然后再依据虚假模态识别准则判定并剔除其他可能存在的虚假模态,从而可以高效可靠地确定结构真实的固有振动模态. 数值算例表明,提出的方法可以有效地消除动力拓扑优化中虚假模态可能产生的不利影响,并保证优化解的可靠性.      

PROPAGATION CHARACTERISTICS OF HIGH ORDER LONGITUDINAL MODES IN STEEL STRANDS AND THEIR APPLICATIONS

Propagation characteristics of high order longitudinal modes of ultrasonic guided waves in seven-wire steel strands are investigated theoretically and experimentally. According to these analysis results, proper longitudinal modes are selected for defect detection in steel strands. Dispersion curves for helical and central wires in a 17.80 mm nominal diameter seven-wire steel strand are numerically obtained firstly, and propagation characteristics of high-order longitudinal modes, such as wave structures, attenuation and dispersion, are analyzed. In experiments, the signals of ultrasonic guided wave at different high frequencies are excited and received at one end of a steel strand by using the same single piezoelectric transducer. The identification of longitudinal modes in the received signals is achieved based on short time Fourier transform. Furthermore, appropriate L（0, 5） mode at 2.54 MHz is chosen for detecting an artificial defect in a helical wire of the steel strand. Results show that high order longitudinal modes in a high frequency range with low dispersion and attenuation whose energy propagates mainly in the center of the wires can be used for defect detection in long range steel strands.     

The initialization of nonlinear waves using an adjustment scheme

A procedure for initializing nonlinear free-surface simulations is developed and validated. Numerical simulations of nonlinear progressive waves are prone to developing spurious high-frequency standing waves unless the flow field is given sufficient time to adjust. An adjustment procedure is developed that allows nonlinear free-surface simulations to be initialized with linear solutions. The adjustment scheme allows the natural development of nonlinear self-wave (locked modes) and wave–wave (free modes) interactions. The implementation of the adjustment procedure is illustrated using a high-order spectral method. Comparisons are made to fully-nonlinear Stokes waves and Schrodinger theory.     

Selective enhancement of oblique waves caused by finite amplitude second mode in supersonic boundary layer

Nonlinear interactions of the two-dimensional(2D) second mode with oblique modes are studied numerically in a Mach 6.0 flat-plate boundary layer, focusing on its selective enhancement effect on amplification of different oblique waves. Evolution of oblique modes with various frequencies and spanwise wavenumbers in the presence of 2D second mode is simulated successively, using a modified parabolized stability equation(PSE) method, which is able to simulate interaction of two modes with different frequencies efficiently. Numerical results show that oblique modes in a broad band of frequencies and spanwise wavenumbers can be enhanced by the finite amplitude 2D second mode instability wave. The enhancement effect is accomplished by interaction of the 2D second mode, the oblique mode, and a forced mode with difference frequency. Two types of oblique modes are found to be more amplified, i.e., oblique modes with frequency close to that of the 2D second mode and low-frequency first mode oblique waves. Each of them may correspond to one type of transition routes found in transition experiments. The spanwise wavenumber of the oblique wave preferred by the nonlinear interaction is also determined by numerical simulations.     

Lamb wave mode decomposition using concentric ring and circular piezoelectric transducers

Lamb waves using surface-bonded piezoelectric transducers (PZTs) have been widely used for nondestructive testing (NDT). However, the identification of individual Lamb wave modes and the subsequent data interpretation are often difficult due to the dispersive and multimodal natures of Lamb waves. To tackle the problem, a Lamb wave mode decomposition technique using concentric ring and circular PZTs is proposed. Its advantages over the conventional approaches are that (1) PZTs need to be placed only a single surface of a specimen and (2) mode decomposition can be performed at any desired frequency without changing the PZT size and/or spacing configuration. The proposed mode decomposition technique is formulated by solving 3D Lamb wave propagation equations considering the PZT size and shape, and this technique requires a specially designed dual PZT composed of concentric ring and circular PZTs. The effectiveness of the proposed technique for the Lamb wave mode decomposition is investigated through numerical simulation and experimental tests performed on an aluminum plate.     

超声速平面剪切层声辐射涡模态数值分析

对Mc = 1.2二维超声速空间发展平面自由剪切层, 进行了扰动模态及流动结构的数值分析. 采用时空三阶改进MacCormack格式, 差分求解可压缩扰动Navier-Stokes方程, 直接数值模拟入口不同基频谐波扰动的非线性演化特征. 采用空间线性稳定性理论证明, 计算所促发的扰动波是声辐射涡模态. 扰动参数及特征函数分析显示, 声辐射涡模态是弱色散的快／慢两种外部模态, 在扰动对流Mach数为超声速一侧呈膨胀／压缩状辐射. 单频受迫扰动可无相差地促发多模态混合扰动波, 而在自然扰动条件下, 剪切层的稳定性受慢模态主导.     

基于裂纹附加模态Timoshenko梁的裂纹损伤识别

将梁中横向裂纹等效为无质量扭转弹簧，并忽略其对梁剪切变形的影响，得到的具有任意裂纹数目Timoshenko 梁自振模态的统一显示解析表达式．将裂纹梁的自振模态分为基本模态和裂纹附加模态，利用最小二乘拟合，建立了利用裂纹附加模态函数的梁裂纹损伤识别方法．通过数值模拟开展了简支单裂纹梁以及悬臂和固支双裂纹梁等的裂纹损伤识别，考察了测量误差对损伤识别的影响，数值结果表明本文所提出的裂纹损伤识别方法对裂纹位置的识别精度高于对裂纹损伤程度的识别精度；随着测量误差的增加，裂纹位置及裂纹损伤程度的识别误差增加，但仍在可接受的范围内，故该裂纹损伤识别方法在实际工程中具有一定的应用价值．     

Prediction on dispersion in elastoplastic unsaturated granular media

This study focuses on the propagation of the plane wave in the elastoplastic unsaturated granular media, and the wave equations and dispersion equations are derived for the media under the framework of Cosserat theory. Due to symmetry, five different wave modes are considered and predicted for the elastoplastic unsaturated granular media based on the Cosserat theory, including two longitudinal waves, one rotational longitudinal wave and two coupled transverse–rotational transverse waves. The correspondence is discussed between these Cosserat wave modes and the classical wave modes. Based on the dispersion equations, the dispersion behaviors are obtained for the five Cosserat wave modes. The results indicated that the different stress-strain stages,including the elastic, hardening and softening stages, have obvious effect on the dispersion behaviors of the Cosserat wave modes.     

Dynamic modeling and modal truncation approach for a high-speed rotating elastic beam

Summary  In the classical finite element analysis of beams, the nonlinear terms of deformation are ignored due to the linearization of deformation based on the assumptions of structural dynamics. Since the number of generalized coordinates is large in flexible bodies when using the finite element method (FEM), the modal truncation approach (MTA) is usually used for improving computational efficiency, and only lower-order transverse modes are chosen. In this paper, dynamic modeling and application of the MTA to a high-speed rotating beam are studied. The foreshortening displacement is included in the longitudinal displacement, therefore the dynamic modeling takes account of the effect of geometric nonlinearity. Equations of a rotating beam are obtained and the FEM and MTA are used for discretization. The applicability of the MTA to a high-speed rotating elastic beam is verified. The comparison of the results obtained by the FEM and MTA shows that in the case of a high-speed rotation, the centrifugal force can excite high-order transverse modes. Since using lower-order transverse modes for modal truncation obviously can cause error, addition of more transverse modes may improve the result. Furthermore, a coupling effect between axial and transverse displacements is revealed. It is shown that in the case of a sudden change of the axial displacement, the inclusion of the axial modes can significantly improve the response. Received 10 April 2001; accepted for publication 26 March 2002 This work was supported by the National Science Foundation of China (19832040) and the National Education Ministry of China (2000024818), for which the authors are grateful.     

薄板弯曲分析的高阶高效无网格法

与传统有限元法相比,无网格法具有节点形函数高度光滑、易于形成高阶近似等优势,更适合于以薄板弯曲问题为代表的高阶偏微分方程的数值求解。然而,高阶无网格法的形函数是非多项式的有理函数,导致弱形式的区域积分难以得到精确计算,通常采用的高阶高斯积分方法需使用大量积分点,计算效率低且精度不高。本文针对薄板弯曲问题的高阶（三阶）无网格法分析,首次发展了与该高阶近似相一致的曲率光顺方案,并基于背景三角形积分单元建立了相应的数值积分格式,大幅度减少了所需的积分点数目。所发展方法的关键在于计算刚度阵所需的形函数的二阶导数由形函数及其一阶导数通过散度定理确定,而非对形函数直接求导获得。数值结果表明,基于标准的高斯积分方案的高阶无网格法精度不高,不能精确再现纯弯曲和线性弯曲模式,且得到的弯矩场分布存在严重的虚假数值振荡。而本文所建议的基于曲率光顺方案的高阶无网格法能够方便高效地求解薄板弯曲问题,尤其是它能精确反映纯弯曲和线性弯曲模式。与标准的高斯积分方法和目前主流的常曲率光顺方法相比,本文方法在计算效率、精度、弯矩分布等方面均展现出显著优势,因而具有较好的应用价值。     

Dynamics of cantilever plates and its hybrid vibration control

Applying Lagrange–Germain’s theory of elastic thin plates and Hamiltonian formulation, the dynamics of cantilever plates and the problem of its vibration control are studied, and a general solution is finally given. Based on Hamiltonian and Lagrangian density function, we can obtain the flexural wave equation of the plate and the relationship between the transverse and the longitudinal eigenvalues.Based on eigenfunction expansion, dispersion equations of propagation mode of cantilever plates are deduced. By satisfying the boundary conditions of cantilever plates, the natural frequencies of the cantilever plate structure can be given.Then, analytic solution of the problem in plate structure is obtained. An hybrid wave/mode control approach, which is based on both independent modal space control and wave control methods, is described and adopted to analyze the active vibration control of cantilever plates. The low-order(controlled by modal control) and the high-order(controlled by wave control) frequency response of plates are both improved. The control spillover is avoided and the robustness of the system is also improved. Finally, simulation results are analyzed and discussed.     

Waves, normal modes and frequencies in periodic and near-periodic rods. Part I

Longitudinal wave motions and localized normal modes in a rod system with periodically-alternating material properties are investigated in this paper. The energy injected into the rod system is shown either to be transported through the whole rod system in pass-bands or to be trapped near the excitation source in stop-bands. For this one-dimensional continuous model, the full power of linear system theory is utilized and a new transfer matrix method is proposed to get closed-form normal mode solutions. Localized normal modes in stop-bands in perfectly-periodic rods with asymmetric bays are identified. It is shown that for this strongly-coupled elastic system, a single small disorder may produce one or two additional modes in each stop-band, these modes are localized around the disordered bay. By understanding basic behavior of such a system, it is hoped ultimately that some insights can be achieved where closed-form results are not possible.     

Non-linear modal analysis of a rotating beam

The free non-linear vibration of a rotating beam has been considered in this paper. The von Karman strain-displacement relations are implemented. Non-linear equations of motion are obtained by Hamilton’s principle. Results are obtained by applying the method of multiple scales to a set of discretized ordinary differential equations which obtained by using the Galerkin discretization method. This set contains coupling between transverse and axial displacements as quadratic and cubic geometric non-linearities. Non-linear normal modes and non-linear natural frequencies with or without internal resonance are observed. In the internal resonance case, the internal resonance between two transverse modes and between one transverse and one axial mode are explored. Obtained results in this study are compared with those obtained from literature. The stability and some dynamic characteristics of the non-linear normal modes such as the phase portrait, Poincare section and power spectrum diagrams have been inspected. It is shown that, for the first internal resonance case, the beam has one stable or degenerate uncoupled mode and either: (a) one stable coupled mode, (b) one unstable coupled mode, (c) two stable and one unstable coupled modes, (d) three stable coupled modes, and (e) one stable coupled mode. On the other hand, for the second internal resonance case, the beam has one stable or unstable or degenerate uncoupled mode and either: (a) two stable coupled modes, (b) two unstable coupled modes, and (c) one stable coupled mode depending on the parameters.     

Lamb wave scattering by a symmetric pair of surface-breaking cracks in a plate

The scattering problem of a Lamb wave incident on a symmetric pair of surface-breaking transverse cracks in a plate is considered. The Lamb wave is assumed to be obliquely incident on the crack plane. Since the cracks are part-through, the scattered field will contain reflected as well as transmitted waves. The energy of the incoming wave is partitioned into reflected and transmitted wave modes. Energy coefficients of the reflected and transmitted waves are calculated as a function of incident frequency and crack depth. The incidence angle of the incoming wave is also treated as a parameter. Both the reflected and transmitted wave fields are considered as linear superpositions of all real and complex wave modes in the plate. Decomposition of modes is achieved with the help of an orthogonality condition based on the principle of reciprocal work. Continuity of displacement and stress fields is imposed at the crack plane. Energy coefficients for reflection and transmission are obtained from the mode amplitudes. Energy coefficients are shown to be a strong function of incident frequency and crack depth. Experiments are conducted with a PZT transducer network interacting with a symmetric pair of machined cracks in an aluminum plate. Trends predicted by the analysis are reflected in the experimental results.     

Effect of initial deformations on wave front propagation in anisotropic plates

The effects of large amplitudes and initial deformations on shock waves and acceleration waves propagating in fiber-reinforced laminated plates are investigated. Three cases are discussed, namely the large amplitude shock under initial in-plane deformations, small amplitude waves under in-plane deformations, and small amplitude waves propagating in a plate with large deflection. It is found that the in-plane force has a substantial effect on the transverse shear mode but little effects on other modes. The large initial deflection, however, is found to have considerable effects on all modes. A general procedure for constructing the wave surfaces is also presented.     

Internal waves from a body accelerating in a thermocline

Many papers study the steady wave system around bodies moving in thermoclines but little attention has been given to unsteady wave systems. This paper concentrates on the unsteady wave systems around accelerating bodies in thermoclines. The wave shapes are calculated using a theory derived from a dispersion relation based on an exp-tanh density profile. All modes of oscillation can be determined and it is shown that for the lowest mode both oblique and transverse waves occur whereas for the higher modes the presence of transverse waves depends on the background conditions and on the speed of the body. Cauchy-Poisson impulsive start waves are included. The theoretical wave shapes compare quite well with those calculated using finite-difference formulations of the full Navier-Stokes equations when a body accelerates from rest.It is also shown how the dispersion relation =N sin together with the WKB approximation can produce the same plan-view wave forms as those obtained using the thermocline wave dispersion relation given by [17, 30].     

深海采矿系统中悬臂式立管涡激振动分析

不同于传统的海洋立管, 深海采矿系统中的垂直提升管道可以被视为一个底部无约束的柔性悬臂式立管, 工作过程中同样面临涡激振动和柔性变形问题. 本文采用一种无网格离散涡方法和有限元耦合的准三维时域求解数值模型, 系统性地研究了不同流速下悬臂式立管的涡激振动问题. 结果表明: 悬臂式立管的横向振动模态阶数随折合速度增加而增大, 在一定折合速度范围内主导振动模态保持不变; 当主导模态转变时, 对应的横向振幅会发生突降, 但是当新的高阶模态被激发后, 立管振幅随来流速度增加而再次逐渐增大; 在相同的振动模态下, 立管底部位移均方根值随折合速度线性增加, 主导振动频率在模态转变时会出现跳跃现象; 特别地, 本文讨论了三阶主导模态下悬臂式立管的振动响应, 无约束的立管底部呈现出较大的振动能量, 且振幅的驻波特征随折合速度增加而逐渐增强; 本文比较了两端铰支立管与悬臂式立管的涡激振动响应特征, 两者在振幅和主导振动频率两方面均表现出了相同的变化趋势.      

Nonlinear coupled vibration of electrostatically actuated clamped–clamped microbeams under higher-order modes excitation

Nonlinear modal interactions have recently become the focus of intense research in micro-resonators for their use to improve oscillator performance and probe the frontiers of fundamental physics. Understanding and controlling nonlinear coupling between vibrational modes is critical for the development of advanced micromechanical devices. This article aims to theoretically investigate the influence of antisymmetry mode on nonlinear dynamic characteristics of electrically actuated microbeam via considering nonlinear modal interactions. Under higher-order modes excitation, two nonlinear coupled flexural modes to describe microbeam-based resonators are obtained by using Hamilton’s principle and Galerkin method. Then, the Method of Multiple Scales is applied to determine the response and stability of the system for small amplitude vibration. Through Hopf bifurcation analysis, the bifurcation sets for antisymmetry mode vibration are theoretically derived, and the mechanism of energy transfer between antisymmetry mode and symmetry mode is detailed studied. The pseudo-trajectory processing method is introduced to investigate the influence of external drive on amplitude and bifurcation behavior. Results show that nonlinear modal interactions can transit vibration energy from one mode to nearby mode. In what follows, an effective way is proposed to suppress midpoint displacement of the microbeam and to reduce the possibility of large deflection. The quantitative relationship between vibrational modes is also obtained. The displacement of one mode can be predicted by detecting another mode, which shows great potential of developing parameter design in MEMS. Finally, numerical simulations are provided to illustrate the effectiveness of the theoretical results.     

低速冲击下复合材料合板的响应过程模拟

运用一种精度较高的高阶位移模型分析了复合材料层合板在低速击下的响应过程,该位移模型能够同时考虑层间正应力和横向剪应力,另外,采用修正的Hertz接触定律与Newmark积分方法相结合,建立了冲击接触的有限元模型,本文的数值模拟计算结果与解析结果相比较证明了该方法的有效和精确性。