CN112100863A - 用于桩基工程的桩基负摩擦力解析解分析方法 - Google Patents
用于桩基工程的桩基负摩擦力解析解分析方法 Download PDFInfo
- Publication number
- CN112100863A CN112100863A CN202011054936.4A CN202011054936A CN112100863A CN 112100863 A CN112100863 A CN 112100863A CN 202011054936 A CN202011054936 A CN 202011054936A CN 112100863 A CN112100863 A CN 112100863A
- Authority
- CN
- China
- Prior art keywords
- pile
- stress
- soil
- load
- length
- 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
Images
Landscapes
- Piles And Underground Anchors (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种用于桩基工程的桩基负摩擦力解析解分析方法,该分析方法包括如下步骤:首先利用传统方法计算群桩势能,引入传统理论荷载传递曲线,并对其进行变分,本发明采用了在桩底为零而沿桩长线性变化的规律,这与实际比较吻合,桩-土滑落超过应力极限之后会产生桩-土滑落,本发明中土为弹塑性,采用初应力法。首先求出弹性解;其次,把超过应力极限的过量应力转换成附加载荷进行新的循环,使得过量应力重新分布到桩上没有超过应力极限的区域。该迭代反复执行,直到所有区域的应力都没有超过应力极限。
Description
技术领域
本发明属于岩土工程技术领域,具体涉及一种用于桩基工程的桩基负摩擦力解析解计算方法。
背景技术
桩基工程对于沿海地区至关重要。收集已发表的科研成果,发现国内对桩基工程的理论研究还处于起步阶段,现有技术关于桩土之间的负摩擦力机理尚不清楚,同时现有的通用有限元计算软件无法反映土的软化,对于桩的负摩阻力随桩深变化规律还没有理论解析方法。
发明内容
为了解决现有技术问题,本发明的目的在于克服已有技术存在的不足,提供一种用于桩基工程的桩基负摩擦力解析解分析方法,能计算用于桩基工程的桩基负摩擦力解析解,以解决了背景技术中提出的问题,作为桩土之间的负摩擦力机理的重要参考。
为达到上述发明创造目的,本发明构思如下:
本发明利用传统方法计算群桩势能,引入传统理论荷载传递曲线,并对其进行变分,采用了在桩底为零而沿桩长线性变化的规律,这与实际比较吻合,桩-土滑落超过应力极限之后会产生桩-土滑落。土为弹塑性,采用初应力法。首先求出弹性解;其次,把超过应力极限的过量应力转换成附加载荷进行新的循环,使得过量应力重新分布到桩上没有超过应力极限的区域。该迭代反复执行,直到所有区域的应力都没有超过应力极限。
根据上述发明构思,本发明采用如下技术方案:
一种用于桩基工程的桩基负摩擦力解析解分析方法,所述计算方法包括如下步骤:
A.确定装群势能:
上述式中第一项为应变能,后三项为外力所作的功,包括沿桩长方向即桩轴线方向的应力{τz}和桩底的应力{σl}以及桩头荷载{pt},{wz}、{wl}和{wt}为它们对应的位移;式中:
{pt}为桩头荷载,
{τz}为沿桩长方向即桩轴线方向的应力,
{σl}为桩底的应力,
{Gl}为土壤沿桩长方向即桩轴线方向的线性变化的刚度,到达桩底时等于{Gb},
{Gb}为土壤在桩底的刚度;
B.构建土壤的载荷-位移关系:
理论载荷传递曲线为:
{τz}=[k]{ωz}
{σl}=[kb]{ωl}
[k]和[kb]为土壤对于基桩和桩底的刚度阵,其中系数随深度线性变化;刚度阵中的距离桩轴线的平均半径rm、土壤运动产生切应力,可由下式计算:
rm=αρl(1-υs)
υs为土壤泊松比,l为桩长,ρ为土壤不均匀性因子而α为:
C.采用基本变分公式:
若在深度z的土壤自由沉降量为{sz},而在桩底为{sl},则响应的纯粹土壤运动
{ω}={sz}+{ωz}
{ωb}={sl}+{ωl}
即{wz}={w}-{sz}和{wt}={wb}-{sl},整理得:
于是
因为桩-土之间位移相容,所以在它们界面处位移相同,用有限项级数表示
其中βij为待定系数;由最小势能原理
导出:
上式写成矩阵形式:
[h]{β}={p}+{Qs}
向量{Qs}为:
{Qs}表示由于土壤运动而作用在桩基上的力,一旦土壤运动规律给定就可确定;采用在桩底为零而沿桩长线性变化的规律,这与实际比较吻合,均质土壤由于地面超载而产生的随深度线性变化的运动;
D.求解桩头响应解:
{β}=[f]{p}+{δs}
其中,{f}=[h]-1,{δs}=[h-1]{Qs};因为{β},{f}和{δs}的行和列根据每根桩的情况进行扩充,上述方程转换成桩头沉降{wt}和桩头荷载{pt}的表达式:
{wt}=[ft]{pt}+{δt}
其中,[ft]和{δt}为[f]和{δs}经过扩充得到的矩阵和向量;对[ft]进行求逆:
{Pt}=[kt]{wt}-{Qt}
其中,{kt}=[ft]-1,{Qt}=[ft]-1{δt};用桩顶刚度阵组装,从而给出桩顶-桩群系统的完全解答;考虑刚性桩头从而没有桩头垂直载荷的特殊情况,桩头荷载之和为零,统一的桩头沉降解:
其中,{1}为单位向量;一旦桩头沉降量解出,桩头载荷和待定系数用相应公式确定;沿桩长任意点的位移和切应力用以上所述相应公式计算;桩上的负向摩擦力利用切应力沿桩长积分求出;
E.进行桩-土滑落分析
超过应力极限之后会产生桩-土滑落,土为弹塑性,采用初应力法;首先求出弹性解;其次,把超过应力极限的过量应力转换成附加载荷进行新的循环,使得过量应力重新分布到桩上没有超过应力极限的区域;该迭代反复执行,直到所有区域的应力都没有超过应力极限。
本发明与现有技术相比较,具有如下显而易见的突出实质性特点和显著优点:
本发明采用用于桩基工程的桩基负摩擦力解析分析方法,首先利用传统方法计算群桩势能,引入传统理论荷载传递曲线,并对其进行变分,本发明采用了在桩底为零而沿桩长线性变化的规律,这与实际比较吻合,桩-土滑落超过应力极限之后会产生桩-土滑落,本发明中土为弹塑性,采用初应力法。首先求出弹性解;其次,把超过应力极限的过量应力转换成附加载荷进行新的循环,使得过量应力重新分布到桩上没有超过应力极限的区域。该迭代反复执行,直到所有区域的应力都没有超过应力极限。
附图说明
图1为本发明桩上应力示意图。其中图1(a)为桩上应力图;图1(b)为土壤剪切模量剖面示意图。
具体实施方式
以下结合具体的实施例子对上述方案做进一步说明,本发明的优选实施例详述如下:
在本实施例中,如图1所示,一种用于桩基工程的桩基负摩擦力解析解分析方法,所述计算方法包括如下步骤:
A.确定装群势能:
式中:
{pt}为桩头荷载,
{τz}为沿桩长方向(桩轴线方向)的应力,
{σl}为桩底的应力,
{Gl}为土壤沿桩长方向(桩轴线方向)线性变化的刚度,到达桩底时等于{Gb},
{Gb}为土壤在桩底的刚度;
式中,第一项为应变能,后三项为外力所作的功,包括沿桩长方向即桩轴线方向的应力{τz}和桩底的应力{σl}以及桩头荷载{pt},{wz}、{wl}和{wt}为它们对应的位移;
B.构建土壤的载荷-位移关系:
理论载荷传递曲线为:
{τz}=[k]{ωz}
{σl}=[kb]{ωl}
[k]和[kb]为土壤对于基桩和桩底的刚度阵,其中系数随深度线性变化;刚度阵中的距离桩轴线的平均半径rm、土壤运动产生切应力,可由下式计算:
rm=αρl(1-υs)
υs为土壤泊松比,l为桩长,ρ为土壤不均匀性因子而α为:
C.采用基本变分公式:
若在深度z的土壤自由沉降量为{sz},而在桩底为{sl},则响应的纯粹土壤运动
{ω}={sz}+{ωz}
{ωb}={sl}+{ωl}
即{wz}={w}-{sz}和{wt}={wb}-{sl},整理得:
于是
因为桩-土之间位移相容,所以在它们界面处位移相同,用有限项级数表示
其中βij为待定系数;由最小势能原理
导出:
上式写成矩阵形式:
[h]{β}={p}+{Qs}
向量{Qs}为:
{Qs}表示由于土壤运动而作用在桩基上的力,一旦土壤运动规律给定就可确定;采用在桩底为零而沿桩长线性变化的规律,这与实际比较吻合,均质土壤由于地面超载而产生的随深度线性变化的运动;
D.求解桩头响应解:
{β}=[f]{p}+{δs}
其中,{f}=[h]-1,{δs}=[h-1]{Qs};因为{β},{f}和{δs}的行和列根据每根桩的情况进行扩充,上述方程转换成桩头沉降{wt}和桩头荷载{pt}的表达式:
{wt}=[ft]{pt}+{δt}
其中,[ft]和{δt}为[f]和{δs}经过扩充得到的矩阵和向量;对[ft]进行求逆:
{Pt}=[kt]{wt}-{Qt}
其中,{kt}=[ft]-1,{Qt}=[ft]-1{δt};用桩顶刚度阵组装,从而给出桩顶-桩群系统的完全解答;考虑刚性桩头从而没有桩头垂直载荷的特殊情况,桩头荷载之和为零,统一的桩头沉降解:
其中,{1}为单位向量;一旦桩头沉降量解出,桩头载荷和待定系数用相应公式确定;沿桩长任意点的位移和切应力用以上所述相应公式计算;桩上的负向摩擦力利用切应力沿桩长积分求出;
E.进行桩-土滑落分析
超过应力极限之后会产生桩-土滑落,土为弹塑性,采用初应力法;首先求出弹性解;其次,把超过应力极限的过量应力转换成附加载荷进行新的循环,使得过量应力重新分布到桩上没有超过应力极限的区域;该迭代反复执行,直到所有区域的应力都没有超过应力极限。
综上所述,本发明用于桩基工程的桩基负摩擦力解析解分析方法,该分析方法包括如下步骤:首先利用传统方法计算群桩势能,引入传统理论荷载传递曲线,并对其进行变分,本发明采用了在桩底为零而沿桩长线性变化的规律,这与实际比较吻合,桩-土滑落超过应力极限之后会产生桩-土滑落,本发明中土为弹塑性,采用初应力法。首先求出弹性解;其次,把超过应力极限的过量应力转换成附加载荷进行新的循环,使得过量应力重新分布到桩上没有超过应力极限的区域。该迭代反复执行,直到所有区域的应力都没有超过应力极限。
上面对本发明实施例结合附图进行了说明,但本发明不限于上述实施例,还可以根据本发明的发明创造的目的做出多种变化,凡依据本发明技术方案的精神实质和原理下做的改变、修饰、替代、组合或简化,均应为等效的置换方式,只要符合本发明的发明目的,只要不背离本发明用于桩基工程的桩基负摩擦力解析解分析方法的技术原理和发明构思,都属于本发明的保护范围。
Claims (1)
1.一种用于桩基工程的桩基负摩擦力解析解分析方法,其特征在于:所述计算方法包括如下步骤:
A.确定装群势能:
式中:
{pt}为桩头荷载,
{τz}为沿桩长方向(桩轴线方向)的应力,
{σl}为桩底的应力,
{Gl}为土壤沿桩长方向即桩轴线方向的线性变化的刚度,到达桩底时等于{Gb},
{Gb}为土壤在桩底的刚度;
第一项为应变能,后三项为外力所作的功,包括沿桩长方向即桩轴线方向的应力{τz}和桩底的应力{σl}以及桩头荷载{pt},{wz}、{wl}和{wt}为它们对应的位移;
B.构建土壤的载荷-位移关系:
理论载荷传递曲线为:
{τz}=[k]{ωz}
{σl}=[kb]{ωl}
[k]和[kb]为土壤对于基桩和桩底的刚度阵,其中系数随深度线性变化;刚度阵中的距离桩轴线的平均半径rm、土壤运动产生切应力,可由下式计算:
rm=αρl(1-υs)
υs为土壤泊松比,l为桩长,ρ为土壤不均匀性因子而α为:
C.采用基本变分公式:
若在深度z的土壤自由沉降量为{sz},而在桩底为{sl},则响应的纯粹土壤运动
{ω}={sz}+{ωz}
{ωb}={sl}+{ωl}
即{wz}={w}-{sz}和{wt}={wb}-{sl},整理得:
于是
因为桩-土之间位移相容,所以在它们界面处位移相同,用有限项级数表示
其中βij为待定系数;由最小势能原理
导出:
上式写成矩阵形式:
[h]{β}={p}+{Qs}
向量{Qs}为:
{Qs}表示由于土壤运动而作用在桩基上的力,一旦土壤运动规律给定就可确定;采用在桩底为零而沿桩长线性变化的规律,这与实际比较吻合,均质土壤由于地面超载而产生的随深度线性变化的运动;
D.求解桩头响应解:
{β}=[f]{p}+{δs}
其中,{f}=[h]-1,{δs}=[h-1]{Qs};因为{β},{f}和{δs}的行和列根据每根桩的情况进行扩充,上述方程转换成桩头沉降{wt}和桩头荷载{pt}的表达式:
{wt}=[ft]{pt}+{δt}
其中,[ft]和{δt}为[f]和{δs}经过扩充得到的矩阵和向量;对[ft]进行求逆:
{Pt}=[kt]{wt}-{Qt}
其中,{kt}=[ft]-1,{Qt}=[ft]-1{δt};用桩顶刚度阵组装,从而给出桩顶-桩群系统的完全解答;考虑刚性桩头从而没有桩头垂直载荷的特殊情况,桩头荷载之和为零,统一的桩头沉降解:
其中,{1}为单位向量;一旦桩头沉降量解出,桩头载荷和待定系数用相应公式确定;沿桩长任意点的位移和切应力用以上所述相应公式计算;桩上的负向摩擦力利用切应力沿桩长积分求出;
E.进行桩-土滑落分析
超过应力极限之后会产生桩-土滑落,土为弹塑性,采用初应力法;首先求出弹性解;其次,把超过应力极限的过量应力转换成附加载荷进行新的循环,使得过量应力重新分布到桩上没有超过应力极限的区域;该迭代反复执行,直到所有区域的应力都没有超过应力极限。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010180938 | 2020-03-16 | ||
CN2020101809381 | 2020-03-16 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112100863A true CN112100863A (zh) | 2020-12-18 |
CN112100863B CN112100863B (zh) | 2022-06-14 |
Family
ID=73783790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011054936.4A Active CN112100863B (zh) | 2020-03-16 | 2020-09-29 | 用于桩基工程的桩基负摩擦力解析解分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100863B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113158291A (zh) * | 2021-03-01 | 2021-07-23 | 中国铁路设计集团有限公司 | 一种桥墩在地震作用下摩擦桩桩长的计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102808429A (zh) * | 2012-08-31 | 2012-12-05 | 机械工业勘察设计研究院 | 基于徐变应变分离的桩基浸水载荷试验方法 |
CN105839678A (zh) * | 2016-04-11 | 2016-08-10 | 湖北工业大学 | 一种改进型桩基竖向承载力试验反向自平衡法及试验装置 |
CN106650131A (zh) * | 2016-12-28 | 2017-05-10 | 中国能源建设集团江苏省电力设计院有限公司 | 组合受荷桩承载力传递矩阵通解方法 |
CN107577852A (zh) * | 2017-08-11 | 2018-01-12 | 上海岩土工程勘察设计研究院有限公司 | 用于超长桩荷载传递机理研究的荷载传递函数计算方法 |
-
2020
- 2020-09-29 CN CN202011054936.4A patent/CN112100863B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102808429A (zh) * | 2012-08-31 | 2012-12-05 | 机械工业勘察设计研究院 | 基于徐变应变分离的桩基浸水载荷试验方法 |
CN105839678A (zh) * | 2016-04-11 | 2016-08-10 | 湖北工业大学 | 一种改进型桩基竖向承载力试验反向自平衡法及试验装置 |
CN106650131A (zh) * | 2016-12-28 | 2017-05-10 | 中国能源建设集团江苏省电力设计院有限公司 | 组合受荷桩承载力传递矩阵通解方法 |
CN107577852A (zh) * | 2017-08-11 | 2018-01-12 | 上海岩土工程勘察设计研究院有限公司 | 用于超长桩荷载传递机理研究的荷载传递函数计算方法 |
Non-Patent Citations (3)
Title |
---|
W.Y.SHEN 等: "群桩-承台基础共同作用的变分分析法", 《地下空间》 * |
喻君: "改进的荷载传递法在桩基沉降计算中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
李林 等: "考虑沉桩效应的群桩非线性荷载-沉降解析", 《岩土力学》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113158291A (zh) * | 2021-03-01 | 2021-07-23 | 中国铁路设计集团有限公司 | 一种桥墩在地震作用下摩擦桩桩长的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112100863B (zh) | 2022-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Borja et al. | Elastoplastic consolidation at finite strain part 2: finite element implementation and numerical examples | |
Venanzi et al. | The effect of soil–foundation–structure interaction on the wind-induced response of tall buildings | |
CN112100863B (zh) | 用于桩基工程的桩基负摩擦力解析解分析方法 | |
CN103161234A (zh) | 工程结构多级设防烈度下的性能化抗震设计方法 | |
CN102787677B (zh) | 工程结构抗震烈度预测方法 | |
CN111460547A (zh) | 考虑加筋效应的扩底抗拔群桩的变形分析方法 | |
Zhang et al. | Seismic behaviour of high-rise frame-core tube structures considering dynamic soil–structure interaction | |
Tso et al. | Dynamic analysis of plane coupled shear walls | |
CN115855686A (zh) | 一种测算水泥搅拌桩复合地基抗剪强度的方法及装置 | |
CN108897966B (zh) | 基于弹塑性修正的屈曲约束支撑结构的等效弹性分析方法 | |
CN104792560A (zh) | 条石砌筑石墙灰缝力学计算模型、试验装置及试验方法 | |
CN116167229A (zh) | 一种高效的海上风电桩周土抗力模型算法及装置 | |
Fernández-Escobar et al. | Lateral deformation and buckling analysis of piles including shear effects: Numerical analysis | |
CN115859418A (zh) | 一种海上风电大直径单桩水平动力响应分析方法 | |
Zhang et al. | The application of the fem in the prediction of seismic subsidence of loess ground | |
CN110990910A (zh) | 时程激励下线性耗能结构响应的快速迭代法 | |
CN204903188U (zh) | 条石砌筑石墙灰缝力学计算模型及试验装置 | |
Baltaji et al. | Non-linear time domain site response and soil structure interaction analyses for nuclear facilities using moose | |
CN117312796A (zh) | 非等长桩体复合地基竖向承载特性的积分方程通解方法 | |
Shukla et al. | Time-dependent settlement analysis of a geosynthetic-reinforced soil | |
Deka | Different analysis methods of piled rafts | |
CN116720372A (zh) | 一种桩基竖向承载性状分析方法及系统 | |
Leti et al. | Use of Incremental Dynamic Analysis for the seismic performance assessment of RC framed structures | |
Chehat et al. | Modelling of non-linear seismic ground response using elasto-plastic constitutive framework within a finite element soil column model | |
Ghorbani et al. | Estimation of seismic demand in tall buildings using mode-acceleration method |
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 |