共查询到20条相似文献，搜索用时 62 毫秒
1.
Mahdi Taiebat Boris Jeremić Yannis F. Dafalias Amir M. Kaynia Zhao Cheng 《Soil Dynamics and Earthquake Engineering》2010
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. 相似文献
2.
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. 相似文献
3.
4.
两相饱和多孔介质的动力响应问题在地震工程领域具有重要的研究意义,由于涉及到固相和液相的动力耦合,使得该问题的求解尤为复杂。本文利用Comsol在求解多场耦合问题上的优点,针对Biot饱和多孔介质u-U耦合形式下的波动方程特征,经过一系列微分算子运算和矩阵变换得到导数形式下的波动方程,基于Comsol Multiphysics提供的广义偏微分方程模式对变形后的波动方程进行求解,并把改进后的无限元边界应用到无限域动力问题的模拟中。通过与饱和多孔介质动力响应的解析解进行对比,验证模型求解技术的可行性和正确性,并在此基础上讨论饱和土地基中空沟隔振效果与饱和土体参数孔隙率、泊松比的关系。通过研究分析,可以为饱和土地基中空沟隔振设计提供一些有价值的参考。 相似文献
5.
Uplift mechanism for a shallow-buried structure in liquefiable sand subjected to seismic load: centrifuge model test and DEM modeling 总被引：2，自引：1，他引：1
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. 相似文献
6.
Amongst different existing techniques developed for absorbing boundaries in transient dynamic analysis, the paraxial approximation represents an elegant framework. This approach is extended to the case of saturated porous media, which is of vital interest in many earthquake engineering problems. With respect to the frequency content of the dynamic process, several formulations exist for modelling the two-phase behaviou of soil or rock materials. The existence of a shear wave and two dilatational waves, which are standard features of porous media, is verified for the adopted formulation. The paraxial approximation is presented in the case of the u – p Biot two-phase dynamic formulation, suitable for seismic analyses and implemented by using a finite element approach. The efficiency of these boundaries is further verified through same numerical examples. 相似文献
7.
8.
9.
M. G. Markov 《Izvestiya Physics of the Solid Earth》2009,45(9):769-776
The problem about the refraction of elastic waves at the interface of two half-spaces (rocks) filled with immiscible fluids
is solved with the use of modified boundary conditions of the dynamics of saturated porous media [Nagy and Nayfeh, 1995].
The solution is obtained within the framework of the Frenkel’-Biot theory with allowance for the surface tension at the interface
of fluids for a porous medium formed by two half-spaces, which differ only in the properties of the fluids filling them. Practically
important cases of the reflection from the fluid-gas (water-air) and fluid-fluid (oil-water) interfaces are considered in
detail. The calculations are performed for both harmonic waves and pulses. The possibility to determine in principle the structural
factor characterizing the pore-space geometry of rocks from measured dynamic parameters of reflected waves is shown. 相似文献
10.
针对饱和多孔介质中热弹性波的传播特性问题,基于多孔介质理论和广义的热弹性模型,研究平面S波在饱和多孔热弹性介质边界上的反射问题。以考虑流-固耦合的饱和多孔介质波动方程和热-弹耦合的广义热弹性基本方程出发,建立饱和多孔介质的热-流-固耦合弹性波动模型。通过引入势函数并考虑自由透水和绝热的边界条件,经过理论推导最终给出在饱和多孔热弹性介质边界上的四种反射波的振幅反射率的理论表达式。在此基础上进行数值计算,分别讨论平面S波的入射频率、入射角和热膨胀系数等参数对四种反射波的振幅反射率的影响情况。结果表明:各反射波的振幅反射率分别随频率和热膨胀系数的增大而增大,同时也受到平面S波入射角变化的影响。该结论对于土动力学的理论研究及其相关的工程勘探具有一定的指导意义。 相似文献
11.
饱和砂土液化是由地震引起的一种最常见的工程地质现象,也是造成重大地震灾害的主要原因之一。由于成因的复杂性和所造成灾害的严重性,饱和砂土液化一直是土动力学和岩土地震工程研究领域的重要课题。针对饱和砂土液化问题,基于开源地震工程数值计算平台OpenSees,对材料库中的4种砂土本构模型进行数值计算。采用二维u-p单元模拟土颗粒位移和孔隙水压力,分析和对比4种模型在循环动力荷载作用下的加速度、超孔隙水压力、位移、剪应力-剪应变和平均有效应力路径方面的响应结果。研究结果表明:(1)砂土对输入加速度表现出一定的放大效应,对于不同的模型,该放大效应存在一些差异;(2) Stress Density模型在循环动力荷载作用下易产生永久变形;(3)在循环动力荷载作用下,PDMY模型和CycLiqCPSP模型的强度逐渐降低,直到完全消失;(4) Stress Density模型和Manzari Dafalias模型在循环动力荷载下表现出明显的剪胀效应。研究成果对砂土液化的数值模拟问题具有重要的理论价值,可为饱和砂土的液化模拟和砂土本构模型的选取提供参考。 相似文献
12.
针对含弱渗透性覆盖层的饱和砂土地基进行一组离心机振动台试验,并采用OpenSees对试验模型进行数值模拟。通过模型试验与数值模拟结果对比讨论OpenSees对于饱和砂土地基地震液化模拟的精度;采用水平方向的Arias强度表示传入某一位置的地震动强度,并以液化时水平方向Arias强度作为该土层的抗液化强度;采用OpenSees计算不同地震动输入时饱和砂土的反应,以此检验Arias强度作为抗液化强度的准确性。结果表明,引起饱和砂土液化所需要的地震动强度随深度增加而增加;当传入的地震强度达到砂土发生液化所需要的地震强度时,该层砂土将会发生液化。 相似文献
13.
14.
15.
16.
Study of pore pressure variation during liquefaction using two constitutive models for sand 总被引：1，自引：0，他引：1
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. 相似文献
17.
18.
Computer simulations are used to calculate the elastic properties of model cemented sandstones composed of two or more mineral phases. Two idealized models are considered – a grain‐overlap clay/quartz mix and a pore‐lining clay/quartz mix. Unlike experimental data, the numerical data exhibit little noise yet cover a wide range of quartz/cement ratios and porosities. The results of the computations are in good agreement with experimental data for clay‐bearing consolidated sandstones. The effective modulus of solid mineral mixtures is found to be relatively insensitive to microstructural detail. It is shown that the Hashin–Shtrikman average is a good estimate for the modulus of the solid mineral mixtures. The distribution of the cement phase is found to have little effect on the computed modulus–porosity relationships. Numerical data for dry and saturated states confirm that Gassmann's equations remain valid for porous materials composed of multiple solid constituents. As noted previously, the Krief relationship successfully describes the porosity dependence of the dry shear modulus, and a recent empirical relationship provides a good estimate for the dry‐rock Poisson's ratio. From the numerical computations, a new empirical model, which requires only a knowledge of system mineralogy, is proposed for the modulus–porosity relationship of isotropic dry or fluid‐saturated porous materials composed of multiple solid constituents. Comparisons with experimental data for clean and shaly sandstones and computations for more complex, three‐mineral (quartz/dolomite/clay) systems show good agreement with the proposed model over a very wide range of porosities. 相似文献
19.
M. D. Sharma 《Journal of Seismology》2017,21(2):425-434
Propagation of surface waves is studied at the pervious boundary of a porous solid saturated with a mixture of two immiscible fluids. An approach, based on continuum mixture theory, is used to derive a secular equation for the propagation of harmonic waves at the stress-free plane surface of this non-dissipative medium. Numerical analysis shows that this secular equation may not represent the propagation of true surface wave in the porous aggregate. Then, this equation is solved numerically for the propagation of pseudo Rayleigh wave or the leaky surface waves. To ensure the existence of pseudo Rayleigh wave, capillary effect between two (wetting and non-wetting) pore-fluids is related to the partial saturation. Effects of porosity and partial saturation coupled with capillary effect are observed on the phase velocity of pseudo Rayleigh waves in sandstone saturated with water-CO2 mixture. 相似文献
20.
黄土具有极强的水敏性和动力易损性，黄土地区多次强震都引起过液化、滑坡等地质灾害，造成了严重的人员伤亡和财产损失，因此振动作用下高含水率黄土的液化问题不容忽视。在大量已有研究的基础上，以宁夏党家岔滑坡为例，研究振动作用下高含水率黄土的液化问题。现场调查发现高含水率滑带土并未达到完全饱和状态（饱和度达95%左右），在新鲜的芯样断面发现有明显的"流态化"液化破坏特征。借助室内试验和数值模拟技术，对党家岔滑坡非饱和黄土层的液化性能及液化发生机理进行分析。结果表明：（1）非饱和黄土层液化发生机理可概括为：地震作用下饱和黄土层孔隙水压力激增，高含水率非饱和黄土层孔压增长响应滞后，随着孔隙水压朝上部消散，地下水向上渗流，当平均有效应力接近0时，高含水率非饱和黄土层发生液化；（2）振动过程中不同饱和度黄土孔隙水压力增长响应具有滞后性，借鉴饱和黄土液化时孔压比的判别标准和Seed简化判别法，初步证实党家岔滑坡高含水率非饱和黄土层可发生振动液化，斜坡前缘和中部土体发生液化的初始饱和度范围分别为68.3%~100%和73.8%~100%，斜坡后缘土体不发生液化。 相似文献