py曲线法的介绍及案例分析PPT课件

1、2021年12月15日1. 1.引言引言第1页/共36页p-y曲线法：2021年12月15日第2页/共36页2 2. .研究现状研究现状2021年12月15日第3页/共36页2021年12月15日p-y曲线法的概念最早由Mcclelland和Focht提出来。他们认为试桩的实测反力与变位的关系曲线与同时进行的土的固结不排水三轴试验应力应变曲线存在相关关系，于是提出了一种求解桩非线性横向阻力的方法。国内国内对p-y曲线法的研究较晚，韩理安教授根据众多的现场试桩资料提出了p-y曲线的，采用相似理论的计算提出了一种构造p-y曲线的简便方法，该法目前已纳入港口工程桩基规范。第4页/共36页2021年1

2、2月15日从目前p-y曲线的研究方法来看，可分为三类：2.土体为，长度3 m，外径10mm，壁厚5mm的 钢管桩模型试验3.土体为，钢管桩模型试验4.土体为，钢管桩模型试验最为可靠，但其成本较高，此外，现场试验具有地域局限性。具有成本低、条件易控制等优点，不足之处在于土体为重塑土，无法反应现场原状土的结构性，此外，模型与实际的一致性目前还难以评价。 优点在于重复性好、条件易控制，最大问题在于土体的本构模型及桩土界面的模拟。第5页/共36页2021年12月15日土体一般分为沙土和黏土，那么，p-y曲线的研究也应分别考虑。沙土的研究1.Reese等基于埋入砂土中一组钢管桩试验研究和分析，从土楔体极

3、限平衡理论推出了砂土的极限土抗力，建立了；2.美国API规范对上述公式进行了修正，修正后的方法计算位移精度较高，对于控制侧向位移为主而桩自身强度较高的水平承载桩，新法尤佳；3.李雨润、王成雷等利用振动台试验对砂土液化条件下的桩土动力相互作用p-y曲线进行了分析并提出了相应的修正计算方法。黏土的研究4.河海大学新统一法；3.Sullivn法；5.同济大学法。2.Reese法；1.Matlock法；五种方法的p-y曲线构造关键在于函数分段依据和极限土抗力的确定。关键第6页/共36页2021年12月15日3 3. .非线性非线性有有 限元分析限元分析第7页/共36页NewmarkNewmark法法N

4、ewmark法的基本原理 把桩划分成若干段，将沿每段桩侧土体的横向抗力变换为一等效的弹簧支座在该段，弹簧的弹性系数则根据该处桩侧土的特性而定，从而所有问题都化为解算支承在一系列弹簧支座上的连续梁，用Newmark数值计算法求出梁的内力及位移。如下图所示。2021年12月15日第8页/共36页Newmark法公式推导 设桩入土长度为h，单元总数为n，单元长度L=h/n .由于K(z)=m (x+xn)n，x0为地面处与土的黏聚力有关的系数，则第i个弹簧的弹性系数为:1/21/2( )iiiKBK x dx2021年12月15日由式（1）积分可得：(1)111000/ 21nnm BKxxn110

5、0/ 21nnjnm BKxhxhn11001/ 21/ 2,1,2,3.11nnjim BKxixiinn(2)(4)(3)第9页/共36页2021年12月15日 基于以上Newmark法弹簧设置计算方法，本文将其推广用于p-y曲线方法的非线性弹簧的设置。其求解基本步骤如下：1.根据土性参数大概取定各层土的m值，通过Newmark方法计算各个相应位置弹簧的弹性系数，将此值作为弹簧弹性系数的初值；基于marc有限元软件计算梁弹簧结构体系的位移、转角，将此位移、转角作为梁的初始位移和转角;2.根据第1步计算的位移、转角，采用有限元插值函数求出任一入土深度x处的水平位移 ;1234izijzjYN

6、 YNN YN3.由X、Y值，计算桩上入土深度x处的水平土抗力 ，从而可得任意点处的地基系数 ; ,/K xP x YY4.设桩入土长度为h，单元总数为n，单元长度为 ，则可得第i个弹簧的弹性系数为(1/2)(1/2)( )iiiKK x dx第10页/共36页2021年12月15日有限元有限元分析分析法法 以上公式均是针对单桩结构物的公式分析方法，仅仅能够分析单根桩结构。那么，对于多桩结构物，我们也同样可以依据弹簧的设置来类似考虑，分析过程中但要考虑群桩效应的影响，即引入群桩效应系数。基于p-y曲线方法的有限元分析过程实际上是一个不断改变弹簧弹性系数的。2021年12月15日第11页/共36

7、页2021年12月15日具体分析过程如下：1.用mark有限元软件建立桩土作用的梁一弹簧结构体系分析模型；2.将模型初始弹簧弹性系数设为Newmark法所得出的弹簧系数，通过mark计算得出桩的位移和转角；3.将第2步位移作为初始位移Y0，通过式5可以得出下一次计算所需的弹簧弹性系数，再通过mark计算桩的位 移 、转角、内力等；5.通过mark结果文件得出桩结构的位移、转角和应力等。4.检验 (为一较小正数)是否成立，若不成立，则Y0 =Y1，返回第3步直到收敛条件满足为止， 执行下一步骤；2102YYEPS局限性：对于比较重要的工程，尤其当土层较多时参数也较多，不可能对所有参数进行反演，

8、该种方法就失灵了，因此，采用敏感性分析就可以避免此种问题。第12页/共36页2021年12月15日敏感性敏感性 分析分析 敏感性分析是指从定量分析的角度研究有关因素发生某种变化对某一个或一组关键指标影响程度的一种不确定分析技术。其实质是通过逐一改变相关变量数值的方法来解释关键指标受这些因素变动影响大小的规律。2021年12月15日 因此只需对桩位移、内力等敏感性较大的参数进行反演即可，本试验选取各参数对泥面位移的敏感性程度作为评判点。结论：参数对泥面位移的敏感性大小不仅与泥面位移对土体参数的偏导数有 关 ，而且还与参数的变异性有关，求解泥面位移对土体参数的偏导 数之后，即可得到土体参数对泥面位

9、移的敏感性的大小。第13页/共36页2021年12月15日实例分析实例分析桩的直径为1m，断面刚度EI=1.332467KPa，泥面以上桩长8.055m，泥面以下桩长35.6m。土有两层：表层为淤泥质亚黏土， ； ；第二层土为砂土， 。总共划分304个单元，其中泥面以上划分64个单元，共设置241个弹簧。32m=18kpa=18KN/muc厚，50=0.080333.6m=26=18KN/m厚，某码头预应力管柱桩试桩资料：第14页/共36页2021年12月15日弯矩沿桩长的变化规律图如下：第15页/共36页2021年12月15日计算结果与实测值对比表水平荷载（水平荷载（KN）比较项目比较项目实

10、测值实测值文文12本文本文40泥面位移（mm）6.3606.3046.335最大弯矩（KNm）372.0384.4385.8最大弯矩位置（m）2.6702.60060泥面位移（mm）11.7010.7910.936最大弯矩（KNm）635.8601.8621.8最大弯矩位置（m）3.1103.050 本文计算的位移值与实测值接近，最大弯矩相差稍大一些。总的来说，泥面位移及最大弯矩均与实测值吻合较好，这说明在横向荷载作用下考虑桩土之间的非线性相互作用的计算是有意义的，亦说明本文方法是合理可行的。第16页/共36页2021年12月15日 在进行敏感性分析时所有参数的变异系数均取为0.1，此时采用其

11、中某一参数值变化而引起泥面位移的改变作为敏感性大小程度的判据，计算结果见下表。黏土参数黏土参数泥面位移对参数的偏导数泥面位移对参数的偏导数参数标准值参数标准值泥面位移标准差值（泥面位移标准差值（mm）-5.62x10-21.80 x1001.01x10-1-1.45x10-21.80 x1002.61x10-25.47x10-28.00 x10-34.37x10-2-2.60 x10-24.50 x10-21.17x10-2-1.94x10-21.80 x1003.50 x10-3-6.35x10-22.60 x1001.65x10-4()uc kpakpa（）bkpa（）0( )()uckp

12、akpa（）b0( )kpa（） 可以看出：黏性土的不排水抗剪强度 和砂土的内摩擦角够对泥面位移的敏感性较大，砂土的重度对泥面位移的敏感性最小。第17页/共36页2021年12月15日4 4. .有限元有限元模型分析模型分析第18页/共36页2021年12月15日1.土体参数：内摩擦角 ，密度 ，弹性模量 ，泊松比 ，采用50mx50mx50m的立方体土体，位移边界，底面和四周边界加x、y、z向约束。上边界为地表面，为自由边界。0=313=2000kg/msE =100MPa=0.32.桩参数：为混凝土桩，密度 ，弹性模量 ，泊松比 ，桩半径 。 3=2400kg/m4sE =3 10 MPa

13、=0.2r=1m3.边界条件：底面和四周边界加x、y、z向约束。上边界为地表面，采用自由边界。 第19页/共36页2021年12月15日分析步设置分两个分析步模拟：第一步施加重力，获得初始地应力场，并通过地应力平衡将位移场清零；第二步施加水平推力。模型图如下：第20页/共36页2021年12月15日 分别对混凝土桩施加密度为2400kg/m3的重力，土体施加密度为2000kg/m3的重力 后，初始应力分布云图如下所示。初始应力分布云图S11应力分布云图第21页/共36页2021年12月15日S11应力分布云图组合体应力分布云图在X正向施加水平推力F=4KN后，应力分布云图如下。第22页/共36

14、页2021年12月15日S11应力分布云图混凝土桩应力分布云图第23页/共36页2021年12月15日在X正向施加水平推力F=4KN后，应变分布云图如下。S11应力分布云图组合体应变分布云图第24页/共36页2021年12月15日S11应力分布云图混凝土桩应变分布云图第25页/共36页2021年12月15日混凝土桩基地下部分节点的应力应变数据如下表所示。位置位置结点号结点号LE11LE11（10-6mm）LE11LE11（10-6mm）LE11LE11（10-6mm）LE11LE11（10-6mm）S11S11( (kPakPa) )S11S11( (kPakPa) )S11S11( (kPa

15、kPa) )S11S11( (kPakPa) )=1.0=0.6=0.8=1.2=1.0=0.6=0.8=1.20.00 4.00 631.00 631.00 631.00 631.00 0.00 0.00 0.00 0.00 -5.00 5.00 274.00 164.40 219.20 328.80 53.60 32.16 42.88 75.04 -10.00 8.00 41.90 25.14 33.52 50.28 175.60 105.36 140.48 245.84 -15.00 10.00 -19.18 -11.51 -15.34 -23.02 150.30 90.18 120.2

16、4 210.42 -20.00 12.00 -48.25 -28.95 -38.60 -57.90 116.20 69.72 92.96 162.68 -25.00 14.00 -90.40 -54.24 -72.32 -108.48 57.80 34.68 46.24 80.92 -30.00 16.00 -132.00 -79.20 -105.60 -158.40 15.80 9.48 12.64 22.12 第26页/共36页2021年12月15日混凝土桩x向水平位移分布曲线图（纵坐标表示距地面深度，横坐标表示水平位移）如下：混凝土桩x向水平位移分布曲线图第27页/共36页2021年12

17、月15日混凝土桩应力分布曲线图（纵坐标表示距地面深度，横坐标表示应力大小）如下：混凝土桩应力分布曲线图第28页/共36页2021年12月15日混凝土桩应力-位移分布曲线图（纵坐标表示应力，横坐标表示水平位移）如下：混凝土桩应力分布曲线图第29页/共36页2021年12月15日2021年12月15日5 5. .存在存在的问题的问题第30页/共36页2021年12月15日1. 基于少数现场试验结果，加上脱离实际场地条件作的纯数学 上的调整而引入的经验参数，理论基础不够严密；2.将桩土间的相互作用用一系列离散的一维非线性弹簧来模拟，忽 视了系统的空间连续性；3.无法考虑桩体抗弯刚度、桩截面形状、桩顶约束情况等一些 影响因素。P-y曲线法先天存在的不足之处：2021年12月15日第31页/共36页2021年12月15日2021年12月15日1. 没有区分桩体施工中是否存在挤土效应；2. 现场试验没有考虑桩体实际还承受着竖向荷载；3. 土力学室内试验应力路径没有与现场桩侧土体所处的真实应力状态对应。P-y曲线法研究中普遍被忽视的还有以下问题：第32页/共36页2021年12月15日6 6. .结语结语第