波浪与地震荷载共同作用下桩的动力响应
土木工程系, 上海 200240
Dynamic Response of Pile at Waterwave Load and Seismic Load
Department of Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
通讯作者: 周香莲,女,副教授;E-mail:zhouxl@sjtu.edu.cn.
责任编辑: 陈晓燕
收稿日期: 2020-05-21
基金资助: |
|
Received: 2020-05-21
作者简介 About authors
刘晨晨(1995-),男,安徽省合肥市人,硕士生,从事海洋土力学的研究
海洋风机单桩基础会受到波浪荷载的循环作用,同时会面临地震荷载的威胁,因此,海洋单桩基础所处的环境较为复杂.然而,目前大多理论研究主要考虑波浪荷载或地震荷载作用下的桩基动力响应问题.本文采用有限元方法建立了波浪-地震-桩-海床耦合模型,数值分析方法基于Abaqus隐式动力分析.采用Morison方程模拟波浪对于桩基础的作用,海床采用Mohr-Coulomb弹塑性模型,桩为弹性介质,地震波通过加速度的形式添加在模型底部.主要研究埋置在海床中桩的动力响应问题,如:加速度、位移、弯矩以及剪力等.结果表明:地震荷载对于海上风机单桩基础具有重要影响.在地震荷载作用下,桩的加速度以及位移响应均有一定程度的放大.土体性质以及桩的参数对于海上风机单桩基础设计至关重要.
关键词:
The pile foundation of offshore wind turbines is not only subjected to the cyclic action of wave load, but also the threat of seismic load. Therefore, the environment of pile in the ocean is complex. However, most theoretical investigations have often focused on the seabed response at water wave load or seismic load respectively. In this paper, a wave-seismic-pile-seabed coupling model is constructed by using the finite element method. The numerical analysis is based on the implicit dynamic analysis in Abaqus. Morison’s equations are used to simulate the effect of water wave on pile foundation. The Mohr-Coulomb model is adopted to simulate the seabed while the pile is considered as an elastic medium. The earthquake is applied on the bottom of the model as acceleration. The dynamic response of the pile embedded in seabed are studied, such as acceleration, displacement, shear force, and bending moment. The results show that seismic load has an important influence on the single pile foundation of offshore wind turbine. Under the action of earthquake load, the acceleration and displacement response of the pile are amplified to a certain extent. The properties of the soil and the parameters of the pile are crucial to the design of the single pile foundation for offshore wind turbines.
Keywords:
本文引用格式
刘晨晨, 张琪, 李明广, 周香莲, 黎蔚杰.
LIU Chenchen, ZHANG Qi, LI Mingguang, ZHOU Xianglian, LI Weijie.
随着传统能源的大量消耗,越来越多的国家开始可再生能源的发展与使用.风能是一种清洁高效的能源,在沿海地区的储备极其丰富.基于成本与技术的考虑,大部分海上风机基础采用的是单桩基础.对于我国而言,海岸线绵长,海上风能资源十分丰富,但大多处于地震烈度较高的地区,海上风机不仅会受到波浪荷载的循环作用,同时也面临着地震破坏的风险.因此,在海上风机单桩基础设计过程中,为确保结构安全性,需充分考虑波浪荷载与地震荷载等复杂的海洋环境荷载对于桩基的影响.
对于海上风机单桩基础受波浪荷载作用时的动力响应及其与海床相互作用的问题,学者们进行了大量研究.其中,Lu等[4]提出了一种耦合模型,用于研究波浪荷载作用下多孔海床和海上单桩的动力响应问题,他们将波浪力考虑为一种外荷载,并在线性波浪理论的背景下,通过波函数展开法进行计算.Sui等[5]采用完全非线性Boussinesq方程模拟了波浪与桩的相互作用,发现单桩基础对于波浪运动以及孔隙水压力响应具有重要影响.李琪等[6]基于COMSOL有限元软件计算分析台风环境中风机结构承受的风荷载、波浪荷载和流荷载,建立典型大直径单桩与导管架基础风机结构简化模型,分析了台风极端工况下风机桩基泥面和轮毂点的转角与位移动态响应.
本文基于Abaqus有限元软件建立了一个三维单桩基础有限元模型,分别考虑了地震荷载单独作用下以及地震荷载与波浪荷载共同作用下单桩的动力响应问题.其中地震荷载以Kobe地震波水平加速度时程曲线作为地震输入,作用在桩上的波浪力采用Morison方程进行简化模拟.并进一步分析讨论了桩的弹性模量、埋深以及土体弹性模量对于单桩动力响应的影响.
1 数值方法
1.1 隐式动力分析
为研究波浪荷载与地震荷载共同作用下桩基动力响应问题,本文采用隐式动力分析方法.隐式积分方法基本假设如下:
式中:u为位移向量,
将式(3)、(4)代入t+Δt时刻平衡方程,则有:
式中:
1.2 波浪理论
本文中波浪运动采用线性波浪理论进行描述,其表达式如下:
式中: ϕ为速度势;H为波高;k为波数,k=2π/Lw,Lw 为波长;T为波浪周期;z竖向坐标;d为水深;x为水平坐标;ω为波浪圆频率,ω=2π/T.
波面方程表达式如下:
式中:η为波面高度.根据式(8),波浪的水平方向速度uw 和加速度
1.3 Morison方程
Morison等[7]提出Morison方程,用于计算波浪荷载作用在结构物上产生的力.目前,与波长相比尺度较小的细长柱体波浪力计算在工程中仍采用Morison方程.一般认为,桩直径与波长的比值,即D/Lw<0.2(D为桩身直径)的柱体为小尺度结构物,适用于Morison方程.根据Morison方程理论,作用在结构物的波浪荷载包括拖曳力和惯性力两个部分.拖曳力及惯性力在时域中表示为
式中:f为单位长度上作用的波浪力;fD、fI分别为单位长度上作用的拖曳力和惯性力;ρ为海水密度;CD为拖曳力系数;CM为惯性力系数.
将式(9)、(10)代入式(12)并积分可得到水平波浪力表达式如下:
式中:g为重力加速度,
本文中,波浪力作为集中力施加在桩上.
2 计算模型
2.1 模型设置
本文采用通用有限元软件Abaqus建立桩-土模型,并对地震荷载与波浪荷载共同作用下桩的动力响应规律进行研究.为更加真实地模拟地震荷载与波浪荷载作用下桩-土动力响应,本文假定土体为Mohr-Coulomb模型,桩为弹性材料.
桩-土之间接触面采用主从接触,刚度较大的桩设为主面,刚度较小的土设为从面.法向接触采用硬接触,切向接触采用库伦摩擦模型,具体表述如下:
式中:τc为接触面发生滑移时的临界切应力;μ为摩擦因数;P为接触面之间的接触压力.
在桩-土动力相互作用中,土体与桩基础中均存在阻尼.本文仅考虑土体的阻尼,忽略桩的阻尼.土体材料阻尼采用瑞利阻尼理论:
式中:α及 β均为阻尼系数,
ω1、ω2分别为土体一阶、二阶自振频率,ξ为阻尼比,取值5%.
2.2 边界条件
本文中,竖向边界允许产生水平位移,约束竖向位移,并对竖向边界采用图1所示的约束压力,以模拟土体真实环境,且在地震过程中,约束压力假定不变.地震波通过水平加速度的形式添加在模型底部x方向,约束压力计算式为
式中:K0为静止土压力系数;ρs为土体密度;zs为土体深度;σn为摩擦角.
图1
3 模型验证
表 1 Peiris模型参数
Tab. 1
土层 | 厚度/ m | 密度/ (kg·m-3) | 弹性模量/ MPa | 泊松比 | 摩擦角/ (°) | 黏聚力/ kPa |
---|---|---|---|---|---|---|
1 | 16 | 1631 | 10 | 0.4 | 0 | 39 |
2 | 6 | 1835 | 15 | 0.4 | 0 | 59 |
3 | 2 | 1886 | 21 | 0.4 | 0 | 83 |
4 | 2 | 1937 | 63 | 0.3 | 35 | 0 |
5 | 7 | 1937 | 248 | 0.3 | 50 | 0 |
图2
4 数值结果和讨论
表2 模型参数
Tab.2
物理量 | 数值 | |
---|---|---|
海水密度/(kg·m-3) | 1000 | |
波高/m | 1.5 | |
波浪周期/s | 5 | |
水深/m | 10 | |
海床厚度/m | 30 | |
海床尺寸/(m×m) | 50×50 | |
土体密度/(kg·m-3) | 1850 | |
土体泊松比 | 0.35 | |
土体弹性模量/MPa | 40 | |
土体摩擦角/(°) | 30 | |
桩长/m | 30 | |
埋深/m | 20 | |
桩直径/m | 2 | |
桩密度/(kg·m-3) | 2300 | |
桩泊松比 | 0.15 | |
弹性模量/GPa | 20 |
图3
4.1 不同荷载组合时桩顶动力响应
图4
图4
不同荷载条件下桩顶动力响应
Fig.4
Dynamic response of head of pile under different load conditions
4.2 桩的弹性模量
基于表格2中模型参数,本节考虑Ep=2,20,100 GPa时桩的动力响应.图5表示具有不同弹性模量的桩受到波浪荷载与地震荷载共同作用时桩顶动力响应.桩顶加速度动力响应受桩的弹性模量影响较大,弹性模量越大,加速度越小.桩顶水平位移受弹性模量影响更为显著,当Ep=2 GPa时,在地震与波浪共同作用下,桩顶水平位移明显大于Ep为20和100 GPa时的桩顶水平位移,且Ep=20,100 GPa时桩顶水平位移相差不大.值得注意的是,在荷载作用在桩的初期,尤其是Ep=2 GPa时,波浪荷载起主要作用.
图5
图5
不同弹性模量时桩的桩顶动力响应
Fig.5
Dynamic response of head of pile with different elastic modulus of pile
为了更加直观地表示桩的弹性模量对于动力响应的影响,图6绘制出了t=9.06 s时的桩身相对加速度变化以及t=6.48 s时的桩身相对位移变化.图中:h为桩身高度,桩顶为原点,向下为正;ar和ur分别为桩身相对水平加速度和相对水平位移.可以看出,埋置部分的桩受桩的弹性模量影响较小,而露出部分影响较大.
图6
图6
不同弹性模量时桩沿桩身动力响应
Fig.6
Dynamic response along pile with different elastic modulus of pile
4.3 桩的埋深
图7
图7
不同埋深时的桩顶动力响应
Fig.7
Dynamic response of head of pile at different buried depths
图8为加速度响应与水平位移响应沿桩身的变化情况.可以看出,埋置部分的桩的加速度与位移变化规律同露出部分有显著差别,埋置部分的桩的相对加速度随着埋深的增加而增加,而露出部分则完全相反.对于相对位移,埋置部分的桩的相对位移也随着埋深增加而增加,露出部分随着埋深增加而减小.由此表明,海床土体对于桩会有约束效应,能够在地震发生时对于桩产生一定的保护作用.埋深越大时,这种保护作用越大.
图8
4.4 土的弹性模量
图9
图9
土的弹性模量不同时的桩顶动力响应
Fig.9
Dynamic response of head of pile with different elastic modulus of soil
图10
图10
土的弹性模量不同时沿桩身动力响应
Fig.10
Dynamic response along pile with different elastic modulus of soil
5 结论
(1) 相较于地震荷载对于桩的影响,波浪荷载影响相对较小.在地震荷载作用下,桩身对于地震加速度具有一定的放大效应.
(2) 在其他条件相同的情况下,随着桩的弹性模量的增加,加速度响应与位移响应均随着弹性模量的增加而减小.埋置部分的桩受桩弹性模量的影响较小,而露出部分影响较大.
(3) 桩顶加速度随着土体弹性模量的增加而减小.土体弹性模量越大,桩的水平位移越小.土体对于桩的约束作用随着土体弹性模量增加而增强.
参考文献
桩基础在阻尼与分层弹性地基场地土波动影响下的横向地震反应分析
[J]. ,
The lateral seismic response analysis for the pile foundation under the influence of damping and layer elasto-foundation bed site soil wave motion
[J].
Vertical earthquake response of megawatt-sized wind turbine with soil-structure interaction effects
[J]. ,
水平环境荷载与地震动联合作用下海上风机单桩基础动力响应模型试验
[J]. ,
Dynamic model test on monopile for offshore wind turbine under jointed lateral environmental load and seismic load
[J].
Dynamic response of a porous seabed and an offshore pile to linear water waves
[C]// ,
Three-dimensional numerical model for wave-induced seabed response around mono-pile
[J]. ,DOI:10.1080/17445302.2015.1051312 URL [本文引用: 1]
台风环境中典型海上风机结构的动力响应数值分析
[J]. ,
Numerical analysis on dynamic response of typical structure of offshore wind turbines under typhoon
[J].
The force exerted by surface waves on piles
[J]. ,DOI:10.2118/950149-G URL [本文引用: 1]
Soil-pile interaction of pile embedded in deep-layered marine sediment under seismic excitation
[J]. ,DOI:10.2749/101686614X13854694314720 URL [本文引用: 1]
/
〈 | 〉 |