复合地基三维动力固结过程的数值仿真.docxVIP

文章编号：

复合地基三维动力固结过程的数值仿真

张萌,王鹏林

（天津市同盟信息系统集成有限公司，天津300191）

摘要：基于有限元方法建立了复合地基固结过程的动力学模型，该模型包括桩-土相互作用、

土壤-孔隙流体材料等非线性特性。通过子模型和嵌入区域技术降低了计算规模、提高了收敛

性，建立了针对此类问题的基本分析流程和方法。

关键词：动力固结；桩-土相互作用；有限元方法；ABAQUS软件

中图分类号：TU473.1文献标识码：A

NumericalSimulationofDynamicConsolidation

ProcessofThree-DimensionalCompositeFoundation

ZHANGMENG,WANGPENGLIN

(TEAMInformationSystemIntegrationCo.,Ltd,Tianjin300191,China)

Abstract:BasedonFEA(FiniteElementAnalysis)methodtocreatedynamicsmodeloftransient

consolidationprocessofcompositefoundation.Thismodelincludessomenon-linear

characteristics,likepile-soilinteraction,soil-porefluidmaterialandsoon.Usesubmodeland

embeddedregiontechnology,toreducethesizeofcalculateandimprovetheconvergenceofthis

model.Analysismethodandprocessaimatthisproblemisestablished.

Keywords:TransientConsolidation;Pile-SoilInteraction;FiniteElementMethod;ABAQUS

Software

0引言

桩-土相互作用、岩土材料塑性本构、多孔介质的渗

随着我国基础建设能力的不断发展，许多与岩流场、半无限计算区域等方面的影响都需要详细计

土工程相关的问题需要考虑更周全的影响因素，进算。但是在有限元模型中同时包含这些非线性因素

行更合理的分析预测。例如复合地基沉降中涉及的将会引起收敛困难。所以在实际应用中，受到计算

2

资源的限制，往往会将其中一些条件忽略或者简化

???????（1-1）

u

f

为线性。而这些忽略或简化的条件对计算结果的影

式中：??——土骨架所承担的有效应力

?——系数，通常由试验确定。在ABAQUS

响却难以评估。

软件中假定??s，即饱和度。

u——孔隙水压力

f

因此通过一个典型复合地基固结过程的算例，

??

?111000?T

建立了针对此类问题的多步计算策略。可以在较小

因此孔隙水压力边界不精确，得到的应力结果的计算规模下，将前面提到的种种非线性因素纳入

将出现误差。在其他半无限计算域的处理方法中（例有限元模型当中。

1规划计算策略

1.1问题描述本例中路基的长度相对其宽度足

如：施加粘滞边界或边界元方法）无法精确的对应

力波或孔隙水压力消散过程进行模拟并且其中很多

够长，且相关的地质条件沿中线不变化。但是地基

采用CFG桩-盖板结构加强，导致路基沿中线的应

参数不易确定。

力分布不均匀，不能简化成二维平面应变模型，必

须建立如图1所示的三维模型。为了简化计算，只

1.3桩-土相互作用为了精确计算沉降量和地

取典型段的四分之一部分加以研究。

基承载力必须考虑桩-土之间的接触作用。有限元方

法处理接触时，需要接触面两侧的网格尺寸尽量一

致。接触单元的尺寸差异过大会导致计算误差增大，

收敛性下降。

由图1可见增大土体边界后，CFG桩与土体的

几何模型尺寸差异很大。如果让二者网格尺寸一致，

图1复合地基几何模型

则会导致网格数量激增或者圆桩的几何特征无法真

1.2半无限计算区域岩土工程问题中土体的边实反映。如果在土体的局部进行网格加密，则会导

界是趋于无限远的。使用有限元方法计算时，计算致土体内部网格尺寸差异过大或者网格过渡区域单

区域大小必须有所限制，因此需要将土体边界人为元质量下降，同样不可取。

截断。对于本例中土质较软的地基，可以使用增大1.4子模型技术为了在有限元模