[发明专利]一种管土动态耦合作用分析方法

权利要求书

1.一种管土动态耦合作用分析方法，采用了同时考虑土体刚度作用和土体阻尼作用的管土动态耦合作用分析模型：

a)土体刚度作用：

f_k＝k_soil·y

b)土体阻尼作用：

fd=csoil·y·]]>

式中：f_k—土体对管道结构的刚度作用力；

k_soil—土体支撑刚度；

y—管道结构的位移；

f_c—土体对管道结构的阻尼作用力；

c_k—土体阻尼系数；

—管道结构的速度；

根据管道结构嵌入土体的最大深度确定土体非线性迟滞环，其中包括管道结构回弹阶段、管道结构局部分离阶段和管道结构再次回压阶段；

(1)管道结构回弹阶段：

P1=a(y1D)b(Su0+Sugy1)D]]>

P1-2=P1+y-y11k0-y-y1(1+υ)P1]]>

P₂＝-φP₁

y2=y1-(1-υ)P1k01+φυ-φ]]>

式中：y₁—管道结构嵌入土体的最大深度；

P₁—最大土体作用力；

a—与管沟形状和管道结构表面粗糙度有关的参数，由试验获得；

b—与管沟形状和管道结构表面粗糙度有关的参数，由试验获得；

S_u0—海底泥面土体不排水抗剪强度；

S_ug—海底土体不排水抗剪强度垂向增量；

D—管道结构外径；

P_1-2—管道结构回弹阶段土体作用力；

υ—土体特性参数，由试验获得；

k₀—土体特性参数，由试验获得；

P₂—土体最大吸附力；

φ—土体特性参数，由试验获得；

y₂—土体最大吸附力出现位置；

(2)管道结构局部分离阶段：

P2-3=P22+P24[3(2y-(y2+y3)(y2-y3))-(2y-(y2+y3)(y2-y3))3]]]>

式中：P_2-3—管道结构局部分离阶段土体作用力；

y₃—土体作用力消失点位置；

—土体特性参数，由试验获得；

(3)管道结构再次回压阶段：

P3-1=P12+P14[3(2y-(y1+y3)(y1-y3))-(2y-(y1+y3)(y1-y3))3]]]>

式中：P_3-1—管道结构再次回压时土体作用力；

由于土体非线性迟滞环与管道结构嵌入土体的最大深度有关，根据管道结构的运动情况，对管道结构嵌入土体的最大深度进行实时更新；

又由于海床土体对管道结构的支撑刚度与管道结构所处的位置和运动方向有关，故：

当管道结构上升远离土体，且管道结构位置处于y₁和y₂之间时，土体支撑刚度为：

ksoil=k12=(υ-φ)k01+υ;]]>

当管道结构上升远离土体，且管道结构位置处于y₂和y₃之间时，土体支撑刚度为：

ksoil=k23=-φk0(υ-φ)ψ(1+υ)(1+φ);]]>

当管道结构下降压缩土体，且管道结构位置处于y₁和y₃之间时，土体支撑刚度为：

ksoil=k31=k0(υ-φ)(1+ψ)(1+υ)(1+φ);]]>

其它状态下，土体支撑刚度均为零；

土体阻尼系数可以根据下式确定：

csoil=ksoilωEDissipated2πEElastic;]]>

EDissipated=|∫y3y2P2-3dy|+∫y3y1P3-1dy-∫y2y1P1-2dy]]>

EElastic=|k12(y1-y2)2|+|k23(y2-y3)2|+|k31(y1-y3)2|]]>

式中：c_soil—土体阻尼系数；

k_soil—线性化的土体支撑刚度；

ω—响应模态的圆频率；

E_Dissipated—土体吸附的弹性势能；

E_Elastic—土体阻尼消耗的能量。