共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
饱和多孔介质的动力响应研究在众多工程领域具有重要意义。充分考虑孔隙率的变化规律与影响因素,有利于合理揭示饱和多孔介质的相关力学行为。为此,将动态孔隙率模型与用于表征饱和多孔介质动力特性的u-U-p型方程结合,构建相应的非线性力学模型,利用Comsol Multiphysis PDE求取相应的数值解,以此研究不同透水条件下,受谐波载荷激励的二维饱和土体的孔隙率、变形量及孔隙水压力的变化规律。结果表明：孔隙率的变化与土骨架的体应变及孔隙水压力直接相关,土体压缩过程中,孔隙率相应减小,土骨架与孔隙流体的相互作用增强,土体运动时所受阻力增大,其无量纲竖向位移小于孔隙率被视为常数时的情况,在此条件下,由于土体的变形量减小,其孔隙水压力也相对减小。故充分考虑动态孔隙率,有利于更加精确地研究等饱土体和多孔介质的相关力学行为。此外,土体上表面透水条件下,孔隙流体可以从土体表面自由排出,土骨架承受的载荷更大,与不透水条件相比,土体孔隙率、竖向位移、孔隙水压力等变化更为显著。 相似文献
3.
In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
4.
This paper presents a numerical model for the elasto‐plastic electro‐osmosis consolidation of unsaturated clays experiencing large strains, by considering electro‐osmosis and hydro‐mechanical flows in a deformable multiphase porous medium. The coupled governing equations involving the pore water flow, pore gas flow, electric flow and mechanical deformation in unsaturated clays are derived within the framework of averaging theory and solved numerically using finite elements. The displacements of the solid phase, the pressure of the water phase, the pressure of the gas phase and the electric potential are taken as the primary unknowns in the proposed model. The nonlinear variation of transport parameters during electro‐osmosis consolidation are incorporated into the model using empirical expressions that strongly depend on the degree of water saturation, whereas the Barcelona Basic Model is employed to simulate the elasto‐plastic mechanical behaviour of unsaturated clays. The accuracy of the proposed model is evaluated by validating it against two well‐known numerical examples, involving electro‐osmosis and unsaturated soil behaviour respectively. Two further examples are then investigated to study the capability of the computational algorithm in modelling multiphase flow in electro‐osmosis consolidation. Finally, the effects of gas generation at the anode, the deformation characteristics, the degree of saturation and the time dependent evolution of the excess pore pressure are discussed. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
5.
General coupling extended multiscale FEM for elasto‐plastic consolidation analysis of heterogeneous saturated porous media 
下载免费PDF全文
下载免费PDF全文 This paper presents a general coupling extended multiscale FEM (GCEMs) for solving the coupling problem of elasto‐plastic consolidation of heterogeneous saturated porous media. In the GCEMs, the numerical multiscale base functions for the solid skeleton and fluid phase of the coupling system are all constructed on the basis of the equivalent stiffness matrix of the unit cell, which not only contain the interaction between the solid and fluid phases but also consider the time effect. Furthermore, in order to improve the computational accuracy for two‐dimensional problems, a multi‐node coarse element strategy for the GCEMs is proposed, and a two‐scale iteration algorithm for the elasto‐plastic consolidation analysis is developed. Some one‐dimensional and two‐dimensional homogeneous and heterogeneous numerical examples are carried out to validate the proposed method through the comparison with the coupling multiscale FEM and standard FEM. Numerical results show that the newly developed GCEMs can almost preserve the same convergent property as the standard FEM and also possesses the advantages of high computational efficiency. In addition, the GCEMs can be easily applied to other coupling multifield and multiphase transient problems. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
6.
This paper presents a fully coupled finite element formulation for partially saturated soil as a triphasic porous material, which has been developed for the simulation of shield tunnelling with heading face support using compressed air. While for many numerical simulations in geotechnics use of a two‐phase soil model is sufficient, the simulation of compressed air support demands the use of a three‐phase model with the consideration of air as a separate phase. A multiphase model for soft soils is developed, in which the individual constituents of the soil—the soil skeleton, the fluid and the gaseous phase—and their interactions are considered. The triphasic model is formulated within the framework of the theory of porous media, based upon balance equations and constitutive relations for the soil constituents and their mixture. An elasto‐plastic, cam–clay type model is extended to partially saturated soil conditions by incorporating capillary pressure according to the Barcelona basic model. The hydraulic properties of the soil are described via DARCY 's law and the soil–water characteristic curve after VAN GENUCHTEN . Water is modelled as an incompressible and air as a compressible phase. The model is validated by means of selected benchmark problems. The applicability of the model to geotechnical problems is demonstrated by results from the simulation of a compressed air intervention in shield tunnelling. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
7.
Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
8.
This study aims at determining the macroscopic strength of porous materials having a Drucker–Prager solid phase at microscale and two populations of saturated pores with different pressures at both micro and meso scales. To this end, and taking account of the available results by Maghous et al. (2009), we first derive a closed‐form expression of approximate criterion for a dry porous medium whose matrix obeys to a general elliptic criterion. The methodology to formulate this criterion is based on limit analysis of a hollow sphere subjected to a uniform strain rate boundary conditions. The obtained results are then implemented in a two‐step homogenization procedure, which interestingly delivers analytical expression of the macroscopic criterion for dry double porous media whose solid phase at microscale obeys to a Drucker–Prager criterion. After a brief discussion of the results, we propose an extension to double porous saturated media, allowing therefore to quantify the simultaneous effects of the different pore pressures applied on each voids population. The results are discussed in terms of the existence or not of effective stresses. Finally, they are assessed by comparing them to recently available results. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
9.
Extreme waves caused by tsunamis and storm surges can lead to soil failures in the near‐shore region, which may have severe impact on coastal environments and communities. Multiphase flows in deformable porous media involve several coupled processes and multiple time scales, which are challenging for numerical simulations. The objective of this study is to investigate the roles of the various processes and their interactions in multiphase flows in unsaturated soils under external wave loading, via theoretical time‐scale analysis and numerical simulations. A coupled geomechanics–multiphase flow model based on conservation laws is used. Theoretical analysis based on coupled and decoupled models demonstrates that transient and steady‐state responses are governed by pore pressure diffusion and saturation front propagation, respectively, and that the two processes are essentially decoupled. Numerical simulations suggest that the compressibility of the pore fluids and the deformation of the soil skeleton are important when the transient responses of the media are of concern, while the steady‐state responses are not sensitive to these factors. The responses obtained from the fully coupled numerical simulations are explained by a simplified time‐scale analysis based on coupled and decoupled models. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
10.
The empirical Darcy's law of water transport in porous media, Fick's law of chemical diffusion, and Fourier's law of thermal transport have been widely used in geophysics/geochemistry for over 150 years. However, the strong couplings between water, temperature, and chemicals in a membrane porous medium have made these laws inapplicable and present a significant hurdle to the understanding of multiphase flow in such a material. Extensive experiments over the past century have observed chemical osmosis and thermal osmosis, but a model for understanding their underlying physicochemical basis has remained unavailable, because of the highly cross‐disciplinary and multiscale‐multiphase nature of the coupling. Based on the fundamental principles of nonequilibrium thermodynamics and mixture coupling theory, a rigorously theoretical and mathematical framework is proposed and a general model accounting for all of the coupled influences is developed. This leads to a simple and robust mathematical matrix for studying multiphase couplings in a membrane porous medium when all chemical components are electrically neutral. 相似文献
11.
基于Biot理论,研究了饱和土中带有衬砌的圆形隧洞在移动环形荷载作用下的动力响应。假定衬砌为弹性体,土体为饱和多孔介质,引入两类势函数来表示土体、孔隙水和衬砌的位移,使隧洞的控制方程解耦。结合边界条件及连续条件,通过傅立叶变换得到频率-波数域中衬砌和土体的应力、位移和孔隙水压力解答,最后用傅立叶积分逆变换得到时-空域中的数值解。计算并比较了3种隧洞模型（弹性土体隧洞、饱和土体隧洞和饱和土衬砌隧洞）的动力响应分析。数值分析结果说明：（1）移动荷载速度对3种隧洞动力响应均具有较大影响;（2）弹性土体隧洞和饱和土体隧洞的动力响应具有明显区别,所以在富水地区的隧洞动力响应中土体应该视为饱和土体;（3）衬砌对隧洞动力响应有较大影响,故隧洞的动力分析中不能忽略衬砌作用。 相似文献
12.
We describe a new approach for simulation of multiphase flows through heterogeneous porous media, such as oil reservoirs.
The method, which is based on the wavelet transformation of the spatial distribution of the single-phase permeabilities, incorporates
in the upscaled computational grid all the relevant data on the permeability, porosity, and other important properties of
a porous medium at all the length scales. The upscaling method generates a nonuniform computational grid which preserves the resolved structure
of the geological model in the near-well zones as well as in the high-permeability sectors and upscales the rest of the geological
model. As such, the method is a multiscale one that preserves all the important information across all the relevant length
scales. Using a robust front-detection method which eliminates the numerical dispersion by a high-order total variation diminishing
method (suitable for the type of nonuniform upscaled grid that we generate), we obtain highly accurate results with a greatly
reduced computational cost. The speed-up in the computations is up to over three orders of magnitude, depending on the degree
of heterogeneity of the model. To demonstrate the accuracy and efficiency of our methods, five distinct models (including
one with fractures) of heterogeneous porous media are considered, and two-phase flows in the models are studied, with and
without the capillary pressure. 相似文献
13.
This paper describes a mathematical model and adaptive numerical simulation of the time‐dependent multiphase, multicomponent
flow which occurs when a gas venting process is used to remove a volatile contaminant from a porous medium. The numerical
simulation is adaptive in both space and time and involves the use of a finite element spatial discretization and the SPRINT2D
software [6] for time integration. Results are presented which demonstrate the high quality of the simulation, both in terms
of the length scales of the features that are resolved and the efficiency of the solutions relative to those obtained on fixed
grids.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
14.
In this paper, a mathematical model is presented for the analysis of dynamic fracture propagation in the saturated porous media. The solid behavior incorporates a discrete cohesive fracture model, coupled with the flow in porous media through the fracture network. The double‐nodded zero‐thickness cohesive interface element is employed for the mixed mode fracture behavior in tension and contact behavior in compression. The crack is automatically detected and propagated perpendicular to the maximum effective stress. The spatial discretization is continuously updated during the crack propagation. Numerical examples from the hydraulic fracturing test and the concrete gravity dam show the capability of the model to simulate dynamic fracture propagation. The comparison is performed between the quasi‐static and fully dynamic solutions, and the performance of two analyses is investigated on the values of crack length and crack mouth opening. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
15.
分数导数模型描述的饱和土桩纵向振动分析 总被引：7，自引：0，他引：7
土体具有黏弹性性质,为了更好地考虑饱和土体固相土骨架的黏弹性性质,在分数导数理论和多孔介质理论的基础上,将土体视为液固饱和两相介质,并利用分数导数模型来描述饱和土固相土骨架的应力-应变关系,建立分数导数模型描述的饱和土的控制方程。在三维轴对称情况下,利用势函数和分离变量法研究了分数导数模型描述的饱和土中桩的振动问题。分析了模型参数对饱和土中桩的竖向振动的影响。研究表明,分数导数模型描述的饱和土的控制方程应用范围更广,模型参数对桩的竖向振动有较大的影响。 相似文献
16.
In this paper, a series of multimaterial benchmark problems in saturated and partially saturated two‐phase and three‐phase deforming porous media are addressed. To solve the process of fluid flow in partially saturated porous media, a fully coupled three‐phase formulation is developed on the basis of available experimental relations for updating saturation and permeabilities during the analysis. The well‐known element free Galerkin mesh‐free method is adopted. The partition of unity property of MLS shape functions allows for the field variables to be extrinsically enriched by appropriate functions that introduce existing discontinuities in the solution field. Enrichment of the main unknowns including solid displacement, water phase pressure, and gas phase pressure are accounted for, and a suitable enrichment strategy for different discontinuity types are discussed. In the case of weak discontinuity, the enrichment technique previously used by Krongauz and Belytschko [Int. J. Numer. Meth. Engng., 1998; 41:1215–1233] is selected. As these functions possess discontinuity in their first derivatives, they can be used for modeling material interfaces, generating only minor oscillations in derivative fields (strain and pressure gradients for multiphase porous media), as opposed to unenriched and constrained mesh‐free methods. Different problems of multimaterial poro‐elasticity including fully saturated, partially saturated one, and two‐phase flows under the assumption of fully coupled extended formulation of Biot are examined. As a further development, problems involved with both material interface and impermeable discontinuities, where no fluid exchange is permitted across the discontinuity, are considered and numerically discussed. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
17.
Three‐dimensional time‐harmonic fundamental solutions for a fluid‐saturated poroelastic half‐space with partially permeable free surface 下载免费PDF全文
下载免费PDF全文 By virtue of a pair of scalar potentials for the displacement of the solid skeleton and the pore fluid pressure field of a saturated poroelastic medium, an alternative solution method to the Helmholtz decomposition is developed for the wave propagation problems in the framework of Biot's theory. As an application, a comprehensive solution for three‐dimensional response of an isotropic poroelastic half‐space with a partially permeable hydraulic free surface under an arbitrarily distributed time‐harmonic internal force field and fluid sources is developed. The Green's functions for the poroelastic fields, corresponding to point, ring, and disk loads, are reduced to semi‐infinite complex‐valued integrals that can be evaluated numerically by an appropriate quadrature scheme. Analytical and numerical comparisons are made with existing elastic and poroelastic solutions to illustrate the quality and features of the solution. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
18.
This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow. 相似文献
19.
Mathematical simulation of non‐isothermal multiphase flow in deformable unsaturated porous media is a complicated issue because of the need to employ multiple partial differential equations, the need to take into account mass and energy transfer between phases and because of the non‐linear nature of the governing partial differential equations. In this paper, an analytical solution for analyzing a fully coupled problem is presented for the one‐dimensional case where the coefficients of the system of equations are assumed to be constant for the entire domain. A major issue is the non‐linearity of the governing equations, which is not considered in the analytical solution. In order to introduce the non‐linearity of the equations, an iterative discretized procedure is used. The domain of the problem is divided into identical time–space elements that cover the time–space domain. A separate system of equations is defined for each element in the local coordinate system, the initial and boundary conditions for each element are obtained from the adjacent elements and the coefficients of the system of equations are considered to be constant in each step. There are seven governing differential equations that should be solved simultaneously: the equilibrium of the solid skeleton, mass conservation of fluids (water, water vapor and gas) and energy conservation of phases (solid, liquid and gas). The water vapor is not in equilibrium with water and different phases do not have the same temperature. The governing equations that have been solved seem to be the most comprehensive in this field. Three examples are presented for analyzing heat and mass transfer in a semi‐infinite column of unsaturated soil. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
20.
An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献