高频纵向导波在钢杆中传播特性的研究

计算了钢杆中纵向轴对称导波模态的衰减频散曲线和群 速度频散曲线. 分析了1~3MHz范围内高频纵向轴对称超声导波在钢杆中的传播特 性. 理论分析表明，各高阶纵向模态都存在一个衰减最小值，在此衰减最小值所对应频率下 的高阶模态能传播较远距离，可用于钢杆导波检测. 建立实验系统，采用轴对称同端激 励接收的方法，根据第1次端面回波做出群速度和端面回波幅值随频率变化曲线，实验结 果与理论分析基本吻合. 表明考虑材料衰减特性的钢杆频散曲线可以作为实验指导依据.     

分数阶粘弹性地基上薄板波动和振动特性分析

本文研究了粘弹性地基上薄板的波动和振动问题．主要讨论了基于分数导数理论的粘弹性地基模型上 薄板弯曲波的传播特性以及固有频率对地基的依赖特性．推导了三种经典粘弹性地基模型的复模量．并利用分 数导数的性质得到分数阶粘弹性地基上 Kirchhoff板中弯曲波的传播速度、衰减系数以及自由振动的复固有频 率．数值算例表明粘弹性地基对弯曲波传播特性存在显著影响,不同粘弹性模型所对应的色散和衰减特性也存 在较大差别．分数阶导数可以实现相邻整数阶导数之间的光滑过渡．利用分数导数的本构关系可以更加真实地 描述粘弹性地基的历史依赖行为,更准确地表现出粘弹性地基板中弯曲波的色散和衰减特性．     

Inhomogeneous wave propagation in partially saturated soils

In the present study, inhomogeneous plane harmonic waves propagation in dissipative partially saturated soils are investigated. The analytical model for the dissipative partially saturated soils is solved in terms of Christoffel equations. These Christoffel equations yields the existence of four wave (three longitudinal and one shear) modes in partially saturated soils. Christoffel equations are further solved to determine the complex velocities and polarizations of four wave modes. Inhomogeneous propagation is considered through a particular specification of complex slowness vector. A finite non-dimensional inhomogeneity parameter is considered to represent the inhomogeneous nature of these four waves. Impact of tortuosity parameter on the movement of pore fluids is considered. Hence, the considered model is capable of describing the wave behavior at high as well as mid and low frequencies. Numerical example is considered to study the effects of inhomogeneity parameter, saturation of water, porosity, permeability, viscosity of fluid phase and wave frequency on the velocity and attenuation of four waves. It is observed that all the waves are dispersive in nature (i.e., frequency dependent).     

On dispersive propagation of surface waves in patchy saturated porous media

Frequency-dependent velocity and attenuation for Rayleigh-wave propagation along a vacuum/patchy saturated porous medium interface are investigated in the low frequency band (0.1–1000 Hz). Conventional patchy saturation models for compressional waves are extended to account for Rayleigh wave propagation along a free surface. The mesoscopic interaction of fluid and solid phases, as a dominant loss mechanism in patchy saturated media, significantly affects Rayleigh-wave propagation and attenuation. Researches on the dispersion characteristics at low frequencies with different gas fractions in patchy saturated media also demonstrate a strong correlation between the Rayleigh-wave mode and the fast compressional wave. Especially, the strongest attenuation with the maximum value of 1/Q$1/Q$ for Rayleigh waves are obtained in the frequency range of 1–200 Hz. Numerical results show that the significant dependence of velocity and attenuation on frequencies and gas fractions presents a distinctive dynamical response of Rayleigh waves in the time domain.     

Wave propagation in transversely isotropic porous piezoelectric materials

Wave propagation in porous piezoelectric material (PPM), having crystal symmetry 6 mm, is studied analytically. Christoffel equation is derived for the propagation of plane harmonic waves in such a medium. The roots of this equation give four complex wave velocities which can propagate in such materials. The phase velocities of propagation and the attenuation quality factors of all these waves are described in terms of complex wave velocities. Phase velocities and attenuation of the waves in PPM depend on the phase direction. Numerical results are computed for the PPM BaTiO₃. The variation of phase velocity and attenuation quality factor with phase direction, porosity and the wave frequency is studied. The effects of anisotropy and piezoelectric coupling are also studied. The phase velocities of two quasi dilatational waves and one quasi shear waves get affected due to piezoelectric coupling while that of type 2 quasi shear wave remain unaffected. The phase velocities of all the four waves show non-dispersive behavior after certain critical high frequency. The phase velocity of all waves decreases with porosity while attenuation of respective waves increases with porosity of the medium. The characteristic curves, including slowness curves, velocity curves, and the attenuation curves, are also studied in this paper.     

A mixture theory analysis for the surface-wave propagation in an unsaturated porous medium

The mixture theory is employed to the analysis of surface-wave propagation in a porous medium saturated by two compressible and viscous fluids (liquid and gas). A linear isothermal dynamic model is implemented which takes into account the interaction between the pore fluids and the solid phase of the porous material through viscous dissipation. In such unsaturated cases, the dispersion equations of Rayleigh and Love waves are derived respectively. Two situations for the Love waves are discussed in detail: (a) an elastic layer lying over an unsaturated porous half-space and (b) an unsaturated porous layer lying over an elastic half-space. The wave analysis indicates that, to the three compressional waves discovered in the unsaturated porous medium, there also correspond three Rayleigh wave modes (R1, R2, and R3 waves) propagating along its free surface. The numerical results demonstrate a significant dependence of wave velocities and attenuation coefficients of the Rayleigh and Love waves on the saturation degree, excitation frequency and intrinsic permeability. The cut-off frequency of the high order mode of Love waves is also found to be dependent on the saturation degree.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

Flutter and resonances of a flag near a free surface

We investigate the effects of a nearby free surface on the stability of a flexible plate in axial flow. Confinement by rigid boundaries is known to affect flag flutter thresholds and fluttering dynamics significantly, and this work considers the effects of a more general confinement involving a deformable free surface. To this end, a local linear stability is proposed for a flag in axial uniform flow and parallel to a free surface, using one-dimensional beam and potential flow models to revisit this classical fluid–structure interaction problem. The physical behaviour of the confining free surface is characterized by the Froude number, corresponding to the ratio of the incoming flow velocity to that of the gravity waves. After presenting the simplified limit of infinite span (i.e. two-dimensional problem), the results are generalized to include finite-span and lateral confinement effects. In both cases, three unstable regimes are identified for varying Froude number. Rigidly-confined flutter is observed for low Froude number, i.e. when the free surface behaves as a rigid wall, and is equivalent to the classical problem of the confined flag. When the flow and wave velocities are comparable, a new instability is observed before the onset of flutter (i.e. at lower reduced flow speed) and results from the resonance of a structural bending wave and one of the fundamental modes of surface gravity waves. Finally, for large Froude number (low effect of gravity), flutter is observed with significant but passive deformation of the free surface in response of the flag’s displacement.     

An explicit formulation for the evolution of nonlinear surface waves interacting with a submerged body

An explicit formulation to study nonlinear waves interacting with a submerged body in an ideal fluid of infinite depth is presented. The formulation allows one to decompose the nonlinear wave–body interaction problem into body and free‐surface problems. After the decomposition, the body problem satisfies a modified body boundary condition in an unbounded fluid domain, while the free‐surface problem satisfies modified nonlinear free‐surface boundary conditions. It is then shown that the nonlinear free‐surface problem can be further reduced to a closed system of two nonlinear evolution equations expanded in infinite series for the free‐surface elevation and the velocity potential at the free surface. For numerical experiments, the body problem is solved using a distribution of singularities along the body surface and the system of evolution equations, truncated at third order in wave steepness, is then solved using a pseudo‐spectral method based on the fast Fourier transform. A circular cylinder translating steadily near the free surface is considered and it is found that our numerical solutions show excellent agreement with the fully nonlinear solution using a boundary integral method. We further validate our solutions for a submerged circular cylinder oscillating vertically or fixed under incoming nonlinear waves with other analytical and numerical results. Copyright © 2007 John Wiley & Sons, Ltd.     

变深度浅水域中非定常船波

以Green—Naghdi(G—N)方程为基础，采用波动方程／有限元法计算船舶经过变深度浅水域时非定常波浪特性．把运动船舶对水面的扰动作为移动压强直接加在Green-Naghdi方程里，以描述运动船体和水面的相互作用．以Series60 CB＝0．6船为算例，给出自由面坡高，波浪阻力在船舶经过一个水下凸包时变化规律，并与浅水方程的结果进行了比较．计算结果表明，当船舶经过凸包时，波浪阻力先增加，后减少，并逐渐趋于正常．同时发现，当船速小于临界速度时(Fr＝√gh＜1．0)，G—N方程给出的船后尾波比浅水方程的结果明显，波浪阻力也比浅水方程的结果有所提高，频率散射必须考虑．当船速大于临界速度时(Fr＝√gh＞1．0)，G—N方程的计算结果与浅水方程差别不大，频率散射的影响可以忽略．     

Oscillatory line source for water waves in shear flow

The linearized water-wave radiation problem for the oscillating 2D submerged source in an inviscid shear flow with a free surface is investigated analytically. The vorticity is uniform, with zero velocity at the free surface. Then there will be at most two emitted waves, and no Doppler effects. Exact far-field waves are derived, with radiation conditions applied at infinity. An upstream wave will always exist, whereas the downstream wave exists only when the angular frequency of oscillation exceeds the vorticity. The wave radiation problem is solved also for oscillating vortex and dipoles. The amplitudes and energy fluxes are calculated.     

缺陷岩体中弹性波传播的多次散射效应分析

弹性波在岩体中传播时与岩体缺陷相互作用形成复杂的传播图案。为研究缺陷对弹性波多次散射作用的影响,建立了双椭圆缺陷模型,基于Green函数基本解,采用边界积分的计算方法,得到了反映缺陷界面条件的刚度矩阵,分析了弹性波在双椭圆缺陷间的多次散射效应。结果表明：与单椭圆缺陷模型相比,双缺陷的相互作用使得弹性波频散和衰减效应增强,定量给出了缺陷的影响区域,从而明确了多次散射效应的尺度界限。进一步探讨了弹性波传播的多尺度效应,结果表明频散的Rayleigh峰、Mie峰和衰减的峰值频率同椭圆长轴和入射波波长两个尺度密切相关,存在明确的定量关系。相应的数值模拟结果表明,弹性波和缺陷相互作用在缺陷界面上诱发界面波,该界面波也存在频率相关性,影响了弹性波宏观传播的频散和衰减特征。     

各向异性粘弹性多孔截止中平面波的传播

This study discusses wave propagation in perhaps the most general model of a poroelastic medium. The medium is considered as a viscoelastic, anisotropic and porous solid frame such that its pores of anisotropic permeability are filled with a viscous fluid. The anisotropy considered is of general type, and the attenuating waves in the medium are treated as the inhomogeneous waves. The complex slowness vector is resolved to define the phase velocity, homogeneous attenuation, inhomogeneous attenuation, and angle of attenuation for each of the four attenuating waves in the medium. A non-dimensional parameter measures the deviation of an inhomogeneous wave from its homogeneous version. An numerical model of a North-Sea sandstone is used to analyze the effects of the propagation direction, inhomogeneity parameter, frequency regime, anisotropy symmetry, anelasticity of the frame, and viscosity of the pore-fluid on the propagation characteristics of waves in such a medium.     

Two dimensional wave problems in rotating elastic media

The rotation of an elastic medium makes it act anisotropically and dispersively. The eigenvectors for plane wave propagation are in general complex and thus the waves are elliptically polarized. In general the waves are neither pure shear nor pure compressional waves, and their speeds depend on the ratio of rotational frequency of the medium and the angular frequency of the wave.The class of problems discussed here involves waves propagating perpendicularly to the axis of rotation and in particular we discuss plane strain modes. The reflection and refraction of plane waves is considered.The plane waves are used to construct a general solution in cylindrical coordinates. The solution is given in terms of Bessel functions. The cylindrical solution is applied to scattering by circular cylinders. The problem of free oscillations is mentioned briefly.     

A linear dynamic model for a saturated porous medium

A linear isothermal dynamic model for a porous medium saturated by a Newtonian fluid is developed in the paper. In contrast to the mixture theory, the assumption of phase separation is avoided by introducing a single constitutive energy function for the porous medium. An important advantage of the proposed model is it can account for the couplings between the solid skeleton and the pore fluid. The mass and momentum balance equations are obtained according to the generalized mixture theory. Constitutive relations for the stress, the pore pressure are derived from the total free energy accounting for inter-phase interaction. In order to describe the momentum interaction between the fluid and the solid, a frequency independent Biot-type drag force model is introduced. A temporal variable porosity model with relaxation accounting for additional attenuation is introduced for the first time. The details of parameter estimation are discussed in the paper. It is demonstrated that all the material parameters in our model can be estimated from directly measurable phenomenological parameters. In terms of the equations of motion in the frequency domain, the wave velocities and the attenuations for the two P waves and one S wave are calculated. The influences of the porosity relaxation coefficient on the velocities and attenuation coefficients of the three waves of the porous medium are discussed in a numerical example.     

Turbulent mass transport and attenuation in Stokes waves

The motion of turbulent Stokes waves on a finite constant depth fluid with a rough bed is considered. First and second order turbulent boundary layer equations are solved numerically for a range of roughness parameters, and from the solutions are calculated the mass transport velocity profiles and attenuation coefficients. A new mechanism of turbulent mass transport is found which predicts a reduction and reversal of drift velocity in shallow water in agreement with experimental observations under turbulent conditions. This transpires because the second order Stokes wave motion, in a turbulent boundary layer, can directly influence the mass transport velocity by mode coupling interactions between different second order Fourier modes of oscillation. It is also found that the Euler contribution due to the radiation stress of the first order motion is reduced to half of it's corresponding laminar value as a consequence of the velocity squared stress law. The attenuation is found to be of inverse algebraic type with the reciprocal wave height varying linearly with either distance or time. The severe wave height restriction applicable to the Longuet-Higgins [4] solution is shown not to apply to progressive waves on a finite constant depth of fluid. The existence of sand bars on sloping beaches exposed to turbulent waves is predicted.     

流体饱和标准线性粘弹性多孔介质中的平面波

研究了流体饱和不可压标准线性粘弹性多孔介质中平面波的传播和反射问题．在固相骨架小变形的假定下,得到了粘弹性多孔介质中波动方程的一般解,讨论了弥散关系和波的衰减特性．结果表明：在流体饱和不可压粘弹性多孔介质中,仅存在一个耦合纵波和一个耦合横波,纵波和横波的波速、衰减率等取决于孔隙流体与固相骨架间的相互作用以及固相骨架本身的粘性．同时,研究了半空间自由边界上入射波(纵波、横波)的反射问题。得到了非均匀反射波的波速、反射系数、衰减率等的表达式及其相关的数值结果．     

基于时间反转法的Lamb波检测技术的研究进展

Lamb波具有传播距离远、衰减小等特点,已被广泛应用于大型板类结构的损伤检测。时间反转是实现超声波聚焦的有效方法之一。本文回顾并总结了近20年来时间反转方法在Lamb波检测领域中的研究进展。针对Lamb波固有的频散和多模态等传播特性,详细论述了时间反转法在板结构Lamb波检测中的应用。最后对基于时间反转法的免基准信号Lamb波损伤检测方法进行了总结,分析其存在的问题并对其应用前景进行了展望。     

Investigation of the far and near fields in the problem of exponentially stratified flow over a bottom irregularity

The three-dimensional problem of finite-depth stratified flow over a small bottom irregularity is considered in mixed Euler-Lagrange variables. The Brunt-Väisälä frequency is assumed to be constant and small, and the free surface condition is replaced by the rigid roof condition. Investigation of the far field showed that the principal wave perturbations lie within an angle which for large values of the internal Froude number is much less than theKelvin angle, while the wave amplitude at infinity is of the order of l/r, where r is the polar radius. The ring perturbations are exponentially damped. As distinct from point source models, the model in question does not lead to divergence of the integrals on the flow axis [1-3]. Appproximate expressions for the radial and ring waves in terms of certain universai functions were obtained for investigating the near and far fields when the bottom irregularity is hemispherical. For the radial waves a law of similarity was obtained for which the characteristic dimension in the direction of the flow axis is the ratio of the flow velocity to the Brunt-Väisälä frequency, and the characteristic dimension in a direction perpendicular to the flow axis the depth of the fluid. In the first approximation the ring perturbations do not depend on the Brunt-Väisälä frequency. It is shown that in the near field the zone of intense wave perturbations is of the order of the fluid depth and not of the dimensions of the obstacle as for Kelvin ship waves on the surface of a homogeneous fluid.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 86–94, September–October, 1987.