输电工程管理中的微型桩杆塔基础有限元数值模拟_陈红发.pdf

98电力与电子技术Power&Electronical Technology电子技术与软件工程Electronic Technology&Software Engineering微型桩是近年来逐渐发展起来的一种新型桩基础，又称树根桩。微型桩采用二次注浆，一般其直径小于 400mm，长细比大于 30，承载力较大1。其主要特点有施工机具小，适用于狭窄的施工作业区；施工振动、噪音小，且对土层适应性非常强 2。微型桩的概念最早于上个世纪 50 年代由意大利的 Lizzi提出，1982年由法国Soletanche公司将微型桩技术引入国内，引起国内学术界的广泛重视3。在国内主要应用于古建筑加固纠偏、房屋和防洪堤坝加固、边坡加固、旧房改造、建筑物加固防震抵抗交替荷载的基础、水池底板的抗浮等工程。近年来，微型桩逐渐在杆塔基础中得到应用，龚健等在软土地基中进行大量微型桩单桩和群桩现场原型试验，结果表明微型桩基础在软土地基杆塔基础中的应用具有足够的安全性能4。苏荣臻等对黄土地基微型桩杆塔基础施工工艺进行了现场原型试验5，郑卫锋等在青海地区进行了微型桩现场试验，针对寒冷地区施工问题进行了分析6。除此之外，中国电力科学研究院与各相关单位协同工作，先后在上海、天津等地进一步试点研究并获得了较好的经济和社会效益。运用 Abaqus 程序对单桩进行了抗压、水平力的数值计算分析，研究了微型桩的受力特性和工作特性，以期对微型桩基础的设计计算进行优化，进一步推进微型桩基础在输电线路中的应用。1 微型桩应用中存在的主要问题1.1 问题现状及分析微型桩作为一种新型杆塔基础，已证实其具有一定的适用性和应用价值。但是，当前研究均是在软土和黄土土基中对其性能表现进行试验研究。而宁夏地处我国西北部，气温多变，地质条件以砂土为主，其中风沙地占宁夏总面积的11.5%，表 1 为宁夏地区某 35kV 线路工程地勘报告。如若在宁夏地区使用微型桩，存在以下问题并对其进行分析：（1）对砂土土基的性能进行试验分析。国内学者对于微型桩在砂土中承载能力的相关研究甚少7，而砂土土体性质与软土、黄土土体性质完全不同，需要针对砂土特征进行相应的原型试验和数值模拟，以得到微型桩在砂土中的受力特性。（2）数值模拟的参数选择困难。在实际工程应用中，重要参数的选取一般通过现场原型试验得出或者参照其他工程的经验取值。因此，微型桩若要在砂土中进行有效应用，必须验证微型桩相关参数，可采取数值模拟与试验结果对比的方法验证参数。吴利霞等人对微型桩基础的承载性能进行了数值模拟分析，结果表明在有限元模型的基础上进行安全强度分析是可行的8-12，但是均未考虑砂土土基情况，因此对于砂土土基中有限元模型的可靠性，需要进一步研究验证。1.2 主要影响因素一般情况下，桩基础的承载能力受土基和桩体本身的共同影响。影响土基性能的因素包括弹性模量、泊松比、含水率、内摩擦角等以及桩-土相互作用的影响，如接触面和摩擦系数等，桩体本身性质包括弹性模量、泊松比、桩长、桩体直径等。例如内摩擦角直接影响桩-土之间摩擦系数，弹性模量和泊松比对土体自身的承载力影响较大，图 1 为某输电线路基础、地脚螺栓施工现场。对于有限元模型而言，本构模型的选取和参数的设置影响较大。常用土体本构模型有剑桥模型、Druker-Prager 模型、Mohr-Coulomb 模型、Lade-Duncan 模型等。剑桥模型常用于不排水固结饱和粘土中，Lade-Duncan 模型常用于无粘性土，Mohr-Coulomb模型适用于多种土体，适用性更加广泛13。参数的设置往往直接影响到模型结果的准确性和可靠性，参数选取、初始条件设置、边界条件设置、荷载设置均需引起注意。例如土基渗透性好坏、加载速度快慢会导致孔压上发生不同程度的上升，反映了不同排水状况下的短期或输电工程管理中的微型桩杆塔基础有限元数值模拟陈红发1吕鸿斌1霍历炜1李燕2白钰3（1.宁夏天净元光电力有限公司 宁夏回族自治区银川市 750004）（2.国网宁夏电力经济技术研究院 宁夏回族自治区银川市 750004）（3.国网银川供电公司 宁夏回族自治区银川市 750004）摘要：本文针对微型桩在砂土土基中的受力情况，利用 Abaqus 有限元分析软件，对桩在纵向和水平荷载作用下的受力情况和形变、位移等特征进行了研究。模拟结果分析揭示了桩身在土体中的受力、位移分布情况，并且水平荷载作用下，桩身 1/4 处弯矩最大。将力学模型计算结果与有限元模型结果进行对比，发现其误差在可接受范围内，从而验证了有限元模型的正确性，为工程项目设计和工程项目管理提供理论支持。关键词：微型桩基础；纵向荷载；水平荷载；有限元分析图 1：某输电线路基础、地脚螺栓施工现场 99电力与电子技术Power&Electronical Technology电子技术与软件工程Electronic Technology&Software Engineering表 1：35kV 某线路工程地勘报告序号（杆号）地貌描述岩土特性描述地基土主要力学指标地下水位（m）技经指标名称层底深度（m）岩性描述成因类型重度(kN/m3)粘聚力C(kPa)内摩擦角()承载力特征值fak(kPa)桩极限侧阻力标准值qsik(kPa)桩极限端阻力标准值Qpk(kPa)#1（J1）丘陵荒地黄土状粉土2.6褐黄色，稍湿，稍密，局部可见针状孔隙，垂直节理较发育，土质不均匀，局部砂性较强，偶含砾。Q4eol1510151103013.0丘陵采用干作业挖孔桩基础（机械施工），结合地质分层进行算费。角砾11.0杂色，稍湿，骨架颗粒一般粒径 2 20mm，最大粒径 30mm，次棱角状，骨架颗粒约占总重的 60%。母岩成分为砂岩、灰岩，充填粉细砂。处中密密实状态。Q4al+pl190352001352500卵石13.0杂色，稍湿，骨架颗粒一般粒径 20 50mm，最大粒径 350mm，次圆状，骨架颗粒约占总重的 65%。母岩成分为砂岩、灰岩，充填粉细砂。处中密密实状态。Q4al+pl200403001603000#2（G2）丘陵荒地黄土状粉土2.4褐黄色，稍湿，稍密，局部可见针状孔隙，垂直节理较发育，土质不均匀，局部砂性较强，偶含砾。Q4eol15101511030丘陵采用干作业挖孔桩基础（机械施工），结合地质分层进行算费。卵石6.5杂色，稍湿，骨架颗粒一般粒径 20 50mm，最大粒径 350mm，次圆状，骨架颗粒约占总重的 65%。母岩成分为砂岩、灰岩，充填粉细砂。处中密密实状态。Q4al+pl20040300160300011.0角砾11.0杂色，稍湿，骨架颗粒一般粒径 2 20mm，最大粒径 30mm，次棱角状，骨架颗粒约占总重的 60%。母岩成分为砂岩、灰岩，充填粉细砂。处中密密实状态。Q4al+pl190352001352500100电力与电子技术Power&Electronical Technology电子技术与软件工程Electronic Technology&Software Engineering长期承载力。建立模型时简化条件的选取会影响建模结果。在处理有限元模型问题时，往往需要进行一定的条件简化，这是由于实际问题的复杂度过高，模型建立困难，另外复杂的荷载和边界条件，往往会使模型结果不收敛，得不到有效结果。2 模型计算理论通过有限元方法，能够很好的模拟桩和土的几何形状、本构模型、边界条件等。为验证有限元计算结果的可靠性，需要与实测结果相互比较、相互印证，以便更好的理解桩的工作形状。计算桩基承载力的方法主要主要有 和 两种方法。其中粘土不排水条件下的桩适用于 方法。而 方法是基于有效应力的分析方法，既适合于粘土又适合于砂土，根据本文所选用的土质条件，故选用 方法来做理论分析14,15。其极限摩阻力为：fs=h （1）式中的 h 为作用在桩壁上的水平有效土压力；为桩和土之间的摩擦系数。考虑 h=K0v，则有：fs=K0v=v （2）式中的 v 为竖向土压力，K0为水平土压力系数。其极限端阻力为：fb=cbNc+v b Nq （3）式中 cb为桩端土体有效粘聚力；v b为桩端土体竖向有效应力；Nc和 Nq为承载力系数。（4）Nc=Nq-1cot （5）式中 为控制桩端破坏面性状的角度，0.330.58。3 模型设置及结果分析模型中桩基础采用线弹性模型，即小变形情况下的非线性模型。砂土地基采用 Mohr-Coulomb 模型进行模拟，该模型为理想弹塑性模型，理想弹塑性材料在应力未达屈服时会发生弹性变形，达到屈服时会产生永久塑性变形且不断发展直至破坏13。3.1 竖直受力模型本文假设桩位于正常排水的砂土、无孔隙压力等。考虑到该模型的对称性，采用轴对称模型进行分析，为确保边界对分析区域的影响最低，以桩长 20 米、桩径 0.4m 为例，分析区域为桩端向下 20 米、水平方向 8 米。对于排水条件，模拟应将时间取得足够长，以保证孔压得到充分的消散。该模型简化了实际情况，不考虑孔压即由有效应力承担全部荷载。模型中土体为砂土，选用 Mohr-Coulomb 模型模拟，#3（J2）丘陵荒地黄土状粉土1.4褐黄色，稍湿，稍密，局部可见针状孔隙，垂直节理较发育，土质不均匀，局部砂性较强，偶含砾。Q4eol1510151103013.0丘陵采用干作业挖孔桩基础（机械施工），结合地质分层进行算费。卵石4.5杂色，稍湿，骨架颗粒一般粒径 20 50mm，最大粒径 350mm，次圆状，骨架颗粒约占总重的 65%。母岩成分为砂岩、灰岩，充填粉细砂。处中密密实状态。Q4al+pl20040300160300013.0卵石13.0杂色，稍湿，骨架颗粒一般粒径 20 50mm，最大粒径 350mm，次圆状，骨架颗粒约占总重的 65%。母岩成分为砂岩、灰岩，充填粉细砂。处中密密实状态。Q4al+pl20040300160300013.0图 2：荷载沉降曲线101电力与电子技术Power&Electronical Technology电子技术与软件工程Electronic Technology&Software Engineering土体参数如表 2 所示。桩采用线弹性模型，弹性模量 E=20GPa，泊松比=0.2。桩土摩擦系数为 0.42（即 tan(0.75)）。对于初始应力的设置，假设土体经过一维K0正常固结，则竖向初始应力 v0和水平初始应力 h0为：v0=z （6）h0=K0v0 （7）式中 K0为初始水平压力系数，因为该模型只考虑到竖直方向的作用，不考虑水平方向，故取(1-)=0.538。有初始平均应力，偏应力 q0=1-k0z。图 2 给出了荷载作用下沉降曲线。由于端阻力在总荷载中占比较大，当变形比较大时才会逐渐显现，故该模型中荷载沉降曲线变化平缓，极限荷载通常取沉降为 0.1 倍直径时对应的荷载约为 371.2kN。图 3 给出桩端附近的等效应力分布。取 0.455，有承 载 力 系 数 Nq=15.63，Nc=25.34，端 阻 力 为 Qb=cbNC+v bNqA=(0.125.34+8015.63)0.22=157.9kN。而 桩 的 极 限侧阻力按式（1）得 Qf=142.0kN。因而总的极限承载力为Q=Qf+Qb=299.9kN，与有限元的计算结果相差约 19%左右。不同时刻的桩侧摩阻力分布如图 4 所示。桩侧摩阻力自上而下的发展过程：上部的外摩擦阻力先到极限，荷载向下传递，下部外摩阻力开始变化，随后保持不变，直至最终外摩阻力呈线性分布。由于桩顶受压后侧向膨胀，导致桩端顶部土体的水平应力变大，相应的摩阻力也变大，致使桩端顶部的外摩阻力有局部的扩大趋势。由于桩端以下的破坏区向上滑动变形，迫使桩侧土体上移，造成桩端摩阻力最大值随着时间有上移的趋势。3.2 水平受力模型水平荷载作用下的桩-土体系可按照桩、土相对刚度的不同划分为两类工作状态和破坏机理，即刚性短桩和弹性长桩，前者因转动或平移破坏，后者因挠曲破坏。微型桩桩身刚度比较低，长细比大于30，柔性性能较好，可视为具有弹性长桩的破坏状态。由