首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The Mechanical Coupling of Fluid-Filled Granular Material Under Shear   总被引：1，自引：0，他引：1  
The coupled mechanics of fluid-filled granular media controls the physics of many Earth systems, for example saturated soils, fault gouge, and landslide shear zones. It is well established that when the pore fluid pressure rises, the shear resistance of fluid-filled granular systems decreases, and, as a result, catastrophic events such as soil liquefaction, earthquakes, and accelerating landslides may be triggered. Alternatively, when the pore pressure drops, the shear resistance of these geosystems increases. Despite the great importance of the coupled mechanics of grain–fluid systems, the basic physics that controls this coupling is far from understood. Fundamental questions that must be addressed include: what are the processes that control pore fluid pressurization and depressurization in response to deformation of the granular skeleton? and how do variations of pore pressure affect the mechanical strength of the grains skeleton? To answer these questions, a formulation for the pore fluid pressure and flow has been developed from mass and momentum conservation, and is coupled with a granular dynamics algorithm that solves the grain dynamics, to form a fully coupled model. The pore fluid formulation reveals that the evolution of pore pressure obeys viscoelastic rheology in response to pore space variations. Under undrained conditions elastic-like behavior dominates and leads to a linear relationship between pore pressure and overall volumetric strain. Viscous-like behavior dominates under well-drained conditions and leads to a linear relationship between pore pressure and volumetric strain rate. Numerical simulations reveal the possibility of liquefaction under drained and initially over-compacted conditions, which were often believed to be resistant to liquefaction. Under such conditions liquefaction occurs during short compactive phases that punctuate the overall dilative trend. In addition, the previously recognized generation of elevated pore pressure under undrained compactive conditions is observed. Simulations also show that during liquefaction events stress chains are detached, the external load becomes completely supported by the pressurized pore fluid, and shear resistance vanishes.  相似文献   

2.
In this paper, a novel coupled pore-scale model of pore-fluid interacting with discrete particles is presented for modeling liquefaction of saturated granular soil. A microscale idealization of the solid phase is achieved using the discrete element method (DEM) while the fluid phase is modeled at a pore-scale using the lattice Boltzmann method (LBM). The fluid forces applied on the particles are calculated based on the momentum exchange between the fluid and particles. The presented model is based on a first principles formulation in which pore-pressure develops due to actual changes in pore space as particles׳ rearrangement occurs during shaking. The proposed approach is used to model the response of a saturated soil deposit subjected to low and large amplitude seismic excitations. Results of conducted simulations show that at low amplitude shaking, the input motion propagates following the theory of wave propagation in elastic solids. The deposit response to the strong input motion indicates that liquefaction took place and it was due to reduction in void space during shaking that led to buildup in pore-fluid pressure. Soil liquefaction was associated with soil stiffness degradation and significant loss of interparticle contacts. Simulation results also indicate that the level of shaking-induced shear strains and associated volumetric strains play a major role in the onset of liquefaction and the rate of pore-pressure buildup.  相似文献   

3.
松散的饱和砂土在液化之前可作为固体看待,但当过量的孔隙水压达到初始侧应力时它成为液体,之后又恢复其强度。这个过程不能各自单独处理,而应视为一个从固态到液态,或从液态到固态的连续变化过程。因此,整个伴随着地基液化———土体流动现象全过程应作为一个在固态和液态之间相变的系列过程。本文导出了一个简单的基本方程式,可以表达松散饱和砂土在液化和地基侧向滑移现象期间的固-液相转换。该基本方程式可用作桩基系统的动力分析,其适用性通过与弹塑性基本方程的比较得到验证。  相似文献   

4.
Numerical analyses of liquefiable sand are presented in this paper. Liquefaction phenomenon is an undrained response of saturated sandy soils when they are subjected to static or dynamic loads. A fully coupled dynamic computer code is developed to predict the liquefaction potential of a saturated sandy layer. Coupled dynamic field equations of extended Biot's theory with u–P formulation are used to determine the responses of pore fluid and soil skeleton. Generalized Newmark method is employed for integration in time. The soil behavior is modelled by two constitutive models; a critical state two-surface plasticity model, and a densification model. A class ‘B’ analysis of a centrifuge experiment is performed to simulate the dynamic response of level ground sites. The results of the numerical analyses demonstrate the capability of the critical sate two-surface plasticity model in producing pore pressures that are consistent with observations of the behavior of liquefiable sand in the centrifuge test.  相似文献   

5.
Accurate prediction of the liquefaction of saturated soils is based on strong coupling between the pore fluid phase and soil skeleton. A practical numerical method for large strain dynamic analysis of saturated soils is presented. The u–p formulation is used for the governing equations that describe the coupled problem in terms of soil skeleton displacement and excess pore pressure. A mixed finite element and finite difference scheme related to large strain analysis of saturated soils based on the updated Lagrangian method is given. The equilibrium equation of fluid-saturated soils is spatially discretized by the finite element method, whereas terms associated with excess pore pressure in the continuity equation are spatially discretized by the finite difference method. An effective cyclic elasto-plastic constitutive model is adopted to simulate the non-linear behavior of saturated soils under dynamic loading. Several numerical examples that include a saturated soil column and caisson-type quay wall are presented to verify the accuracy of the method and its usefulness and applicability to solutions of large strain liquefaction analysis of saturated soils in practical problems.  相似文献   

6.
The 1995 Hyogoken–Nambu earthquake caused severe liquefaction over wide areas of reclaimed land. Furthermore, the liquefaction induced large ground displacement in horizontal directions, which caused serious damage to foundations of structures. However, few analyses of steel pipe piles based on field investigation have so far been conducted to identify the causes and process of such damage. The authors conducted a soil–pile-structure interaction analysis by applying a multi-lumped-mass-spring model to a steel pipe pile foundation structure to evaluate the causes and process of its damage. The damage process analyzed in the time domain corresponded well with the results of detailed field investigation. It was found that a large bending moment beyond the ultimate plastic moment of the pile foundation structure was induced mainly by the large ground displacement caused by liquefaction before lateral spreading of the ground and that the displacement appeared during the accumulating process of the excess pore water pressure.  相似文献   

7.
The paper presents a mathematical model for the deformation of soil under irregular cyclic loading in the simple-shear conditions. The model includes the possible change in the effective pressure in saturated soil due to the cyclic shearing, the reciprocal influence of the effective pressure on the response of the soil to the shear loading, and the pore pressure dissipation due to the seepage of the pore fluid. The hysteresis curves for the strain–stress relationship are constructed in such a way that they produce both the required backbone curve and the required damping ratio as functions of the strain amplitude. At the same time, the approach enables the constitutive functions involved in the model to be specified in various ways depending on the soil under study. The constitutive functions can be calibrated independently of each other from the conventional cyclic shear tests. The constitutive model is incorporated in the boundary value problem for the dynamic site response analysis of level ground. A numerical solution is presented for the dynamic deformation and liquefaction of soil at the Port Island site during the 1995 Hyogoken-Nambu earthquake.  相似文献   

8.
The seismic risk is fairly high in Hong Kong even though it is located in an intreplate area with low to moderate seismicity. This is because of its high seismic vulnerability due to the presence of many steep loose fill slopes with a marginal static factor of safety, and a high consequence ‘value’ as a result of the dense population and intense economic activity in Hong Kong. In order to investigate the seismic stability and potential flow liquefaction of loose fill slopes, dynamic centrifuge tests in uni-axial and bi-axial directions were performed on saturated model embankments made of loose completely decomposed granite (CDG). Three windowed sinusoidal waves with peak shaking amplitudes ranging from 0.08 g to 0.3 g (prototype scale) were adopted. During the strong uni-axial shaking of 0.3 g, the measured maximum excess pore pressure ratios ranged from 0.70 to 0.85 and a relatively small crest settlement of 5.8 mm (0.22 m prototype) was measured. No soil liquefaction or flow slides were observed. Comparing the results between the strong uni-axial and bi-axial shaking, the maximum pore pressure ratios measured from the bi-axial test varied from 0.75 to 0.87, which were marginally larger than those obtained from the uni-axial test. Although the measured crest settlement during the bi-axial shaking was about 27% larger than that of the uni-axial test, soil liquefaction and flow slide did not occur. These test results suggest that loose CDG fill slopes are likely to be stable under the proposed design PGA ranging from 0.08 to 0.11 g in Hong Kong.  相似文献   

9.
To predict the earthquake response of saturated porous media it is essential to correctly simulate the generation, redistribution, and dissipation of excess pore water pressure during and after earthquake shaking. To this end, a reliable numerical tool requires a dynamic, fully coupled formulation for solid–fluid interaction and a versatile constitutive model. Presented in this paper is a 3D finite element framework that has been developed and utilized for this purpose. The framework employs fully coupled dynamic field equations with a u–p–U formulation for simulation of pore fluid and solid skeleton interaction and a SANISAND constitutive model for response of solid skeleton. After a detailed verification and validation of the formulation and implementation of the developed numerical tool, it is employed in the seismic response of saturated porous media. The study includes examination of the mechanism of propagation of the earthquake-induced shear waves and liquefaction phenomenon in uniform and layered profiles of saturated sand deposits.  相似文献   

10.
The interaction between fluid and sediment particles is widely involved in hydraulic engineering problems. In the current study, an explicit incompressible mesh-free method in the framework of the Moving Particle Semi-implicit(MPS) method is proposed to simulate the interaction between the two phases in submerged conditions. The proposed method solves two sets of the continuity and momentum equations, respectively, for the fluid phase and the sediment phase according to the mixture theory. In th...  相似文献   

11.
This study considers the effects of heat transfer and fluid flow on the thernal, hydrologic, and mechanical response of a fault surface during seismic failure. Numerical modeling techniques are used to account for the coupling of the thermal, fluid-pressure, and stress fields. Results indicate that during an earthquake the failure surface is heated to a tempeature required for the thermal expansion of pore fluids to balance the rate of fluid loss due to flow and the fluid-volume changes due to pore dilatation. Once this condition is established, the pore fluids pressurize and the shear strength decreases rapidly to a value sufficient to maintain the thermal pressurization of pore fluids at near-lithostatic values. If the initial fluid pressure is hydrostatic, the final temperature attained on the failure surface will increase with depth, because a greater pressure increase can occur before a near-lithostatic pressure is reached. The rate at which thermal pressurization proceeds depends primarily on the hydraulic characteristics of the surrounding porous medium, the coefficient of friction on the fault surface, and the slip velocity. If either the permeability exceeds 10–15 m2 or the porous medium compressibility exceeds 10–8 Pa–1, then frictional melting may occur on the fault surface before thermal pressurization becomes significant. If the coefficient of friction is less than 10–1 and if the slip velocity is less than 10–2 msec–1, then it is doubtful that either thermal pressurization or frictional melting on the fault surface could cause a reduction in the dynamic shear strength of a fault during an earthquake event.  相似文献   

12.
The destructive 1999 Chi–Chi earthquake (Mw 7.5) was the largest inland earthquake in Taiwan in the 20th century. Several observations witness the non-linear seismic soil response in sediments during the earthquake. In fact, large settlements as well as evidence of liquefaction attested by sand boils and unusual wet ground surface were observed at some sites. In this paper, we present a seismic response simulation performed with CyberQuake software on a site located within the Chang-Hwa Coastal Industrial Park during the 1999 Chi–Chi earthquake in Taiwan. A non-linear multi-kinematic dynamic constitutive model is implemented in the software. Computed NS, EW and UP ground accelerations obtained with this model under undrained and two-phase assumptions, are in good agreement with the corresponding accelerations recorded at seismic station TCU117, either for peak location, amplitudes or frequency content. In these simulations, liquefaction occurs between depths 1.3 and 11.3 m, which correspond to the observed range attested by in place penetration tests and other liquefaction analyses. Moreover, the computed shear wave velocity profile is very close to post-earthquake shear wave velocity profile derived from correlations with CPT and SPT data. Finally, it is shown that in non-linear computations, even though a 1D geometry is considered, it is necessary to take into account the three components of the input motion.  相似文献   

13.
Compaction or densification of loose saturated soils has been the most popular method of reducing earthquake related liquefaction potential. Such compaction of a foundation soil is only economical when limited in extent, leading to a case of an ‘island’ of improved ground (surrounded by unimproved ground). The behavior of the densified sand surrounded by liquefied loose sand during and following earthquakes is of great importance in order to design the compacted area rationally and optimize both safety and economy. This problem is studied herein by means of dynamic centrifuge model tests. The results of three heavily-instrumented dynamic centrifuge tests on saturated models of side-by-side loose and dense sand profiles are discussed. The test results suggest the following concerns as relates to ‘islands’ of densified soil: (1) there is a potential strength degradation in the densified zone as a result of pore pressure increase due to migration of pore fluid into the island from the adjacent loose liquefied ground; (2) there is a potential for lateral deformation (sliding) within the densified island as the surrounding loose soil liquefies.  相似文献   

14.
A dominant mechanism for residual trapping of a nonwetting fluid in porous media during imbibition is snap-off or the disconnection of a continuous stream of the nonwetting fluid when it passes through pore constrictions and when a criterion based on capillary pressure imbalance is met. While quasi-static criteria for Roof snap-off have been defined for pores based on the imbalance between capillary pressure across the front/tail meniscus and local capillary pressure at the pore throat, and expressed in terms of pore body to pore throat ratio for simplification, we extended the previous quasi-static snap-off criterion by considering the local capillary pressure imbalance between the pore body and the pore throat for both circular and noncircular pores when the wetting film exists. We then used the criterion to analyze results from computational fluid dynamics (CFD) simulations of multi-phase flow with supercritical CO2 as the nonwetting fluid and water as the wetting fluid. The extended criterion successfully described most situations we modeled. Furthermore, we compared fluid interface shape for a noncircular 3D pore predicted by the minimum surface energy (MSE) theory against 3D CFD simulations. While the fluid interface shape at the pore throat for 3D simulation was consistent with the shape predicted by MSE theory, the shape could not be successfully predicted by the MSE theory at the upstream and downstream pore body. Moreover, film flow existed for the noncircular pore at the downstream pore body.  相似文献   

15.
In this paper, a two-dimensional integrated numerical model is developed to examine the influences of cross-anisotropic soil behaviour on the wave-induced residual liquefaction in the vicinity of a pipeline buried in a porous seabed. In the wave model, the RANS (Reynolds Averaged Navier–Stokes) equation is used to govern the wave motion. In the seabed model, the residual soil response in the vicinity of an embedded pipeline is considered with the 2-D elasto-plastic solution, where the phase-resolved shear stress is used as a source for the build-up of residual pore pressure. Classical Biot׳s consolidation equation is used for linking the solid-pore fluid interaction. The validation of the proposed integrated numerical model is conducted by the comparisons with the previous experimental data. Numerical examples show that the pore pressures can accumulate to a large value, thus resulting in a larger area of liquefaction potential in the given anisotropic soil compared to that with isotropic solution. The influences of anisotropic parameters on the wave-induced residual soil response in the vicinity of pipeline are significant. A high rate of pore pressure accumulation and dissipation is observed and the liquefaction potential develops faster as the anisotropic parameters increase. Finally, a simplified approximation based on a detailed parametric investigations is proposed for the evaluation of maximum liquefaction depth (zL) in engineering application.  相似文献   

16.
For single-phase flow through a network model of a porous medium, we report (1) solutions of the Navier–Stokes equation for the flow, (2) micro-particle imaging velocimetry (PIV) measurements of local flow velocity vectors in the “pores throats” and “pore bodies,” and (3) comparisons of the computed and measured velocity vectors. A “two-dimensional” network of cylindrical pores and parallelepiped connecting throats was constructed and used for the measurements. All pore bodies had the same dimensions, but three-different (square cross-section) pore-throat sizes were randomly distributed throughout the network. An unstructured computational grid for flow through an identical network was developed and used to compute the local pressure gradients and flow vectors for several different (macroscopic) flow rates. Numerical solution results were compared with the experimental data, and good agreement was found. Cross-over from Darcy flow to inertial flow was observed in the computational results, and the permeability and inertia coefficients of the network were estimated. The development of inertial flow was seen as a “two-step” process: (1) recirculation zones appeared in more and more pore bodies as the flow rate was increased, and (2) the strengths of individual recirculation zones increased with flow rate. Because each pore-throat and pore-body dimension is known, in this approach an experimental (and/or computed) local Reynolds number is known for every location in the porous medium at which the velocity has been measured (and/or computed).  相似文献   

17.
Based on a centrifuge model test and distinct element method(DEM), this study provides new insights into the uplift response of a shallow-buried structure and the liquefaction mechanism for saturated sand around the structure under seismic action. In the centrifuge test, a high-speed microscopic camera was installed in the structure model, by which the movements of particles around the structure were monitored. Then, a two-dimensional digital image processing technology was used to analyze the microstructure of saturated sand during the shaking event. Herein, a numerical simulation of the centrifuge experiment was conducted using a two-phase(solid and fl uid) fully coupled distinct element code. This code incorporates a particle-fl uid coupling model by means of a "fi xed coarse-grid" fl uid scheme in PFC3D(Particle Flow Code in Three Dimensions), with the modeling parameters partially calibrated based on earlier studies. The physical and numerical models both indicate the uplifts of the shallow-buried structure and the sharp rise in excess pore pressure. The corresponding micro-scale responses and explanations are provided. Overall, the uplift response of an underground structure and the occurrence of liquefaction in saturated sand are predicted successfully by DEM modeling. However, the dynamic responses during the shaking cannot be modeled accurately due to the restricted computer power.  相似文献   

18.
Seismic stability, liquefaction, and deformation of earth structures are critical issues in geotechnical earthquake engineering practice. At present, the equivalent linear approach is considered the ‘state of practice’ in common use. More recently, dynamic analyses incorporating nonlinear, effective-stress-based soil models have been used more frequently in engineering applications. This paper describes a bounding surface hypoplasticity model for sand [Wang ZL. Bounding surface hypoplasticity model for granular soils and its applications. PhD Dissertation for the University of California at Davis, U.M.I. Dissertation Information Service, Order No. 9110679; 1990; Wang ZL, Dafalias YF, Shen CK. Bounding surface hypoplasticity model for sand. ASCE, J Eng Mech 1990;116(5):983–1001; Wang ZL, Makdisi FI. Implementing a bounding surface hypoplasticity model for sand into the FLAC program. In: Proceedings of the international symposium on numerical modeling in geomechanics. Minnesota, USA; 1999. p. 483–90] incorporated into a two-dimensional finite difference analysis program [Itasca Consulting Group, Inc. FLAC (Fast Lagrangian Analysis of Continua), Version 4. Minneapolis, MN; 2000] to perform nonlinear, effective-stress analyses of soil structures. The soil properties needed to support such analyses are generally similar to those currently used for equivalent linear and approximate effective-stress analyses. The advantages of using a nonlinear approach are illustrated by comparison with results from the equivalent linear approach for a rockfill dam. The earthquake performance of a waterfront slope and an earth dam were evaluated to demonstrate the model's ability to simulate pore-pressure generation and liquefaction in cohesionless soils.  相似文献   

19.
饱和黄土在不同外荷载作用下其液化机理具有显著差异.为研究饱和黄土动态液化和静态液化机理的差异性,基于室内动三轴试验和静三轴试验,研究岷县永光饱和黄土动态液化后的动应力与轴向动应变关系、动孔隙水压力比与轴向动应变关系,分析其静态液化后的偏应力与轴向应变关系、孔隙水压力比与轴向应变关系,并结合液化前、后的SEM试验结果,研...  相似文献   

20.
The model presented in the companion paper is validated in both the linear and nonlinear cases under steady-state single frequency harmonic and transient ground motions. The crest acceleration responses of the Santa Felicia earth dam subjected to the 1971 San Fernando earthquake and of the Long Valley earth dam subjected to the strongest of the 1980 Mammoth Lake earthquakes are computed and compared with the motions recorded at the site. Acceleration time histories for the solid and fluid phases in both horizontal and vertical directions, as well as stress-strain and pore water pressure-strain time histories for points along the height of the dam are presented. The ability of the model to simulate the occurrence of liquefaction in a dam is also demonstrated.  相似文献   

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

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