共查询到20条相似文献，搜索用时 31 毫秒
1.
Lilja R. Bjarnadóttir Denise C. Rüther Monica C. M. Winsborrow Karin Andreassen 《Boreas: An International Journal of Quaternary Research》2013,42(1):84-107
A marine geophysical study reveals a complex deglaciation pattern in the Kveithola trough, W Barents Sea. The data set includes multibeam swath bathymetry and sub‐bottom sediment profiler (chirp) data acquired for the whole extent of a palaeo, marine‐terminating ice stream, along with high‐resolution single‐channel seismic data from chosen profiles. The multibeam data show a geomorphic landform assemblage characteristic of ice streams. The results of a combination of seismic and chirp unit stratigraphy reveal that the seabed geomorphology is governed by a deeper‐lying reflector. The reflector dominates surface expressions of several subglacial and ice‐marginal units, each connected to a separate episode of ice‐margin stillstand/advance. Analysis of the combined data set has resulted in a conceptual model of the ice‐stream retreat. The model depicts complex deglaciation of a small, confined ice‐stream system through episodic retreat. It describes the formation of several generations of grounding‐zone systems, characterized by high meltwater discharges and the deposition of fine‐grained grounding‐line fans. The inferred style of grounding‐zone deposition in Kveithola deviates from that of other accounts, and is suggested to be intermediate in the previously described continuum between morainal banks and grounding‐line wedges. The results of this paper have implications for grounding‐zone theory and should be of interest to modellers of grounding‐line dynamics and ice‐stream retreat. 相似文献
2.
The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter‐layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat‐type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic–perfectly plastic. Condition of slip at the interfaces are determined by a Mohr–Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14 :87–104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub‐vertical slopes. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
3.
In this paper, a numerical model to predict flow‐induced shear failure along pre‐existing fractures is presented. The framework is based on a discrete fracture representation embedded in a continuum describing the damaged matrix. A finite volume method is used to compute both flow and mechanical equilibrium, whereas specifically tailored basis functions are used to account for the physics at discontinuities. The failure criterion is based on a maximum shear strength limit, which changes with varying compressive stress on the fracture manifold. The displacements along fracture manifolds are obtained such that force balance is achieved under conditions, where shear stress of the failing fracture segment is constrained to the maximum shear strength at the segment. Simultaneously, the fluid pressure is computed independently of the shear slip. A relaxation model approach is used to obtain the maximum shear limit on the fracture manifold, which leads to grid convergence. 相似文献
4.
A three‐dimensional phenomenological model is developed to describe the long‐term creep of gypsum rock materials. The approach is based on the framework of continuum damage mechanics where coupling with viscoelasticity is adopted. Specifically, a local damage model based on the concept of yield surface is proposed and deeply investigated. Among the many possibilities, we choose in this work its coupling with a generalized Kelvin–Voigt rheological model to formulate the whole behavior. Long‐term as well as short‐term relaxation processes can be integrated in the model by means of as many as necessary viscoelastic processes. The numerical discretization is described for an easy integration within a finite element procedure. Finally, a set of numerical simulations is given to show the possibilities of the presented model. It shows good agreement with some experimental results found in the literature. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
5.
The purpose of this paper is to present a physically based plasticity model for non‐coaxial granular materials. The model, which we shall call the double slip and rotation rate model (DSR2 model), is a pair of kinematic equations governing the velocity field. The model is based on a discrete micro‐analysis of the kinematics of particles in contact, and is formulated by introducing a quantity called the averaged micro‐pure rotation rate (APR) into the unified plasticity model which was proposed by one of the authors. Our macro–micro mechanical analysis shows that the APR is a non‐linear function of, among other quantities, the macro‐rotation rate of the major principal axis of stress taken in the opposite sense. The requirement of energy dissipation used in the double‐sliding free‐rotating model appears to be unduly restrictive as a constitutive assumption in continuum models. In the DSR2 model the APR tensor and the spin tensor are directly linked with non‐coaxiality of the stress and deformation rate tensors. We also propose a simplified plasticity model based on the DSR2 model for a class of dilatant materials, and analyse its material stability. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
6.
A numerical investigation was made of the relationships between fracture initiation, growth, stress field and boundary conditions. Two-dimensional plane strain continuum models were used in which fractures appeared as zones of strain localization developed through application of a strain softening Mohr–Coulomb constitutive model. R and R′ fractures developed first, followed by Y fractures at larger strains. The models showed that equal development of conjugate R and R′ fractures is easily changed to favor one or the other set by minor variations in model initial conditions. Strength loss in fractures caused stress field rotations in regions bounded by fractures, altering the orientation of subsequent fractures. The amount and sense of stress field rotation is dependent on the strength loss during displacement on the fractures, the orientation of fractures, and on the boundary conditions. Y oriented fractures could be explained on the basis of a Mohr–Coulomb failure criterion provided that stress field rotation is accounted for. Monitoring of fracture slip activity showed that, under conditions of constant boundary velocity, slip was discontinuous in time, alternating on fractures throughout the model. 相似文献
7.
The objective of this paper is to formulate and validate an accurate MPM approach for the numerical simulation of the large displacement of membranes containing soil. In the proposed approach, the membrane is discretised by a surface mesh that allows accurate simulation of membrane stresses. The membrane is free to move through a three‐dimensional grid for a continuum consisting of tetrahedral elements. The approach is applied to model a geocontainer being released from a split barge, taking into account the frictional contact between the geotextile and the barge. No‐slip contact is assumed between the geotextile and the soil inside. The effect of geocontainer interaction is investigated by dropping a second container. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
8.
A numerical model to predict landslide movements along pre‐existing slip surfaces from rainfall data is presented. The model comprises: a transient seepage finite‐element analysis to compute the variations of pore water pressures due to rainfall; a limit equilibrium stability analysis to compute the factors of safety along the slip surface associated with transient pore pressure conditions; an empirical relationship between the factor of safety and the rate of displacement of the slide along the slip surface; an optimization algorithm for the calibration of analyses and relationships based on available monitoring data. The model is validated with reference to a well‐monitored active slide in central Italy, characterized by very slow movements occurring within a narrow band of weathered bedrock overlaid by a clayey silt colluvial cover. The model is conveniently divided and presented in two parts: a groundwater model and a kinematic model. In the first part, monthly recorded rainfall data are used as time‐dependent flow boundary conditions of the transient seepage analysis, while piezometric levels are used to calibrate the analysis by minimizing the errors between monitoring data and computed pore pressures. In the second part, measured inclinometric movements are used to calibrate the empirical relationship between the rate of displacement along the slip surface and the factor of safety, whose variation with time is computed by a time‐dependent stability analysis. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
9.
A Lagrangian numerical approach for the simulation of rapid landslide runouts is presented and discussed. The simulation approach is based on the so‐called Particle Finite Element Method. The moving soil mass is assumed to obey a rigid‐viscoplastic, non‐dilatant Drucker–Prager constitutive law, which is cast in the form of a regularized, pressure‐sensitive Bingham model. Unlike in classical formulations of computational fluid mechanics, where no‐slip boundary conditions are assumed, basal slip boundary conditions are introduced to account for the specific nature of the landslide‐basal surface interface. The basal slip conditions are formulated in the form of modified Navier boundary conditions, with a pressure‐sensitive threshold. A special mixed Eulerian–Lagrangian formulation is used for the elements on the basal interface to accommodate the new slip conditions into the Particle Finite Element Method framework. To avoid inconsistencies in the presence of complex shapes of the basal surface, the no‐flux condition through the basal surface is relaxed using a penalty approach. The proposed model is validated by simulating both laboratory tests and a real large‐scale problem, and the critical role of the basal slip is elucidated. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
10.
11.
Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution 
下载免费PDF全文
下载免费PDF全文 We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
12.
Masoud K. Darabi Rashid K. Abu Al‐Rub Eyad A. Masad Dallas N. Little 《国际地质力学数值与分析法杂志》2012,36(7):817-854
Based on the continuum damage mechanics, a general and comprehensive thermodynamic‐based framework for coupling the temperature‐dependent viscoelastic, viscoplastic, and viscodamage behaviors of bituminous materials is presented. This general framework derives systematically Schapery‐type nonlinear viscoelasticity, Perzyna‐type viscoplasticity, and a viscodamage model analogous to the Perzyna‐type viscoplasticity. The resulting constitutive equations are implemented in the well‐known finite element code Abaqus via the user material subroutine UMAT. A systematic procedure for identifying the model parameters is discussed. Finally, the model is validated by comparing the model predictions with a comprehensive set of experimental data on hot mix asphalt that include creep‐recovery, creep, uniaxial constant strain rate, and repeated creep‐recovery tests in both tension and compression over a range of temperatures, stress levels, and strain rates. Comparisons between model predictions and experimental measurements show that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
13.
Earthquakes are recognized as resulting from a stick–slip frictional instability along faults. Based on the node‐to‐point contact element strategy (an arbitrarily shaped contact element strategy applied with the static‐explicit algorithm for modelling non‐linear frictional contact problems proposed by authors), a finite element code for modelling the 3‐D non‐linear friction contact between deformable bodies has been developed and extended here to analyse the non‐linear stick–slip frictional instability between deformable rocks with a rate‐ and state‐dependent friction law. A typical fault bend model is taken as an application example to be analysed here. The variations of the normal contact force, the frictional force, the transition of stick–slip instable state and the related relative slip velocity along the fault between the deformable rocks and the stress evolution in the total bodies during the different stages are investigated, respectively. The calculated results demonstrate the usefulness of this code for simulating the non‐linear frictional instability between deformable rocks. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
14.
本文将岩石视为颗粒体材料,采用两种模型对巷道围岩的应力、应变及破坏区的分布规律进行了数值模拟。第一种模型是连续介质模型,其中考虑了颗粒、界面及基体。第二种模型是将第一种模型中的基体去掉。研究结果表明:当基体强度参数降低较少时,巷道围岩中的环向和径向应力在传统结果附近波动;当基体强度参数降低较多时,两种应力的波动幅度提高,而且,基体位置的应力向其周围的颗粒或界面转移。第二种模型结果的波动幅度更大。随着基体强度参数的降低,巷道围岩中的应变集中区向深部转移,形成相互交织的滑移线网,滑移线网的位置主要位于基体和界面中,这与第二种模型的结果有明显的差异(多个环向的应变集中区)。 相似文献
15.
A novel conceptual model of the mechanics of sands is developed within an elastic–plastic framework. Central to this model is the realization that volume changes in anisotropic granular materials occur as a result of two fundamentally different mechanisms. The first is purely kinematic, dilative, and is the result of the changes in anisotropic fabric. There is also a second volume change in granular media that occurs as a direct response to changes in stress as in a standard elastic/plastic continuum. The inclusion of the two sources of volume change results in three important datum states. When subjected to isotropic strains, the resulting stress state in granular materials is not isotropic but lies upon the kinematic normal consolidation line. There exists a state at which the fabric‐induced volumetric strain rate becomes equal to the stress‐induced volumetric strain rate making the total plastic volumetric strain rate equal to zero. Granular response changes from contractive to dilative at this phase transformation line. The third datum state is the one in which the stress‐induced volumetric strain rate is zero. The sand, however, continues to dilate at this state with the difference between stress and dilation ratio a constant as predicted by Taylor's stress–dilatancy rule. These predictions are shown in accordance with experimental data from a series of drained tests and undrained on Ottawa sand. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
16.
17.
In this paper we assess two competing tectonic models for the development of the Isa Superbasin (ca 1725–1590 Ma) in the Western Fold Belt of the Mt Isa terrane. In the ‘episodic rift‐sag’ tectonic model the basin architecture is envisaged as similar to that of a Basin and Range province characterised by widespread half‐graben development. According to this model, the Isa Superbasin evolved during three stages of the Mt Isa Rift Event. Stage I involved intracontinental extension, half‐graben development, the emergence of fault scarps and tilt‐blocks, and bimodal volcanism. Stage II involved episodic rifting and sag during intervening periods of tectonic quiescence. Stage III was dominated by thermal relaxation of the lithosphere with transient episodes of extension. Sedimentation was controlled by the development of arrays of half‐grabens bounded by intrabasinal transverse or transfer faults. The competing ‘strike‐slip’ model was developed for the Gun Supersequence stratigraphic interval of the Isa Superbasin (during stage II and the beginning of stage III). According to this model, sinistral movements along north‐northeast‐orientated strike‐slip faults took place, with oblique movements along northwest‐orientated faults. This resulted in the deposition of southeast‐thickening ramp sequences with local sub‐basin depocentres forming to the west and north of north‐northeast‐ and northwest‐trending faults, respectively. It is proposed that dilation zones focused magmatism (e.g. Sybella Granite) and transfer of strike‐slip movement resulted in transient uplift along the western margin of the Mt Gordon Arch. Our analysis supports the ‘episodic rift‐sag’ model. We find that the inferred architecture for the strike‐slip model correlates poorly with the observed structural elements. Interpretation is made difficult because there has been significant modification and reorientation of fault geometry during the Isan Orogeny and these effects need to be removed before any assertion as to the basin structure is made. Strike‐slip faulting does not explain the regional‐scale pattern of basin subsidence. The ‘episodic rift‐sag’ model explains the macroscopic geometry of the Isa Superbasin and is consistent with the detailed sedimentological analysis of basin facies architecture, and the structural history and geometry. 相似文献
18.
Frédéric Dufour Grégory Legrain Gilles Pijaudier‐Cabot Antonio Huerta 《国际地质力学数值与分析法杂志》2012,36(16):1813-1830
Damage models are capable of representing crack initiation and mimicking crack propagation within a continuum framework. Thus, in principle, they do not describe crack openings. In durability analyses of concrete structures however, transfer properties are a key issue controlled by crack propagation and crack opening. We extend here a one‐dimensional approach for estimating a crack opening from a continuum‐based finite element calculation to two‐dimensional cases. The technique operates in the case of mode I cracking described in a continuum setting by a nonlocal isotropic damage model. We used the global tracking method to compute the idealized crack location as a post‐treatment procedure. The original one‐dimensional problem devised in Dufour et al. [4] is recovered as profiles of deformation orthogonal to the idealized crack direction are computed. An estimate of the crack opening and an error indicator are computed by comparing finite element deformation profiles and theoretical profiles corresponding to a displacement discontinuity. Two estimates have been considered: In the strong approach, the maxima of the profiles are assumed to be equal; in the weak approach, the integrals of each profile are set equal. Two‐dimensional numerical calculations show that the weak estimates perform better than do the strong ones. Error indicators, defined as the distance between the numerical and theoretical profiles, are less than a few percentages. In the case of a three‐point bending, test results are in good agreement with experimental data, with an error lower than 10% for widely opened crack (> 40µm). Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
19.
This paper investigates the possibility of interpreting progressive shear failure in hard soils and soft rocks as the result of shear propagation of a pre‐existing natural defect. This is done through the application of the principles of fracture mechanics, a slip‐weakening model (SWM) being used to simulate the non‐linear zone at the tips of the discontinuity. A numerical implementation of the SWM in a computation method based on the boundary element technique of the displacement discontinuity method (DDM) is presented. The crack and the non‐linear zone at the advancing tip are represented through a set of elements, where the displacement discontinuity (DD) in the tangential direction is determined on the basis of a friction law. A residual friction angle is assumed on the crack elements. Shear resistance decreases on elements in the non‐linear zone from a peak value at the tip, which is characteristic of intact material, to the residual value. The simulation of a uniaxial compressive test in plane strain conditions is carried out to exemplify the numerical methodology. The results emphasize the role played by the critical DD on the mechanical behaviour of the specimen. A validation of the model is shown through the back analysis of some experimental observations. The results of this back analysis show that a non‐linear fracture mechanics approach seems very promising to simulate experimental results, in particular with regards to the shear band evolution pattern. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
20.
We propose a double slip non-coaxial plastic model within the framework of a Cosserat continuum theory. In a Cosserat continuum,
a material point possesses the degrees of freedom of an infinitesimal rigid body: two translations and one rotation in 2D.
We formulate the plastic model into viscous-plastic constitutive relationships and illustrate the viscous-plastic behaviour
of the model by means of numerical solution of a simple shear problem. 相似文献