一种全频散波方程的反问题

基于一种全频散波方程研究了对于谐波和波包的反问题。首先根据Mindlin理论建立了描述无耗散微结构线性固体中波传播模型一一一种全频散波方程,并讨论了其频散特性。然后基于该全频散波方程,提出了利用四种不同谐波的频率和相应波数确定波方程四个未知系数的反问题,并用严格的数学理论论证了此反问题。研究证明,通过测量同一种无耗散微结构线性固体中传播的四种不同谐波的频率和相应波数,在正常频散和反常频散情况下可唯一地确定波方程的未知系数,即材料的未知参数。     

Properties of surface waves in an elastic cylinder filled with a liquid

The properties of harmonic surface waves in an elastic cylinder filled with a liquid are studied. The case of elastic material for which the shear wave velocity is higher than the sound velocity in a liquid is considered. The wave motion is described based on the complete system of equations of the dynamic theory of elasticity and the equation of motion of an ideal compressible liquid. The asymptotic analysis of the dispersion equation in the region of large wave numbers and qualitative analysis of the dispersion spectrum showed that in such a waveguiding system there exist two surface waves, the Stoneley and the Rayleigh waves. The lowest normal wave forms the Stoneley wave on the internal surface of the cylinder. In this waveguide phase, velocities of all normal waves, except for the lowest one, have the velocity of sound in the liquid as their limit. Therefore, the Rayleigh wave on the external surface of the cylinder is formed by all normal waves in the range of frequencies and wave numbers in which phase velocities of normal waves of the composite waveguide and the lowest normal wave of the elastic hollow cylinder coincide.     

Cumulative solutions of nonlinear longitudinal vibration in isotropic solid bars

Based on the strain invariant relationship and taking the high-order elastic energy into account, a nonlinear wave equation is derived, in which the excitation, linear damping, and the other nonlinear terms are regarded as the first-order correction to the linear wave equation. To solve the equation, the biggest challenge is that the secular terms exist not only in the fundamental wave equation but also in the harmonic wave equation （unlike the Duffing oscillator, where they exist only in the fundamental wave equation）. In order to overcome this difficulty and to obtain a steady periodic solution by the perturbation technique, the following procedures are taken： （i） for the fundamental wave equation, the secular term is eliminated and therefore a frequency response equation is obtained; （ii） for the harmonics, the cumulative solutions are sought by the Lagrange variation parameter method. It is shown by the results obtained that the second- and higher-order harmonic waves exist in a vibrating bar, of which the amplitude increases linearly with the distance from the source when its length is much more than the wavelength; the shift of the resonant peak and the amplitudes of the harmonic waves depend closely on nonlinear coefficients; there are similarities to a certain extent among the amplitudes of the odd- （or even-） order harmonics, based on which the nonlinear coefficients can be determined by varying the strain and measuring the amplitudes of the harmonic waves in different locations.     

A numerical study on the propagation of Rayleigh and guided waves in cortical bone according to Mindlin's Form II gradient elastic theory

Cortical bone is a multiscale heterogeneous natural material characterized by microstructural effects. Thus guided waves propagating in cortical bone undergo dispersion due to both material microstructure and bone geometry. However, above 0.8 MHz, ultrasound propagates rather as a dispersive surface Rayleigh wave than a dispersive guided wave because at those frequencies, the corresponding wavelengths are smaller than the thickness of cortical bone. Classical elasticity, although it has been largely used for wave propagation modeling in bones, is not able to support dispersion in bulk and Rayleigh waves. This is possible with the use of Mindlin's Form-II gradient elastic theory, which introduces in its equation of motion intrinsic parameters that correlate microstructure with the macrostructure. In this work, the boundary element method in conjunction with the reassigned smoothed pseudo Wigner-Ville transform are employed for the numerical determination of time-frequency diagrams corresponding to the dispersion curves of Rayleigh and guided waves propagating in a cortical bone. A composite material model for the determination of the internal length scale parameters imposed by Mindlin's elastic theory is exploited. The obtained results demonstrate the dispersive nature of Rayleigh wave propagating along the complex structure of bone as well as how microstructure affects guided waves.     

一维非线性声波传播特性

针对一维非线性声波的传播问题进行了有限元仿真和实验研究. 首先推导了一维非线性声波方程的有限元形式, 含有高阶矩阵的非线性项导致声波具有波形畸变、谐波滋生、基频信号能量向高次谐波传递等非线性特性. 编制有限元程序对一维非线性声波进行了计算并对仿真得到的畸变非线性声波信号进行处理, 分析其传播性质和物理意义. 为验证有限元计算结果, 开展了水中的非线性声波传播的实验研究, 得到了不同输入信号幅度激励下和不同传播距离的畸变非线性声波信号. 然后对基波和二次谐波的传播性质进行详细讨论, 分析了二次谐波幅度与传播距离和输入信号幅度的变化关系及其意义, 拟合出二次谐波幅度随传播距离变化的方程并阐述了拟合方程的物理意义. 结果表明, 数值仿真信号及其频谱均与实验结果有较好的一致性, 证实计算方法和结果的正确性, 并提出了具有一定物理意义的二次谐波随传播距离变化的简单数学关系. 最后还对固体中的非线性声波传播性质进行了初步探讨. 本研究工作可为流体介质中的非线性声传播问题提供理论和实验依据.     

The problem of surface wave propagation in a reduced Cosserat medium

This article continues a series of publications devoted to the study of waves in the framework of the asymmetric theory of elasticity, where the deformed state of the medium is characterized by independent vectors of translation and rotation. The problem of acoustic Rayleigh wave propagation in half space is considered within a model of the reduced Cosserat medium. A general analytic solution of this problem is obtained. The analysis of this solution is compared with the corresponding solution for a classical elastic medium and full linear Cosserat medium. It is shown that the Rayleigh wave is characterized by a range of forbidden frequencies, where this wave cannot propagate. The dispersion curve consists of two branches. One of them has a cut-off frequency and cut-off wavenumber.     

Nonlinear,nondispersive wave equations: Lagrangian and Hamiltonian functions in the hodograph transformation

The hodograph transformation is generally used in order to associate a system of linear partial differential equations to a system of nonlinear (quasilinear) differential equations by interchanging dependent and independent variables. Here we consider the case when the nonlinear differential system can be derived from a Lagrangian density and revisit the hodograph transformation within the formalism of the Lagrangian-Hamiltonian continuous dynamical systems.Restricting to the case of nondissipative, nondispersive one-dimensional waves, we show that the hodograph transformation leads to a linear partial differential equation for an unknown function that plays the role of the Lagrangian in the hodograph variables. We then define the corresponding hodograph Hamiltonian and show that it turns out to coincide with the wave amplitude. i.e., with the unknown function of the independent variables to be solved for in the initial nonlinear wave equation.     

微结构固体中的孤立波及其存在条件

根据Mindlin理论,考虑宏观尺度非线性效应、二次和三次微尺度非线性效应以及微尺度频散效应,建立了描述一维微结构固体中纵波传播的一种新模型.用动力系统定性分析理论,分析了微结构固体中孤立波的存在条件及其几何特征,证明了当介质参数和孤立波传播速度满足适当条件时,在二次微尺度非线性效应的影响下微结构固体中可以形成一种非对称孤立波,在三次微尺度非线性效应的影响下微结构固体中可以形成一种对称孤立波.最后,用数值方法进一步验证了上述结论.     

Flexural edge waves in semi-infinite elastic plates

This paper presents a detailed analysis of the dispersion for flexural edge waves in semi-infinite isotropic elastic plates. A solution to the dynamic equations of motion is constructed by the superposition of two partial solutions, each providing zero shear stresses at the plate faces. A dispersion equation is expressed via the determinant of an infinite system of linear algebraic equations. The system is reduced to a finite one by taking into account the asymptotic behaviour of unknown coefficients. The accuracy of the solution is confirmed by a good agreement with the available experimental data and by a proper satisfaction of the prescribed boundary conditions.A detailed analysis of dispersion properties for the edge wave and corresponding displacements at various frequencies is carried out. In addition to the well-known results it is shown that the plate height does not influence the existence of the edge wave at high frequencies and, as the frequency increases, the phase velocity of the edge wave in a semi-infinite plate asymptotically approaches the velocity of an edge wave in a right-angled wedge. The performed analysis allows evaluating the plate theories such as the Kirchhoff theory or other refined plate theories developed for modeling edge waves in semi-infinite elastic plates at low frequencies.     

Generation of Elastic Waves by a Magnetic Field Associated with Reversible Rotations in Triaxial Magnets

In the first-order anharmonicity approximation of Hooke's law, the amplitudes of elastic waves, absorption coefficients, wave numbers of the fundamental wave and of the first harmonic generated by an alternating magnetic field with a preset orientation relative to the basis axes of a crystal having arbitrary dimensions are calculated for multidomain magnets with rigidly fixed domain boundaries in terms of concentrations of magnet phases and magnetic structure parameters with allowance for the wave equation and angular momenta.     

Investigation of Lamb elastic waves in anisotropic multilayered composites applying the Green’s matrix

This article presents a numerical study of dispersion characteristics of some symmetric and antisymmetric composites modelled as multilayered packets of layers with arbitrary anisotropy of each layer. The authors introduce a subsidiary boundary problem of three-dimensional elasticity theory for the system of partial differential equations describing the harmonic oscillations of the composite caused by a surface load. The problem reduces to a boundary problem for ordinary differential equations by employing the Fourier transform. An algorithm of constructing the Fourier transform of the Green’s matrix of the given boundary problem is presented. The wave numbers of Lamb waves propagating in composites, their phase velocity surfaces and group wave surfaces are presented through the poles of the transform of the Green’s matrix. The authors obtain the dispersion curves for different directions and frequencies and investigate the dispersion curves and surfaces of wave numbers, phase velocities and group wave surfaces for various composites. The numerical results are then compared with the results obtained by applying other methods.     

VIBRATION EIGENFREQUENCY ANALYSIS OF A SINGLE-LINK FLEXIBLE MANIPULATOR

We analyse the vibration eigenfrequencies of a flexible slewing beam with a payload attached at one end. A wave propagation method (WPM) is used. There are four types of waves which propagate along a beam—two dispersive waves travelling in opposite directions, and two evanescent waves near the endpoints. We add a fifth time-harmonic function corresponding to oscillation of the beam at the payload end. We show that the large frequencies are asymptotically identical to those for the clamped–free beam, independent of the payload. For small eigenfrequencies, we incorporate WPM with a perturbation iteration procedure, the results of which agree well with “exact” values which result from solving a transcendental equation cited elsewhere in the literature.     

粘性物质中正激波稳定性分析

用线性稳定性理论,分析粘性物质中的正激波稳定性问题.粘性物质中任意强度的一维激波,其稳定性问题可归结为处理复数范围内的特征值问题,该特征值问题由两个一阶常微分方程及一个二阶常微分方程构成.这些常微分方程的系数依赖于流动的基本流场的物理量及其梯度.所获得的特征值问题由一个四阶精度的有限差分离散求解.分析考虑物质粘性的金属铝中的正激波稳定性,可以看出,正激波运动是稳定的,并且激波速度对波前和波后的小扰动量的衰减有相反的作用,而物质粘性有致稳的作用.     

Nonlinear gravitational waves on the surface of viscous fluid: Lagrangian approach

The problems of the asymptotic theory of weakly nonlinear surface waves in viscous fluid are discussed. For standing waves on deep water, the solutions obtained in the first- and second-order approximations in a small parameter—wave steepness—are analyzed. The evolution equation for the amplitude of wave packet envelope is obtained where the inverse Reynolds number is equal to the squared steepness. It is shown that this is a nonlinear Schrödinger equation with linear dissipation.     

Leaky surface-elastic-waves on flat and highly corrugated surfaces for isotropic non-dissipative media

The dispersion relation for Rayleigh waves on a grating-surface of a semi-infinite, isotropic, nondissipative, elastic medium is solved numerically, with complex wave vector k or complex frequency, in the radiative region (above the bulk transverse sound-line) and within the first frequency-gap on the Brillouin zone boundary created by the grating periodicity. The acoustic attenuation, found from the imaginary part of k, agrees well with experiment. A dispersive branch, with corresponding complex-solutions for the flat surface, between the bulk transverse and longitudinal sound-lines, representing a new leaky mode or surface resonance, accounts for the principal peak in the attenuation.     

Detection of surface acoustic waves by scanning tunneling microscopy

A new method for the investigation of ultrasonic waves on surfaces of solids based on scanning tunneling microscopy is presented. A sinusoidal high frequency signal is added to the tip voltage. Hence the tunneling current contains a component whose frequency is the difference of the frequencies of the acoustic wave field and the ac tip voltage. Amplitude and phase of this component carry the full information about the wave field.     

Measuring of viscoelastic properties of homogeneous soft solid using transient elastography: an inverse problem approach

Two main questions are at the center of this paper. The first one concerns the choice of a rheological model in the frequency range of transient elastography, sonoelasticity or NMR elastography for soft solids (20-1000 Hz). Transient elastography experiments based on plane shear waves that propagate in an Agar-gelatin phantom or in bovine muscles enable one to quantify their viscoelastic properties. The comparison of these experimental results to the prediction of the two simplest rheological models indicate clearly that Voigt's model is the better. The second question studied in the paper deals with the feasibility of quantitative viscosity mapping using inverse problem algorithm. In the ideal situation where plane shear waves propagate in a sample, a simple inverse problem based on the Helmholtz equation correctly retrieves both elasticity and viscosity. In a more realistic situation with nonplane shear waves, this simple approach fails. Nevertheless, it is shown that quantitative viscosity mapping is still possible if one uses an appropriate inverse problem that fully takes into account diffraction in solids.     

Stability of a Composite Wave of Two Viscous Shock Waves for the Full Compressible Navier-Stokes Equation

In this paper we investigate the asymptotic stability of a composite wave consisting of two viscous shock waves for the full compressible Navier-Stokes equation. By introducing a new linear diffusion wave special to this case, we successfully prove that if the strengths of the viscous shock waves are suitably small with same order and also the initial perturbations which are not necessarily of zero integral are suitably small, the unique global solution in time to the full compressible Navier-Stokes equation exists and asymptotically tends toward the corresponding composite wave whose shifts (in space) of two viscous shock waves are uniquely determined by the initial perturbations. We then apply the idea to study a half space problem for the full compressible Navier-Stokes equation and obtain a similar result. Research is supported in part by NSFC Grant No. 10471138, NSFC-NSAF Grant No. 10676037 and 973 project of China, Grant No. 2006CB805902, in part by Japan Society for the Promotion of Science, the Invitation Fellowship for Research in Japan (Short-Term). Research is supported in part by Grant-in-Aid for Scientific Research (B) 19340037, Japan.