Loading [MathJax]/jax/output/SVG/fonts/TeX/Size1/Regular/Main.js
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

上限法简化模型分析鱼雷锚水平承载力

瑜璐, 杨庆, 张金利, 杨钢

瑜璐, 杨庆, 张金利, 杨钢. 上限法简化模型分析鱼雷锚水平承载力[J]. 岩土工程学报, 2020, 42(4): 773-781. DOI: 10.11779/CJGE202004022
引用本文: 瑜璐, 杨庆, 张金利, 杨钢. 上限法简化模型分析鱼雷锚水平承载力[J]. 岩土工程学报, 2020, 42(4): 773-781. DOI: 10.11779/CJGE202004022
YU Lu, YANG Qing, ZHANG Jin-li, YANG Gang. Simplified model for upper bound method to analyze horizontal bearing capacity of torpedo anchors[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 773-781. DOI: 10.11779/CJGE202004022
Citation: YU Lu, YANG Qing, ZHANG Jin-li, YANG Gang. Simplified model for upper bound method to analyze horizontal bearing capacity of torpedo anchors[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 773-781. DOI: 10.11779/CJGE202004022

上限法简化模型分析鱼雷锚水平承载力  English Version

基金项目: 

国家重点研发计划项目 2016YFE0200100

国家自然科学基金重点项目 51639002

详细信息
    作者简介:

    瑜璐(1983—),女,博士研究生,主要从事深海锚固系统研究。E-mail: yulu@mail.dlut.edu.cn

    通讯作者:

    杨庆, E-mail: qyang@dlut.edu.cn

  • 中图分类号: TU431

Simplified model for upper bound method to analyze horizontal bearing capacity of torpedo anchors

  • 摘要: 鱼雷锚的承载力问题是工业界关注的首要问题。基于塑性极限分析理论上限法提出一种简化模型用于计算鱼雷锚水平承载力。通过有限元方法得到考虑锚翼长宽比(Lf/Dw)和土体强度(su)等因素的平均水平承载力系数Np,avg=f(su,Lf/Dw),将其用于上限法简化模型。并采用有限元方法和de’Sousa研究结果验证上限法简化模型的合理性。研究进一步发现:①旋转中心位置低于锚翼长时,减小锚翼长度或宽度,旋转中心位置靠向锚尖;②鱼雷锚的锚翼长度大于1/2锚长时,增加锚翼宽度比增加锚翼长度更能显著提高承载力,该结论为优化锚型设计提供参考依据。
    Abstract: The bearing capacity of torpedo anchors is the major concern for the marine industry. In this study, a simple model is proposed to calculate their horizontal bearing capacity based on the upper-bound limit analysis of plasticity. The horizontal bearing factor,Np,avg=f(su,Lf/Dw), is obtained considering the influences of aspect ratio of torpedo anchor fluke(Lf/Dw) and soil strength (su) by using the finite element method (FEM), and it is applied in the simple model. The analytic results of the simple model are verified through the comparative study with FEM and de’Sousa results. It is found that: (1) When the rotation center position is lower than the fluke of torpedo anchors, the rotation center position is close to the anchor tip by reducing fluke width or fluke length. (2) When the fluke length is larger than 1/2 length of torpedo anchors, compared with that by increasing the fluke length, the bearing capacity is obviously improved by increasing the fluke width. It is helpful for the optimization design of torpedo anchor shape.
  • 近年随着海洋工业从浅海向深海发展,新型适用于深海系泊的鱼雷锚,因其安装简单方便、成本低等诸多优势深受海洋工业界的青睐。自2000年鱼雷锚首次应用于深海浮式平台系泊系统,到目前已有上千根鱼雷锚系泊采油平台和船支(主要应用于巴西坎普斯盆,挪威西海岸附近,墨西哥湾等海域)。鱼雷锚的承载力影响上部采油平台的生产能力,因此鱼雷锚的承载力问题是其研究的首要问题。研究人员希望通过试验和数值方法构建具有广泛指导意义的承载力包络线(VH曲线),预测鱼雷锚的承载力。但因鱼雷锚是深埋锚,锚尖埋深达1.5~2.5倍锚长(Lieng等[1],De Araujo等[2],Brandão等[3]),试验方法很难获得相对于锚端的水平拉拔承载力,因此更难获得承载力的VH曲线。O'Beirne等[4]通过数值方法获得小比尺模型锚的水平承载力,进一步构建承载力包络线(VH曲线),提出水平承载力是1.1倍的竖向承载力(Hmax=1.1Vmax)。但是该关系仅限于O'Beirne等研究的土体强度和锚型。Fu等[5]通过离心机试验和数值方法对比分析了3种锚型在高岭土中的承载力,并构建了新型承载力包络线,θ>45°时(θ为锚链与锚端水平面夹角)倾斜拉拔的竖向分力等于竖向承载力(Vi=Vmax),θ<25°时倾斜拉拔的水平分力等于水平承载力(Hi=Hmax),25°<θ<45°时承载力下限边界Vi/Vmax=1.71.4Hi/Hmax。Fu等水平承载力计算是基于Brom对桩的水平承载力研究,首先假设桩周围土体极限反力分布,再根据整体受力平衡来计算极限承载力[6]。但是因为鱼雷锚形状不规则锚周围土体反力分布很难准确假定。

    到目前为止,水平承载力的研究主要集中于数值方法,解析方法研究相对较少。本文基于塑性极限分析上限法提出一种简化模型用于推导鱼雷锚水平承载力。将有限元方法构建的平均水平承载力系数Np,avg=f(su,Lf/Dw)Lf/Dw锚翼长宽比;su不排水抗剪强度)用于简化模型,推导不同锚型和土体强度下鱼雷锚水平承载力的解析解,并与ABAQUS数值解以及前人研究结果比较分析,进一步验证简化模型合理性。

    Murff等[7]率先将塑性极限分析上限法用于三维桩水平承载力研究。根据桩的动荷载边界条件构建虚拟速度场,外力做功耗散功率等于桩周土体内能耗散功率,利用能量平衡原理计算载荷,不必解复杂的偏微分方程。虚拟速度场满足几何方程,体积不变和位移速度边界条件。因此更适合解决形状不规则的鱼雷锚的水平承载力问题。

    鱼雷锚锚链通常连接在锚端部,在水平荷载作用下围绕旋转中心发生旋转(图1)。假设土体为满足相关流动法则的理想刚塑性Tresca材料,锚端的旋转速度v0,锚侧速度分布随深度(z)成线性变化vi=(1z/L0)v0,(其中L0旋转中心与锚端距离)。

    图  1  鱼雷锚水平荷载受力图
    Figure  1.  Horizontal loads on a torpedo anchor
    ˙ES=˙Einter+˙Eshear+˙Efric=FHv0=v0D0L0|1zL0|(Np, avgsu+γz)dz (1)

    式中 ˙Einter土体内塑性变形耗散功率(软化,硬化[8-10],剪胀,剪缩等);˙Eshear速度不连续面剪切功率;˙Efric界面摩擦所消耗功率;FH水平外荷载;v0锚端速度;Np, avg平均水平承载力系数;γ土浮重度;su不排水抗剪强度;L0旋转中心距离锚顶部的距离;L锚的长度;D0锚的有效直径,当无锚翼时D0=DD为锚杆直径);当有锚翼时D0=Dw=D+2Wf(β=0)D0=D+2Wfsinβ(β0)Wf为锚翼宽度,β为拉拔力施加方向与锚翼夹角,如图2所示)。为减小鱼雷锚安装过程中锚翼端部的贯入阻力,通常将锚翼底部和顶部设计成斜面,因此锚翼底部和顶部的有效直径Dw(i)=D+2Wfsinβh(i)/Lf(i),h(i)为锚翼中某点i高度距离锚翼顶部或底部的距离,Lf(i)是锚翼上斜面的竖向投影高度Lf1或者锚翼下斜面的竖向投影高度Lf3

    图  2  鱼雷锚水平荷载加载方向示意图
    Figure  2.  Horizontal loading diagram of a torpedo anchor

    锚的单位面积水平向承载力计算公式pmax=Npsu+γzNp水平承载力系数)。根据Randoph等[11]对桩的水平承载力的研究,深埋桩土体围绕桩发生平面流动,近似为平面应变状态,土重对桩水平承载力影响较小,可忽略土重影响。因此对于深埋鱼雷锚可忽略土重的影响,可简化为

    ˙ES=FHv0=v0D0L0|1zL0|Np, avgsudz (2)

    鱼雷锚旋转中心位置存在3种情况,如图1,可能位于锚翼中部A点(Lf1<L0<(LfLf3));旋转中心位于锚翼下斜面B点((LfLf3)<L0<Lf);旋转中心低于锚翼C点(L0>Lf)。

    (1)当旋转中心位于A点时(Lf1<L0<(LfLf3)

    ˙ES=v0Dw(i)Lf10(1zL0)Np, avgsudz+v0DwL0Lf1(1zL0)Np, avgsudz+v0DwLfLf3L0(zL01)Np, avgsudz+v0Dw(i)LfLfLf3(zL01)Np,avgsudz+v0DLLf(zL01)Np,avgsudz (3)

    (2)当旋转中心位于B点时((LfLf3)<L0<Lf

    ˙ES=v0Dw(i)Lf10(1zL0)Np, avgsudz+v0DwLfLf3Lf1(1zL0)Np, avgsudz+v0Dw(i)L0LfLf3(1zL0)Np, avgsudz+v0Dw(i)LfL0(zL01)Np, avgsudz+v0DLLf(zL01)Np, avgsudz (4)

    (3)当旋转中心位于C点时(L0>Lf

    ˙ES=v0Dw(i)Lf10(1zL0)Np, avgsudz+v0DwLfLf3Lf1(1zL0)Np, avgsudz+v0Dw(i)LfLfLf3(1zL0)Np, avgsudz+v0DL0Lf(1zL0)Np, avgsudz+v0DLL0(zL01)Np, avgsudz (5)

    当桩发生旋转,旋转中心到桩底的距离小于旋转中心到桩底部边缘的距离因此桩底部形成一个球形破坏面。

    简化球形破坏面耗散功率公式为[12]

    ˙ET=(π/4)v0suD20 (6)

    鱼雷锚锚尖是圆锥或椭圆形,在旋转时就不会在锚尖形成球形成破坏面。反而因为鱼雷锚是深埋锚,锚杆顶部会出现球形破坏面。

    鱼雷锚水平荷载作用下周围土体内部总损耗率:

    ˙E=˙ES+˙ET (7)

    鱼雷锚在水平荷载下发生旋转,旋转中心的位置将影响水平承载力的大小。通过优化旋转中心的位置得到上限法的最小水平承载力。根据数学求函数极值的方法,解公式(8)得到最优旋转中心位置。研究人员[13]通过刚性桩的横向拉拔试验研究发现,刚性桩绕旋转中心转动,失稳后旋转中心在桩上的位置基本不变。本文假设鱼雷锚在水平荷载作用力下失稳后旋转中心位置也是不变化。

    FH=d(˙Es/v0)dL0=0 (8)

    鱼雷锚无锚翼部分近似于桩,Randolph通过上下限方法研究桩在完全摩擦时水平承载力系数为11.94。锚翼部分近似于DEPLA(动力深埋平板锚dynamically embedded plate anchors)锚。O’loughlin等[14]通过离心机试验研究DEPLA锚的承载力系数14.2~15.8,并指出DEPLA锚的破坏面沿着交叉锚翼向外扩展的正交破坏面,因此其承载力系数大于Martin等[15]给出深埋极薄圆型平板锚承载力系数精确解13.11,以及Wang等[16]给出平板锚厚度是宽度的1/20时,平板锚的承载力系数14.33。鱼雷锚水平荷载的运动趋势与DEPLA锚略有不同,DEPLA锚在拉拔荷载下发生平动,鱼雷锚水平荷载下发生旋转,因此鱼雷锚的水平承载力系数应低于DEPLA锚的承载力系数而高于桩的水平承载力系数。

    鱼雷锚包括有锚翼和无锚翼部分,锚翼的长宽影响水平承载力系数,因此鱼雷锚很难获得单位面积水平承载力系数。本文通过有限元结果反推全锚的平均水平承载力系数,考虑土体强度和Lf/DwNp,avg=f(su,Lf/Dw)函数影响,其推导流程如下图3

    图  3  水平承载力系数计算流程图
    Figure  3.  Flow chart for calculating horizontal capacity factor

    众所周知鱼雷锚锚翼的长度(Lf)和宽度(Wf)对水平承载力影响很大。存在锚翼长度与有效直径比(Lf/Dw)相同但锚翼侧面积(Af)不同和锚翼的Af相同但Lf/Dw不同。如表1,W-0和N-4鱼雷锚的Lf/Dw基本相同,但是Af相差4.3倍。W-1和N-1鱼雷锚的Af相同,但是Lf/Dw却不同。因此针对此情况,本文设计了两类锚型:W型-鱼雷锚:保持锚翼长度不变,改变锚翼宽度,达到不同的Lf/Dw。N型-鱼雷锚:保持锚翼宽度不变,改变锚翼长度,从而实现Lf/Dw的变化。见图4图5

    表  1  鱼雷锚尺寸表
    Table  1.  Sizes of torpedo anchor
     锚型符号锚长锚直径锚尖长锚翼厚锚翼宽锚翼长长径比面积
    L/mD/mLt/mt/mWf/mLf/mLf1/mLf2/mLf3/mLf/DwAf/m2
    W型-鱼雷锚W-0120.961.820.11.915.920.483.841.61.2474.57
    W-1120.961.820.11.445.920.483.841.61.5456.22
    W-2120.961.820.10.835.920.483.841.62.2632.40
    W-3120.961.820.10.455.920.483.841.63.1817.57
    W-4120.961.820.10.395.920.483.841.63.4015.23
    N型-鱼雷锚N-0120.961.820.10.909.390.487.311.63.4060.12
    N-1120.961.820.10.908.850.486.771.63.2156.23
    N-2120.961.820.10.906.150.484.971.62.2336.79
    N-3120.961.820.10.904.350.484.071.61.5823.83
    N-4120.961.820.10.903.450.483.171.61.2517.35
    验证锚型T-0120.961.820.11.165.920.483.841.61.8145.29
    T-1120.961.820.10.655.920.483.841.62.6325.38
    T-2120.961.820.10.907.050.480.681.62.5543.27
    T-3120.961.820.10.905.250.482.271.61.9030.31
    注:面积(Af)是锚翼侧面积之和。
    下载: 导出CSV 
    | 显示表格
    图  4  鱼雷锚锚型示意图
    Figure  4.  Schematic diagram of torpedo anchors
    图  5  锚翼部分侧面与Lf/Dw关系曲线
    Figure  5.  Relationship between fluke area of torpedo anchor and Lf/Dw

    所有锚型预埋在同样深度,锚尖埋深为1.5L,锚顶部埋深0.5L。土体被考虑为各向同性的纯黏土,满足Tresca屈服准则,浮重度为6 kN/m3。基于ABAQUS软件建立鱼雷锚预埋于土中3D数值模型(图6),土体模型尺寸长和宽为20.83D,深为27.08D,弹性模量为500Su,泊松比为0.49。采用C3D8R八面体单元对模型进行剖分,通过网格密度对承载力影响验证,确定网格加密区网单元尺寸为D/9.6,土体总网格数量在140000~230000。锚–土界面采用完全摩擦。

    图  6  1/2有限元模型
    Figure  6.  1/2 model of finite element model

    (1)均匀土体,土体强度su=10 kPa;

    (2)正常固结土,土体强度su=kz kPa (k =0.1,1,2,3,4,5,6),土体强度分布如图7所示。

    图  7  土体强度曲线
    Figure  7.  Curves of soil strength

    在水平面上,圆桩是完全对称将不受拉拔方向的影响。而鱼雷锚因带有锚翼不是完全对称结构,拉拔力施加方向与锚翼夹角(β),将导致有锚翼部分有效直径不同。如图2,当β=0°时,有效直径最大Dw=D+2Wf,当β=45°时,有效直径最小Dw=D+2Wfsin45°。De Sousa等[17]指出β=0°与β=45°的水平承载力差值是2.61%。图8(a),(b)所示也进一步证明,在均匀土体和正常固结土中锚翼拉拔夹角度对水平承载力略有影响。拉拔夹角β=45°和β=0°的W型–鱼雷锚或N型–鱼雷锚的水平承载力差值1.5%~5.5%。本次研究针对鱼雷锚有效直径等于最大有效直径,即β=0°的工况。

    图  8  β=45°和β=0°鱼雷锚水平承载力
    Figure  8.  Horizontal loads on torpedo anchors under β=45°and β=

    图8(a),(b)W-1鱼雷锚和N-1鱼雷锚侧面积相等,W-1鱼雷锚的水平承载力大于N-1鱼雷锚的水平承载力,可知侧面积相等时,锚翼宽度越宽水平承载力越大;W-3鱼雷锚和N-4鱼雷锚的Lf/Dw分别是3.18和1.25且W-3鱼雷锚的侧面积Af=17.57略大于N-4鱼雷锚的侧面积Af=17.04,但是N-4鱼雷锚的水平承载力大于W-3鱼雷锚的水平承载力。因为N-4鱼雷锚的等效直径Dwβ=0°时Dw=2.76β=45°时Dw=2.23)大于W-3鱼雷锚的Dwβ=0°时Dw=1.86β=45°时Dw=1.6)。由此可见,将锚翼宽度增加一倍,即使锚翼面积减少3%,也可以获得较高的水平承载力。

    图9(a),(b)为均匀土体中鱼雷锚在水平荷载作用下的位移云图和应力云图。由图9(a)可知,鱼雷锚在水平荷载作用下围绕旋转中心(位移为0)发生旋转,锚体位移随着深度基本呈线性变化,由此可见前文假设锚侧速度分布随深度成线性变化是合理的。

    图  9  均匀土体的应力云图和位移云图
    Figure  9.  Stress patterns and displacements of uniform soils

    图10(a),在均匀土中(su=10 kPa),W型-鱼雷锚保持锚翼长度不变,增加或减少锚翼宽度,水平承载力系数Nk=0p,avg随着Lf/Dw增加没有明显变化,则水平承载力系数取均值Nk=0p,avg=12.97。N型-鱼雷锚保持锚翼宽度不变,改变锚翼长度,Nk=0p,avg随着Lf/Dw增加逐渐减小,水平承载力系数与Lf/Dw关系见下式:

    图  10  水平承载力系数或旋转中心与Lf/Dw关系曲线
    Figure  10.  Relationship between horizontal capacity factor or rotation center and Lf/Dw
    Nk=0p,avg=14.22.16e(Lf/Dw4.73)28 (9)

    为使水平承载力系数Nk=0p,avgLf/Dw关系在工程实践中具有广泛通用性,并避免过高预测水平承载力,取两水平承载力系数的最小值:

    Nk=0p,avg=Min{12.9714.22.16e(Lf/Dw4.73)28 (10)

    表2比较了在均匀土体(su=10 kPa)中N-0~N-4鱼雷锚和W-0~W-4型鱼雷锚的数值解和简化模型解析解。当N型-鱼雷锚Lf/Dw较小时和W型-鱼雷锚Lf/Dw较大时,计算的水平承载力偏保守,解析解比数值解低5.7%左右。

    表  2  鱼雷锚水平承载力数值解与解析解对比
    Table  2.  Comparison of horizontal capacity between FEM results and analytical solutions
     锚型su=10 kPasu=1z kPa
    数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%
    W-02032.702054.42-1.071806.921847.862.27
    W-11766.031755.14-0.621556.131563.070.45
    W-21311.341309.16-0.171154.111172.851.62
    W-31000.55975.32-2.52882.76900.091.96
    W-4958.70917.55-4.29843.42847.720.51
    N-01420.921416.60-0.301283.301265.55-1.38
    N-11410.731404.80-0.421278.231269.89-0.65
    N-21395.331377.29-1.291232.931235.570.21
    N-31285.381215.69-5.421069.061068.05-0.09
    N-41153.281087.87-5.67940.64950.431.04
    下载: 导出CSV 
    | 显示表格

    在正常固结土(su=1z kPa)中,如图10(a),W型-鱼雷锚与N型-鱼雷锚的水平承载力系数Np,avg基本吻合并且随着Lf/Dw增加略有降低。土体强度su=1z kPa时,水平承载力系数Nk=1p,avg

    Nk=1p,avg=10+0.61+10(6.4+2LfDw) (11)

    表2中,比较在正常固结土(su=1z kPa)中N型-鱼雷锚和W型-鱼雷锚的数值解和简化模型解析解,两者差值均小于2.3%。

    图10(b)可知,对于N型-鱼雷锚,当旋转中心位置低于锚翼长时(L0>Lf),随着Lf/Dw增加,旋转中心靠近锚顶。当旋转中心位置高于锚翼长时(L0<Lf),随着Lf/Dw增加,旋转中心向锚尖靠近。也就是锚翼宽度不变,锚翼越短和越长旋转中心都向锚尖靠近。对于W型-鱼雷锚,无论旋转中心位置高于锚翼长还是低于锚翼长,都随着Lf/Dw增加,旋转中心越靠近锚尖。即锚翼长度不变,减小锚翼宽度,旋转中心向锚尖方向移动。

    图11所示N-1鱼雷锚的水平承载力与土体强度关系曲线。N-1鱼雷锚的水平承载力随着拉拔位移(U)增加逐渐增大,当U/D=0.2时水平承载力基本达到稳定,并且水平承载力随着土体强度增加而增大。图12所示,当土体强度梯度0.01<k<0.1时,水平承载力系数随着k值增大急剧增加;当土体强度梯度0.1<k<1时,随着k值增长水平承载力系数的增长逐渐缓慢;当土体强度梯度k >1时,随着k值增长水平承载力系数反而逐渐减小。下式表达鱼雷锚水平承载力系数与土体强度关系。值得注意,N-1鱼雷锚的旋转中心位置没有随着土体强度梯度的增加发生变化。

    图  11  水平承载力与U/D关系曲线
    Figure  11.  Relationship between horizontal capacity and U/D
    图  12  水平承载力系数和旋转中心与土强度梯度关系
    Figure  12.  Relationship between soil strength gradient and horizontal capacity factor or rotation center
    Nk=ip, avg={Nk=1p,avgk1.50.00052+k1.5(0.01<k1)Nk=1p,avg1.6e(k6)23.38(k>1)  (12)

    综上所述,得到水平承载力系数与土体强度和Lf/Dw的关系式如下:

    Np,avg={Min{12.9714.22.16e(Lf/Dw4.73)28(k=0)[10+0.61+10(6.4+2LfDw)]k1.50.00052+k1.5(0.01<k1)[10+0.61+10(6.4+2LfDw)]1.6e(k6)23.38(k>1)  (13)

    将上文中T-1~T-4型鱼雷锚用于验证推导的水平承载力系数Np, avg的可行性。分别将T-1~T-4鱼雷锚在均匀土体(su=10 kPa, 50 kPa)和正常固结土体(su=1z kPa, 5z kPa中的数值解和简化模型解析解进行比较,结果如表3所示,T-1~T-4数值结果和简化模型的解析解差值均小于3.2%,由此可见Np, avg用于水平承载力计算是可行的。

    表  3  验证锚水平承载力的数值解与简化模型解析解比较
    Table  3.  Comparison of horizontal capacity between FEM results and analytical solutions of simplified models of torpedo anchors
    土体强度su=10 kPasu=50 kPasu=1z kPasu=5z kPa
    名称单位数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%
    T-11563.491575.160.757817.467769.07-0.621381.691387.070.396178.646142.38-0.59
    T-21168.061178.480.895840.325800.86-0.681025.471052.092.594673.034652.62-0.44
    T-31401.441414.520.937007.226973.86-0.481270.511277.190.525576.405651.631.35
    T-41350.361327.22-1.716751.796541.69-3.111150.691164.111.175303.295152.94-2.83
    下载: 导出CSV 
    | 显示表格

    De Souse等[17]应用数值方法研究了锚翼宽度,土体强度等因素对鱼雷锚承载力的影响。为了进一步论证本文研究的水平承载力系数和上限法简化模型预测鱼雷锚水平承载力的合理性。用上限法简化模型计算De’Sousa的工况水平承载,并与De Sousa的结果进行对比,如表4所示,当Lf/Dw从3.83增长到7.93时De Sousa的结果与简化模型解析解的差值小于2.4%。鱼雷锚的Lf/Dw=3.83,土体强度梯度k从1.5增长到6时,De’Sousa的结果与简化模型解析解的差值随着k值增加从0.68%增加到5.24%。由此可见,该简化锚型预测较低土体强度梯度k时差值更小,而海洋土体强度梯度范围常在k=0.5~2之间,很少出现土体强度梯度k>3的工况。

    表  4  De Sousa等[17]的结果与简化模型解析解比较
    Table  4.  Comparison between de’sousa FEM results and analytical solutions of simplified model
    名称单位锚直径/m锚翼宽/m锚翼长/mLf/Dw 旋转中心/m数值解/MN解析解/MN差值/%
    su=3z kPa1.10.1511.17.9312.766.76.862.39
    1.10.4511.15.5511.968.918.950.45
    1.10.911.13.8310.6111.7211.841.02
    su=1.5z kPa1.10.911.13.8310.615.95.940.68
    su=6z kPa1.10.911.13.8310.6120.9819.885.24
    下载: 导出CSV 
    | 显示表格

    通过上述比较分析,上限法简化模型对水平承载力预测是合理的。将该简化模型预测的水平承载力(Hmax)与API公式预测的竖向承载力(Vmax)带入前人研究的V-H包络线(FicosθHmax)m+(FisinθVmax)n=1θ是锚链与锚端水平面夹角,m,n是系数),可预测任意倾角θ的倾斜承载力Fi。鱼雷锚现场安装过程中根据贯入深度和土体强度就可预测倾斜极限承载力。再依据API基础设计提供的竖向和水平向承载力安全系数2.0和1.5,可得到任意θ角的设计承载力[18]。上述方法节省了有限元计算的时间和经济成本,适用于工程实践。

    图13可见,当N-型鱼雷锚的锚翼侧面积(Af)小于37 m2即锚翼长度小于1/2锚长时,锚的水平承载力随着锚翼长度增加快速增大。当锚翼侧面积(Af)大于37 m2即锚翼长度大于1/2锚长时,增大锚翼长度将不能显著提高锚的水平承载力。而W-型鱼雷锚的水平承载力一直随着锚翼侧面积增大成线性增大。即锚翼长度不变,增加锚翼宽度可以有效提高水平承载力。该结论为锚型优化提供参考依据。

    图  13  鱼雷锚水平承载力与锚翼侧面积关系
    Figure  13.  Relationship between lateral area of fluke of torpedo anchor and horizontal capacity

    本文基于塑性极限分析上限法建立鱼雷锚水平承载力解析解的简化模型。用于海洋工业界预测鱼雷锚水平承载力。

    (1)拉拔角度对水平承载力略有影响,锚翼面积相等时锚翼越宽水平承载力越大。

    (2)均匀土体,鱼雷锚锚翼宽度不变改变锚翼长度,水平承载力系随着Lf/Dw增加而减小。鱼雷锚锚翼长度不变,改变锚翼宽度,水平承载力系数随着Lf/Dw增加没有明显变化,Nk=0p,avg=12.97

    (3)正常固结土,水平承载力系数随Lf/Dw增大略有降低。随着土体强度梯度k增加,水平承载力系数先增大后减小,但是旋转中心位置不随k值变化。

    (4)当锚翼宽度不变旋转中心小于锚翼长度时,锚翼越长旋转中心都越靠近锚尖。当旋转中心小于锚翼长度时,锚翼越短旋转中心都越靠近锚尖。锚翼长度不变,减小锚翼宽度旋转中心靠近锚尖。

    (5)鱼雷锚的锚翼长度大于1/2锚长时,增加锚翼宽度比增加锚翼长度更能明显提高承载力。

  • 图  1   鱼雷锚水平荷载受力图

    Figure  1.   Horizontal loads on a torpedo anchor

    图  2   鱼雷锚水平荷载加载方向示意图

    Figure  2.   Horizontal loading diagram of a torpedo anchor

    图  3   水平承载力系数计算流程图

    Figure  3.   Flow chart for calculating horizontal capacity factor

    图  4   鱼雷锚锚型示意图

    Figure  4.   Schematic diagram of torpedo anchors

    图  5   锚翼部分侧面与Lf/Dw关系曲线

    Figure  5.   Relationship between fluke area of torpedo anchor and Lf/Dw

    图  6   1/2有限元模型

    Figure  6.   1/2 model of finite element model

    图  7   土体强度曲线

    Figure  7.   Curves of soil strength

    图  8   β=45°和β=0°鱼雷锚水平承载力

    Figure  8.   Horizontal loads on torpedo anchors under β=45°and β=

    图  9   均匀土体的应力云图和位移云图

    Figure  9.   Stress patterns and displacements of uniform soils

    图  10   水平承载力系数或旋转中心与Lf/Dw关系曲线

    Figure  10.   Relationship between horizontal capacity factor or rotation center and Lf/Dw

    图  11   水平承载力与U/D关系曲线

    Figure  11.   Relationship between horizontal capacity and U/D

    图  12   水平承载力系数和旋转中心与土强度梯度关系

    Figure  12.   Relationship between soil strength gradient and horizontal capacity factor or rotation center

    图  13   鱼雷锚水平承载力与锚翼侧面积关系

    Figure  13.   Relationship between lateral area of fluke of torpedo anchor and horizontal capacity

    表  1   鱼雷锚尺寸表

    Table  1   Sizes of torpedo anchor

     锚型符号锚长锚直径锚尖长锚翼厚锚翼宽锚翼长长径比面积
    L/mD/mLt/mt/mWf/mLf/mLf1/mLf2/mLf3/mLf/DwAf/m2
    W型-鱼雷锚W-0120.961.820.11.915.920.483.841.61.2474.57
    W-1120.961.820.11.445.920.483.841.61.5456.22
    W-2120.961.820.10.835.920.483.841.62.2632.40
    W-3120.961.820.10.455.920.483.841.63.1817.57
    W-4120.961.820.10.395.920.483.841.63.4015.23
    N型-鱼雷锚N-0120.961.820.10.909.390.487.311.63.4060.12
    N-1120.961.820.10.908.850.486.771.63.2156.23
    N-2120.961.820.10.906.150.484.971.62.2336.79
    N-3120.961.820.10.904.350.484.071.61.5823.83
    N-4120.961.820.10.903.450.483.171.61.2517.35
    验证锚型T-0120.961.820.11.165.920.483.841.61.8145.29
    T-1120.961.820.10.655.920.483.841.62.6325.38
    T-2120.961.820.10.907.050.480.681.62.5543.27
    T-3120.961.820.10.905.250.482.271.61.9030.31
    注:面积(Af)是锚翼侧面积之和。
    下载: 导出CSV

    表  2   鱼雷锚水平承载力数值解与解析解对比

    Table  2   Comparison of horizontal capacity between FEM results and analytical solutions

     锚型su=10 kPasu=1z kPa
    数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%
    W-02032.702054.42-1.071806.921847.862.27
    W-11766.031755.14-0.621556.131563.070.45
    W-21311.341309.16-0.171154.111172.851.62
    W-31000.55975.32-2.52882.76900.091.96
    W-4958.70917.55-4.29843.42847.720.51
    N-01420.921416.60-0.301283.301265.55-1.38
    N-11410.731404.80-0.421278.231269.89-0.65
    N-21395.331377.29-1.291232.931235.570.21
    N-31285.381215.69-5.421069.061068.05-0.09
    N-41153.281087.87-5.67940.64950.431.04
    下载: 导出CSV

    表  3   验证锚水平承载力的数值解与简化模型解析解比较

    Table  3   Comparison of horizontal capacity between FEM results and analytical solutions of simplified models of torpedo anchors

    土体强度su=10 kPasu=50 kPasu=1z kPasu=5z kPa
    名称单位数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%数值解/kN解析解/kN差值/%
    T-11563.491575.160.757817.467769.07-0.621381.691387.070.396178.646142.38-0.59
    T-21168.061178.480.895840.325800.86-0.681025.471052.092.594673.034652.62-0.44
    T-31401.441414.520.937007.226973.86-0.481270.511277.190.525576.405651.631.35
    T-41350.361327.22-1.716751.796541.69-3.111150.691164.111.175303.295152.94-2.83
    下载: 导出CSV

    表  4   De Sousa等[17]的结果与简化模型解析解比较

    Table  4   Comparison between de’sousa FEM results and analytical solutions of simplified model

    名称单位锚直径/m锚翼宽/m锚翼长/mLf/Dw 旋转中心/m数值解/MN解析解/MN差值/%
    su=3z kPa1.10.1511.17.9312.766.76.862.39
    1.10.4511.15.5511.968.918.950.45
    1.10.911.13.8310.6111.7211.841.02
    su=1.5z kPa1.10.911.13.8310.615.95.940.68
    su=6z kPa1.10.911.13.8310.6120.9819.885.24
    下载: 导出CSV
  • [1]

    LIENG J T, TJELTA T I, SKAUGSET K. Installation of two prototype deep penetrating anchors at the GjØa Field in the North Sea[C]//Proceedings of the Offshore Technology Conference, 2010, Houston.

    [2]

    DE ARAUJO J B, MACHADO R D, DE MEDEIROS JUNIOR C J. High holding power torpedo pile – results for the first long term application[C]//Proceedings of the 23rd International Conference on Offshore Mechanics and Arctic Engineering, 2004, New York.

    [3]

    BRANDÃO F E N, HENRIQUES C C D, DE ARAUJO J B, et al. Albacora Leste field development – FPSO P-50 mooring system concept and installation[C]//Proceedings of the 38th Offshore Technology Conference, 2006, Houston.

    [4]

    O'BEIRNE C, O'LOUGHLIN C D, WANG D, et al. Capacity of dynamically installed anchors as assessed through field testing and three-dimensional large-deformation finite element analyses[J]. Canadian Geotechnical Journal, 2015, 52(5): 548-562. doi: 10.1139/cgj-2014-0209

    [5]

    FU Y, ZHANG X, LI Y, et al. Holding capacity of dynamically installed anchors in normally consolidated clay under inclined loading[J]. Canadian Geotechnical Journal, 2017, 54(9): 1257-1271. doi: 10.1139/cgj-2016-0537

    [6]

    BROMS B B. Lateral resistance of piles in cohesive soils[J]. Journal of the Soil Mechanics and Foundations Division, ASCE, 1964, 90(2): 27-64. doi: 10.1061/JSFEAQ.0000611

    [7]

    MURFF J D, HAMILTON J M. P-ultimate for undrained analysis of laterally loaded piles[J]. Journal of Geotechnical Engineering, 1993, 119(1): 91-107. doi: 10.1061/(ASCE)0733-9410(1993)119:1(91)

    [8]

    YAO Y P, HOU W, ZHOU A N. UH model: Three-dimensional unified hardening model for overconsolidated clays[J]. Géotechnique, 2009, 59(5): 451-469. doi: 10.1680/geot.2007.00029

    [9]

    YAO Y P, ZHOU A N. Non-isothermal unified hardening model: a thermo-elasto-plastic model for clays[J]. Géotechnique, 2013, 63(15): 1328-1345. doi: 10.1680/geot.13.P.035

    [10]

    YAO Y P, SUN D A, MATSUOKA H. A unified constitutive model for both clay and sand with hardening parameter independent on stress path[J]. Computers and Geotechnics, 2008, 35(2): 210-222. doi: 10.1016/j.compgeo.2007.04.003

    [11]

    RANDOLPH M F, HOULSBY G T. The limiting pressure on a circular pile loaded laterally in cohesive soil[J]. Géotechnique, 1984, 34(4): 613-623. doi: 10.1680/geot.1984.34.4.613

    [12]

    AUBENY C P, MOON S K, MURFF J D. Lateral undrained resistance of suction caisson anchors[J]. Int J Offshore Polar Engineering, 2001, 11(2): 95-103.

    [13] 崔新壮, 丁桦, 金青, 等. 刚性桩横向承载力的三维极限分析[J]. 岩石力学与工程学报, 2006, 25(3): 641-646. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200603035.htm

    CUI Xin-jian, DING Hua, JIN Qing, et al. Three- dimensional limit analysis of lateral bearing capacity of rigid plies[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(3): 641-646. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200603035.htm

    [14]

    O'LOUGHLIN C D, BLAKE A P, RICHARDSON M D. Installation and capacity of dynamically embedded plate anchors as assessed through centrifuge tests[J]. Ocean Engineering, 2014, 88(4): 204-213.

    [15]

    MARTIN C M, RANDOLPH M F. Application of the lower bound and the upper bound theorems of plasticity to collapse of circular foundations[C]//Proceedings of the 10th International Conference of the International Association for Computer Methods and Advances in Geomechanics, 2001, Tucson: 1417-1428.

    [16]

    WANG D, YUXIA HU, RANDOLPH , et al. Three- dimensional large deformation finite-element analysis of plate anchors in uniform clay[J]. Journal of Geotechnical & Geoenvironmental Engineering, 2010, 136(2): 355-365.

    [17]

    DE SOUSA, J R M, DE AGUIAR C S, ELLWANGER G B, et al. Undrained load capacity of torpedo anchors embedded in cohesive soils[J]. Journal of Offshore Mechanics and Arctic Engineering, 2011, 133(2): 1-12.

    [18]

    American Petroleum Institute. Geotechnical and Foundation Design Consideration[M]. 1st ed. Washington D C: API Recommended Practice 2 GEO, 2011.

  • 期刊类型引用(2)

    1. 吴学震,姜杰,李大勇,蒋宇静. 深水注浆锚设计原理及其承载特性试验研究. 岩土力学. 2022(10): 2707-2716 . 百度学术
    2. 刘金龙,祝磊,姚军,肖赟. 基于鱼雷锚施工的海洋组合式锚泊基础. 安徽建筑. 2021(05): 92-94 . 百度学术

    其他类型引用(5)

图(13)  /  表(4)
计量
  • 文章访问数:  281
  • HTML全文浏览量:  53
  • PDF下载量:  132
  • 被引次数: 7
出版历程
  • 收稿日期:  2019-03-20
  • 网络出版日期:  2022-12-07
  • 刊出日期:  2020-03-31

目录

/

返回文章
返回