粘弹性损伤混凝土介质中波的传播

混凝土的损伤会引起弹性和粘性性质的改变,利用弹性和粘性的双参数来描述损伤,建立双参数损伤理论。根据不同的损伤程度对混凝土介质进行分区处理,建立基本方程式,考虑连续性条件和边界条件,求解波动方程。比较是否考虑粘性时的波幅和波传播时间与损伤的关系,表明粘性是损伤混凝土介质中波传播所必须要考虑的问题。分析了粘弹性混凝土介质中损伤区域长度、损伤度等对波传播的影响,给出了它们的关系曲线,可为波的反分析提供依据。     

The Peculiarities of Linear Wave Propagation in Double Porous Media

The features of propagation of longitudinal and transverse waves (LW and TW) in fractured porous medium (FPM) saturated with liquid are investigated by methods of multiphase mechanics. The mathematical model of FPM accounting for inequality of velocities and pressures of liquid in pores and fractures, liquid mass exchange and nonstationary interaction forces is developed. Processes of monochromatic wave propagation are studied. The dispersion relation is obtained and the effect of model parameters on wave propagation is analysed. It is established that one transverse and three longitudinal waves propagate in FPM saturated with liquid. The fastest LW is a deformational wave and the two others are filtrational. Filtrational waves attenuate much stronger than deformational and transverse waves. Distinction of velocities and pressures in liquid in various pore systems provides an explanation for the existence of the two filtrational waves in porous medium with two different characteristic sizes of pores.     

Seismic Wave Propagation in Composite Elastic Media

It has been known since the time of Biot–Gassman theory (Biot, J Acoust Soc Am 28:168–178, 1956, Gassmann, Naturf Ges Zurich 96:1–24, 1951) that additional seismic waves are predicted by a multicomponent theory. It is shown in this article that if the second or third phase is also an elastic medium then multiple p and s waves are predicted. Futhermore, since viscous dissipation no longer appears as an attenuation mechanism and the media are perfectly elastic, these waves propagate without attenuation. As well, these additional elastic waves contain information about the coupling of the elastic solids at the pore scale. Attempts to model such a medium as a single elastic solid causes this additional information to be misinterpreted. In the limit as the shear modulus of one of the solids tends to zero, it is shown that the equations of motion become identical to the equations of motion for a fluid filled porous medium when the viscosity of the fluid becomes zero. In this limit, an additional dilatational wave is predicted, which moves the fluid though the porous matrix much similar to a heart pumping blood through a body. This allows for a connection with studies which have been done on fluid-filled porous media (Spanos, 2002).     

不可压饱和粘弹性多孔介质固结问题的Gurtin型变分原理

基于多孔介质理论,在固相骨架和孔隙流体微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,利用卷积积分的性质,本文首先建立了以固相骨架位移、孔隙流体相对速度和孔隙流体压力为宗量的流体饱和粘弹性多孔介质固结问题的一个Gurtin型变分原理.其次,利用Lagrange乘子法解除相关的变分约束条件,建立了流体饱和粘弹性多孔介质固结问题的若干广义Gurtin型变分原理,包括第三类的Hu-Washizu型变分原理.最后,简单讨论了等价初边值问题的相应变分原理.这些Gurtin型变分原理的建立不仅丰富了饱和粘弹性多孔介质的相关理论,而且为相关数值模拟方法,如有限元法、无网格法等的建立奠定了理论基础.     

The Generalized Reflection and Transmission Matrix Method for Wave Propagation in Stratified Fluid-Saturated Porous Media

A systematic and efficient algorithm, the generalized reflection and transmission matrix method, has been developed for wave propagation in stratified fluid-saturated poroelastic half-space. The proposed method has the advantage of computational efficiency and numerical stability for high frequencies and large layer thickness. A wide class of seismic sources, ranging from a single-body force to double couples, is introduced by utilizing the moment tensor concept. In order to validate the proposed algorithm, we applied our formulation to calculate wave fields in a homogeneous poroelastic half-space. It is shown that the numerical results computed with the present approach agree well with those computed with the analytical solution. Numerical examples for a two-layer model subjected to various sources such as double couple, dipole, and explosive sources are provided. From the waveforms of surface displacements, the arrivals of transmitted and converted PS and SP waves at the interface of the two-layer model can be clearly observed. As expected, it is impossible to observe the arrivals of transmitted $S$ and transmitted and converted SP waves from the waveforms induced by fluid withdrawal.     

Propagation and Evolution of Wave Fronts in Two-Phase Porous Media

An investigation is made into the propagation and evolution of wave fronts in a porous medium which is intended to contain two phases: the porous solid, referred to as the skeleton, and the fluid within the interconnected pores formed by the skeleton. In particular, the microscopic density of each real material is assumed to be unchangeable, while the macroscopic density of each phase may change, associated with the volume fractions. A two-phase porous medium model is concisely introduced based on the work by de Boer. Propagation conditions and amplitude evolution of the discontinuity waves are presented by use of the idea of surfaces of discontinuity, where the wave front is treated as a surface of discontinuity. It is demonstrated that the saturation condition entails certain restrictions between the amplitudes of the longitudinal waves in the solid and fluid phases. Two propagation velocities are attained upon examining the existence of the discontinuity waves. It is found that a completely coupled longitudinal wave and a pure transverse wave are realizable in the two-phase porous medium. The discontinuity strength of the pore-pressure may be determined by the amplitude of the coupled longitudinal wave. In the case of homogeneous weak discontinuities, explicit evolution equations of the amplitudes for two types of discontinuity waves are derived.     

Propagation of Elastic Waves in Periodically Inhomogeneous Media

The structural theory of waves and vibrations in periodically inhomogeneous media is set out. Relevant research results are presented. Emphasis is on the principles of the theory, surface waves, and other surface effects     

损伤介质中波传播特性的研究

本文研究了弹性波在损伤介质中传播的基本理论,建立了基本动力学方程式,并根据其边界条件和连续性条件求得末端接收波位移的时域曲线。接收波则综合携带了损伤介质内部特性信息,可通过这些信息研究损失介质损坏的程度及位置。通过实例计算,定量地分析了结构内部损伤度、损伤位置和损伤区域大小对接收波幅值和波到达时间的影响。通过分析可知损伤对波传播的影响较大;并且可以根据给定的曲线图,直接由测得的接收波峰值及到达时间探测到结构内部损伤大小与分布。这为无损检测的定量化反分析提供了可靠的依据。     

Wave Propagation along Axes of Symmetry in Linearly Elastic Media with Initial Stress

The proofs of the Fedorov–Stippes and Fedorov's theorems, which hold for linearly elastic homogeneous bodies in natural configurations, remain valid for any linearly elastic medium with initial stress provided that Hooke's tensor be replaced by a suitable elasticity tensor.     

Attenuated Wave Field in Fluid-Saturated Porous Medium with Excitations of Multiple Sources

This research addresses the investigation of an elastic wave field in a homogeneous and isotropic porous medium which is fully saturated by a Newtonian viscous fluid. A new methodology is developed for describing the wave field in the medium excited by multiple energy sources. To quantify the relative displacements between the fluid and solid of the medium, the governing equations of the elastic wave propagation are derived in the form of displacements specially. The velocities and attenuation of the waves are considered as functions of viscosity and frequency. Making use of the Hankel function and the moving-coordinate method, a model of the wave motion with multiple cylindrical wave sources is built. Making use of the model established in this research, the relative displacement between the fluid and the solid can be quantified, and the wave field in the porous media can then be determined with the given energy sources. Numerical simulations of cylindrical waves from multiple energy sources propagating in the porous medium saturated by viscous fluid are performed for demonstrating the practicability of the model developed.     

A Computational Method for Free-Surface Viscous Flows Through Porous Media

In this paper we have investigated free-surface unsteady viscous flow through porous media using a finite difference technique via primitive formulation. A transformation has been used to convert the irregular flow domain into a rectangular domain. Two illustrations have been presented in case of a standing wave problem and a cavity flow problem with a free surface. The method is found satisfactory for numerical solution for Reynolds number upto 200.     

非均匀吸收介质中地震波的广义传播矩阵解法

本文在复频域内,通过应用混合变量粘弹性波方程和线性常微分方程组的指数矩阵解法,给出了一种计算非均匀吸收介质中地震波传播的广义传播矩阵解法。该方法适用于各种粘弹性模型,可模拟任意震源及所产生的各种体波、面波,数值结果表明具有很高的计算精度。     

Application of a Direct Stiffness Method to Wave Propagation in Multiphase Poroelastic Media

In this paper, it is demonstrated how a direct stiffness method for wavepropagation in multilayered saturated poroelastic media, based on integraltransform techniques, can be modified to account for a small amount of gasin the pores. Unsaturated media with small gas fractions can be treatedusing Smeulders extension of Biots poroelastic theory. The effect of thepresence of gas bubbles on the fluid bulk modulus and the dispersioncharacteristics of a water-saturated sand of Mol is demonstrated. Thedirect stiffness method is illustrated with a numerical example wheretransient wave propagation in a dry, saturated and unsaturated halfspaceis considered.     

Rayleigh–Taylor Instability of Rotating Oldroydian Viscoelastic Fluids in Porous Medium in Presence of a Variable Magnetic Field

The instability of the plane interface between two oldroydian viscoelastic superposed fluids in the presence of uniform rotation and variable magnetic field in porous medium is considered. For potentially stable configuration, the system is found to be stable for disturbances of all wave numbers. The magnetic field succeeds in stabilizing certain wave-number range, which were unstable in the absence of magnetic field and rotation for the potentially unstable configuration. Sub cases of magnetic free and rotation free configurations are also considered, separately.     

可变形多孔介质渗透系数的测定方法

在Ｂｉｏｔ理论基础上给出可变形多孔介质耦合渗流基本方程；求出小试件一维定常耦合渗流问题的解答；表明在一维流固耦合情况下试件内部压力梯度有明显的非均匀性。因此通过实验确定可变形多孔介质渗透系数在数学上可归结为微分方程的反问题，传统的测试渗透系数的方法需要改进。介绍了可变形多孔介质渗透系数的测试原理和测试方法。对粒状多孔材料实验的结果表明，传统实验方法得到的渗透系数误差较大。     

Validation of the Slow Compressional Wave in Porous Media: Comparison of Experiments and Numerical Simulations

We perform numerical simulation of ultrasonic experiments on poroelastic samples, in which Biot's slow compressional wave had been observed. The simulation is performed using OASES modeling code, which allows to compute elastic wave fields in layered poroelastic media. Modeled were the experiments of Plona (1980), Rasolofosaon (1988), and our own measurements. In all the three situations, a good agreement between experiment and simulations has been observed. This further confirms the fact that Biot's theory of poroelasticity, on which the simulations were based, adequately describes the behavior of the porous materials under investigations at ultrasonic frequencies.     

Propagation of Transient Acoustic Waves in Layered Porous Media: Fractional Equations for the Scattering Operators

Acoustic waves scattering from a rigid air-saturated porous medium is studied in the time domain. The medium is one dimensional and its physical parameters are depth dependent, i.e., the medium is layered. The loss and dispersion properties of the medium are due to the fluid-structure interaction induced by wave propagation. They are modeled by generalized susceptibility functions which express the memory effects in the propagation process. The wave equation is then a fractional telegraphists equation. The two relevant quantities are the scattering operators—transmission and reflection operators—which give the scattered fields from the incident wave. They are obtained from Volterra equations which are fractional equations for the scattering operators.     

油水两相渗流压力脉动的实验及数值分析

本文用高精度的压差传感器和高速度的数据采集设备同时测量出油水两相流体在多孔介质流动时产生的较大压力降和微小压力脉动,应用自行研制的大型数据处理软件对恒流速水驱油的两相渗流压力脉动实验数据进行了分析,发现不同阶段压力脉动具有明显的频谱特性和时间相关特性的不同,在第三阶段（油为主,水增加阶段）谱能增加最大,时间正相关程度最强。     

Time discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium

A time-discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium is proposed. As compared with the existing discontinuous Galerkin finite element methods, the distinct feature of the proposed method is that the continuity of the displacement vector at each discrete time instant is automatically ensured, whereas the discontinuity of the velocity vector at the discrete time levels still remains. The computational cost is then obviously reduced, particularly, for material non-linear problems. Both the implicit and explicit algorithms to solve the derived formulations for material non-linear problems are developed. Numerical results show a good performance of the present method in eliminating spurious numerical oscillations and providing with much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain. The project supported by the National Natural Science Foundation of China (19832010, 50278012, 10272027) and the National Key Basic Research and Development Program (973 Program, 2002CB412709)     

A Mathematical Model for Hysteretic Two-Phase Flow in Porous Media

We develop a mathematical model for hysteretic two-phase flow (of oil and water) in waterwet porous media. To account for relative permeability hysteresis, an irreversible trapping-coalescence process is described. According to this process, oil ganglia are created (during imbibition) and released (during drainage) at different rates, leading to history-dependent saturations of trapped and connected oil. As a result, the relative permeability to oil, modelled as a unique function of the connected oil saturation, is subject to saturation history. A saturation history is reflected by history parameters, that is by both the saturation state (of connected and trapped oil) at the most recent flow reversal and the most recent water saturation at which the flow was a primary drainage. Disregarding capillary diffusion, the flow is described by a hyperbolic equation with the connected oil saturation as unknown. This equation contains functional relationships which depend on the flow mode (drainage or imbibition) and the history parameters. The solution consists of continuous waves (expansion waves and constant states), shock waves (possibly connecting different modes) and stationary discontinuities (connecting different saturation histories). The entropy condition for travelling waves is generalized to include admissible shock waves which coincide with flow reversals. It turns out that saturation history generally has a strong influence on both the type and the speed of the waves from which the solution is constructed.