海上风电单桩基础周围土体地震液化分析.pdf

第4 7 卷第4 期 2 0 1 7 年4 月 中国海洋大学学报 P E R l0 D I C A L0 Fo C E A NU N I V E R S I T YO FC H I N A 4 7 ( 4 ) ：0 9 3 0 9 9 A p r ，2 0 1 7 海上风电单桩基础周围土体地震液化分析。 刘红军1 ，王超2 ( 1 中国海洋大学海洋环境与生态教育部重点实验室，山东青岛2 6 6 1 0 0 ；2 中国海洋大学环境科学与工程学院，山东青岛2 6 6 1 0 0 ) 摘要：基于室内土工动三轴土体液化试验，采用有限元分析计算与剪应力比较法相结合的综合判别方法，对于黄河三 角洲地区海上风电单桩基础周围土体是否会液化进行了判断，并且通过对比不同位置处土体单元受力情况，总结出单桩基 础的存在对于其周围土体地震荷载作用下动力响应的影响。研究发现，土体液化只会在一定深度内发生，1 2m 以下土体 由于上覆有效应力较大基本不会发生液化。单桩基础对于桩内土体约束作用十分明显，同时会剧烈加大桩侧土体的受力。 本文通过对土体施加x 轴方向上的加速度的方式模拟地震作用，结果显示，位于单桩基础x 轴方向上的土体单元相较于 y 轴方向上等距离的土体单元，受到的剪应力振荡幅度要大的多。 关键词： 海上风电；单桩基础；A B A Q U S ；地震荷载；抗液化剪应力 中图法分类号：P 7 3 6 1 文献标志码： A 文章编号： 1 6 7 2 5 1 7 4 ( 2 0 1 7 ) 0 4 0 9 3 0 7 D O I ： 1 0 1 6 4 4 1 j c n k ih d x b 2 0 1 6 0 1 1 8 引用格式：刘红军，王超海上风电单桩基础周围土体地震液化分析 J 中国海洋大学学报( 自然科学舨) ，2 0 1 7 ，4 7 ( 4 ) ：9 3 9 9 L I UH o n g - J u n ，W A N GC h a o S t u d yo ns e is m icr e s p o n s ea n dliq u e f a c t io no fs o ila r o u n dp ilef o u n d a t io no fo f f s h o r ew in d p o w e r J P e r io d ic a lo fO c e a nU n iv e r s it yo fC h in a ，2 0 1 7 ，4 7 ( 4 ) ：9 3 9 9 风能是是一种清洁、可再生的能源，尤其是海上风 力发电具有风速大且稳定、不占用陆地资源、建设周期 短、运营成本低等优势，越来越受到重视。相关统计显 示，到2 0 1 3 年底为止，欧洲地区已有6 9 座海上风电场 20 8 0 台海上风机运营，总装机容量超过66 0 0M W ，此 外另有超过30 0 0M W 的海上风机正在建设中 1 。与 此相对的，中国目前已有海上风电总装机容量约为4 3 0 M W 。据最新公布的数据显示2 0 1 4 - - 2 0 1 6 年中国计划 要建成总装机容量10 5 3M W 的海上风机 2 | 。 国家发改委会曾经颁布过风能资源评价技术规 定，年有效风能密度能达到2 0 0W m 2 以上则属于风能 丰富区。据其测定，黄河三角洲的风能年平均有效小 时数和年平均风功率密度均符合标准。尤其是滨州 市，海岸线长达2 4 0k m ，是中国重要的风能地区之 一L 3 。目前黄河三角洲地区风电项目已在规划中。海 上风电项目基础型式丰富，其中应用最为广泛的是单 桩基础。单桩基础施工周期短，施工工艺较为简单且 价格较低，而且建成后占用海床面积最小，具有承载力 高、沉降量小且均匀的特点。目前欧洲建成和规划的 海上风电机组的6 0 左右都是采用单桩基础。东海大 桥海上风电场( 中国第一个建成的海上风电项目) 采用 了直径4 8I n ，长度5 0m 的钢管桩 4 。黄河三角洲水 深小于2 0m 的滩涂和浅海区域广阔，适用目前应用较 广的海上风电单桩基础。 许多学者已经就海上风电单桩基础展开了研究并 取得了一定的成果。祁德庆等对非线性波浪作用下风 机单桩基础结构响应进行研究，分析出结构顶部的位 移频率谱曲线 5 。荣维栋等采用三种地震波和两种水 深的波浪荷载任意组合的方式，通过A N A Y S 有限元 模拟完成了海上风电单桩基础在波浪地震联合作用下 的动力特性分析 6 。N A la t i等对带三脚架和钢支架 的单桩基础在地震作用下的动力响应的研究证明了地 震荷载对于桩基础受力影响显著，并对于海上工作承 台的地震危险性进行了风险评估 7 。袁宇等综合利用 A N A Y S 有限元模拟及p - y 曲线法，在充分考虑桩土相 互作用的前提下得到了海上风机全斜桩承台与全直桩 承台在地震动力荷载作用下的动力响应 8 。I A n a s t a s o p o u lo s 等通过施加地震荷载和环境荷载研究了不 * 基金项目：国家自然科学基金项目“波浪作用下海上风电场桩基土体液化响应研究一以黄河三角洲为例”( 4 1 5 7 2 2 4 7 ) ；山东省科技攻关资助项目 “黄河三角洲海上风电场桩基设计关键技术研究”( 2 0 1 4 G G X l0 4 0 0 7 ) 资助 S u p p o r t e db yN a t io n a lN a t u r eS c ie n c eF u n d sR e s e a r c ho nL iq u e f a c t io nR e p o n s eo fS o ilB o d yA r o u n dP ileF o u n d a t io no fO f f s h o r eW in d T u r b in eU n d e rW a v eC o n d it io n s - - T a k in gY e llo wR iv e r D e lt aa st h eC a s e ( 4 1 5 7 2 2 4 7 ) ，R e s e a r c ho nK e yT e c h n iq u e sinP ileF o u n d a t io n D e s ig no fO f f s h r eW in dF a r minY e llo wR iv e rD e lt a ( 2 0 1 4 G G X l0 4 0 0 7 ) 收稿日期：2 0 1 6 - 0 4 0 7 ；修订日期：2 0 1 6 0 6 2 7 作者简介：刘红军( 1 9 6 6 一) ，男，教授，博导，主要从事海洋工程地质方面的科研工作。E - m a il：h o n g j u n O U C e d u c a 中国海洋大学学报 同深度的单桩基础的动力响应的区别 9 。目前已有的 研究对于海上风机单桩基础动力荷载作用下的分析大 都集中在对于桩身内力及位移分析方面，对于桩周土 的研究尚不多见。黄河三角洲地区地处饱和粉土带， 饱和细粒土在地震这种频率较高的动力荷载作用下很 容易发生较大规模的液化现象。本文重点分析桩周地 基土在地震荷载作用下的动力响应，进行液化判别。 研究成果可以为海上风电场建设提供一定参考。 1 有限元模型 本文首先对于地震荷载作用下的海上风电单桩基 础周围土体的动力响应进行分析计算，选用大型有限 元计算软件为A B A Q U S 。A B A Q U S 是一款功能强大 的有限元软件，拥有多种多样的模型库，可以模拟很多 类型的材料，其中包括金属、橡胶、复合材料、土壤和岩 石等，尤其是在岩土工程领域内，拥有多种岩土本构模 型可供选择，弹性或是弹塑性模型计算都可以很好的 完成。 1 1 物理参数及模型设置 单桩基础为钢管桩，直径6m ，厚度6 0in t o ，长4 0m ， 入土2 0m 深。土体模型直径6 0m ，深度2 5m 。数值 模拟中所需土体参数来自于青东、垦东地区现场所取 原状样在室内土工实验室2 贝4 定。黄河口地区地层条件 十分复杂，为了减少数值模拟过程的难度，本文将土体 划分为两层，分别采用黄河三角洲地区比较有代表性 粉质黏土和粉土。因为土体液化深度有限，为了保证 两层不同土质均有机会参与液化反应，所以上层粉质 黏土层设置为4m ，以下直到2 5I ll均为粉土层。模型 参数见表1 。表1 中lD 为材料密度，E 为弹性模量，p 为 泊松比，C 为粘聚力，p 为内摩擦角。 表1 模型部件及物理参数 T a b le1P a r tin s t a n c e sa n dp h y s ic a lp a r a m e t e r s 部件P a r t s几何尺寸D im e n s io n s 材料属性P r o p e r t ie s 钢管桩厚6。mm，入土深度矗一N t o o l n il a K m L = 4 0 m 直径6 0m ,P 2 0 直径D = 6 0m 深度H 一 粉质黏土4m 的圆柱体，以其圆心 S ilt yc la y 为圆心半径3m 处带有 圆环状钢管桩槽 线弹性材料_ 0 78 0 0 k g m a ，E 一2 1 0G P a ， “一O 3 弹塑性材料，采用M - C 模型，p 一19 8 0k g m 3 E 一5 6 2M P a ，“一 0 2 7 ，c 一1 1 5k P a ，。一 5 1 。 1 2 模型建立 本文建立的模型包括钢管桩和地基土两部分，模 型参数上文已经列出。在有限元模拟过程中必须注意 的还有两个部件之间的接触和相互作用。桩端直接与 土体绑定，桩体表面与土体在动力荷载作用下发生相 对位移时的相互作用遵循M o h r - C o u lo m b 摩擦定律， 摩擦系数取0 3 0 。大直径钢管桩桩内土体对于桩身的 稳定性及承载性能均有较大影响，在进行分析时必须 予以考虑。在本文的模型中，设置桩内土体与桩外土 体高度相同，且对于桩身内侧与桩内土体的接触的与 桩身外侧与土体的相同。网格划分后的模型见图1 。 图1网格剖分示意图 F ig 1 M e s hg e n e r a t io ns k e t c h 1 3 地应力平衡 自然界的海床，不论何种土质，都已经经过了在重 力作用下的很长时期的固结沉降，达到了稳定状态，本 文把这种稳定状态称为初始地应力状态。在这种状态 下，土体内本身具有恰好与重力相平衡的应力。这种 内部应力可以保证土体在不受外力作用时不发生大规 模的沉降。而我们建立的有限元模型是没有内部应力 的。如果直接在建立的模型上加载动力荷载进行计 算，那么土体必然会同时受重力影响发生固结沉降，导 致最后的计算结果中位移与应力都远远大于实际值。 为了有效模拟实际情况，必须在加载动力荷载之前通 过初始条件的设置使土体达到初始地应力状态。但是 知道土体内部应力的真实数值是很困难的，为此采取 的方法为：施加动力荷载前，对整个模型施加重力荷 载，得到的土体内部应力作为初始条件写入模型中，再 施加重力以平衡。经过此过程后，土体在重力荷载下 的位移小于1 0 m ，符合工程中的精度要求。 1 4 地震荷载 本文研究针对黄河三角洲地区，该地区抗震设防 烈度为7 度，根据水利水电工程地质勘察规范 ( G B 5 0 4 8 7 - 2 0 0 8 ) 中规定，对于抗震设防烈度为7 度的 4 期 刘红军，等：海上风电单桩基础周围土体地震液化分析 地区，计算地震作用时地震加速度时程的最大值取0 1 倍的重力加速度，本文约取为1m s 。 1 9 7 6 年天津宁河地震记录以及1 9 4 0 年美国加利 福尼亚州帝谷地震在E IC e n t r o 台站的记录等两列地 震波因为记录较为准确因而被广泛地应用于地震分析 中。为了减少偶然性的计算失误，本文同时采用此两 列地震波。宁河天津地震波记录的最大加速度值为 1 4 5 8m s 2 ，所以应将此地震波的加速度时程曲线同步 除去1 4 5 8 而后得到新的加速度时程曲线并将其加载 到土体上。采用地震波的前1 5S 。波峰1m s 2 出现 在7 6 4S 。同理，E I 地震波亦需要调幅处理，波峰 1m s 2 出现在7 1 1S 。 在利用A B A Q U S 软件计算地震力作用时，在土体 底面z 方向施加随时间变化而变化的加速度，地震加 速度时程曲线如图2 所示。 1 0 1 o 8 M c “打o 0 6J 0 4J 0 2J 时间T im e s 图2 地震加速度时程曲线 F ig 2 T h et im e - p a t hc u r v e so fs e is m ica c c e le r a t io n 2 动三轴试验及剪应力比较法 A B A Q U S 有限元模拟软件拥有前文所述的诸多 优点，但是在进行海床土体流固耦合计算的时候， A B A Q U S 无法输出土体孔隙水压力是一大不足。本 文采用A B A Q U S 有限元计算与剪应力比较法相结合 的方法来进行液化判定。 本文采用胜利油田青东、垦东地区标准地层勘察 项目中在该地区所取原状样进行室内动三轴试验所得 的结果来计算土体抗液化剪应力。 试验采用W 3 Z 昏2 0 型微机控制电液伺服土动三 轴试验机进行测定土样抗液化强度的应力控制振动三 轴试验，激振波形采用正弦波，试样尺寸采用伞一 3 9 5m m X8 0m m ，激振频率选用1H z 。试验中采用3 个固结压力分别为9 5 、1 1 5 和1 2 5K P a 。1 个固结比为 K 。一1 0 。 目前工程上多以初始液化或轴向应变值e 。一5 作为破坏标准，即动三轴试验的终点。在本次试验中， 采用5 轴向应变作为液化破坏标准。 2 1 抗液化剪应力计算 动三轴试验的结果以动剪应力比( a d 2 a 7 。) 和破坏 振次lo g N f 的关系曲线表示( 见图3 、4 ) o0 4 器 O 3 6 上0 3 2 曼o2 8 I - + 割q 0 2 4 幅O 2 1 0 010 0 0 振次V ib r a t io nt im e s 图3 粉质粘土土样动强度拟合曲线 F ig 3D y n a m ics t r e n g t hf it t in gc u r v eo fs ilt yc la y 0 11 01 0 010 0 0 振次V ib r a t io nt im e s 图4 粉土土样动强度拟合曲线 F ig 4 D y n a m ics t r e n g t hf it t in gc u r v eo fs ilt ys o il 天然埋藏状态下土层的抗液化剪应力可由下式来 确定： r C r a7 y 羲。 式中：r 为土体的抗液化剪应力( k P a ) ；C ，为应力校正 系数；仃7 y 为上覆土层自重有效压力( k P a ) ；鳓为破坏动 应力( k P a ) ；d 7 。为土体围压( k P a ) 。 依据上述公式，由动三轴试验数据可以得到黄河 h g nNnaT河宁津天 111 0 8 6 4 2 1 O O O O r O O O O O 1 0 0 0 0 0 中国海洋大学学报 口三角洲地区不同土体深度处抗液化剪应力计算结果 如表2 所示。因为本文所研究土体均为海床土，因而 在计算过程中均取其土体的浮容重。 2 2 最大剪应力计算 S e e d 剪应力比较法 1 0 中剪应力是通过公式将地 震作用转化为等效循环剪应力。本文通过A B A Q U S 来模拟土体地震作用下动力响应进而输出土体单元受 到的最大剪应力，以此与抗液化剪应力对比，通过两者 的大小关系来对土体进行液化判别。 T r e s c a 屈服条件又称最大切应力不变条件，是根 据最大剪应力来判定屈服的。所以不难得出土体单元 的T r e s c a 应力值与土体剪应力值有函数关系如下： 口t = 口1 一= 2 t o 。 式中O t r e s e a 为土体的T r e s c a 应力值。 如前文所述，A B A Q U S 无法输出土体孔隙水压 力。但是A B A Q U S 丰富的模型库和强大的计算及后 处理功能使得其计算的土体应力十分可靠。本文通过 输出土体单元在整个地震作用过程中的T r e s c a 应力时 程曲线。由此可以得到土体单元的剪应力时程变化曲 线，取其最高点即可得到土体单元在整个过程中的最 大剪应力Z - m a x 。因为地震作用无规则，而且最大加速度 作用时间很短，因而取等效剪应力k 为液化判别指标。 S e e d 简化法在提出时对于两者的大小给出了倍数关 系：r a y = o 6 5 ，通过0 6 5 这一修正系数将随机振动 转换为等效均匀循环振动，此后数十年的研究及应用 中均采用这一数值。本文沿用此修正系数对数值计算结 果进行处理。将同一土体单元的等效剪应力r a v 与其抗 液化剪应力进行比较就可以对土体单元是否液化做出判 断。根据所有液化的土体单元即可确定液化的区域。 3 输出结果分析 根据有限元软件输出结果，取各个土体单元在整 个数值模拟过程中最大的T r e s c a 应力，计算出相应的 ，并与之前已经计算得出的抗液化剪应力列表对比， 液化判别结果如表2 所示。 表2 抗液化剪应力计算及液化判别结果 T a b le2R e s u lt so fa n t i- liq u e f a c t io ns h e a rs t r e s s c a lc u la t io na n dL iq u e f a c t io nd e t e r m in a t io n 3 1 单桩基础对土体受力的影响 通过计算输出结果可以明显看出地震荷载的作用 下，模型周边的地基土与围绕桩周围的土体单元的动 力响应是完全不同的。尤其是单桩基础深入地基土达 到2 0m ，必然会对周边土体在动力加载过程中的表现 有很大影响。同时，单桩基础的X 轴方向上的土体与 y 轴上的也有很大区别，这说明单桩基础对于周边土体 既有固定约束作用，又有动力挤压作用，即便士体单元 与单桩基础距离相同但是方向不同受到的综合作用也 是不同的。 3 1 1 选取研究关注的土体单元为了更好的表现 不同位置处土体单元的区别，进而得到单桩基础的存 在对于土体在地震荷载作用下动力响应的作用，本文 选取了几个处于不同位置的土体单元，并标记如下。A 点为桩中心处；B 点在单桩基础的X 轴方向上并且距 离单桩基础5m ；C 点在单桩基础的X 轴方向上并且距 离单桩基础3 0m 。D 、E 、F 点水平相对位置与A 、B 、C 点同时埋深为5m 。与之相对的，G 、H 、N 、M 点相对位 置与以上诸点相同，只是均位于单桩基础y 轴方向上。 如图5 所示。 4 期 刘红军等：海上风电单桩基础周围土体地震液化分析 ( a ) 工轴方向 ( a ) X a x isd ir e c t io n ( b ) r 轴方向 ( b ) _ Fa x isd ir e c t io n 图5 选取的土体单元位置示意图 F ig 5T h elo c a t io n so ft h ec h o o s e ns o ile le m e n t s 将所选取的不同土体单元的T r e s c a 应力时程数据 输出并画成曲线，进行单桩基础不同相对位置处的土 体单元受力情况对比，从而得到单桩基础的存在对于 土体的影响。 3 1 20m 深度不同位置处土体单元受力情况对比 在0m 深度处，分析数据可以得到曲线如图6 、7 、8 。 0 3691 2151 8 时间T im e s 图6A 、B 、C 点土体单元的T r e s c a 应力时程曲线 F ig 6 T h et im e - p a t hc u r v e so fT r e s c as t r e s so fp o in tA 、B 、C A 、B 、C 点同在深度0m 的海床面上，由图6 可以 看出，A 点由于处于钢管桩中，钢管桩对于桩内的土体 单元约束作用很强，桩内土体单元受力情况相对而言 较为稳定。距离单桩基础5m 处B 点的土体单元距离 4 5 4 0 3 5 重3 0 j2 5 当 2 0 R 倒1 5 1 0 5 O 03691 21 51 8 时间T im e s 图7B 、G 点土体单元的T r e s c a 应力时程曲线 F ig 7T h et im e - p a t hc u r v e so fT r e s c as t r e s so fp o in tB 、G 03691 21 51 8 时间T im e s 图8c 、H 点土体单元的T r e s c a 应力时程曲线 F ig 8 T h et im e - p a t hc u r v e so fT r e s c as t r e s so fp o in tC 、H 单桩基础有一定距离，由于单桩基础的撞击作用而受 力变化幅度较大。在距离单桩基础较远的3 0m 处的C 点，受到的更多的是地震通过土体自下而上传导过来 的动力作用，与B 点相比较而言和单桩基础相互作用 的部分减弱。这也解释了为何靠近单桩基础外侧的土 体单元更容易发生液化。 B 、G 两点分别为X 、y 轴方向上距离单桩基础5m 的土体单元，C 、H 两点分别为X 、y 轴方向上距离单桩 基础3 0m 的土体单元。通过图7 与8 可以明显看出， 距离单桩基础距离相同的情况下，由于输入的地震作 用为X 轴方向上的水平加速度，因而单桩基础与在其 x 轴方向上的土体单元之间的动力相互作用更强，即 单桩基础因物理、力学性质的区别，其存在与否会对周 围土体在动力荷载作用下的受力情况产生较大的影 响。 3 1 35m 深度不同位置处土体单元受力情况对比 在土体的5m 深度处，分析数据可以得到如下曲线： 根据图9 、1 0 、1 1 可以看出，5m 深度处土体单元的 受力的趋势及大小关系和前文分析的0m 海床面处大 同小异，只是因为深度处于5m 处，环境对于研究对象 的约束更加严格，因而单桩基础和土体表现出的活动 如 筋 加 ” 加 5 o sso皇R倒 钙 ” 如 筋 加 ：2 m 5 o mdIsojR目 中国海洋大学学报 8 0 7 0 6 0 室5 0 耄 曼4 0 R3 0 垲 2 0 1 0 0 O3 6 91 2 1 51 8 时间T im e s 图9D 、E 、F 点土体单元的T r e s c a 应力时程曲线 F ig 9T h et im e - p a t hc u r v e so fT r e s c as t r e s so fp o in tD 、E 、F 8 0 7 0 6 0 星5 0 呈 羔4 0 R3 0 型 2 0 lO 0 03691 21 51 8 时间T im e s 图1 0E 、M 点土体单元的T r e s c a 应力时程曲线 F ig 1 0 T h et im e - p a t hc u r v e so fT r e s e as t r e s so fp o in tE 、M 03691 21 51 8 时间T im e s 图1 1F 、N 点土体单元的T r e s e a 应力时程曲线 F ig 1 1 T h et im e - p a t hc u r v e so fT r e s c as t r e s so fp o in tF 、N 性都较差，同时因无法通过运动的方式排解从而造成 了应力的积聚。总体而言，51 T I 深度处的土体单元受到 的剪应力在数值上来说是更大的。 4 结论 综上所述，正因为钢管桩在物理、力学性质上与土 体存在巨大区别，所以单桩基础周围土体在动力荷载 作用下的动力响应与自由海床相比较有很大的不同。 根据前述判别结果进行分析，可以得到如下结论： ( 1 ) 地震造成的液化受土体深度的因素影响，0 6m 土层内，土体单元的6 6 7 发生了液化；在6 1 2I n 深 度区间，有1 6 9 7 的土体单元发生了液化；1 2r n 以下 土体单元未发生液化。总体上来看，随着土体深度的 影响，土体所受到的上覆土层压力逐渐增大，导致抗液 化剪应力增大，液化的土体单元数量减少非常明显。 只有深度小于1 21 T I 的土体才会随深度的不同或多或 少有土体单元发生液化；当土体深度大于1 21 T I 时，土 体基本不发生液化，这是因为自重压密的影响。 ( 2 ) 单桩基础桩内的土体单元受力较小，不容易发生液 化，而在单桩基础外部的土体单元则明显受到所处位 置的影响。发生液化的土体单元主要分布在沿X 轴经 过单桩基础的位置，在单桩基础y 方向上的土体尤其 是远离单桩基础的土体单元液化程度明显偏低。 ( 3 ) 单桩基础对于桩内土体具有较强的约束作用，相当 于室内试验中额外添加的围压，所以桩内的土体受到 约束固定作用较为明显；同时严重远离单桩基础的土 体受单桩影响较小，此二者均不易发生液化。而在其 他区域的土体单元，既受到单桩基础的固定约束作用， 又受到振动过程中单桩基础的撞击动力作用，处于单 桩基础y 轴方向上的土体单元受到的固定作用大于动 力作用，液化程度明显较弱；而X 方向上在振动过程中 与单桩基础发生了强烈的相互撞击，因而受到的动力 作用远远大于固定作用，大部分发生了液化。 参考文献： E I - I E w E A - T h eE u r o p e a no f f s h o r ew in din d u s t r y - k e yt r e n d sa n ds t a t is t ic s2 0 1 3 R 3 I s 1 ：E W E A ，2 0 1 4 E 2 3 张海龙中国新能源发展研究 D 长春：吉林大学，2 0 1 4 Z h a n gH a i- lo n g R e s e a r c ho nt h eN e wE n e r g yD e v d o p m e n tinC I li n a E D C h a n g e h u n ：J ilinU n iv e r s it y ，2 0 1 4 3 3 刘清志，范延磊黄河三角洲风电利用分析科学研究 J 农业工 程技术，2 0 1 0 ，2 3 ( 8 ) ：3 0 3 2 L iuQ in g - z h i，F a nY a n - le iA n a ly s iso fw in dp o w e ru t iliz a t io nin t h eY e llo wR iv e rD e lt a E J A g r ic u lt u r eE n g in e e r in gT e c h n o lo g y ， 2 0 1 0 ，2 3 ( 8 ) ：3 0 - 3 2 4 3 刘星，潘培永浅析东海风电基础承台施工技术难点对策及实施 效果l- J 中国水运，2 0 1 2 ，1 2 ( 8 ) ：2 3 1 2 3 4 L iuX in g ，P a nP e i- y o n g A n a ly s iso nt h ec o n s t r u c t io nt e c h n o lo g y o ft h eE a s tC h in aS e aw illdp o w e rb a s ep la t f o r ma n dit sim p le m e n t a t io ne f f e c t J C h in aW a t e rT r a n s p o r t ，2 0 1 2 ，1 2 ( 8 ) ：2 3 1 2 3 4 5 祁德庆，崩兴峰不规则波作用下风机单桩基础结构响应分析口 力学季刊，2 0 1 2 ，2 2 ( 4 ) ：6 5 7 6 6 4 Q iD e - q in g ，B e n gX in g - f e n g A n a ly s iso nt h eo f f s h o r ew in dt u r b in g u n d e rw a v elo a d J C h in e s eQ u a r t e r lyo fM e c h a n ic s ，2 0 1 2 ，2 2 ( 4 ) ：6 5 7 6 6 4 加 如 加 如 如 加 o 芒：o扫R趟 4 期刘红军，等：海上风电单桩基础周围土体地震液化分析 9 9 6 3 E 7 3 8 荣维栋，李洪斌近海单桩风机在波浪地震联合作用下的动力特性 分析 J 江苏科技大学学报( 自然科学版) ，2 0 1 5 ，2 9 ( 1 ) ：2 7 3 2 R o n gW e i- d o n g ，L iH o n g - h in A n a ly s iso nt h ed y n a m icr e s p o n s eo f o f f s h o r ew in dt u r b in eu n d e rt h ej o in ta c t io no fw a v ea n ds e is m ic _ J J o u r n a l o fJ ia n g s uU n iv e r s it yo fS c ie n c ea n dT e c h n o lo g y ( N a t u r a lS c ie n c eE d it io n ) ，2 0 1 5 ，2 9 ( 1 ) ：2 7 3 2 N a t a leM a r l，G iu s e p p eF a illa ，F e lic eA r e n a S e is m ica n a ly s iso f o f f s h o r ew in dt u r b in e so nb o t t o m - f ix e ds u p p o r ts t r u c t u r e s J E n g in e e r in gS t r u c t u r e ，2 0 1 4 ，8 6 ( A 3 7 3 ) ：1 - 2 2 袁宇，黄挺，郑金海考虑桩土非线性作用的近海风机高桩平台地 震响应分析 J 水道港口，2 0 1 6 ，3 7 ( 1 ) ：8 1 8 8 Y a nY u ，H u a n gT in g ，Z h e n gJ in - h a iS e is m icr e s p o n s ea n a ly s iso f h ig h - p ilep la t f o r mo fo f f s h o r e 、i1 1 dt u r b in ec o n s id e r in gn o n lin e a r p ile - s o ilin t e r a c t io n J J o u r n a l o fW a t e r w a ya n dH a r b o r ，2 0 1 6 ， 3 7 ( 1 ) ：8 1 8 8 9 3A n a s t a s o p o u lo sI ，T h e o f ilo uMH y b r idf o u n d a t io nf o ro f f s h o r e w in dt u r b in e s ：E n v ir o n m e n t a la n ds e is m iclo a d in g E J S o ilD y n a m ic sa n dE a r t h q u a k eE n g in e e r in g ，2 0 1 6 ，8 0 ( 2 5 ) ：1 9 2 2 0 9 E 1 0 S E E DH B ，I D R I S SI LS im p lif ie dp r o c e d u r ef o re v a lu a t in gs o il liq u e f a c t io np o t e n t ia l J J o u r n a l o ft h eS o ilM e c h a n ic sa n dF o u n d a t io n sD iv is io n ，A S C E ，1 9 7 1 ，9 7 ( 9 ) ：1 2 4 9 1 2 7 3 S t u d yo nS e is m icR e s p o n s ea n dL iq u e f a c t io no fS o il A r o u n dP ile F o u n d a t io no fO 仟s h o r eW in dP o w e r L I UH