岩体体积节理数(Jv)的现场测量方法评价

采用横观各向同性的流体饱和孔隙介质地层模型, 模拟了井中由点声源激发的声场. 计算出了井孔声场函数的全部极点, 发现最低阶泄漏模式波极点的虚部与正常模式波极点的虚部大小属于同一数量级, 因而是不可忽略的. 在算例中还发现, 伪瑞利波相速度的低频极限值不一定等于地 层横波速度. 计算了由极点留数对应的斯通利波、伪瑞利波和两种泄漏模式波的激 发频谱和时域波形, 以及由割线积分给出的激发频谱和时域波形. 通过与全波的比较, 发现在一些情况下仅由准横波割线积分不能获得准横波首波, 必须将泄漏模式波与准横波割线积分波叠加后才能得到准横波首波, 泄漏模式对全波的贡献不能忽略.     

白云岩受压声学特性及其在裂缝研究中的应用

通过实验手段研究了白云岩在三轴和单轴加载过程中的声波波速、幅度、频谱特性的变化特征. 结果发现:(1)白云岩在受压变形过程中,声波波速的变化较好地反映了岩石内部裂缝的闭合、产生、扩展和贯穿等变化. 总体而言,横波波速比纵波能更好地预测裂缝的产生,而纵波对裂缝的非稳定发展或贯穿有较敏感的反映. (2)随着轴压的增大,纵、横波波形的末端都有散射波的出现,且横波波形末端的"鱼尾"状散射波信号更为明显(此时应力点为极限强度的60%左右),预示着岩石内部裂缝的产生和稳定扩展. (3)频谱曲线也很好地反映了岩石内部结构变形情况. 随着岩石受力的增大,频谱曲线上的振幅都呈现增大趋势,标志着岩石的压实阶段;当频谱曲线上低频段较高频段活跃时,标志着裂缝的产生;更有甚者,低频还会取代高频成为主频. (4)在岩石受压裂缝闭合阶段,首波振幅和频谱主振幅都呈现上升趋势;而在裂缝不稳定扩展阶段,主振幅表现出比首波振幅低的上升趋势;振幅曲线达到峰值后都存在一突降拐点,预示着岩石的贯通破坏. 研究对于地层岩石的动态长期监测和工程岩体的稳定性评价都具有重要的理论参考价值.      

孔隙地层中声波诱导的电磁场测井响应模拟

声电效应测井过程中,当充流体井中发射的声波传播至井外时,在孔隙地层中引起含带电离 子的孔隙流体相对于固体骨架流动,导致电磁场. 在忽略转换电磁场对声场影响的条件下,先采用Biot-Rosenbaum模型计算孔隙地层中的 声场,再计算声波诱导的电场和磁场. 计算结果表明,有伴随纵波波群和横波波群的电磁场、有伴随斯通利波的电磁场,还有微弱的 临界折射电磁波. 计算还表明: 接收器偏离井轴越远,电场和磁场的强度越大; 在井壁处,磁场强度的波形与轴向电场强度的波形相似.     

关于匀质弹性半空间上刚体垂直振动的轴对称混合边值问题

本文继,Awojobi等之后对匀质弹性半空间上刚体垂直振动轴对称混合边值问题进行了新的研究。抛弃由Reissner所建议的用消除驻波确定自由波项系数的方法,直接由弹性力学方程与有关问题的混合边界条件处理了自由波项,使著名的瑞利极点化为可去奇点,从而建立了与过去不同的表征有关边值问题的第二类弗里德荷姆积分方程,其积分核是非奇异的并在寻常积分定义下有意义。对所获得的方程进行了详细的解答,并与大型动力基础实测资料进行了对比,从而作出了半空间内介质的阻尼损耗不可忽略的结论。     

区域脉冲载荷下二维Lamb问题的精确求解

采用积分变换方法,并利用两类积分公式克服反变换求解的困难,求得了区域脉冲载荷下一个二维Lamb问题的代数形式的精确解.基于该分段函数形式的代数结果,纵波、横波、Rayleigh波等应力波成分在弹性表面的激发和传播过程得到详尽分析,其中很多结论是已有的解析积分结果或者数值计算结果不曾得到的.     

线状脉冲激光激发双层圆柱的波传播

该文建立了一个物理模型来预测均匀、各向同性双层固体圆柱在激光烧蚀或热弹激发下的弹性波传播规律.在建模中忽略了激光源的光学穿透和材料的热扩散.采用空间和时间上的双重傅氏变换对理论模型进行求解,得到其位移场的变换解.通过对铝壳铜芯和铜壳铝芯的数值计算,得到不同内外径比下瞬态位移响应的时域和频域数值波形.用传播路径分析发现结果中存在界面反射波,这可以表明所建立的理沦模型和计算方法是有效的,同时其余波的传播也被分析和解释.     

铁木辛柯梁中的卸载弯曲波及二次断裂

半无限长梁承受恒定弯矩作用后, 如果自由端的初始弯矩突然释放, 将在梁中激发出一列卸载弯曲应力波. 采用铁木辛柯梁和瑞利梁来研究突然卸载所激发出的弯曲波的传播特征. 利用拉普拉斯变换方法进行分析, 首先推导出铁木辛柯梁和瑞利梁中的卸载弯曲波的像函数解析解, 采用数值反变换方法给出了时域上波传播的响应解, 并研究了梁中各点的横向位移、弯矩和剪力随时间的变化规律. 计算结果表明: 与简化的欧拉梁不同, 旋转惯性的引入使铁木辛柯梁和瑞利梁中的弯曲波传播具有强烈的局部化效应, 特别是梁中各点经历的弯矩变化, 和其距离自由端的位置相关, 不同时刻的弯矩峰值大小不同;瑞利梁中离自由端不同距离各点的峰值弯矩先增大后降低, 最后达到一个渐近值;铁木辛柯梁中各点的峰值弯矩总体上随着时间单调增大到同一个渐近值, 该渐近值与欧拉梁中的峰值弯矩值相同, 均为1.43.切应力效应的引入进一步降低了铁木辛柯梁中卸载弯曲波的波速, 同时也使得铁木辛柯梁中弯矩峰值的最大值小于瑞利梁中的最大值. 对于脆性细长梁的纯弯曲断裂, 铁木辛柯梁可以较好地预测二次断裂的发生位置, 相应的碎片尺寸约为7倍梁横截面厚度.      

激光热弹激励流-固泄漏Lamb波的有限元模拟

根据变分原理,得到热弹体运动方程和热传导方程相对应的有限元方程. 通过数值积分方法求解有限元方程,得到脉冲激光线源在水/铝、空气/铝这两种流-固界面上热弹激励的泄漏Lamb波瞬态波形. 计算结果表明,泄漏Lamb波不但存在于液-固界面,而且存在于气-固界面;和Lamb波相反,泄漏Lamb波的S_0模态是反对称的,而A_0模态是对称的;但由于这两种流-固界面的性质差异导致泄漏Lamb波的波形和幅度不同.      

模式自适应连续小波去除趋势项方法在爆破振动信号分析中的应用

为了更精确提取爆破振动信号峰值速度、能量等重要特征,必须对爆破振动加速度信号时域积分中的趋势项予以去除。通过对实测爆破振动加速度信号进行梯形数值积分,提出以时域积分后的爆破振动速度信号来构造模式自适应小波基的方法,并用此方法去除时域积分后爆破振动速度信号中的趋势项,然后对去除趋势项后的爆破振动速度信号进行能量特征分析。结果表明:模式自适应连续小波法成功去除了时域积分后爆破振动速度信号中的趋势项;与建立在传统Fourier变换基础上的频谱分析相比,小波变换的能量分析具有更精细的频率分辨率,更适合于对频率分辨率要求更高的爆破振动信号进行分析;各频率区间范围划分越宽,爆破振动加速度信号与速度信号各频率区间内能量分布的相关程度越高,反之,相关程度越低。     

滚波演化中聚合过程的数值模拟研究

滚波是一种重力作用下自由液面失稳诱发的水面波动现象, 通常可分为具有相对稳定波形和波速的周期性滚波与波形和波速不断变化的不规则滚波(自然滚波). 不规则滚波的相互作用和发展演化过程十分复杂, 至今对其认识尚不成熟. 本文采用基于雷诺平均Navier-Stokes方程的立面二维数值模型, 对不规则滚波发展过程中的吸收聚合和追赶聚合现象进行了数值模拟研究. 分析了两种聚合模式的演化过程, 给出了滚波聚合过程中完整的波形、波速、速度剖面以及湍流黏性等重要信息. 结果表明滚波的聚合过程是不规则滚波演化和增长的重要机制, 在特定条件下滚波增长由自然增长模式转变为以吸收聚合和追赶聚合为主的增长模式. 滚波聚合过程中, 依次经历后波追赶、爬升、与前波合并、内部流场调制等多个步骤, 最终形成一个具有更大波长和波高的滚波. 本文发现了在3个滚波间距较近的情况下, 会发生二重聚合现象, 即后两个滚波首先聚合, 然后与前波进一步聚合形成一个新的滚波.      

Comments on the relation between surface wave theory and the theory of reflection

The sextic Stroh formalism, previously extensively used in the analysis of subsonic phenomena, has been used for the analysis of reflection phenomena and leaky surface waves in the first transsonic range of velocities. In particular the behaviour of the reflection problem at the limiting velocity is studied. It is shown that when the condition of free surface can be satisfied without the inhomogeneous partial wave, a situation which would appear to be the natural limiting case of a surface wave of infinite penetration, the body wave alone satisfies the condition of free surface. This result illuminates the Barnett-Lothe existence theorem for subsonic surface waves. The close connection between the reflection problem and the leaky surface wave problem becomes very apparent in the formalism used. It is shown that for a point on a branch of leaky waves where the solution is undamped, the conditions for simple reflection, i.e. reflection only involving the two body waves, are also present. In the vicinity of such a point reflection is accompanied by resonance excitation of leaky waves. The paper concludes with some explicit calculations for transversely isotropic solids.     

Volume and surface stability of uniformly compressed transversely isotropic linearly elastic materials

The stability loss of a transversely isotropic linearly elastic medium is studied. The medium is uniformly compressed in both horizontal directions, and the initial stress in the vertical direction is equal to zero. The standard analysis based on the Hadamard condition is used. The bifurcation equation divides into two parts, and therefore, two kinds of buckling modes are possible. The critical initial compression is found, but the buckling modes remain indefinite (as the wave length so the relation between the wave numbers is arbitrary). The stability loss of a compressed half-space with a free surface is studied. Only one kind of buckling mode localized near the free surface is possible, and as for an entire space, the buckling mode and the wave length are indefinite. In these problems, linear as well as non-linear approaches are used. In the linear approach, the pre-buckling deformations are ignored. It is shown that for some values of parameters, the linear approach leads to qualitatively incorrect results. The stability loss of an uniformly compressed plate lying on a soft elastic half-space is studied. By using the non-linear post-critical analysis, it is shown that the buckling mode is a chessboard-like one.     

MEASUREMENT AND ANALYSIS OF LASER GENERATED RAYLEIGH AND LAMB WAVES CONSIDERING ITS PULSE DURATION

In this article, laser generated Rayleigh and Lamb waves are studied by taking into account its pulse duration. The physical model and theoretical solution are presented to predict the corresponding waveforms for aluminum samples under the ablation generation regime.The waveforms of the excited Rayleigh and Lamb waves by laser with selected pulse duration were measured by laser interferometer and analyzed theoretically, and the agreement between measurement and analysis is demonstrated for the validation of the theoretical model and solution.The broadening of the Rayleigh wave and the disappearing of high order Lamb wave modes can be found with the increase of the pulse duration by the laser ultrasonic technique.     

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.     

Transient thermoelastic waves in an anisotropic hollow cylinder due to localized heating

This paper presents a theoretical study of transient ultrasonic guided waves generated by concentrated heating of the outer surface of an infinite anisotropic hollow circular cylinder. Generalized thermoelastic theory proposed by Lord and Shulman is adopted to model the dynamic thermoelastic behavior of the cylinder. The concentrated heat source model used is to represent heating due to a pulsed laser beam, which is focused on the outer surface of the cylinder. A semi-analytical finite element (SAFE) method is employed to evaluate guided wave modes in the cylinder. Using integral transform techniques, the modal wave forms are obtained in frequency and wave number domains. Time histories of the propagating modes are then calculated by applying inverse Fourier transformation in the time domain. Numerical results showing the dispersion curves for the group velocities of the propagating modes and transient radial displacements are presented. For this purpose it is assumed that the cylinder is made of transversely isotropic silicon nitride (Si₃N₄). Attention is focused on the propagation characteristics of longitudinal and flexural modes separately.     

Stress Evaluation Using Ultrasonic Interference Spectrum of Leaky Lamb Waves

This paper presents a nondestructive stress evaluation technique using the ultrasonic interference spectrum of leaky Lamb waves. By using a specific pitch-catch ultrasonic setup, the symmetric and antisymmetric modes of Lamb waves in a finite plate are decoupled, leading to simple relationships between the modal frequency spacing of two adjacent modes in the interference spectrum and the acoustic wave velocities that are functions of stress. As a result, the stress in the plate can be determined by measuring the modal frequency spacing instead of the relative flight times to calculate the acoustic wave velocity. Extensive experiments were carried out to verify the viability and robustness of the new technique using a simple testing system. It has been demonstrated that the new technique is about 25 times more accurate than existing flight-time approaches using the same testing system. The experimental results agree well with the results obtained by other ultrasonic methods using expensive equipment.     

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.     

Scattering coefficients for a sphere in a visco-acoustic medium for arbitrary partial wave order

Analytical solutions are reported for the scattering coefficients of a solid elastic sphere suspended in a viscous fluid for arbitrary partial wave order. Expressions are derived for incident compressional and shear wave modes, taking into account the viscosity of the surrounding fluid and resultant wave mode conversion. The long compressional wavelength limit is employed to simplify the derivation, whereas no restriction is placed on the shear wavelength in the fluid compared to the particle dimension. The analytical approximations are compared with numerical results obtained from matrix inversion of the boundary equations and agree within the validity domain of the solutions.