共查询到20条相似文献，搜索用时 15 毫秒
1.
A constitutive model for mechanical and chemo‐mechanical compaction in sedimentary basins and finite element analysis 
下载免费PDF全文
下载免费PDF全文 Compaction and associated fluid flow are fundamental processes in sedimentary basin deformation. Purely mechanical compaction originates mainly from pore fluid expulsion and rearrangement of solid particles during burial, while chemo‐mechanical compaction results from Intergranular Pressure‐Solution (IPS) and represents a major mechanism of deformation in sedimentary basins during diagenesis. The aim of the present contribution is to provide a comprehensive 3D framework for constitutive and numerical modeling of purely mechanical and chemo‐mechanical compaction in sedimentary basins. Extending the concepts that have been previously proposed for the modeling of purely mechanical compaction in finite poroplasticity, deformation by IPS is addressed herein by means of additional viscoplastic terms in the state equations of the porous material. The finite element model integrates the poroplastic and poroviscoplastic components of deformation at large strains. The corresponding implementation allows for numerical simulation of sediments accretion/erosion periods by progressive activation/deactivation of the gravity forces within a fictitious closed material system. Validation of the numerical approach is assessed by means of comparison with closed‐form solutions derived in the context of a simplified compaction model. The last part of the paper presents the results of numerical basin simulation performed in one dimensional setting, demonstrating the ability of the modeling to capture the main features in elastoplastic and viscoplastic compaction. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
2.
The paper deals with the modeling of some aspects, such as the formulation of constitutive equations for sediment material or finite element approach for basin analysis, related to mechanical compaction in sedimentary basins. In addition to compaction due to gravity forces and pore‐pressure dissipation, particular emphasis is given to the study of deformation induced by tectonic sequences. The numerical model relies upon the implementation of a comprehensive constitutive model for the sediment material formulated within the framework of finite poroplasticity. The theoretical model accounts for both hydromechanical and elasticity–plasticity coupling due to the effects of irreversible large strains. From the numerical viewpoint, a finite element procedure specifically devised for dealing with sedimentary basins as open systems allows to simulate within a two‐dimensional setting the process of sediment accretion or erosion. Several basin simulations are presented. The main objective is to analyze the behavior of a sedimentary basin during the different phases of its life cycle: accretion phase, pore‐pressure dissipation phase and compressive/extensional tectonic motions. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
3.
This paper is the first part of a work that aims at developing a mechanical model for the behaviour of propellant‐like materials under high confining pressure and strain rate. The behaviour of a typical material is investigated experimentally. Several microstructural deformation processes are identified and correlated with loading conditions. The resulting behaviour is complex, non‐linear, and characterized by multiple couplings. The general structure of a relevant model is sought using a thermodynamic framework. A viscoelastic‐viscoplastic‐compaction model structure is derived under suitable simplifying assumptions, in the framework of finite, though moderate, strains. Model development, identification and numerical applications are given in the companion paper. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
4.
A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
5.
6.
A new model for calculating the resistance to penetration into geological or geologically derived materials is proposed. We assume steady‐state flow of the target material over the penetrator. The target medium is described by a rate dependent constitutive equation that accounts for combined effects of strain rate and compaction on yielding. The wedge‐shaped penetrator is considered to be rigid. The influence of the characteristics of the penetrator/target interface, impact velocity, target mechanical properties and nose geometry on the resistance to penetration is investigated. It is found that for low to intermediate impact velocities, accounting for friction results in a blunter optimal wedge geometry for optimal penetration performance. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
7.
This paper completes the study presented in the accompanying paper, and demonstrates a numerical algorithm for parameter prediction from the piezocone test (CPTU) data. This part deals with a development of neural network (NN) models which are able to map multi‐variable input data onto typical geotechnical characteristics and constitutive parameters of the modified Cam clay model, which has been applied in this study. The identification procedure is designed for the coupled hydro‐mechanical boundary value problem in normally‐and lightly overconsolidated clayey soils including partially drained conditions that may appear during cone penetration. The NN models are trained with pseudo‐experimental measurements derived with the aid of the numerical model of the piezocone test, presented in the accompanying paper. Different input configurations containing CPTU measurements and some complementary data are studied with respect to the accuracy of predicted parameter values. Finally, the performance of the developed NN predictors is tested with field CPTU data which are derived from three well‐documented characterization sites, and the obtained predictions are compared with benchmark laboratory results. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
8.
9.
An adaptive substepping explicit integration scheme is developed for a porosity‐dependent hydro‐mechanical model for unsaturated soils. The model is referred to as the modified σ –Θ model in this paper, which features the employment of the subloading surface plasticity and the stress–saturation approach. On numerical aspects, convex/nonconvex subloading surfaces in the σ –Θ space may result in incorrect loading–unloading decisions during the integration. A new loading–unloading decision method is developed here to solve the problem and then embedded into the explicit integration scheme for the modified σ –Θ model. In addition, to enhance the accuracy of the explicit integration, local errors from both hydraulic and mechanical components are included in the error control for each substep. A drift correction method is also developed to ensure the state point lies on the subloading surface in the σ –Θ space within a set error level. The performance of the loading–unloading decision method for the modified σ –Θ model is discussed through comparing it with the conventional loading–unloading decision method. The importance of involving the hydraulic component in the error control is also demonstrated. The accuracy and efficiency of the proposed adaptive substepping explicit integration scheme for the modified p–Θ model are also studied via several numerical examples. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
10.
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. 相似文献
11.
This paper introduces an unconventional constitutive model for soils, which deals with a unified thermo‐mechanical modelling for unsaturated soils. The relevant temperature and suction effects are studied in light of elasto‐plasticity. A generalized effective stress framework is adopted, which includes a number of intrinsic thermo‐hydro‐mechanical connections, to represent the stress state in the soil. Two coupled constitutive aspects are used to fully describe the non‐isothermal behaviour. The mechanical constitutive part is built on the concepts of bounding surface theory and multi‐mechanism plasticity, whereas water retention characteristics are described using elasto‐plasticity to reproduce the hysteretic response and the effect of temperature and dry density on retention properties. The theoretical formulation is supported by comparisons with experimental results on two compacted clays. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
12.
Discrete element method has been widely adopted to simulate processes that are challenging to continuum-based approaches. However, its computational efficiency can be greatly compromised when large number of particles are required to model regions of less interest to researchers. Due to this, the application of DEM to boundary value problems has been limited. This paper introduces a three-dimensional discrete element–finite difference coupling method, in which the discrete–continuum interactions are modeled in local coordinate systems where the force and displacement compatibilities between the coupled subdomains are considered. The method is validated using a model dynamic compaction test on sand. The comparison between the numerical and physical test results shows that the coupling method can effectively simulate the dynamic compaction process. The responses of the DEM model show that dynamic stress propagation (compaction mechanism) and tamper penetration (bearing capacity mechanism) play very different roles in soil deformations. Under impact loading, the soil undergoes a transient weakening process induced by dynamic stress propagation, which makes the soil easier to densify under bearing capacity mechanism. The distribution of tamping energy between the two mechanisms can influence the compaction efficiency, and allocating higher compaction energy to bearing capacity mechanism could improve the efficiency of dynamic compaction.
相似文献13.
Deep geological repository involving a multibarrier system constitutes one of the most promising options for isolating high‐level radioactive waste from the human environment. To certify the efficiency of waste isolation, it is essential to understand the behaviour of confining geomaterial under a variety of environmental conditions. To this end, results from a near‐to‐real experiment, the full‐scale engineered barriers in situ experiment, are studied by means of a thermo–hydro–mechanical finite element approach, including a consistent thermoplastic constitutive model for unsaturated soils. Laboratory tests are simulated to calibrate model parameters. The results of the numerical simulations are compared with sensor measurements and show the ability of the model to reproduce the main behavioural features of the system. The influence of the hysteretic and temperature‐dependent retention of water on the mechanical response is exhibited. Finally, those results are interpreted in the light of thermoplasticity of unsaturated soils, which reveals the highly coupled and non‐linear characters of the processes encountered. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
14.
Numerical modelling of bit–rock fracture mechanisms in percussive drilling with a continuum approach
Timo Saksala 《国际地质力学数值与分析法杂志》2011,35(13):1483-1505
This paper presents a numerical method for continuum modelling of the dynamic bit–rock interaction process in percussive drilling. The method includes a constitutive model based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. The interaction between the drill bit and rock is modelled using contact mechanics by treating the bit as a rigid body. As the bit–rock interaction in percussive drilling is a transient event, the method is implemented in explicit dynamics FEM. The rock strength heterogeneity is characterized at the mesoscopic level statistically using the Weibull distribution. The bit–rock interaction is simulated under axisymmetric conditions using cylindrical and hemispherical buttons. The choice of the quite complex constitutive model accounting, e.g. for plastic compaction, viscoplastic shear and tensile failure along with induced damage and rate dependency is justified by numerical simulations. Moreover, the quasi‐static and dynamic cases are compared in plane strain simulations. Finally, some results clarifying the discrepancy of opinions found in the literature concerning the side (lateral) crack formation are obtained. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
15.
In this paper, several mechanical deformation curves of limestone are reviewed, and the effects of temperature, confining pressure, and fluid are discussed. Generally, Mohr–Coulomb is used for limestone brittle fracture. The characteristic of low temperature cataclastic flow and the conditions and constitutive equations of intracrystal plastic deformation such as dislocation creep, diffusion creep, and superplastic flow are discussed in detail. Specifically, from the macroscopic and microscopic view, inelastic compression deformation (shear-enhanced compaction) of large porosity limestone is elaborated. Compared with other mechanics models and strength equations, the dual porosity (macroporosity and microporosity) model is superior and more consistent with experimental data. Previous research has suffered from a shortage of high temperature and high pressure limestone research; we propose several suggestions to avoid this problem in the future: (1) fluid-rock interaction research; (2) mutual transition between natural conditions and laboratory research; (3) the uniform strength criterion for shear-enhanced compaction deformation; (4) test equipment; and (5) superplastic flow mechanism research. 相似文献
16.
砂岩压实作用研究主要包括压实作用机制、影响因素和响应特征三方面内容,重点阐述了砂岩压实作用模型与数值模拟研究、压实作用物理模拟与数值模拟融合研究两方面的新进展,结果表明:1）压实作用模型包括机械压实作用和化学压实作用模型,机械压实作用模型重点介绍了石英砂岩的粒间体积-深度曲线;而化学压实作用模型是与压溶作用相关的溶解作用和胶结作用,着重介绍了压溶成因的石英胶结作用。2）压实作用物理模拟与数值模拟的融合研究,主要是在压实作用机理认识基础上,通过压实模拟实验,建立基于实验数据的压实作用模型,进而开展数值模拟研究。最后指出压实作用研究趋势将朝着砂岩组分更加复杂的压实作用模型、压实作用与流体流动耦合、垂向-侧向多动力学机制、砂岩-泥岩协同机制、构造-沉积-成岩"三相"耦合的压实作用数值模拟等方面发展。 相似文献
17.
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. 相似文献
18.
Z. Liu N. Boukpeti X. Li F. Collin J.‐P. Radu T. Hueckel R. Charlier 《国际地质力学数值与分析法杂志》2005,29(9):919-940
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. 相似文献
19.
This paper presents a three‐dimensional elastoplastic constitutive model for predicting the hydraulic and mechanical behaviour of unsaturated soils. It is based on experimental results obtained from a series of controlled‐suction triaxial tests on unsaturated compacted clay with different initial densities. Hydraulic hysteresis in the water‐retention behaviour is modelled as an elastoplastic process, with the elastic part modelled by a series of scanning curves and the elastoplastic part modelled by the main drying and wetting curves. The effect of void ratio on the water‐retention behaviour is studied using data obtained from controlled‐suction wetting–drying cyclic tests on unsaturated compacted clay with different initial densities. The effect of the degree of saturation on the stress–strain‐strength behaviour and the effect of void ratio on the water‐retention behaviour are considered in the model, as is the effect of suction on the hydraulic and mechanical behaviour. The initial density dependency of the compacted soil behaviour is modelled by experimental relationships between the initial density and the corresponding yield stress and, thereby, between the initial density and the normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure and yield surfaces in the deviatoric stress plane are given by the Matsuoka–Nakai criterion. Model predictions of the stress–strain and water‐retention behaviour are compared with those obtained from triaxial tests with different initial densities under isotropic compression, triaxial compression and triaxial extension, with or without variation in suction. The comparisons indicate that the model accurately predicts the hydraulic and mechanical behaviour of unsaturated compacted soils with different initial densities using the same material constant. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
20.
Kousik Deb 《国际地质力学数值与分析法杂志》2008,32(10):1267-1288
The present study pertains to the development of a mechanical model for predicting the behavior of granular bed‐stone column‐reinforced soft ground. The granular layer that has been placed over the stone column‐reinforced soft soil has been idealized by the Pasternak shear layer. The saturated soft soil has been idealized by the Kelvin–Voigt model to represent its time‐dependent behavior and the stone columns are idealized by stiffer Winkler springs. The nonlinear behavior of the granular fill has been incorporated in this study by assuming a hyperbolic variation of shear stress with shear strain as in one reported literature. Similarly, for soft soil it has also been assumed that load‐settlement variation is hyperbolic in nature. The effect of consolidation of the soft soil due to inclusion of the stone columns has also been included in the model. Plane‐strain conditions are considered for the loading and foundation soil system. The numerical solutions are obtained by a finite difference scheme and the results are presented in a non‐dimensional form. Parametric studies for a uniformly loaded strip footing have been carried out to show the effects of various parameters on the total as well as differential settlement and stress concentration ratio. It has been observed that the presence of granular bed on the top of the stone columns helps to transfer stress from soil to stone columns and reduces maximum as well as differential settlement. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献