共查询到20条相似文献，搜索用时 595 毫秒
1.
岩石与水相互作用的研究不仅在理论上而且在工程实践上具有重大的作用,大量的注水采油工程、核废料贮存等就是利用岩石的各向异性和其水力学特征,同时,水又时常引起岩石和土体滑坡坍塌等自然灾害。在微观分析基础上,从理论上建立了各向异性损伤岩石的能量原理及相应的各向异性Biot系数基本方程;并以红砂岩为例,运用两种不同试验方法获取各向同性Biot系数;以三轴注水、排水等试验,论证了在应力-应变全过程曲线中,Biot系数具有各向异性损伤特征。这种特性对注水和注浆等岩土工程具有重要的指导意义。 相似文献
2.
基于Leclaire对饱和双相孔隙弹性介质Biot模型的扩展,研究含有两种不同固相组分的三相多孔弹性介质中体波的传播特性。以饱和冻土为例,分析了各相体积分数、颗粒形状,接触参数等因素对波动方程中惯性参数、黏性参数、刚度参数的影响;对该三相介质模型进行了退化,分析了孔隙中只含液态水或固态冰时体波的特性;以饱和冻土为例,通过数值计算,探讨了饱和冻土中体波的相速度和衰减系数与胶结参数、接触参数、频率、饱和度、孔隙率等参数的关系。结果表明：与一般的饱和土不同,饱和冻土中存在5种体波,即3种纵波和2种横波;5种体波均具有弥散性和衰减性,且P1波、S1波弥散性和衰减性远小于P2、P3、S2波;胶结参数、饱和度、孔隙率对5种体波的传播特性影响显著,接触参数对传播特性影响较小。 相似文献
3.
Parametric study of the smooth‐joint contact model on the mechanical behavior of jointed rock 
下载免费PDF全文
下载免费PDF全文 The smooth‐joint contact model based on distinct element method has been widely used to represent discontinuity in the simulation of fractured rock mass, but there is rare efficient guidance for the selection of proper parameters of smooth‐joint contact model, which is the basement for using this model properly. In this paper, the effect of smooth joint parameters on the macroscopic properties and failure mechanism of jointed rock under triaxial compression test is investigated. The numerical results reveal that the friction coefficient of smooth joint plays a dominant role in controlling mechanical behaviors. The stiffness of smooth joint has a relative small influence on the mechanical behaviors. Poisson ratio decreases with the reduction of normal stiffness but increases with the reduction of shear stiffness. The reduction of smooth joint strength, which is determined by normal strength, cohesion, and friction angle of smooth joint, contributes to the breakage of bonded smooth joint and ultimately decreases the strength of the specimen. We proposed a detailed calibration process for smooth‐joint contact model according to the relationship between smooth‐joint parameters and mechanical properties. By following this process, the numerical results are validated against corresponding experimental results and good agreement between them can be found in stress‐strain curves and failure modes of different joint orientations. Further analyses from the microperspective are performed by looking at transmission of contact force, the nature and distribution of microcracks, and the particle displacement to show the failure process and failure modes. 相似文献
4.
Xuhui Zhang Jia Lin Xiaobing Lu Lele Liu Changling Liu Mingyao Li Yewang Su 《国际地质力学数值与分析法杂志》2018,42(7):931-942
Gas hydrate‐bearing sediments (GHBSs) have been considered as a potential energy resource. In this paper, the mechanical properties of GHBS are firstly investigated by the integrated test apparatus for synthesis of GHBS using silty sand as skeleton. Triaxial tests indicate an obvious transition of stress‐strain relationship from strain hardening under low hydrate saturation and strain softening under high hydrate saturation. The hypoplastic models coupled with Drucker‐Prager criterion and the Mohr‐Coulomb criterion are proposed to analyze the stress‐strain relationship of GHBS with considering the effective porosity because of the hydrate filling in the pores of GHBS. The strain hardening and softening behaviors are well predicted with less material parameters compared with the classical models. Compared with the test results, the proposed hypoplastic models are verified to be capable of capturing the salient features of the mechanical behaviors of GHBS under the conditions of little temperature change and no hydrate dissociation. 相似文献
5.
A hierarchical thermo‐hydro‐plastic constitutive model for unsaturated soils and its numerical implementation 下载免费PDF全文
下载免费PDF全文 《国际地质力学数值与分析法杂志》2018,42(15):1785-1805
Unsaturated soils are highly heterogeneous 3‐phase porous media. Variations of temperature, the degree of saturation, and density have dramatic impacts on the hydro‐mechanical behavior of unsaturated soils. To model all these features, we present a thermo‐hydro‐plastic model in which the hydro‐mechanical hardening and thermal softening are incorporated in a hierarchical fashion for unsaturated soils. This novel constitutive model can capture heterogeneities in density, suction, the degree of saturation, and temperature. Specifically, this constitutive model has 2 ingredients: (1) it has a “mesoscale” mechanical state variable—porosity and 3 environmental state variables—suction, the degree of saturation, and temperature; (2) both temperature and mechanical effects on water retention properties are taken into account. The return mapping algorithm is applied to implement this model at Gauss point assuming an infinitesimal strain. At each time step, the return mapping is conducted only in principal elastic strain space, assuming no return mapping in suction and temperature. The numerical results obtained by this constitutive model are compared with the experimental results. It shows that the proposed model can simulate the thermo‐hydro‐mechanical behavior of unsaturated soils with satisfaction. We also conduct shear band analysis of an unsaturated soil specimen under plane strain condition to demonstrate the impact of temperature variation on shear banding triggered by initial material heterogeneities. 相似文献
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.
高放废物地质处置库处于温度?渗流?应力（THM）多场耦合环境中,对高放废物处置库进行安全评估时,需进行多场耦合分析。然而,高放废物处置库开挖引起硐壁附近围岩应力重分布,产生损伤,导致围岩热学参数（T）、渗流参数（H）和力学参数（M）发生变化,且在空间上分布不均匀,这将会对运营期处置库THM耦合演化过程产生显著影响。通过分析高放废物处置库温度?渗流?应力三场的耦合原理和处置库围岩损伤的分布和演化规律,定义了损伤变量和损伤演化准则,并将损伤变量与热学参数、渗流参数、力学参数以及多场耦合参数（Biot系数、Biot模量和温度排水系数）建立联系,将围岩损伤与温度?渗流?应力建立联系,形成了一个弹塑性损伤温度?渗流?应力多场耦合数值模型,然后利用建立的模型对瑞士Mont Terri高放废物地质处置库围岩加热试验进行模拟,对比了模拟值和试验值,比较了考虑开挖损伤和不考虑开挖损伤对高放废物地质处置库温度?渗流?应力的影响,并分析了在多场耦合作用下开挖损伤的演化规律。 相似文献
8.
The Barcelona basic model cannot predict the mechanical behaviour of unsaturated expansive soils, whereas the Barcelona expansive model (BExM) can only predict the stress–strain behaviour of unsaturated expansive soils without the water‐retention behaviour being incorporated. Moreover, the micro‐parameters and the coupling function between micro‐structural and macro‐structural strains in the BExM are difficult to determine. Experimental data show that the compression curves for non‐expansive soils under constant suctions are shifted towards higher void ratios with increasing suction, whereas the opposite is true for expansive soils. According to the observed water‐retention behaviour of unsaturated expansive soils, the air‐entry value increases with density, and the relationship between the degree of saturation and void ratio is linear at constant suction. According to the above observation, an elastoplastic constitutive model is developed for predicting the hydraulic and mechanical behaviour of unsaturated expansive soils, based on the existing hydro‐mechanical model for non‐expansive unsaturated soil. The model takes into consideration the effect of degree of saturation on the mechanical behaviour and that of void ratio on the water‐retention behaviour. The concept of equivalent void ratio curve is introduced to distinguish the plastic potential curve from the yield curve. The model predictions are compared with the test results of an unsaturated expansive soil, including swelling tests under constant net stress, isotropic compression tests and triaxial shear tests under constant suction. The comparison indicates that the model offers great potential for quantitatively predicting the hydraulic and mechanical behaviour of unsaturated expansive soils. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
9.
A unified mechanical and retention model for saturated and unsaturated soil behaviour 总被引：1，自引：1，他引：0
The coupled mechanical and water retention elasto-plastic constitutive model of Wheeler, Sharma and Buisson (the Glasgow coupled model, GCM) predicts unique unsaturated isotropic normal compression and unsaturated critical state planar surfaces for specific volume and degree of saturation when soil states are at the intersection of mechanical (M) and wetting retention (WR) yield surfaces. Experimental results from tests performed by Sivakumar on unsaturated samples of compacted speswhite kaolin confirm the existence and form of these unique surfaces. The GCM provides consistent representation of transitions between saturated and unsaturated conditions, including the influence of retention hysteresis and the effect of plastic volumetric strains on retention behaviour, and it gives unique expressions to predict saturation and de-saturation conditions (air-exclusion and air-entry points, respectively). Mechanical behaviour is modelled consistently across these transitions, including appropriate variation of mechanical yield stress under both saturated and unsaturated conditions. The expressions defining the unsaturated isotropic normal compression planar surfaces for specific volume and degree of saturation are central to the development of a relatively straightforward methodology for determining values of all GCM parameters (soil constants and initial state) from a limited number of laboratory tests. This methodology is demonstrated by application to the experimental data of Sivakumar. Comparison of model simulations with experimental results for the full set of Sivakumar’s isotropic loading stages demonstrates that the model is able to predict accurately the variation of both specific volume and degree of saturation during isotropic stress paths under saturated and unsaturated conditions. 相似文献
10.
Empirical Correlations for Predicting Strength Properties of Rocks from P-Wave Velocity Under Different Degrees of Saturation 总被引：5，自引：3，他引：2
Determination of P-wave velocity (V p), which is closely related to intact rock properties both in laboratory and in situ conditions, is a non-destructive, easy and less complicated procedure. Due to these advantages, there is an increasing trend to predict the physico-mechanical properties of rocks from V p. By considering that no attempt on the estimation of mechanical properties of rocks from V p under different degrees of saturation has been made, in this study, it was aimed to correlate strength properties (uniaxial compressive and tensile strengths) with V p of various rock types under different degrees of saturation. For this purpose, fourteen different rock types were collected from several parts of Turkey and a comprehensive laboratory testing program was conducted. Experimental results indicated that strength and deformability properties of the rocks decreased with increasing degree of saturation, while V p showed increasing and decreasing trends depending on degree of saturation. Simple regression analysis results indicated that although prediction of the strength properties of rocks directly from V p at different degrees of saturation was possible, the equations developed would yield some under- or over-predictions. In the second stage of statistical analyses, a series of different prediction relationships were developed by using independent variables such as V p, degree of saturation and effective clay content (ECC). The statistical tests suggested that the resultant multivariate equations had very high prediction performances and were very useful tools to estimate the strength properties from V p determined at any degree of saturation. In addition, the comparisons between the theoretical Gassmann–Biot velocities, which were calculated at different degrees of saturation, and the experimental results suggest that the theoretical Gassmann–Biot velocities show inconsistencies with the experimental results obtained from the investigated rock types with high ECC. Therefore, it was concluded that the use of theoretical velocities is not suitable for rock types with high ECC. 相似文献
11.
Alexander H.D. Cheng 《国际地质力学数值与分析法杂志》2020,44(4):455-483
This paper extends the Biot theory of poroelasticity from the saturated to unsaturated case. The Biot phenomenological model uses parameters that are easily observable, such as the deformation of porous frame, total stress, pore pressure, and fluid specific discharge. Such model is preferred for engineering applications. At this macroscopic level, the extension of Biot theory from saturated to unsaturated is straightforward. The constitutive constants, however, are combined properties of solid, pore space, and fluids. In the unsaturated case, the constants are functions of the degree of saturation. Their measurements and tabulation over a range of saturation is generally not feasible for practical applications. In this work, a Biot-Willis type analysis is performed for the unsaturated case to provide a theory that the bulk material constants can be evaluated using laboratory measurable micromechanical constants under saturated condition, plus a capillary pressure curve (saturation versus suction pressure) typically available for unsaturated porous medium, without the need of measurement at each state of saturation. In particular, it is demonstrated that the surface energy contained in the meniscus interface manifests as a “capillary modulus,” given by the negative inverse slope of the capillary pressure curve. A rigorous analysis based on the thermodynamic variational energy approach is also conducted to lend theoretical support to the phenomenological approach. The presented model can bring a closure to the practical engineering modeling of the deformation of partially saturated porous medium that lacks the information of material constants over the range of saturation. 相似文献
12.
红层软岩内部孔隙具有随机、多样化的分布特点,孔隙结构的变化是影响其宏观力学性能的关键所在。由SEM扫描电镜获取岩样不同饱水时间下的细观结构图像,根据盒维数法计算出孔隙的分形维数,发现随着饱水时间的加长,孔隙的分形维数呈现增大趋势。同时对孔隙的数量、大小进行统计分析,得出不同饱水时间下岩样内部孔隙的分布特征。基于多重分形理论,采用统计矩的方法对孔隙结构进行定量表征。结果表明,孔隙结构的分配函数与q阶次趋于线性关系,验证了该结构的自相似性与无标度性,由广义分形维数D(0)>D(1)>D(2)说明了孔隙具有多重分形特征,由多重分形谱参数分析了孔隙结构的不规则性与复杂程度,更好地表征孔隙大小各异的分布情况。结合孔隙结构的多重分形特征与岩样抗压强度,建立起孔隙结构变化与其力学性能的关联性,对研究红层软岩的损伤过程具有一定的指导意义。 相似文献
13.
根据非均质岩石参数赋值方法具有随机性的特点,以随机参数的变异系数定义岩石非均质度,推导了Weibull分布参数赋值方法的非均质度计算方法,并与形状参数进行了对比分析;推导了基于矿物细胞元参数赋值方法的岩石非均质度的计算方法;利用上述2种赋值方法研究了岩石非均质度对其力学性能的影响。研究结果表明,利用变异系数定义岩石非均质度是合理的,可用于描述不同非均质岩石模型的非均质程度;非均质岩石的力学特性由细观单元力学特性和其非均质特性共同决定,岩石的非均质性对其极限强度具有弱化影响;对于不同的非均质岩石模型,可建立相同的线性函数关系表示非均质岩石极限强度受其非均质度的影响。 相似文献
14.
长期处于变低温环境中的岩石热力学效应主要表现为材料的力学特性参数、热力系数和饱和冻结状态孔隙冰胀系数等随温度改变而发生变化,这种变化对岩体的变形及强度特性均有显著影响。根据变物性参数的非线性热弹性理论,建立了考虑热力系数和冰胀力系数的本构方程,给出了单轴压缩条件下热力效应系数随温度变化的分析格式和确定方法;借助于花岗岩不同低温干燥和饱和冻结状态的单轴压缩试验资料,获得了花岗岩低温热力效应与温度的关系,并探讨了其对花岗岩抗压强度的影响特性。分析表明,在变低温环境下花岗岩热力系数和冰胀力系数均随相对温差的增大而逐渐减小,热力系数降低速率小于冰胀系数降低速率,在同一温差下热力系数大于冰胀系数;热应力和冰胀应力与温度改变量呈非线性增长关系;花岗岩抗压强度在两种状态下均呈增大趋势,但主要以热应力为主,热力系数引起的试样轴向应力增量总是大于冰胀力引起的轴向应力增量。所给方法及研究成果可用于长期处于变低温状态下材料的力学性能研究及其工程应用。 相似文献
15.
16.
Ahmad Zalooli David Martín Freire-Lista Mashalah Khamehchiyan Mohammad Reza Nikudel Rafael Fort Shahram Ghasemi 《Environmental Earth Sciences》2018,77(8):315
Durability of building stones is an important issue in sustainable development. Crystallization of soluble salts is recognized as one of the most destructive weathering agents of building stones. For this reason, durability of Ghaleh-khargushi rhyodacite and Gorid andesite from Iran was investigated against sodium sulfate crystallization aging test. Petrographic and physico-mechanical properties and pore size distribution of these stones were examined before and after the aging test. The characteristics of the microcracks were quantified with fluorescence-impregnated thin sections. Durability and physico-mechanical characteristics of Ghaleh-khargushi rhyodacite are mainly influenced by preferentially oriented preexisting microcracks. Stress induced by salt crystallization led to the widening of preexisting microcracks in Ghaleh-khargushi rhyodacite, as confirmed by the pore size distributions before and after the aging test. The preexisting microcracks of Gorid andesite were attributed to the mechanical stress induced by contraction of lava during cooling. The number of transcrystalline microcracks was significantly increased after the aging test. The degree of plagioclase microcracking was proportional to its size. Durability of the studied stones depends on initial physico-mechanical properties, pore size distribution, and orientation of microcracks. Initial effective porosity is found to be a good indicator of the stones’ durability. Salt crystallization resulted in an increase in the effective porosity with a parallel decrease in the wave velocities. Surface microroughness parameters increased with the development of salt crystallization-induced microcracking. Gorid andesite showed higher quality and durability than Ghaleh-khargushi rhyodacite. 相似文献
17.
The knowledge about thermo-mechanical properties of granite is still limited to some extent. Individual measurements are necessary to obtain reliable properties for specific granite types. A reliable numerical model of thermal cracking behaviours of granite exposed to extreme high temperatures (e.g. 800–1000 °C) is missing. In this study, the impact of temperature up to 1000 °C on physical, mechanical, and thermal properties as well as thermo-mechanical coupled behaviour of Eibenstock granite was investigated by laboratory testing and numerical simulations. The physical properties including mineral composition, density, P-wave velocity, and open porosity are measured to be temperature dependent. Uniaxial compression and Brazilian tests were carried out to measure uniaxial compressive strength (UCS), Young’s modulus, stress–strain relationship, and tensile strength of Eibenstock granite before and after thermal treatment, respectively. Thermal properties including specific heat, thermal conductivity, thermal diffusivity, and linear thermal expansion coefficient are also measured and found to be temperature dependent, especially the expansion coefficient which shows a steep increase around 573 °C as well as at 870 °C. The numerical simulation code FLAC3D was used to develop a numerical scheme to simulate the thermal-induced damage of granite at high temperatures. Statistical methods combined with real mineral composition were used to characterize the heterogeneity of granite. The numerical model is featured with reliable temperature-dependent parameters obtained from laboratory tests. It can well reproduce the laboratory results in form of thermal-induced micro- and macrocracks, as well as the stress–strain behaviour and the final failure pattern of Eibenstock granite after elevated temperatures up to 1000 °C. The simulation results also reveal that the thermal-induced microcracks are randomly distributed across the whole sample. Although most thermal-induced damages are tensile failures, shear failure begins to develop quickly after 500 °C. The obvious UCS reduction in granite due to heating is mainly caused by the increase in shear failure. The simulation also shows that the dominant impact of α–β quartz transition is widening pre-existing cracks rather than the formation of new microcracks.
相似文献18.
ABSTRACTSkempton coefficient B is commonly used to evaluate the saturation of a sample before the shearing stage in a triaxial test. For unsaturated soils, measuring the degree of saturation Sr at the end of the consolidation stage is very difficult whereas measuring Skempton coefficient B is very simple. Using this coefficient to obtain Sr of a sample in a triaxial apparatus without complex procedures might be an improvement, especially as regards testing procedures for unsaturated soils. Different authors have proposed a relationship between these two parameters based on the hypothesis that a variation of B coefficient with Sr is mostly due to the compressibility of air in the pores. However, the presence of air as a fluid phase gives rise to suction after the equilibrium is reached inside the sample leading to an increase in stiffness. Based on this observation, this article presents a new relationship that enables Sr to be evaluated with a given Skempton coefficient taking into account not only the initial void ratio and the isotropic elastic modulus but also the suction and the elastoplastic behaviour of the soil. The results of the proposed approach are then confronted with the original relationship and compared to experimental data. 相似文献
19.
Three-dimensional poromechanical back analysis of the pulse test accounting for transverse isotropy 总被引：3，自引：2，他引：1
The pulse test is usually considered as an efficient test for measuring the permeability of weakly permeable porous media. It is generally interpreted by comparison of experimental and theoretical curves obtained by a simplified 1D analytical solution. In a previous work, we proposed an improvement for test interpretation, through fully coupled numerical modeling in 2D axisymmetric configurations. The numerical modeling took into account the poromechanical coupling and exhibited 2D effects that cannot be assessed with classical interpretation. The modeling was coupled with an inverse method permitting the assessment of intrinsic permeability, mechanical parameters such as Young??s modulus, and poromechanical coupling parameters such as Biot coefficient. The interpretation of the tests on argillites showed anisotropy of the material and revealed the fundamental importance of this anisotropy for further works. The present paper aims at taking into account the transverse isotropy of the rock samples in the poromechanical back analysis of the test. A transverse isotropic poroelastic constitutive law has been developed and implemented in Code_Aster (Edf), and 3D modeling is required to account for all possible orientations of the samples. The poromechanical modeling is still coupled with an inverse method, also implemented in the finite element code. The method was applied to laboratory tests on Meuse/Haute-Marne argillites, on cylindrical samples with axis either parallel or perpendicular to the isotropy planes, and the results were compared with isotropic interpretation. The 3D transverse isotropic modeling provides more consistent values of the parameters than the 2D modeling, particularly for the mechanical parameters. 相似文献
20.
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. 相似文献