基于BESO算法的高桩承台式水平轴风力机支撑结构优化.pdf

收稿日期:修回日期:录用日期:基金项目:上海市教育委员会科研创新计划自然科学重大项目(E ),国家自然科学基金(,),上海市曙光计划(S G ),上海交通大学深蓝计划(S L P T ),湖南省自然科学基金(J J )作者简介:占玲玉(),硕士生,从事海上风机结构优化、减振研究通信作者:周岱,教授,博士生导师;E m a i l:z h o u d a i s j t u e d u c n 文章编号:()D O I:/j c n k i j s j t u 基于B E S O算法的高桩承台式水平轴风力机支撑结构优化占玲玉,何文君,周岱,韩兆龙,朱宏博,张凯,涂佳黄(上海交通大学 船舶海洋与建筑工程学院,上海 ;上海交通大学 海洋工程国家重点实验室,上海 ;湘潭大学 土木工程与力学学院,湖南 湘潭 )摘要:研究大型水平轴风力机塔筒的合理支撑结构形式对风力发电结构系统安全至关重要针对高桩承台式水平轴风力机,基于反比例型删除率的双向渐进结构优化(B E S O)算法,拓扑优化风力机塔筒结构,根据计算流体动力学方法和桩土相互作用原理,建立考虑桩土作用的风力机结构风致动力响应数值模型,比较分析优化前后风力机结构动力响应特性,验证结构优化方法的可靠性经拓扑优化得到新的高桩承台水平轴风力机支撑结构,相较于初始结构,优化后的支撑结构风致动力响应显著降低,成果可为高桩承台式水平轴风力机支撑结构设计优化提供技术参考关键词:水平轴风力机;高桩承台式风力机基础;双向渐进结构优化算法;风致动力响应中图分类号:P ;O 文献标志码:AS u p p o r tS t r u c t u r eO p t i m i z a t i o no fH i g h P i l eC a pS u p p o r t e dH o r i z o n t a lA x i sW i n dT u r b i n eS y s t e mB a s e do nB E S OA l g o r i t h mZHANL i n g y u,HEW e n j u n,ZHO UD a i,HANZ h a o l o n g,ZHUH o n g b o,ZHANGK a i,T UJ i a h u a n g(S c h o o l o fN a v a lA r c h i t e c t u r e,O c e a na n dC i v i lE n g i n e e r i n g,S h a n g h a i J i a oT o n gU n i v e r s i t y,S h a n g h a i ,C h i n a;S t a t eK e yL a b o r a t o r yo fO c e a nE n g i n e e r i n g,S h a n g h a i J i a oT o n gU n i v e r s i t y,S h a n g h a i ,C h i n a;C o l l e g eo fC i v i lE n g i n e e r i n ga n dM e c h a n i c s,X i a n g t a nU n i v e r s i t y,X i a n g t a n ,H u n a n,C h i n a)A b s t r a c t:T h es t u d yo fr e l i a b l es u p p o r ts t r u c t u r ei so fg r e a ts i g n i f i c a n c et ot h es a f e t yo fl a r g e s c a l e dh o r i z o n t a l a x i sw i n dt u r b i n e(HAWT)s y s t e m I nt h i sp a p e r,f o rc a p s u p p o r t e dHAWT w i t hh i g hp i l e,t h eb i d i r e c t i o n a le v o l u t i o n a ls t r u c t u r eo p t i m i z a t i o n(B E S O)a l g o r i t h m b a s e do ni n v e r s e l yp r o p o r t i o n a ld e l e t i o nr a t ew a su s e dt oo p t i m i z e i t ss u p p o r ts t r u c t u r e U s i n gc o m p u t a t i o n a l f l u i dd y n a m i c s(C F D)a n dt h ep r i n c i p l eo fp i l e s o i l i n t e r a c t i o n,t h e f i n i t ee l e m e n tm o d e lo fHAWTw a se s t a b l i s h e d,w h e r et h ew i n dl o a da n dp i l e s o i l i n t e r a c t i o nw e r et a k e ni n t oc o n s i d e r a t i o n T h er e l i a b i l i t yo ft h es t r u c t u r a lo p t i m i z a t i o nm e t h o dw a sv e r i f i e d t h r o u g h t h e c o m p a r i s o no f t h ed y n a m i c r e s p o n s e c h a r a c t e r i s t i c sb e t w e e n t h e i n i t i a l a n dt h eo p t i m i z e dm o d e l T h er e s u l t ss h o wt h a tt h ec u r r e n tB E S Oa l g o r i t h mc a ne f f e c t i v e l yg e n e r a t ean o v e l第 卷 第期 年月上 海 交 通 大 学 学 报J OUR NA LO FS HANGHA I J I A OT ON GUN I V E R S I T YV o l N o A u g 上海交通大学学报第 卷s u p p o r ts t r u c t u r e f o r mf o rh i g h p i l eHAWT,w h o s ed y n a m i c r e s p o n s e i s s i g n i f i c a n t l y r e d u c e d T h e r e s u l t sc a np r o v i d eu s e f u l r e f e r e n c e s f o rHAWT sd e s i g n s K e yw o r d s:h o r i z o n t a la x i sw i n dt u r b i n e(HAWT);h i g h p i l ec a pf o u n d a t i o n;b i d i r e c t i o n a lp r o g r e s s i v es t r u c t u r a l o p t i m i z a t i o na l g o r i t h m(B E S O);w i n d i n d u c e dd y n a m i c r e s p o n s e风力机大型化可有效提高风力发电量和风能利用率,是海上风力发电的世界发展趋势研发建造包括、和 MW在内的大型风力机是我国海洋绿色能源开发利用的战略重点然而随着风力机转子、塔筒尺寸的不断增加,大型风力机结构系统柔性增加,抗风能力明显降低,风力机风致动力响应更加复杂,系统风致动力灾变风险更加凸显研究并优化大型水平轴风力机支撑结构对保证系统运行安全十分重要风力机大型化导致塔筒顶部大型转子和发电机质量增加、支撑塔筒高度增大常用的圆锥塔筒无法满足风力机结构的强度和刚度要求迄今,学者们对大型水平轴风力机支撑结构优化已有一定研究 G e n t i l s等基于遗传算法,考虑多 准 则 约 束,对 美 国 国 家 可 再 生 能 源 实 验 室(N R E L)提 出 的 MW单 桩 式 风 力 机(N R E LMW)塔架和基础进行优化分析并计算风力机支撑结构优化尺寸参数;O e s t等基于疲劳损伤荷载、准静态分析和动态分析,对海上风力机支撑结构进行梯度优化设计;葛旭等结合有限元分析与进化策略,对海上三桩单立柱风力机支撑结构主尺度参数进行多参数同步优化设计;T i a n等对导管架式海上风力机的导管架基础进行拓扑优化,结合尺寸和形状优化实现结构重分布,得到质量更轻的导管架结构上述研究仅为针对特定风力机结构形式的参数优化,未改变其支撑结构形式,虽对风力机结构稳定性有一定提高,但设计方法具有局限性相较于参数优化,拓扑优化方法具有更多的设计自由度,是优化大型水平轴支撑结构的新途径拓扑优化方法应用于包括结构优化在内的广泛领域,其双向渐进结构优化(B E S O)算法优化效果较好王章骏等采用B E S O算法对吸能管进行拓扑优化,提高了带隔板薄壁方管的耐撞性能,实现了吸能管轻量化的设计目标;刘昆等基于B E S O算法对船体梁腹板开孔进行设计优化,提出适用于船体结构孔型的B E S O算法;B E S O还在罗马小体育宫、卡塔尔国家会议中心 等建筑工程设计中得到应用目前,大型水平轴风力机塔筒支撑结构拓扑优化还鲜有研究基础是风力机系统的关键结构之一,高桩承台基础的结构安全性能高,适合我国海床地质条件,施工工 艺 成 熟,在 我 国 近 海 风 电 场 建 设 中 多 有 应用 但不考虑桩土相互作用会导致风力机位移响应被严重低估 ,是结构安全设计的重要隐患考虑桩土相互作用的风力机支撑结构性能,开展支撑结构优化分析基于反比例型删除率的B E S O算法,提出大型高桩承台式水平轴风力机支撑结构拓扑的优化算法,运用有限元软件AN S Y S开展大型水平轴风力机支撑结构寻优分析,获得一种风力机优化支撑结构,显著降低了风力机风致动力响应 风力机结构建模与分析 风力机支撑结构建模以应用较广的N R E LMW 风力机为背景,进行考虑桩土作用的近海坐底式高桩承台大型水平轴风力机塔筒结构优化研究高桩承台式水平轴风力机初始模型如图(a)所示,其支撑系统包括上部塔筒结构、承台和群桩水平轴风力机塔筒为圆筒式支撑结构,如图(b)所示,塔筒底部外径Dm,塔筒高度H m 塔筒选用Q 钢材,密度 g/m,弹性模量E N/m,泊松比 ,风力机额定风速为 m/s,具体参数如表所示图水平轴风力机塔筒结构初始模型F i g H o r i z o n t a l a x i sw i n dt u r b i n e i n i t i a lm o d e l第期占玲玉,等:基于B E S O算法的高桩承台式水平轴风力机支撑结构优化 表MW水平轴风力机塔架参数T a b MWh o r i z o n t a l a x i sw i n dt u r b i n ep a r a m e t e r s参数数值参数数值塔筒底部外径,D/m额定输出功率/MW塔筒顶部外径/m 叶片长度,L/m 塔架高度,H/m 额定风速/(ms)转子质量/k g 额定叶尖速度/(ms)机舱质量/k g 塔架质量/k g 图基础与土体模型F i g F o u n d a t i o ng e o m e t r i cm o d e l a n ds o i lp r o f i l e s图风场计算域和边界条件F i g W i n df i e l dc o m p u t a t i o nd o m a i na n db o u n d a r yc o n d i t i o n参照上海东海大桥海上大型风力发电工程项目采用高桩承台式 基 础,混 凝 土 承 台 直 径D m,承台厚度Sm,设置根直径dm、壁厚t mm的 全 直 钢 管 基 桩,桩 基 长 度H m,入土深度h m承台材料选用C 标号混凝土,弹性模量E G P a,密度为 k g/m,泊松比 ,基础模型如图(a)所示基于J u n g等 和C h e n等 确定土体模型和力学参数:土体浮重度、弹性模量E、泊松比、内摩擦角和不排水抗剪强度su、黏土泊松比u、的最大有效应力下的轴向应变、土层深度z、膨胀角、黏聚力c、地基反力模量k,如图(b)所示 风力机支撑结构风荷载对高桩承台大型水平轴风力机结构系统,基于计算流体动力学(C F D)方法并运用S t a r C CM计算软件 数值模拟获得风力机叶片气动荷载,风力机风轮直径为d,选取长方体风场计算域,计算域长度、宽度和高度分别为 d、d和d,如图(a)和(b)所示计算域划分共 万个网格计算域 上海交通大学学报第 卷中,风力机叶片圆柱体旋转域在竖平面的直径为 d,旋转域高度即顺风向尺度为 d;旋转域之外的计算域等其他部分为静止域对旋转域和静止域设置交界面,采用滑移网格处理风轮运行旋转问题,如图(c)和(d)所示叶片表面采用棱柱层网格,延伸因子为 ,总厚度为 m,无量纲化壁面距离为y对大型风力机支撑结构采用锥筒结构,塔筒高度为 m将塔筒简化为受弯曲变形、轴向压缩和扭转变形的悬臂梁风力机运行时,需考虑塔筒沿高度的风压变化沿塔筒高度l处的风荷载 为f(l)CdDlV(l)()式中:Cd为风阻力系数,取值 ;为空气密度;Dl为塔筒高度l处的塔筒截面直径;V(l)为沿塔架高度l处的瞬时风速 桩土相互作用美国石油协会(A P I)规范推荐的土体py曲线,即非线性弹簧模型广泛应用于海洋工程桩基设计引入用于分析风力机支撑结构的桩土相互作用,该模型考虑了土体非线性、土层变化和土体种类等因素的影响,软土中各层py曲线为ppuyyc/,yyc pu XXRyyc yc,ycy yc puXXR,y yc()式中:p为深度X处单位桩长的极限水平土抗力标准值;yc为实验室中对未扰动工试样做不排水压缩试验时,其应力达到最大应力一半时的变形;y为泥面以下深度X处桩的侧向水平变形;XR为泥面以下到土抗力减少区域底部的深度;pu为深度X处单位桩长的极限水平土抗力标准值对于砂土情况,py曲线为pA put a n hK XA puy()式中:K为地基反力初始模量;A为计入静力荷载和循环荷载条件的参数基于A P I规范的py曲线法,计算循环荷载作用情 况 下 的 桩 侧 土 抗 力,使 用AN S Y S软 件 的C o m b i n e 单元模拟py曲线,该单元是具有非线性功能的单向拉压单元,常用于桩土相互作用模拟在搭建有限元模型时,根据计算得到的土体py曲线参数,每间隔m设置一个C o m b i n e 土弹簧单元在桩土相互作用模拟过程中,为简化计算未考虑桩土循环弱化作用及桩端设置竖向简支约束 风力机支撑结构系统振动响应采用N e wm a r k 法 中的平均常加速度法对风力机支撑结构系统风致振动进行瞬态分析,其中控制N e wm a r k 法精度和稳定性的 参数取值为 ,采用瑞利阻尼模拟风力机结构阻尼,刚度阻尼系数取值为 风力机结构系统动力激励下的振动方程 为M u(t)C u(t)K u(t)FI(t)FC(t)FK(t)()式中:u(t)、u(t)、u(t)分别为水平轴风力机的加速度、速度和位移矩阵;M、C、K分别为大型水平轴风力机结构的质量、阻尼和刚度矩阵;FI(t)、FC(t)和FK(t)分别为风力机的惯性力、阻尼力和恢复力矩阵 大型水平轴风力机支撑拓扑优化 B E S O算法的数学模型针对大型水平轴风力机支撑结构优化问题,建立以结构静刚度最大为目标、以结构体积为约束条件的数学力学模型,表示为m i nCTCT(x)uTK u()s t nixiviV,xixm i n,()式中:CT为结构平均柔顺度;V为目标体积比;vi为第i个单元的体积;n为迭代次数;xi为设计变量;xm i n和分 别 代 表 某 一 单 元 在 设 计 域 内 缺 失 或存在 单元灵敏度过滤及更新为避免产生奇异单元刚度矩阵,算法以xm i n代替表示空单元,在该软删单元的B E S O算法 中所有单元的材料弹性模量可根据材料插值模型设置为E(xi)Exqi()式中:q为惩罚因子单元灵敏度表示增删单元对结构平均柔顺度的影响程度,单元灵敏度大小为单元灵敏度的相对排序,表示为iqCTxiuTiKiui,xixq m i nuTiKiui,xixm i n()式中:Ki和ui分别为单元刚度矩阵和单元位移向量为抑制优化过程中出现实体单元与空洞单元交第期占玲玉,等:基于B E S O算法的高桩承台式水平轴风力机支撑结构优化 替出现的棋盘格等数值不稳定现象,通过权重系数引入节点灵敏度节点概念:jMiwiei()式中:j为节点j的灵敏度;M为对节点j有影响的单元总数;wi为权重系数;ei为单元i的灵敏度在迭代过程中,根据节点灵敏度对单元灵敏度进行过滤并更新单元灵敏度,更新过程为niJj(rm i nri j)niJj(rm i nri j)()inin i()式中:ri j为j节点到i单元的中心的距离;rm i n为以i单元中心为原点的过滤圆半径;ni为第n次迭代的单元灵敏度,将式()计算得到的i赋予ni,进行下一次迭代 基于B E S O算法的结构拓扑优化鉴于水平轴风力机结构系统对称性,将支撑结构三维优化问题转化为二维平面优化,以风力机塔筒支撑形式和支撑点位置为优化对象将风力机塔筒主轴高度方向和塔筒半径方向围成的平面作为初始设计域,采用AN S Y S软件中的节点P l a n e 单元对计算模型进行网格划分设置底部为固定端约束,选取最不利气动荷载作为外部风荷载采用反比例型删除率优化模型 对风力机支撑结构进行拓扑优化:m a xVVVim i nm a xm i n()式中:为结构材料的删除率;m a x为给定的最大删除率,取 ;m i n为给定的最小删除率,取 ;Vi为优化过程中每一步迭代过程的体积比;V值取 在迭代过程中,根据删除率对体积比进行更新Vk Vk(),k,()式中:Vk为第k次迭代的体积比:Vk 为第k次迭代的体积比当迭代结果满足根据体积约束和收敛准则 时,终止迭代为提高迭代效率,设置收敛准则,表示当前迭代步之前 次迭代中平均柔顺度的相对变化量,当目标函数变化随着迭代变化足够小时,则判定迭代结束:c hi(Cki Ck i)iCki()式中:Ck为第k次迭代的平均柔顺度;为收敛因子,取值为 根据上述计算过程实现基于B E S O算法的水平轴风力机支撑结构优化,流程如图所示图拓扑优化流程图F i g F l o wc h a r to f t o p o l o g yo p t i m i z a t i o n 优化前后模型受力特性 基于B E S O算法的支撑结构优化运用B E S O双向渐进结构优化算法,对MW大型水平轴风力机塔筒结构(高度为 m)进行拓扑优化,如图所示为防止风力机塔筒顶部发生过大位移,通常对风力机塔筒受弯矩作用较大的底部区域进行加固,所选取的底部拓扑区域为宽度a m、长度b m所围成的平面,以此作为初始设计域,采用节点P l a n e 单元进行网格划分在设计域区域施加外部荷载,包括风力机机舱、转子质量等纵向荷载Fz,风力机叶片运动产生的推力、沿塔筒高度分布的梯度风荷载等在设计域顶部产生的荷载Fx y和弯矩Mx y底部边界条件为固定端约束,拓扑优化的结构目标体积分数为 图塔架结构底部受力及拓扑设计区域F i g T o w e rs t r u c t u r ef o r c ed i a g r a m a n dt o p o l o g yd e s i g na r e a 上海交通大学学报第 卷分析拓扑过程中体积比即单元删除量变化和平均柔顺度相对迭代次数的变化历程,如图所示拓扑迭代初期,拓扑单元较多,单元删除量较大,拓扑迭代速率较快随着迭代进行,单元体积比逐渐减小,迭代 次后,结构单元体积比开始趋近于目标体积约束条件,拓扑形态开始稳定,最终得到反比例型删除率的优化模型柔顺度量值为 Nm对初始设计域,通过上述二维拓扑找形得到带图塔筒支撑结构拓扑优化F i g T o p o l o g yo p t i m i z a t i o no f t o w e r s u p p o r t s t r u c t u r e s图本文结果及其分析对比F i g P o w e r c o m p a r i s o nb e t w e e nC F Dm o d e l a n dr e f e r e n c e图迭代过程结构平均柔顺度和体积分数变化F i g A v e r a g e f l e x i b i l i t ya n dv o l u m e f r a c t i o no fo p t i m i z e ds t r u c t u r e有侧向支撑的塔筒形状如图(a)所示,其侧向支撑高度为 m、宽度为m基于支撑结构的对称性,考虑基础为四桩布局,将二维拓扑结果扩展到三维空间,获得的四支撑结构优化形式如图(b)所示根据拓扑优化结果,构建高桩承台式水平轴风力机支撑结构模型 优化后风力机支撑结构特性基于C F D数值方法计算优化后风力机支撑结构的风致水平推力和风力机功率为考虑最不利荷载工况,在进行C F D数值计算时,设速度入口风速等于切出风速为 m/s,与边界元法(B EM)、致动盘(A D)及有限体积法(F VM)中的结果进行对比分析,如图所示可知,本文风力机功率、风力机旋转时叶片对塔筒的水平推力(T k N)均与相关文献结果接近分别计算优化后支撑结构和原塔筒结构的静力特性,注意到风力机塔筒经受周期性脉动风荷载作用,可将风力机旋转多圈并稳定后叶片对塔筒的水平推力用于计算最不利工况下的塔筒位移,运用式第期占玲玉,等:基于B E S O算法的高桩承台式水平轴风力机支撑结构优化 ()计算沿风机塔筒高度的风荷载风力机支撑结构顺风向(x向)、横风向(y向)和竖向(z向)最大动位移均发生于塔筒顶部,其中顺风向位移明显大于横风向位移针对 MW大型水平轴风力机,分别计算初始支撑结构和优化后支撑结构在额定风速下的塔顶位移初始结构塔顶顺风向位移为 mm,竖向位移为 mm,而优化后结构塔顶的顺风向位移为 mm,竖向位移为 mm优化后支撑结构的位移明显减少,顺风向和竖向位移相对初始支撑结构减少了 和,如图所示图初始结构与优化结构位移响应云图F i g C o m p a r i s o no fd i s p l a c e m e n to f i n i t i a l s t r u c t u r ea n do p t i m i z e ds t r u c t u r e图 初始和优化结构的塔顶顺风向速度时程和加速度时程F i g C o m p a r i s o no fv e l o c i t ya n da c c e l e r a t i o nt i m eh i s t o r yc u r v ea t t h e t o po f t o w e r 桩土作用优化后支撑结构风致动力响应分析考虑桩土作用,分析MW风力机初始和优化后结构在顺风向和竖向的动力响应设来流风速为 m/s,风力机叶片转速为 r a d/s,风力机运动周期为 s 风力机旋转 圈后气动性能基本稳定,相应的风力机旋转时长为 s 初始和优化后结构的塔筒顶端顺风向位移(Ux)时程曲线如图 所示计算发现在t s时刻,风力机结构振幅较大因此分析该时刻优化前后支撑结构的顺风向和竖向位移,得到初始和优化后塔顶顺风向最大位移分别为 和 mm,减少 ;竖向位移分别为 mm和 mm,最大竖向位移减少了 由此可见,在保持与初始结构体积相同情况下,优化后风力机支撑结构位移明显减少,结构刚度大于初始结构图 初始结构和优化结构的塔顶顺风向位移时程曲线F i g C o m p a r i s o no fd i s p l a c e m e n t t i m eh i s t o r yc u r v eo ft w om o d e l sa t t h e t o po f t o w e r考虑桩土作用的初始和优化后结构在 s时段顺风向的速度(vx)响应时程和加速度(ax)响应时程如图 所示可见优化后结构的速度和加速度 上海交通大学学报第 卷响应小于初始结构,风振特性得到显著改善,在风力机旋转多圈如 s进入稳态后改善效果更明显,如塔顶加速度由初始结构的 m/s减小为优化后结构的 m/s,减幅达 考虑桩土作用的支撑结构振动效应明显减小,具有较好的减振效果结论针对大型高桩承台式水平轴风力机支撑结构系统,运用C F D方法计算风致效应,考虑桩土作用效应,运用基于反比例函数的双向渐进结构优化算法对风力机塔筒支撑结构进行拓扑优化,开展风力机结构系统的静动力响应特性分析,得到结论如下:()运用基于反比例型删除率的B E S O算法并考虑桩土作用影响,开展水平轴风力机塔筒支撑结构拓扑优化,得到风力机支撑结构的新形式()优化后支撑结构风致效应明显减小,在保持结构体积不变情况下,优化后支撑结构整体刚度明显提高,塔顶顺风向和竖向静力位移分别减小 和,动力分析下塔顶在顺风向和竖向最大位移分别减少 和 运用B E S O算法进行水平轴风力机支撑结构的二维平面优化,后续将对风力机支撑结构开展三维优化研究参考文献:A GU I L A R VA R GA SS,T E L L E SE S T E V E SG R,MA A I R AP M,e ta lW i n dp o w e rg e n e r a t i o n:Ar e v i e wa n dar e s e a r c ha g e n d aJ J o u r n a lo fC l e a n e rP r o d u c t i o n,:KA EWN I AMP,C A O MS,A L KAY EM NF,e t a lR e c e n t a d v a n c e s i nd a m a g ed e t e c t i o no fw i n dt u r b i n eb l a d e s:As t a t e o f t h e a r t r e v i e wJ R e n e w a b l ea n dS u s t a i n a b l eE n e r g yR e v i e w s,:孟珣基于动力特性的海上风力发电支撑结构优化技术研究 D青岛:中国海洋大学,ME N GX u n O p t i m u mt e c h n o l o g yo ns u p p o r ts t r u c t u r e s o f o f f s h o r ew i n d t u r b i n eb a s e do nd y n a m i cp r o p e r t i e sDQ i n g d a o:O c e a n U n i v e r s i t yo fC h i n a,G E N T I L S T,WAN G L,KO L I O S AI n t e g r a t e ds t r u c t u r a l o p t i m i s a t i o no fo f f s h o r ew i n dt u r b i n es u p p o r ts t r u c t u r e sb a s e do nf i n i t ee l e m e n ta n a l y s i sa n dg e n e t i ca l g o r i t h mJA p p l i e dE n e r g y,:O E S TJ,S AN D A L K,S CHA F H I R TS,e ta lO ng r a d i e n t b a s e do p t i m i z a t i o no fj a c k e ts t r u c t u r e sf o ro f f s h o r ew i n dt u r b i n e sJW i n dE n e r g y,():葛旭,徐业鹏,黄丹基于进化策略的海上风电支撑结构多参数同步优化设计 J可再生能源,():G EX u,X U Y e p e n g,HUAN G D a n M u l t i p a r a m e t e r s i m u l t a n e o u s l yo p t i m a l d e s i g no f s u p p o r t s t r u c t u r ef o ro f f s h o r ew i n dt u r b i n eb a s e do ne v o l u t i o n a r ys t r a t e g yJ R e n e w a b l eE n e r g yR e s o u r c e s,():T I ANXJ,S UNXY,L I UGJ,e t a lO p t i m i z a t i o nd e s i g no ft h ej a c k e ts u p p o r ts t r u c t u r ef o ro f f s h o r ew i n d t u r b i n eu s i n g t o p o l o g yo p t i m i z a t i o nm e t h o dJO c e a nE n g i n e e r i n g,:王章骏,许平,邢杰,等基于改进B E S O算法的带隔板薄壁方管耐撞性拓扑优化 J铁道科学与工程学报,():WAN GZ h a n g j u n,X UP i n g,X I NGJ i e,e t a l C r a s h w o r t h i n e s st o p o l o g yo p t i m i z a t i o no ft h et h i n w a l l e ds q u a r et u b e w i t h d i a p h r a g m s b a s e d o n i m p r o v e dB E S Oa l g o r i t h mJ J o u r n a lo fR a i l w a yS c i e n c ea n dE n g i n e e r i n g,():刘昆,郭德松,王仁华,等基于改进B E S O方法的多工 况 船 体 开 孔 孔 形 优 化 J中 国 海 洋 平 台,():L I U K u n,GUO D e s o n g,WAN G R e n h u a,e ta lS h a p eo p t i m i z a t i o no fh u l l s t r u c t u r eh o l eo p e n i n gu n d e r m u l t i p l ec o n d i t i o n s b a s e d o ni m p r o v e d B E S Om e t h o dJ C h i n aO f f s h o r eP l a t f o r m,():YANX,B A OD W,Z HOU YF,e t a l D e t a i l c o n t r o ls t r a t e g i e sf o rt o p o l o g yo p t i m i z a t i o ni na r c h i t e c t u r a l d e s i g na n dd e v e l o p m e n tJ F r o n t i e r s o fA r c h i t e c t u r a lR e s e a r c h,():周豪,江滨矶崎新 引领世界建筑创新潮流的先锋J中国勘察设计,():Z HOU H a o,J I AN GB i n A r a t aI s o z a k i i sap i o n e e ri nt h ew o r l do f a r c h i t e c t u r a l i n n o v a t i o nJ C h i n aE n g i n e e r i n g&C o n s u l t i n g,():孙明明,李昕,李炜海上风力机高桩承台基础反应特性研究 J太阳能学报,():S UN M i n g m i n g,L IX i n,L IW e i S t u d yo nr e s p o n s eo fo f f s h o r ew i n d t u r b i n ew i t hh i g h p i l e c a p f o u n d a t i o nJ A c t aE n e r g i a eS o l a r i sS i n i c a,():F UT A IM M,HA I GH SK,MA D A B HU S H IGSP C o m p a r i s o no f t h ed y n a m i c r e s p o n s e s o fm o n o p i l e sa n dg r a v i t yb a s ef o u n d a t i o n sf o ro f f s h o r ew i n dt u r b i n e si ns a n du s i n gc e n t r i f u g e m o d e l l i n gJS o i l sa n dF o u n d a t i o n s,():A B H I NAVKA,S AHAN N o n l i n e a rd y n a m i c a l b e h a v i o u ro f j a c k e t s u p p o r t e do f f s h o r ew i n dt u r b i n e s i nl o o s es a n dJM a r i n eS t r u c t u r e s,:第期占玲玉,等:基于B E S O算法的高桩承台式水平轴风力机支撑结构优化 刘晨晨,张琪,李明广,等波浪与地震荷载共同作用下桩的动力响应 J上海交通大学学报,():L I UC h e n c h e n,Z HAN G Q i,L IM i n g g u a n g,e ta lD y n a m i cr e s p o n s eo fp i l ea tw a t e r w a v e l o a da n ds e i s m i c l o a dJ J o u r n a l o fS h a n g h a i J i a oT o n gU n i v e r s i t y,():J ONKMANJ,B UT T E R F I E L DS,MU S I A L W,e ta lD e f i n i t i o no fa MWr e f e r e n c ew i n dt u r b i n ef o ro f f s h