共查询到20条相似文献，搜索用时 31 毫秒
1.
2.
Rohit Tiwari Tanusree Chakraborty Vasant Matsagar 《Geotechnical and Geological Engineering》2017,35(4):1491-1512
The present work deals with the three-dimensional nonlinear finite element (FE) analyses of the tunnel in soil subjected to internal blast loading. The analyses are performed using the coupled Eulerian–Lagrangian analysis tool in FE software Abaqus/Explicit. The soil and reinforced concrete lining are modeled using the Lagrangian elements. The explosive Trinitrotoluene (TNT) is modeled using the Eulerian elements. The stress–strain response of soil, concrete, and reinforcement are simulated using strain rate dependent Drucker–Prager plasticity, concrete damaged plasticity and Johnson–Cook (J–C) plasticity models, respectively. The pressure–volume relationship of the TNT explosive is simulated using the Jones-Wilkins-Lee equation of state. Parametric sensitivity studies have been performed for different (1) tunnel lining thicknesses, (2) explosive charge weights and (3) angles of internal friction of soil. It is observed from the results that blast induced pressure on the tunnel lining increases with the increase in charge weight. Both the lining and the surrounding soil undergo significant deformation. The deformation of the tunnel lining increases with increasing charge weight and decreases with increasing lining thickness and increasing the angle of internal friction of soil. Blast-induced velocity in soil attenuates with increasing distance from the source of the blast. 相似文献
3.
The migration of contaminant through soil is usually modeled using the advection‐dispersion equation and assumes that the porous media is stationary without introducing a constitutive equation to represent soil structure. Consequently, time‐dependent deformation induced by soil consolidation or physical remediation is not considered, despite the need to consider these variables during planning for the remediation of contaminated ground, the prediction of contaminated groundwater movement, and the design of engineered landfills. This study focuses on the numerical modeling of solute transfer during consolidation as a first step to resolve some of these issues. We combine a coupling theory‐based mass conservation law for soil‐fluid‐solute phases with finite element modeling to simulate solute transfer during deformation and groundwater convection. We also assessed the sensitivity of solute transfer to the initial boundary conditions. The modeling shows the migration of solute toward the ground surface as a result of ground settlement and the dissipation of excess pore water pressure. The form of solute transport is dependent on the ground conditions, including factors such as the loading schedule, contamination depth, and water content. The results indicate that an understanding of the interaction between coupling phases is essential in predicting solute transfer in ground deformation and could provide an appropriate approach to ground management for soil remediation. 相似文献
4.
Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
5.
An efficient finite–discrete element method applicable for the analysis of quasi‐static nonlinear soil–structure interaction problems involving large deformations in three‐dimensional space was presented in this paper. The present method differs from previous approaches in that the use of very fine mesh and small time steps was not needed to stabilize the calculation. The domain involving the large displacement was modeled using discrete elements, whereas the rest of the domain was modeled using finite elements. Forces acting on the discrete and finite elements were related by introducing interface elements at the boundary of the two domains. To improve the stability of the developed method, we used explicit time integration with different damping schemes applied to each domain to relax the system and to reach stability condition. With appropriate damping schemes, a relatively coarse finite element mesh can be used, resulting in significant savings in the computation time. The proposed algorithm was validated using three different benchmark problems, and the numerical results were compared with existing analytical and numerical solutions. The algorithm performance in solving practical soil–structure interaction problems was also investigated by simulating a large‐scale soft ground tunneling problem involving soil loss near an existing lining. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
6.
B. A. Cheeseman S. Wolf C. F. Yen R. Skaggs 《Fragblast: International Journal for Blasting and Fragmentation》2006,10(1):1-8
Shallow buried explosives pose a significant threat to lightweight vehicles and their onboard personnel. To date, designers of lightweight vehicles are limited in their knowledge of what occurs during the blast. The high intensity, short term loading imparted by the explosion is enormously complex and can be significantly affected by a number of parameters including the size, shape, type, detonation point and depth of burial (DOB) of the explosive and the type, density and water content of the soil. Recent advancements in numerical simulations have enabled the complex blast event to be accurately modelled by coupling Eulerian and Lagrangian analyses: the former is well suited to modelling the blast and while the latter, the structural response. Further validation of the modelling technique is considered in the current paper, which details simulations performed utilising the coupled Eulerian-Lagrangian analysis to study the blast output of explosives buried in saturated sand. These experiments varied explosive charge size, its depth of burial, the target stand-off (SO) distance and the dimensions of the target plate. The investigation concludes with a discussion of the accuracy of the numerical simulations when compared with the experimental observations. 相似文献
7.
In this paper we present a study of the performance of sparse iterative solvers regarding the resolution of three-dimensional and non-linear problems encountered in soil/structure interaction. It is composed of two parts. In the first one, we present briefly iterative methods and preconditioners used in this study, then we analyse their performance on three soil/structure interaction problems: a shallow foundation under a vertical loading, a single pile subjected to a lateral loading and the construction of a lined tunnel in a soft soil. Tests are performed assuming an elastic–perfectly plastic constitutive law for the soil material with a non-associated Mohr–Coulomb flow rule. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
8.
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. 相似文献
9.
Three porous media flow problems, in which the fluid mechanical interactions are critical, are studied in a mesoscopic–microscopic coupling system. In this system, fluid flow in the pore space is explicitly modeled at mesoscopic level by the lattice Boltzmann method, the geometrical representation and the mechanical behavior of the solid skeleton are modeled at microscopic level by the particulate distinct element method (DEM), and the interfacial interaction between the fluid and the solids is resolved by an immersed boundary scheme. In the first benchmark problem, the well‐known and frequently utilized Ergun equation is validated in periodic particle and periodic pore models. In the second problem, the upward seepage problem is simulated over three stages: The settlement of the column of sphere under gravity loading is measured to illustrate the accuracy of the DEM scheme; the system is solved to hydrostatic state with pore space filled with fluid, showing that the buoyancy effect is captured correctly in the mesoscopic–microscopic coupling system; then, the flow with constant rate is supplied at the bottom of the column; the swelling of the ground surface and pore pressure development from the numerical simulation are compared with the predictions of the macroscopic consolidation theory. In the third problem, the fluid‐flow‐induced collapse of a sand arch inside a perforation cavity is tested to illustrate a more practical application of the developed system. Through comparing simulation results with analytical solutions, empirical law and physical laboratory observations, it is demonstrated that the developed lattice Boltzmann–distinct element coupling system is a powerful fundamental research tool for investigating hydromechanical physics in porous media flow. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
10.
Analysis of large deformation of geomaterials subjected to time‐varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator‐split arbitrary Lagrangian–Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid–fluid coupling and strong material non‐linearity. Each time step of the operator‐split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one‐dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
11.
An objective of this paper is to demonstrate that the small strain model developed by the authors can be incorporated into the conventional kinematic hardening plasticity framework to predict pre‐failure defor mations. The constitutive model described in this paper is constituted by three elliptical yield surfaces in triaxial stress space. Two inner surfaces are rotated ellipses of the same shape, representing the boundaries of the linear elastic and small strain regions, while the third surface is the modified Cam clay large‐scale yield surface. Within the linear elastic region, the soil behaviour is elastic with cross‐coupling between the shear and volumetric stress–strain components. Within the small strain region, the soil behaviour is elasto‐plastic, described by the kinematic hardening rule with an infinite number of loading surfaces defined by the incremental energy criterion. Within the large‐scale yield surface, the soil behaviour is elasto‐plastic, described by kinematic and isotropic hardening of the small strain region boundary. Since the yield surfaces have different shapes, the uniqueness of the plastic loading condition imposes a restriction on the ratio between their semi‐diameters. The model requires 12 parameters, which can be determined from a single consolidated undrained triaxial compression test. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
12.
Philip Cardiff 《国际地质力学数值与分析法杂志》2015,39(13):1410-1430
Accurate prediction of the interactions between the nonlinear soil skeleton and the pore fluid under loading plays a vital role in many geotechnical applications. It is therefore important to develop a numerical method that can effectively capture this nonlinear soil‐pore fluid coupling effect. This paper presents the implementation of a new finite volume method code of poro‐elasto‐plasticity soil model. The model is formulated on the basis of Biot's consolidation theory and combined with a perfect plasticity Mohr‐Coulomb constitutive relation. The governing equation system is discretized in a segregated manner, namely, those conventional linear and uncoupled terms are treated implicitly, while those nonlinear and coupled terms are treated explicitly by using any available values from previous time or iteration step. The implicit–explicit discretization leads to a linearized and decoupled algebraic system, which is solved using the fixed‐point iteration method. Upon the convergence of the iterative method, fully nonlinear coupled solutions are obtained. Also explored in this paper is the special way of treating traction boundary in finite volume method compared with FEM. Finally, three numerical test cases are simulated to verify the implementation procedure. It is shown in the simulation results that the implemented solver is capable of and efficient at predicting reasonable soil responses with pore pressure coupling under different loading situations. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
13.
The paper presents a constitutive model for simulating the high strain‐rate behavior of sands. Based on the concepts of critical‐state soil mechanics, the bounding surface plasticity theory and the overstress theory of viscoplasticity, the constitutive model simulates the high strain‐rate behavior of sands under uniaxial, triaxial and multi‐axial loading conditions. The model parameters are determined for Ottawa and Fontainebleau sands, and the performance of the model under extreme transient loading conditions is demonstrated through simulations of split Hopkinson pressure bar tests up to a strain rate of 2000/s. The constitutive model is implemented in a finite‐element analysis software Abaqus to analyze underground tunnels in sandy soil subjected to internal blast loads. Parametric studies are conducted to examine the effect of relative density and type of sand and of the depth of tunnel on the variation of stresses and deformations in the soil adjacent to the tunnels. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
14.
W. L. Fourney U. Leiste R. Bonenberger D. J. Goodings 《Fragblast: International Journal for Blasting and Fragmentation》2005,9(4):205-217
This paper investigates the various mechanisms and parameters that are responsible for delivering impulse to a vehicle that is unfortunate enough to detonate a buried mine. Small scale tests are used to examine the effects of air blast or ejected sand in imparting impulse to a plate that is located above the surface of the saturated soil that contains the explosive. Parameters such as confinement, stand off distance, depth of burial of the explosive, density of the soil, and saturation level of the soil are also examined. 相似文献
15.
This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
16.
Large deformation soil behavior underpins the operation and performance for a wide range of key geotechnical structures and needs to be properly considered in their modeling, analysis, and design. The material point method (MPM) has gained increasing popularity recently over conventional numerical methods such as finite element method (FEM) in tackling large deformation problems. In this study, we present a novel hierarchical coupling scheme to integrate MPM with discrete element method (DEM) for multiscale modeling of large deformation in geomechanics. The MPM is employed to treat a typical boundary value problem that may experience large deformation, and the DEM is used to derive the nonlinear material response from small strain to finite strain required by MPM for each of its material points. The proposed coupling framework not only inherits the advantages of MPM in tackling large deformation engineering problems over the use of FEM (eg, no need for remeshing to avoid mesh distortion in FEM), but also helps avoid the need for complicated, phenomenological assumptions on constitutive material models for soil exhibiting high nonlinearity at finite strain. The proposed framework lends great convenience for us to relate rich grain-scale information and key micromechanical mechanisms to macroscopic observations of granular soils over all deformation levels, from initial small-strain stage en route to large deformation regime before failure. Several classic geomechanics examples are used to demonstrate the key features the new MPM/DEM framework can offer on large deformation simulations, including biaxial compression test, rigid footing, soil-pipe interaction, and soil column collapse. 相似文献
17.
This paper presents a numerical formulation of a three dimensional embedded beam element for the modeling of piles, which incorporates an explicit interaction surface between soil and pile. The formulation is herein implemented for lateral loading of piles but is able to represent soil–pile interaction phenomena in a general manner for different types of loading conditions or ground movements. The model assumes perfect adherence between beam and soil along the interaction surface. The paper presents a comparison of the results obtained by means of the present formulation and by means of a previously formulated embedded pile element without interaction surface, as well as reference semi‐analytical solutions and a fully 3D finite element (FE) model. It is seen that the proposed embedded element provides a better convergence behavior than a previously formulated embedded element and is able to reproduce key features of a full 3D FE model. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
18.
M. A. Hicks 《国际地质力学数值与分析法杂志》2000,24(5):453-476
An adaptive mesh refinement algorithm has been developed for non‐linear computations in geomechanics, based on a smoothed stress–strain finite element formulation. This uses estimates of error in the incremental shear strain invariant to guide the regeneration of unstructured meshes at regular intervals during loading. Following each mesh‐update, no re‐analysis of previous increments with the new mesh is necessary. Algorithm performance has been investigated by analysing a passive earth pressure problem using a linear elastic‐perfectly plastic Mohr–Coulomb soil model. Perfectly drained behaviour has been considered, as have partially drained situations using hydromechanical coupling, while undrained behaviour has been approximated using time steps close to zero. In all cases, mesh adaptivity has been successful in capturing regions of high strain gradient. The results have been compared with analytical solutions. Accurate computations of limit load and shear band orientation have been obtained for a wide range of material dilation angles. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
19.
A novel three‐dimensional particle‐based technique utilizing the discrete element method is proposed to analyze the seismic response of soil‐foundation‐structure systems. The proposed approach is employed to investigate the response of a single‐degree‐of‐freedom structure on a square spread footing founded on a dry granular deposit. The soil is idealized as a collection of spherical particles using discrete element method. The spread footing is modeled as a rigid block composed of clumped particles, and its motion is described by the resultant forces and moments acting upon it. The structure is modeled as a column made of particles that are either clumped to idealize a rigid structure or bonded to simulate a flexible structure of prescribed stiffness. Analysis is done in a fully coupled scheme in time domain while taking into account the effects of soil nonlinear behavior, the possible separation between foundation base and soil caused by rocking, the possible sliding of the footing, and the dynamic soil‐foundation interaction as well as the dynamic characteristics of the superstructure. High fidelity computational simulations comprising about half a million particles were conducted to examine the ability of the proposed technique to model the response of soil‐foundation‐structure systems. The computational approach is able to capture essential dynamic response patterns. The cyclic moment–rotation relationships at the base center point of the footing showed degradation of rotational stiffness by increasing the level of strain. Permanent deformations under the foundation continued to accumulate with the increase in number of loading cycles. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
20.
Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献