CN114676610A - 一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法 - Google Patents
一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法 Download PDFInfo
- Publication number
- CN114676610A CN114676610A CN202210278842.8A CN202210278842A CN114676610A CN 114676610 A CN114676610 A CN 114676610A CN 202210278842 A CN202210278842 A CN 202210278842A CN 114676610 A CN114676610 A CN 114676610A
- Authority
- CN
- China
- Prior art keywords
- soil body
- underwater structure
- analysis
- model
- area
- 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.)
- Pending
Links
Images
Classifications
-
- 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
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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)
- Piles And Underground Anchors (AREA)
Abstract
本发明公开了一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,本发明考虑了水下结构开挖过程,将土体划分为桩靴下压导致的柱孔区、水下结构的开挖区、影响区三个区域进行分析,依据柱孔扩张理论,将桩靴贯入过程对周围土体的作用力简化为作用在桩靴下压深度处的基底压力和作用在土体侧壁的扩孔张力,进而准确计算实际工程中插桩作业对邻近结构物变形及应力的影响;本发明运用ABAQUS有限元软件,采用静力分析步,计算效率高,耗时短,计算水下结构变形和应力数值可靠,计算结果与离心机试验结果吻合较好,可以有效解决实际工程中的问题。本发明符合工程实际,方法简单明确,易于计算。
Description
技术领域
本发明属于数值模拟技术领域,具体涉及一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法。
背景技术
经济的快速发展往往伴随着能源的大量消耗,而陆上油气资源的开发程度很高,已探明的可开发油气资源越来越少。相比之下,海上油气资源较为丰富,因此大规模开发海上油气资源已成必然趋势。自升式钻井平台属于海洋移动平台,具有定位能力强、作业稳定性好的优势,在海洋油气资源勘探开发领域中的应用越来越广泛。自升式钻井平台的桩靴直径一般在10~20m之间,钻井作业时桩靴的插拔会排开大量土体,土体的移动势必会对周边平台基础或水下防护结构等产生不利影响,危害设施安全,如何评估自升式钻井平台插拔桩作业对邻近结构物的影响成为近年来的研究重点。
目前,常用的研究桩靴插拔对邻近结构物影响的数值分析方法主要是基于ABAQUS数值分析软件的耦合欧拉-拉格朗日(CEL)方法。CEL方法主要用来做流固耦合仿真或者大变形仿真,其在解决桩靴插拔对邻近结构物的影响问题时,存在三个缺陷:①无法准确处理插桩对周边开挖坑槽等孔穴结构的影响。现有CEL方法中一般在欧拉土体上方10m范围内设置空穴,为桩靴插拔过程中土体向上隆起变形提供空间,插桩作业时土体不断向上方空穴流动。但若土体存在开挖结构,其在开挖处需将欧拉土体也设置为空穴,此时土体流动机制改变,甚至会出现时间效应,影响邻近结构物的应力和变形。②由于CEL方法原理是计算过程中空间网格的形状、大小和位置保持不变,土体在网格间不断流动,因此使用CEL大变形算法计算得到的承载力曲线不稳定,得到的邻近结构物变形结果偏大。③CEL算法计算效率低,耗时较长。
针对以上问题,以柱孔扩张理论为依据,基于ABAQUS中的静力分析,提出了一种分析插桩作业对邻近结构物影响的有限元分析方法,此方法可考虑存在开挖坑槽等孔穴存在时桩靴插拔对周围土体及结构物的影响,且采用静力分析方法,计算耗时短,有效解决了计算效率低的问题。
发明内容
本发明的目的旨在提供一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,该方法依据柱孔扩张理论,将桩靴贯入过程对周围土体的作用力简化为作用在桩靴下压深度处的基底压力和作用在土体侧壁的扩孔张力,并使用ABAQUS中的单元生死功能,实现静力法对桩靴连续贯入过程的模拟,进而准确评估实际工程中水下结构受桩靴插拔作用的影响。
本发明所采用的技术方案:
一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,其具体计算步骤为:
S1、根据分析工况,建立对应的有限元分析模型,具体包括:土体模型和水下结构模型,土体模型划分为三个区域,包括:桩靴下压导致的柱孔区、水下结构的开挖区和影响区,桩靴直径即柱孔区的直径为D1,桩靴贯入深度即柱孔区的深度为d,开挖区直径为D2,开挖区的深度为h,影响区的宽度为≥10D1,影响区的深度为≥10d;水下结构模型的外径为D2,水下结构模型的高度为H,水下结构模型的壁厚为t。
S2、建立土体材料属性和水下结构材料属性,并将其分别赋予土体模型与水下结构模型。
S3、根据计算工况,将水下结构模型与土体模型装配至一起,形成装配体,需设置的参数包括:桩靴与水下结构模型之间的距离和水下结构模型的埋置深度。
S4、建立数值分析模型的计算分析步。共建5个分析步,包括:地应力分析步、水下结构安放分析步、水下结构开挖分析步、柱孔区形成分析步和柱孔区扩张分析步。
S5、建立相互作用。首先建立水下结构与土体的接触属性,一般采用面面接触。然后使用相互作用中的单元生死功能,模拟水下结构的安放与内部土体开挖过程及桩靴下压产生柱孔区的过程。具体步骤为:使用单元生死功能将水下结构在地应力分析步设置为无效,并在水下结构安放分析步设置为重激活;将水下结构开挖区在开挖分析步设置为无效,以模拟该区域土体的开挖;将柱孔区在柱孔形成分析步中设置为无效,以模拟该区域土体由于桩靴下压形成柱孔区。
S6、建立模型的边界条件。对土体模型设置合理的位移边界,包括:土体底部全约束、土体侧面放开沿深度方向的位移,其余全约束、土体顶面约束沿深度方向的位移。
S7、建立荷载条件。包括:①在地应力分析步中对全部土体施加重力。②在水下结构安放分析步中对水下结构施加重力。③在柱孔区扩张分析步中对区域内底部土体表面施加等同于基底压力大小的压强。④在柱孔区扩张分析步中对区域内侧向土体表面施加等同于扩孔力大小的压强。
S8、对土体施加地应力场。
S9、对土体模型及水下结构模型进行网格划分。
S10、提交作业进行计算。
本发明具有的优点和积极效果是:
1、考虑了水下结构开挖过程,将土体划分为桩靴下压导致的柱孔区、水下结构的开挖区、影响区三个区域进行分析,依据柱孔扩张理论,将桩靴贯入过程对周围土体的作用力简化为作用在桩靴下压深度处的基底压力和作用在土体侧壁的扩孔张力,进而准确计算实际工程中插桩作业对邻近结构物变形及应力的影响。
2、运用ABAQUS有限元软件,采用静力分析步,计算效率高,耗时短,计算水下结构变形和应力数值可靠,计算结果与离心机试验结果吻合较好,可以有效解决实际工程中的问题。
3、综上所述,本发明符合工程实际,方法简单明确,易于计算。
附图说明
图1是土体模型区域划分示意图。
图2是柱孔扩张示意图。
图3是环肋相对位置示意图。
图4是有限元网格划分示意图。
图5是离心试验钢圆筒位移及应变片测点布置图。
图6是钢圆筒筒身应变离心试验与数模结果对比图。
图7是钢圆筒筒身位移数值计算云图。
对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,可以根据以上附图获得其他的相关附图。
具体实施方式
为了进一步表述本发明的内容、计算效果,结合附图列举如下实例进行说明:
实施例1
S1、根据分析工况,建立对应的有限元分析模型,具体包括:土体模型和钢圆筒模型。并将土体模型划分为三个区域,包括:桩靴下压导致的柱孔区、钢圆筒内部开挖区和影响区。桩靴直径即柱孔区的直径为D1,桩靴贯入深度即柱孔区的深度为d,开挖区直径为D2,开挖区的深度为h,影响区的宽度为≥10D1,影响区的深度为≥10d;水下结构模型的外径为D2,水下结构模型的高度为H,水下结构模型的壁厚为t。
在该步骤中,算例参数如下,桩靴直径D1=18m,桩靴贯入深度d=3m;邻近水下结构为大直径薄壁钢圆筒,其直径D2=28m,高度H=18m,钢圆筒筒壁和横肋为Q345B型钢材,钢圆筒上环肋厚度为25mm,环肋宽度为260mm,环肋相对位置图如图3所示,钢圆筒边缘距离桩靴边缘为0.8D1;土体为均质黏土,钢圆筒筒内土体开挖深度h=9.2m,具体的钢圆筒及土体参数如表1所示。
表1计算参数表
S2、建立土体材料属性、钢圆筒材料属性,并将其分别赋予土体模型与钢圆筒模型。
S3、根据计算工况,确定桩靴与钢圆筒距离为0.8D1,钢圆筒埋置深度H为18m,将钢圆筒模型与土体模型装配至一起,形成装配体。
S4、建立相对应的计算分析步(静力、通用),包括:地应力分析步、钢圆筒安放分析步、钢圆筒开挖分析步、柱孔区形成分析步、柱孔区扩张分析步。
S5、建立相互作用。首先建立钢圆筒与土体的接触属性,采用面面接触,法向行为设置为“硬”接触,切向行为罚函数为0.4。然后使用相互作用中的单元生死功能,模拟钢圆筒的安放与内部土体开挖过程及桩靴下压产生柱孔区的过程。具体步骤为:使用单元生死功能将钢圆筒在地应力分析步设置为无效,并在钢圆筒安放分析步设置为重激活;将钢圆筒内部开挖区在开挖分析步设置为无效,以模拟该区域土体的开挖;将柱孔区在柱孔形成分析步中设置为无效,以模拟该区域土体由于桩靴下压形成柱孔区。
S6、建立模型的边界条件。对土体模型设置合理的位移边界,包括:对土体侧面进行位移约束,1方向和2方向设置为0;对土体底面设置位移约束,1方向、2方向、3方向均设置为0;对土体顶面设置位移约束,3方向设置为0。
S7、建立荷载条件。包括:①在地应力分析步中对全部土体施加重力,竖向分量设为10。②在钢圆筒安放分析步中对钢圆筒施加重力,竖向分量设为10。③在柱孔区扩张分析步中对区域内底部土体表面施加等同于基底压力大小的压强,具体数值为220kPa。④在柱孔区扩张分析步中对区域内侧向土体表面施加等同于扩孔力大小的压强,荷载分布形式沿深度方向为矩形分布,最大值大小由图2中柱孔扩张理论推导求得。根据柱孔扩张理论解,土体满足下述平衡微分方程:
桩靴在黏土中的贯入过程为不排水过程,土体材料服从Tresca屈服准则,即:
σr-σθ=2K (2)
式中:K为Tresca常数,K=su。
联立式(1)和式(2),并代入相应边界条件后积分可得:
S8、对土体施加地应力场以平衡地应力,由于土体高度为40m,因此土体底部地应力值为-380kPa。
S9、对土体模型及钢圆筒模型进行网格划分。土体及钢圆筒模型均采用实体单元模拟,网格属性为六面体单元C3D8R,网格加密区内网格尺寸为0.05D1。桩靴一倍直径范围外受影响较小,因此网格选用单精度,网格尺寸为0.05D1-0.2D1,具体的网格划分方式见图4。
S10、提交作业进行计算。
实施例2
本实施例针对算例工况开展了筒靴间距为0.8D1时桩靴贯入对邻近钢圆筒结构影响的离心机模型试验。试验用离心机为TK-C500型土工离心机。试验中按1:100的缩尺参照算例工况建立桩靴及钢圆筒模型。离心机试验能较为准确地测量筒壁应变及位移。试验时,为测量钢圆筒变形,在钢圆筒顶面布置四个位移测点F1-F4,测点数据由激光位移传感器采集,具体布置方式如图5(a)所示;在钢圆筒外壁沿环向方向每隔45°设置一组应变片,依次编号为1-8号,应变片为BX120-2AA型,具体布置方式如示意图5(b)所示。具体的土体制备方法及试验加载过程可参见文献:张海洋,刘润,贾沼霖.自升式平台插桩对邻近平台桩基础的影响研究[J].岩土工程学报,2021,43(05):867-876。
图6为使用本方法计算得到的受桩靴插拔影响邻近钢圆筒筒身应变结果与离心机试验结果对比图。由图6可以看出本发明提出的方法与离心机试验结果吻合较好,能较准确地考虑桩靴插拔对邻近结构物应力应变的影响。本方法计算得到的钢圆筒筒身位移云图如图7所示,其位移在钢圆筒顶部筒靴连线位置处达到最大,位移变形模式与离心机试验结果相同。具体数值与离心机试验结果对比如表2所示,由表2可知,数值计算位移结果与离心机试验位移结果相差1.70%,可见本发明提出的方法能较好地模拟桩靴插拔对邻近结构物位移的影响。
表2钢圆筒筒身位移对比
以上对本发明做了示例性的描述,应该说明的是,在不脱离本发明的核心的情况下,任何简单的变形、修改或者其他本领域技术人员能够不花费创造性劳动的等同替换均落入本发明的保护范围。
Claims (5)
1.一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,其特征在于:
S1、根据分析工况,建立对应的有限元分析模型,具体包括:土体模型和水下结构模型,土体模型划分为三个区域,包括:桩靴下压导致的柱孔区、水下结构的开挖区和影响区,桩靴直径即柱孔区的直径为D1,桩靴贯入深度即柱孔区的深度为d,开挖区直径为D2,开挖区的深度为h,影响区的宽度为≥10D1,影响区的深度为≥10d;水下结构模型的外径为D2,水下结构模型的高度为H,水下结构模型的壁厚为t;
S2、建立土体材料属性和水下结构材料属性,并将其分别赋予土体模型与水下结构模型;
S3、根据计算工况,将水下结构模型与土体模型装配至一起,形成装配体;
S4、建立数值分析模型的计算分析步;
S5、建立相互作用,首先建立水下结构与土体的接触属性,然后使用相互作用中的单元生死功能,模拟水下结构的安放与内部土体开挖过程及桩靴下压产生柱孔区的过程;
S6、建立模型的边界条件,对土体模型设置合理的位移边界;
S7、建立荷载条件,包括:①在地应力分析步中对全部土体施加重力,②在水下结构安放分析步中对水下结构施加重力,③在柱孔区扩张分析步中对区域内底部土体表面施加等同于基底压力大小的压强,④在柱孔区扩张分析步中对区域内侧向土体表面施加等同于扩孔力大小的压强;
S8、对土体施加地应力场;
S9、对土体模型及水下结构模型进行网格划分;
S10、提交作业进行计算。
2.根据权利要求1所述的一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,其特征在于:在S3中,需设置的参数包括:桩靴与水下结构模型之间的距离和水下结构模型的埋置深度。
3.根据权利要求1所述的一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,其特征在于:在S4中,共建5个分析步,包括:地应力分析步、水下结构安放分析步、水下结构开挖分析步、柱孔区形成分析步和柱孔区扩张分析步。
4.根据权利要求1所述的一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,其特征在于:在S5中,具体步骤为:使用单元生死功能将水下结构在地应力分析步设置为无效,并在水下结构安放分析步设置为重激活;将水下结构开挖区在开挖分析步设置为无效,以模拟该区域土体的开挖;将柱孔区在柱孔形成分析步中设置为无效,以模拟该区域土体由于桩靴下压形成柱孔区。
5.根据权利要求1所述的一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法,其特征在于:在S6中,位移边界包括:土体底部全约束、土体侧面放开沿深度方向的位移,其余全约束、土体顶面约束沿深度方向的位移。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210278842.8A CN114676610A (zh) | 2022-03-21 | 2022-03-21 | 一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210278842.8A CN114676610A (zh) | 2022-03-21 | 2022-03-21 | 一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114676610A true CN114676610A (zh) | 2022-06-28 |
Family
ID=82074910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210278842.8A Pending CN114676610A (zh) | 2022-03-21 | 2022-03-21 | 一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114676610A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117077250A (zh) * | 2023-07-17 | 2023-11-17 | 浙江大学 | 一种用于离心模型实验的考虑静动力相似模型桩设计方法 |
-
2022
- 2022-03-21 CN CN202210278842.8A patent/CN114676610A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117077250A (zh) * | 2023-07-17 | 2023-11-17 | 浙江大学 | 一种用于离心模型实验的考虑静动力相似模型桩设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Centrifuge modeling of lateral bearing behavior of offshore wind turbine with suction bucket foundation in sand | |
CN107894311B (zh) | 土石坝地震破坏的模型试验方法 | |
Khosravi et al. | Arch action over an excavated pit on a stable scarp investigated by physical model tests | |
CN111595703A (zh) | 一种基于模型试验的节理岩质边坡爆破失稳规律的测试方法 | |
CN114676610A (zh) | 一种适用于模拟插桩作业对邻近结构物影响的拟静力计算方法 | |
Chen et al. | Analytical solution for kinematic response of offshore piles under vertically propagating S-waves | |
Wang et al. | Undrained bearing capacity of suction caissons for offshore wind turbine foundations by numerical limit analysis | |
CN108681523B (zh) | 小净距隧道深、浅埋临界深度计算方法 | |
CN114372314A (zh) | 预测降压降水引起地面沉降的方法 | |
CN105756107A (zh) | 验证支护结构与土体共同作用的离心试验模型及制作方法 | |
CN110321569B (zh) | 一种适用于桩靴插拔对临近桩基影响的数值模拟方法 | |
CN113011062A (zh) | 一种变电站地基处理数值仿真分析方法 | |
Yuan et al. | Simulation of tensile cracking in earth structures with an adaptive RPIM-FEM coupled method | |
Zhou et al. | Engineering characteristics and reinforcement program of inclined pre-stressed concrete pipe piles | |
CN103088847B (zh) | 一种地基加固离心模型试验砂桩制作装置及制作方法 | |
Ghazavi et al. | Numerical investigation of load-settlement characteristics of multi-edge shallow foundations | |
CN115048611A (zh) | 一种深水无封底嵌岩钢板桩围堰嵌固深度的计算方法 | |
Zhuo et al. | Stability of long trench in soft soils by bentonite-water slurry | |
Fan et al. | Experimental and numerical study on the damage evolution of random rock joint surface during direct shear under CNL condition | |
CN113006084A (zh) | 一种无支撑基坑支档系统的施工方法及其动态设计方法 | |
CN109635489B (zh) | 一种考虑桩靴脚印卸荷作用对邻近桩基础承载力影响的有限元计算方法 | |
Li et al. | Load-displacement behaviour of tapered piles: theoretical modelling and analysis | |
Dey et al. | Numerical Study of the Effect of Pile Driving on the Position of a Neutral Plane | |
Moghaddam et al. | Mechanism of wave-induced liquefaction around suction caissons | |
CN114329827A (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 |