钝体表面附着物对低速水流压电俘能器性能影响研究.pdf

55(7):1463-1472piezoelectric energy harvester performance in low velocity water flow.Chinese Journal of Theoretical and Applied Mechanies,2023.Sun Weipeng,Liu Chenhan,Yu Xiaobin,Hu Shen,ZhongKexin,ZhaoDaoli.Effect of attachment forbluff body surfaceon1463-1472引用格式：孙维鹏，刘宸涵，郁小彬，胡坤，钟可欣，赵道利.钝体表面附着物对低速水流压电俘2023,55(7):流体力学Chinese Journal of Theoretical and Applied MechanicsJul.,20232023年7 月Vol.55,No.7力第55卷第7 期报学学钝体表面附着物对低速水流压电俘能器性能影响研究孙维鹏*刘宸涵*郁小彬胡坤*钅钟可欣*赵道利*,2)*（西安理工大学省部共建西北旱区生态水利国家重点实验室，西安7 1 0 0 48)+（华东天荒坪抽水蓄能有限责任公司，杭州31 0 0 1 2)摘要低速水流能可作为可再生能源的重要补充，适用于低速水流能开发的压电俘能器已在流场中表现出了优异的性能.而在钝体的表面增加附着物可影响压电俘能器的能量转换.通过改变附着物的形状、凸起高度和凹陷深度探究其对驰振式压电俘能器输出特性的影响.利用拓展的哈密顿原理建立俘能器的机械控制方程，引入高斯定理建立电场一位移控制方程，并基于准稳态假设计算驰振的水动力及力矩，进而得到压电俘能器机电耦合分段参数模型.用伽辽金法离散悬臂梁的位移，并在此基础上解耦控制方程，从而得到输出功率及悬臂梁位移近似解.通过水槽实验获得压电俘能器输出功率的实验值，并验证了数学模型的准确性.结果表明：0.51 m/s流速下，-2 mm椭球型凸起的压电俘能器RMS功率为1.41 1 mW，与无附着物的椭圆柱相比增幅为6 9.8 8%；当钝体包裹6 mm凸起的圆柱型附着物时，悬臂梁末端的振幅为3.0 7 mm，相比于无附着物的情况降低了8 4.8 3%.三维数值模拟的结果表明凸起、凹陷为2 mm附着物的应用会使在钝体上形成的压差升高，进而加剧流致振动，且会影响从钝体两侧脱落涡的大小及强度.关键词低低速水流能，附着物，压电俘能器，分段参数模型，驰振中图分类号：TN712+.5文献标识码：Adoi:10.6052/0459-1879-23-065EFFECTOFATTACHMENTFORBLUFFBODYSURFACEONPIEZOELECTRICENERGYHARVESTERPERFORMANCEINLOWVELOCITYWATERFLOWI)Sun Weipeng*JLiu Chenhan*Yu Xiaobin tHu Shen*Zhong Kexin*Zhao Daoli*,2)(State Key Laboratory of Eco-hydraulics in Northwest Arid Region,Xian University of Technology,Xian 710048,China)(East China Tianhuangping Pumped Storage Power,Co.,Ltd.,Hangzhou 310012,China)AbstractLow velocity water flow energy is an important complement to renewable energy and the piezoelectricenergy harvesters developed by using this energy have shown excellent performance in flow fields.The energyconversion of piezoelectric energy harvesters can be affected by the added attachments on bluff body surface.Therefore,the output characteristics of the piezoelectric energy harvester are investigated by varying the shape,theconvexity height and the concavity depth of the attachment.Mechanical governing equation of the piezoelectric energyharvester is developed by using the extended Hamiltons principle,the electric field-displacement governing equation is2023-03-01收稿，2 0 2 3-0 6-0 5录用，2 0 2 3-0 6-0 6 网络版发表.1）国家自然科学基金资助项目（52 1 7 90 8 9,52 2 0 91 1 5）.2)通讯作者：赵道利，教授，主要研究方向为压电能量采集、水轮机结构分析.E-mail:力14642023年第55卷报学学obtained by applying Gauss law,the quasi-steady-state assumption is employed to calculate the hydrodynamic force andmoment of galloping,and then the electromechanical coupling piecewise distributed parameter model of thepiezoelectric energy harvester is established.The displacement of the cantilever beam is discrete by Galerkin procedure,and the governing equations are decoupled on this basis to obtain the approximate solution of harvester output powerand bluff body displacement.Experimental solutions of the output power are obtained by water channel experiment andthe accuracy of model is further verified.The results showed that at the flow velocity of 0.51 m/s,the RMS power of thepiezoelectric energy harvester wrapped ellipsoidal attachment with-2 mm convexity is 1.411 mW,an increase of69.88%compared to the elliptic cylinder without attachment,and the amplitude of cantilever beam is 3.07 mm when thebluff body is equipped on the cylindrical attachment with 6 mm convexity,a reduction of 84.83%compared to the onewithout the attachment.The numerical simulation results show that the application of the convex and concaveattachment for 2 mm increases the pressure differential across bluff body,which intensifies vibration,it also affects thesize and strength of shedding vortices from both sides of the bluff body.Key words low velocity water flow energy,attachment,piezoelectric energy harvester,piecewise distributed parametermodel,galloping引言近年来，能源危机愈演愈烈，为了应对新环境下全球能源格局发生的深刻调整2 ,高效开发利用可再生能源迫在眉睫3-6 .水能作为占比最高的可再生能源，如何高效利用一直备受关注.天然河道中的高速水流能可通过水轮发电机组加以利用7-8 ,这一技术目前已相对成熟；而湖泊与小型河流中的低速水流能虽可作为新能源的重要补充，但因其分布广、能量低而存在利用困难的问题.驰振9 作为一种能够发生在低速水流中的现象，有流速响应范围宽和无锁定限制的特点.因此基于驰振设计且已在流场中展现出游优异性能的驰振式压电俘能器为低速水流能的利用提供了解决方案1 0 .利用驰振式压电俘能器将低速水流能量转化为电能，可以为微机电系统1-1 2 和无线传感器1 3-1 4 供电.目前已经存在一些关于利用低速水流能的压电俘能器研究.曹东兴等1 5 设计一种利用永磁铁增强流致振动的压电俘能器.实验结果表明：水流速度为0.5m/s时，输出功率的均方根值为1 2 0 W，与无附加磁体的压电俘能器相比提高了57.8%.单小彪等1 6 提出一种在水中完全淹没的尾流激振式压电俘能器.结果表明，当钝体的直径为1 5mm时，在上游距其中心2 0 mm处放置直径5mm扰流柱能够使压电俘能器输出2 3.35V的电压，与没有扰流柱相比输出电压提升了51%.Liu等1 7 基于驰振的触发方式，提出一种新的时间变流型方法.该方法解释了传统的时间平均流型方法难以解释的横向驰振机制.Lian等1 8 研究钝体截面为正方形、三角形和圆形的压电俘能器在水中的性能.实验表明：截面为三角形的钝体，其振幅比大于1.8；圆柱型钝体的能量转化率最高为8 0%.燕翔等1 9 实验研究T型截面钝体的驰振特性，并进一步探究其在不同外力条件及不同截面宽高比和阻尼比条件下的振动响应差异。He等2 0 通过研究证实改变钝体的几何形状，可以实现耦合或独立的驰振和涡激振动.Zhao等2 1 采用圆柱形扰流柱，对钝体截面形状为圆形、D形和倒D形的压电俘能器进行水槽实验测试.结果表明，3种压电俘能器的采集功率分别达到57 0.3，59 6.4和1 0 7 4W/m3，是无尾流激励情况的43.2，2 5.3和31倍.张军等2 2 对三棱柱的流致振动进行了实验研究,并评估其能量转化效率.Song等2 3 研究圆柱型钝体压电俘能器在不同直径、质量和阻尼比条件下的输出性能.研究证实在一定范围内，降低质量和增加钝体直径都可以提高压电俘能器的输出.以上研究主要通过改变钝体的形状来优化压电俘能器，而在钝体表面添加小型附着物的研究还很少.存在于钝体表面的附着物能够影响压电俘能器系统的振动与输出2 4-2 5,目前已存在一些将附着物与压电俘能器相结合的应用：Zhu等2 6 提出一种由圆柱体和五角星组成的新型压电俘能器，双向流固耦合仿真的结果表明其功率密度可达到8 8 5W/m3.Hu等2 7 通过调整圆形、三角形和矩形杆件在圆柱形钝体上的位置来提升设备的风能量采集效率.实验证实当安装的三棱柱夹角为6 0 时，能量采集器的性能最佳.Zhu等2 8 设计了在圆柱体表面安装对称鳍状条带的新钝体，数值模拟的结果表明收集到的最大功率超过6 0 W/m.目前关于附着物在压电俘能器1465第7 期孙维鹏等：钝体表面附着物对低速水流压电俘能器性能影响研究上的研究都只局限于其形状及安装形式，缺少在钝体上进行切除特定形状的尝试.附着物的凸起高度和凹陷深度值对压电俘能器性能的影响有待研究.综上所述，本文主要研究不同凸起高度/凹陷深度和形状的附着物对驰振式压电俘能器输出特性的影响.引入水的附加质量力改进现有的分段参数模型.基于准稳态假设，建立驰振的水动力模型并计算出驰振力.水槽实验验证了数学模型的准确性，进一步通过模型分析压电俘能器的输出功率和悬臂梁振幅与凸起高度/凹陷深度之间的关系.最后基于数值模拟从流体动力学的角度阐述超表面对系统影响的根本原因.1数学模型本节利用拓展的哈密顿原理推导出压电俘能器的机械控制方程，引入高斯定理建立电场一位移控制方程，并基于准稳态假设计算驰振的水动力及力矩，利用伽辽金法离散悬臂梁的位移，并在此基础上解耦控制方程，得出输出功率和悬臂梁位移的近似解，1.1俘能器机械控制方程本文研究的压电俘能器由钝体、悬臂梁和MFC压电片组成，如图1所示.悬臂梁的上端固定，中部黏贴有压电片，下端与钝体连接.x轴为悬臂梁的长度方向，轴为水流方向，受到水流冲击的钝体在z方向上发生振动，考虑到振幅较大，选择脆性低、韧性高的宏观纤维复合材料(macrofibercomposites,fixedMFCpiezoelectricsheetcantileverbeambluff bodydirectionofwaterflowvibratingdirectionof mass block图1压电俘能器Fig.1PiezoelectricenergyharvesterMFC)作为压电材料.椭圆柱钝体的表面装配了人工设计的附着物，其形状采用椭球型和圆柱型两种，各自的最大高度被定义为特征值，分别为Eell和Ecy，如图2 所示.椭球型和圆柱型附着物的特征值均设置为6，4，2,-1和-2 mm,其中正值为凸起，负值为凹陷.基于拓展的哈密顿原理 2 9-3 0 ,压电俘能器的控制方程可表示为04v(x,t)0v(x,t)05v(x,t)EI+Ca+Cs1十0 x40 xox4ot22v(x,t)ds(x-l1)ds(x-2)m十90t2dxdxds(x-l)Ftipo(x-l)-Mip(1)dxEel=4 mmE.=-2mm(a)Ee=4,-2 mm 的椭球型附着物结构(a)Structure of ellipsoidal attachments with Ee=4,-2 mmE=2mmCyEey=-2 mm(b)Ecy=2,-2 mm的圆柱型附着物结构(b)Structure of cylindrical attachments with Eey=2,-2 mm图2 附着物的结构图Fig.2Structural diagram of the attachments力1466报学2023年第55卷学式中，v(x,t)表示t时刻悬臂梁x处的位置在z方向上的位移；c.和c，分别表示悬臂梁的气动阻尼系数和黏性应变系数；表示克罗内克符号，其值在x=1时取1,反之取0;9 p为压电耦合系数，9 p=-e31bp(z1+z2)/2,e31为压电应力系数,e31=E,d31,d31为压电层的应变系数；V(t)压电材料的电压;Fip和Mip分别表示悬臂梁的驰振力和力矩.1.2电场一位移控制方程引入高斯定理 3 1ddVDndAC3dA(2)二didtRAA式中，D表示电位移矢量；n表示梁平面法向量；D3是电位移分量,表示为D(s,t)=e31e11(s,t)+833E3,其中33表示恒定应变下的介电常数，8 11是压电层的平均应变分量，表示为1Cz2 0v(x,t)11(x,z,t)=zdzZ2-Z1Jz10.x2系统电场一位移控制方程可表示为V()V)12(s,t)(3)ds0tos2PdtR式中C,是压电俘能器的电容,表示为3b,(2-1)Cp=hp1.3驰振的水动力计算为了建立驰振的水动力模型，引入准稳态假设.而假设中的驰振力会受到攻角的影响.攻角可表示为(l,t)+xi(l,t)Q=tan(4)U而攻角与驰振系数C,存在如下关系C,=ai tan +a tan(5)式中a和a3可通过实验获得.考虑到附加质量力对驰振力的影响，钝体单位长度的驰振力可表示为 3 2 wU2btipCy元pwbipCMF,=i(l,t)+xi(l,t)(6)24式中CM表示附加质量系数.悬臂梁所受到的驰振水动力和力矩可由下式计算Fip=(7)Mtif01.4驰振模型的机电解耦应用伽辽金法 3 3-3 4 离散悬臂梁的位移v(x,t),将其划分为空间变量和时间变量8v(x,t)=(8)i=1式中,Q(t)和i(t)分别为悬臂梁的第i阶模态坐标和模态.驰振一阶模态的水动力可以根据)=d(D)Fip+d()Mtip来计算.依据式（7)，得出一阶模态的水动力(）f(t)=k1+k3(9)k2OU式中,ki,kz和k3分别为a1puU2bipk1=2tipCMPu元b22()lip+$()s(lptipk2=+tip43a3puU2bipk3=0()+x0(D)dx+2x(0(1)+x0(0)dx)Abdelkefi等 3 5 的研究证实，悬臂梁的振动主要受一阶模态支配,Bibo等 3 6 也发现一阶模态满足被捕获的条件.因此本文采用一阶模态来简化悬臂梁的机电耦合模型Q(t)+2gwQ(t)+wQ(t)+0,V(t)=f(t)V()(10)C,V(t)+=0pQ(t)R式中，（表示阻尼系数；表示压电俘能器的一阶固有频率；，是压电耦合项，表示为p=9pb(l2)-(l);(t)表示一阶振型的水动力.解耦的控制方程 3 7 由式(9)和式(10)得250+Ce0()+0(0)=1+k21a?Q(t)+k1+k31+k2(1i)式中，和C。分别为压电俘能器的修正频率和阻尼,表示为1467孙维鹏等：钝体表面附着物对低速水流压电俘能器性能影响研究第7 期w2+Cp?Ce1+k2(12)RC.=CRa2+1进一步可得位移模态坐标与电压之间的关系Vo=CeRu?Qg(13)模态坐标的振幅为4(x-26w-Ce)Q0(14)-3ya2式中，x和分别为线性和立方项的水动力阻尼，并且表示为X=岁，=艺k1.-k3T3由式(14)确定模态振幅Q。后,压电俘能器的采集功率可表示为P=(15)R2实验设置本研究的实验部分在图3 所示的U型水槽中进行.水从地下水库通过泵与进水管被引入矩形薄壁堰中，薄壁堰上游的蜂窝网可以降低水的流，从而减少水流不稳定对上游水位探针读数的影响；水流离开水堰流经第一段长直段，在这里被充分稳定以满足实验段对稳定水流的需求；随后水流流入实验段，实验段的中部安装了压电俘能器；最后水流由尾水口流回地下水库.实验段的长度为6.2 m，其断面宽度为0.4m.实验中的流速是通过泵阀门及尾水口上方的电动阀门控制的，转动泵阀门能够调节进入水槽内的流量；调整电动阀门的高度可控制水槽内的水位，然后利用流量计算公式得出流速压电俘能器的钝体、悬臂梁及MFC压电片的参数如表1所示.本实验的流速范围为0.2 6 waterlevelgaugeNo.1inlet pipepiezoelectric energypumpvalveharvesterwaterlevelgaugeNo.2testingsystem图3 U型水槽Fig.3U-shape water channel表1钝体、悬臂梁和压电片的物理参数Table 1Physical parameters of bluff body,cantilever beam andpiezoelectric sheetParameterValuelength,width and thickness of cantilever beam/mm100,20,0.8Youngs modulus of cantilever beam/GPa110density of cantilever beam/(gcm-3)8.9length,width and thickness of MFC/mm28,14,0.3Youngs modulus of MFC/GPa33.336density of MFC/(gcm-3)5.44strain coefficient of MFC/(pC-N-l)-170dielectric constant of MFC/(nFml)14.04length and width of bluff body/mm100,50mass of bluff body/g200height of attachment/mm6,4,2,-1,20.51m/s，取0.0 5m/s为一流速步长.实验中利用示波器测试出与压电片串联的负载电阻的电压并将数据储存在计算机，然后计算出电压的均方根值(RMSvoltage).实验的负载电阻值为3.510 52.3模型验证为了验证所建立数学模型的准确性，分别绘制0.41,0.46和0.51m/s流速下椭球型和圆柱型附着物在2 mm凹陷条件下电压与预测值的对照，如图4所示.图中的实线为实验测得的电压曲线，是从单次exp-modexp-modexp-mod40200-20-400.41m/s-0.46m/s0.51m/s(a)凸起的2 mm椭球型附着物(a)Convex 2 mm ellipsoidall attachmentexp-modexp-modexp-mod40200-20-400.41 m/s-0.46m/s0.51m/s(b）凸起的2 mm圆柱型附着物(b)Convex 2 mm cylindrical attachment图4模型计算和实验结果的时域曲线对比Fig.4Comparison of time histories of the model calculations andexperimental results力1468报学2023年第55卷学140s的实验电压中选取相对稳定的4s绘制的；虚线为电压的数学模型计算值.图中实线和虚线基本重合，二者之间的微小差异是由水流的轻微波动引起的.所建立数学模型的计算结果与实验结果十分吻合，由此验证了数学模型的准确性.4压电俘能器输出特性和悬臂梁振幅的分析结合实验值和数学模型的近似解，分析椭球型和圆柱型附着物在不同的特征值和流速条件下对驰振式压电俘能器RMS输出功率和钝体振幅的影响.4.1流速对压电俘能器输出特性影响的分析椭球型附着物的压电俘能器RMS输出功率随流速的变化如图5(a)所示，可以看出在所有特征值下，压电俘能器的输出均随着流速的增加而增加.这是由于钝体的振幅随流速而增大，从而带动压电片发生更大程度的形变进而产生更高的电压，符合驰振的特性.此外附着物的特征值对压电俘能器输出1.8E=01.5Eel=6 mm-1 mm4 mm-2 mm1.22 mm0.90.60.300.260.310.360.410.460.51velocity/(ms-l)(a)椭球型附着物(a)Ellipsoidal attachment1.6E=0Eey=6 mm-1 mm1.24 mm-2 mm2mm0.80.400.260.310.360.410.460.51velocity/(ms-l)(b)圆柱型附着物(b)Cylindrical attachment图5流速对压电俘能器RMS输出功率的影响Fig.5Effect of flow velocity on the RMS output power of thepiezoelectricenergy harvester功率的影响也不同.在U=0.51m/s的流速下,Eel=-2mm的压电俘能器RMS输出功率为1.411mW相比于椭圆柱增幅为6 9.8 8%；而部分特征值的附着物则抑制压电俘能器的输出：当Eell=6,4mm时，压电俘能器的RMS输出功率在各个流速下始终低于未包裹附着物的椭圆柱，其中U=0.51m/s时，6mm特征值的椭球型附着物使输出功率降低了80.02%.图5(b)展示了包裹圆柱型附着物对压电俘能器性能的影响.与椭球型附着物相比，圆柱型附着物的施加使得压电俘能器的出输出整体偏低.当Eey=4,6mm时，压电俘能器在所有流速下的输出均被大幅度抑制，其中Ecy=6mm的压电俘能器在U=0.51m/s的输出为0.0 59 mW，降低了9 2.9%.此外随着流速的不断增加,Eey=-1,-2mm在流速超过0.41m/s后能够提高输出功率,其中Ecy=-2mm在0.51m/s的流速下能够产生1.3 0 8 mW的RMS功率，4.2特征值对压电俘能器输出特性的影响4.1中比较了不同特征值下流速对于压电俘能器RMS输出功率的影响，接下来结合数学模型的计算结果，进一步分析特征值的变化对于压电俘能器性能的影响.图6(a)表明随着特征值从-2 mm逐渐增大到6 mm,椭球型附着物压电俘能器的输出先减小再增大后减小的趋势，最大值分别在Eel1=-2,2mm处.这说明当特征值绝对值在0 附近，即凸起或凹陷相对较小时，附着物的施加能够增强钝体的振动，提高压电俘能器的输出功率，并且最大值的出现证明实现这种增强效果的附着物存在“最优特征值”.在0.51m/s流速下，未包裹附着物的原始钝体对应的压电俘能器RMS输出功率为0.8 3 mW,而分别应用和-2 mm特征值的椭球型与圆柱型附着物均有最大输出，分别对应于1.41和1.3 1mW，相比于原始俘能器分别提高了6 9.8 8%和57.8 3%.此外当特征值大于2 mm后，压电俘能器的输出产生不同程度的降低，这表明较大高度的附着物会抑制钝体的振动.对于圆柱型附着物如图6(b)所示，随着附着物高度从最大凹陷增加到最大凸起，压电俘能器的RMS输出功率在Ecy0时又呈现先增大后减小的趋势.这与椭球型附着物的趋势基本一致，说1469期孙维鹏等：钝体表面附着物对低速水流压电俘能器性能影响研究第2.4.molexp0.51 m/s2.0-.molexp0.46m/s.-molexp0.41 m/s1.6molexp0.36m/s1.20.8咖0.4-2-10246eigenvalue/mm(a)椭球型附着物(a)Ellipsoidal attachment2.0molexp0.51 m/smolexp0.46 m/s1.6molexp0.41 m/smolexp0.36m/s1.20.80.4012-10246eigenvalue/mm(b)圆柱型附着物(b)Cylindrical attachment图6 不同形状附着物的压电俘能器RMS输出功率随特征值的变化Fig.6Variation of RMS output power for energy harvesters versuseigenvalues with different shaped attachment applied明两种形状附着物的特征值变化对压电俘能器输出的影响规律是基本相同的.此外，附着物的恰当应用能够增加压电俘能器的输出.在椭球型附着物特征值为-2 mm、流速为0.51m/s的条件下，俘能器的最大RMS输出功率为1.411mW，对应的输出功率密度p可被计算为11998.3W/m3，与部分研究的俘能器输出功率密度对比如表2 所示表2 与同类研究对比Table 2Comparison with similar studiesAutherFlowvelocity/(ms-l)Maximumpp/(W-m-3)Zhao et al,210.651074Song et al,380.3175.9Cao et al,150.51142.86Song et al,230.35301.75Cao et al,390.53295.24ourresearch0.5111998.294.3特征值对悬臂梁振幅影响的分析为了进一步分析附着物特征值的变化对压电俘能器输出的影响，通过机电解耦得到悬臂梁的底端振幅，如图7 所示.通过与图6 中相应特征值条件下RMS输出功率的对比可以发现，悬臂梁的振幅与压电俘能器输出功率的变化趋势基本相同，二者呈正相关.此外当U=0.51m/s时，没有施加附着物的悬臂梁振幅为2 0.2 5mm相同流速下钝体包裹了Eey=6mm附着物的悬臂梁振幅则为3.0 7 mm，降低了84.83%，此时悬臂梁几乎静止,钝体的振动得到大幅度抑制.40.0.51m/s.0.46m/s0.41 m/s（0.36m/s3020100一2-10246eigenvalue/mm图7 圆柱型附着物悬臂梁的振幅随特征值的变化Fig.7Variation of amplitude with eigenvalue of cantilever beam forcylindrical attachment4.4包裹附着物的钝体扰流数值模拟适当特征值的附着物对压电俘能器功率输出具有增幅效应，并且无论是椭球型还是圆柱型附着物，凹陷或者凸起高度为2 mm时，增幅效应均最大.因此选取 Eell=-2 mm,Eey=-2 mm,Eell=2 mm,Eey=2mm和E=0的5个钝体进行三维数值模拟研究.D为钝体的迎水面宽度，计算域设置为2 0 D6D2D，钝体距上游面6 D.围绕钝体的周围设置加密的四面体网格，其他区域采用六面体网格.对网格进行无关性验证，均衡了效率和精度，最终网格数目控制在6 0 万左右，整体网格质量控制在0.5以上.为了定性地研究附着物对钝体流致振动的影响，并考虑到计算效率及收敛性，将整个计算域的网格设为固定网格，端流模型采用SSTk-w，设置0.51m/s流速入口和环境出口，模拟时间步长为0.0 0 12 5s，计算步数为8 0 0 0 步.图8 绘制了5个钝体的压力和涡量云图，从左侧流入的水流在受到钝体阻碍后流速降低并在其上游形成高压区；随后从侧边绕过加速流向下游并产力1470报学学2023年第55卷Eel=-2 mmEe,=-2 mmvortexcorepressure/Paregion/s-11.85010241017.650101-3.2001012101-1.405102El=2 mm-2.4901023.5751020-4.6601025.745102-2101-6.830102-7.9151029.000102Eey=2 mm-4.101E=0(a)压力云图(b)三维涡量云图(a)Pressurecontour(b)Three-dimensional vortical contour图8 U=0.51m/s时,Eell=2,-2mm,Eey=2,-2mm和E=0的云图Fig.8Contour of Eel=2,-2 mm,Ecy=2,-2 mm and E=0 with U=0.51 m/s生低压区，加速过程中水流受黏性的影响形成脱落涡.对比图8(a)所示的压力云图发现：包裹附着物的钝体周围产生的高、低压区之间的压力差明显大于未包裹附着物情况，这也对应更大的升力系数，如表3所示.图8(b)所示的三维涡量云图也表明包裹附着物后，从钝体两侧脱落的尾涡明显更大.绝对值为2mm附着物的存在，使在钝体上形成的压差升高，加剧流致振动，最终提高了压电俘能器的输出，这与试验及模型得出的结果是一致的.此外进一步对比包裹两种类型附着物钝体的脱落尾涡，发现在相同特征值下，包裹椭球型附着物的钝体尾流涡街要比包裹圆柱型附着物的情况略微紊乱，推测是因为附着物的结构干扰了流体分层，这种表3 升力系数Table3LiftcoefficientBluff bodyLiftcoefficientEel=-2 mm1.125Eey=-2 mm1.082Eell=2 mm1.037Eey=2 mm0.976E=00.716干扰作用在椭球型附着物上要大于圆柱型附着物5总结本文采用数学模型、实验及数值模拟相结合的方式研究了附着物的形状和附着物的特征值对驰振式压电俘能器输出的影响.基于拓展的哈密顿原理建立压电俘能器的机械控制方程，利用高斯定理建立电场-位移控制方程，通过准稳态假设计算出驰振力，在利用伽辽金法离散悬臂梁位移的基础上，实现对机电耦合控制模型的化简和解耦，最终得到系统输出及悬臂梁振动位移的近似解.模型计算的准确性得到了水槽实验的验证.利用数值模拟研究进一步分析了附着物在钝体的振动过程中发挥的作用.本文主要结论如下.(1)在U=0.51m/s和Eell=-2mm的条件下,压电俘能器最大RMS输出功率为1.411mW，相比于E=0的情况增加了6 9.8 8%.(2）随着附着物的特征值从-2 mm增加到6mm，压电俘能器的输出功率均出现先减小再增大后减小的趋势，凸起和凹陷分别在-2 和2 mm取得最大值.对于增强压电俘能器的输出而言，附着物的凸起和凹陷分别存在“最优特征值”1471期第孙维鹏等：钝体表面附着物对低速水流压电俘能器性能影响研究(3)当U=0.51m/s时,Eey=6mm条件下悬臂梁末端振幅为3.0 7 mm，与E=0的情况相比降低了84.83%，附着物的减振效果为大跨度柔性结构的减振提供了参考（4）数值模拟的结果表明：特征值绝对值为2mm的附着物能够令在钝体上形成的压差升高、升力变大，引起钝体愈发剧烈的流致振动；并且附着物的应用会干扰钝体两侧的流体分层，进而对脱落尾涡的大小产生影响.参考文献1赵宏图.国际能源转型进程中的能源安全新挑战.国家安全研究，2022,6(6):88-104,167(Zhao Hongtu.New challenges confrontingenergy security during the ongoing international energy transitions.China Security Studies,2022,6(6):88-104,167(in Chinese)2李瑞婷，胡滨，腰正君等.新谈科普展览设计中叙事知识链的构建-以“同构能源格局 常设展览为例.自然科学博物馆研究，2022,5(5):82-91(Li Ruiting,Hu Bin,Yao Zhengjun,et al.A newdiscussion on the construction of narrative knowledge chain in thedesign of popular science exhibitions:Taking the permanent exhibi-tion of Isomorphic energy pattern as an example.Journal of NaturalScience Museum Research,2022,5(5):82-91(in Chinese)3曹蕾蕾，武建华，樊浩等.考虑可制造性约束的声子晶体多目标拓扑优化.力学学报，2 0 2 2,54(4):113 6-1144(CaoLeilei,WuJian-hua,Fan Hao,et al.Multi-objective topology optimization ofphononic crystals considering manufacturing constraint.ChineseJournal of Theoretical and Applied Mechanics,2022,54(4):1136-1144(in Chinese)4种照辉，姜信洁，何则.国际能源贸易依赖网络特征及替代关系研究：化石能源与可再生能源.地理研究，2 0 2 2，41(12)：3 2 14-3 2 2 8(Zhong Zhaohui,Jiang Xinjie,He Ze.Research on the network de-pendence characteristics