Dynamic responses of calcareous foundation reinforced by microbially induced calcite precipitation
-
摘要: 微生物加固技术是一种有效提高土体强度、抑制土体发生液化破坏的绿色环保加固技术。采用温控微生物加固法对南海钙质砂地基模型进行微生物加固处理,并开展了一系列模型振动台试验研究,系统讨论了微生物加固程度和土体深度对微生物加固钙质砂地基的动应力应变关系、剪切模量、剪切波速以及动强度等动力特性的影响。试验研究表明微生物加固对钙质砂地基的动力学特性影响十分显著,具体表现在:随着微生物加固程度的提高,动剪应变显著降低,滞回圈骨干曲线的斜率逐渐增大,滞回圈面积和土体能量的耗散逐渐减小;动剪切模量随着加固程度的提高而增大,但增大幅度逐渐降低;深部土体的剪切模量衰减要高于上部土体;剪切波速值随着微生物加固程度的提高而显著提高,振动后的土体剪切波速值要高于振动前。动三轴试验获得的统一动强度方程,通过折减系数进行修正处理,可在一定程度上模拟振动台试验所获得的微生物加固钙质砂动强度发展规律。Abstract: The microbially induced calcite precipitation (MICP) is a green and ecofriendly technique to efficiently improve soil strength and mitigate liquefaction potential of soil. In this study, the temperature-controlled MICP method is used to reinforce the foundation model made with the calcareous sand from the South China Sea. A series of shaking table tests are performed to investigate the effects of biocementation level and soil depth on the dynamic stress-strain relationship, shear modulus, shear wave velocity and dynamic strength of the MICP-treated calcareous sand foundation. The test results show that the MICP can affect the dynamic response of the calcareous sand foundation significantly, indicating that with the increase of the biotreatment level, the dynamic shear strain decreases significantly; the area of stress-strain hysteresis ring and the energy dissipation decrease; the dynamic shear modulus increases with the decrease of increment amplitude; the shear modulus degradation at a larger depth is higher than that at a lower depth. The shear wave velocity increases with the biocementation level and becomes higher after shaking. The dynamic strength equation for the biotreated soil models from the triaxial cyclic loading tests, multiplying a reduction coefficient, can be used to simulate the dynamic soil strength of the MICP-treated calcareous sand foundation in the shaking table tests.
-
Keywords:
- biotreatment /
- calcareous sand /
- dynamic response /
- dynamic strength /
- shaking table test
-
0. 引言
大力发展可再生能源是加快生态文明建设的重要支撑,其中海上风机得到了快速发展[1]。但中国海底地震频繁,海底地震及其引发的次生灾害会对海工结构造成毁灭性的损害[2]。因此,如何保障海上风机的地震安全性是海上风机发展面临的重要问题。
海上风机地震响应分析是保障其地震安全性的重要环节。其中自由场响应作为海上风机地震响应分析的输入,是研究地震波散射的基础。Thomson[3]最先给出了层状介质中波传播问题的传递矩阵解,但国内外学者借此研究时多考虑场地为干土情形。在构建自由场时,大多直接在底部输入地震动,忽略自由场的作用[4],或只考虑底部、侧面自由场中的一个[5]。杜修力等[6]分析指出,不完整考虑自由场效应将无法得到正确结果。同时直接将海床表面的海域场地地震记录作用在基岩底部计算自由场并不合理。
其次,人工边界的实施可以将无限域波动散射问题转化为有限域问题来解决。目前人工边界主要有透射边界[7]、黏性边界[8]和黏弹性边界[9]等。
在此基础上,考虑海水-海床-风机间的相互作用分析。其中海水-结构属于流固耦合,多简化为Morison方程、辐射波浪理论和附加质量法;海床-结构属于土-结相互作用,多简化为集总参数法、子结构法和整体分析直接法,求解十分不便。冯玉涛等[10]基于Morison方程实现了考虑动水效应的桩-土-结构地震响应分析。魏凯等[11]提出了动水附加质量简化有限元计算方法,并表示Morison方程存在一定的适用范围。Zuo等[12]利用ABAQUS软件研究了在土-结相互作用下风力机在不同运行状态下的动力响应。Haciefendioglu[13]建立了海水-单桩-土耦合的三维数值有限元模型,分析了场地条件对海上风机的影响。Lee等[14]采用非线性p-y、t-z曲线和附加质量法考虑了桩-土作用和海水对风机的影响。Kim等[15]通过建立非线性弹簧单元模拟了桩土相互作用。Francesca等[16]利用简化模型计算了考虑桩土相互作用的风机结构地震响应。Yang等[17]发现土-结相互作用显著影响海上风机塔顶位移和泥面弯矩。Wang等[18]的工作印证了这一结论,指出分析海上风机地震响应时应考虑气动-伺服-水动-结构弹性耦合及土-结相互作用的影响。目前国内外大多只考虑海水-海床、海床-结构和海水-结构三种耦合中的一个或两个,且是单向耦合,并做简化处理,或同时考虑三种耦合做全显式分析或做全隐式分析,计算效率有待提高。
基于此,本文拟将场地采用显式有限元分析,并统一为广义饱和多孔介质实现耦合,结构采用隐式有限元分析,显隐式通过重叠一层单元实现场地-结构每时步耦合,结合人工边界条件,实现海水-饱和海床-风机相互作用高效分区分析,并探讨海水深度、海床波速、入射角度对海上风机地震响应的影响。
1. 基本理论
1.1 问题描述
海水-饱和海床-风机相互作用分析问题本质上属于地震波散射问题,涉及结构与周围半无限域环境的动力耦合作用。采用数值方法分析无限域波动散射问题时,一般将无限域截取为有限计算区域,并通过施加人工边界来考虑无限域的影响,其具体求解过程为:自由场分析(结构不存在时的波动分析),为散射问题分析提供波动输入;人工边界的实施;环境介质-结构相互作用分析(海水-海床-风机耦合),如图 1所示。
1.2 自由场分析
自由场模型一般假设为水平成层半空间模型,该假设对风机结构地震响应影响较小。如图 2(a)所示,当场地为海水-海床-露头基岩时,一般将控制点地震动响应折减一半作为基岩半空间的地震输入,即
I2s(t)=12R1s(t) 。 (1) 式中:R1s为#1点实测的露头基岩地震动的地震响应;I2s指#2点基岩底部的地震输入;t为时间。
但在海域地震安全性评价及海洋工程结构抗震规范中,一般给定海床表面的加速度谱或设计谱,因此需要通过反演得到基岩输入处的地震动(入射波)。本文假设在如图 2(b)所示水平成层半空间模型下,根据已知水-土界面处的控制点A的地震运动,反演得到基岩半空间控制点B输入处的地震激励,即
H(ω)=RAp(ω)IBp(ω) 。 (2) 式中:ω为频率;IBp为基岩底部控制点B处的脉冲波输入;RAp为在基岩底部输入脉冲波后海床表面控制点A处得到的脉冲波响应,该响应可根据水的声波方程、Boit饱和多孔介质运动方程、弹性介质的波动方程,以及各界面边界条件,采用传递矩阵方法得到。两者比值得到由B点至A点的传递函数H,并基于H,由海床表面控制点A处的地震响应RAs反演得到基岩底部控制点B处的地震输入IBs:
IBs(ω)=RAs(ω)H(ω) 。 (3) 同样,若已知水平成层半空间模型任意控制点i处地震动,可反演得到基岩输入处B地震动输入。
得到基岩输入处的入射波后,可得到透射边界区的自由场响应uf和Uf,分别输入到海水-饱和海床-风机系统的左侧、右侧、前侧、后侧和底部。
1.3 海域场地分析
场地分析包含海水、饱和海床和基岩,桩基础含在场地内,与基岩分析相同。该计算区分为人工边界结点区和内部结点区:人工边界结点采用多次透射公式计算;内部结点根据广义饱和多孔介质理论,将基岩(基础)视为孔隙率为0的饱和介质,而海水视为孔隙率为1的饱和介质,海床为孔隙率介于0~1的饱和介质。详细推导见参考文献[19,20]。
(1)内部结点
内部结点分为一般内部点和界面点,一般内部点是指具有同一孔隙率的广义饱和介质内部点,界面点指两种不同孔隙率的广义饱和介质界面处点。
广义饱和多孔介质基本微分方程如下:
固相平衡方程:
LTsσ′−(1−β)LTwP+b(˙U−˙u)=(1−β)ρs¨u 。 (4) 液相平衡方程:
−βLTwP+b(˙u−˙U)=βρw¨U 。 (5) 相容方程(考虑初始孔压和初始体应变为零时):
−βP=Ew[βew+(1−β)es] 。 (6) 式中:Ls,Lw为微分算子矩阵;σ′为有效应力;b=β2/k,k为流体渗透系数,β为孔隙率;P为孔隙水压;u,U分别是固相和液相位移;ρs,ρw分别是固相和液相密度;es,ew分别为固相和液相体应变;Ew为流体体变模量。
依据微分方程(4)、(5)和相容方程(6)以及边界条件,利用伽辽金法建立广义饱和多孔介质的有限元空间离散方程,得到任意内部结点i的固相、液相运动平衡方程:
¨uiMsi+Fsi+Tsi−Ssi=0, (7) ¨UiMwi+Fwi+Twi−Swi=0 。 (8) 式中:ui,Ui分别为结点i的固、液相位移;Msi,Mwi分别为结点i的固、液相质量;Fsi,Fwi分别为结点i的固、液相本构力;Tsi,Twi分别为结点i的固、液相黏性阻力;Ssi,Swi分别为结点i的固、液相界面力。
对方程(7),(8)采用显式中心差分法进行时间积分,得到一般内部结点i的固、液相位移递推公式:
u(p+1)i=2upi−u(p−1)i−(Δt)2msi(Fspi+Tspi−Sspi) , (9) U(p+1)i=2Upi−U(p−1)i−(Δt)2mwi(Fwpi+Twpi−Swpi) 。 (10) 式中:上标p−1,p,p+1分别为结点i的3个相邻时刻;Δt为时间步长。
当结点i在同一饱和介质内部时,根据位移和应力连续条件,方程(9),(10)中的Sspi、Swpi均为0,当在不同饱和介质交界面上时,方程(9),(10)变为
u(p+1)i=2upi−u(p−1)i−(Δt)2msi(Fspi+Tspi−SspNi−SspTi) , (11) U(p+1)i=2Upi−U(p−1)i−(Δt)2mwi(Fwpi+Twpi−SwpNi−SwpTi) , (12) u(p+1)j=2upj−u(p−1)j−(Δt)2msj(Fspj+Tspj−¯SspNj−¯SspTj) , (13) U(p+1)j=2Upj−U(p−1)j−(Δt)2mwj(Fwpj+Twpj−¯SwpNj−¯SwpTj) 。 (14) 其中:方程(11),(12)为结点i的固液相位移递推公式;方程(13),(14)为与结点i在不同饱和介质交界面处形成结点对的结点j的固液相位移递推公式;SspNi,SspTi分别为p时刻作用在结点i固相的界面法向力和界面切向力;SwpNi,SwpTi为p时刻作用在结点i液相的界面法向力和界面切向力,但SwpTi=0;¯SspNj,¯SspTj分别为p时刻作用在结点j固相的界面法向力和界面切向力;¯SwpNj,¯SwpTj为p时刻作用在结点j液相的界面法向力和界面切向力,但¯SwpTj=0。
根据界面连续条件,推导可得p时刻作用在结点i固、液相的界面力为
SspNi=A22B1−A12B2A22A11−A12A21 , (15) SwpNi=A21B1−A11B2A21A12−A11A22 。 (16) 其中,
A11=(1msi+1msj)(Δt)2 , (17) A12=(1−βjβi)(Δt)2msj , (18) A21=(βi−βj)(Δt)2msj , (19) A22=(βimwi+β2jβimwj+(βj−βi)2βimsj)(Δt)2 , (20) B1=ni[ni⋅(ˆu(p+1)j−ˆu(p+1)i)] , (21) B2=nini⋅[βj(ˆU(p+1)j−ˆu(p+1)j)−βi(ˆU(p+1)i−ˆu(p+1)j)]。 (22) 式中:下标j为与界面结点i形成结点对的另一饱和介质的界面结点,而ˆu(p+1)i,ˆu(p+1)j,ˆU(p+1)i和ˆU(p+1)j分别为界面结点i和结点j不考虑界面力时的固、液相位移,即方程(9)、(10)代入界面力为0后求得的p+1时刻的位移;ni为界面结点i所在界面法向向量。求得SspNi,SspNi后,可进一步由界面连续条件解得¯SspNj,¯SwpNj。再根据固相位移连续条件,求得结点i界面切向力为
SspTi=(ˆu(p+1)j+Δu(p+1)Nj−ˆu(p+1)i−Δu(p+1)Ni)msimsj(msi+msj)(Δt)2 。 (23) 其中,
Δu(p+1)Ni=(Δt)2msiSspNi , (24) Δu(p+1)Nj=(Δt)2msjSspNj 。 (25) 通过界面连续条件求得结点j固相切向界面力¯SspTk,并按方程(11)~(14)求得界面点的位移响应。
(2)人工边界点
由于透射边界不依赖于具体的波动方程,具有普适性和使用简便的特点,本文人工边界条件采用多次透射公式(multi-transmitting formula,简称MTF)。
边界节点的总位移场u和U通过波场分解分为散射场us,Us和自由场uf,Uf,即
u=us+uf , (26) U=Us+Uf 。 (27) 其中,自由场位移uf,Uf可根据2.2节求得。利用多次透射技术,模拟外行波向外散射,可得边界节点i的散射场位移为
u(p+1)is=N∑k=1(−1)k+1CNku(p+1−k)ks , (28) U(p+1)is=N∑k=1(−1)k+1CNkU(p+1−k)ks 。 (29) 其中CNk为二项系数:
CNk=N!k!(N−k)! 。 (30) 节点k为沿边界节点i处的法线方向指向计算区内部的点,其散射位移可由式(26),(27)式求得,N为透射阶数。
由式(28),(29)获得边界点散射场位移后us,Us,便可由式(26),(27)求得边界点的总位移u,U。
1.4 风机结构分析
海上风机结构采用Newmark隐式积分法计算。从海域场地计算区获得风机计算区边界结点p+1时刻反应,则其计算区内p+1时刻动力平衡方程为
(K+1β(Δt)2M+γβΔtC)up+1=Fp+1+M⋅(1β(Δt)2up+1βΔt˙up+(12β−1)¨up)+C(γβΔtup+(γβ−1)˙up+(γ2β−1)Δt¨up) 。 (31) 式中:γ,β为Newmark积分格式参数。
1.5 耦合作用分析
海水-饱和海床-风机耦合作用分析通过分区并行方法实现,即海水、海床和风机基础采用广义饱和介质理论统一框架描述,作为显式区通过集中质量显式有限元进行分析;风机上部结构作为隐式计算区采用隐式有限元进行分析;风机上部结构与海域场地之间的耦合分析通过设置一层显-隐式重叠单元来实现,如图 3所示。
具体的交互方式如下:若计算区p及以前时刻响应已知,则根据2.3节方法可将显式计算区响应更新至p+1时刻,而隐式计算区边界点是显式计算区内部点,可将显式计算区内部点p+1时刻响应赋给隐式计算区边界点;再根据式(36)求得隐式计算区内部点p+1时刻响应,显式计算区边界点是隐式计算区内部点,将隐式计算区内部点p+1时刻响应赋给显式计算区边界点,得到显式计算区所有点p+1时刻响应。循环以上步骤,则可得到系统各时刻响应。
2. 近海单桩式风机地震响应分析
2.1 模型与地震动输入
(1)场地与风机模型
本文以江苏如东北部海域为例,水平成层场地,计算区尺寸取40 m×40 m×70 m,其中海水层厚20 m,饱和海床层厚40 m,基岩层厚10 m;选用额定功率为5 MW的风机,风机基础采用大直径单桩基础,场地和基础介质参数见表 1所示。
表 1 环境介质参数表Table 1. Parameters of environmental media材料 剪切模量/
(108 Pa)泊松比 密度/
(kg·m-3)阻尼比 剪切波速/
(m·s-1)海水 0 0.05 1000 — 0 海床 2.50 0.30 2000 0.03 445 基岩 48 0.20 2643 0.05 1348 单桩 808 0.30 7850 0.05 3208 风机模型参数采用文献[21]中的数据,具休如下:塔架高度为77.6 m,桩露水面高度为10 m,轮毂高度为90 m,桩入海床深度为60 m,塔筒顶部,底部直径分别为3.87,6 m,塔筒顶部,底部壁厚分别为0.019,0.027 m,单桩直径,壁厚分别为6,0.06 m,机舱技师为240 t,转子质量+轮毂质量为110 t,塔筒弹性模量为2.1 GPa,塔筒泊松比为0.30,塔筒密度为8500 kg/m3,转子直径为126 m。
采用ANSYS软件建模,其中机舱和塔筒采用BEAM189单元,叶片采用SHELL63单元,基础采用SOLID185单元建模。对结构整体进行模态分析,控制塔筒反应的前六阶固有频率分别为1.1736,1.2037,2.1105,3.1873,3.5426,4.2909 Hz。场地和基础通过自编程序建模,采用六面体实体单元,单元尺寸为1 m×1 m×1 m,以SV波垂直入射,时间步距Δt =4×10-5 s,满足波动模拟的精度要求。近海单桩式风机模型如图 4所示。
为考虑海水深度、海床波速、入射角度3个因素对海上风机地震响应的影响,通过控制变量设计了12个工况,如表 2所示。其中工况1~工况4控制变量为海水深度;工况5~工况8控制变量为海床波速;工况9~工况12控制变量为入射角度。
表 2 分析工况表Table 2. Working conditions工况 海水深度/m 海床波速/(m·s-1) 入射角度/
(°)工况1 10 445 0 工况2 15 445 0 工况3 20 445 0 工况4 25 445 0 工况5 20 445 0 工况6 20 753 0 工况7 20 953 0 工况8 20 1128 0 工况9 20 445 0 工况10 20 445 15 工况11 20 445 25 工况12 20 445 35 (2)地震动输入
根据江苏如东海域地震安评成果,此场区为Ⅲ类场地,地震基本烈度为Ⅶ度,设计地震分组为第三组,场地特征周期为0.65 s,场地水平地震动峰值加速度为0.1725g,竖向地震动峰值加速度为0.115g。
根据《海上固定平台规划、设计和推荐作法—荷载抗力系数设计法(增补1):SY/T10009—2002》[22],地震反应谱采用加速度反应谱曲线,临界阻尼比为5%,土类型为C类,已有位于海床表面处的目标加速度反应谱,拟根据NGA-West地震动数据库,选取3条与场地特性符合的地震波进行时程分析,并根据2.2节方法反演地震波至基岩底部,如图 5所示。
假设地震波以SV波垂直入射。第一条地震波为1981年Corinth台站记录的Corinth_Greece近海岸地震波,地震波持续时间41.29 s。第二条地震波为1989年Coyote Lake Dam - Southwest Abutment台站记录的Loma Prieta近海岸地震波,地震波持续时间39.99 s。第三条地震波为2007年Joetsu Yasuzukaku Yasuzuka台站记录的Chuetsu-oki_Japan海底地震波,地震波持续时间59.99 s。3条地震记录水平向加速度时程图和频谱图如图 6所示。
2.2 结果分析
(1)海水深度影响
控制海水深度从10 m变化至15,20,25 m,分析监控点A和B(图 4所示)的响应,探究海水深度对海上风机地震响应的影响。表 3给出了各水深工况下监控点处的位移、速度、加速度和弯矩的最大值。
表 3 不同水深工况下地震反应结果Table 3. Seismic response results under different water depths地震记录 监控点 工况
Hw/m绝对位移峰值/
m绝对速度峰值/
(m·s-1)绝对加速度峰值/
(m·s-2)绝对最大弯矩/
(N·m)地震记录
1A 10 0.0690 0.5379 12.4667 2.7735E8 15 0.0640 0.5564 11.5243 2.6007E8 20 0.0820 0.8348 19.2107 4.1018E8 25 0.0633 0.4518 12.1952 2.8194E8 B 10 0.1000 0.9037 17.8810 1.5733E9 15 0.0908 0.9365 17.2672 1.3769E9 20 0.1046 1.6624 17.6343 2.8165E9 25 0.0995 0.7379 13.8724 1.4435E9 地震记录
2A 10 0.0948 0.5798 12.5576 2.8056E8 15 0.0941 0.5866 12.6743 2.7385E8 20 0.0936 0.9776 21.4524 4.5699E8 25 0.0728 0.5185 13.3567 3.0869E8 B 10 0.1359 0.9042 16.7664 1.6644E9 15 0.1360 0.9331 21.2555 1.3756E9 20 0.1345 2.1096 19.5376 3.3482E9 25 0.1131 0.8231 14.1074 1.5988E9 地震记录
3A 10 0.0922 0.7312 18.5752 3.4591E8 15 0.0941 0.7848 16.6196 3.4940E8 20 0.0834 1.0505 22.7070 4.7975E9 25 0.0671 0.6256 13.8432 3.0876E8 B 10 0.1244 1.1889 29.6767 1.9616E9 15 0.1263 1.2035 26.7635 1.9556E9 20 0.1201 1.9692 26.3036 3.4581E9 25 0.0993 1.0606 20.7832 1.8110E9 由表 3结果可知,风机塔筒顶部点和底部点绝对位移地震记录1输入时,随水深增加先减小再增大,在20 m时到达最大,然后减小;在记录2、3输入时,整体上随水深增加而减小。风机塔筒顶部点和底部点绝对速度及弯矩峰值在水深从10 m变化到15 m时略微增大或减小,在20 m时显著增大,20到25 m时又急剧减小,10,15,25 m水深时差别不大。塔筒顶部点的绝对加速度峰值随水深的变化规律与弯矩一致,但塔筒顶部点的加速度峰值在水深为10,15,20 m时变化不是很明显,但在30 m时明显较小。海水深度对风机地震反应的影响规律受地震输入影响较小,可能是3条记录的频谱特性差异不大,频谱主要集中在5 Hz范围内(如图 6);但不同地震记录下风机反应大小还是有差异,如地震记录3输入下的绝对加速度峰值和弯矩峰值要明显大于地震记录1和2输入时的反应值。海水深度的变化,一方面改变了海域自由场(地震输入),另一方面,改变了海域场地-风机整体体系的自振特性,从而影响风机结构的反应;风机塔筒顶部点和底部点的位移、速度、加速度和弯矩在海水深度为20 m时达到了峰值,应该是在水深为20 m时海域场地-风机整体体系的自振频率更接近地震输入的卓越频率。
(2)海床波速影响
在不同海床波速下的塔筒顶部结点A位移响应和塔筒底部结点B弯矩响应如图 7~9所示。可以发现,随着海床剪切波速增大,位移时程图 7(a)、8(a)、9(a)的峰值先增大后减小,地震记录1输入时,在海床剪切波速为753 m/s时达到最大,而记录2和3输入时,在953 m/s时到达最大,由位移频谱图 7(b),8(b),图 9(b)可看出,与地震记录频谱有关。随海床剪切波速增加,最大弯矩先增大后减小,在753 m/s时达到最大。对比位移时程图和弯矩时程图可以发现,弯矩时程图高频明显,因为弯矩与加速度有关,与实际“加速度信号对高频敏感,位移信号对低频敏感”相符,并且可以发现随着波速慢慢增大,塔底倾覆弯矩先小幅增大再逐渐减小,并且当波速为953 m/s,风机塔底倾覆弯矩与其余3种工况方向相反。
(3)入射角度影响
随着入射角度的改变,水平向位移和竖向位移均会发生变化。从图 10~12可见,在不同地震输入下,随着入射角度的增加,风机地震响应变化大体一致,其中竖向位移和加速度反应越来越大,水平向位移和加速度反应总体上越来越小,最大弯矩越来越小。水平向位移变化略有不同,其中,在地震记录1下,随着入射角度增大,水平向位移逐渐减小,而在地震记录2和3下,随着入射角度增大,水平向位移先大幅减小再小幅增大。可见,当改变入射角度时,风机地震响应结果变化受地震输入的影响较小。随着入射角度的增加,出现上述变化规律是因为随着入射角度在xz平面下的增大,x方向分量逐渐减小,z方向分量逐渐增大。故x方向位移和加速度以及受x方向力影响的弯矩均有所减小,而z方向的加速度有所增大。
3. 结论
本文基于广义饱和多孔介质统一计算框架,提出了一套海水-海床-近海单桩式风机耦合作用分区分析方法,并分析了海水深度、海床波速和入射角度对海上风机塔筒地震响应的影响,得到如下3点结论。
(1)海水深度的变化,一方面改变了海域自由场(地震输入),另一方面,改变了海域场地-风机整体体系的自振特性,从而影响风机结构的反应。因此,海水深度对风机地震响应影响较大;当海水增加到一定深度时(如本文算例中20 m时),海域场地-风机整体体系的自振频率与输入地震波频率接近时,风机反应最大。
(2)海床剪切波速对风机地震响应的影响机理本质上与海水深度相同。相对来说,海床剪切波速对风机倾覆弯矩的影响要比塔顶位移大。
(3)随着地震波入射角度的增大,水平向位移和加速度以及倾覆弯矩均有不同程度的减小,竖向位移和加速度均有不同程度的增大。
本文算例仅为线性情形,未考虑海床的非线性以及风机结构的非线性;非线性情形时海水深度、剪切波速、入射角度等因素对风机反应的影响规律应会不同,如斜入射会增加塔顶的竖向加速度,若考虑P-Δ效应,则会加大结构水平向反应,需进一步研究。
-
-
[1] WANG X Z, JIAO Y Y, WANG R, et al. Engineering characteristics of the calcareous sand in Nansha Islands, South China Sea[J]. Engineering Geology, 2011, 120(1/2/3/4): 40-47. http://www.onacademic.com/detail/journal_1000034075073310_47e1.html
[2] XIAO Y, WANG L, JIANG X, et al. Acoustic emission and force drop in grain crushing of carbonate sands[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2019, 145(9): 04019057. doi: 10.1061/(ASCE)GT.1943-5606.0002141
[3] SINGH S C, CARTON H, TAPPONNIER P, et al. Seismic evidence for broken oceanic crust in the 2004 sumatra earthquake epicentral region[J]. Nature Geoscience, 2008, 1(11): 777-781. doi: 10.1038/ngeo336
[4] 胡进军, 郝彦春, 谢礼立. 潜在地震对我国南海开发和建设影响的初步考虑[J]. 地震工程学报, 2014, 36(3): 616-621. https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201403031.htm HU Jinjun, HAO Yanchun, XIE Lili. Effects of potential earthquakes on construction and development in South China Sea region[J]. China Earthquake Engineering Journal, 2014, 36(3): 616-621. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201403031.htm
[5] VAHDANI S, PYKE R, SIRIPRUSANEN UJTRN. Liquefaction of Calcareous Sands and Lateral Spreading Experienced in Guam as A Result of the 1993 Guam Earthquake[C]// 5th US- Japan workshop on Earthquake resistant design of lifeline facilities and countermeasures against soil liquefaction, Snowbird Utah, 1994: 117-123.
[6] SEED HB, IDRISS IM. Soil Moduli and Damping Factors for Dynamic Response Analysis[R]. Berkeley: Earthquake Engineering Research Center (EERC), 1970.
[7] 高冉, 叶剑红. 中国南海吹填岛礁钙质砂动力特性试验研究[J]. 岩土力学, 2019, 40(10): 3897-3908, 3919. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910024.htm GAO Ran, YE Jianhong. Experimental investigation on the dynamic characteristics of calcareous sand from the reclaimed coral reef islands in the South China Sea[J]. Rock and Soil Mechanics, 2019, 40(10): 3897-3908, 3919. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910024.htm
[8] 马维嘉, 陈国兴, 秦悠, 等. 初始主应力方向角对饱和珊瑚砂液化特性影响的试验[J]. 岩土工程学报, 2020, 42(3): 592-600. doi: 10.11779/CJGE202003022 MA Weijia, CHEN Guoxing, QIN You, et al. Experimental studies on effects of initial major stress direction angles on liquefaction characteristics of saturated coral sand[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(3): 592-600. (in Chinese) doi: 10.11779/CJGE202003022
[9] ZHANG Y L, DING X M, CHEN Z X, et al. Seismic responses of slopes with different angles in coral sand[J]. Journal of Mountain Science, 2021, 18(9): 2475-2485. doi: 10.1007/s11629-020-6546-9
[10] MURFF J D. Pile capacity in calcareous sands: state if the art[J]. Journal of Geotechnical Engineering, 1987, 113(5): 490-507. doi: 10.1061/(ASCE)0733-9410(1987)113:5(490)
[11] WU Q, DING X M, CHEN Z X, et al. Shaking table tests on seismic responses of pile-soil-superstructure in coral sand[J]. Journal of Earthquake Engineering, 2022, 26(7): 3461-3487. doi: 10.1080/13632469.2020.1803160
[12] LIU H X, ZHANG J M, ZHANG X D, et al. Seismic performance of block-type quay walls with liquefiable calcareous sand backfill[J]. Soil Dynamics and Earthquake Engineering, 2020, 132: 106092. doi: 10.1016/j.soildyn.2020.106092
[13] 何稼, 楚剑, 刘汉龙, 等. 微生物岩土技术的研究进展[J]. 岩土工程学报, 2016, 38(4): 643-653. doi: 10.11779/CJGE201604008 HE Jia, CHU Jian, LIU Hanlong, et al. Research advances in biogeotechnologies[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(4): 643-653. (in Chinese) doi: 10.11779/CJGE201604008
[14] 彭劼, 温智力, 刘志明, 等. 微生物诱导碳酸钙沉积加固有机质黏土的试验研究[J]. 岩土工程学报, 2019, 41(4): 733-740. doi: 10.11779/CJGE201904017 PENG Jie, WEN Zhili, LIU Zhiming, et al. Experimental research on MICP-treated organic clay[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(4): 733-740. (in Chinese) doi: 10.11779/CJGE201904017
[15] 桂跃, 吴承坤, 刘颖伸, 等. 利用微生物技术改良泥炭土工程性质试验研究[J]. 岩土工程学报, 2020, 42(2): 269-278. doi: 10.11779/CJGE202002008 GUI Yue, WU Chengkun, LIU Yingshen, et al. Improving engineering properties of peaty soil by biogeotechnology[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(2): 269-278. (in Chinese) doi: 10.11779/CJGE202002008
[16] 欧孝夺, 莫鹏, 江杰, 等. 生石灰与微生物共同固化过湿性铝尾黏土试验研究[J]. 岩土工程学报, 2020, 42(4): 624-631. doi: 10.11779/CJGE202004004 OU Xiaoduo, MO Peng, JIANG Jie, et al. Experimental study on solidification of bauxite tailing clay with quicklime and microorganism[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 624-631. (in Chinese) doi: 10.11779/CJGE202004004
[17] 谢约翰, 唐朝生, 尹黎阳, 等. 纤维加筋微生物固化砂土的力学特性[J]. 岩土工程学报, 2019, 41(4): 675-682. doi: 10.11779/CJGE201904010 XIE Yuehan, TANG Chaosheng, YIN Liyang, et al. Mechanical behavior of microbial-induced calcite precipitation(MICP)- treated soil with fiber reinforcement[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(4): 675-682. (in Chinese) doi: 10.11779/CJGE201904010
[18] LIU Hanlong, XIAO Peng, XIAO Yang, et al. State-of-the-art review of biogeotechnology and its engineering applications[J]. Journal of Civil and Environmental Engineering, 2019, 41(1): 1-14. http://en.cnki.com.cn/Article_en/CJFDTotal-JIAN201901001.htm
[19] XIAO Y, HE X, EVANS T M, et al. Unconfined compressive and splitting tensile strength of basalt fiber–reinforced biocemented sand[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2019, 145(9): 04019048. http://www.researchgate.net/publication/334837170_Unconfined_Compressive_and_Splitting_Tensile_Strength_of_Basalt_Fiber-Reinforced_Biocemented_Sand
[20] XIAO Y, ZHANG Z, STUEDLEIN A W, et al. Liquefaction modeling for biocemented calcareous sand[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2021, 147(12): 04021149. http://www.nstl.gov.cn/paper_detail.html?id=5abfd949bb8b9a4d39182baeb0a7309f
[21] 刘汉龙, 肖鹏, 肖杨, 等. MICP胶结钙质砂动力特性试验研究[J]. 岩土工程学报, 2018, 40(1): 38-45. doi: 10.11779/CJGE201801002 LIU Hanlong, XIAO Peng, XIAO Yang, et al. Dynamic behaviors of MICP-treated calcareous sand in cyclic tests[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(1): 38-45. (in Chinese) doi: 10.11779/CJGE201801002
[22] FATTAHI S M, SOROUSH A, HUANG N. Biocementation control of sand against wind erosion[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2020, 146(6): 04020045. http://www.researchgate.net/publication/340607105_Biocementation_Control_of_Sand_against_Wind_Erosion
[23] JIANG N J, SOGA K. The applicability of microbially induced calcite precipitation (MICP) for internal erosion control in gravel–sand mixtures[J]. Geotechnique, 2017, 67: 42-55. http://www.onacademic.com/detail/journal_1000039580441110_860b.html
[24] CLARÀ SARACHO A, HAIGH S K, EHSAN JORAT M. Flume study on the effects of microbial induced calcium carbonate precipitation (MICP) on the erosional behaviour of fine sand[J]. Géotechnique, 2021, 71(12): 1135-1149. http://www.researchgate.net/publication/342933967_Flume_study_on_the_effects_of_microbial_induced_calcium_carbonate_precipitation_MICP_on_the_erosional_behaviour_of_fine_sand
[25] XIAO P, LIU H L, XIAO Y, et al. Liquefaction resistance of bio-cemented calcareous sand[J]. Soil Dynamics and Earthquake Engineering, 2018, 107: 9-19. http://www.researchgate.net/profile/Yang_Xiao27/publication/322537057_Liquefaction_resistance_of_bio-cemented_calcareous_sand/links/5a5ea9b0a6fdcc68fa99305a/Liquefaction-resistance-of-bio-cemented-calcareous-sand.pdf
[26] XIAO P, LIU H L, STUEDLEIN A W, et al. Effect of relative density and biocementation on cyclic response of calcareous sand[J]. Canadian Geotechnical Journal, 2019, 56(12): 1849-1862. doi: 10.1139/cgj-2018-0573
[27] ZHANG X L, CHEN Y M, LIU H L, et al. Performance evaluation of a MICP-treated calcareous sandy foundation using shake table tests[J]. Soil Dynamics and Earthquake Engineering, 2020, 129: 105959. http://www.sciencedirect.com/science/article/pii/S0267726119306281
[28] 肖鹏. 微生物温控加固钙质砂动力与液化特性研究[D]. 重庆: 重庆大学, 2020. XIAO Peng. Study on Dynamic and Liquefaction Charcteristics of Temperature Controlled MICP-Treated Calcareous Sand[D]. Chongqing: Chongqing University, 2020. (in Chinese)
[29] PRASAD S, TOWHATA I, CHANDRADHARA G P, et al. Shaking table tests in earthquake geotechnical engineering[J]. Current Science, 2004, 87: 1398-1404. http://www.nstl.gov.cn/paper_detail.html?id=37706bab806aa8423f61c9fca70631ae
[30] TSUKAMOTO Y, ISHIHARA K, SAWADA S, et al. Settlement of rigid circular foundations during seismic shaking in shaking table tests[J]. International Journal of Geomechanics, 2012, 12(4): 462-470.
[31] DARBY K M, BOULANGER R W, DEJONG J T, et al. Progressive changes in liquefaction and cone penetration resistance across multiple shaking events in centrifuge tests[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2019, 145(3): 04018112. http://www.xueshufan.com/publication/2908166530
[32] TURAN A, HINCHBERGER S D, EL NAGGAR H. Design and commissioning of a laminar soil container for use on small shaking tables[J]. Soil Dynamics and Earthquake Engineering, 2009, 29(2): 404-414. http://www.sciencedirect.com/science/article/pii/S0267726108000584
[33] KOGA Y, MATSUO O. Shaking table tests of embankments resting on liquefiable sandy ground[J]. Soils and Foundations, 1990, 30(4): 162-174. http://www.onacademic.com/detail/journal_1000039882878110_82e5.html
[34] MONTOYA B M, DEJONG J T, BOULANGER R W. Dynamic response of liquefiable sand improved by microbial-induced calcite precipitation[J]. Géotechnique, 2013, 63(4): 302-312. http://smartsearch.nstl.gov.cn/paper_detail.html?id=f94d6a520c012ed5cd3e532fa707b176
[35] AFACAN K B, BRANDENBERG S J, STEWART J P. Centrifuge modeling studies of site response in soft clay over wide strain range[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2014, 140(2): 04013003.
[36] FENG K, MONTOYA BM. Influence of confinement and cementation level on the behavior of microbial-induced calcite precipitated sands under monotonic drained loading[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 142(1): 04015057. http://smartsearch.nstl.gov.cn/paper_detail.html?id=b35ac2ea3a19d5c35aa8fa22eea2c772
[37] CONLEE C T, GALLAGHER P M, BOULANGER R W, et al. Centrifuge modeling for liquefaction mitigation using colloidal silica stabilizer[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2012, 138(11): 1334-1345.
[38] DIHORU L, BHATTACHARYA S, MOCCIA F, et al. Dynamic testing of free field response in stratified granular deposits[J]. Soil Dynamics and Earthquake Engineering, 2016, 84: 157-168. http://www.researchgate.net/profile/Luiza_Dihoru/publication/298710335_Dynamic_testing_of_free_field_response_in_stratified_granular_deposits/links/573301d708aea45ee838d3a8.pdf
[39] HARDIN B O, DRNEVICH V P. Shear modulus and damping in soils: design equations and curves[J]. Journal of the Soil Mechanics and Foundations Division, 1972, 98(7): 667-692. http://www.researchgate.net/profile/Bobby_Hardin/publication/247509918_Shear_Modulus_and_Damping_in_Soils_Design_Equations_and_Curves/links/5762f0ab08aecb4f6fee03a0/Shear-Modulus-and-Damping-in-Soils-Design-Equations-and-Curves.pdf
[40] WANG K, BRENNAN A. Behaviour of saturated fibre- reinforced sand in centrifuge model tests[J]. Soil Dynamics and Earthquake Engineering, 2019, 125: 105749. http://www.xueshufan.com/publication/2953448460
[41] ELGAMAL A, YANG Z H, LAI T, et al. Dynamic response of saturated dense sand in laminated centrifuge container[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2005, 131(5): 598-609. http://www.researchgate.net/profile/Bruce_Kutter/publication/238179978_Dynamic_Response_of_Saturated_Dense_Sand_in_Laminated_Centrifuge_Container/links/55598dd608aeaaff3bf99f8b.pdf
[42] 张鑫磊, 陈育民, 张喆, 等. 微生物灌浆加固可液化钙质砂地基的振动台试验研究[J]. 岩土工程学报, 2020, 42(6): 1023-1031. doi: 10.11779/CJGE202006005 ZHANG Xinlei, CHEN Yumin, ZHANG Zhe, et al. Performance evaluation of liquefaction resistance of a MICP-treated calcareous sandy foundation using shake table tests[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(6): 1023-1031. (in Chinese) doi: 10.11779/CJGE202006005
[43] 刘汉龙, 马国梁, 肖杨, 等. 微生物加固岛礁地基现场试验研究[J]. 地基处理, 2019, 1(1): 26-31. https://www.cnki.com.cn/Article/CJFDTOTAL-DJCL201901007.htm LIU Hanlong, MA Guoliang, XIAO Yang, et al. In situ experimental research on calcareous foundation stabilization using MICP technique on the reclaimed coral reef islands[J]. Journal of Ground Improvement, 2019, 1(1): 26-31. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DJCL201901007.htm
[44] VAN PAASSEN L A, GHOSE R, VAN DER LINDEN T J M, et al. Quantifying biomediated ground improvement by ureolysis: large-scale biogrout experiment[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2010, 136(12): 1721-1728. http://www.researchgate.net/profile/Mark_Van_Loosdrecht/publication/221932480_Quantifying_biomediated_ground_improvement_by_ureolysis_Large-scale_biogrout_experiment/links/5640dda608ae24cd3e40b8c4.pdf
-
期刊类型引用(19)
1. 梁贤伟,郭子亮,袁梓文,穆保岗. 基于MRA-BP神经网络的桩基托换桥墩沉降预测模型研究. 特种结构. 2025(01): 82-86+97 . 百度学术
2. 杜策,张振,张力,周振,叶怀文. 桩基托换对邻近深基坑历史建筑的变形影响控制研究. 施工技术(中英文). 2024(23): 23-28 . 百度学术
3. 李旭,郑少强,佟庆伟,绳虎. 多种支护模式下基坑监测数据处理与分析. 山西建筑. 2023(05): 80-84 . 百度学术
4. 徐现启,马炜,邓才兵,孙书淳. 岩质顺层破碎采石场开挖边坡稳定性分析. 云南水力发电. 2023(02): 17-24 . 百度学术
5. 田书广,高佳豪,周喻,李想,裴博文,刘汝辉,田峰,李坤. 不断流情况下明挖基坑穿河施工变形特性研究. 建筑结构. 2023(S2): 2024-2032 . 百度学术
6. 孙若翔. 复杂环境下城市建筑深基坑变形与数值模拟分析. 江西建材. 2022(01): 37-40 . 百度学术
7. 臧一平,卜成,邵苏文,孙晋晶. 地铁盾构隧道下穿桥梁桩基被动托换加固研究. 山西建筑. 2022(06): 139-143 . 百度学术
8. 高成雷. 城区明挖基坑地下连续墙支撑体系稳定性研究及监测分析. 建筑技术开发. 2022(10): 148-151 . 百度学术
9. 阮小勇. 桥梁桩基的门架式托换结构及工作特性研究. 科技创新与应用. 2022(31): 82-88 . 百度学术
10. 陈治雄. 福州火车站南广场深基坑支护和土方开挖施工关键技术. 建筑技术开发. 2022(19): 139-142 . 百度学术
11. 杨正华,段军朝,郭庆军,鄢玉胜. 异型变宽连续梁桩基托换变形监测分析. 城市轨道交通研究. 2022(12): 106-111 . 百度学术
12. 李红,王初生. 隧道桩基主动托换加固效果仿真评价方法. 计算机仿真. 2021(04): 188-192 . 百度学术
13. 朱会强,张明,张波,张建设,郭波锋. 黄河中下游冲洪积地层PRC管桩深基坑支护. 河南科学. 2021(06): 964-970 . 百度学术
14. 陈祉阳,王文柱. 地铁车站明挖基坑下穿互通立交桥基础的施工关键技术及变形影响分析. 建筑结构. 2021(S1): 1952-1956 . 百度学术
15. 吴毅彬,郑伟. 某城际铁路车站天桥主动托换施工方案. 厦门理工学院学报. 2021(03): 62-69 . 百度学术
16. 闫强,廉向东,凌建明. 边坡开挖支护时序有限元分析. 交通运输工程学报. 2020(03): 61-71 . 百度学术
17. 郭在旭,向俐蓉,刘会娟,孙晓. 藏东南砂层路堑分层开挖下钻孔桩响应研究. 高速铁路技术. 2020(05): 51-56 . 百度学术
18. 陈伟志,蒋关鲁,刘勇,陈虹羽. 川藏铁路钻孔桩加固斜坡路基的振动台试验研究. 岩石力学与工程学报. 2020(12): 2540-2556 . 百度学术
19. 付艳军. 车站基坑开挖对邻近立交桥的影响特性研究. 云南水力发电. 2019(06): 117-121 . 百度学术
其他类型引用(5)
-
其他相关附件