一种群桩动力响应和动力阻抗的分析方法和系统

摘要

本发明公开了一种群桩动力响应和动力阻抗的分析方法和系统，本发明考虑波浪荷载对群桩的影响进行了群桩动力稳定问题的研究，采用改进的Vlasov地基模型来计算地基反力，利用相互作用因子法和矩阵传递法相结合的方式建立主动桩和被动桩的动力稳定方程，得到了邻桩之间的动力相互作用因子和群桩阻抗，并对群桩的稳定性进行了参数分析。通过研究发现，波浪荷载的存在使得群桩的动力响应明显增大；群桩的动力阻抗和相互作用因子主要受土体的参数影响，但波浪荷载的存在会影响土体的一些参数；增大表层土体的弹性模量可以有效提高群桩阻抗；本发明中改进的Vlasov地基模型能够准确符合工程实际。

说明书

技术领域

本发明属于岩土工程技术领域，具体涉及一种基于相互作用因子叠加法计算群桩动力响应和动力阻抗的计算方法。

背景技术

目前，对于单桩在简谐荷载作用和冲击荷载作用下的动力稳定性问题的研究，在实际的工程中大多数的桩基都是以群桩的形式而出现的，常见的有2x2，3x3，6x6群桩等，在一些大型的结构中桩数可能会更多。相比于单桩来说，群桩的动力分析要复杂许多，不同于单桩仅需要考虑自身所承受的荷载作用，群桩中的桩基还要考虑其他相邻的桩基对其自身的影响作用，即群桩效应。考虑群桩效应就要涉及到桩-土-桩之间的相互作用问题，而在本发明中主要研究的是动力问题，因此主要是桩土之间的动力相互作用，分析土层中的桩土动力相互作用是进一步研究桩基动力响应的基础。

发明内容

为了解决现有技术问题，本发明的目的在于克服已有技术存在的不足，提供一种群桩动力响应和动力阻抗的分析方法和系统，基于相互作用因子叠加法计算群桩动力响应和动力阻抗的计算方法。相互作用因子叠加法，推导过程简单，计算量小，是目前来说一种比较适合的计算群桩动力响应和动力阻抗的方法。之前的一些研究中普遍采用Winkler地基模型来模拟土体的抗力，但由于Winkler模型较为简单，简化太多，缺点也比较明显，因此在本发明的计算中考虑了土体之间的连续特性，采用了改进的Vlasov地基模型来计算土体的地基反力；以解决了上述背景技术中提出的问题。

为达到上述发明创造目的，本发明采用如下技术方案：

一种群桩动力响应和动力阻抗的分析方法，采用改进的Vlasov地基模型来计算地基反力；利用相互作用因子法和矩阵传递法相结合的方式，建立主动桩和被动桩的动力稳定方程，得到了邻桩之间的动力相互作用因子和群桩阻抗，并对群桩的稳定性进行了参数分析，得到群桩动力响应和动力阻抗。

优选地，本发明群桩动力响应和动力阻抗的分析方法，包括如下步骤：

(1)参数选取

桩-土-桩之间的动力相互作用是分析群桩动力响应的重要部分，通过对群桩之间动力相互作用的分析，得到主动桩-土-被动桩之间的关系，继续分析群桩的动力响应，而动力相互作用的分析首先要从主动桩开始；主动桩动力分析模型如下：

设定，N

其中：

g为重力加速度，取9.8m/s

根据模型，得到土层的运动平衡方程为：

其中，k

根据群桩涉及到桩-土-桩之间的动力相互作用，描述桩-土之间的相互作用，并基于连续介质模型推导而来的VLasov地基模型来模拟土体的抗力，具体的计算公式如下所示：

q(x)＝k

式中

h(z)为竖向位移的衰减函数该函数的参数，采用Vallabhan和Das，通过使用另一个新的参数γ把位移函数和衰减函数相互联系起来，得到位移函数和衰减函数的准确表达式，称为改进的Vlasov地基模型；将采用改进的Vlasov地基模型来进行地基反力的计算，根据Vallabhan和Das，基于桩基横向位移下的地基模型参数为：

式中，η为lamé常数，

G为土体的剪切模量，

γ为衰减参数通过迭代法计算，

K

K

其中

r为柱坐标中的变量；

地基土反力q(x)公式为：

q(x)＝k

土体的阻尼c

其中ρ

(2)模型方程的建立

通过桩身运动平衡方程得到桩身的稳态振动方程的一般形式为：

考虑桩基部分嵌固到土体中，桩身处于水中的部分要承受波浪荷载的作用，没有土体的约束作用，将其桩身分为两部分，桩身处于土体中的部分振动方程如上式所示，而桩身露出土体的部分振动方程则如下式所示：

将桩身的位移U

桩身深入土体部分：

桩身处于水中部分：

式中：

其中：

h

然后通过求解上述的高阶振动微分方程式可以得到如下形式的通解：

其中：

A

上式的通解为：

U′

其中：

A′

(3)分析桩身露出土体的部分，即分析承受波浪荷载的桩身部分：

对于桩身露出土体的部分，将其视为一个单元层，类似于土层的划分，将其视为一层，对于该部分的截面转角

有如下的关系式：

将其整理为矩阵的形式如下式所示：

取桩顶的z＝0，得到如下所示的式子：

然后桩身处于水中的部分与土层的交界处，令z＝h

通过对矩阵的变换后，将桩身露出土体部分的桩顶的位移和水土交界处的位移相联系起来，如下所示：

假设露出土体部分的桩长为L

桩身处于土体中的部分的研究：根据桩身处于土体中的部分涉及到了土体的约束作用和土体的分层问题，具体计算步骤如下：

桩身处于土体部分的位移U

此时桩顶部的位移就变为了水-土交界面处的位移，桩底位移即为实际的桩底位移；土层单元内的剪力，弯矩与桩身水平位移的关系如下：

将上述公式整理为矩阵形式如下：

令：

假设此时的桩顶处即土体的表面处z＝0，可以得到：

同样的在桩基下部取z＝h

如果将土体分为多层，根据土体的连续性原则，则得到，u

使用传递矩阵法来将土层之间的位移，剪力，转角和弯矩通过参数传递矩阵相互连接起来，如下式所示：

式中L

令：

则上式可以表示为如下的形式：

假设桩底的边界条件为固定端，桩顶为自由端，则：

整理上式可以得到：

[K

整理上式可以得到：

U

最后计算桩顶总位移和总的转角时，将处于土体部分的桩顶位移U

上式即为结合两部分的桩身动力响应后得到的最终桩顶的位移，转角，剪力和弯矩；

根据单桩水平阻抗的定义，得到如下所示的单桩阻抗计算公式；

其中，阻抗R

(4)群桩模型建立：

4-1)群桩的模型分析：

设定ψ(s,θ)为土体应力波的衰减函数，f

式中，

此处s为桩间距，θ为桩-桩之间的夹角；V

式中，V

主动桩振动引起的应力波传出时的位移为U

U

假设被动桩的位移为U

桩身处于水中部分的振动平衡方程：

桩身处于土体部分的振动平衡方程：

被动桩相比于主动桩波浪荷载f

上式的计算过程如下所示：

首先令

则上式表示为如下的形式：

式中，

上式的通解为：

式中,

在土层单元内，每个桩基截面的截面转角

式中，

其中

根据传递矩阵将每一土层的位移，转角，剪力，弯矩联系起来，如下式所示，整理后的传递关系矩阵为：

式中，

可以将上式表示为：

根据模型假设边界条件为桩顶固定，所以：

然后将边界条件代入到上式中得到：

其中

[μ

根据相互作用因子的定义，得到：

群桩的水平相互作用因子为：

群桩的摇摆相互作用因子：

群桩的桩顶总位移，转角参数与单桩计算方法相同具体为：

4-2)群桩阻抗分析：

群桩水平阻抗的计算具体计算如下所示，假设群桩桩数为n，群桩的水平位移u

设主动桩j对被动桩i的影响因子为χ

式中，R

群桩的水平动力阻抗为：

K

一种群桩动力响应和动力阻抗的分析系统，包括：

存储子系统，由于存储计算机程序；

信息处理子系统：用于执行计算机程序时实现本发明群桩动力响应和动力阻抗的分析方法的步骤。

本发明与现有技术相比较，具有如下显而易见的突出实质性特点和显著优点：

1.本发明利用相互作用因子法和矩阵传递法相结合的方式建立主动桩和被动桩的动力稳定方程，得到了邻桩之间的动力相互作用因子和群桩阻抗，并对群桩的稳定性进行了参数分析，本发明中改进的Vlasov地基模型能够准确符合工程实际；

2.本发明采用相互作用因子叠加法，推导过程简单，计算量小，是目前来说一种比较适合的计算群桩动力响应和动力阻抗的方法；在本发明的计算中考虑了土体之间的连续特性，采用了改进的Vlasov地基模型来计算土体的地基反力，计算开销显著降低；

3.本发明方法能降低成本，能够准确符合工程实际，适合推广使用。

附图说明

图1为本发明中主动桩A模型图。

图2为本发明中所述Vlasov地基模型桩土相互作用模型图。

图3为本发明中迭代计算流程图。

图4为本发明中群桩模型图。

图5为本发明中主动桩和被动桩位置图。

图6为本发明中2x2群桩示意图。

图7为本发明中群桩相互作用因子-实部。

图8为本发明中群桩相互作用因子-虚部。

图9为本发明中群桩阻抗实部刚度。

图10为本发明中群桩阻抗虚部刚度。

图11为本发明中群桩阻抗实部刚度随土层弹性模量比的变化图。

图12为本发明中群桩阻抗虚部随土层弹性模量比变化图。

图13为本发明中不同地基模型群桩阻抗实部变化图。

图14为本发明中不同地基模型群桩阻抗虚部变化图。

图15为本发明中水平动力相互因子随a

图16为本发明中水平动力相互因子随a

图17为本发明中波高对群桩位移响应u的影响图。

图18为本发明中波长对群桩位移响应u的影响图。

具体实施方式

以下结合具体的实施例子对上述方案做进一步说明，本发明的优选实施例详述如下：

实施例一：

在本实施例中，参见图1-2，一种群桩动力响应和动力阻抗的分析方法，采用改进的Vlasov地基模型来计算地基反力；利用相互作用因子法和矩阵传递法相结合的方式，建立主动桩和被动桩的动力稳定方程，得到了邻桩之间的动力相互作用因子和群桩阻抗，并对群桩的稳定性进行了参数分析，得到群桩动力响应和动力阻抗。

本实施例对群桩的稳定性进行了参数分析，本发明中改进的Vlasov地基模型能够准确符合工程实际。

实施例二：

本实施例与实施例一基本相同，特别之处在于：

在本实施例中，一种群桩动力响应和动力阻抗的分析方法，包括如下步骤：

(1)参数选取

桩-土-桩之间的动力相互作用是分析群桩动力响应的重要部分，通过对群桩之间动力相互作用的分析，得到主动桩-土-被动桩之间的关系，继续分析群桩的动力响应，而动力相互作用的分析首先要从主动桩开始；如下图1所示，图1所示为主动桩A的模型示意图；

主动桩动力分析模型如下：

设定，N

其中：

g为重力加速度，取9.8m/s

根据模型，得到土层的运动平衡方程为：

其中，k

根据群桩涉及到桩-土-桩之间的动力相互作用，描述桩-土之间的相互作用，并基于连续介质模型推导而来的VLasov地基模型来模拟土体的抗力，具体的计算公式如下所示：

q(x)＝k

式中

h(z)为竖向位移的衰减函数该函数的参数，采用Vallabhan和Das，通过使用另一个新的参数γ把位移函数和衰减函数相互联系起来，得到位移函数和衰减函数的准确表达式，称为改进的Vlasov地基模型；将采用改进的Vlasov地基模型来进行地基反力的计算，根据Vallabhan和Das，基于桩基横向位移下的地基模型参数为：

式中，η为lamé常数，

G为土体的剪切模量，

γ为衰减参数通过迭代法计算，迭代过程见图3，

K

K

其中

r为柱坐标中的变量；

地基土反力q(x)公式为：

q(x)＝k

土体的阻尼c

其中ρ

(2)模型方程的建立

通过桩身运动平衡方程得到桩身的稳态振动方程的一般形式为：

考虑桩基部分嵌固到土体中，桩身处于水中的部分要承受波浪荷载的作用，没有土体的约束作用，将其桩身分为两部分，桩身处于土体中的部分振动方程如上式所示，而桩身露出土体的部分振动方程则如下式所示：

将桩身的位移U

桩身深入土体部分：

桩身处于水中部分：

式中：

其中：

h

然后通过求解上述的高阶振动微分方程式可以得到如下形式的通解：

其中：

A

上式的通解为：

U′

其中：

A′

(3)分析桩身露出土体的部分，即分析承受波浪荷载的桩身部分：

对于桩身露出土体的部分，将其视为一个单元层，类似于土层的划分，将其视为一层，对于该部分的截面转角

有如下的关系式：

将其整理为矩阵的形式如下式所示：

取桩顶的z＝0，得到如下所示的式子：

然后桩身处于水中的部分与土层的交界处，令z＝h

通过对矩阵的变换后，将桩身露出土体部分的桩顶的位移和水土交界处的位移相联系起来，如下所示：

假设露出土体部分的桩长为L

桩身处于土体中的部分的研究：由于桩身处于土体中的部分涉及到了土体的约束作用，和土体的分层问题，因此相比于桩身处于水中部分的位移，它的计算过程要稍微复杂一些，具体计算步骤如下：

桩身处于土体部分的位移U

此时桩顶部的位移就变为了水-土交界面处的位移，桩底位移即为实际的桩底位移；土层单元内的剪力，弯矩与桩身水平位移的关系如下：

将上述公式整理为矩阵形式如下：

令：

假设此时的桩顶处即土体的表面处z＝0，可以得到：

同样的在桩基下部取z＝h

如果将土体分为多层，根据土体的连续性原则，则得到，u

使用传递矩阵法来将土层之间的位移，剪力，转角和弯矩通过参数传递矩阵相互连接起来，如下式所示：

式中L

令：

则上式可以表示为如下的形式：

假设桩底的边界条件为固定端，桩顶为自由端，则：

整理上式可以得到：

[K

整理上式可以得到：

U

最后计算桩顶总位移和总的转角时，将处于土体部分的桩顶位移U

上式即为结合两部分的桩身动力响应后得到的最终桩顶的位移，转角，剪力和弯矩；

根据单桩水平阻抗的定义，得到如下所示的单桩阻抗计算公式；

其中，阻抗R

(4)群桩模型建立：

4-1)群桩的模型分析如图4所示：

图中，设定ψ(s,θ)为土体应力波的衰减函数，f

式中，

此处s为桩间距，θ为桩-桩之间的夹角；V

式中，V

主动桩振动引起的应力波传出时的位移为U

U

假设被动桩的位移为U

桩身处于水中部分的振动平衡方程：

桩身处于土体部分的振动平衡方程：

被动桩相比于主动桩波浪荷载f

上式的计算过程如下所示：

首先令

则上式表示为如下的形式：

式中，

上式的通解为：

式中,

在土层单元内，每个桩基截面的截面转角

式中，

其中

根据传递矩阵将每一土层的位移，转角，剪力，弯矩联系起来，如下式所示，整理后的传递关系矩阵为：

式中，

可以将上式表示为：

根据模型假设边界条件为桩顶固定，所以：

然后将边界条件代入到上式中得到：

其中

[μ

根据相互作用因子的定义，得到：

群桩的水平相互作用因子为：

群桩的摇摆相互作用因子：

群桩的桩顶总位移，转角参数与单桩计算方法相同具体为：

4-2)群桩阻抗分析：

群桩水平阻抗的计算具体计算如下所示，假设群桩桩数为n，群桩的水平位移u

设主动桩j对被动桩i的影响因子为χ

式中，R

群桩的水平动力阻抗为：

K

本实施例方法实现了计算相互作用群桩动力响应和动力阻抗的解析，目前国内外关于这方面的研究主要为数值方法，数值方法的适用范围比较广泛，但是计算过程复杂，计算量大，在一些复杂的结构中还存在一些困难而且计算太慢；不方便在实际工程设计中的计算分析。基于此本发明考虑波浪荷载对群桩的影响进行了群桩动力稳定问题的研究，采用改进的Vlasov地基模型来计算地基反力，利用相互作用因子法和矩阵传递法相结合的方式建立主动桩和被动桩的动力稳定方程，得到了邻桩之间的动力相互作用因子和群桩阻抗，并对群桩的稳定性进行了参数分析。通过研究发现，波浪荷载的存在使得群桩的动力响应明显增大；群桩的动力阻抗和相互作用因子主要受土体的参数影响，但波浪荷载的存在会影响土体的一些参数；增大表层土体的弹性模量可以有效提高群桩阻抗；本实施例方法中改进的Vlasov地基模型能够准确符合工程实际。

实施例三：

本实施例与实施例一基本相同，特别之处在于：

在本实施例中，本实施例具体涉及一种基于相互作用因子叠加法计算群桩动力响应和动力阻抗的计算方法，如图6所示的2x2群桩来进行算例验证，波浪荷载方向如图所示，桩身的参数为：桩长l＝37.6m，桩身深入土体中的长度取为18.2m,桩径d＝1.0m，横向桩距为5m，纵向桩距为6m，E＝30GPa，ρ

表1：土体参数表

模拟成果图如图7～18，本实施例算例首先分析了桩1与其他三根桩之间的相互作用因子，通过与黄茂松学者等的结果进行了对比验证，两个结果吻合相对较好，满足精度要求。如图7，图8所示，图中文献解即为黄茂松等的解。图7为通过计算所得的群桩相互作用因子的实部部分，图8为计算所得的群桩相互因子的虚部部分，通过对比发现，本实施例的结果略高于黄茂松等的结果，这是由于本实施例中考虑了土体之间的连续性影响所致。从上面两幅图中还可以看到，源桩1对于不同位置处的桩之间的相互作用因子是不相同的，桩1-3之间的相互作用因子和桩1-4之间的相互因子趋势基本相同，桩1-2位置最近，两桩之间的相互作用因子曲线则明显区别于桩1-3和桩1-4之间的相互因子曲线。

本实施例算例分析还进行了群桩之间的动力阻抗研究，将所得的结果进行了无量纲化处理，并与kaynia等的结果进行了对比验证，同时还对比了在理论计算中是否考虑轴力对群桩阻抗的影响，具体如图9，图10所示。从图中看到，当桩间距与桩径之比较小时，随a

参数分析

桩基参数为：桩长L＝55m，桩径d＝1.6m，桩身深入土中部分长度为30m，E＝3×10

假设土层第二层和第三层的弹性模量不变，改变表层土体的模量使其为：E

从两幅图中看到，随着a

之前的一些研究中普遍采用Winkler模型来模拟土体反力，在发明中为了更好的考虑土体的连续特性，从而来深入的研究群桩的桩土特性，所以采用了基于连续介质模型的改进的Vlasov地基模型来计算桩侧的土体反力。针对两种不同的地基模型，下图13分析了Winkler模型和Vlasov地基模型下群桩的的阻抗随a

首先，从上面两幅图中发现，两种模型下实部阻抗都是先随着a

群桩中不同位置的桩，相对应的桩基之间的相互作用因子也不相同，这在实施例验证部分已经得到了证实，接下来的部分将会具体研究桩间距与桩径的比值对水平动力相互作用因子的影响，具体如下图15，图16所示：

从两幅图中看到，随着桩间距和桩径比值的增大，动力相互作用因子曲线的波动变化开始增大，当s/d＝3时，动力相互因子随a

波浪荷载作为一种动力荷载施加在桩身上，对于群桩的动力响应也会产生一定的影响，下面将要分析波浪荷载的波长和波高对于群桩的位移响应影响。如下图17，图18所示：

从图中看出，随着波高H的不断增大，桩身的位移响应呈线性增加，相比于单桩来说增加的幅度要小；波长L对位移响应的影响与波高H有所不同，位移响应u随着波长L的增加，不再是单纯的线性的增加，而是呈非线性的增加，且增加的速度越来越快，与单桩时波长对位移的影响基本相同，由于存在群桩效应的影响，波浪荷载作用下群桩的位移响应随波长增加的幅度要比单桩时小许多。

本实施例考虑波浪荷载对群桩的影响进行了群桩动力稳定问题的研究，采用改进的Vlasov地基模型来计算地基反力，利用相互作用因子法和矩阵传递法相结合的方式建立主动桩和被动桩的动力稳定方程，得到了邻桩之间的动力相互作用因子和群桩阻抗，并对群桩的稳定性进行了参数分析。

实施例四：

本实施例与上述实施例基本相同，特别之处在于：

在本实施例中，一种群桩动力响应和动力阻抗的分析系统，包括：

存储子系统，由于存储计算机程序；

信息处理子系统：用于执行计算机程序时实现所述群桩动力响应和动力阻抗的分析方法的步骤。本实施例通过计算中考虑了土体之间的连续特性，采用了改进的Vlasov地基模型来计算土体的地基反力，计算开销显著降低，解决成本，适合在工程实际中推广应用。

上面对本发明实施例结合附图进行了说明，但本发明不限于上述实施例，还可以根据本发明的发明创造的目的做出多种变化，凡依据本发明技术方案的精神实质和原理下做的改变、修饰、替代、组合或简化，均应为等效的置换方式，只要符合本发明的发明目的，只要不背离本发明的技术原理和发明构思，都属于本发明的保护范围。