首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
In a companion Part I of this paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.735 ), a coupled hydro‐mechanical (HM) formulation for geomaterials with discontinuities based on the finite element method (FEM) with double‐node, zero‐thickness interface elements was developed and presented. This Part II paper includes the numerical solution of basic practical problems using both the staggered and the fully coupled approaches. A first group of simulations, based on the classical consolidation problem with an added vertical discontinuity, is used to compare both the approaches in terms of accuracy and convergence. The monolithic or fully coupled scheme is also used in an application example studying the influence of a horizontal joint in the performance of a reservoir subject to fluid extraction. Results include a comparison with other numerical solutions from the literature and a sensitivity analysis of the mechanical parameters of the discontinuity. Some simulations are also run using both a full non‐symmetric and a simplified symmetric Jacobian matrix. On top of verifying the model developed and its capability to reflect the conductivity changes of the interface with aperture changes, the results presented also lead to interesting observations of the numerical performance of the methods implemented. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

2.
An advanced elasto‐plastic constitutive model for frictional materials, whose incremental version is presented in a companion paper (Int. J. Numer. Anal. Meth. Geomech., 2002; 26 :647), is implemented in a user‐defined material module. The general calculation strategy inside this module is presented and discussed, including the initial intersection of the yield surface and the techniques for updating of stresses and hardening modulus. Several integration schemes are implemented in the module and their capabilities in relation to the advanced, three‐dimensional constitutive model are evaluated. The forward Euler, modified Euler, and Runge–Kutta–Dormand–Prince integration schemes are explained in detail, compared, and evaluated in view of error tolerances and computational efficiency. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

3.
The paper describes an approach to the simulation of arbitrary laboratory tests with homogeneous stress and strain fields applicable to most constitutive models common in geomechanics. The method by Bardet and Choucair [Int. J. Numer. Anal. Meth. Geomech. 15(1):1–19, 1991] is generalized for an arbitrary number of controlling and controlled variables. The approach is illustrated for time‐dependent thermo–hydro–mechanical constitutive models. The purpose of the method is to define an endless variety of different laboratory tests declaratively by means of two matrices E and S , which define the constraints on the controlled quantities such as stress, strains, suction, or temperature. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

4.
A constitutive model for simulation of the behavior of unsaturated interfaces is presented here. The model is an extension of an existing critical state compatible interface model for dry and saturated interfaces that was already proposed by one of the authors [Lashkari, A. 2013. Int. J. Numer. Anal. Meth. Geomech. 37 (8): 904–931]. For a proper simulation of the behavior of partially saturated interfaces, the extended model is formulated in terms of two pairs of work conjugate stress–strain‐like variables. The modified model simulations are compared with the existing data of dry, unsaturated, and saturated interfaces. For each interface type, it is shown that the proposed model can capture the essential elements of the behavior using a unique set of parameters. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

5.
The governing equations for an elasto‐plastic constitutive model for frictional materials such as soil, rock, and concrete are presented, and the incremental form is indicated in preparation for implementation of the model in a user‐defined module for finite element calculations. This isotropic, work‐hardening and ‐softening model employs a single yield surface, it incorporates non‐associated plastic flow, and its capability of capturing the behaviour of different types of frictional materials under various three‐dimensional conditions has been demonstrated by comparison with measured behaviour, as presented in the literature. The incrementalization procedure is indicated and the resulting equations for the single hardening model are presented together with parameters for a dense sand. Following the implementation of the model, these parameters are used for evaluation of different integration schemes as presented in a companion paper by Jakobsen and Lade (Int. J. Numer. Anal. Meth. Geomech. 2002; 26 :661). Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

6.
  相似文献   

7.
The objective of this work was to evaluate the theoretical capabilities of the non-destructive impact–response method in estimating the length and cross-sectional area of intact piles. Three-dimensional (3-D) axisymmetric finite element models were developed to simulate the testing. The results obtained were compared to one-dimensional solutions to evaluate the importance of 3-D effects. Extensive parametric studies were then performed on piles without defects. In each parametric study, the results from the direct use of time histories of displacements or velocities, the mobility function and the Fourier transform of the recorded displacements (impact-echo method) were compared in order to assess their relative advantages and disadvantages. The effects of the relative stiffness of the surrounding soil to that of the pile and of the embedment depth were also investigated for all three methods. In a companion paper the use of these procedures to detect defects such as bulbs (increases in the cross-sectional area of the pile) or necks (decreases in area) is studied. © 1997 by John Wiley & Sons, Ltd. Int. J. Numer. Anal. Meth. Geomech., vol. 21, 255–275 (1997)  相似文献   

8.
In this paper, a series of multimaterial benchmark problems in saturated and partially saturated two‐phase and three‐phase deforming porous media are addressed. To solve the process of fluid flow in partially saturated porous media, a fully coupled three‐phase formulation is developed on the basis of available experimental relations for updating saturation and permeabilities during the analysis. The well‐known element free Galerkin mesh‐free method is adopted. The partition of unity property of MLS shape functions allows for the field variables to be extrinsically enriched by appropriate functions that introduce existing discontinuities in the solution field. Enrichment of the main unknowns including solid displacement, water phase pressure, and gas phase pressure are accounted for, and a suitable enrichment strategy for different discontinuity types are discussed. In the case of weak discontinuity, the enrichment technique previously used by Krongauz and Belytschko [Int. J. Numer. Meth. Engng., 1998; 41:1215–1233] is selected. As these functions possess discontinuity in their first derivatives, they can be used for modeling material interfaces, generating only minor oscillations in derivative fields (strain and pressure gradients for multiphase porous media), as opposed to unenriched and constrained mesh‐free methods. Different problems of multimaterial poro‐elasticity including fully saturated, partially saturated one, and two‐phase flows under the assumption of fully coupled extended formulation of Biot are examined. As a further development, problems involved with both material interface and impermeable discontinuities, where no fluid exchange is permitted across the discontinuity, are considered and numerically discussed. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

9.

A series of nonlinear dynamic FE simulations have been performed to investigate the seismic performance of a flexible propped retaining wall in a saturated clay. The simulations have been carried out considering different acceleration time histories at the bedrock and two different inelastic soil models: the classical elastoplastic Modified Cam Clay model and the advanced hypoplastic model for clays proposed by Ma?ín (Int J Numer Anal Methods Geomech 29:311–336, 2005) equipped with the intergranular strains extension. The results of the simulations highlight the major role played by the choice of the constitutive model for the soil on the predicted seismic response, in terms of predicted wall displacements and structural loads. In particular, the results show that a key role is played by the model ability to correctly reproduce soil dilatancy as a function of the current stress state and loading history. This has a major impact on the inelastic volumetric deformations accumulated during the undrained seismic shearing and on the development of excess pore water pressures around the excavations.

  相似文献   

10.
Earlier solutions (Bull. Seismol. Soc. Amer. 1985; 75 :1135–1154; Bull. Seismol. Soc. Amer. 1992; 82 :1018–1040) of deformations caused by the movements of a shear or tensile fault in an isotropic half-space for finite rectangular sources of strain nucleus have been extended for a transversely isotropic half-space. Results of integrating previous solutions (Int. J. Numer. Anal. Meth. Geomech. 2001; 25 (10): 1175–1193) of deformations due to a shear or tensile fault in a transversely isotropic half-space for point sources of strain nucleus over the fault plane are presented. In addition, a boundary element (BEM) model (POLY3D:A three-dimensional, polygonal element, displacement discontinuity boundary element computer program with applications to fractures, faults, and cavities in the Earth's crust. M.S. Thesis, Stanford University, Department of Geology, 1993; 62) is given. Different from similar researches (e.g. Thomas), the Akaike's view on Bayesian statistics (Akaike Information Criterion Statistics. D. Reidel Publication: Dordrecht, 1986) is applied for inverting deformations due to a fault to obtain displacement discontinuities on the fault plane. An example is given for checking displacements predicted by proposed analytical expressions. Another example is generated for the use of proposed BEM model. It demonstrates the effectiveness of this model in exploring displacement behaviours of a fault. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

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

12.
The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

13.
This paper generalizes the finite strain Coulomb solution of Vrakas and Anagnostou (Int J Numer Anal Meth Geomech 2014; 38(11): 1131–1148) for the classic tunnel mechanics problem of the ground response curve to elastoplastic grounds satisfying a non‐linear Mohr's failure criterion. A linear (Coulomb‐type) plastic potential function is used, leading to a non‐associated flow law, and edge plastic flow is considered in the plastic zone. The solution for a general non‐linear Mohr's failure criterion is semi‐analytical in that it requires the evaluation of definite integrals. In the special case of the Hoek–Brown criterion, however, these integrals are calculated analytically, resulting in a rigorous closed‐form series solution. The applicability of the derived solution is illustrated through the example of the Yacambú‐Quibor tunnel, where very large deformations were observed when crossing of weak graphitic phyllites. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

14.
The residual mechanism of wave‐induced pore water pressure accumulation in marine sediments is re‐examined. An analytical approximation is derived using a linear relation for pore pressure generation in cyclic loading, and mistakes in previous solutions (Int. J. Numer. Anal. Methods Geomech. 2001; 25 :885–907; J. Offshore Mech. Arctic Eng. (ASME) 1989; 111 (1):1–11) are corrected. A numerical scheme is then employed to solve the case with a non‐linear relation for pore pressure generation. Both analytical and numerical solutions are verified with experimental data (Laboratory and field investigation of wave– sediment interaction. Joseph H. Defrees Hydraulics Laboratory, School of Civil and Environmental Engineering, Cornell University, Ithaca, NY, 1983), and provide a better prediction of pore pressure accumulation than the previous solution (J. Offshore Mech. Arctic Eng. (ASME) 1989; 111 (1):1–11). The parametric study concludes that the pore pressure accumulation and use of full non‐linear relation of pore pressure become more important under the following conditions: (1) large wave amplitude, (2) longer wave period, (3) shallow water, (4) shallow soil and (5) softer soils with a low consolidation coefficient. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

15.
We reply to Prof. Linkov comments on our article entitled “Explicit versus implicit front advancing schemes for the simulation of hydraulic fracture growth” (Int. J. Numer. Anal. Methods Geomech., 2019, 43 (6), 1300–1315). We present additional results indicating that both the implicit and explicit front advancement schemes are robust even in the case of a large stress contrast.  相似文献   

16.
This paper analyses the problem of a fluid‐driven fracture propagating in an impermeable, linear elastic rock with finite toughness. The fracture is driven by injection of an incompressible viscous fluid with power‐law rheology. The relation between the fracture opening and the internal fluid pressure and the fracture propagation in mobile equilibrium are described by equations of linear elastic fracture mechanics (LEFM), and the flow of fluid inside the fracture is governed by the lubrication theory. It is shown that for shear‐thinning fracturing fluids, the fracture propagation regime evolves in time from the toughness‐ to the viscosity‐dominated regime. In the former, dissipation in the viscous fluid flow is negligible compared to the dissipation in extending the fracture in the rock, and in the later, the opposite holds. Corresponding self‐similar asymptotic solutions are given by the zero‐viscosity and zero‐toughness (J. Numer. Anal. Meth. Geomech. 2002; 26 :579–604) solutions, respectively. A transient solution in terms of the crack length, the fracture opening, and the net fluid pressure, which describes the fracture evolution from the early‐time (toughness‐dominated) to the large‐time (viscosity‐dominated) asymptote is presented and some of the implications for the practical range of parameters are discussed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

17.
A simple analytical procedure that applies classical beam-column theory for evaluating passive rockbolt roof reinforcement is presented in this paper. The analytical model is derived from first principles and is capable of modelling any number of reinforcing bolts. Each rockbolt is modelled as a linear spring and the model allows for non-uniform bolt spacing. In this study the rock beam is assumed to be isotropic and linearly elastic for the sake of simplicity. However, the analytical model can be extended to include anisotropic rockmass as well as inelastic material behaviour. The solution to the coupled set of governing equations is obtained by using a simple numerical solution procedure. The results from the analytical model indicate that the critical buckling load of a rock beam is strongly influenced by the ambient rock modulus. For salt-rock excavations the rock modulus typically declines with time due to various phenomena, and a diminished modulus could seriously compromise roof stability. The other main conclusion of this study is that rockbolts loose their effectiveness in restraining a roof beam once its critical buckling load is approached. In such a situation, increasing bolt stiffness does not improve its reinforcing action on a roof beam but it enhances the possibility of bolt failure due to anchor pull-out. © 1997 by John Wiley & Sons, Ltd. Int. J. Numer. Anal. Meth. Geomech., vol. 21, 241–253 (1997)  相似文献   

18.
The application of Pastor–Zienkiewicz constitutive model for sands to dynamic consolidation problems is presented in this paper. This model is implemented in a coupled code formulated in terms of displacements for both solid and fluid phases (u?w formulation), which is firstly compared with u?pw formulation for some simple examples. Its range of validity, previously established for elastic problems and harmonic loading, is explored. Once the suitability of the u?w formulation has been ascertained for this kind of dynamic problems in soils, one‐ and two‐dimensional (plane strain) dynamic consolidation numerical examples are provided, aiming to give some light into the physics of this ground improvement technique. A ‘wave of dryness’, observed at the soil surface during the impact in field cases, is numerically reproduced and justified. Some hints on the influence of the loading zone size are also given. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

19.
Research investigations on three-dimensional (3-D) rectangular hydraulic fracture configurations with varying degrees of fluid lag are reported. This paper demonstrates that a 3-D fracture model coupled with fluid lag (a small region of reduced pressure) at the fracture tip can predict very large excess pressure measurements for hydraulic fracture processes. Predictions of fracture propagation based on critical stress intensity factors are extremely sensitive to the pressure profile at the tip of a propagating fracture. This strong sensitivity to the pressure profile at the tip of a hydraulic fracture is more strongly pronounced in 3-D models versus 2-D models because 3-D fractures are clamped at the top and bottom, and pressures in the 3-D fractures that are far removed from the fracture tip have little effect on the stress intensity factor at the fracture tip. This rationale for the excess pressure mechanism is in marked contrast to the crack tip process damage zone assumptions and attendant high rock fracture toughness value hypotheses advanced in the literature. A comparison with field data is presented to illustrate the proposed fracture fluid pressure sensitivity phenomenon. This paper does not attempt to calculate the length of the fluid lag region in a propagating fracture but instead attempts to show that the pressure profile at the tip of the propagating fracture plays a major role in fracture propagation, and this role is magnified in 3-D models. © 1997 by John Wiley & Sons, Ltd. Int. J. Numer. Anal. Meth. Geomech., vol. 21, 229–240 (1997).  相似文献   

20.
This work presents analytical solutions for displacements caused by three‐dimensional point loads in a transversely isotropic full space, in which transversely isotropic planes are inclined with respect to the horizontal loading surface. In the derivation, the triple Fourier transforms are employed toyield integral expressions of Green's displacement; then, the triple inverse Fourier transforms and residue calculus are performed to integrate the contours. The solutions herein indicate that the displacements are governed by (1) the rotation of the transversely isotropic planes (?), (2) the type and degree of material anisotropy (E/E′, ν/ν′, G/G′), (3) the geometric position (r, φ, ξ) and (4) the types of loading (Px, Py, Pz). The solutions are identical to those of Liao and Wang (Int. J. Numer. Anal. Methods Geomechanics 1998; 22 (6):425–447) if the full space is homogeneous and linearly elastic and the transversely isotropic planes are parallel to the horizontal surface. Additionally, a series of parametric study is conducted to demonstrate the presented solutions, and to elucidate the effect of the aforementioned factors on the displacements. The results demonstrate that the displacements in the infinite isotropic/transversely isotropic rocks, subjected to three‐dimensional point loads could be easily determined using the proposed solutions. Also, these solutions could realistically imitate the actual stratum of loading situations in numerous areas of engineering. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

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

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