共查询到20条相似文献，搜索用时 218 毫秒
1.
This paper develops a semi-analytical solution for the transient response of an unsaturated single-layer poroviscoelastic medium with two immiscible fluids by using the Laplace transformation and the state-space method. Using the elastic–viscoelastic correspondence principle, we first introduce the Kelvin–Voigt model into Zienkiewicz’s unsaturated poroelastic model. The vibrational response for unsaturated porous material can be obtained by combining these two models and assuming that the wetting and non-wetting fluids are compressible, the solid skeleton and solid particles are viscoelastic, and the inertial and mechanical couplings are taken into account. The Laplace transformation and state-space method are used to solve the basic equations with the associated initial and boundary conditions, and the analytical solution in the Laplace domain is developed. To evaluate the responses in the time domain, Durbin’s numerical inverse Laplace transform method is used to obtain the semi-analytical solution. There are three compressional waves in porous media with two immiscible fluids. Moreover, to observe the three compressional waves clearly, we assume the two immiscible fluids are water and oil. Finally, several examples are provided to show the validity of the semi-analytical solution and to assess the influences of the viscosity coefficients and dynamic permeability coefficients on the behavior of the three compressional waves. 相似文献
2.
单层不可压缩饱和多孔介质一维瞬态响应精确解 总被引：1，自引：0，他引：1
基于Biot理论,考虑惯性、黏滞和机械耦合作用,假定固体颗粒和流体均不可压缩,得到了表面任意竖向荷载作用下单层饱和多孔介质一维瞬态响应的精确解。导出了以固体骨架位移表示的无量纲控制方程,并将边界条件齐次化。求解对应无黏滞耦合作用的特征值问题,得到一组满足齐次边界条件、关于空间坐标的正交函数基。利用变异系数法和基函数的正交性,得到一系列相互解耦的、关于时间的二阶常微分方程及相应的初始条件,并采用状态空间法求解常微分方程,得到位移分量。对整体平衡方程关于空间坐标积分,根据边界条件可确定总应力,并进而求得孔隙压力。通过算例验证所得解法的正确性 相似文献
3.
4.
在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。 相似文献
5.
In this study, a linear poroviscoelastic model based on the Biot theory is proposed to analyze the dynamic response of partially saturated soil. Both the flow-dependent and flow-independent poroviscoelastic behaviors are described in the proposed model. The compressibilities of both the constituents, including the skeleton material and porefluid, are considered, and the effective skeleton stress is determined using the linear viscoelasticity law based on the generalized spring-dashpot model. The soil surface is subjected to two types of harmonic loading, namely, vertical compressive loading and lateral shear loading. The influences of the relaxation time, saturation degree, and soil permeability on the surface and interior dynamic responses of the soil layer are investigated. It is revealed that the amplification factor, effective stress, and pore pressure decrease as the viscous damping increases, indicating that neglecting the viscoelastic property of solid skeleton could overestimate the induced dynamic responses. Furthermore, in comparison with the viscous damping of the solid skeleton, the viscous coupling involving the viscous resisting forces between the solid skeleton and pore fluid has limited effect on the dynamic behavior of the soil layer. 相似文献
6.
研究了横观各向同性饱和土地基在地表动力荷载作用下的三维瞬态响应。基于饱和多孔介质的三维Biot波动理论,利用Laplace变换,建立圆柱坐标系下横观各向同性饱和土的波动方程;解耦波动方程后,根据算子理论,并借助Fourier展开和Hankel变换技术,得到瞬态荷载作用下,饱和土介质的土骨架位移和应力、孔隙水相对位移和孔隙水压力的一般解;利用一般解,给出横观各向同性饱和地基在地表集中荷载激励下的瞬态Lamb问题的解答。数值算例结果表明,采用各向同性饱和介质的动力学模型,不能准确描述具有明显各向异性特性的饱和土地基的瞬态动力特性。 相似文献
7.
Zhendong Shan Daosheng Ling Zhinan Xie Liping Jing Yongqiang Li 《国际地质力学数值与分析法杂志》2019,43(13):2184-2199
One-dimensional transient wave propagation in a saturated single-layer porous medium with a fluid surface layer is studied in this paper. An analytical solution for a special case with a dynamic permeability coefficient kf → ∞ and a semianalytical solution for a general case with an arbitrary dynamic permeability coefficient are presented. The eigenfunction expansion and precise time step integration methods are employed. The solution is presented in series form, and thus, the long-term dynamic responses of saturated porous media with small permeability coefficients can be easily computed. We first transform the nonhomogeneous boundary conditions into homogeneous boundary conditions, and then we obtain the eigenvalues and orthogonal eigenfunctions of the fluid–solid system. Finally, the solutions in the time domain are developed. As the model is one dimensional, geometric attenuation is absent, and only the attenuation in the saturated porous medium is considered. We can apply this model to analyse the influences of different seabed types on the propagation of acoustic waves in the fluid layer, which is very important in ocean acoustics and ocean seismic. This solution can also be employed to validate the accuracies of various numerical methods. 相似文献
8.
基于Zienkiewicz提出的非饱和多孔介质波动理论,考虑两相流体和固体颗粒的压缩性以及惯性、黏滞和机械耦合作用,采用半解析的方法获得了一类典型边界条件下单层非饱和多孔介质一维瞬态响应解。首先推导出无量纲化后以位移表示的控制方程,并将其写成矩阵形式;然后,将边界条件齐次化,求解控制方程所对应的特征值问题,得到了满足齐次边界条件的特征值和相对应的特征函数。根据变异系数法并利用特征函数的正交性,得到了一系列仅黏滞耦合的关于时间的二阶常微分方程及相应的初始条件。在此基础上,运用精细时程积分法给出了常微分方程组的数值解。最后,通过若干算例验证了结果的正确性并探讨了单层非饱和多孔介质一维瞬态动力响应的特点。该方法可推广应用于其他典型的边界条件。 相似文献
9.
Biot's equations of wave propagation through fluid-saturated porous elastic media are discretized spatially using the finite element method in conjunction with Galerkin's procedure. Laplace transformation of the discretized equations is used to suppress the time variable. Introducing Laplace transforms of constituent velocities at nodal points as additional variables, the quadratic set of equations in the Laplace transform parameter is reduced to a linear form. The solution in the Laplace transform space is inverted, term by term, to get the complete time history of the solid and fluid displacements and velocities. Since the solution is exact in the time domain, the error in the calculated response is entirely due to the spatial approximation. The procedure is applied to one-dimensional wave propagation in a linear elastic material and in a fluid-saturated elastic soil layer with ‘weak’, ‘strong’ as well as ‘moderate’ coupling. With refinement of the spatial mesh, convergence to the exact solution is established. The procedure can provide a useful benchmark for validation of approximate temporal discretization schemes and estimation of errors due to spatial discretization. 相似文献
10.
孔隙度和岩石压缩系数等方面的研究在评价油藏弹性产能和动态地质储量方面有重要的应用价值.结合国家地震前兆台网中心数据库8口井的数字化水位资料等, 研究了承压井含水层介质在不排水状态下的孔隙度、固体骨架的体积压缩系数和含水层内水的体积压缩系数.结果表明, 孔隙度与固体骨架的体积压缩系数和含水层内水的体积压缩系数间存在幂函数关系.在第1象限内, 各井含水层介质固体骨架的体积压缩系数随着孔隙度的增大而增大; 含水层内水的体积压缩系数随着孔隙度的增大而减小.固体骨架和含水层内水的体积压缩系数间满足一元二次多项式关系, 且含水层内水的体积压缩系数要比固体骨架的体积压缩系数大, 水更易压缩.另外, 灰岩骨架的压缩系数相对小于砂岩骨架的压缩系数. 相似文献
11.
饱和多孔介质的动力响应研究在众多工程领域具有重要意义。充分考虑孔隙率的变化规律与影响因素,有利于合理揭示饱和多孔介质的相关力学行为。为此,将动态孔隙率模型与用于表征饱和多孔介质动力特性的u-U-p型方程结合,构建相应的非线性力学模型,利用Comsol Multiphysis PDE求取相应的数值解,以此研究不同透水条件下,受谐波载荷激励的二维饱和土体的孔隙率、变形量及孔隙水压力的变化规律。结果表明：孔隙率的变化与土骨架的体应变及孔隙水压力直接相关,土体压缩过程中,孔隙率相应减小,土骨架与孔隙流体的相互作用增强,土体运动时所受阻力增大,其无量纲竖向位移小于孔隙率被视为常数时的情况,在此条件下,由于土体的变形量减小,其孔隙水压力也相对减小。故充分考虑动态孔隙率,有利于更加精确地研究等饱土体和多孔介质的相关力学行为。此外,土体上表面透水条件下,孔隙流体可以从土体表面自由排出,土骨架承受的载荷更大,与不透水条件相比,土体孔隙率、竖向位移、孔隙水压力等变化更为显著。 相似文献
12.
Unsaturated soils are considered as porous continua, composed of porous skeleton with its pores filled by water and air. The governing partial differential equations (PDE) are derived based on the mechanics for isothermal and infinitesimal evolution of unsaturated porous media in terms of skeleton displacement vector, liquid, and gas scalar pressures. Meanwhile, isotropic linear elastic behavior and liquid retention curve are presented in terms of net stress and capillary pressure as constitutive relations. Later, an explicit 3D Laplace transform domain fundamental solution is obtained for governing PDE and then closed‐form analytical transient 3D fundamental solution is presented by means of analytical inverse Laplace transform technique. Finally, a numerical example is presented to validate the assumptions used to derive the analytical solution by comparing them with the numerically inverted ones. The transient fundamental solutions represent important features of the elastic wave propagation theory in the unsaturated soils. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
13.
Chongbin Zhao Lynn B. Reid Klaus Regenauer-Lieb Thomas Poulet 《Computational Geosciences》2012,16(3):735-755
In dealing with chemical-dissolution-front propagation problems in fluid-saturated porous media, the chemical dissolution front represented by the porosity of the medium may have a very steep slope (i.e., a very large porosity gradient) at the dissolution front, depending on the mineral dissolution ratio that is defined as the equilibrium concentration of the dissolved minerals in the pore-fluid to the solid molar density of the dissolvable minerals in the solid matrix. When the mineral dissolution ratio approaches zero, the theoretical value of the porosity gradient tends to infinity at the chemical dissolution front. Even for a very small value of the mineral dissolution ratio, which is very common in geochemical systems, the porosity gradient can be large enough to cause the solution hard to converge when the conventional finite element method is used to solve a chemical dissolution problem in a fluid-saturated porous medium where the pore-fluid is compressible. To improve the convergent speed of solution, a porosity-gradient replacement approach, in which the term involving porosity-gradient computation is replaced by a new term consisting of pore-fluid density-gradient and pressure-gradient computation, is first proposed and then incorporated into the finite element method in this study. Through comparing the numerical results obtained from the proposed approach with the theoretical solutions for a benchmark problem, it has been demonstrated that not only can the solution divergence be avoid, but also the accurate simulation results can be obtained when the proposed porosity-gradient replacement approach is used to solve chemical-dissolution-front propagation problems in fluid-saturated porous media including pore-fluid compressibility. 相似文献
14.
The authors propose a semi‐analytical approach to studying wave propagation in multilayered poroviscoelastic grounds due to transient loads. The theoretical development is based on the exact stiffness matrix method for the Biot theory coupled with a matrix conditioning technique. It is developed in the wavenumber frequency domain after a Fourier transform on the surface space variables and the time variable. The usual methods yield a poorly conditioned numerical system. This is due in particular to the presence of mismatched exponential terms. In this article, increasing exponential terms are eliminated and only decreasing exponential terms remain. Consequently, the method can be applied to a large field of configurations without restriction concerning high frequencies, large Fourier transform parameters or large layer thicknesses. Validation and efficiency of the method are discussed. Effects of layering show that the layer impedance influence on solid and fluid displacements. Moreover, this approach can be of interest for the validation of numerical tools. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
15.
The complete solution is presented for the transient effects of pumping fluid from a point sink embedded in a saturated, porous elastic half-space. It is assumed that the medium is homogeneous and isotropic with respect to its elastic properties and homogeneous but anisotropic with respect to the flow of pore fluid. The soil skeleton is modelled as a linear elastic material obeying Hooke's law, while the pore fluid is assumed to be incompressible with its flow governed by Darcy's law. The solution has been evaluated for a particular value of Poisson's ratio of the solid skeleton, i.e. 0.25, and the results have been presented graphically in the form of isochrones of excess pore pressure and surface profile for the half-space. The solutions presented may have application in practical problems such as dewatering operations in compressible soil and rock masses. 相似文献
16.
An analytical solution to the two-dimensional wave propagation in fluid-saturated half-space subjected to a strip load with vertical harmonic oscillation at the surface is presented. The basic equations have been derived on the basis of Biot's linear theory of poro-elasticity and then solved using Fourier complex transform for the horizontal direction. The importance of a number of soil characteristics including compressibility, degree of saturation and soil permeability has been examined. It is shown that the effect of pore fluid is dominant only for fully saturated soils with incompressible solid grains and low permeability. For partially saturated, compressible or very permeable soils, the stresses would be mainly transferred to solid part and there will be considerable reduction in pore pressure amplitude. 相似文献
17.
Modeling of unsaturated granular flows by a two-layer approach 总被引：2，自引：1，他引：1
Flows of partially saturated grain-fluid mixtures over complex curved topography are commonly observed in nature. However, comprehensive understanding of the physics behind them is to date out of reach. To investigate their dynamic process, a two-layer approach is proposed, in which the fluid-saturated granular layer is overlaid by the pure granular material. More specifically, the lower layer is described by a two-phase mixture theory of density preserving solid and fluid constituents. For the upper layer, the single-phase granular mass is treated as a frictional Coulomb-like continuum, and the dilation effect and the influence of the interstitial air are ignored. The capillarity effects and grains-size segregation are not considered in both the layers. The lower and upper layers interact at an interface which is a material surface for the fluid phase, but across which the mass exchange for the granular phase may take place. The granular mass exchange across the layer interface is parameterized by an entrainment type postulate. In addition, the classical jump conditions are employed to connect both layers at the interface dynamically. Furthermore, we perform the depth-averaged technique for the saturated grain-fluid mixture lower layer and the pure granular upper layer, respectively, to simplify the governing equations established. It is demonstrated that the resulting model equations can be reduced to most of the existing single-layer pure granular flow models and saturated two-phase single-layer debris flow models. Numerical solutions demonstrate that the present two-layer model can describe flows of partially saturated grain-fluid mixtures and the transition process of a saturated grain-fluid mixture into an under-saturated state. 相似文献
18.
A Laplace transform is used to solve the problem of the steady state and transient response of a pinned head pile embedded into a viscoelastic Winkler soil medium. The pile is modeled as an Euler–Bernoulli beam while the soil medium is modeled using a Winkler subgrade approach. Two analytical solutions are developed to specifically address both steady state and transient loads encountered during dynamic pile testing. After choosing a proper contour integration in the complex plane, inverse integration is evaluated. The steady state solutions are associated to the residues of the integration around the poles while the transient solutions are associated to the integration paths along the contour integration. The derived solutions are applied to a case history for which results of dynamic pile tests are available. Dynamic pile flexion is generated by delivering eccentric impact using a dynamic loading test module. Validity of the proposed solution is discussed basing on geotechnical campaign and recorded pile head bending moment and rotation rate. 相似文献
19.
横观各向同性饱和土体三维粘弹性动力分析 总被引：2，自引：2，他引：0
采用针对横观各向同性饱和土体u-w形式三维粘弹性动力方程,考虑土骨架的粘弹性性质且基于粘弹性理论,通过运用Fourier 展开、Laplace 和Hankel 积分变换方法和引入中间变量,将含有粘弹性参数的六元二阶偏微分运动控制方程组,化为2组各含4个未知变量的常微分方程组,从而给出了柱坐标系下粘弹性横观各向同性饱和土体在非轴对称动力荷载作用下的瞬态反应的土骨架位移分量、孔隙流体相对于土骨架的位移分量瞬态反应一般解。在此基础上,引入初始条件和边界条件,对垂直向和水平向动力荷载作用下半空间边值问题进行了求解。根据动力时域解答的一般解,利用Laplace和Hankel 数值逆变换技术,编制了相应的数值计算程序。并进行了实例验证和弹性、粘弹性解的对比分析。结果表明,在进行横观各向同性饱和土体动力分析时,考虑土骨架的粘弹性是必要的。 相似文献