CN110263483A - 基于flac3d软件的抗滑桩设计推力曲线计算方法 - Google Patents

基于flac3d软件的抗滑桩设计推力曲线计算方法 Download PDF

Info

Publication number
CN110263483A
CN110263483A CN201910588650.5A CN201910588650A CN110263483A CN 110263483 A CN110263483 A CN 110263483A CN 201910588650 A CN201910588650 A CN 201910588650A CN 110263483 A CN110263483 A CN 110263483A
Authority
CN
China
Prior art keywords
section
monitoring cross
thrust
soil body
soil
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
Application number
CN201910588650.5A
Other languages
English (en)
Other versions
CN110263483B (zh
Inventor
雷达
蒋关鲁
邹祖银
付正道
范凯
董青
戢林
贾少敏
冯研
高焱
李保险
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Agricultural University
Original Assignee
Sichuan Agricultural University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan Agricultural University filed Critical Sichuan Agricultural University
Priority to CN201910588650.5A priority Critical patent/CN110263483B/zh
Publication of CN110263483A publication Critical patent/CN110263483A/zh
Application granted granted Critical
Publication of CN110263483B publication Critical patent/CN110263483B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

本发明公开了一种基于FLAC3D软件的抗滑桩设计推力曲线计算方法,包括步骤:建立横向宽度1米的边坡三维模型;计算边坡滑动面的极限内摩擦角和极限粘聚力;划分边坡模型的竖向监测截面;计算监测截面后侧土体产生的弹性力;计算监测截面后侧土体产生的滑坡推力;计算监测截面前侧土体抗力;根据不同的工程设计安全系数,计算各个监测截面的抗滑桩设计推力,绘制边坡的抗滑桩设计推力曲线。本发明利用摩尔‑库仑塑性本构模型、弹性本构模型进行边坡的抗滑桩设计推力计算,能合理反映边坡滑坡推力的变化规律,提高了计算精度。

Description

基于FLAC3D软件的抗滑桩设计推力曲线计算方法
技术领域
本发明涉及岩土工程数值模拟计算技术领域,特别涉及一种基于FLAC3D软件的抗滑桩设计推力曲线计算方法。
背景技术
进行抗滑桩加固边坡的支挡结构设计时,抗滑桩的设计推力计算十分重要,抗滑桩设计推力的计算结果直接影响桩身截面尺寸、纵向配筋、锚固段长度、桩间距等设计指标,也会间接影响工程施工的难易程度、项目预算、施工周期等。基于工程设计目标,采用合理的计算理论,提高抗滑桩设计推力计算的合理性和精度,具有非常实用的工程应用价值和理论意义。
在我国滑坡灾害治理中,多采用传递系数法计算抗滑桩的设计推力,该方法采用刚体极限平衡理论,计算过程忽略边坡几何形态的影响,滑动面倾角变化能显著改变抗滑桩设计推力的计算结果。由于刚体极限平衡理论忽略土体的弹塑性本构特性,不能合理反映边坡土体的受力变形特性和破坏过程。针对传递系数法计算理论的局限性,有必要建立一种以弹塑性本构模型为基础的抗滑桩设计推力曲线计算方法,充分考虑土体的受力变形特性和边坡形态,弥补现有传递系数法计算理论的不足,具有重要的工程实践意义。
针对传递系数法计算边坡的抗滑桩设计推力曲线不能考虑土体的弹塑性本构特性、边坡形态,以及抗滑桩设计推力计算结果过于依赖滑动面倾角变化的特点,基于FLAC3D数值模拟计算软件,以土体的弹塑性本构模型为基础,建立横向宽度为1米的边坡三维模型,根据工程设计安全系数FS,按照公式φF=tan-1((tanφ)/FS)折减滑动面极限内摩擦角φ,按照公式cF=c/FS折减滑动面极限粘聚力c,分别计算抗滑桩桩后土体的滑坡推力、弹性力,以及抗滑桩桩前土体抗力,桩后土体滑坡推力减去弹性力和桩前土体抗力的计算结果,即为监测截面的抗滑桩设计推力,将边坡各个监测截面的抗滑桩设计推力值连成曲线,建立基于FLAC3D软件的抗滑桩设计推力曲线计算方法。
发明内容
本发明针对现有技术的缺陷,提供了一种基于FLAC3D软件的抗滑桩设计推力曲线计算方法,解决现有传递系数法忽略土体的弹塑性本构特性、不考虑边坡几何形态、滑动面倾角变化对滑坡推力影响过大的缺陷。
为了实现以上发明目的,本发明采取的技术方案如下:
一种基于FLAC3D软件的抗滑桩设计推力曲线计算方法,其特征在于,包括:横向宽度为1米的边坡三维模型、滑动面极限内摩擦角和极限粘聚力的计算方法、监测截面后侧土体产生弹性力的计算方法、监测截面后侧土体产生滑坡推力的计算方法、监测截面前侧土体抗力的计算方法和监测截面抗滑桩设计推力的计算方法;
根据所述横向宽度为1米的边坡三维模型,基于二分法的计算原理,利用LOOP循环语句折减滑动面内摩擦角和粘聚力的计算方法计算极限值,利用弹性本构模型计算监测截面后侧土体产生的弹性力,根据工程设计安全系数,利用摩尔-库仑本构模型计算滑动面极限内摩擦角和极限粘聚力折减时的桩后滑坡推力和桩前土体抗力,桩后滑坡推力减去弹性力和桩前土体抗力的计算结果,即为各个监测截面的抗滑桩设计推力,将它们连成曲线,得到边坡的抗滑桩设计推力曲线。
进一步地,所述横向宽度为1米的边坡三维模型具体为:边坡模型采用实体单元网格模拟,几何尺寸和原型保持一致,边坡模型的横向宽度均为1米,在边坡模型的前侧边界、后侧边界、底部边界和横向宽度方向上,设置固定速度和位移的边界条件。
进一步地,所述滑动面极限内摩擦角和极限粘聚力的计算方法,具体为:建立横向宽度为1米的边坡三维模型,以力不平衡比率不大于1×10-5为收敛指标,将土体粘聚力和抗拉强度设为1×1010Pa,采用更改强度参数的摩尔-库仑本构模型生成初始地应力场。设置折减系数的上限值k2和下限值k1,设定折减系数的计算精度,并预定最大计算时步数。根据二分法的计算原理,采用FISH语言编写LOOP循环语句,参照公式折减边坡滑动面的内摩擦角φ0和粘聚力c0,当力不平衡比率不大于1×10-5时,则下一次计算中的下限值替换为当力不平衡比率大于1×10-5时,则下一次计算中上限值替换为当经过数次折减计算,最终折减系数仍为k2时,则应设置更大的上限值;当经过数次折减计算,最终折减系数仍为k1时,则应设置更小的下限值。利用摩尔-库仑本构模型求解相应的地应力场,通过循环折减计算,当上下限值的差值小于计算精度时,计算结束,得到滑动面的极限内摩擦角和极限粘聚力。
进一步地,所述监测截面后侧土体产生弹性力的计算方法,具体为:监测截面沿坡脚方向的岩土体定义为刚体,利用弹性本构模型计算监测截面相邻土体单元网格的水平土压力,再乘以相应监测截面土体单元网格高度,求和汇总得到监测截面后侧土体产生的弹性力。
进一步地,所述监测截面后侧土体产生滑坡推力的计算方法,具体为:根据工程设计安全系数FS,按照公式φF=tan-1((tanφ)/FS)和cF=c/FS,对滑动面的极限内摩擦角φ和极限粘聚力c进行折减,将监测截面沿坡脚方向的岩土体定义为刚体,利用摩尔-库仑本构模型计算监测截面相邻土体单元网格的水平土压力,再乘以相应监测截面土体单元网格高度,求和汇总得到监测截面后侧土体产生的滑坡推力。
进一步地,所述监测截面前侧土体抗力的计算方法,具体为:根据工程设计安全系数FS,按照公式φF=tan-1((tanφ)/FS)和cF=c/FS,对滑动面的极限内摩擦角φ和极限粘聚力c进行折减,将监测截面沿坡顶方向的岩土体定义为刚体,利用摩尔-库仑本构模型计算监测截面相邻土体单元网格的水平土压力,再乘以相应监测截面土体单元网格高度,求和汇总得到监测截面前侧的土体抗力。
进一步地,所述监测截面抗滑桩设计推力的计算方法,具体为:抗滑桩设计推力=滑坡推力-弹性力-土体抗力。
与现有技术相比,本发明的优点在于:基于FLAC3D数值模拟计算软件,利用摩尔-库仑塑性本构模型、弹性本构模型进行边坡的抗滑桩设计推力计算,能合理反映边坡滑坡推力的变化规律,提高计算精度。
附图说明
图1为本发明基于FLAC3D软件的抗滑桩设计推力曲线计算流程图;
图2为本发明横向宽度1米的边坡三维模型图;
图3为本发明滑动面极限内摩擦角和极限粘聚力的计算流程图;
图4为本发明边坡模型的竖向监测截面划分示意图;
图5为本发明计算监测截面后侧土体弹性力或滑坡推力的三维模型图;
图6为本发明计算监测截面前侧土体抗力的三维模型图;
图7为本发明边坡的抗滑桩设计推力曲线图;
图7a为工程设计安全系数FS=1.1时的抗滑桩设计推力曲线图;
图7b为工程设计安全系数FS=1.25时的抗滑桩设计推力曲线图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下根据附图并列举实施例,对本发明做进一步详细说明。
一种基于FLAC3D软件的抗滑桩设计推力曲线计算方法,如图1所示,该方法包括:
步骤1:建立横向宽度1米的边坡三维模型
根据边坡原型几何尺寸,利用FLAC3D数值模拟软件建立实体单元网格的边坡三维模型,横向宽度均为1米,如图2所示,在边坡模型的前侧边界、后侧边界、底部边界、横向宽度方向均固定节点速度和位移。
步骤2:计算边坡滑动面的极限内摩擦角和极限粘聚力
如图3所示,建立横向宽度为1米的边坡三维模型,以力不平衡比率不大于1×10-5为收敛指标,将土体粘聚力和抗拉强度设为1×1010Pa,采用更改强度参数的摩尔-库仑本构模型生成初始地应力场。设置滑动面内摩擦角和粘聚力的折减系数上限值k2(如4)和下限值k1(如0.5),设定折减系数的计算精度为0.001,将计算时步的最大值设为10000步。根据二分法的计算原理,采用FISH语言编写LOOP循环语句,参照公式折减边坡滑动面的内摩擦角φ0和粘聚力c0,当力不平衡比率不大于1×10-5时,则下一次计算中的下限值替换为当力不平衡比率大于1×10-5时,则下一次计算中上限值替换为当经过数次折减计算,最终折减系数仍为k2时,则应设置更大的上限值;当经过数次折减计算,最终折减系数仍为k1时,则应设置更小的下限值。利用摩尔-库仑本构模型求解相应的地应力场,通过循环折减计算,当上下限值的差值小于计算精度0.001时,计算结束,得到滑动面的极限内摩擦角φ和极限粘聚力c。
步骤3:划分边坡模型的监测截面
根据工程设计要求,沿水平方向在边坡模型内部划分一定数量的竖向监测截面,如图4所示,监测截面即抗滑桩的考察桩位。
步骤4:计算监测截面后侧土体产生的弹性力
监测截面沿坡脚方向的岩土体定义为刚体,如图5所示,利用弹性本构模型进行计算,提取监测截面相邻土体单元网格的水平土压力b,再乘以相应监测截面土体单元网格高度h,监测截面后侧土体产生的弹性力为(b1×h1+…+bn×hn)。
步骤5:计算监测截面后侧土体产生的滑坡推力
监测截面沿坡脚方向的岩土体定义为刚体,如图5所示,以力不平衡比率不大于1×10-5为收敛指标,将土体粘聚力和抗拉强度设为1×1010Pa,利用更改强度参数的摩尔-库仑本构模型生成初始地应力场。根据工程设计安全系数FS,按照公式φF=tan-1((tanφ)/FS)和cF=c/FS,对滑动面的极限内摩擦角φ和极限粘聚力c进行折减,利用摩尔-库仑本构模型进行计算,提取监测截面相邻土体单元网格的水平土压力a,再乘以相应监测截面土体单元网格高度h,监测截面后侧土体产生的滑坡推力为(a1×h1+…+an×hn)。
步骤6:计算监测截面前侧土体抗力
监测截面沿坡顶方向的岩土体定义为刚体,如图6所示,以力不平衡比率不大于1×10-5为收敛指标,将土体粘聚力和抗拉强度设为1×1010Pa,利用更改强度参数的摩尔-库仑本构模型生成初始地应力场。根据工程设计安全系数FS,按照公式φF=tan-1((tanφ)/FS)和cF=c/FS,对滑动面的极限内摩擦角φ和极限粘聚力c进行折减,利用摩尔-库仑本构模型进行计算,当力不平衡比率大于1×10-5时,不计入抗滑桩前侧土体抗力,当力不平衡比率不大于1×10-5时,提取监测截面相邻土体单元网格的水平土压力c,再乘以相应监测截面土体单元网格高度h,监测截面前侧的土体抗力为(c1×h1+…+cn×hn)。
步骤7:根据不同的工程设计安全系数,计算监测截面的抗滑桩设计推力,对于每一个监测截面,抗滑桩设计推力=(a1×h1+…+an×hn)-(b1×h1+…+bn×hn)-(c1×h1+…+cn×hn),绘制边坡的抗滑桩设计推力曲线,如图7所示。
基于FLAC3D软件的抗滑桩设计推力曲线计算方法的计算精度和边坡三维模型单元体网格划分的精细程度相关,两者呈正相关性。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (7)

1.一种基于FLAC3D软件的抗滑桩设计推力曲线计算方法,其特征在于,包括:横向宽度为1米的边坡三维模型、滑动面极限内摩擦角和极限粘聚力的计算方法、监测截面后侧土体产生弹性力的计算方法、监测截面后侧土体产生滑坡推力的计算方法、监测截面前侧土体抗力的计算方法和监测截面抗滑桩设计推力的计算方法;
根据所述横向宽度为1米的边坡三维模型,基于二分法的计算原理,利用LOOP循环语句折减滑动面内摩擦角和粘聚力的计算方法计算极限值,利用弹性本构模型计算监测截面后侧土体产生的弹性力,根据工程设计安全系数,利用摩尔-库仑本构模型计算滑动面极限内摩擦角和极限粘聚力折减时的桩后滑坡推力和桩前土体抗力,桩后滑坡推力减去弹性力和桩前土体抗力的计算结果,即为各个监测截面的抗滑桩设计推力,将它们连成曲线,得到边坡的抗滑桩设计推力曲线。
2.根据权利要求1所述的方法,其特征在于:所述横向宽度为1米的边坡三维模型具体为:边坡模型采用实体单元网格模拟,几何尺寸和原型保持一致,边坡模型的横向宽度均为1米,在边坡模型的前侧边界、后侧边界、底部边界和横向宽度方向上,设置固定速度和位移的边界条件。
3.根据权利要求1所述的方法,其特征在于:所述滑动面极限内摩擦角和极限粘聚力的计算方法,具体为:建立横向宽度为1米的边坡三维模型,以力不平衡比率不大于1×10-5为收敛指标,将土体粘聚力和抗拉强度设为1×1010Pa,采用更改强度参数的摩尔-库仑本构模型生成初始地应力场;设置折减系数的上限值k2和下限值k1,设定折减系数的计算精度,并预定最大计算时步数;根据二分法的计算原理,采用FISH语言编写LOOP循环语句,参照公式折减边坡滑动面的内摩擦角φ0和粘聚力c0,当力不平衡比率不大于1×10-5时,则下一次计算中的下限值替换为当力不平衡比率大于1×10-5时,则下一次计算中上限值替换为当经过数次折减计算,最终折减系数仍为k2时,则应设置更大的上限值;当经过数次折减计算,最终折减系数仍为k1时,则应设置更小的下限值;利用摩尔-库仑本构模型求解相应的地应力场,通过循环折减计算,当上下限值的差值小于计算精度时,计算结束,得到滑动面的极限内摩擦角和极限粘聚力。
4.根据权利要求1所述的方法,其特征在于:所述监测截面后侧土体产生弹性力的计算方法,具体为:监测截面沿坡脚方向的岩土体定义为刚体,利用弹性本构模型计算监测截面相邻土体单元网格的水平土压力,再乘以相应监测截面土体单元网格高度,求和汇总得到监测截面后侧土体产生的弹性力。
5.根据权利要求1所述的方法,其特征在于:所述监测截面后侧土体产生滑坡推力的计算方法,具体为:根据工程设计安全系数FS,按照公式φF=tan-1((tanφ)/FS)和cF=c/FS,对滑动面的极限内摩擦角φ和极限粘聚力c进行折减,将监测截面沿坡脚方向的岩土体定义为刚体,利用摩尔-库仑本构模型计算监测截面相邻土体单元网格的水平土压力,再乘以相应监测截面土体单元网格高度,求和汇总得到监测截面后侧土体产生的滑坡推力。
6.根据权利要求1所述的方法,其特征在于:所述监测截面前侧土体抗力的计算方法,具体为:根据工程设计安全系数FS,按照公式φF=tan-1((tanφ)/FS)和cF=c/FS,对滑动面的极限内摩擦角φ和极限粘聚力c进行折减,将监测截面沿坡顶方向的岩土体定义为刚体,利用摩尔-库仑本构模型计算监测截面相邻土体单元网格的水平土压力,再乘以相应监测截面土体单元网格高度,求和汇总得到监测截面前侧的土体抗力。
7.根据权利要求1所述的方法,其特征在于:所述监测截面抗滑桩设计推力的计算方法,具体为:抗滑桩设计推力=滑坡推力-弹性力-土体抗力。
CN201910588650.5A 2019-07-02 2019-07-02 基于flac3d软件的抗滑桩设计推力曲线计算方法 Expired - Fee Related CN110263483B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910588650.5A CN110263483B (zh) 2019-07-02 2019-07-02 基于flac3d软件的抗滑桩设计推力曲线计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910588650.5A CN110263483B (zh) 2019-07-02 2019-07-02 基于flac3d软件的抗滑桩设计推力曲线计算方法

Publications (2)

Publication Number Publication Date
CN110263483A true CN110263483A (zh) 2019-09-20
CN110263483B CN110263483B (zh) 2022-12-13

Family

ID=67923790

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910588650.5A Expired - Fee Related CN110263483B (zh) 2019-07-02 2019-07-02 基于flac3d软件的抗滑桩设计推力曲线计算方法

Country Status (1)

Country Link
CN (1) CN110263483B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110851898A (zh) * 2019-10-24 2020-02-28 中铁二院工程集团有限责任公司 一种基于桥梁墩台变形控制的滑坡抗滑桩设计方法
CN111553013A (zh) * 2020-05-06 2020-08-18 中铁二院工程集团有限责任公司 一种模拟抗滑桩对桥梁墩台作用的简化数值模拟方法
CN112924649A (zh) * 2021-03-17 2021-06-08 四川农业大学 一种滑坡模型水平加速度放大系数测试方法
CN115168953A (zh) * 2022-07-11 2022-10-11 哈尔滨工业大学 一种基于边坡稳定性的山区公路线位优化方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105956281A (zh) * 2016-05-05 2016-09-21 中国人民解放军国防科学技术大学 固体火箭发动机装药设计方法
CN107066739A (zh) * 2017-04-17 2017-08-18 安徽理工大学 一种基于flac3d数值模拟的断层注浆加固效果评价方法
CN107330182A (zh) * 2017-06-28 2017-11-07 西北农林科技大学 基于增湿条件的强度折减计算边坡稳定安全系数的方法
CN108717503A (zh) * 2018-05-25 2018-10-30 中建隧道建设有限公司 一种设计工况下埋入式抗滑桩承担的滑坡推力计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105956281A (zh) * 2016-05-05 2016-09-21 中国人民解放军国防科学技术大学 固体火箭发动机装药设计方法
CN107066739A (zh) * 2017-04-17 2017-08-18 安徽理工大学 一种基于flac3d数值模拟的断层注浆加固效果评价方法
CN107330182A (zh) * 2017-06-28 2017-11-07 西北农林科技大学 基于增湿条件的强度折减计算边坡稳定安全系数的方法
CN108717503A (zh) * 2018-05-25 2018-10-30 中建隧道建设有限公司 一种设计工况下埋入式抗滑桩承担的滑坡推力计算方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
万玲等: "抗滑桩抗滑阻力的优化计算方法", 《重庆大学学报》 *
曹晓川等: "强度折减法在抗滑桩加固土坡稳定分析中的应用", 《山西建筑》 *
李起龙等: "基于多目标综合评价法的边坡抗滑桩桩位优化设计", 《铁道科学与工程学报》 *
程昌明等: "降雨入渗对砂质边坡变形和稳定性的影响", 《金属矿山》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110851898A (zh) * 2019-10-24 2020-02-28 中铁二院工程集团有限责任公司 一种基于桥梁墩台变形控制的滑坡抗滑桩设计方法
CN111553013A (zh) * 2020-05-06 2020-08-18 中铁二院工程集团有限责任公司 一种模拟抗滑桩对桥梁墩台作用的简化数值模拟方法
CN111553013B (zh) * 2020-05-06 2022-06-14 中铁二院工程集团有限责任公司 一种模拟抗滑桩对桥梁墩台作用的简化数值模拟方法
CN112924649A (zh) * 2021-03-17 2021-06-08 四川农业大学 一种滑坡模型水平加速度放大系数测试方法
CN115168953A (zh) * 2022-07-11 2022-10-11 哈尔滨工业大学 一种基于边坡稳定性的山区公路线位优化方法

Also Published As

Publication number Publication date
CN110263483B (zh) 2022-12-13

Similar Documents

Publication Publication Date Title
CN110263483A (zh) 基于flac3d软件的抗滑桩设计推力曲线计算方法
Tang et al. A novel approach for determining landslide pushing force based on landslide-pile interactions
CN107330224B (zh) 一种分条间作用力倾角非假定的边坡稳定性分析条分法
Xie et al. Lateral load bearing capacity of offshore high-piled wharf with batter piles
Andrianopoulos et al. Explicit integration of bounding surface model for the analysis of earthquake soil liquefaction
Jia et al. Three-dimensional modeling of bank erosion and morphological changes in the Shishou bend of the middle Yangtze River
US7966165B2 (en) Soil-water coupled analyzer and soil-water coupled analysis method
Li et al. A novel optimal plane arrangement of stabilizing piles based on soil arching effect and stability limit for 3D colluvial landslides
CN104899380A (zh) 一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法
Galavi Groundwater flow, fully coupled flow deformation and undrained analyses in PLAXIS 2D and 3D
CN111046552A (zh) 桩类抗滑结构内力计算的位移加载方法
CN107563064A (zh) 一种海啸波越浪过程的二维数值模拟方法
Memon A comparison between limit equilibrium and finite element methods for slope stability analysis
Yazd et al. Optimal designing of concrete gravity dam using particle swarm optimization algorithm (PSO)
Dong et al. Limit analysis of vertical anti-pulling screw pile group under inclined loading on 3D elastic-plastic finite element strength reduction method
CN115983068A (zh) 一种非饱和土边坡稳定性测量方法、装置及存储介质
Mirsaidov et al. Assessment of the stress-strain state of earth dams
CN106759083B (zh) 适用于深水卸荷式板桩岸壁结构的概化模型设计方法
CN115169191A (zh) 一种斜直组合支护结构受力计算方法
Farzaneh et al. Active earth pressure induced by strip loads on a backfill
Bürg et al. Application of a second-order implicit material point method
Montesi Soil Structure Interaction Using p–y Analyses for Earth Retaining Structures
JP7103046B2 (ja) 多次元圧密理論に基づく地盤の解析方法
Ai et al. The Analysis of Earth Pressure on Retaining Wall Based on Abaqus
Dysarz Basics of river flow modelling

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20221213