首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Owing to imperfect boundary conditions in laboratory soil tests and the possibility of water diffusion inside the soil specimen in undrained tests, the assumption of uniform stress/strain over the sample is not valid. This study presents a qualitative assessment of the effects of non‐uniformities in stresses and strains, as well as effects of water diffusion within the soil sample on the global results of undrained cyclic simple shear tests. The possible implications of those phenomena on the results of liquefaction strength assessment are also discussed. A state‐of‐the‐art finite element code for transient analysis of multi‐phase systems is used to compare results of the so‐called ‘element tests’ (numerical constitutive experiments assuming uniform stress/strain/pore pressure distribution throughout the sample) with results of actual simulations of undrained cyclic simple shear tests using a finite element mesh and realistic boundary conditions. The finite element simulations are performed under various conditions, covering the entire range of practical situations: (1) perfectly drained soil specimen with constant volume, (2) perfectly undrained specimen, and (3) undrained test with possibility of water diffusion within the sample. The results presented here are restricted to strain‐driven tests performed for a loose uniform fine sand with relative density Dr=40%. Effects of system compliance in undrained laboratory simple shear tests are not investigated here. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

2.
In the present paper the results of experimental and numerical investigations of hygro-thermo-mechanical properties of “Schilfsandstein” from Baden-Württemberg are presented and discussed. The aim of the study was to investigate the influence of extreme environmental conditions causing severe damage that can often be observed on various monuments during restoration and reconstruction process. In the experimental investigations the stone specimen was exposed to extreme hygro-thermal conditions, which caused significant damage (cracking). During the experiment the ultrasonic and acoustic emission analyses were performed. Subsequently, a 3D finite element analysis was carried out using a coupled hygro-thermo-mechanical model for stone. The experimental and numerical results indicate that wetting in combination with either decrease or increase of temperature lead to severe damage. The most critical locations for initialization of cracks are sharp spikes (notches), or acute angles that are inward or outward oriented (corners). Stiff but brittle layers, like late diagenetic iron cementations, might be a weak point for the initialization of crack. Furthermore, the bedding direction (anisotropy) seems to be a preferred direction for crack propagation. The results of the numerical analysis are in relatively good agreement with experimental observations.  相似文献   

3.
A finite element analysis is presented of some direct shear tests carried out on samples of a stiff, heavily overconsolidated clay. The effects of softening are accounted for through a simple material model which views this phenomenon as a consequence of an intrinsic loss of strength of the material taking place when plastic strains increase. After illustrating the main characteristics of the finite element approach, with particular reference to the correct determination of the non-reversible strains, its use in the interpretation of the mentioned shear tests is discussed. Some comments are also presented on the mesh dependence of the numerical results, and on the use of a relation between material parameters and element size which reduces this negative effect. Finally, the solution procedure is applied to the simulation of a compression test in order to investigate the possibility of numerically detecting the ‘spontaneous’ (i.e. not directly induced by the boundary conditions) formation of shear zones within the specimen.  相似文献   

4.
A new constitutive formulation for simulating the behaviour of nearly saturated sands under seismic loads is presented. The formulation is based on combining the Henry's law for dissolution of gas in water, the ideal or perfect gas law and the law of conservation of mass. The effects of transient air dissolution in water on the compressibility of partially saturated soils are also taken into account. The model was calibrated based on numerical simulations of isotropically consolidated cyclic triaxial tests conducted on partially saturated samples of Toyoura sand. A multi‐yield plasticity soil constitutive model implemented in the finite element code DYNAFLOW was used for these numerical simulations. It is shown that the formulation proposed here is able to reasonably predict the soil cyclic undrained behaviour at various degrees of saturation (95% and higher). Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

5.
Numerical modeling of stress-permeability coupling in rough fractures   总被引:2,自引:0,他引:2  
A numerical model is described for coupled flow and mechanical deformation in fractured rock. The mechanical response of rock joints to changes in hydraulic pressure is strongly influenced by the geometric characteristics of the joint surfaces. The concept of this work is to combine straightforward finite element solutions with complex and realistic fracture surface geometry in order to reproduce the non-linear stress-deformation-permeability coupling that is commonly observed in fractures. Building on the numerous studies that have expanded the understanding of the key parameters needed to describe natural rough-walled fractures, new methods have been developed to generate a finite element mesh representing discrete fractures with realistic rough surface geometries embedded in a rock matrix. The finite element code GeoSys/Rockflow was then used to simulate the coupled effects of hydraulic stress, mechanical stress, and surface geometry on the evolving permeability of a single discrete fracture. The modeling concept was experimentally verified against examples from the literature. Modeling results were also compared to a simple interpenetration model.  相似文献   

6.
贾善坡  吴渤  陈卫忠  伍国军  高敏  龚俊 《岩土力学》2014,35(8):2375-2384
以热力学和弹塑性力学理论为基础,分析岩石热-力完全耦合作用及其对力学参数和热特性参数的影响,建立了岩石热-力-损伤耦合模型及其参数演化方程,以ABAQUS软件为平台对其进行二次开发,并通过典型算例验证了岩石热-力完全耦合的重要性。然后以某深埋软岩隧洞为例,研究温度和开挖卸载共同作用下的隧洞围岩力学行为和损伤过程。计算结果表明:温度对岩石的力学性质和损伤演化过程影响显著,开挖损伤和热应力诱发的损伤对围岩热力学参数的影响不可忽略;所提出的力学模型可以有效反映围岩损伤演化、调热圈演化以及热力学参数演化,具有一定的借鉴作用。  相似文献   

7.
A numerical model is proposed to analyse elastic as well as elastoplastic behaviour of stone-column reinforced foundations. The stone-columns are assumed to be dispersed within the in situ soil and a homogenization technique is invoked to establish equivalent material properties for in situ soil and stone-column composite. The difficulties encountered in carrying out elastoplastic analyses of composite materials are overcome by adopting a separate yield function for each of the constituent materials and a sub-iteration procedure within an implicit backward Euler stress integration scheme. In the proposed procedure, equilibrium as well as kinematic conditions implied in the homogenization procedure are satisfied for both elastic as well as elastoplastic stress states. The proposed model is implemented in an axi-symmetric finite element code and numerical prediction is made for the behaviour of model circular footings resting on stone-column reinforced foundations. This prediction indicates good agreement with experimental observation. Finally, a new scheme in which the length of stone-column is variable is proposed and its behaviour is examined through a numerical example. © 1998 John Wiley & Sons, Ltd.  相似文献   

8.
《Computers and Geotechnics》2006,33(6-7):316-329
The aim of the paper is to provide new elements concerning the constitutive behavior of sedimentary rocks and the numerical aspects for basin simulators. A comprehensive model for mechanical compaction of sedimentary basins is developed within finite poroplasticity setting. Particular concern is paid to the effects of large porosity changes on the poromechanical properties of the sediment material. A simplified micromechanics-based approach is used to account for the stiffness increase and hardening induced by large plastic strains.A key challenge for numerical assessment of sedimentary basin evolution is to integrate multiple coupled processes in the context of open material systems. To this end, a numerical approach inspired from the ‘deactivation/reactivation’ method used for the simulation of excavation process and lining placement in tunnel engineering, has been developed. Periods of sediments accretion are simulated by progressive activation of the gravity forces within a fictitious closed system. Fundamental components of the constitutive model developed before (hydromechanical coupling, dependence of poroelastic properties on large plasticity, impact of irreversible porosity changes on the hardening rule, evolution of permeability with porosity) are included into our finite element code.Illustrative examples of basin simulation are performed in the one-dimensional case. Various aspects of the constitutive model are investigated. Their influence on the corresponding basin response is analyzed in terms of compaction law, porosity and fluid pressure profiles.  相似文献   

9.
张翠莲 《岩土力学》2016,37(9):2721-2727
基于连续介质损伤力学框架,通过损伤张量和有效应力来描述节理岩体的力学性能,自主研发了基于损伤力学模型的有限元程序(CD-FEM),用于节理岩体等效力学性能研究。同时,采用Karhunen-Loeve(K-L)展开来分解随机输入场,用混沌多项式来表示随机输出场,采用概率配点法生成配点,再由连续损伤有限元分析系统CD-FEM求解确定性方程组,最终得到输出域的统计数据,从而提出了一种将随机分析与基于连续损伤力学模型的数值分析方法解耦进行节理岩体不确定性力学行为分析的方法。利用该方法,对一典型节理岩体在加载条件下的力学行为进行不确定性分析,并与蒙特卡罗方法进行对比,结果表明,所提方法的计算量大大减少,极大地提高了节理岩体力学性能不确定性分析的效率,可应用于对节理岩体在不同载荷条件下的不确定性进行分析。  相似文献   

10.
Many numerical methods that simulate groundwater flow, particularly the continuous Galerkin finite element method, do not produce velocity information directly. Many algorithms have been proposed to improve the accuracy of velocity fields computed from hydraulic potentials. The differences in the streamlines generated from velocity fields obtained using different algorithms are presented in this report. The superconvergence method employed by FEFLOW, a popular commercial code, and some dual-mesh methods proposed in recent years are selected for comparison. The applications to depict hydrogeologic conditions using streamlines are used, and errors in streamlines are shown to lead to notable errors in boundary conditions, the locations of material interfaces, fluxes and conductivities. Furthermore, the effects of the procedures used in these two types of methods, including velocity integration and local conservation, are analyzed. The method of interpolating velocities across edges using fluxes is shown to be able to eliminate errors associated with refraction points that are not located along material interfaces and streamline ends at no-flow boundaries. Local conservation is shown to be a crucial property of velocity fields and can result in more accurate streamline densities. A case study involving both three-dimensional and two-dimensional cross-sectional models of a coal mine in Inner Mongolia, China, are used to support the conclusions presented.  相似文献   

11.
This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

12.
为了揭示循环扰动和高温影响下砂岩的岩爆机制,开展了循环扰动和高温作用下砂岩单轴压缩试验和CT扫描试验,研究了不同条件下砂岩的力学特性、岩爆倾向性及破坏特征,探讨了岩爆倾向性与破坏特征的关系。研究结果表明:循环扰动和高温对砂岩的力学性能及岩爆倾向性影响效果显著;无循环扰动砂岩的单轴抗压强度、弹性模量及岩爆倾向性随温度的升高呈先增加后降低的趋势,200℃为该类砂岩的阈值温度,受循环扰动砂岩的力学特性及岩爆倾向性随温度的上升而降低;而砂岩的力学特性及岩爆倾向性随循环应力幅值增加而下降;随着循环应力幅值和温度的增加,砂岩破坏模式由劈裂破坏向剪切破坏转变,同时砂岩的岩爆倾向性与裂隙三维分形维数呈良好的负相关关系;此外,高温对砂岩的力学性能、岩爆倾向性及破坏程度影响效果强于循环扰动。研究结果可为高温工程岩爆防治提供理论依据和工程参考价值。  相似文献   

13.
This paper presents parallel and serial viscoelasto‐plastic models to simulate the rate‐independent and the rate‐dependent permanent deformation of stone‐based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto‐plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement‐based incremental formulation for the viscoelastic part and an elastic predictor—plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto‐plastic model employed a viscoelastic predictor—plastic corrector algorithm. The stone‐based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto‐plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto‐plastic matrix models. The simulation results showed that the serial viscoelasto‐plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto‐plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate‐independent and rate‐dependent viscoelasto‐plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto‐plastic behaviors of the stone‐based materials. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

14.
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.  相似文献   

15.
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.  相似文献   

16.
A procedure for numerical approximation to two-dimensional, hydraulically-driven fracture propagation in a poroelastic material is described. The method uses a partitioned solution procedurè to solve a finite element approximation to problems described by the theory of poroelasticity, in conjunction with a finite difference approximation for modelling fluid flow along the fracture. An equilibrium fracture model based on a generalized, Dugdale–Barenblatt concept is used to determine the fracture dimensions. An important feature is that the fracture length is a natural product of the solution algorithm. Two example problems verify the accuracy of the numerical procedure and a third example illustrates a fully-coupled simulation of fracture propagation. Photographs taken from a high-performance engineering workstation provide insight into the nature of the coupling among the physical phenomena.  相似文献   

17.
In recent years, a number of constitutive models have been proposed to describe mathematically the mechanical response of natural clays. Some of these models are characterized by complex formulations, often leading to non‐trivial problems in their numerical integration in finite elements codes. The paper describes a fully implicit stress‐point algorithm for the numerical integration of a single‐surface mixed isotropic–kinematic hardening plasticity model for structured clays. The formulation of the model stems from a compromise between its capability of reproducing the larger number of features characterizing the behaviour of structured clays and the possibility of developing a robust integration algorithm for its implementation in a finite elements code. The model is characterized by an ellipsoid‐shaped yield function, inside which a stress‐dependent reversible stiffness is accounted for by a non‐linear hyperelastic formulation. The isotropic part of the hardening law extends the standard Cam‐Clay one to include plastic strain‐driven softening due to bond degradation, while the kinematic hardening part controls the evolution of the position of the yield surface in the stress space. The proposed algorithm allows the consistent linearization of the constitutive equations guaranteeing the quadratic rate of asymptotic convergence in the global‐level Newton–Raphson iterative procedure. The accuracy and the convergence properties of the proposed algorithm are evaluated with reference to the numerical simulations of single element tests and the analysis of a typical geotechnical boundary value problem. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

18.
Monowire block cutting machines can be used for natural stone block squaring and slab cutting operations. The plants where the cutting operations are performed demand high product quality with minimum operational costs. The major parameters affecting the economy of the operation are the energy consumed and the wear induced on the diamond beads during the cutting operation. An efficient cutting operation can only be maintained by selecting proper cutting parameters. Therefore, cutting parameters should be clearly understood. Experimental studies and numerical modeling methods are significant in terms of identifying the energy consumption occurring during natural stone cutting with monowire. Experimental studies and numerical modeling using discrete element method were performed on Afyon White Marble. Experimental studies have been performed by using a specially designed, fully automatic monowire cutting machine, and numerical analyses were carried out by commercially available software called three-dimensional particle flow code (PFC3D). A discrete element model for the cutting operation was developed, and various numerical models were performed for different peripheral speeds and cutting speeds, while, at the same time, the actual cutting operations were being carried out in the laboratory. Finally, the data obtained from the experimental works were compared with the data from numerical modeling. A comparison indicates that the frictional energy values obtained by means of numerical modeling are in good agreement with the results of the laboratory measurements. This study clearly put forward the influences of effective parameters on monowire cutting operations in natural stone industry. Furthermore, it filled an important space in the literature about the use of monowire block squaring machines.  相似文献   

19.
The use of geosynthetic-encased stone columns as a method for soft soil treatment is extensively used to increase the bearing capacity and reduce the settlement of raft foundations and the foundation of structures like embankments. Pre-strain is an effect occurring in the encasement during stone column installation due to the compaction of the stone material. The present study uses the finite element program Plaxis to perform a numerical analysis of the soft clay bed reinforced by geosynthetic-encased stone columns. An idealization is proposed for simulation of installation of geosynthetic-encased stone columns in soft clay based on the unit-cell concept. In the analyses, initially, the validity of the analysis of the single column-reinforced soil in the unit-cell model was performed through comparison with the group columns. Then, by considering a unit-cell model, the finite element analyses were carried out to evaluate the stiffness of the reinforced ground to estimate the settlement. The results of the analyses show that the improved stiffness of the encased stone column is not only due to the confining pressure offered by the geosynthetic after loading, but the initial strain of the geosynthetic that occurred during installation also contributes to the enhancement of the stiffness of the stone column and the reduction of the settlement.  相似文献   

20.
To reduce the weathering rate of natural building stones, a wide variety of water repellents and consolidants are commercially available. Although a lot of research is performed on these products, it remains difficult to determine which product is appropriate to use for a certain type of building stone. Each type of building stone has its own petrographic characteristics (mineralogy, texture, …), leading to typical physical and technical properties which influence its rate of decay. The localisation of the products inside a stone type is not only depending on the properties of the products themselves, but also linked to the texture and structure of the stone. The impregnation depth of the products strongly influences their efficiency and is therefore a key issue in the determination if a product is functional for a certain type of stone. X-ray micro-CT has recently been introduced as a non-destructive material evaluation technique for engineering and geology purposes. The fact that micro-CT can provide information about the internal structure and properties of natural building stones, is a major advantage in the study of their conservation. Firmly linked with more classical research techniques, this non-destructive technique offers an extra dimension to the cultural heritage research.In this paper, non-destructive X-ray computed micro-tomography (micro-CT) turns out to be a powerful tool as it can visualise the presence of water repellents and consolidants inside the stone and can help to detect the influence that these treatments exert on porosity and pore-size distribution. For the visualization of the product its atomic number and density and the amount of product inside the natural building stone is crucial. Besides the contrast in attenuation, the resolution of the micro-CT also needs to be taken into account to obtain a good contrast between stone and product. By doping with 3-bromopropyltrimethoxysilane, more contrast occurs between the stone material and the conservation product. When the amount of product inside the samples is very low and the product is homogeneously spread inside the sample, the doped product is difficult to distinguish from the stone material due to the resolution of the micro-CT. When the amount of doped product inside the stone material reaches a certain threshold value, the dope will appear in the 2D cross-sections derived from micro-CT, creating a good visualization of the products inside the stone samples.  相似文献   

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

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