共查询到20条相似文献，搜索用时 46 毫秒
1.
Cracks on a natural soil surface provide preferential pathways for water infiltration and contaminant solute transport. Before the mechanical property and permeability of cracked soil can be studied, a crack network model must first be established. Based on statistical analysis of crack geometry from field observations, a new method of representing a 3D crack network was developed. A horizontal plane of a crack network was derived using the Voronoi diagram. Each vertical crack was idealized as an inverted triangular prism. The 3D permeability tensor was determined by modeling the water flow through the crack network. A computer program was developed to generate a 3D crack network automatically and compute the permeability tensor of cracked soil at different depths. The model was verified by comparing the measured permeability and computed permeability of cracked soil. The results showed that the computed permeability was consistent with the measured permeability. 相似文献
2.
A two‐scales numerical analysis is set up in order to upscale the permeability of fractured materials such as concrete. To that aim, we couple finite element (FE) kinematics enhancements (strong discontinuities) representing fine scale cracks to the fine scale permeability tensor. The latter may be split into two parts: the first one is isotropic and corresponds to flows within the porosity of the material; the second one, based upon a set of cracks with different orientations and openings, is anisotropic. For the latter, each crack is a path for mass flow according to the Poiseuille law considering two infinite planes. We show how the upscaling procedure leads both to the definition of macroscopic permeability tensors as well as the flow rate evaluation for components of concrete structures. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
3.
The model proposed in this article relates permeability to porosity measurements that can easily be performed in the laboratory. The pore size distribution (PSD) curve is updated with strains and damage. The updated volumetric fractions of natural pores and cracks are introduced in the expression of permeability. Contrary to classical permeability models based on PSD integrations, the model proposed in this article accounts for possible changes in the porosity modes: one mode for undamaged samples and two modes for cracked samples. The proposed approach also accounts for varying states of damage, as opposed to classical fracture network models, in which the cracks pattern is fixed. The only material parameters that are required to describe the microstructure are the lower and upper bounds of the pores size for both natural pores and cracks. All the other PSD parameters involved in the model are related to macroscopic parameters that can easily be determined in the laboratory, such as the initial void ratio. The framework proposed in this article can be used in any damage constitutive model to determine the permeability of a brittle porous medium. Drained triaxial compression tests have been simulated. Before cracks initiation, permeability decreases while the larger natural pores are getting squeezed. After the occurrence of damage, permeability grows due to the increase of cracks density. The model performs well to represent the influence of the confining pressure on damage evolution and permeability variations. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
4.
In a cracked material, the stress intensity factors (SIFs) at the crack tips, which govern the crack propagation and are associated with the strength of the material, are strongly affected by the crack inclination angle and the orientation with respect to the principal direction of anisotropy. In this paper, a formulation of the boundary element method (BEM), based on the relative displacements of the crack tip, is used to determine the mixed‐mode SIFs of isotropic and anisotropic rocks. Numerical examples of the application of the formulation for different crack inclination angles, crack lengths, and degree of material anisotropy are presented. Furthermore, the BEM formulation combined with the maximum circumferential stress criterion is adopted to predict the crack initiation angles and simulate the crack propagation paths. The propagation path in cracked straight through Brazilian disc specimen is numerically predicted and the results of numerical and experimental data compared with the actual laboratory observations. Good agreement is found between the two approaches. The proposed BEM formulation is therefore suitable to simulate the process of crack propagation. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
5.
For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
6.
This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
7.
This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek−1 is 𝒪(1). These variations are estimated to be less than 5% if Ek−1<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
8.
Reservoir depletion results in rock failure, wellbore instability, hydrocarbon production loss, oil sand production, and ground surface subsidence. Specifically, the compaction of carbonate reservoirs with soft rocks often induces large plastic deformation due to rock pore collapse. On the other hand, following the compaction of reservoirs and failure of rock formations, the porosity and permeability of formations will, in general, decrease. These bring a challenge for reservoir simulations because of high nonlinearity of coupled geomechanics and fluid flow fields. In this work, we present a fully implicit, fully coupled, and fully consistent finite element formulation for coupled geomechanics and fluid flow problems with finite deformation and nonlinear flow models. The Pelessone smooth cap plasticity model, an important material model to capture rock compaction behavior and a challenging material model for implicit numerical formulations, is incorporated in the proposed formulation. Furthermore, a stress-dependent permeability model is taken into account in the formulation. A co-rotational framework is adopted for finite deformation, and an implicit material integrator for cap plasticity models is consistently derived. Furthermore, the coupled field equations are consistently linearized including nonlinear flow models. The physical theories, nonlinear material and flow models, and numerical formulations are the focus of part I of this work. In part II, we verify the proposed numerical framework and demonstrate the performance of our numerical formulation using several numerical examples including a field reservoir with soft rocks undergoing serious compaction. 相似文献
9.
A new anisotropic poroelastic damage model is proposed for saturated brittle porous materials. The model is formulated in the framework of the continuum damage mechanics. A second‐rank symmetric tensor is used to characterize material damage due to oriented microcracks. The classic Biot poroelastic theory is then extended to include poroelastic damage coupling. Both the deterioration of elastic properties and poroelastic coefficients is taken into account. A suitable procedure for determination of model parameters from standard laboratory tests is presented. The validity of the model is tested through comparison between numerical predictions and experimental data in various loading conditions. The overall performance of the model is evaluated. The choice of relevant effective stress for the microcrack propagation criterion in saturated cohesive geomaterials is discussed. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
10.
A thermodynamic framework is proposed to couple the effect of mechanical stress and temperature on crack opening and closure in rocks. The model is based on continuum damage mechanics, with damage defined as the second-order crack density tensor. The free energy of the damaged rock is expressed as a function of deformation, temperature, and damage. The damage criterion captures mode I crack propagation, the reduction in toughness due to heating, and the increase in energy release rate with cumulated damage. Crack closure is modeled through unilateral effects produced on rock stiffness. The model was calibrated and verified against published experimental data. Thermo-mechanical crack opening (resp. closure) was studied by simulating a triaxial compression test (resp. uniaxial extension test), including a thermal loading phase. The degradation of stiffness due to tensile stress and recovery of stiffness induced by both mechanical and thermo-mechanical unilateral effects are well captured. The thermo-mechanical energy release rate increases with thermal dilation and also decreases with ambient temperature. It was observed that there is a temperature threshold, below which the rock behaves elastically. A parametric study also showed that the model can capture hardening and softening during thermo-mechanical closure (for specific sets of parameters). These numerical observations may guide the choice of rock material used in geotechnical design, especially for nuclear waste disposals or compressed-air storage facilities. 相似文献
11.
12.
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. 相似文献
13.
外部荷载作用下的裂隙扩展在空间上一般是非均匀的,引起岩石材料的衍生各向异性。将材料离散成大量随机分布的由力键连结的物质点,基于力键的方向性,且将局部弹性张量离散成一定数量的方向张量,理论推导出力键模量与宏观弹性参数之间的关系。通过考虑力键断裂效应,建立了各向异性弹性损伤本构模型。为了模拟中等孔隙率岩石在常规三轴压缩试验中脆性向延性转变的力学行为,在力键断裂效应中引入损伤抑制函数。通过模拟Tennessee大理岩和Indiana石灰岩的常规三轴实验,并与试验数据对比,验证了模型的合理性和有效性。 相似文献
14.
Fluid flow during permeation grouting of fine sands with a microcement-based grout is studied by assuming that the heterogeneous medium composed of the initial granular skeleton, filtered cement and the interstitial fluid phase can be replaced by a continuous equivalent medium at the macroscopic level. Consequently, the method of Homogenization of Periodic Structures (HPS) is used to identify the effective permeability tensor evolution under the effect of cement filtration. The expression of the macroscopic permeability tensor derived through the HPS procedure is shown to depend on the permeating fluid viscosity and the geometrical arrangement of the sand grains and cement deposit within the microstructure. Numerical computations are made using various two-dimensional and three-dimensional microstructures, and the model results are confronted with grouting experiments performed on small scale columns in the laboratory. 相似文献
15.
The paper examines ion (chloride) transport equations in porous media (concrete) integrated over a representative elementary volume, that is to say, averaging over the macroscopic level the phenomena that occur really at the pore scale. There are three basic variables to be used: chloride concentration, moisture and temperature. The diffusion process is examined, in addition to other phenomena such as convection (the motion of dissolved substances caused by flow of water in a pore solution of partially saturated media) or chloride binding (the capacity of free chloride of being chemically bound, particularly with C3A to form Friedel salts). Contrary to other approaches, such effects are not considered by means of apparent diffusion coefficients but by developing the complete set of time‐dependent equations for both the chloride concentration within the pore solution and the moisture content within the pore space. Once the general model is described, the system of equations can be solved numerically by means of a two‐dimensional finite element formulation. The main objective is to reproduce results of experimental tests by means of a priori parameter estimation, according to the characteristics of materials and external environment conditions, thereby superseding the well‐known best fit a posteriori through Fick's second equation. While the introduction of hygrometric conditions and convection phenomena appears to be of high significance, other factors like temperature, surface concentration, chloride binding or equivalent hydration time are analysed too. The proposed model can reproduce bidimensional complex geometries, for example, cracked concrete cover, as well as variable surface condition. An application case is developed through a realistic model of the geometry of a crack. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
16.
17.
Nathan Benkemoun Hayder Al Khazraji Philippe Poullain Marta Choinska Abdelhafid Khelidj 《国际地质力学数值与分析法杂志》2018,42(3):449-468
In this paper, mesoscale hydromechanical simulations are performed to study (1) fracture features and (2) crack‐gas permeability coupling evolution in the context of the tensile splitting test. The mesostructure is based on a 2‐phase 3‐D representation of heterogeneous materials, such as concrete, where stiff aggregates are embedded into a mortar matrix. To take into account these heterogeneities without any mesh adaptation, a weak discontinuity is introduced into the strain field. In addition, a strong discontinuity is also added to take into account microcracking. This mechanical model is cast into the framework of the enhanced finite element method. Concerning the coupling with gas permeability, a double‐porosity method is used to simulate the flow through the cracks and the porosity. The apparent gas permeability is afterwards evaluated by a homogenization method. On the basis of finite element simulations, influence of aggregate size on ultimate crack opening, macroscopic ultimate tensile stress, total dissipated energy, and gas permeability evolution is numerically investigated. Furthermore, gas permeability evolution is also compared with experimental results from the literature. In addition, in the spirit of a sequential multiscale approach, macroscale gas permeability equations are identified from the hydromechanical results coming from the mesoscale computations. These equations lead to a relation between macroscale gas permeability evolution and crack opening. Besides, we show how the aggregate size influences the percolation threshold and that after this threshold, a cubic relation between macroscale gas permeability and crack opening is obtained. 相似文献
18.
19.
20.
Yazhou Zou 《Acta Geotechnica》2012,7(2):129-137
A macroscopic model for predicting the relative hydraulic permeability of unsaturated soils is proposed. In this model, pores in unsaturated soils are considered to be parallel flow tubes. The water flow in the pores is assumed to take place in the water film on the inside wall of the flow tubes. The viscosity of pore water is considered to be different from the viscosity of pure water and variable with the variation of degree of saturation. The values of tortuosity factor and pore shape factor of unsaturated soils are estimated theoretically. The theoretical model is verified using experimental data for 32 different soils. For application in engineering practice, the value of viscosity of pore water in different soils is proposed. 相似文献