苍梧县城区防洪排涝工程桩基础非线性有限元分析.pdfVIP

苍梧县城区防洪排涝工程桩基础非线性有限元分析 简介： 介绍采用 ANSYS 有限元分析软件分析苍梧县城区防洪堤边坡 及桩基整体稳定性。 关键字：非线性 有限元分析 桩基础 稳定性 1、前言 　　苍梧县城（龙圩镇）位于广西东部、珠江流域西江水系的浔江 与下小河交汇处，距梧州市约 12Km。城区防洪堤全长约4.5Km，堤基 处于一级阶地前沿缘 15m 高程线附近，阶地岸坡曲折，坡度变化较大， 地势相对较为挟窄。防洪堤基础岩土共分七层，在某些地段有缺失， 自上而下的覆盖层分别为①人工填土，②褐灰色粉质粘土，③黄色粘 土，④灰黑色粉质粘土，⑤砂卵石层，⑥花岗岩残坡积层，⑦花岗岩。 我们根据地质部门提供的岩土力学参数，应用圆弧法对岸坡稳定进行 了大量的计算工作，分析结果表明，在未采取任何工程措施情况下， 不少堤段的岸坡稳定安全系数小于 1，计算的结果为不稳定。为此， 我们设计采用挡水结构加抗滑桩工程措施。鉴于目前对抗滑桩的稳定 计算尚未有成熟的计算方法[1]。本文介绍应用ANSYS 有限元分析软 件，对桩与土的相互作用进行有限元分析，以判定防洪堤的稳定性。 　　2、计算资料 　　主要计算资料有： 　　a) 选取最危险滑动堤段断面，断面的土层及结构如图 1。 　　b) 土层及桩的物理力学性能参数（表一） 　　表一土层及桩的物理力学性能参数表 参数土层 E （MPa）μ（kN/m3）（kN/m3）C （kPa）φ（度）① 4.20.0917.014.56.06.5②5.00.1518.7214.314.511.5③ 7.30.319.315.126.015.5④5.120.218.914.811.510.0 花岗岩 26000.622.022.0 混凝土 260000.16725.025.0 　　桩基相关参数表 岩土名称类别①填土②粉质粘土③粘土④粉质粘土弱风化花岗岩人 孔工挖桩桩周土摩擦力 qs915302070 桩端土承载力 q 　　 　　图 1 分析模型剖面图 　　c).结构荷载 　　 (1)结构及土体自重； 　　 (2)路面活动荷载 6.86kPa； 　　 (3)水压力（计算到最高洪水及最低洪水位）； 　　 (4)上部结构传到桩台的力（表二）。 　　表二上部结构传到桩台的力 工况 M （kN.m）V （kN）H （kN）洪水期 16532134472972 退水期 18858173613340　　3、计算范围及简化 　　a)计算范围：以挡水结构中心线为基准，向河床方向取 24.7m，向街道方向取 25.3m，地基深度取至-8.9m 高程。 　　b)不考虑渗透动水压力，在浸润线以上取地基土的湿容重γ湿， 浸润线以下取饱和容重。 　　c)边界条件：左右边界水平约束，下端固定。 　　4、力学模型 　　a).土与桩共同作用下的非线性有限元分析 　　本报告旨在对桩与土的共同作用的力学性能进行分析，必须考 虑土的非线性及边坡的滑移可能出现的大变形，因此对防洪堤进行非 线性分析时要考虑土材料的非线性及几何非线性。采用 Drucker-Prager 屈服准则。该准则对库仑——莫尔准则予以近似， 由此来修正Von Mises 屈服准则，该准则适用于混凝土、岩石、土壤 等材料。相应的等效应力的表达式为： 　　=3+[ 　　其中： {S}为偏应力，[M]为 Mises 准则中的[M]， 　　=，Φ为内摩擦角； 　　材料的屈服参数定义为： 　　其中： C 为粘聚力。 　　为了使计算结果更准确地反映结构的实际变形，进行有限元分 析时，在土和桩之间采用双弹簧接触单元以反映桩与土的相互作用， 桩和土分别剖分为 6 节点三角形单元。计算过程中逐级加载，按照 Newton-Raphson 迭代算法和线性搜索技术以保证计算的收敛性。 　　b).平面应变分析模型 　　为了得到较为精确的平面问题的解，对每根桩的刚度进行等效 简化，把原来圆截面桩等效为同厚的矩形截面桩，则等效后转化为等 截面的柱体（图 2），可以按平面应变问题来分析。设原来桩的弹性 模量 E0，惯性矩为 I0，等效后分别为E1，I1；则 　　 　　图 2 桩刚度沿长度等效简化图 　　E0I0=E1I1 　　= 　　则等效后桩的等效弹性模量为：E1=0.147E0 　　5、计算结果分析 　　 (a)桩的典型截面应力