软黏土层一维有限应变固结的超静孔压消散研究

根据土力学固结理论计算分析软黏土层固结过程的超静孔隙水压力值,确定软黏土体固结过程的强度增长,对排水固结法处理软土地基至关重要。软黏土层固结过程中土体变形较大时,有限应变固结理论和小应变固结理论计算分析软黏土固结所得结果差异较大。利用非线性有限元法及程序,通过对软黏土层固结工程算例的计算结果分析,研究了有限应变固结理论和小应变固结理论计算分析软黏土层一维固结超静孔压值消散的差异;探讨了软黏土体一维固结过程中,几何非线性、土体渗透性变化和压缩性变化对超静孔隙水压力消散的影响。研究结果表明,当土体的变形较大时,有限应变固结理论计算出的超静孔压要比小应变固结理论得到的值消散的更快。考虑土体固结过程中渗透性的变化时,超静孔压消散变慢;可用软黏土渗透性变化指数ck 反映渗透性变化对超静孔压消散的影响,渗透性变化指数ck值越小、超静孔压消散越慢。固结过程中软黏土压缩性的大小及变化也影响超静孔压的消散,可用软黏土的压缩指数cc反映固结过程中压缩性的大小及变化对超静孔压消散的影响,软黏土的压缩指数cc越小,固结过程软黏土层中的超静孔压消散越快。     

考虑温度影响的改进西原流变模型及一维固结解

随着我国热能源开发,热力管道工程大规模应用,由此带来的岩土体应力、变形、渗流和温度耦合作用问题已经成为岩土工程领域的研究热点。在长期荷载作用下,岩土体的变形和强度随时间而变化,尤其是软土,具有明显的流变特征。传统的各流变模型并未考虑升温方式和温度对其各参数的影响。为了能更好地反映流变固结特性,考虑软黏土的黏弹塑性,在西原模型基础上,引入温度膨胀系数,建立了热力耦合作用下的改进西原模型,给出了其应力-应变关系,并利用Laplace变换和逆变换求解了一维热固结方程,得到了其解析解。计算结果表明:固结压力一定时,温度升高加快土体固结;温度一定时,土体中孔压随固结压力增大而减小;改进西原模型的弹性模量越大,孔压消散速率越大,土体固结越快;土颗粒受温度影响而产生的热膨胀对于孔压消散的作用非常小;黏滞系数越小,孔压消散更快,土体固结越快;试验结果与考虑土体的黏弹塑性的热力耦合改进西原模型计算结果吻合,改进西原模型可以较好地描述土体的热力耦合特性。     

基于相似性原理超固结土不排水扩张半解析解

在实际工程中,土体往往因卸载、再加载等复杂应力路径而处于超固结状态,而现有的圆孔扩张问题的计算模型往往不能反映超固结土中剪胀、软化等一些特殊性质。为了解决这一问题,基于相似性原理和统一硬化（UH）模型,结合相关联的流动法则和大变形理论,采用相似求解技术求解了超固结土不排水扩张问题的半解析解答。通过理想化算例分析了圆孔扩张挤土产生的应力和孔压响应,并通过分析不同超固结比OCR的土体应力路径的变化规律,讨论了UH模型的适用性。结果表明：对于轻超固结土,空腔周围土体孔压在塑性区沿径向单调递减,随着OCR增大,塑性区内孔压分布呈现出“S”形的趋势,孔壁附近的孔压逐渐减小,孔壁周围甚至出现负孔压。随着OCR增大,压力?扩张曲线收敛变慢。在扩孔过程中正常固结土一直处于剪缩硬化阶段。而对于超固结土,土体则经历了临界状态→剪胀硬化阶段→临界（特征）状态→剪缩硬化阶段。该研究成果不仅丰富了相似求解技术的应用,而且为超固结土中桩基承载力、隧道围岩变形预测和原位测试参数等岩土工程问题的计算提供了理论依据。     

黄河水下三角洲硬壳层特征及其液化过程研究

为揭示黄河口水下三角洲硬壳层的土性特征和风浪作用下的液化破坏状况,选择典型研究区,在现场利用普氏贯入仪测试硬壳层的强度特征,原位取1m原状样进行室内土工试验;利用重锤锤击荷载板的方式模拟波浪对硬壳层的动力作用,通过孔压探头和普氏贯入仪实时监测土体内孔压和振动前后强度的变化;通过理论计算研究硬壳层在不同风浪等级下的液化深度。通过研究发现,（1）硬壳层土体基本上处于超固结状态,且超固结比随深度增加而减小,大王北和刁口地区硬壳层强度约是新滩和广利港的两倍,离河口近的地区强度的变异系数比远的地区要大;（2）根据孔压随振次的变化关系,现场土体在振动过程中,孔压的增长经历了4个阶段,即初始阶段、增长阶段、稳定阶段和衰减阶段,且表层土体达到液化,深层的未液化;（3）大王北硬壳层在6~10级风浪下的平均液化深度仅为7~11cm,新滩和广利港硬壳层在6~10级风浪作用下的液化深度达32~42cm。强度高的硬壳层液化深度小,低的液化深度大,这种液化深度的差异性造成了地貌上的凹凸不平。     

冲击载荷作用下饱和软粘土孔压增长和消解规律

通过大量室内试验,研究饱和软粘土在冲击荷载作用下的变形和孔压增长规律,讨论孔压增长的内在机理和不同围压、不同冲击能作用下孔压量的变化。讨论土体的再固结变形规律,提出再固结体积压缩系数的确定方法。最后分析冲击荷载作用后的强度变化。这些讨论为动静结合法处理软基提供了依据。     

循环荷载作用下超固结软黏土软化-孔压模型研究

随着循环次数的增加,循环荷载作用下正常固结软黏土孔隙水压力逐渐增加,而孔压的增加将导致土体发生软化现象。对于超固结软黏土,循环荷载作用下在循环初期将产生负孔压,而此时也发生循环软化现象,这显然与有效应力原理是相矛盾的。以往的研究往往通过建立软化指数与循环次数的关系来描述土体的循环软化特性,从而不能反映土体的循环过程中残余孔压变化对土体循环软化的影响。通过对杭州饱和软黏土进行应力控制的循环三轴试验,对循环荷载作用下超固结软黏土的软化特性及孔压发展规律及两者关系进行了研究。在试验的基础上建立了超固结软黏土循环软化-孔压模型,该模型反映了残余孔压增长对超固结软黏土循环软化特性的影响规律。     

现场软土欠固结程度孔压消散试验评价方法

原位测试技术分析软土的应力历史可以避免取样及室内土工试验对土样的扰动,其结果能真实地反映现场土体的工程特性。以往基于孔压静力触探（CPTU）测试技术的超固结比（OCR）计算方法主要是针对超固结土取得的,缺乏对现场处于欠固结状态土体的考虑,具有一定的局限性。在已有研究成果的基础上,提出采用不完全孔压消散曲线的末段及时间平方根倒数外推法计算原位初始孔压,如果初始孔压大于静水压力,表明原位土层中存在固有孔压,为欠固结土。在此基础上,通过计算固结状态参数,可进一步对软土层的欠固结程度进行定量评价。工程应用表明,该方法不受软土性质和地域限制,具有普遍适用性,是合理可行的。在缺乏室内固结试验成果的情况下,可高效快捷地判别软土层的固结状态。研究成果对合理评价欠固结软土的工程特性具有一定的指导意义。     

波浪作用下饱和砂质海床土体与管线相互作用规律研究

针对波浪引起的饱和砂质海床土体和管线相互作用问题,将Biot动力固结理论与笔者课题组提出的砂土液化变形弹塑性本构模型相耦合,较为合理地再现了简谐波浪作用下较浅饱和砂质海床中管线周围可液化海床土体的超静孔隙水压力瞬态累积变化规律与液化过程。数值计算结果与Sumer等的试验规律一致。结果表明：由于管线的存在,改变了饱和砂质海床液化区域的空间分布。液化首先由管线下部土体开始产生,随着波浪荷载的持续作用,液化区域沿着管线外壁向上演化;同时海床表层土体产生液化并向深层发展,最终管线周围土体都发生液化,这是导致空管上浮的主要原因。当饱和砂质海床中存在管线时,管线附近海床土体液化深度明显变深。超静孔压累积和渗透力变化的耦合作用是导致饱和砂质海床土体产生液化的原因。与将海床土体视为饱和弹性多孔介质相比,可考虑液化全过程的弹塑性动力分析能更为合理地揭示实际波浪作用下饱和砂质海床土体的渗流场和应力场的瞬态时空演变规律。     

考虑固结的透水管桩沉桩全过程有限元模拟

运用透水管桩技术加快沉桩后桩周土体内超静孔隙水释放,进而消除沉桩施工中超静孔压的不利影响。基于有限元数值计算法,利用ABAQUS有限元软件建立透水管桩模型,实现透水管桩贯入过程以及桩周土体固结过程模拟。对比CEM圆柱孔扩张理论验证数值计算结果;阐述透水管桩贯入过程中位移场以及超静孔压场变化规律;对比分析静压桩和透水管桩桩周土体固结性状,结果表明透水管桩能加速超静孔压消散,短期内实现桩基承载力的快速提升。     

基于应变空间多机构CG模型的地震液化大变形数值分析

砂土液化导致的地基侧向大变形是地震中许多重要的工程设施和建筑物破坏的主要原因之一。简要介绍了可进行液化大变形分析的散粒体材料本构模型--应变空间多机构CG模型,基于FLIP ROSE程序平台,建立了预测和研究倾斜地基砂土液化导致侧向大变形的二维有限元数值分析方法。采用该模型对相同工况的土工动态离心模型试验进行了模拟,通过对比超孔隙水压力、剪切波水平加速度以及地基侧向位移发现,数值预测与试验结果吻合良好,从而验证了该有限元数值分析模型的可靠性。最后利用该数值分析模型预测了倾斜率不同的地基受到相同剪切波作用时,倾斜地基不同深度产生的侧向位移。预测结果显示,随着地基深度的减小,倾斜率对于地震液化导致倾斜地基侧向大变形的影响越来越显著。     

Assessment of soil liquefaction potential based on numerical method- a case study of an urban railway project

Paying special attention to geotechnical hazards such as liquefaction in huge civil projects like urban railways especially in susceptible regions to liquefaction is of great importance. A number of approaches to evaluate the potential for initiation of liquefaction, such as Seed and Idriss simplified method have been developed over the years. Although simplified methods are available in calculating the liquefaction potential of a soil deposit and shear stresses induced at any point in the ground due to earthquake loading, these methods cannot be applied to all earthquakes with the same accuracy, also they lack the potential to predict the pore pressure developed in the soil. Therefore, it is necessary to carry out a ground response analysis to obtain pore pressures and shear stresses in the soil due to earthquake loading. Using soil historical, geological and compositional criteria, a zone of the corridor of Tabriz urban railway line 2 susceptible to liquefaction was recognized. Then, using numerical analysis and cyclic stress method using QUAKE/W finite element code, soil liquefaction potential in susceptible zone was evaluated based on design earthquake.     

CONE模型与地基动力液化的非线性有效应力分析

在分析建筑物地基的地震（或振动）液化时,为了简化分析土壤-地基的动力相互作用,引入Cone模型,在总应力分析的方法基础上,结合了Martin等提出的孔隙水压力发展模型,并且考虑土壤发生液化的动力非线性过程中材料刚度退化的影响,提出了可以直接求解各区域的孔隙水压力的有效应力简化分析法,并进行液化可能性判别。     

循环荷载作用下海床结构粉质土的液化渗流机理定性研究

通过对黄河三角洲海床粉质土的基本性质和循环荷载作用下的喷水冒砂现象的分析,认识了海床粉质土的结构性。利用岩土破损力学原理,建立了海床结构粉质土的二元介质模型。定性分析了结构带土体在加荷、卸荷、液化渗流过程中的变形破损和在液化渗流过程中有效应力和超孔隙水压力的相互转化,指出在液化深度内,土体恢复后的固结强度应提高。这与试验结果相一致。用二元介质模型对解释试验中液化渗流差异是结构带内出现喷水冒砂和结构体上出现遍地冒水现象。     

考虑波致孔压累积与应力耦合效应的埋置管线周围海床响应特征与机制分析

波浪导致的海床液化是埋置管线失稳的一个重要因素。超静孔隙水压力累积引起的液化深度较深,对海底管线稳定性的影响较大,因此波浪作用下管线-海床系统的累积响应特征一直受到研究者的重点关注。本文基于考虑孔压累积与海床应力耦合发展的数值计算模型,对非线性行进波作用下含埋置管线的海床累积响应特征进行了模拟计算,并与非耦合模型的计算结果进行了对比分析,结果表明,当考虑孔压累积与海床应力的耦合效应时,管线附近累积孔压在水平方向上的不均匀分布会导致海床循环剪应力的增大,从而会极大地促进管线周围海床累积孔压的发展,增大管线的影响范围; 忽略孔压累积与海床应力的耦合效应,会在一定程度上低估管线周围海床的液化深度,不利于管线的安全。     

粉质海床波浪响应的数值模拟及试验研究

针对波浪荷载作用下粉质海床的动力响应问题,结合波浪槽模型试验研究的结果,提出了一种粉质海床在波浪作用下超静孔隙水压力增长的数值分析方法。计算结果表明：该方法能较好地拟合试验数据,揭示粉质海床孔隙水压力发展的机理,能运用于海床液化势的评价分析。     

Effects of driving and subsequent consolidation on behaviour of driven piles

The effects and simulation of driving of structures (piles) into saturated soil media are discussed, and procedures for numerical simulation of driving are proposed. Consolidation caused by changes in stresses and in pore water pressures in the soil mass due to the driving is solved by using a finite element procedure. The changes in stresses and pore water pressures due to driving are obtained on the basis of the cavity expansion approach, and are introduced in the finite element procedure as initial conditions. Stresses and deformations around a pile as consolidation proceeds are plotted and related to the quantities such as a wall friction, and point and total loads relevant to analysis and design.     

Seismic response of pile foundations in liquefiable soil: parametric study

The performance of pile foundations in liquefiable soil subjected to earthquake loading is a very complex process. The strength and stiffness of the soil decrease due to the increase in pore pressure. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. This paper presents the response of vertical piles in liquefiable soil under seismic loads. A finite difference model, known as fast Lagrangian analysis of continua, is used to study the pile behavior considering a nonlinear constitutive model for soil liquefaction and pile?Csoil interaction. The maximum lateral displacement and maximum pile bending moment are obtained for different pile diameters, earthquake predominant frequencies, Arias intensities, and peak accelerations. It is found that the maximum lateral displacement and the maximum pile bending moment increase when the predominant earthquake frequency value decreases for a given peak acceleration value.     

Liquefaction potential of coastal slopes induced by solitary waves

Tsunami runup and drawdown can cause liquefaction failure of coastal fine sand slopes due to the generation of high excess pore pressure and the reduction of the effective over burden pressure during the drawdown. The region immediately seaward of the initial shoreline is the most susceptible to tsunami-induced liquefaction failure because the water level drops significantly below the still water level during the set down phase of the drawdown. The objective of this work is to develop and validate a numerical model to assess the potential for tsunami-induced liquefaction failure of coastal sandy slopes. The transient pressure distribution acting on the slope due to wave runup and drawdown is computed by solving for the hybrid Boussinesq—nonlinear shallow water equations using a finite volume method. The subsurface pore water pressure and deformation fields are solved simultaneously using a finite element method. Two different soil constitutive models have been examined: a linear elastic model and a non-associative Mohr–Coulomb model. The numerical methods are validated by comparing the results with analytical models, and with experimental measurements from a large-scale laboratory study of breaking solitary waves over a planar fine sand beach. Good comparisons were observed from both the analytical and experimental validation studies. Numerical case studies are shown for a full-scale simulation of a 10-m solitary wave over a 1:15 and 1:5 sloped fine sand beach. The results show that the soil near the bed surface, particularly along the seepage face, is at risk to liquefaction failure. The depth of the seepage face increases and the width of the seepage face decreases with increasing bed slope. The rate of bed surface loading and unloading due to wave runup and drawdown, respectively, also increases with increasing bed slope. Consequently, the case with the steeper slope is more susceptible to liquefaction failure due to the higher hydraulic gradient. The analysis also suggests that the results are strongly influenced by the soil permeability and relative compressibility between the pore fluid and solid skeleton, and that a coupled solid/fluid formulation is needed for the soil solver. Finally, the results show the drawdown pore pressure response is strongly influenced by nonlinear material behavior for the full-scale simulation.     

Numerical simulation of mitigation for liquefaction-induced soil deformations in a sandy ground improved by cement grouting

This paper presents a numerical study of mitigation for liquefaction during earthquake loading. Analyses are carried out using an effective stress based, fully coupled, hybrid, finite element-finite differences approach. The sandy soil behavior is described by means of a cyclic elastoplastic constitutive model, which was developed within the framework of a nonlinear kinematic hardening rule. In theory, the philosophies of mitigation for liquefaction can be summarized as two main concepts, i.e. prevention of excess pore water pressure generation and reduction of liquefaction-induced deformations. This paper is primarily concerned with the latter approach to liquefaction mitigation. Firstly, the numerical method and the analytical procedure are briefly outlined. Subsequently, a case-history study, which includes a liquefaction mitigation technique of cement grouting for ground improvement of a sluice gate, is conducted to illustrate the effectiveness of liquefaction countermeasures. Special emphasis is given to the computed results of excess pore water pressures, displacements, and accelerations during the seismic excitation. Generally, the distinctive patterns of seismic response are accurately reproduced by the numerical simulation. The proposed numerical method is thus considered to capture the fundamental aspects of the problems investigated, and yields results for design purposes. From the results in the case, excess pore water pressures eventually reach fully liquefied state under the input earthquake loading and this cannot be prevented. However, liquefaction-induced lateral spreading of the foundation soils can be effectively reduced by the liquefaction mitigation techniques. An erratum to this article can be found at     

Stresses around Pressure Tunnels with Semi-permeable Liners

Summary. The liner of a pressure tunnel needs to be designed such that it can withstand the loads from the ground, the internal pressure, and minimize the development of significant pore pressures at the liner-ground interface. Pore pressures behind the liner reduce the effective stresses in the ground immediately in contact with the liner and can ultimately produce loss of support from the ground. Deformations and loads of the liner are intimately connected to the interplay that exists between liner, ground, and pore pressures in the ground. A closed-form analytical solution has been derived that accounts for the inter-relation between liner, ground, and pore pressures. Elastic response of the liner and ground, and plane strain conditions at any cross-section of the tunnel are assumed. The solution shows that stresses in the ground depend on the following dimensionless factors: relative stiffness of the ground and liner, ground Poisson’s ratio, surface slope angle, coefficient of earth pressure at rest, relative tunnel depth, and magnitude of the pore pressure behind the liner relative to the internal pressure. The minimum ground effective tangential stresses at the ground-liner interface increase with the relative stiffness of the liner, with the coefficient of earth pressure at rest, and with tunnel depth. They decrease with increasing surface slope angle and pore pressures behind the liner. As leakage through the liner increases, the pore pressures in the ground increase. This results in a decrease of effective radial and tangential stresses in the ground while displacements and loads of the liner are relatively less affected.