砂土液化变形的有限元-无网格耦合方法

通过构造Biot固结理论u-p方程的无网格伽辽金-有限元耦合方法,对砂土液化变形问题进行了数值模拟。对于饱和砂土,采用Oka等提出的弹塑性本构模型,同时采用更新的Lagrange计算格式推导了控制方程。耦合方法能够发挥有限元和无网格各自的优点,既避免了由于单元变形扭曲而引起的计算中断,也可节约计算时间,算例验证了该方法在地震液化问题中的有效性。     

基于数字图像的砂土压缩变形模式的试验研究

基于数字图像处理技术,利用图像分析程序包PhotoInfor分析了侧限压缩条件下砂土各加荷阶段土体位移的产生、发展与演化全过程,发现砂土的位移场类似于抛物线形状,位移等值线呈轴对称分布,位移量最大的区域总是发生在上透水石的下底面中心处。单调加荷条件下,位移增量主要集中在局部变形带附近,其范围有继续延展和贯通的趋势,顶部、中部测点的土体位移与整体压缩量基本呈同步增长;顶部的土体位移量始终大于中部,加荷的中后期差别更加明显。与粉质黏土位移场进行对比试验分析,砂土内部的应力传递、变形模式以及荷载的影响范围均表现出明显的差异。上述研究为深入探讨基底压力的分布规律,利用宏观试验现象研究砂土细观结构变化的定量模型提供了新的思路。     

土体应变局部化现象的理论解析

引起土体失稳的应变局部化现象是在特定应力状态下,土体本构产生的分叉特性。基于有限变形理论推导了应变局部化产生的三维解析解。基于应变局部化的理论解析,分析了轴对称和平面应变条件下应变局部化现象在弹塑性硬化阶段的存在性以及剪切带的方向性。 理论分析表明,在轴对称条件下,土体应变局部化产生于土体应力-应变的软化阶段,而平面应变条件下,土体应变局部化一般出现在应力-应变的硬化阶段,其剪切带方向角的理论预测与Arthur等[1]建议值较为一致。     

砂土的渐进破坏及其数值模拟

紧密砂土在排水剪切过程中由于材料的应变软化特性,会出现变形局部化并伴随剪切带的产生。对此类问题进行数值模拟时会遇到很多困难,比如如何模拟紧砂的应变软化特性,如何处理有限元数值分析过程中刚度矩阵的非正定导致计算结果发散问题,以及如何处理应力不连续的问题,拟针对以上问题进行了探讨。     

饱和砂土液化后大变形机制的离散元细观分析

采用颗粒流软件模拟了饱和砂土在不排水条件下的循环剪切试验,研究了不同因素对液化的影响,并进一步分析了饱和砂土液化后宏观变形的基本规律。在此基础上,从孔隙分布角度解释了砂土液化后的大变形的细观物理机制。通过自编程序对颗粒排列和孔隙分布的演化过程进行定量描述,给出孔隙率标准差作为液化后体积收缩势的度量,并研究了孔隙率标准差与液化后大变形的关系。离散元细观数值模拟再现了室内试验中的宏观现象,证实了室内试验中饱和砂土液化后的有限剪切大变形是客观真实的材料响应。土体体积收缩势的累积所导致的孔隙均匀化以及土颗粒间自由空隙增大正是饱和砂土液化后循环剪应变逐渐增大的细观机制。孔隙率标准差作为孔隙均匀化的量化指标,与循环剪应变各周次幅值有良好的相关性。     

饱和砂土动静力特性的网格自适应数值模拟

在有限元数值计算平台上,建立了一套基于后验误差评估和Delaunay三角剖分的网格自适应方案,针对饱和砂土静力受压和地震液化的特性进行模拟。验证了超收敛单元片回归（SPR）误差评估中原用于四边形单元的双线性回归函数在用于三角形单元网格时的适用性和可靠性;在饱和砂土动、静力算例中,网格自适应计算获得的变形、应变、超孔压比等的变化规律与常规有限元结果趋势一致。随着网格的再生成,参考点的位移和全域的平均相对误差逼近精确值。对于初始网格,讨论了合理的自适应程度并应用于地震液化的自适应数值模拟中,也对Delaunay三角剖分实施了一些改进。最终证明该自适应方案在提高计算效率的同时,亦可以保证计算所需的精准度。     

大变形固结理论最终沉降量分析

大变形理论的力学机制较为严格,但由于其固有的复杂性和非线性,目前对大小变形理论计算的沉降量差别仍存在争议。基于连续介质力学基本原理,采用土力学符号约定推导了T.L.描述的大变形固结增量有限元方程组及其总应力分析形式。结合算例,分析了材料线性化和材料几何双线性化简化算法的影响,探讨了较合适的荷载增量步长,讨论了大、小变形法计算的最终沉降量差别。结果表明,荷载增量步长对两种线性化法计算结果较为敏感;大变形法计算的沉降量偏大于小变形结果,这可能是由几何非线性问题的“拉压刚度不同”效应引起;初应力矩阵对最终沉降量有较大影响等。     

减饱和砂土缓倾场地的液化性状分析

减饱和法是一种通过减小砂土地基的饱和度,从而提高地基抗液化强度的新方法。基于减饱和砂土中流体模量同步更新的改进算法对减饱和砂土离心机振动台试验进行了数值分析,并与单一流体模量的简化算法进行了对比分析。结果表明：由于改进算法中考虑了因孔压变化引起的等效流体模量的变化,其计算结果更接近试验结果,而简化算法低估了减饱和砂土的孔压积累。基于改进算法开展了不同饱和度、倾斜角度的缓倾场地上液化变形的数值模拟研究,分析发现超孔隙水压力增长的速度及其峰值随着饱和度的增加而增大,饱和度从100%降低至96.4%,同一深度处的超孔压峰值降低约20%～65%,加速度响应的峰值也有明显的降低;沿地基深度0.75 m到9.00 m,侧向位移减少约20%～50%,表明饱和度的降低对抑制倾斜场地上可液化砂土层的侧向变形有显著效果,随着地基深度的减小,饱和度对于侧向位移的影响越来越明显。     

强间断分析方法在土工结构物渐进破坏过程中的应用

对于超固结黏土和密实砂土等软化材料或非关联塑性材料组成的地基、边坡及挡土墙墙后土体,在其破坏过程中,会产生应变局部化现象,使得控制方程的类型发生改变,从而导致出现数值解不惟一和解的网格相关性等现象。为了克服这些数值困难,基于强间断分析方法,及单元内嵌不连续面的有限元模型,对地基、土坡、墙后土体的渐进破坏过程进行了数值模拟。计算结果表明,单元内嵌不连续面模型可以有效地模拟土工结构失稳破坏过程,并且能够明显地改善采用常规有限元方法所产生的网格尺寸相关性问题。这一方法可作为传统极限平衡法进行稳定分析、承载力分析的有益补充。     

基于自由场大型振动台试验的饱和砂土固-液相变特征研究

基于饱和砂土自由场地大型振动台试验,利用加速度记录的线性方法反演得到模型地基土动剪应力-动剪应变,并根据流体动力学相关理论,引入非牛顿流体中的表观黏度、零剪切黏度等概念,对场地液化后的剪切稀化特征进行探讨,研究了饱和砂土在地震荷载下的固-液相变特征。研究结果表明:大震激励时,饱和砂土层上部土体孔压比达到1,饱和砂土液化,从动剪应力-动剪应变曲线变化中可知液化后的砂土动剪切模量下降明显,这表明土体逐渐软化;从动剪应力、动剪应变推导得出的饱和砂土剪应力-剪应变率变化趋势与非牛顿流体的流变曲线相近,且饱和砂土上部液化土层的表观黏度大幅下降,液化后砂土表现出"剪切稀化"的假塑性流体特性。     

A two-scale constitutive framework for modelling localised deformation in saturated dilative hardening material

Undrained deformation of dilative sand generates negative excess pore pressure. It enhances the strength, which is called dilative hardening. This increased suction is not permanent. The heterogeneity at the grain scale triggers localisations causing local volume changes. The negative hydraulic gradient drives fluid into dilating shear zones. It loosens the soil and diminishes the shear strength. It is essential to understand the mechanism behind this internal drainage and to capture it numerically. The purpose of this paper is to develop a macroscopic constitutive relationship for the undrained deformation of saturated dense sand in the presence of a locally fully or partially drained shear band. Separate constitutive relations are generated for the band and intact material. Both time and scale dependence during pore fluid diffusion in saturated sand are captured, eliminating the mesh dependency for finite element implementations. The model is applied to the Gauss points that satisfy the bifurcation criterion. The proposed method is calibrated to recreate the undrained macroscopic response bestowed by an extra-small mesh. The microscopic behaviours inside and outside shear band predicted by this model are qualitatively in good agreement with individual material point behaviours inside and outside the shear band in the extra-small mesh. Depending on the loading rate and the shear band thickness, the response inside the band can be fully or partially drained, which governs the ultimate global strength. The calibrated model is exploited to simulate an upscaled biaxial compression test with semipermeable boundaries.     

软土地基上路堤填筑的破坏性状分析

在软土地区建设高速公路、铁路等的填筑路堤越来越多,由于软土的不排水抗剪强度较低,软土地区路堤的填筑十分困难,而软土地基的渗透性很低,通常认为路堤的填筑是在不排水条件下进行。经典的设计方法是基于极限平衡法,无法考虑变形等因素,因此,能综合考虑稳定及变形等因素的有限元方法成为分析此类问题的一种有效手段。依据连云港铁路路堤在软土地基上的填筑试验资料,采用大变形有限元方法及Mohr-Coulomb模型,采用现场试验得到的土体强度值,对平面应变条件下路堤填筑的破坏性状及极限填筑高度等进行了分析,并对比了有限元分析结果及现场试验结果,结果表明,采用大变形有限元方法能够较准确的得到填筑路堤的极限高度,并有效地分析填筑路堤的破坏性状。     

饱和悬浮塑料砂流动变形可视环剪试验研究

为了分析剪切条件下零有效应力状态砂土的流动变形规律,对传统环剪仪进行了试样可视化改造,研制了透明环形剪切盒,通过对剪切盒膨胀性能分析及与标准砂的剪切试验对比,验证了环剪装置改造的合理性。通过分析不排水条件下饱和悬浮塑料砂的剪应力-应变曲线,发现其剪切强度具有应变软化的特性。通过分析环剪仪中饱和悬浮塑料砂试样的有效应力,证实了可视环剪试验中的土体基本处于零有效应力状态。开展了饱和悬浮塑料砂的流动变形可视环剪试验研究,结果表明：饱和悬浮塑料砂在固结不排水条件下剪切变形不连续,直接在剪切面发生断裂;在不固结不排水的条件下饱和悬浮塑料砂的剪切变形表现出流动形态,且与剪切速率有关：在低剪切速率下,剪切变形仅在剪切面处形成具有曲线轨迹的流动变形而在其他区域不发生变形;而在高剪切速率下剪切变形为整体的倾斜变形,符合黏性流体的流动变形特征。     

Kinematic modelling of shear band localization using discrete finite elements

Modelling shear band is an important problem in analysing failure of earth structures in soil mechanics. Shear banding is the result of localization of deformation in soil masses. Most finite element schemes are unable to model discrete shear band formation and propagation due to the difficulties in modelling strain and displacement discontinuities. In this paper, a framework to generate shear band elements automatically and continuously is developed. The propagating shear band is modelled using discrete shear band elements by splitting the original finite element mesh. The location or orientation of the shear band is not predetermined in the original finite element mesh. Based on the elasto‐perfect plasticity with an associated flow rule, empirical bifurcation and location criteria are proposed which make band propagation as realistic as possible. Using the Mohr–Coulomb material model, various results from numerical simulations of biaxial tests and passive earth pressure problems have shown that the proposed framework is able to display actual patterns of shear banding in geomaterials. In the numerical examples, the occurrence of multiple shear bands in biaxial test and in the passive earth pressure problem is confirmed by field and laboratory observations. The effects of mesh density and mesh alignment on the shear band patterns and limit loads are also investigated. Copyright © 2003 John Wiley & Sons, Ltd.     

Effects of boundary conditions and partial drainage on cyclic simple shear test results—a numerical study

Owing to imperfect boundary conditions in laboratory soil tests and the possibility of water diffusion inside the soil specimen in undrained tests, the assumption of uniform stress/strain over the sample is not valid. This study presents a qualitative assessment of the effects of non‐uniformities in stresses and strains, as well as effects of water diffusion within the soil sample on the global results of undrained cyclic simple shear tests. The possible implications of those phenomena on the results of liquefaction strength assessment are also discussed. A state‐of‐the‐art finite element code for transient analysis of multi‐phase systems is used to compare results of the so‐called ‘element tests’ (numerical constitutive experiments assuming uniform stress/strain/pore pressure distribution throughout the sample) with results of actual simulations of undrained cyclic simple shear tests using a finite element mesh and realistic boundary conditions. The finite element simulations are performed under various conditions, covering the entire range of practical situations: (1) perfectly drained soil specimen with constant volume, (2) perfectly undrained specimen, and (3) undrained test with possibility of water diffusion within the sample. The results presented here are restricted to strain‐driven tests performed for a loose uniform fine sand with relative density D_r=40%. Effects of system compliance in undrained laboratory simple shear tests are not investigated here. Copyright © 2004 John Wiley & Sons, Ltd.     

Undrained shear band in water saturated granular media: A critical revisiting with numerical examples

By taking into account the continuity of stresses and pore pressure across the boundaries of a shear band, it is theoretically shown that incipient localization may take place in saturated loose contractive sand in the regime of decreasing deviator stress under undrained conditions. The undrained shear band orientation primarily depends on the Poisson's ratio and the dilatancy characteristics of the material. However, under strict isochoric constraint, localization is precluded in dense dilative sand with deviator stress increasing only. If any local volume change takes place with a certain mechanism, undrained shear band would be possible in dense sand, most likely in the regime of decreasing friction angle after the peak friction angle is mobilized. Numerical examples are given to demonstrate the variation of shear band inclinations with void ratios and initial consolidation pressures. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of Lode angle‐dependent failure criteria on shear localization analysis in sand

Geotechnical experiments show that Lode angle‐dependent constitutive formulations are appropriate to describe the failure of geomaterials. In the present study, we have adopted one such class of failure criteria along with a versatile constitutive relationship to theoretically analyze the effects of Lode angle on localized shear deformation or shear band formation in loose sand for both drained and undrained conditions. We determine the variation in the possible stress states for shear localization due to the introduction of Lode angle by considering the localized deformation as a bifurcation problem. Further, similar bifurcation analysis is performed for the stress states along a specific loading path, namely, plane strain compression at the constitutive level. In addition, the plane strain compression tests have been simulated as a boundary value finite element problem to see how Lode angle affects the post‐localization response. Results show that the inclusion of a Lode angle parameter within the failure criterion has considerable effects on the onset, plastic strain, and propagation of shear localization in loose sand specimens. For drained condition, we notice early inception of shear localization and multiple band formation when the Lode angle‐dependent failure criterion is used. Undrained localization characteristics, however, found to be independent of Lode angle consideration.     

常剪应力路径下含气砂土的三轴试验

天然气水合物完全分解时,产生的气体使得能源土孔隙压力急速增加,有效应力减小,进而引起土体液化破坏。此时深海能源土斜坡的应力状态与静力液化失稳过程可简化为含气土在常剪应力排水(或不排水)应力路径下的破坏问题。以此为背景,提出了制备含气砂土试样的改进充气管法,并开展了含气砂土的常剪应力路径三轴试验。22组试验结果表明:同一孔隙比的含气密砂在不同围压与常剪应力下具有相同的失稳线;含气砂土试样失稳时的应力比和体变均随初始相对密实度的增大而增大;含气密砂在常剪应力路径下饱和度对失稳特征影响的规律性在排水与不排水条件下均不明显,但在不排水条件下含气砂土的孔压(或体变)对变形的敏感性降低;含气密砂在常剪应力路径到达失稳点之后,排水条件下是瞬变的液化鼓胀破坏,不排水条件下是渐变的剪切破坏。     

吸力式沉箱基础极限拉拔承载力的数值分析

作为新型的深水海洋基础型式,吸力式沉箱基础被广泛地用于系泊深水海洋设施中,从而承受巨大的倾斜上拔荷载。在上拔荷载水平分量与竖向分量的共同作用下,吸力式沉箱的承载特性及其工作性能评价是海洋工程设计与建设中的关键技术问题之一。然而现有的理论分析与试验研究并不能满足工程实践的需要,因此,对吸力式沉箱基础的极限承载力分析建立了有限元数值计算方法。当沉箱基础在快速拔出过程中,正常固结黏土处于完全不排水状态,沉箱基础发生整体破坏时表现出反向地基承载力失稳模式,此时沉箱基础所发挥的极限承载能力往往最大。为此,在数值计算中直接假定沉箱基础及其周围土体处于完全不排水状态,针对不同的沉箱长径比,分别确定了在竖向上拔荷载和水平拉拔的单独作用下沉箱基础极限承载力。对比发现：竖向上拔极限承载力有限元解能够较好地与理论计算结果相符合,而水平极限承载力解与理论计算结果存在一定的差异。