共查询到20条相似文献，搜索用时 31 毫秒
1.
A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
2.
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. 相似文献
3.
A viscoplastic subloading soil model for rate‐dependent cyclic anisotropic structured behaviour 
下载免费PDF全文
下载免费PDF全文 This paper presents a new purely viscoplastic soil model based on the subloading surface concept with a mobile centre of homothety, enabling the occurrence of viscoplastic strains inside the yield surface and avoiding the abrupt change in stiffness of the traditional overstress viscoplastic models. This is required for overconsolidated soils. The model is formulated to reproduce the soil rate‐dependent behaviour under cyclic loading (changes in loading direction) and incorporates both initial and induced anisotropy, as well as destructuring. The model shows good qualitative response to some imposed three‐dimensional stress paths under quasi‐inviscid (elastoplastic) behaviour. Some of the main time‐dependent aspects of soil behaviour that the model is capable of reproducing were also illustrated. The capability of the model to adequately reproduce the results from an undrained triaxial test performed on stiff overconsolidated clays from the Lisbon region (Formação de Benfica), with an unloading–reloading deviatoric stress cycle at constant mean stress, that incorporates a series of staggered fast loading and creep stages, was evaluated. The model was able to reproduce well the main observed aspects of the time‐dependent stress–strain response and pore pressure evolution of a stiff overconsolidated clay under complex loading. The revised and generalised viscoplastic subloading surface concept is viable and can be applied to a consistent extension to viscoplasticity, including in the interior of the yield surface, of existing elastoplastic models formulated for soils and other materials. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
4.
This paper presents an elasto‐plastic model for unsaturated compacted soils and experimental results obtained from a series of suction‐controlled triaxial tests on unsaturated compacted clay with different initial densities. The initial density dependency of the compacted soil behaviour is modelled by establishing experimental relationships between the initial density and the corresponding yield stress and thereby between the initial density and the location and slope of normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure surface and the yield surface in the deviatoric plane are given by the extended SMP criterion. A considerable number of the isotropic compression, triaxial compression and extension tests on unsaturated compacted clay with different initial densities were performed using a suction‐controllable triaxial apparatus, to measure the stress–strain–volume change in different stress paths and wetting paths. The model has well‐predicting capabilities to reproduce the mechanical behaviour of specimens compacted under different conditions not only in isotropic compression but also in triaxial compression and triaxial extension. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
5.
Damage‐viscoplastic consistency model for rock fracture in heterogeneous rocks under dynamic loading
Timo Saksala 《国际地质力学数值与分析法杂志》2010,34(10):1041-1062
This paper presents a damage‐viscoplastic consistency model for numerical simulation of brittle fracture in heterogeneous rocks. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept with separate damage variables in tension and compression. This approach does not suffer from ill‐posedness, caused by strain softening, of the underlying boundary/initial value problem since viscoplasticity provides the regularization by introducing a length scale effect under dynamic loading conditions. The model uses the Mohr–Coulomb yield criterion with the Rankine criterion as a tensile cut‐off. The damage law in compression is calibrated via the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition occurring in confined compression at a certain level of confinement. The heterogeneity of rock is accounted for by the statistical approach based on the Weibull distribution. Numerical simulations of confined compression test in plane strain conditions demonstrate a good agreement with the experiments at both the material point and structural levels as the fracture modes are realistically predicted. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
6.
A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
7.
A multiphase coupled elasto‐viscoplastic finite element analysis formulation, based on the theory of porous media, is used to describe the rainfall infiltration process into a one‐dimensional soil column. Using this framework, we have numerically analyzed the generation of pore water pressure and deformations when rainfall is applied to the soil. A parametric study, including rainfall intensity, soil–water characteristic curves, and permeability, is carried out to observe their influence on the changes in pore water pressure and volumetric strain. From the numerical results, it is shown that the generation of pore water pressure and volumetric strain is mainly controlled by material parameters α and n′ that describe the soil–water characteristic curve. A comparison with the laboratory results shows that the proposed method can describe very well the characteristics observed during the experiments of one‐dimensional water infiltration into a layered unsaturated soil column. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
8.
Rate‐dependent behaviour of chalk and other porous rocks has undergone widespread study in geomechanics due to its implications on the performance of engineering structures. We present a rate‐dependent constitutive model for chalk and other porous rocks with several new features. The model formulation is based on a viscoplastic rate‐lines approach in which the axial strain rate depends on the proximity of the stress point to an elliptical reference surface. A non‐associated viscoplastic potential surface and an axial scaling algorithm are used to determine the viscoplastic strain components. The model predicts that axial yields stress varies as a power function of applied axial strain rate, as shown by published laboratory data. Comparisons with published experimental data indicate that the model is capable of reproducing observed rate‐dependent behaviour of chalk under a variety of loading conditions. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
An already available non‐associated elastic–viscoplastic constitutive model with anisotropic strain hardening is modified in order to describe both the constitutive parameter dependency on relative density and the spatio‐temporal evolution of strain localization. To achieve this latter goal, two distinct but similar approaches are introduced: one inspired by the gradient theory and one by the non‐local theory. A one‐dimensional case concerning a simple shear test for a non‐homogeneous infinitely long dense sand specimen is numerically discussed and a finite difference scheme is employed for this purpose. The results obtained by following the two different approaches are critically analysed and compared. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
10.
Piles may be subjected to lateral soil pressures as a result of lateral soil movements from nearby construction‐related activities such as embankment construction or excavation operations. Three‐dimensional finite element analyses have been carried out to investigate the response of a single pile when subjected to lateral soil movements. The pile and the soil were modelled using 20‐node quadrilateral brick elements with reduced integration. For compatibility between the soil–pile interface elements, 27‐node quadrilateral brick elements with reduced integration were used to model the soil around the pile adjacent to the soil–pile interface. A Mohr–Coulomb elastic–plastic constitutive model with large‐strain mode was assumed for the soil. The analyses indicate that the behaviour of the pile was significantly influenced by the pile flexibility, the magnitude of soil movement, the pile head boundary conditions, the shape of the soil movement profile and the thickness of the moving soil mass. Reasonable agreement is found between some existing published solutions and those developed herein. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
11.
A three‐phase soil model is proposed to simulate stress wave propagation in soil mass to blast loading. The soil is modelled as a three‐phase mass that includes the solid particles, water and air. It is considered as a structure that the solid particles form a skeleton and their voids are filled with water and air. The equation of state (EOS) of the soil is derived. The elastic–plastic theory is adopted to model the constitutive relation of the soil skeleton. The damage of the soil skeleton is also modelled. The Drucker–Prager strength model including the strain rate effect is used to describe the strength of the soil skeleton. The model is implemented into a hydrocode Autodyn. The recorded results obtained by explosion tests in soil are used to validate the proposed model. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
12.
The aim of this paper is to present a methodology for identifying the soil parameters controlling the delayed behaviour from laboratory and in situ pressuremeter tests by using an elasto‐viscoplastic model (EVP‐MCC) based on Perzyna's overstress theory and on the elasto‐plastic Modified Cam Clay model. The influence of both the model parameters and the soil permeability was studied under the loading condition of pressuremeter tests by coupling the proposed model equations with Biot's consolidation theory. On the basis of the parametric study, a methodology for identifying model parameters and soil permeability by inverse analysis from three levels of constant strain rate pressuremeter tests was then proposed and applied on tests performed on natural Saint‐Herblain clay. The methodology was validated by comparing the optimized values of soil parameters and the values of the same parameters obtained from laboratory test results, and also by using the identified parameters to simulate other tests on the same samples. The analysis of the drainage condition and the strain rate effect during a pressuremeter test demonstrated the coupled influence of consolidation and viscous effects on the test results. The numerical results also showed that the inverse analysis procedure could successfully determine the parameters controlling the time‐dependent soil behaviour. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
13.
Timo Saksala 《国际地质力学数值与分析法杂志》2013,37(3):309-324
This paper considers numerical modelling of rock fracture induced by dynamic bit–rock interaction in percussive drilling. The work presented here extends the author's earlier research on the topic from the axisymmetric case to 3D case. The numerical method for modelling rock fracture includes a constitutive model for rock and a contact mechanics‐based technique to simulate the bit–rock interaction. The constitutive model is based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. This model is improved here from its earlier state by calibrating the softening laws using fracture energies GIc and GIIc in tension and compression, respectively. Moreover, the viscosity modulus in tension is calibrated based on the dynamic Brazilian disc test. With these enhancements, the developed method is applied to 3D case of the bit–rock interaction problem assuming one symmetry plane. Single impact with single and multiple‐button bits is simulated. In the latter case, an initial borehole is modelled in order to simulate the usual in‐situ drilling conditions. The different failure types observed in the experiments as well as the interaction between the buttons resulting in chipping are realistically captured in the simulations. Copyright © 2012 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.
Timo Saksala 《国际地质力学数值与分析法杂志》2010,34(13):1362-1386
This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
16.
17.
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. 相似文献
18.
A closed form analytical solution for two‐dimensional plane strain consolidation of unsaturated soil stratum 下载免费PDF全文
下载免费PDF全文 This paper discusses the excess pore‐air and pore‐water pressure dissipations and the average degree of consolidation in the 2D plane strain consolidation of an unsaturated soil stratum using eigenfunction expansion and Laplace transformation techniques. In this study, the application of a constant external loading on a soil surface is assumed to immediately generate uniformly or linearly distributed initial excess pore pressures. The general solutions consisting of eigenfunctions and eigenvalues are first proposed. The Laplace transform is then applied to convert the time variable t in partial differential equations into the Laplace complex argument s. Once the domain is obtained, a simplified set of equations with variable s can be achieved. The final analytical solutions can be computed by taking a Laplace inverse. The proposed equations predict the two‐dimensional consolidation behaviour of an unsaturated soil stratum capturing the uniformly and linearly distributed initial excess pore pressures. This study investigates the effects of isotropic and anisotropic permeability conditions on variations of excess pore pressures and the average degree of consolidation. Additionally, isochrones of excess pore pressures along vertical and horizontal directions are presented. It is found that the initial distribution of pore pressures, varying with depth, results in considerable effects on the pore‐water pressure dissipation rate whilst it has insignificant effects on the excess pore‐air pressure dissipation rate. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
19.
D. W. Smith 《国际地质力学数值与分析法杂志》2000,24(8):693-722
Contaminant migration through soil is usually modelled mathematically using the dispersion–advection equation. This type of model finds application when planning the remediation of contaminated land, predicting the movement of polluted groundwater and designing engineered landfills. Usually the analysis assumes that the porous media through which the contaminant migrates is stationary. However, the construction of landfills on clay soils means that the soil beneath the landfill will undergo time‐dependent deformation as the soil consolidates. To date, there are no published data on the effect a deforming porous media may have on contaminant transport beneath a landfill; indeed, there appears to be no theory of contaminant migration through a deforming soil. In this paper, a one‐dimensional theory of contaminant migration through a saturated deforming porous media is developed based on a small and large strain analysis of a consolidating soil and conservation of contaminant mass. By selection of suitable parameters, the new transport equation reduces to the familiar one‐dimensional dispersion–advection equation for a saturated soil with linear, reversible, equilibrium controlled sorption of the contaminant onto the soil skeleton. Analytic solutions to a quasi‐steady‐state contaminant transport problem for a deforming media are presented, and a preliminary assessment made of the potential importance of soil deformation on the results of a contaminant migration analysis. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
20.
Dynamic two‐phase interaction of soil can be modelled by a displacement‐based, two‐phase formulation. The finite element method together with a semi‐implicit Euler–Cromer time‐stepping scheme renders a discrete equation that can be solved by recursion. By experience, it is found that the CFL stability condition for undrained wave propagation is not sufficient for the considered two‐phase formulation to be numerically stable at low values of permeability. Because the stability analysis of the two‐phase formulation is onerous, an analysis is performed on a simplified two‐phase formulation that is derived by assuming an incompressible pore fluid. The deformation of saturated porous media is now captured in a single, second‐order partial differential equation, where the energy dissipation associated with the flow of the fluid relative to the soil skeleton is represented by a damping term. The paper focuses on the different options to discretize the damping term and its effect on the stability criterion. Based on the eigenvalue analyses of a single element, it is observed that in addition to the CFL stability condition, the influence of the permeability must be included. This paper introduces a permeability‐dependent stability criterion. The findings are illustrated and validated with an example for the dynamic response of a sand deposit. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献