一种水平力下钻井船插桩对邻近群桩影响的数值分析方法与流程

本发明涉及海洋工程技术领域，具体为一种水平力下钻井船插桩对邻近群桩影响的数值分析方法。

背景技术：

自升式钻井船在钻井过程中，往往要在导管架固定平台桩基础附近，将大直径的钻井船桩靴基础贯入海底泥面以下一定深度，以提供钻井船作业时需要的承载力。这种大尺寸桩靴在贯入土层时，不可避免对邻近平台的桩基稳定性产生不利影响，甚至危害平台的安全。因此，有必要对钻井船桩靴贯入对邻近平台桩基产生的影响进行分析评价。

评价钻井船桩靴贯入对邻近平台桩基影响的实质是能够客观分析插桩挤土导致的桩身位移与相应的弯矩。然而迄今为止，国内外关于钻井船插桩对邻近平台桩基影响的研究较少，已有的一些研究主要集中在模型试验，由于1g条件下的模型试验，有一定局限性，而离心模型试验的复杂性和过高的试验成本，故已有的这些试验研究只能模拟较为简单土层中的插桩对邻近桩的影响程度。

近年来，随着计算机硬件设施不断提高和有限元分析方法不断完善，为利用数值模拟方法解决钻井船插桩对邻近桩影响的问题提供了可能。但是通过常规有限元计算，分析插桩导致的邻近桩桩身应力、应变及位移的分布，可以避免桩靴贯入土体导致的非线性大变形引起的计算网格畸变，但无法模拟插桩的整个连续过程中不断累积的挤土效应对邻近桩的影响，计算结果与实际相比有很大差异。因此，采用非常规有限元数值模拟方法解决钻井船插桩对邻近桩影响问题是十分必要的。

目前，在解决非线性大变形问题的数值模拟方法有耦合欧拉-拉格朗日方法(cel方法)、任意拉格朗日-欧拉方法(ale方法)与网格重划分法。采用ale方法分析插桩对邻近平台桩基的影响，但分析得到的位移明显大于离心模型试验结果，有关研究表明采用cel方法对软黏土场地上钻井船插桩对桩头固定条件下邻近桩的影响离心模型试验进行数值仿真模拟，得到了与离心模型试验较为一致的弯矩分析结果。因此采用cel有限元方法模拟水平力作用下钻井船插桩对邻近的群桩影响，具有十分重要的工程意义。

技术实现要素：

本发明的目的在于提供一种水平力下钻井船插桩对邻近群桩影响的数值分析方法，以解决上述背景技术中提出的问题。

为实现上述目的，本发明提供如下技术方案：一种水平力下钻井船插桩对邻近群桩影响的数值分析方法，包括以下几个步骤：

s1、依据工程地质，建立符合实际的数值分析模型；

s2、在群桩桩头施加固定的水平荷载，然后模拟插桩过程，输出插桩前后桩体和土体的水平位移云图及桩体的水平位移曲线；

s3、根据桩体参数及api规范，建立对应的非线性地基梁模型；

s4、把s2得到的桩身位移输入到非线性地基梁模型中，输出插桩前后弯矩沿桩身的分布曲线；

s5、对比插桩前后的模拟结果，对水平力作用下钻井船插桩对邻近群桩影响程度进行评价。

作为本发明一种优选的技术方案，步骤s1中，所述依据工程地质，建立符合实际的数值分析模型，具体包括以下几个步骤：

s11、考虑到桩靴插入土层时桩靴周围会产生隆起的土体，为模拟该部分，故在土体表面以上设置一定深度的空穴单元；插桩深度影响范围内的土体设置成为欧拉材料，插桩深度影响范围以外的土体、桩靴及受影响桩均采用拉格朗日实体建模；桩靴贯入土层过程中的自身变形很小，故将其设为刚体；桩靴初始位置位于空穴底部的泥面中心处；

s12、根据土层参数、桩靴及受影响桩的力学参数，分别给土体、桩靴及桩赋值属性；

s13、在设置接触时，欧拉材料与拉格朗日体之间的接触采用“通用接触”的接触算法，拉格朗日体之间的接触采用“主面—从面”的接触算法，选择刚度大的桩体表面作为主面，刚度小的土体表面作为从面；

s14、根据实际情况，对欧拉土体四周施加水平两个方向速度为0的约束，土体底面施加竖直方向速度为0的约束，并且设置欧拉吸收边界，以消除由于计算模型尺寸选取导致的边界效应，减小计算结果的振荡；对拉格朗日土体按照相同的方式施加位移为0的约束；

s15、划分网格时，土体、空穴和桩均采用六面体单元进行划分，桩靴采用四面体单元进行划分。

作为本发明一种优选的技术方案，所述的步骤s11中，欧拉土体的深度应大于桩靴的插入深度，特别对于插桩深度范围内土层均为黏土层的情况，优先使用以下规则设置cel有限元模型中欧拉土体的范围：

①若插桩深度小于0.75倍桩靴直径，将欧拉土体的范围设置为泥面以下1～1.25倍桩靴直径；

②若插桩深度大于0.75倍桩靴直径，则将插桩深度与插桩深度以下0.5倍桩靴直径范围内的土体均设置成欧拉土体。

作为本发明一种优选的技术方案，所述的步骤s12中，根据土层参数、桩靴及受影响桩的力学参数，分别给土体、桩靴及桩赋值属性，具体如下：

由于在利用cel有限元方法模拟钻井船插桩时，利用理想弹塑性模型来描述土层的应力应变关系，遵循摩尔-库仑强度准则，故土体采用摩尔-库伦本构模型，包括四个基本参数：弹性模量、泊松比、粘聚力和内摩擦角；

对于受影响桩，桩的材料取为弹性，并按广义虎克定律确定应力-应变关系；

设定桩体为钢管桩，为了提高计算效率，将受影响桩设置为实体桩体，计算中，按钢管桩与实体桩抗弯刚度相等的原则，由下式确定计算模型中实心桩的等效弹性模量，式中：e1、e2分别为钢管桩和等效实体桩的弹性模量；i1、i2为钢管桩和等效实体桩的惯性矩。

作为本发明一种优选的技术方案，所述的步骤s15中，为了减少桩靴贯入阻力的振荡对结果的影响，在靠近桩靴部分设置为细网格区域，以外区域设置为粗网格区域，细网格尺寸取为0.05倍桩靴直径，细网格范围取为2倍桩靴直径。

作为本发明一种优选的技术方案，步骤s2中，所述在群桩桩头施加固定的水平荷载，然后模拟插桩过程，输出插桩前后桩体和土体的水平位移云图及桩体的水平位移曲线，具体包括以下几个步骤：

s21、在模型中群桩桩头建立参考点并与桩头顶面耦合，若桩头有承台时，也可在承台表面建立参考点并于表面耦合；

s22、根据桩头或承台受到的水平力大小，在参考点上加上相应的水平荷载；

s23、在桩靴底面建立参考点，并把桩靴约束为刚体；

s24、建立分析步，分别为重力分析步、加载分析步、插桩分析步，并分别设定各个分析步的时间，其中插桩分析步时间要根据插桩深度及插桩速度确定；

s25、通过在边界条件中修改桩靴在插桩分析步中的竖向位移，控制插桩深度；

s26、建立作业并提交计算，在结果中输出插桩前后的桩体和土体的水平位移云图及桩体的水平位移曲线。

作为本发明一种优选的技术方案，步骤s22及s25中，所述的力和位移的施加应避免瞬间施加，一般要为施加的力和位移建立一个平滑的幅值曲线，让力和位移的施加更加平稳，以消除激振力产生的震颤。

作为本发明一种优选的技术方案，步骤s3中，所述的根据桩体参数及api规范，建立对应的非线性地基梁模型，由于通过cel有限元方法进行数值分析，可以直接确定受影响桩的桩身位移，为了依据桩身位移，按土层的p-y关系确定相应的桩身弯矩，这里采用非线性地基梁有限元模型计算插桩导致的受影响桩的桩身弯矩，具体包括以下几个步骤：

s31、按实际桩的截面尺寸确定梁单元的截面特征参数；

s32、沿桩长设置若干非线性p-y弹簧，模拟桩-土之间相互作用；

s33、参照api规范，确定土层的p-y数据，并输入到地基梁中；

s34、根据桩头的约束情况，在地基梁中施加相应的约束。

作为本发明一种优选的技术方案，步骤s33中，所述的参照api规范，确定土层的p-y数据，对于软粘土，硬粘土和砂性土中桩的p-y曲线api规范分别给出了相应的计算公式：

(1)砂性土的p-y曲线

砂性土分为浅层破坏形式和深层破坏形式，浅层破坏和深层破坏的极限土抗力转折点深度xr按计算，式中：c1、c2、c3为系数，d为外径；

当x＜xr时，是浅层破坏；当x＞xr时，是深层破坏；

砂性土的桩侧的极限土抗力随深度不同而变化，浅层破坏和深层破坏极限土抗力为p，p取pux＝(c1x+c2d)γx和pud＝c3dγx两式中的较小值，式中：γ为土体有效重度；x为土表面下计算深度；

给定深度为x的砂性土p-y曲线用表示，式中：pu为深度x处的极限土抗力，a为考虑循环荷载或者静荷载状态的系数，k为地基反力系数，其中，循环荷载：a＝0.9，短期静载：

(2)软粘土中的p-y曲线

软粘土中的p-y曲线法一般假定pu沿深度分布，采用下列两则公式计算，

x＜xr，

x＞xr，pu＝9su；

式中：b为桩侧土计算宽度，这里一般采用桩径，γ为土的有效重度，su为粘土不排水抗剪强度，j为常数，一般粘土取0.5，稍硬粘土取0.25，x为从地面算起的任一深度；

联列以上两列公式，求得地表下土的强度降低的塑性变形深度xr：

作为本发明一种优选的技术方案，经验关系深度xr不可应用于土的强度变化不稳定的情况下，当土体重度和抗剪强度随深度变化时，xr的最小值约为桩径的2.5倍。

与现有技术相比，本发明的有益效果是：本发明一种水平力下钻井船插桩对邻近群桩影响的数值分析方法，依据地质条件，利用abaqus数值模拟软件建立数值分析模型并建立对应的非线性地基梁模型，通过模拟钻井船插桩过程，进而对比插桩前后桩身位移曲线、桩身弯矩曲线，更能直观合理的评价在水平力作用下钻井船插桩对邻近群桩的影响，且操作简单，易于实际应用。

附图说明

图1为系数c1，c2，c3的取值图；

图2为相对密度图；

图3为本发明钻井船插桩对邻近群桩影响的数值分析方法流程图；

图4为本发明实例群桩分析模型图；

图5为本发明实例非线性地基梁模型图；

图6为本发明实例插桩前桩体及土体水平位移云图；

图7为本发明实例插桩后桩体及土体水平位移云图；

图8为本发明实例插桩进行前，前桩桩体弯矩云图；

图9为本发明实例插桩完成后，前桩桩体弯矩云图；

图10为本发明实例插桩进行前，后桩桩体弯矩云图；

图11为本发明实例插桩完成后，后桩桩体弯矩云图；

图12为本发明实例插桩前、后桩体水平位移对比图；

图13为本发明实例插桩前、后桩体弯矩对比图。

具体实施方式

以下对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

实施例：如图1-13所示，本发明提供一种技术方案：一种水平力下钻井船插桩对邻近群桩影响的数值分析方法，包括以下几个步骤：

s1、依据工程地质，建立符合实际的数值分析模型；

s2、在群桩桩头施加固定的水平荷载，然后模拟插桩过程，输出插桩前后桩体和土体的水平位移云图及桩体的水平位移曲线；

s3、根据桩体参数及api规范，建立对应的非线性地基梁模型；

s4、把s2得到的桩身位移输入到非线性地基梁模型中，输出插桩前后弯矩沿桩身的分布曲线；

s5、对比插桩前后的模拟结果，对水平力作用下钻井船插桩对邻近群桩影响程度进行评价。

步骤s1中，所述依据工程地质，建立符合实际的数值分析模型，具体包括以下几个步骤：

s11、考虑到桩靴插入土层时桩靴周围会产生隆起的土体，为模拟该部分，故在土体表面以上设置一定深度的空穴单元；插桩深度影响范围内的土体设置成为欧拉材料，插桩深度影响范围以外的土体、桩靴及受影响桩均采用拉格朗日实体建模；桩靴贯入土层过程中的自身变形很小，故将其设为刚体；桩靴初始位置位于空穴底部的泥面中心处；

s12、根据土层参数、桩靴及受影响桩的力学参数，分别给土体、桩靴及桩赋值属性；

s13、在设置接触时，欧拉材料与拉格朗日体之间的接触采用“通用接触”的接触算法，拉格朗日体之间的接触采用“主面—从面”的接触算法，选择刚度大的桩体表面作为主面，刚度小的土体表面作为从面；

s14、根据实际情况，对欧拉土体四周施加水平两个方向速度为0的约束，土体底面施加竖直方向速度为0的约束，并且设置欧拉吸收边界，以消除由于计算模型尺寸选取导致的边界效应，减小计算结果的振荡；对拉格朗日土体按照相同的方式施加位移为0的约束；

s15、划分网格时，土体、空穴和桩均采用六面体单元进行划分，桩靴采用四面体单元进行划分。

所述的步骤s11中，欧拉土体的深度应大于桩靴的插入深度，特别对于插桩深度范围内土层均为黏土层的情况，优先使用以下规则设置cel有限元模型中欧拉土体的范围：

①若插桩深度小于0.75倍桩靴直径，将欧拉土体的范围设置为泥面以下1～1.25倍桩靴直径；

②若插桩深度大于0.75倍桩靴直径，则将插桩深度与插桩深度以下0.5倍桩靴直径范围内的土体均设置成欧拉土体。

所述的步骤s12中，根据土层参数、桩靴及受影响桩的力学参数，分别给土体、桩靴及桩赋值属性，具体如下：

由于在利用cel有限元方法模拟钻井船插桩时，利用理想弹塑性模型来描述土层的应力应变关系，遵循摩尔-库仑强度准则，故土体采用摩尔-库伦本构模型，包括四个基本参数：弹性模量、泊松比、粘聚力和内摩擦角；

对于受影响桩，桩的材料取为弹性，并按广义虎克定律确定应力-应变关系；

设定桩体为钢管桩，为了提高计算效率，将受影响桩设置为实体桩体，计算中，按钢管桩与实体桩抗弯刚度相等的原则，由下式确定计算模型中实心桩的等效弹性模量，式中：e1、e2分别为钢管桩和等效实体桩的弹性模量；i1、i2为钢管桩和等效实体桩的惯性矩。

所述的步骤s15中，为了减少桩靴贯入阻力的振荡对结果的影响，在靠近桩靴部分设置为细网格区域，以外区域设置为粗网格区域，细网格尺寸取为0.05倍桩靴直径，细网格范围取为2倍桩靴直径；

结合图3对s1进行解释说明；

1)、依据工程地质，建立符合实际的数值分析模型；具体如下：

1.1)本实例中钻井船桩靴直径为8m，邻近受影响桩为箱形截面，截面长、宽分别为2m、1.25m，壁厚0.098m，桩长为36m；土层水平向为34m；土层深32.5m，土层为砂土层；桩靴初始位于土层水平方向中心位置；受影响群桩边缘距离桩靴边缘分别为2m、5.125m，其中群桩组两单桩中心距离为3.125m；桩靴最大插桩深度为4m；

1.2)本实例利用软件abaqus建立cel有限元模型；考虑到土层厚度为32.5m，计算模型总高度取36.5m，长度取34m，宽度取12.5m；从计算模型顶面至顶面以下4m为空穴单元，该单元初始没有任何材料属性，以模拟桩靴插入土层时桩靴周围隆起的土体进入空穴单元的情况，从顶面以下4～36.5m为土体单元，考虑到插桩深度为4m，将计算模型中4～11m的土体设为欧拉体，11～36.5m的土体设为拉格朗日体；

1.3)桩靴及受影响桩均采用拉格朗日实体建模，桩靴贯入土层过程中的自身变形很小，故将其设为刚体，桩靴初始位置位于空穴底部的泥面中心处，在距离桩靴边缘2m处设置一群桩组，群桩组中两桩中心距离为3.125m；

1.4)由于cel有限元数值分析采用显式算法，因此，为了提高计算效率，将受影响桩设置为实心桩体，对于受影响桩，桩的材料取为弹性，并按广义虎克定律确定应力-应变关系，计算中，按钢管桩与实体桩抗弯刚度相等的原则，由确定计算模型中实心桩的等效弹性模量，其中e1、e2分别为钢管桩和等效实体桩的弹性模量，i1、i2为钢管桩和等效实体桩的惯性矩；结果为32.6gpa，泊松比取0.3；对于砂土层，理想弹塑性模型中的弹性模量取4.5倍的cpt试验测得的锥尖阻力(qc＝8.89mpa),即40mpa，泊松比取0.3；内摩擦角取30°，黏聚力取0；

1.5)按广义接触条件描述欧拉体与拉格朗日体之间的接触，这种接触算法基于加强的侵入边界方法，采用罚函数和有限滑动接触离散方式，通过指定接触面中的主面与从面，自动计算和跟踪拉格朗日体和欧拉体之间的分界面，使拉格朗日体可以占据欧拉体中的空穴单元，是cel有限元数值分析模型中常用的接触条件，受影响桩与欧拉土体之间的接触以及桩靴与拉格朗日土体之间的接触运均采用这种接触；对于受影响桩和拉格朗日土体之间的接触，采用拉格朗日算法中的面-面接触，且选择刚度大的桩体表面作为主面，土体表面作为从面，桩-土间的摩擦系数取为0.4；

1.6)计算模型的边界条件设置如下：①将欧拉土体及空穴3个方向边界节点的水平速度设置为0，并且设置为欧拉吸收边界，以防止欧拉土体穿越计算模型的边界，消除由于计算模型尺寸选取可能导致的边界效应；②将拉格朗日土体区域的侧面和底面位移设置为0；③通过控制桩靴的竖向位移，使得桩靴以0.5m/s的速率插入到泥面以下4m的位置；

步骤s2中，所述在群桩桩头施加固定的水平荷载，然后模拟插桩过程，输出插桩前后桩体和土体的水平位移云图及桩体的水平位移曲线，具体包括以下几个步骤：

s21、在模型中群桩桩头建立参考点并与桩头顶面耦合，若桩头有承台时，也可在承台表面建立参考点并于表面耦合；

s22、根据桩头或承台受到的水平力大小，在参考点上加上相应的水平荷载；

s23、在桩靴底面建立参考点，并把桩靴约束为刚体；

s24、建立分析步，分别为重力分析步、加载分析步、插桩分析步，并分别设定各个分析步的时间，其中插桩分析步时间要根据插桩深度及插桩速度确定；

s25、通过在边界条件中修改桩靴在插桩分析步中的竖向位移，控制插桩深度；

s26、建立作业并提交计算，在结果中输出插桩前后的桩体和土体的水平位移云图及桩体的水平位移曲线。

步骤s22及s25中，所述的力和位移的施加应避免瞬间施加，一般要为施加的力和位移建立一个平滑的幅值曲线，让力和位移的施加更加平稳，以消除激振力产生的震颤。

结合图3对s2进行解释说明；

2)、在群桩桩头施加固定的水平荷载，然后模拟插桩过程，输出插桩前后桩体和土体的水平位移云图及桩体的水平位移曲线；具体如下：

2.1)、在模型中群桩桩头建立参考点并与桩头顶面耦合；

2.2)、在参考点上加上相应的水平荷载870000n；

2.3)、在桩靴底面建立参考点，并把桩靴约束为刚体；

2.4)、建立分析步，分别为重力分析步、加载分析步、插桩分析步，并分别设定各个分析步的时间，其中插桩分析步时间要根据插桩深度及插桩速度确定，此案例插桩时间取11.3s；

2.5)、通过在边界条件中修改桩靴在插桩分析步中的竖向位移，控制插桩深度，由于插桩深度是桩靴最大直径顶面处入土深度，而参考点设置在桩尖，故此处竖向位移为5.65；

2.6)、建立作业并提交计算，在结果中输出插桩前后的桩体和土体的水平位移云图及桩体的水平位移曲线；

步骤s3中，所述的根据桩体参数及api规范，建立对应的非线性地基梁模型，由于通过cel有限元方法进行数值分析，可以直接确定受影响桩的桩身位移，为了依据桩身位移，按土层的p-y关系确定相应的桩身弯矩，这里采用非线性地基梁有限元模型计算插桩导致的受影响桩的桩身弯矩，具体包括以下几个步骤：

s31、按实际桩的截面尺寸确定梁单元的截面特征参数；

s32、沿桩长设置若干非线性p-y弹簧，模拟桩-土之间相互作用；

s33、参照api规范，确定土层的p-y数据，并输入到地基梁中；

s34、根据桩头的约束情况，在地基梁中施加相应的约束。

步骤s33中，所述的参照api规范，确定土层的p-y数据，对于软粘土，硬粘土和砂性土中桩的p-y曲线api规范分别给出了相应的计算公式：

(1)砂性土的p-y曲线

砂性土分为浅层破坏形式和深层破坏形式，浅层破坏和深层破坏的极限土抗力转折点深度xr按计算，式中：c1、c2、c3为系数，d为外径；如图1所示；φ为沙性土内摩擦角；

当x＜xr时，是浅层破坏；当x＞xr时，是深层破坏；

砂性土的桩侧的极限土抗力随深度不同而变化，浅层破坏和深层破坏极限土抗力为p，p取pux＝(c1x+c2d)γx和pud＝c3dγx两式中的较小值，式中：γ为土体有效重度；x为土表面下计算深度；

给定深度为x的砂性土p-y曲线用表示，式中：pu为深度x处的极限土抗力，a为考虑循环荷载或者静荷载状态的系数，k为地基反力系数，如图2所示，循环荷载：a＝0.9，短期静载：

(2)软粘土中的p-y曲线

软粘土中的p-y曲线法一般假定pu沿深度分布，采用下列两则公式计算，

x＜xr，

x＞xr，pu＝9su；

式中：b为桩侧土计算宽度，这里一般采用桩径，γ为土的有效重度，su为粘土不排水抗剪强度，j为常数，一般粘土取0.5，稍硬粘土取0.25，x为从地面算起的任一深度；

联列以上两列公式，求得地表下土的强度降低的塑性变形深度xr：

经验关系深度xr不可应用于土的强度变化不稳定的情况下，当土体重度和抗剪强度随深度变化时，xr的最小值约为桩径的2.5倍。

结合图3对s3进行解释说明；

3)、根据桩体参数及api规范，建立对应的非线性地基梁模型；具体包括：

3.1)、按实际桩的截面尺寸确定梁单元的截面特征参数，桩截面长、宽分别为2m、1.25m，壁厚0.098m，桩长为36m，弹性模量为71gpa，泊松比为0.3；

3.2)、沿桩长设置若干非线性p-y弹簧，模拟桩-土之间相互作用；

3.3)、参照api规范，确定土层的p-y数据，并输入到地基梁中；

p-y数据计算如下表：

结合图3对s4进行解释说明；

4)、把第二步得到的桩身位移输入到非线性地基梁模型中，输出插桩前后弯矩沿桩身的分布曲线，具体包括：

4.1)、通过边界条件管理器，在桩的两端结点及弹簧结点上输入插桩前相应的水平位移；

4.2)、建立作业并提交计算。输出插桩前弯矩沿桩身的分布曲线；

4.3)、通过边界条件管理器，在桩的两端结点及弹簧结点上输入插桩后相应的水平位移；

4.4)、建立作业并提交计算，输出插桩后弯矩沿桩身的分布曲线；

结合图3对s5进行解释说明；

5)、对比插桩前后的模拟结果，对水平力作用下钻井船插桩对邻近群桩影响程度进行评价；

由图12可以看出前桩插桩前桩头位移为0.285m，插桩后桩头位移为0.615m，则插桩后前桩的桩头位移增加了0.33m，约为该桩插桩前桩头位移的1.16倍；后桩插桩前桩头位移为0.253m，插桩后桩头位移为0.717m，则插桩后后桩的桩头位移增加了0.464m，约为该桩插桩前桩头位移的1.83倍，所以，插桩过程对后桩(即靠近桩靴侧的桩体)影响较大；

由图13可以看出插桩前后，前桩与后桩的最大弯矩位置一致，插桩前、后最大弯矩位置分别为距离土面10.5m、16.5m，即插桩后最大弯矩位置下移了6.5m，约为插桩深度的1.6倍，前桩插桩前、后最大弯矩分别为16.5mn*m、21.2mn*m，则该桩插桩后最大弯矩约为插桩前的1.28倍：后桩插桩前、后最大弯矩分别为16.3mn*m、24mn*m，则该桩插桩后最大弯矩约为插桩前的1.47倍。所以，插桩过程对后桩(即靠近桩靴侧的桩体)影响较大；

综上，本发明能直观合理的评价在水平力作用下钻井船插桩对邻近群桩的影响程度，并且操作简单，易于实际应用。

尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。