CN112784332A - 基于K-m法的加固桩板墙桩前地基的桩身内力计算方法 - Google Patents
基于K-m法的加固桩板墙桩前地基的桩身内力计算方法 Download PDFInfo
- Publication number
- CN112784332A CN112784332A CN202011310962.9A CN202011310962A CN112784332A CN 112784332 A CN112784332 A CN 112784332A CN 202011310962 A CN202011310962 A CN 202011310962A CN 112784332 A CN112784332 A CN 112784332A
- Authority
- CN
- China
- Prior art keywords
- pile
- section
- foundation
- pile body
- node
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 93
- 238000004364 calculation method Methods 0.000 title claims abstract description 69
- 238000006073 displacement reaction Methods 0.000 claims abstract description 59
- 238000005452 bending Methods 0.000 claims abstract description 39
- 230000003014 reinforcing effect Effects 0.000 claims abstract description 35
- 230000002787 reinforcement Effects 0.000 claims abstract description 8
- 238000006243 chemical reaction Methods 0.000 claims description 55
- 239000002689 soil Substances 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000010008 shearing Methods 0.000 description 8
- 239000011435 rock Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000004873 anchoring Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Bulkheads Adapted To Foundation Construction (AREA)
Abstract
本发明公开了一种基于K‑m法的加固桩板墙桩前地基的桩身内力计算方法,包括以下步骤:确定桩长和截面尺寸;确定加固端和未加固端的相关参数;确定嵌固段桩顶处作用的弯矩M0和剪力Q0;对于加固段采用有限差分法的K法进行迭代计算求得各节点系数;对于未加固段采用有限差分法的m法进行迭代计算求得各节点系数;计算出全桩各节点处的位移;通过迭代计算确定所有节点的内力值。本发明针对桩板墙工程中在桩前地基采用旋喷桩加固桩前地基,对板桩墙桩前滑动面下的嵌固段内力分为加固段和未加固段两部分进行计算,对于加固段采用有限差分法的K法进行计算,对于未加固段采用有限差分法的m法进行计算,使地基内力的计算结果更贴合于实际的计算结果。
Description
技术领域
本发明属于岩土工程领域,具体涉及一种基于K-m法的加固桩板墙桩前地基的桩身内力计算方法。
背景技术
桩板墙或排桩挡土墙桩基水平承载力计算中,锚固点以下的桩身变位和内力,应根据锚固点处的弯矩、剪力和地基的弹性抗力,采用地基系数法进行计算,根据岩土条件可选用“K法”或“m法"。当地基为较完整的岩层时,认为侧向地基系数是常数,不随深度而变化,此常数通常以“Kp”表示,相应的计算方法称为“K法”:当地基为土层或严重风化破碎岩层的情形时,认为侧向地基系数是随深度而呈直线规律的变化,地基系数随深度变化的比例系数通常以符号“m”表示,相应的计算方法称为“m法”。地基土水平抗力的比例系数m是个重要的计算参数,它的取值直接影响到桩基水平承载力计算结果。
工程中经常遇到地层为土层,且设置桩板墙或排桩挡土墙较为合理的工点,但是天然土层地基的地基比例系数较小,且横向承载力低,设置桩板墙或排桩挡土墙后,锚固点及以下的桩身变位和地基的横向抗力不能满足挡墙的正常使用状态或极限承载状态。为了提高桩前地基土的横向抗力、控制桩体位移,桩前土体可采用旋喷桩或搅拌桩等进行土体加固。为了提高设计的准确性、合理性,有必要对桩前地基加固措施和计算方法进行系统化研究。
发明内容
本发明的目的在于:针对于传统的板桩墙桩前地基桩身内力的计算方法,未考虑加固对土体变化的影响,故提供一种基于K-m法的加固桩板墙桩前地基的桩身内力计算方法。
本发明的技术方案是:一种基于K-m法的加固桩板墙桩前地基的桩身内力计算方法,包括以下步骤:
A确定桩长和截面尺寸
确定桩板墙结构尺寸以及桩的结构尺寸;
B.确定加固端和未加固端的相关参数
确定板桩墙嵌固端加固段和未加固段的长度,并确定加固段、未加固段交界面的位置及参数,确定地基反力系数K和地基反力系数与深度变化的比例系数m;
C.确定嵌固段桩顶处作用的弯矩M0和剪力Q0
利用悬臂段受力计算求得嵌固段桩顶处作用的弯矩M0和剪力Q0;
D.对于加固段采用有限差分法的K法进行迭代计算求得各节点系数
将加固段h1的桩体视作桩顶受水平集中荷载的侧荷桩,地基比例系数按K法考虑,得到加固段的挠曲微分方程,进而得到各节点系数;
E.对于未加固段采用有限差分法的m法进行迭代计算求得各节点系数
将未加固段h2的桩体视为桩顶受水平集中荷载的侧荷桩,地基比例系数按照m法考虑,考虑桩身的连续性,得到未加固段的挠曲微分方程,进而得到各节点系数;
F.计算出全桩各节点处的位移
确定得到矩阵方程,计算求出全桩各节点处的位移;
G.通过迭代计算确定所有节点的内力值
根据各内力与位移值的关系,通过迭代计算求出所有节点的内力值。
更进一步的,步骤D中加固段的挠曲微分方程如下:
B0为桩的计算宽度;K为地基反力系数。
更进一步的,步骤D中的各节点系数具体过程如下:
在h1段桩身进行有限差分计算时,任意节点n处的控制差分方程如下:
y0=a0y1-b0y2+c0
根据所述公式,通过递推法,得到含有参数yi、yi-1和yi-2的公式,将上述公式全部代入挠曲微分方程当中,整理得到公式如下:
(an-2an-1-bn-2-4an-1+6+B0Kh4/EI)yi+(an-2bn-1-4bn-1+4)yn+1-4cn-1+an-2cn-1+cn-2=0
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
根据公式和边界条件,以桩顶为起始点至交界面,通过迭代计算可求得加固段各节点的系数。
更进一步的,步骤E中未加固段的挠曲微分方程如下:
其中,E为桩体的弹性模量;I为桩体的横截面惯性矩;y为桩的侧向位移;Z为桩的埋深;B0为桩的计算宽度;K为地基反力系数; m为地基反力系数与深度变化的比例系数;h为桩身长度;h1为加固段长度;A为交界面处的地基系数,A=mh1。
更进一步的,步骤E中各节点系数具体过程如下:
h2段桩身任意节点n处的控制差分方程如下:
y0=a0y1-b0y2
根据所述公式,通过递推法,得到含有参数yn、yn-1和yn-2的公式,将上述公式全部代入挠曲微分方程当中,整理得到公式如下:
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;m为地基反力系数与深度变化的比例系数h为桩身长度;A为交界面处的地基系数,A=mh1;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
根据公式和边界条件所述,以桩底为起始点至交界面通过迭代计算求得未加固段各节点系数。
更进一步的,步骤F中矩阵方程具体如下:
[a]{y}=[c]
{y}=[yA-2 yA-1 yA yA+1 yA+2 y'A-2 y'A-1 y'A y'A+1 y'A+2]T
[c]=[0 0 0 0 cA cA-1 cA-2 0 0 0]T。
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
更进一步的,通过迭代计算就求出所有节点的内力值,
其中,加固段桩身各节点处的弯矩Mn、剪力Qn和地基土反力Pn,计算公式如下:
其中,Mn为桩身各节点处弯矩值;Qn为桩身各节点处剪力值; Pn为桩身各节点处地基土反力值;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度;yn为节点n处桩的侧向位移;Z为桩的埋深;B0为桩体的计算宽度;m为地基反力系数与深度变化的比例系数。
更进一步的,通过迭代计算就求出所有节点的内力值,
其中,未加固段桩身各节点处的弯矩Mn、剪力Qn和地基土反力 Pn,计算公式如下:
其中,Mn为桩身各节点处弯矩值;Qn为桩身各节点处剪力值;Pn为桩身各节点处地基土反力值;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度;yn为节点n处桩的侧向位移;Z为桩的埋深;B0为桩体的计算宽度;m为地基反力系数与深度变化的比例系数;A为交界面处的地基系数,A=mh1,其中,h1为加固段桩长。
本发明的有益效果如下:
本发明针对桩板墙工程中在桩前地基采用旋喷桩加固桩前地基, 对板桩墙桩前滑动面下的嵌固段内力分为加固段和未加固段两部分进行计算,对于加固段采用有限差分法的K法进行计算,对于未加固段采用有限差分法的m法进行计算,区别于传统的计算方法,充分考虑了桩前地基加固后对土体产生的影响,从而使地基内力的计算结果更贴合于实际的计算结果。
本发明结合旋喷桩加固板桩墙桩前地基的特点,具有参数易取、便于计算、计算结果精确等特点。具有很强的实用性和可靠性。
附图说明
图1是本发明中桩板墙嵌固段分析模型图;
图2是本发明中有限差分法计算简图;
图3是本发明中加固段桩身弯矩计算图;
图4是本发明中未加固段桩身弯矩计算图;
具体实施方式
以下,参照附图和实施例对本发明进行详细说明:
如图1-4所示,一种基于k-m法的加固桩板墙桩前地基的桩身内力计算方法,包括以下步骤:
A、确定桩板墙结构尺寸,桩的结构尺寸;
B、确定板桩墙嵌固端加固段和未加固段的长度,确定交界面位置及参数;确定地基反力系数K和地基反力系数与深度变化的比例系数m。
C、由悬臂段受力计算求得嵌固段桩顶处作用的弯矩M0和剪力 Q0。
将加固段h1的桩体视作桩顶受水平集中荷载的侧荷桩,地基比例系数按“K”法考虑时,其挠曲微分方程为:
其中,E为桩体的弹性模量;I为桩体的横截面惯性矩;y为桩的侧向位移;Z为桩的埋深;B0为桩的计算宽度;K为地基反力系数。
h1段桩身任意节点n处的控制差分方程如下:
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;M0为嵌固段桩顶处作用的弯矩;Q0嵌固段桩顶处作用的剪力。
在h1段桩身进行有限差分计算时,任意节点n处的控制差分方程如下:
y0=a0y1-b0y2+c0
根据上述公式,通过递推法,得到含有参数yn、yn-1和yn-2的公式,将上述公式全部代入挠曲微分方程当中,整理得到公式如下:
(ai-2ai-1-bi-2-4ai-1+6+B0Kh4/EI)yi+(ai-2bi-1-4bi-1+4)yi+1-4ci-1+ai-2ci-1+ci-2=0
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
根据所述公式和边界条件,以桩顶为起始点至交界面,通过迭代计算可求得加固段各节点的系数。
将未加固段h2的桩体视为桩顶受水平集中荷载的侧荷桩,地基比例系数按照“m”法考虑,考虑桩身的连续性,其挠曲微分方程如下:
式中,E为桩体的弹性模量;I为桩体的横截面惯性矩;y为桩的侧向位移;Z为桩的埋深;B0为桩的计算宽度;K为地基反力系数;h 为桩身长度;h1为加固段长度;A为交界面处的地基系数,A=mh1。
h2段桩身任意节点n处的控制差分方程如下:
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B0为桩的计算宽度;K为地基反力系数;m为地基反力系数与深度变化的比例系数h为桩身长度;h1为加固段长度; h2为未加固段长度;A为交界面处的地基系数,A=mh1。
在h1段桩身进行有限差分计算时,任意节点n处的控制差分方程如下:
y0=a0y1-b0y2
根据上述公式,通过递推法,得到含有参数yn、yn-1和yn-2的公式,将上述公式全部代入挠曲微分方程当中,整理得到公式如下:
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B0为桩的计算宽度;K为地基反力系数;m为地基反力系数与深度变化的比例系数;h为桩身长度;h1为加固段长度;h2为未加固段长度;A为交界面处的地基系数,A=mh1。
根据上述公式和边界条件所述,以桩底为起始点至交界面通过迭代计算求得未加固段各节点系数。
确定矩阵方程,求解得出交界面处及附近节点的位移,计算求出全桩各节点处的位移,矩阵方程如下所示:
[a]{y}=[c]
{y}=[yA-2 yA-1 yA yA+1 yA+2 y'A-2 y'A-1 y'A y'A+1 y'A+2]T
[c]=[0 0 0 0 cA cA-1 cA-2 0 0 0]T
[c]=[0 0 0 0 cA cA-1 cA-2 0 0 0]T。
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
求解上述方程可算出交界面处及附近节点的位移,并根据各内力与位移值的关系,通过迭代计算就可求出所有节点的内力值。
加固段桩身各节点处的弯矩Mn、剪力Qn和地基土反力Pn,计算公式如下:
其中,Mn为桩身各节点处弯矩值;Qn为桩身各节点处剪力值; Pn为桩身各节点处地基土反力值;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度;yn为节点n处桩的侧向位移;Z为桩的埋深;B0为桩体的计算宽度;m为地基反力系数与深度变化的比例系数。
未加固段桩身各节点处的弯矩Mn、剪力Qn和地基土反力Pn,计算公式如下:
其中,Mn为桩身各节点处弯矩值;Qn为桩身各节点处剪力值; Pn为桩身各节点处地基土反力值;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度;yn为节点n处桩的侧向位移;Z为桩的埋深;B0为桩体的计算宽度;m为地基反力系数与深度变化的比例系数;A为交界面处的地基系数,A=mh1,其中,h1为加固段桩长。
实际工程当中,抗滑桩的计算理论目前主要有静力平衡法、布鲁姆法和弹性地基梁法三种计算理论等。
优选的,弹性地基梁法计算简便且计算结果与实际结果吻合度很高,因此目前的桩板墙设计主要采用此法。
弹性地基梁法假定地基土抗力与位移成正比。
具体解法有三种:一、用数学方法计算桩体的弹性挠曲线微分方程;二、用有限差分法计算桩体的挠曲线微分方程;三、采用有限单元法,根据力的平衡和位移条件求解桩体内力和位移。
优选的,当采用有限差分法等数值分析的方法求出数值解,然后采用泰勒级数展开法把微分方程的基本方程组和边界条件改写成差分方程,然后进行求解,对于求解大变形的几何非线性问题。
步骤A中,桩板墙结构尺寸包括桩长L,桩的结构尺寸包括桩的宽度B,桩的计算宽度B0。
步骤B中,嵌固段的长度为h1,未嵌固段的长度为h2。所述物理学参数包括路K为地基反力系数;m为地基反力系数与深度变化的比例系数。
步骤B中,需要对嵌固段的长度进行换算,即:桩身的嵌固段长度与变形系数的乘积。引入桩的变形系数如下所示:
“K”法计算桩的变形系数:
“m”法计算桩的变形系数:
式中,β,α为桩的变形系数;K为地基反力系数;m为地基反力系数与深度变化的比例系数;E桩体的弹性模量;B0为桩体受荷面的计算宽度;I为桩体的横截面惯性矩;
假设桩全长为l,平均分成n等份,每段长为h,其中,桩顶节点为3,桩底节点为n+3。为了计算方便,设置多个虚拟节点辅助计算。
步骤E中,对于未加固段桩身内力计算中,本文采用“m”法进行计算土的抗力集度,公式为:p=-(A+mz)B0y。
其中:p为土的抗力集度;B0为桩体受荷面的计算宽度;y为桩的侧向位移;z为桩的埋深;m为地基反力系数与深度变化的比例系数;A为交界面处的地基系数,A=mh1,其中,h1为加固段桩长。
实施例一
确定桩板墙长度如下:
其中,确定桩的宽度B=3m;桩的计算宽度B0=2m;加固段长h1=9m,未加固段长h2=9m;加固段与未加固段相交处称为交界面。
桩的抗弯刚度EI=5.05×105KN/m;地基反力系数K=198000;地基反力系数与深度变化的比例系数m=2200;
由悬臂段受力计算得嵌固段桩顶处作用有弯矩M0和剪力Q0。
将加固段h1的桩体视作桩顶受水平集中荷载的侧荷桩,地基比例系数按K法考虑时,其挠曲微分方程为:
其中,E为桩体的弹性模量;I为桩体的横截面惯性矩;y为桩的侧向位移;Z为桩的埋深;B0为桩的计算宽度;K为地基反力系数。
在h1段桩身进行有限差分计算时,任意节点n处的控制差分方程为:
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B0为桩的计算宽度;K为地基反力系数;h为桩身长度。
由桩顶边界条件:
其中,yn为节点n处桩的侧向位移;z为桩的埋深;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度;M0为嵌固段桩顶处作用的弯矩;Q0嵌固段桩顶处作用的剪力。
整理得到:
式中, b-1=1、B0为桩的计算宽度;K为地基反力系数;h为桩身长度;yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;M0为嵌固段桩顶处作用的弯矩;Q0嵌固段桩顶处作用的剪力。
整理后得到下述公式:
y0=a0y1-b0y2+c0
则有yn=anyn+1-biyn+2+cn
则有yn-1=an-1yi-bn-1yn+1+cn-1
则有yn-2=an-2yn-1-bn-2yn+cn-2
将由递推得到的公式,代入挠曲微分方程整理后,得到下述方程:
(ai-2ai-1-bi-2-4ai-1+6+B0Kh4/EI)yi+(ai-2bi-1-4bi-1+4)yi+1-4ci-1+ai-2ci-1+ci-2=0
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
根据上述公式和边界条件,以桩顶为起始点至交界面,通过迭代计算可求得加固段各节点的系数。
将未加固段h2的桩体视为桩顶受水平集中荷载的侧荷桩,地基比例系数按照m法考虑,考虑桩身的连续性,其挠曲微分方程表示为:
其中,E为桩体的弹性模量;I为桩体的横截面惯性矩;y为桩的侧向位移;Z为桩的埋深;B0为桩的计算宽度;K为地基反力系数; m为地基反力系数与深度变化的比例系数;h为桩身长度;h1为加固段长度;A为交界面处的地基系数,A=mh1。
未加固段h2桩身任意节点n处的控制差分方程为:
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;n为节点编号;m为地基反力系数与深度变化的比例系数; h1为加固段长度;h2为未加固段长度;A为交界面处的地基系数, A=mh1。
桩底按照自由端考虑,即桩底处的弯矩和剪力均为0,得到下述方程:
当n=0时:y-1=2y0-y1
当n=0时:y-2=2y-1-2y1+y2
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度。
边界条件认为桩底为自由端,即桩底弯矩和剪力为0,
y0=a0y1-b0y2,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I 为桩体的横截面惯性矩;B0为桩的计算宽度;K为地基反力系数;h 为桩身长度;h1为加固段长度;h2为未加固段长度;A为交界面处的地基系数,A=mh1。
则有yn=anyn+1-bnyn+2
则有yn-1=an-1yn-bn-1yn+1
则有yn-2=an-2yn-1-bn-2yn
由递推得到的公式,代入挠曲微分方程整理后,得到下述方程:
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;m为地基反力系数与深度变化的比例系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
根据上述公式和边界条件,以桩底为起始点至交界面通过迭代计算求得未加固段各节点系数。
由于桩身是连续的,在交界面处满足桩身位移、转角、弯矩和剪力的连续条件。由连续条件得到10阶矩阵方程如下所示:
[a]{y}=[c]
{y}=[yA-2 yA-1 yA yA+1 yA+2 y'A-2 y'A-1 y'A y'A+1 y'A+2]T
[c]=[0 0 0 0 cA cA-1 cA-2 0 0 0]T
[c]=[0 0 0 0 cA cA-1 cA-2 0 0 0]T。
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
解上述矩阵方程可算出交界面处及附近共10个节点的位移,进而得到全桩各节点处的位移,根据各内力与位移值的关系,通过迭代计算就可求出所有节点的内力值。
加固段桩身各节点处的弯矩Mn、剪力Qn和地基土反力Pn,计算公式如下:
其中,Mn为桩身各节点处弯矩值;Qn为桩身各节点处剪力值; Pn为桩身各节点处地基土反力值;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度;yn为节点n处桩的侧向位移;K为地基反力系数;B0为桩体的计算宽度;m为地基反力系数与深度变化的比例系数。
未加固段桩身各节点处的弯矩Mn、剪力Qn和地基土反力Pn,,计算公式如下:
其中,Mn为桩身各节点处弯矩值;Qn为桩身各节点处剪力值;Pn为桩身各节点处地基土反力值;E为桩体的弹性模量;I为桩体的横截面惯性矩;h为桩身长度;yn为节点n处桩的侧向位移;Z为桩的埋深;B0为桩体的计算宽度;m为地基反力系数与深度变化的比例系数;A为交界面处的地基系数,A=mh1;h1为加固段长度。
为验证上述内力计算方法的合理性,使用进行有限元计算,由于将加固段类比于硬质岩,为对比加固段分别以K法和m法计算的结果,假设加固段地基比例系数为70MN/m4,未加固段取地基比例系数为5MN/m4进行计算。
如图3、图4所示,经过迭代计算后得到的桩身内力值变化如下,可以看出区别于传统的计算方法,充分考虑了桩前地基加固后对土体产生的影响,从而使地基内力的计算结果更贴合于实际的计算结果。对板桩墙桩前滑动面下的嵌固段内力分为加固段和未加固段两部分进行分步计算的方法,具有一定的实际工程指导价值。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于K-m法的加固桩板墙桩前地基的桩身内力计算方法,其特征在于:包括以下步骤:
(A)确定桩长和截面尺寸
确定桩板墙结构尺寸以及桩的结构尺寸;
(B)确定加固端和未加固端的相关参数
确定板桩墙嵌固端加固段和未加固段的长度,并确定加固段、未加固段交界面的位置及参数,确定地基反力系数K和地基反力系数与深度变化的比例系数m;
(C)确定嵌固段桩顶处作用的弯矩M0和剪力Q0
利用悬臂段受力计算求得嵌固段桩顶处作用的弯矩M0和剪力Q0;
(D)对于加固段采用有限差分法的K法进行迭代计算求得各节点系数
将加固段h1的桩体视作桩顶受水平集中荷载的侧荷桩,地基比例系数按K法考虑,得到加固段的挠曲微分方程,进而得到各节点系数;
(E)对于未加固段采用有限差分法的m法进行迭代计算求得各节点系数
将未加固段h2的桩体视为桩顶受水平集中荷载的侧荷桩,地基比例系数按照m法考虑,考虑桩身的连续性,得到未加固段的挠曲微分方程,进而得到各节点系数;
(F)计算出全桩各节点处的位移
确定得到矩阵方程,计算求出全桩各节点处的位移;
(G)通过迭代计算确定所有节点的内力值
根据各内力与位移值的关系,通过迭代计算求出所有节点的内力值。
3.根据权利要求2所述的基于K-m法的加固桩板墙桩前地基的桩身内力计算方法,其特征在于:步骤(D)中的各节点系数具体过程如下:
在h1段桩身进行有限差分计算时,任意节点n处的控制差分方程如下:
y0=a0y1-b0y2+c0
根据所述公式,通过递推法,得到含有参数yn、yn-1和yn-2的公式,将上述公式全部代入挠曲微分方程当中,整理得到公式如下:
(ai-2ai-1-bi-2-4ai-1+6+B0Kh4/EI)yi+(ai-2bi-1-4bi-1+4)yi+1-4ci-1+ai-2ci-1+ci-2=0
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
其中,yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度,
根据公式和边界条件,以桩顶为起始点至交界面,通过迭代计算可求得加固段各节点的系数。
5.根据权利要求4所述的基于K-m法的加固桩板墙桩前地基的桩身内力计算方法,其特征在于:步骤(E)中各节点系数具体过程如下:
在h2段桩身进行有限差分计算时,任意节点n处的控制差分方程如下:
y0=a0y1-b0y2
其中,y0为桩底的侧向位移;yn为节点n处桩的侧向位移;d0为桩的截面直径;M0为嵌固段桩顶处作用的弯矩;Q0为嵌固段桩顶处作用的弯矩;E为桩体的弹性模量;I为桩体的横截面惯性矩;m为地基反力系数与深度变化的比例系数;h为桩身长度;h1为加固段长度;h2为未加固段长度;A为交界面处的地基系数,A=mh1,
根据所述公式,通过递推法,得到含有参数yn、yn-1和yn-2的公式,将上述公式全部代入挠曲微分方程当中,整理得到公式如下:
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
yn为节点n处桩的侧向位移;E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;h2为未加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度,
根据公式和边界条件所述,以桩底为起始点至交界面通过迭代计算求得未加固段各节点系数。
6.根据权利要求5所述的基于K-m法的加固桩板墙桩前地基的桩身内力计算方法,其特征在于:步骤(F)中矩阵方程具体如下:
[a]{y}=[c]
{y}=[yA-2 yA-1 yA yA+1 yA+2 y'A-2 y'A-1 y'A y'A+1 y'A+2]T
[c]=[0 0 0 0 cA cA-1 cA-2 0 0 0]T,
其中,d(n)=a(n-2)*a(n-1)-b(n-2)-4*a(n-1)+6+B0*m*K*h^5/EI;
a(n)=(a(n-2)*b(n-1)-4*b(n-1)+4)/d(n);
b(n)=1/d(n);
c(n)=(B*h^4/(EI*h1)*(q1*h1+(q2-q1)*k*h)+4*c(n-1)-a(n-2)*c(n-1)-c(n-2))/d(n);
其中,E为桩体的弹性模量;I为桩体的横截面惯性矩;B为桩的实际宽度;B0为桩的计算宽度;K为地基反力系数;h为桩身长度;h1为加固段长度;q1为未加固段顶部荷载集度;q2为加固段底部荷载集度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011310962.9A CN112784332B (zh) | 2020-11-20 | 2020-11-20 | 基于K-m法的加固桩板墙桩前地基的桩身内力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011310962.9A CN112784332B (zh) | 2020-11-20 | 2020-11-20 | 基于K-m法的加固桩板墙桩前地基的桩身内力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112784332A true CN112784332A (zh) | 2021-05-11 |
CN112784332B CN112784332B (zh) | 2022-05-13 |
Family
ID=75750557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011310962.9A Active CN112784332B (zh) | 2020-11-20 | 2020-11-20 | 基于K-m法的加固桩板墙桩前地基的桩身内力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112784332B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378283A (zh) * | 2021-07-01 | 2021-09-10 | 中铁二院工程集团有限责任公司 | 锚固点增强型桩结构计算方法 |
CN117973165A (zh) * | 2024-04-02 | 2024-05-03 | 湖南大学 | 一种深厚软土区水平受荷桩桩体内力及位移计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015135045A (ja) * | 2013-12-18 | 2015-07-27 | 国立大学法人豊橋技術科学大学 | 長尺物の設置安定化法および長尺物の設置方法 |
CN108959752A (zh) * | 2018-06-26 | 2018-12-07 | 湘潭大学 | 一种适用于计算钢管三排桩桩身位移与桩身内力的计算方法 |
CN111046552A (zh) * | 2019-12-10 | 2020-04-21 | 中铁西北科学研究院有限公司 | 桩类抗滑结构内力计算的位移加载方法 |
US20200157763A1 (en) * | 2018-03-19 | 2020-05-21 | Dalian University Of Technology | Measuring device and method for horizontal dynamic impedance of specified foundation depth based on differential response analysis of pulse excitation |
-
2020
- 2020-11-20 CN CN202011310962.9A patent/CN112784332B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015135045A (ja) * | 2013-12-18 | 2015-07-27 | 国立大学法人豊橋技術科学大学 | 長尺物の設置安定化法および長尺物の設置方法 |
US20200157763A1 (en) * | 2018-03-19 | 2020-05-21 | Dalian University Of Technology | Measuring device and method for horizontal dynamic impedance of specified foundation depth based on differential response analysis of pulse excitation |
CN108959752A (zh) * | 2018-06-26 | 2018-12-07 | 湘潭大学 | 一种适用于计算钢管三排桩桩身位移与桩身内力的计算方法 |
CN111046552A (zh) * | 2019-12-10 | 2020-04-21 | 中铁西北科学研究院有限公司 | 桩类抗滑结构内力计算的位移加载方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378283A (zh) * | 2021-07-01 | 2021-09-10 | 中铁二院工程集团有限责任公司 | 锚固点增强型桩结构计算方法 |
CN117973165A (zh) * | 2024-04-02 | 2024-05-03 | 湖南大学 | 一种深厚软土区水平受荷桩桩体内力及位移计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112784332B (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112784332B (zh) | 基于K-m法的加固桩板墙桩前地基的桩身内力计算方法 | |
Osman et al. | Ground movement predictions for braced excavations in undrained clay | |
CN106777520A (zh) | 一种埋入式抗滑桩的设计方法 | |
CN109783924A (zh) | 地下水对地铁盾构隧道施工影响的数值模拟分析方法 | |
Ma et al. | Calculation of earth pressure on culvert underlying flexible subgrade | |
CN108763833B (zh) | 一种考虑土抗力突变的基坑支护桩挠度的计算方法 | |
CN106021963A (zh) | 一种用于基坑设计的计算方法 | |
CN113420482B (zh) | 一种基于结构内力监测值的管片荷载正交数值反演方法 | |
Feng et al. | Field monitoring of geogrid-reinforced and pile-supported embankment at bridge approach | |
Aicha et al. | Numerical modelling of the behavior of tunnel in soft surrounding rock: A case study of Djebel El-Ouahch Tunnel, Algeria | |
CN105117510B (zh) | 基于截面不均匀收缩效应的主梁挠曲变形预测方法 | |
Li et al. | Distributions of earth pressure and soil resistance on full buried single-row anti-sliding piles in loess slopes in northern Shaanxi based on in-situ model testing | |
Lian et al. | Theoretical analysis and optimization of frame protection to control shallow slope stability and soil erosion | |
Chen et al. | Failure analyses of a reinforced embankment by strength reduction and limit equilibrium methods considering hardening of soft clay | |
CN116992525A (zh) | 考虑自重、桩顶荷载和开挖作用的桩基竖向响应计算方法 | |
Ooi et al. | Numerical study of an integral abutment bridge supported on drilled shafts | |
Gorska et al. | An analysis of excavation support safety based on experimental studies | |
Chen et al. | Soil arching effect associated with ground movement and stress transfer adjacent to braced excavation in clayey ground | |
CN116484471A (zh) | 一种圆形地连墙墙体等效模型的等效参数确定方法 | |
CN115828605A (zh) | 一种顶端与沉陷道路刚性连接钢管桩竖向承载力计算方法 | |
Basile | A practical method for the non-linear analysis of piled rafts | |
CN111985021B (zh) | 一种盾构开挖面的遍布节理流固耦合的安全度分析方法 | |
Gorska et al. | Analysis of displacement of excavation based on inclinometer measurements | |
CN107476354A (zh) | 采用等代墩基模式预测超长群桩沉降量的计算方法 | |
CN113987645A (zh) | 一种山区斜坡直-斜组合桩基简化的内力计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |