基于仿生尼龙丝的圆柱绕流被动控制.pdf

文章编号：1672-9897(2023)04-0066-10doi：10.11729/syltlx20230019基于仿生尼龙丝的圆柱绕流被动控制陈文礼1,2,3，林隆瀚1,2,3，邓质1,2,3，高东来1,2,3,*1.哈尔滨工业大学结构工程灾变与控制教育部重点实验室，哈尔滨1500902.哈尔滨工业大学土木工程智能防灾减灾工业和信息化部重点实验室，哈尔滨1500903.哈尔滨工业大学智慧基础设施研究中心，哈尔滨150090摘要：通过风洞实验研究了鸟类羽毛仿生尼龙丝对圆柱绕流场的控制效果及机理。仿生尼龙丝布置于圆柱的前驻点，雷诺数为 2.67104，特征变量为尼龙丝长度与圆柱直径之比 L/D。采用表面压力测量系统获取圆柱表面压力系数以分析圆柱周围的空气动力。利用高速粒子图像测速（PIV）系统获取圆柱的二维流场信息，并将获取的信息进行本征正交分解（POD），得到流场瞬时特性和时均特性。实验结果表明：L/D1.0 时，尼龙丝可以显著降低圆柱尾流场的湍动能和雷诺应力，并抑制圆柱的升力系数和阻力系数；当 L/D 足够大时，尼龙丝可以抑制剪切层之间的相互接触，从而改变圆柱的卡门涡街旋涡脱落模式。关键词：圆柱尾流；空气动力；仿生；流动控制；尼龙丝；风洞实验；粒子图像测速；本征正交分解中图分类号：O351.2文献标识码：APassive control on flow past a circular cylinder withbionic nylon wiresCHEN Wenli1,2,3，LIN Longhan1,2,3，DENG Zhi1,2,3，GAO Donglai1,2,3,*1.Key Lab of Structures Dynamic Behavior and Control of the Ministry of Education,Harbin Institute of Technology,Harbin150090,China2.Key Lab of Smart Prevention and Mitigation of Civil Engineering Disasters of the Ministry of Industry andInformation Technology,Harbin Institute of Technology,Harbin150090,China3.Lab of Intelligent Civil Infrastructure,Harbin Institute of Technology,Harbin150090,ChinaAbstract:We investigated the control effectiveness and mechanism of the control of thecircularcylinderflowfieldusingbionicnylonwiresinspiredbybirdfeathersbywindtunneltests.Inthisexperiment,ataReynoldsnumberof2.67104,thebionicnylonfilamentwasarrangedatthefrontstationofthecylinderandthelengthratioL/Dbetweennylonwirelengthandcylinderdiameterwasusedasthecharacteristicparameter.Thesurfacepressuremeasurementsystem was used to obtain the pressure coefficients around the cylinder to analyze theaerodynamic forces acting on the cylinder.The two-dimensional flow field information of thecylinderwasobtainedbyahigh-speedParticleImageVelocimetry(PIV)measurementsystem,and the Proper Orthogonal Decomposition(POD)was used to obtain the instantaneous andtime-averagedcharacteristicsoftheflowfield.TheresultsshowthatatL/D1.0,thenylonwirescansignificantlyreducetheturbulentkineticenergyand Reynolds stress of the cylindrical wake field and suppress the lift and drag coefficientdistributions around the circular cylinder.And at high values,nylon wires can inhibit the收稿日期：2023-02-27；修回日期：2023-03-26；录用日期：2023-07-06基金项目：国家自然科学基金优秀青年科学基金项目（51722805）；国家自然科学基金项目（51978222，51578188）*通信作者E-mail：引用格式：陈文礼,林隆瀚,邓质,等.基于仿生尼龙丝的圆柱绕流被动控制 J.实验流体力学,2023,37(4):66-75.CHENWL,LINLH,DENGZ,etal.PassivecontrolonflowpastacircularcylinderwithbionicnylonwiresJ.JournalofExperimentsinFluidMechanics,2023,37(4):66-75.第37卷第4期实验流体力学Vol.37，No.42023年8月Journal of Experiments in Fluid MechanicsAug.，2023interactionbetweenshearlayersandthuschangethevonKrmnvortexsheddingpatternofthecylinder.Keywords:cylinderwake;aerodynamics;bionic;flowcontrol;nylonwires;windtunneltest;PIV(ParticleImageVelocimetry);POD(ProperOrthogonalDecomposition)0 引言钝体绕流是工程界关注的经典课题，在土木工程、海洋工程、航空航天工程等领域都是热点研究方向。当钝体浸没于流体中运动，且钝体表面旋涡脱落频率与钝体固有频率接近时，会引起钝体的涡激共振。涡激共振会导致桥梁、高层建筑等细长型结构产生剧烈振动，并可能对结构造成重大破坏。目前，工程领域常见的振动控制措施包括结构措施、机械措施和气动措施等1-3。此外，研究者还采用不同方法改变工程结构流动形式，进行流动控制，从而达到减阻、减震、降噪等目的，控制效果较好，具有较强的实用价值。Choi 等4根据有无能量消耗将常见流动控制方法分为主动控制方法和被动控制方法。被动控制方法一般通过改变物体几何形状实现，而主动控制方法一般通过加入外界干扰改变初始流场性质来实现。定常吸吹气是近年来主动控制方法研究热点之一。Gao5和 Yu6等通过研究发现，该方法可有效提高尾流稳定性。尽管主动控制方法可通过调整外部能量条件改变控制参数，但对能量的依赖使其应用受到限制；相比而言，被动控制方法成本低、结构简单、维护方便，更易应用于实际工程。多种被动控制方法已被应用于工程实践中，如在结构表面设置凸起7、凹坑8、导流板9-10、开缝11、开槽12-13和螺纹14等。Chauhan 等15通过研究发现，方柱后侧的分流板使方柱的阻力系数最大下降了 23%。Sahu 等16通过实验研究发现，圆柱后驻点的刚性分流板可以有效控制圆柱结构的涡激振动。这些被动控制方法无需输入额外能量，能够有效控制工程结构的非稳态空气动力。鸟类羽毛可以有效地降低阻力和噪声。受鸟类羽毛启发，仿生被动控制方法在工程领域得到了极大关注。Bocanegra 等17研究了一种可以缓解流动分离且不会明显增大湍动能（Turbulence KineticEnergy,TKE）的仿生涂层。Brcker18、Talboys19、Kunze20、Geyer21、Kamps22等使用自适应柔性涂层在圆柱上进行了一系列实验，流动可视化结果表明：从毛发上产生的行波可以缩短分离流长度，增大旋涡脱落的频率，从而有效控制噪声辐射和旋涡脱落现象。Deng 等23研究发现，小迎角下 S833 机翼后缘的柔性尾翼可以抑制上下层剪切层的相互作用，并降低 TKE 和雷诺应力。Deng 等24对布置不同数量、长度和附着角柔性丝线的圆柱进行了数值模拟，结果表明：当 L/D（丝线长度与圆柱直径之比）为 0.5 的柔性丝线设置于圆柱的前后驻点时，平均阻力和波动升力分别降低了 10.8%和 34.6%。基于数值模拟结果，对布置不同数量、长度和附着角的薄膜的圆柱进行了实验研究25，结果表明：在雷诺数 Re=1.07 104条件下，L/D=2.0、附着角为 060和135180的薄膜控制效果最佳。Deng 等24-25的数值模拟和实验结果表明，在圆柱前驻点布置薄膜能够有效控制圆柱绕流场。为改变圆柱的卡门涡街旋涡脱落模式，优化圆柱气动力，降低湍流度，本文采用高速粒子图像测速技术（Particle Image Velocimetry,PIV）和压力测量技术，对前驻点布置仿生尼龙丝的圆柱绕流场进行研究，基于圆柱的气动力特征、时均流场和瞬态结果，分析仿生尼龙丝的控制效果及控制机理。1 实验装置实验在哈尔滨工业大学土木工程学院风洞与浪槽联合实验室 1 号风洞（SMCWT1）进行。该回流式风洞试验段尺寸为 505 mm 505 mm 1050 mm，以透光性能良好的有机玻璃制成；风速在 425 m/s之间连续可调，湍流度约为 0.36%。本文实验的来流风速 u0控制在 8 m/s。1.1 布置尼龙丝的圆柱模型采用的圆柱模型如图 1 所示。模型展向长度为504.0 mm，直径为 50.0 mm，以聚甲基丙烯酸甲酯3D 打印制成，制造误差为 0.1 mm。图 1 实验模型照片Fig.1 Photo of the test model第 4 期陈文礼等：基于仿生尼龙丝的圆柱绕流被动控制67在不同流场条件下，植物和鸟类羽毛会产生不同的变形26-28，改变周围的流场。尼龙丝也有在不同流场条件下产生不同变形的特性，因此本文采用如图 2（a）所示的尼龙丝仿生植物和鸟类羽毛，以达到控制圆柱绕流场的目的。为将尼龙丝布置于圆柱表面的前驻点上，首先在前驻点处开一狭槽，然后将18 束尼龙丝（每束宽度为 5 mm）两两间隔 5 mm 粘贴于一凹形条上，最后将凹形条嵌入狭槽中，如图 2（b）所示。175 mm1090 mm凹形条尼龙丝（a）尼龙丝照片（b）尼龙丝通过凹形条布置于圆柱表面狭槽图 2 仿生被动控制的实现Fig.2 Implementation of bionic passive flow control实验中，外露于圆柱表面的尼龙丝长度 L 以20 mm 为增量由 10 mm 变化至 90 mm，以无量纲参数 L/D（尼龙丝长度与圆柱直径之比）作为特征变量，则 L/D 以 0.4 为增量从 0.2 变化至 1.8。根据来流风速 u0和圆柱直径 D，本文实验雷诺数为 2.67 104。在控制工况中，圆柱表面布置尼龙丝区域的宽度（沿展向）为 175 mm。本文忽略圆柱两端未布置尼龙丝区域的三维效应。1.2 表面压力测量系统如图 3 所示，在测压环上以 10为间隔均匀布置36 个内径为 1.0 mm 的测压孔。以长度为 650 mm、直径为 1.0 mm 的聚氯乙烯管（PVC）连接测压孔和数字压力测量设备（型号 DSM3400，具有 64 个测量通道的扫描阀）。每个工况的采样频率为 500 Hz，采样时间为 40 s。根据 Irwin 等29的研究结果，本实验中由测压管（长度 650 mm）引起的瞬时压力信号振幅衰减和相位滞后导致的阻尼效应可以忽略。测压环测压环90270 xyxz尼龙丝尼龙丝模型端板0（360）180OO图 3 测压环示意图Fig.3 Schematic diagram of the pressure measurement plane根据 Chen 等30的研究，由表面压力分布估计的圆柱模型表面不稳定空气动力（升力和阻力）系数可以用以下公式积分得到：CD=0.5iCpi cosi（1）CL=0.5iCpi sini（2）Cpi=2(pi p)/u20（3）式中：CD、CL分别为圆柱表面阻力系数和升力系数，Cpi为圆柱表面压力系数；pi为圆柱表面静压，p为来流静压；i为测压孔所在的方位角（如图 3 所示），为相邻测压孔方位角的差值，=10。1.3 PIV 测量系统实验中，采用高速粒子图像测速（PIV）系统获得模型尾流流动特征。该系统由烟雾发生器（Rosco，MiniV）、CCD 高速相机（pco.dimax HS4，最大分辨率 2000 像素 2000 像素）、双脉冲 Nd:YAG 激光器（脉冲能量 425 mJ，波长 532 nm）和数字式延迟发生器组成，如图 4 所示。计算机双脉冲激光器模型激光面CCD高速相机烟雾发生器数字式延迟发生器1 050 mm505 mmu0=8 m/szxy505 mm图 4 PIV 测量实验示意图Fig.4 Schematic diagram of the PIV measurement experiment68实验流体力学http:/在实验中，以烟雾发生器生成的直径为 1.05.0 m 的小油滴作为示踪粒子，并以中跨截面作为目标平面。为确保收敛性，在每个测试工况中均采集 2000 对图像（分辨率为 1790 像素 1000 像素，采样频率为 200 Hz，采样时间为 15 s）。对图像采用互相关算法（有效重叠度为 50%，询问窗口大小为32 像素 32 像素）进行处理，以获得圆柱尾流的瞬时速度场。参照 Park 等31的方法，本文使用 PIV 系统所获流速的不确定度约在 5%以内。2 实验结果及分析 2.1 对气动力的控制效果根据压力测量结果，绘制了无控工况（L/D=0）和控制工况下圆柱中跨截面的平均压力系数和脉动压力系数分布，如图 5 所示。由图 5（a）可知，在无控工况下，圆柱背风侧存在对称分布的负压区平台，这是在圆柱上下表面分离的剪切层相互作用下形成的，导致圆柱近尾流区域形成了再循环区域。在控制工况下，背风侧负压区平台“上升”，表明圆柱中跨截面的阻力系数有降低趋势；平均压力系数在 L/D=1.4 时达到最小值，其分布的对称性则明显有所下降，这是由于尼龙丝在圆柱前驻点上下两侧分布不均匀。图 5（b）显示了不同工况下瞬时压力系数的均方根（脉动压力系数）分布。在无控工况下，脉动压力系数分布与平均压力系数分布类似，同样呈现出中心对称性。圆柱背风侧的脉动压力系数比其他位置更高，表明脉动风荷载在圆柱表面的变化和波动较大。与无控工况相比，控制工况下的脉动压力系数降低（特别是在背风侧），表明模型表面的非定常脉动风荷载得到显著控制，脉动压力系数在 L/D=1.4 时达到最小值。1.5L/D=0.2L/D=0.6L/D=1.8L/D=1.0L/D=0L/D=1.4L/D=0.2L/D=0.6L/D=1.8L/D=1.0L/D=0L/D=1.41.00.500.5平均压力系数脉动压力系数1.01.52.00.50.40.30.20.10060120角度/()180240300360060120角度/()（a）平均压力系数（b）脉动压力系数180240300360图 5 无控和控制工况下的圆柱表面压力分布Fig.5 Pressure distributions around the baseline and controlledcircular cylinder model采用傅里叶变换（FFT）对无控和控制工况下圆柱中跨截面的瞬时升力系数的时程曲线进行分析，得到相应的功率谱，如图 6 所示。在无控工况下，能量集中于功率谱的主导频率 19.300 Hz 上，表明模型2 5002 0001 500幅值1 000500001020频率/Hz3040502 5002 0001 500幅值1 000500001020频率/Hz3040502 5002 0001 500幅值1 000500001020频率/Hz304050L/D=0f =19.627 Hzf =19.321 Hzf =16.460 Hz2 5002 0001 500幅值1 000500001020频率/Hz3040502 5002 0001 500幅值1 000500001020频率/Hz3040502 5002 0001 500幅值1 000500001020频率/Hz304050L/D=1.0L/D=1.4L/D=1.8L/D=0.6L/D=0.2f =16.423 Hzf =19.300 Hz图 6 无控和控制工况下圆柱非稳态升力系数的功率谱分析Fig.6 Frequency spectrum analysis of the nonstationary lift coefficients about the cylindrical model with controlled and baseline cases第 4 期陈文礼等：基于仿生尼龙丝的圆柱绕流被动控制69尾流中有规律的周期性旋涡脱落。该主导频率对应的斯特劳哈尔数（Sr=fD/u0）为 0.121。当 L/D 为0.2 和 0.6 时，功率谱振幅略有下降，主导频率变化很小，表明此时模型周围的气动力几乎不受尼龙丝影响。当 L/D=1.0 时，主导频率下降至 16.423 Hz，功率谱振幅明显下降，表明此时尼龙丝能显著控制不稳定的旋涡脱落模式。为证明仿生被动控制方法对降低圆柱空气动力系数的有效性，比较了无控和控制工况下圆柱中跨截面的平均阻力系数和脉动升力系数。图 7 显示了圆柱中跨截面的平均阻力系数和脉动升力系数的分布。从整体趋势可以看出，随着 L/D 增大，平均阻力系数和脉动升力系数均减小，L/D=1.4 时，两者比 无 控 工 况 下 分 别 减 小 约 40.3%和 84.6%，表 明L/D=1.4 的控制工况对圆柱的升力和阻力均有明显的控制效果。1.21.1平均阻力系数脉动升力系数1.00.90.80.70.60.60.50.40.30.20.1000.20.40.60.81.0L/D（a）平均阻力系数（b）脉动升力系数1.21.41.61.82.000.20.40.60.81.0L/D平均阻力系数脉动升力系数1.21.41.61.82.0图 7 无控和控制工况下的圆柱升力系数和阻力系数Fig.7 Drag and lift coefficients acting on the baseline and controlledcircular cylinder model 2.2 对尾流结构的控制效果本征正交分解（POD）是一种处理离散数据的方法，可以将高阶模型简化为低阶模型，简化复杂的随机过程，从而获得所需的时空信息。其基本思想是将复杂的随机变量分解为与空间特征有关的系数矩阵和与时间特征有关的基函数。其中，基函数的确定原则32是：在每一次分解过程中，使最低阶模态含有最多的能量。改进的 snapshot POD 算法由 Sirovich32引入PIV 测量，其有效性被 Deri33和 Prothin34等证实。下文介绍流体力学领域中 snapshot POD 方法的相关理论。POD 模态参数计算采用如下公式：CAi=iAi（4）Ni=1i=u u（5）C=1nuTu（6）u（x,y,t）=u（x,y）+u（x,y,t）（7）u式中：A 为含有模态系数的 POD 模态矩阵；为特征值；C 为自相关矩阵；N 为时间序列长度；为脉动流速向量。根据 Sirovich32的理论，脉动流速可分解为一系列正交基 i与模态系数 ai乘积的和：u(x,y,t)Ni=1ai(t)i(x,y)（8）式（8）为 POD 模态分解的核心。通过模态分解，将二维随机变量脉动速度矩阵分解为 2 个一维变量，实现降维，简化计算，得到流场主要模态阶数的运动特征分布，进而得到代表主要流动能量不同阶数的流动模式信息。通过分析 PIV 系统获得的 2000 对流场照片，采用 snapshot POD 算法，提取了圆柱流动特征和时间演变信息。采用 POD 模态分析得到的前 15 阶模态对流场进行重构，得到一系列瞬时涡量图（图 8），选择其中 1 个周期（T）的 4 幅瞬时涡量图进行分析。在无控工况下，旋涡在圆柱上下两侧被交替拉伸，导致尾流场中出现周期性的交替脱落，如图 8（a）所示。当 L/D 从 0.2 增大至 1.0 时，旋涡脱落模式几乎没有太大变化，仅旋涡形成长度随 L/D 增大而略有下降，如图 8（b）（d）所示。当 L/D=1.4 时，尼龙丝诱导产生的旋涡结构改变了之前的旋涡脱落模式，尾流场中上下两侧剪切层的接触点被推迟，剪切流被拉长，旋涡脱落结构变得稳定，如图 8（e）所示。当L/D=1.8 时，尼龙丝诱导产生的旋涡结构和剪切层之间的相互作用压缩了剪切层，抑制了旋涡脱落，如图 8（f）所示。上述结果表明，尼龙丝对增大尾流场的稳定性具有显著效果。根据 Chen 等30的研究，湍动能（TKE）是一个与动态风荷载有关的流动特征参数，可用于评估模70实验流体力学http:/型表面压力脉动特性和气动力幅值。图 9 显示了无控和控制工况下的湍动能和流线时均分布（图中红色数字为 TKE 峰值）。可以发现，在无控工况下，圆柱上下两侧流动分离，在圆柱背风侧尾流处形成了1.00.500.5y/D1.001.02.0 x/D3.04.001.02.0 x/D3.04.001.02.0 x/D3.04.001.02.0 x/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dx/Dt0+0.25T t0+0.50Tt0+0.75Tt0+T3.04.01.00.500.5y/D1.001.02.03.04.001.02.03.04.001.02.03.04.001.02.03.04.01.00.500.51.001.02.03.04.001.02.03.04.001.02.03.04.001.02.03.04.01.00.500.5y/Dy/D1.001.02.03.04.001.02.03.04.001.02.03.04.001.02.03.04.01.00.500.5y/D1.001.02.03.04.001.02.03.04.001.02.03.04.001.02.03.04.01.00.500.5y/D1.001.02.03.04.001.02.03.04.001.02.03.04.001.02.03.04.01.0Vorticity:0.500.51.0（a）L/D=0（b）L/D=0.2（c）L/D=0.6（d）L/D=1.0（e）L/D=1.4（f）L/D=1.8图 8 PIV 系统获得的瞬时涡量结构演变图Fig.8 Instantaneous flow structures evolution measured by PIV measurement system第 4 期陈文礼等：基于仿生尼龙丝的圆柱绕流被动控制71再循环区域，其中存在一对关于中心线对称的旋涡，旋涡形成长度约为 x/D=1.6。这对旋涡是导致圆柱前后表面压力分布不均匀的原因。随着与圆柱背风侧表面距离增大，湍动能的值先增大后减小。在控制工况下，随着 L/D 增大，再循环区的面积几乎没有改变，旋涡形成长度都约为 x/D=1.5；而湍动能峰值先减小后增大，L/D=1.4 时，湍动能峰值达到最小值，比无控工况下减小了约 50%。上述结果表明，本文的仿生被动控制方法可以有效减小尾流流场的湍流度。此外，当 L/D 为 1.4 和 1.8 时，尾流旋涡的不对称性是由圆柱上下两侧尼龙丝变形不均匀所导致的，如图 10 所示。TKE0.050.100.150.200.251.00.500.5y/D1.001.5x/D（a）L/D=00.232 10.197 70.154 10.128 20.173 13.04.51.00.500.5y/D1.0（b）L/D=0.21.00.500.5y/D1.0（c）L/D=0.61.00.500.5y/D1.0（d）L/D=1.01.00.500.5y/D1.0（e）L/D=1.41.00.500.5y/D1.0（f）L/D=1.801.5x/D3.04.501.5x/D3.04.501.5x/D3.04.501.5x/D3.04.501.5x/D3.04.50.231 80.232 10.197 70.154 10.128 20.173 10.231 8图 9 无控和控制工况下的湍动能分布和时均尾流场Fig.9 TKE distributions and time-averaged wake field about the baseline and controlled circular cylinder model（a）L/D=0（b）L/D=0.2（c）L/D=0.6（d）L/D=1.0（e）L/D=1.4（f）L/D=1.8图 10 无控和控制工况下的尼龙丝变形Fig.10 The nylon wires deformation of the baseline and controlled circular cylinder model72实验流体力学http:/图 11 显示了无控工况和控制工况下的雷诺剪切应力（Reynolds Shear Stress,RSS）时均分布。雷诺剪切应力被认为是衡量流场中速度变化幅度的一个流场特征，可用于表示不同流速流层间的动量交换程度。在无控工况下，雷诺剪切应力呈明显的 4 峰值中心对称结构，雷诺剪切应力峰值较大，且应力变化梯度也较大。在控制工况下，随着 L/D 增大，4 个峰值先减小后增大，L/D=1.4 时，回流区的 2 个应力峰几乎消失，这是由于尼龙丝诱导产生的旋涡结构抑制了圆柱近壁面的速度波动；远离圆柱壁面的 2 个峰值均比无控工况下减小了约 30%，表明本文的仿生被动控制方法对抑制圆柱尾流区的动量交换效应具有明显效果。图 12 展示了无控工况和控制工况下顺流向雷诺正应力（Reynolds Normal Stress,RNS）时均分布。在无控工况下，顺流向雷诺正应力呈明显的“峰值对1.00.500.5y/Dy/Dy/Dy/Dy/Dy/D1.01.301.56.1756.553（a）L/D=03.04.51.00.500.51.001.56.9166.373（b）L/D=0.23.04.51.00.500.51.001.54.6415.823（c）L/D=0.63.04.51.00.500.51.001.55.2385.425x/Dx/Dx/Dx/Dx/Dx/D（d）L/D=1.03.04.51.00.500.51.001.53.9174.864（e）L/D=1.43.04.51.00.500.51.001.55.4065.573（f）L/D=1.83.04.56420246Reynolds Shear Stress/(m2s2):图 11 无控和控制工况下的圆柱尾流场的雷诺剪切应力分布Fig.11 RSS concentrations in the wake field about the baseline and controlled circular cylinder model1.00.500.51.001.53.2213.2773.04.51.00.500.51.03.2093.3321.00.500.51.02.7873.0271.00.500.51.02.9312.9021.00.500.51.02.5082.7271.00.500.51.02.9283.1450.5 1.0 1.5 2.0 2.5 3.0 3.5 4.001.53.04.501.53.04.501.53.04.501.53.04.501.53.04.5Reynolds Normal Stress/(m2s2):（a）L/D=0（b）L/D=0.2（c）L/D=0.6x/Dx/Dx/Dx/Dx/Dx/D（d）L/D=1.0（e）L/D=1.4（f）L/D=1.8y/Dy/Dy/Dy/Dy/Dy/D图 12 无控和控制工况下圆柱的顺流向雷诺正应力分布Fig.12 Pulsating velocity distributions in streamwise direction of the baseline and controlled circular cylinder model第 4 期陈文礼等：基于仿生尼龙丝的圆柱绕流被动控制73峰值”中心对称结构。在控制工况下，随着 L/D 增大，2 个峰值先减小后增大，L/D=1.4 时，峰值比无控工况下减小了约 16%，雷诺正应力“轮廓”则呈顺流向“先压缩后拉伸”变化，L/D 为 1.4 和 1.8 时的雷诺正应力不对称分布也是由圆柱上下两侧尼龙丝变形不均匀所导致的。3 结论本文采用压力测量系统和 PIV 系统对布置尼龙丝的圆柱气动力特征和绕流场结构进行了风洞实验测量，实验结果表明 L/D 对控制效果具有显著影响。主要结论如下：1）L/D 对表面压力分布具有明显影响。在高L/D 值时，脉动升力和平均阻力均被明显抑制，且圆柱的主导频率与 L/D 密切相关。2）L/D 对被动控制机理存在影响。低 L/D 值下，尼龙丝对圆柱旋涡脱落模式几乎没有影响；高L/D 值下，尼龙丝诱导产生的旋涡与剪切层之间的相互作用阻碍了圆柱上下两侧剪切层之间的相互作用，从而抑制了传统的卡门涡街旋涡脱落模式。3）在高 L/D 值下，尼龙丝抑制了雷诺应力和湍动能分布。此外，尾流中顺流向脉动速度峰值下降，在高 L/D 值下其分布区域被压缩。成对的旋涡结构在尼龙丝附近产生，然后向圆柱尾流区移动，在高L/D 值下，尼龙丝诱导产生的旋涡结构延迟了圆柱上下两侧剪切层之间的接触。来自圆柱的分离流和尼龙丝诱导产生的旋涡结构对旋涡脱落产生了显著的抑制。这些现象表明，尼龙丝有效增大了圆柱尾流的稳定性，并抑制了作用于圆柱周围的气动力。参考文献：陈政清.桥梁风工程M.北京:人民交通出版社,2005.CHEN Z Q.Bridge wind engineeringM.Beijing:ChinaCommunicationsPress,2005.1刘庆宽,孙一飞,张磊杰,等.凹痕对斜拉桥斜拉索气动性能影响研究J.工程力学,2019,36(S1):272277.LIU Q K,SUN Y F,ZHANG L J,et al.Study on theinfluenceofdentonaerodynamicperformanceofstaycablesof cable-stayed bridgeJ.Engineering Mechanics,2019,36(S1):272277.doi：10.6052/j.issn.1000-4750.2018.05.S0532董欣,丁洁民,邹云峰,等.倒角化处理对于矩形高层建筑风荷载特性的影响机理研究J.工程力学,2021,38(6):151162,208.DONGX,DINGJM,ZOUYF,etal.Effectofroundedcorners on wind load characteristics of rectangular tallbuildingsJ.Engineering Mechanics,2021,38(6):151 162,208.3doi：10.6052/j.issn.1000-4750.2020.07.0451CHOIH,JEONWP,KIMJ.ControlofflowoverabluffbodyJ.Annual Review of Fluid Mechanics,2008,40:113139.doi：10.1146/annurev.fluid.39.050905.1101494GAO D L,MENG H,HUANG Y W,et al.Active flowcontrolofthedynamicwakebehindasquarecylinderusingcombined jets at the front and rear stagnation pointsJ.PhysicsofFluids,2021,33(4):047101.doi：10.1063/5.00431915YUHY,CHENWL,XUZH,etal.Wakestabilizationbehind a cylinder by secondary flow over the leewardsurfaceJ.PhysicsofFluids,2022,34(5):055110.doi：10.1063/5.00907976BEARMAN P,BRANKOVI M.Experimental studies ofpassive control of vortex-induced vibrationJ.EuropeanJournalofMechanics-B,2004,23(1):915.doi：10.1016/j.euromechflu.2003.06.0027ZHOUB,WANGXK,GUOW,etal.Controlofflowpasta dimpled circular cylinderJ.Experimental Thermal andFluidScience,2015,69:1926.doi：10.1016/j.expthermflusci.2015.07.0208ZHOUX,WANGJJ,HUY.Experimentalinvestigationonthe flow around a circular cylinder with upstream splitterplateJ.JournalofVisualization,2019,22(4):683695.doi：10.1007/s12650-019-00560-x9GAODL,HUANGYW,CHENWL,etal.ControlofcircularcylinderflowviabilateralsplitterplatesJ.PhysicsofFluids,2019,31(5):057105.doi：10.1063/1.509730910GAODL,CHENWL,LIH,etal.Flowaroundacircularcylinder with slitJ.Experimental Thermal and FluidScience,2017,82:287301.doi：10.1016/j.expthermflusci.2016.11.02511KONWM,LEUNGYC,CHENJJJ.FlowpastV-groovecircularcylindersJ.AIAAJournal,1987,25(6):806811.doi：10.2514/3.970412LIMHC,LEESJ.Flowcontrolofcircularcylinderswithlongitudinal grooved surfacesJ.AIAA Journal,2002,40(10):20