首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 272 毫秒
1.
This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

2.
A numerical model is presented to describe the evolution of fracture aperture (and related permeability) mediated by the competing chemical processes of pressure solution and free‐face dissolution/precipitation; pressure (dis)solution and precipitation effect net‐reduction in aperture and free‐face dissolution effects net‐increase. These processes are incorporated to examine coupled thermo‐hydro‐mechano‐chemo responses during a flow‐through experiment, and applied to reckon the effect of forced fluid injection within rock fractures at geothermal and petroleum sites. The model accommodates advection‐dominant transport systems by employing the Lagrangian–Eulerian method. This enables changes in aperture and solute concentration within a fracture to be followed with time for arbitrary driving effective stresses, fluid and rock temperatures, and fluid flow rates. This allows a systematic evaluation of evolving linked mechanical and chemical processes. Changes in fracture aperture and solute concentration tracked within a well‐constrained flow‐through test completed on a natural fracture in novaculite (Earth Planet. Sci. Lett. 2006, in press) are compared with the distributed parameter model. These results show relatively good agreement, excepting an enigmatic abrupt reduction in fracture aperture in the early experimental period, suggesting that other mechanisms such as mechanical creep and clogging induced by unanticipated local precipitation need to be quantified and incorporated. The model is applied to examine the evolution in fracture permeability for different inlet conditions, including localized (rather than distributed) injection. Predictions show the evolution of preferential flow paths driven by dissolution, and also define the sense of permeability evolution at field scale. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

3.
Material behaviour that exhibits characteristics of creep induced by a spontaneous mineral dissolution enhanced by material damage is studied. It is believed that the characteristic rates of the chemical processes involved determine the time‐rate dependence of the resulting strain. A basic model of a combined chemo‐plastic softening and chemically enhanced deviatoric strain hardening for saturated geomaterials is presented. Chemical softening is postulated to occur as a consequence of the net mass removal resulting from dissolution and precipitation of specific minerals occurring at the damage‐generated inter‐phase interfaces. Closed and open systems are discussed. In the former case, deformation at constant stress results entirely from a local compensation mechanism between the chemical softening and strain hardening. The classical three stages of creep are interpreted in terms of mechanisms of dissolution and precipitation, as well as the variation in the reaction surface areas involved in the mass exchange. In an open system, the above local mechanism is enhanced by the removal of mass via diffusion of species affecting the mass balance. Such a system is addressed via a boundary value problem as shown in an example. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

4.
In this paper, we used a theoretical model for the variation of Eulerian porosity, which takes into account the adsorption process known to be the main mechanism of production or sequestration of gas in many reservoir of coal. This process is classically modeled using Langmuir's isotherm. After implementation in Code_Aster, a fully coupled thermo‐hydro‐mechanical analysis code for structures calculations, we used numerical simulations to investigate the influence of coal's hydro‐mechanical properties (Biot's coefficient, bulk modulus), Langmuir's adsorption parameters, and the initial liquid pressure in rock mass during CO2 injection in coal. These simulations showed that the increase in the values of Langmuir's parameters and Biot's coefficient promotes a reduction in porosity because of the adsorption process when the gas pressure increases. Low values of bulk modulus increase the positive effect (i.e., increase) of hydro‐mechanical coupling on the porosity evolution. The presence of high initial liquid pressure in the rock mass prevents the progression of injected gas pressure when CO2 dissolution in water is taken into account. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
The volumetric compaction due to wetting processes is a phenomenon observed quite often in unsaturated soils. Under certain circumstances, saturation events can result into a sudden and unexpected collapse of the system. These phenomena are usually referred to as wetting‐induced collapses, without providing any detailed theoretical justification for this terminology. In order to predict in a general fashion the occurrence of coupled instabilities induced by saturation processes, a generalization of the theoretical approaches usually employed for saturated geomaterials is here provided. More specifically, this paper addresses the problem of hydro‐mechanical instability in unsaturated soils from an energy standpoint. For this purpose, an extension of the definition of the second‐order work is here suggested for the case of unsaturated porous media. On the basis of some examples of numerical simulations of laboratory tests, coupled hydro‐mechanical instabilities are then interpreted in the light of this second‐order energy measure. Finally, the implications of the theoretical results here presented are commented from a constitutive modelling perspective. Two possible alternative approaches to formulate incremental coupled constitutive relations are indeed discussed, showing how the onset of hydro‐mechanical instabilities can be predicted using an extended form of Hill's stability criterion. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

6.
This paper discusses a series of stress point algorithms for a breakage model for unsaturated granular soils. Such model is characterized by highly nonlinear coupling terms introduced by breakage‐dependent hydro‐mechanical energy potentials. To integrate accurately and efficiently its constitutive equations, specific algorithms have been formulated using a backward Euler scheme. In particular, because implementation and verification of unsaturated soil models often require the use of mixed controls, the incorporation of various hydro‐mechanical conditions has been tackled. First, it is shown that the degree of saturation can be replaced with suction in the constitutive equations through a partial Legendre transformation of the energy potentials, thus changing the thermomechanical state variables and enabling a straightforward implementation of a different control mode. Then, to accommodate more complex control scenarios without redefining the energy potentials, a hybrid strategy has been used, combining the return mapping scheme with linearized constraints. It is shown that this linearization strategy guarantees similar levels of accuracy compared with a conventional strain–suction‐controlled implicit integration. In addition, it is shown that the use of linearized constraints offers the possibility to use the same framework to integrate a variety of control conditions (e.g., net stress and/or water‐content control). The convergence profiles indicate that both schemes preserve the advantages of implicit integration, that is, asymptotic quadratic convergence and unconditional stability. Finally, the performance of the two implicit schemes has been compared with that of an explicit algorithm with automatic sub‐stepping and error control, showing that for the selected breakage model, implicit integration leads to a significant reduction of the computational cost. Such features support the use of the proposed hybrid scheme also in other modeling contexts, especially when strongly nonlinear models have to be implemented and/or validated by using non‐standard hydro‐mechanical control conditions. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

7.
This paper presents an upscaling concept of swelling/shrinking processes of a compacted bentonite/sand mixture, which also applies to swelling of porous media in general. A constitutive approach for highly compacted bentonite/sand mixture is developed accordingly. The concept is based on the diffuse double layer theory and connects microstructural properties of the bentonite as well as chemical properties of the pore fluid with swelling potential. Main factors influencing the swelling potential of bentonite, i.e. variation of water content, dry density, chemical composition of pore fluid, as well as the microstructures and the amount of swelling minerals are taken into account. According to the proposed model, porosity is divided into interparticle and interlayer porosity. Swelling is the potential of interlayer porosity increase, which reveals itself as volume change in the case of free expansion, or turns to be swelling pressure in the case of constrained swelling. The constitutive equations for swelling/shrinking are implemented in the software GeoSys/RockFlow as a new chemo‐hydro‐mechanical model, which is able to simulate isothermal multiphase flow in bentonite. Details of the mathematical and numerical multiphase flow formulations, as well as the code implementation are described. The proposed model is verified using experimental data of tests on a highly compacted bentonite/sand mixture. Comparison of the 1D modelling results with the experimental data evidences the capability of the proposed model to satisfactorily predict free swelling of the material under investigation. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

8.
Damage induced by microcracking affects not only the mechanical behaviour of geomaterials but also their hydraulic properties. Evaluating these impacts is important for many engineering applications, such as the safety assessment of radioactive waste disposal facilities. This paper presents a new constitutive model accounting simultaneously for the impact of damage on hydraulic and mechanical properties of unsaturated poroplastic geomaterials. The hydro‐mechanical coupling is formulated by means of the thermodynamic framework for partially saturated media, extended by taking into account isotropic damage and plasticity. State and complementary laws are governed by the so‐called plastic effective stress and equivalent pore pressure. Assuming a bimodal pore size distribution for cracked porous media, the hydraulic part (water retention curve and hydraulic conductivity) is modelled using phenomenological functions of damage variable. The participation of damage on both mechanical and hydraulic part enables this model to describe bilateral couplings between them. This coupled model is then validated against a number of experimental data obtained from Callovo‐Oxfordian argillite, which is the possible host rock for a radioactive waste disposal in France. Parametric studies are also carried out to check the consistency and to better demonstrate the bilateral couplings in the model. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
The use of the asymptotic limit can greatly simplify the theoretical analysis of chemical dissolution front instabilities in fluid‐saturated rocks and therefore make it possible to obtain mathematical solutions, which often play a crucial role in understanding the propagation behavior of chemical dissolution fronts in chemical dissolution systems. However, there has been a debate in recent years that the asymptotic limit of the acid dissolution capacity (i.e., the acid dissolution capacity number approaching zero) alone cannot lead to a sharp dissolution front of the Stefan type in the acidization dissolution system, in which the dissolvable minerals of carbonate rocks are chemically dissolved by the injected acid flow. The acid dissolution capacity number is commonly defined as the ratio of the volume of the carbonate rock dissolved by an acid to that of the acid. In this paper, we use four different proof methods, including (i) direct use of the fundamental concepts; (ii) use of the mathematical governing equations of an acidization dissolution system; (iii) use of the different time scaling approach; and (iv) use of a moving coordinate system approach, to demonstrate that the asymptotic limit of the acid dissolution capacity can indeed lead to sharp dissolution fronts of the Stefan type in acidization dissolution systems on a much larger time scale (than the dissolution time scale). Our new finding is that on the reaction time scale, the condition of the conventional time derivative of porosity approaching zero alone can ensure that the acidization dissolution front has a sharp shape of the Stefan type. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

10.
Li  Zhenze  Su  Grant  Zheng  Quinn  Nguyen  Thanh Son 《Acta Geotechnica》2020,15(3):635-653

Significant chemical influence on the swelling potential of MX-80 bentonite was observed during swelling tests where specimens were hydrated with highly concentrated brine. The maximum swelling pressure for specimens hydrated with brine was about 30% of the maximum swelling pressure for the same specimens hydrated with de-ionized water. The maximum swelling pressure was attained within tens of hours of brine infiltration and further decreased by half within a year. A fully coupled hydro–mechanical–chemical (HMC) dual-porosity model is proposed in this paper to interpret the swelling behaviour of MX-80 when infiltrated with brine. The dependence of hydraulic and mechanical properties on such factors as porosity, salinity and water content was investigated. A nonlinear elastic constitutive model was proposed to correlate the swelling pressure with the variation in the microporosity. The chemical effects on the mechanical behaviour were coupled at the micropore level. A number of relationships have been developed for MX-80, i.e. micropore permeability as a function of void ratio, water retention characteristics of micropores and macropores, micropore dependence on water content and the diffusion coefficients of the two types of pore structure. The proposed model was successful in reproducing both quantitatively and qualitatively the experimental results from two sets of infiltration experiments on compacted MX-80 bentonite.

  相似文献   

11.
The dissolution of crinoid columnals during tectonic deformation of the Appalachian Plateau was enhanced by stress-induced changes of chemical potential of calcite in solution at the surface of the crinoid columnal. Pressure solution on the outside surface occurred within areas of highest normal stress developed at grain to grain contacts whereas free-face dissolution on the inside bore (the axial canal) of the columnal occurred where parts of the crystal lattice next to the surface were subject to the highest strain, as indicated by the presence of mechanical twins. For free-face dissolution, the chemical potential of the solute is affected more by the larger strain energy associated with mechanical twinning than strain energy from elastic strain prior to mechanical twinning of the crystal lattice. These observations suggest that free-face dissolution as well as pressure solution may contribute significantly to dissolution during rock deformation by diffusive mass transfer.  相似文献   

12.
水对名义无水矿物变形的影响   总被引：1，自引：0，他引：1  
在固体地球中,水虽微量,但对众多地质过程(例如,岩石部分熔融与火山喷发、地震活动等)和岩石的物理化学性质(例如,电导率、滞弹性、地震波性质、相变动力学等)影响重大。更为重要的是,水能通过影响矿物的变形机制来控制岩石的流变强度,进而制约着地球动力学的过程。名义无水矿物(NAMs:Nominal anhydrous minerals)即为分子式中不含氢的矿物,其晶格的容水量远小于正常含水矿物(如,角闪石,蛇纹石等)的容水量,但由于NAMs在固体地球中体积比例甚大,仅上地幔的橄榄石中所能溶解的水可能比全部地表水还多。因此了解水对NAMs(尤其是分别作为地壳和上地幔主要组成矿物的石英和橄榄石)变形的影响对于精确地构建岩石圈强度剖面和深刻理解构造地质学与地球动力学过程至关重要。本文将系统地回顾水对NAMs变形的影响,首先通过回顾水在固体地球内部的存在形式提出了NAMs是固体地球中的重要水库,接着阐述了NAMs中水的存在形式、溶解机制、溶解度影响因素及扩散动力学,最后着重论证了水致弱化在石英和富镁石榴石中最强,然后依次是单斜辉石、长石、橄榄石,瓦德利石和林伍徳石。  相似文献   

13.
Since low permeability rock has the potential to behave as an actual semi-permeable osmotic membrane, chemical osmosis may have an effect on water flow and needs to be accounted for in analysis. Hence, an advanced fully coupled formulation for modelling unsaturated hydro, chemical and mechanical effects is here presented based on modified mixture theory. Firstly, the classical Darcy’s law has been modified to include the chemical potential. Secondly, the mechanical deformation, which affects the permeability, has been coupled with the water and chemical flows by using Helmholtz free energy and the Gibbs–Duhem equation. Finally, an illustrative numerical example with relevance to radioactive waste disposal has been analysed to investigate the performance of the coupled formulation and the relative influence of parameters.  相似文献   

14.
The dependence of rock behavior on the deformation rate is still not well understood. In salt rock, the fundamental mechanisms that drive the accumulation of irreversible deformation, the reduction of stiffness, and the development of hysteresis during cyclic loading are usually attributed to intracrystalline plasticity and diffusion. We hypothesize that at low pressure and low temperature, the rate‐dependent behavior of salt rock is governed by water‐assisted diffusion along grain boundaries. Accordingly, a chemo‐mechanical homogenization framework is proposed in which the representative elementary volume (REV) is viewed as a homogeneous polycrystalline matrix that contains sliding grain‐boundary cracks. The slip is related to the mass of salt ions that diffuse along the crack surface. The relationship between fluid inclusion‐scale and REV‐scale stresses and strains is established by using the Mori–Tanaka homogenization scheme. It is noted from the model that a lower strain rate and a larger number of sliding cracks enhance stiffness reduction and hysteresis. Thinner sliding cracks (i.e., thinner brine films) promote stiffness reduction and accelerate stress redistributions. The larger the volume fraction of the crack inclusions, the larger the REV deformation and the larger the hysteresis. Results presented in this study shed light on the mechanical behavior of salt rock that is pertinent to the design of geological storage facilities that undergo cyclic unloading, which could help optimize the energy production cycle with low carbon emissions.  相似文献   

15.
Effective capabilities of combined chemo‐elasto‐plastic and unsaturated soil models to simulate chemo‐hydro‐mechanical (CHM) behaviour of clays are examined in numerical simulations through selected boundary value problems. The objective is to investigate the feasibility of approaching such complex material behaviour numerically by combining two existing models. The chemo‐mechanical effects are described using the concept of chemical softening consisting of reduction of the pre‐consolidation pressure proposed originally by Hueckel (Can. Geotech. J. 1992; 29 :1071–1086; Int. J. Numer. Anal. Methods Geomech. 1997; 21 :43–72). An additional chemical softening mechanism is considered, consisting in a decrease of cohesion with an increase in contaminant concentration. The influence of partial saturation on the constitutive behaviour is modelled following Barcelona basic model (BBM) formulation (Géotech. 1990; 40 (3):405–430; Can. Geotech. J. 1992; 29 :1013–1032). The equilibrium equations combined with the CHM constitutive relations, and the governing equations for flow of fluids and contaminant transport, are solved numerically using finite element. The emphasis is laid on understanding the role that the individual chemical effects such as chemo‐elastic swelling, or chemo‐plastic consolidation, or finally, chemical loss of cohesion have in the overall response of the soil mass. The numerical problems analysed concern the chemical effects in response to wetting of a clay specimen with an organic liquid in rigid wall consolidometer, during biaxial loading up to failure, and in response to fresh water influx during tunnel excavation in swelling clay. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

16.
Halite single crystals in saturated solution were used to study dissolution precipitation creep (DPC) at conditions where plastic deformation is negligible. Specifically, the free unloaded surfaces of these crystals were investigated by a novel Linnik-based phase shift interference microscope. The method allows observations of the crystal surface in-situ and with an axial resolution in the nanometer scale. Transport phenomena in open systems, temperature gradients, and gradients in strain energy density were found to cause morphological changes on the free crystal surface by dissolution/reprecipitation. We did not find evidence for DPC by applying a homogeneous stress field to the crystal as long as plastic deformation was avoided. These findings suggest that deformation of rocks by DPC in situations where dislocation creep is not activated, but is rather promoted by fluid transport through the rock or by episodic changes of extensive parameters affecting solubility than by homogeneous stress alone.Editorial responsibility: J. Hoefs  相似文献   

17.
A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

18.
In western Vanoise (French Alps), karstic pockets of Triassic‐to‐Jurassic metabauxites embedded in carbonate rocks and containing several generations of metamorphic veins were studied. During blueschist facies metamorphism, a cumulative amount of ～13 vol% of water is inferred to have been produced locally by successive dehydration reactions, and part of this fluid remained in the bauxitic lenses during most of the metamorphic cycle. Field and geochemical evidence show that these rocks have been isolated from large‐scale fluid flow (closed‐system behaviour). It is proposed that the internally derived fluid has promoted the opening of fluid‐filled open spaces (as attested by the euhedral habits of vein minerals) and served as medium for mass transfer from rock to vein. Indeed, the vein infill is obviously the result of chemical interactions, at the millimetre‐to‐centimetre scale, between the rock minerals and the locally produced aqueous fluid. Two vein types can be distinguished based on mineralogical and textural features: (i) some veins are filled with newly formed products of either prograde (chloritoid) or retrograde (chlorite) metamorphic reactions; in this case, fluid‐filled open spaces seem to offer energetically favourable nucleation/growth sites; (ii) the second vein type is infilled with cookeite or pyrophyllite, that were present in the host rock prior to the vein formation. In this closed chemical system, the components for the vein infill minerals have been transferred from rock to vein through the fluid, in a dissolution–transport–precipitation process, possibly stress‐assisted. These different vein generations all contain Al‐rich mineral infills, suggesting that Al was a mobile element (cm scale) during metamorphism. In these HP rocks, fluid flow may have been restricted, and if so mass transfer occurred by diffusion in an almost stagnant fluid. Metamorphic veins can be seen as witnesses of fluid and mass redistribution that partly accommodate the rock deformation (lateral segregation).  相似文献   

19.
A significant increase of the permeability of concrete upon micro‐cracking and a good correlation between the evolution of damage (material stiffness) and permeability are observed experimentally. The present contribution investigates this correlation theoretically, with the help of lattice analyses. Scaling analysis of lattices which contain elastic brittle bonds has shown that the material degradation should be described by the evolution of the material stiffness, or compliance, in a continuum setting (damage models). This result is reviewed and further documented in the first part of the paper. In the second part, hydro‐mechanical problems are considered with the construction of a hydraulic lattice, dual to the mechanical one. We observe that the average permeability upon micro‐cracking is the lattice scale‐independent controlling variable in the hydraulic problem. Additionally, results show that in a continuum poro‐mechanical approach, the evolution of the material permeability ought to be related to the elastic unloading stiffness, described e.g. with the help of continuum damage variables. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

20.
In this paper a new analytical model is proposed to determine the permeability tensor for fractured rock masses based on the superposition principle of liquid dissipation energy. This model relies on the geometrical characteristics of rock fractures and the corresponding fracture network, and demonstrates the coupling effect between fluid flow and stress/deformation. This model empirically considers the effect of pre‐peak shear dilation and shear contraction on the hydraulic behavior of rock fractures and can be used to determine the applicability of the continuum approach to hydro‐mechanical coupling analysis. Results of numerical analysis presented in this paper show that the new model can effectively describe the permeability of fractured rock masses, and can be applied to the coupling analysis of seepage and stress fields. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号