基于虚功率原理的鲹科鱼类波状推进分析.pdf

andAppliedMechanics,2023.55(7):1439-1450Huang Kaijun,Yu Yongliang.Analysis of carangifom undulation based on virtual power principle.Chinese Journal of Theoretical引用格式：黄恺俊，王虚功率原理的修科鱼大波状推进分机2023,55(7):1439-1450流体力学Chinese Journal of Theoretical and Applied MechanicsJul.,20232023年7 月力Vol.55,No.7期第55卷第报学学基于虚功率原理的科鱼类波状推进分析黄恺俊余永亮2)(中国科学院大学生物运动力学实验室，北京1 0 0 0 49)摘要自然界中的大多数鱼类通过波状摆动的方式实现推进，这是动态变形的鱼体和周围流体相互作用的结果，研究推进中流体对鱼体变形的响应不仅可以增强对波状推进的认识，还可以为流动控制提供依据.以科模式推进的仿生二维模型为研究对象，通过数值计算获得鱼体波动产生的流场以及鱼体受到的流体力数据.基于虚功率原理，将鱼体受力分解为4部分，分别是鱼体边界变速运动的瞬时贡献、流场中流体旋转和应变速率相对大小的贡献、壁面剪切应力的类摩阻分量和壁面摩阻分量.结果表明，当鱼体波状摆动产生推力时，鱼体边界变速运动是主要的正推力来源，并且该项8 0%的推力贡献来源于2 0%的鱼尾部分的边界变速运动.鱼尾两侧边界层中的流体旋转和应变速率的相对大小和壁面摩阻对推力都是负贡献.对于低雷诺数的情况，流体旋转和应变速率的相对大小的负贡献低于壁面摩阻的负贡献，而在高雷诺数的情况下，流体旋转和应变速率的相对大小的负贡献强于壁面摩阻的负贡献.壁面类摩阻分量相对于其他3项总是较小的.结合标度律分析，在摆动推进的标度关系中，与雷诺数无关的推力部分是由边界的变速运动、流场中流体旋转和应变速率共同提供，且流体旋转和应变速率也贡献了摆动推力中与雷诺数有关的部分，而这一部分接近类摩阻和摩阻的一半，同时，类摩阻和摩阻还提供了常阻力分量.关键词同虚功率原理，科模式，波状推进，推力分解中图分类号：Q352,R318.01文献标识码：Adoi:10.6052/0459-1879-23-076ANALYSISOFCARANGIFORMUNDULATIONBASEDONVIRTUALPOWERPRINCIPLE1)Huang KaijunYu Yongliang2）(The Lab of Biomechanics for Animal Locomotion,University of Chinese Academy of Sciences,Beijing 100049,China)AbstractMost fish in nature achieve propulsion through undulatory movements,which are the result of the interactionbetween the deforming fish body and the surrounding fluid.To study the response of the fluid can enhance ourunderstanding of undulatory propulsion and flow control.A two-dimensional deforming airfoil is used to model thecarangiform fish.The flow field generated by fish body and the fluid forces acting on the fish body were obtained byusing computational fluid dynamics.Using the principle of virtual power,the thrust on the fish body was decomposedinto four parts,which are the instantaneous contribution of the boundary acceleration,the contribution of the relativemagnitude of fluid rotation and strain rate in the flow field,the wall friction-like component and the wall frictioncomponent.The results show that the instantaneous contribution of the boundary acceleration is the main source of2023-03-08收稿，2 0 2 3-0 5-1 2 录用，2 0 2 3-0 5-1 3网络版发表.1)国家自然科学基金（1 2 1 7 2 355,1 1 6 7 2 2 91）和中央高校基本科研业务费专项资金（E1E42201)资助项目.2)余永亮，教授，主要研究方向为生物运动力学、空气动力学.E-mail:力14402023年第55卷报学学positive thrust.The 80%of the thrust contribution of this term comes from the instantaneous boundary accelerationmovement of the rear 20%of the fish body.The fluid rotation and strain rate in the boundary layer on both sides of thefish tail and the friction contribute to resistance.For high Reynolds number,the negative contribution of the relativemagnitude of fluid rotation and strain rate is stronger than that of wall friction,while for low Reynolds number,it is lowerthan that of wall friction.However,the wall friction-like component is always smaller compared to the other three terms.In the analysis of the scaling law of undulatory propulsion,it was found that there is a component independent of theReynolds number which is provided by the first two parts,while the other component that is dependent on the Reynoldsnumber is provided by the last three parts.Furthermore,the frictional force and the friction-like force provide constantresistance.Keywordsvirtual power principle,carangiform,undulatory propulsion,thrust decomposition引言自然界中的大多数鱼类进化出波状摆动的方式实现推进1-3,其中振幅包络线从前部到后部逐渐增加的波状摆动更有利于产生推力4,这类波动模式通常被称为科模式.科模式与鳗模式、金枪鱼模式都是通过身体和尾鳍的摆动进行推进的，科模式的鱼体前1/3不波动，后体摆动的波幅随体长增大3，常见的斑马鱼、草鱼、鲤鱼等都采用这类推进模式.鳗模式的鱼类全部鱼体都大幅摆动，而金枪鱼模式的鱼体仅尾部往复摆动.科鱼类的游速较高，常被用于水下仿生航行器的仿生设计中.因此，深入研究科模式推进的流体力学机制不仅可以增强对波状推进的认识，还可以为水下航行器的设计提供参考.Videler2和Lighthill5预测了鱼尾摆动后的流场尾迹，描述波状推进的鱼类在尾部会脱落出交错排列的离散尾涡，被称为反卡门涡街.前人通过数字粒子图像测速仪(DPIV)观测到太阳鱼的尾迹中存在交替的涡结构6-8 ,在振荡翼的模型实验中也发现相同的流动9-1 1 ;通过对二维波状翼型的数值模拟,Dong等1 2 、Deng等1 3 和Liu等1 4 也同样观察到这种典型的尾迹结构.与阻力型的卡门涡街不同，反卡门涡街是一种推力型的尾流结构.Yu等1 5 基于波状推进的鱼体产生反卡门涡街的共性，利用涡动力学分析方法，推导并提出了统一描述波状推进的水动力与运动参数、雷诺数的关联关系的标度律.该规律对鳗模式、科模式以及金枪鱼模式都适用.在研究鱼体运动与受力关系方面，近1 0 年来也有不少实验和计算分析的工作.Gazzola等1 6 结合实验和数值模拟的结果，通过流体力中推、阻力的匹配关系,得到巡游速度的标度律.Thekkethil 等1 7 使用统一的运动学模型进行二维数值模拟，该模型涵盖各种鱼体波动和鱼尾刚体俯仰运动，并研究了波动和俯仰运动的等价性以及推力产生的机制.Hu等1 8 对鱼体的低频波状摆动进行数值模拟,获得鱼体阻力与波动参数之间的标度律.Gupta 等1 9则不光给出波动鱼体受力，还研究了输出功率、效率的标度律关系.Gao 等2 0 在Yu等1 5 的工作基础上，基于压差力标度律和摩阻力标度律的匹配关系，也给出巡游速度的标度律.这些工作很好地给出运动与水动力性能的宏观联系，为进一步分析推进的流体力学机制提供数据支持.事实上，通过计算手段获得流场数据只是全面了解流动的第一步，为了找到流动和受力的定量关系，引入气动力理论的方法可以从流动数据中提取物理规律2 1 .目前有两种方法可以建立起物体所受力与流场结构之间的联系.一种是Burgers22,Lighthil23 和Wu24 提出的导数矩的方法,通过流体的涡量一阶矩计算流场动量，并通过动量变化求得物体受力.后来被Wu等2 5 发展到有限域的涡量矩理论.另一种是Quartapelle等2 6 、Howe27和Chang28采用的投影法,通过引入一个势函数，将力分解为附加质量力，兰姆矢量相关的旋涡贡献，以及黏性相关的项.Yu29引入虚速度的概念，通过变分得到了流体力学中的虚功率原理，并指出引入的势函数实际上是物体虚运动引起的虚速度场，虚速度场具有远场衰减的特性，在受力分解中包含物面积分项和流体结构的体积分项，并分析各项的伽利略不变形式。为了研究科模式波状推进的流体动力1441期第黄恺俊等：基于虚功率原理的科鱼类波状推进分析响应，从虚功率原理出发，可以深入分析变形体的边界瞬时流体力效应和流场结构对流体力的贡献本文选取科模式推进的鱼体作为研究对象，研究鱼游波状推进的流体力机理.第1 节介绍鱼体波动模型以及基于虚功率原理的分析方法.第2 节通过研究典型参数下科鱼波状摆动流体力的分解来阐明推力的主要来源，然后讨论不同雷诺数条件下推力与Strouhal数的关系，并结合波状推进的标度律分析，研究标度律公式中各项物理意义.1计算模型与分析方法1.1鱼游模型借鉴前人1 2,1 7,1 9,30 研究科模式鱼类游泳的物理模型，本文采用NACA0012翼型来模拟鱼体，翼型的中线代表了鱼体在静态平衡时的脊柱.鱼体脊柱波动的侧向位移y(x,t)由Videler等31 通过实验提供的行波方程来模拟2元y(x,t)=Am(x)cosOxL(1)其中，入和c分别表示波动波长和波速，波速是频率与波长入的乘积.如图1 所示，鱼体静止时，头部位于x=0,尾部在x=L处.Am(x)表示不同身体位置的最大振幅.对于科模式的波动，通常认为波长入等于一倍体长，身体的振幅包络线Am(x)可由二次函数给出2Am(x)=ao+aix+a2x(2)其中 Am(0)=0.02L,Am(0.2)=0.01L,Am(1)=0.1L.图 2给出包络线示意图（虚线）和鱼体中线的侧向位移波动曲线(实线)，可以看出，最大振幅出现在尾尖处，而最小振幅大概在0.2 3L处.鱼体头部有小幅的波动，而尾部相对于头部有更高的振幅.本文定义鱼体尾部最大振幅为鱼体波动的振幅，即波动振幅A=0.1L.文图1 仿鱼翼型波状摆动示意图Fig.1Schematic of a fish-like undulatory airfoil0.100t=OTt=T/80.075T/40.050t=31/80.0250-0.025-0.050-0.075-0.10000.20.40.60.81.0 x/L图2 一个波动周期不同时刻的鱼体中心线及其包络线Fig.2 The centerline of the fish body at various time instants for oneundulation cycle and their envelops在本文中，由波动鱼模型引起的二维非定常流动,可由以下无量纲化的不可压缩Navier-Stokes(N-S)方程控制Vu=0Qu+(u.)u=一-Vp+vu(3)tPP其中，u是速度场，p是压强场，为流体介质的黏性系数.本文选取特征的参考速度U(来流速度),鱼体的体长L以及流体密度对所有的物理量进行无量纲化，则体现流场惯性力和黏性力比值大小的雷诺数为Re=pUL/.另一个常用的无量纲数是Strouhal数，它表征了流动的非定常性，定义为St=2fA/U.本文采用的数值求解N-S方程的程序是开源程序包OpenFOAM的pimpleFoam求解器，程序的验证和确认见附录，该程序在文献1 5 中已得到了验证.1.2基于虚功率原理的推力分析方法通过数值手段获得鱼体波状摆动的流场信息，便可通过鱼体表面的法向应力pn和切向应力T=un积分，得到鱼体的受力FF=S)(pn+T)ds(4)其中，法向量n指流体的外法向量，和鱼体表面的外法向量方向相反.该公式简洁地给出鱼体与流体的相互作用结果，可用于验证以下基于虚功率原理的推力分析方法.为了建立流动与受力的定量联系，本文引入Yu29提出的虚功率原理的分析方法.假定鱼体在前进方向上产生瞬时的虚运动u，则流场中的虚功率关系可得力14422023年第55卷报学学(uo);QuduiLdV+i-udVxjdt(5)其中，u*是由鱼体虚运动引起的虚速度场，ji代表施加在流体微团上的表面力的张量(称为应力张量).式（5）的物理意义是物体所受水动力沿着虚运动方向所做的虚功率等于流场中的内应力和广义力的虚功率.广义力包含流场中的体积力fi和惯性力dui/dt.进一步地，方程(5)左侧的值就是推力FT，即水动力沿着虚运动方向的分量，所以对方程(5)右侧各项的分析将有助于理解推力的物理成因.Yu29在文中提到，可以假定虚流场是不可压缩势流.为了使势函数唯一，采用了在无穷远处势函数为0 的边界条件2 6.2 8-2 9.因此，虚流场满足的方程和边界条件如下V2=0,on VVon=uon,on Sfish(6)IVdl0,0,on S)其中V是虚流场的速度矢量u*，S f i s h 和S分别表示鱼体表面和无穷远边界.虚流场是无源无旋的，可通过在鱼体表面分布基本解(源汇和偶极子)以满足方程（6)的鱼体和无穷远处的边界条件.将虚流场方程(6)代入虚功率关系(5)，并考虑流动不可压缩且体积力可忽略，则表达式(5）右式的两项分别为ou(nxw).VodS+xjSfish(nxw)uods(7)SfishduduVodVnddS+dtSfishdtduodV.(u.Vu)dV:Pndds-fishdt2pdQdV(8)从式（7)可以看出流场的流体内应力对虚应变做功的功率转换为壁面的切应力T所做的虚功率，一部分是对虚速度场做功的功率，另一部分是对鱼体虚运动做功的功率.前者被称为类摩阻（该说法来自Chang28),后者正好为摩阻在推进方向的大小.从式(8)可以看出流场中惯性力的虚功率可以转换为与鱼体边界加速度有关的表面积分和与流场中对流散度有关的体积分.表面积分表明壁面的瞬时运动即可确定部分虚功率的大小，其值可以看作在黏性流动中鱼体运动状态瞬时所确定的推力贡献.而流体对流的散度正好等于负的2 倍Q（Q 是速度梯度张量的第二不变量，定义为Q=0.5(2I2-DI)，其中2 和D分别是速度梯度分解成的反对称张量和对称张量)，Q0代表流体质点的旋转强于应变速率，Hunt等32 就提出用Q0作为旋涡的识别.因此，式(8)中右边第二项可以看作是流场中流体的旋转和应变速率相对大小对推力的贡献，为了方便，我们称之为Q贡献.因此,本文将鱼体受到的推力FT分解为4个分量du1ndds2pdQdV-Odtu*1f,u(nxw)uds(9)分别表示边界加速度对推力的贡献FT,a、流场中Q对推力的贡献FT,Q、边界的类摩阻分量FT,以及摩阻分量FT.f，其推力系数及各分量的系数定义为FTFT.aFT.0CT.CT.a=T,CT.Q=T二30U2L2(10)FT.vFT,f,CT.f=TCT1ou定义推力系数的周期均值算符一六Cto+T()dt,其中T代表鱼的波动周期.本文基于虚功率原理出发得到的式(9)与Gao等33 从压强Poission方程出发得到的结果是一致的.因为满足方程(6)的虚流场是势流场，所以u增加多少倍，和对应增加多少倍.推力分解式(9)中各项积分后除以u相当于一种归一化的计算，u的取值不会影响每一项最终的计算结果.本文为了方便式(9)的分析，在计算中取u的模为1.2结果与讨论2.1典型问题的推力分解为了考察式(9)受力分解与式(4)在推力分析上的一致性，选择无量纲摆尾幅度为A=0.10，波长为入=1.0，波速为c=2.0，流动的雷诺数为Re=5000作为典型的案例，对比了一个摆动周期里推力系数随时间的变化,如图3(a)所示.结果表明基于虚功率原理的推力曲线与由式（4)直接积分的推力曲线重1443黄恺俊等：基于虚功率原理的科鱼类波状推进分析期第0.5Eq.(4)0.4Eq.(9)0.30.20.10-0.1-0.2-0.300.20.40.60.81.0t/T(a)式(4)以及式(9)分别计算的推力系数CT(a)The time-dependent thrust coefficient C,calculatedby Eq.(4)and Eq.(9)0.50.40.30.20.10-0.1-0.2-0.300.20.40.60.81.0t/T(b)式(9)计算的 Cr.,C.o Cr.w Cr.f.(b)The time-dependent Cr.a,Cr.o,Cr.Cr.calculated by Eq.(9)图3推力系数在一个周期内的变化Fig.3The time-dependent thrust coefficient in an undulation cycle合，具有一致性.进一步地，图3(b)给出了4个分量CT,a,CT,CT,和CT,随着时间变化的曲线,发现边界加速度对推力始终是一个正的贡献，而流场中Q对推力有很大的负贡献，两者相位几乎相反.边界的类摩阻分量相对于其他分量要小得多，几乎可以忽略不计，而摩阻分量在一个周期内始终是阻力.四者的周期均值分别为0.2 47 4，-0.1 1 92，0.0 0 0 3和-0.0 6 0 7，可见CT,v的确可以忽略，推力主要由CT,a,CT,o,CT.f3项组成，其中CT,o接近摩阻CT,f的两倍.接下来我们将对推力贡献最重要的两部分CT.a和CT.o做进一步的分析.对于二维波动模型，可以对边界加速度对推力的贡献FT.a做进一步的推导.分别沿着上、下表面进行积分可以得到dudunddsnddsdtdtdu+dtdu(11)dt其中，下标+和下标_分别代表上下表面的物理量，y+和y-分别代表上下表面的曲线方程，y=oy/x代表曲线对x的导数.根据变形规则，在相同的x坐标，鱼体边界与鱼体脊柱波动一致，则边界加速度的值du/dt是相同的.边界加速度可由波动方程对时间的二阶导数获得，流体在鱼体表面的外法向量可由曲线方程获得，具体关系如下du+du0dtdt0t2(12)n+n_=(-)将方程(1 2)代入方程(1 1),我们可以得到2FT,a=?+一d_)dx(13)00t2基于式（1 3)我们可以看出由边界加速度引起的推力贡献沿着鱼体中线的分布与两部分相关因素相关.一部分是侧向运动的方向加速度（波动加速度)，另一部分是上下表面的虚速度势的差值.根据势流的特点，鱼体表面的虚速度势是由当前时刻的鱼体外形决定的，因此2y2元=-4元 Am(x)F2 cosot2(14)Ap(x,t)=+-图4(a)给出了一个摆尾周期内边界加速度对推力贡献沿着中线的分布，其中实线和点被用来区分摆尾前半周期和后半周期的特定时刻.从图4中可以看出，前半周期和后半周期的摆尾动作产生的FT.a是一致的.从曲线看，头部区域由于变形小，加速度小，对FTa的影响就小，尾部区域的值波动大，大部分时刻都处于正值，这意味着对推力有大贡献.进一步地，为了更好地展示鱼体各个部位变速运动带来的贡献，将鱼体沿中线等分成5个区域，并统计了5个区域分别对推力贡献的周期均值，如图4(b)所示.从数据上可以看出，前3个区域的均值与整体均值之比大约为3%4%，且均为负贡献，处于尾部的第4和第5区对整体正推力有很大的正贡献，特别是加速度最大的第5区，虽然只占鱼体长度的2 0%，但这个区域的贡献占整个贡献的8 3.2 5%.根据式(9)右式第二项可知，流场中的Q分布对推力的贡献依赖两个因素，一个是Q的强度，另一个力14442023年第55卷报学学2.52.01.51.00.50-0.500.20.40.60.81.0 x/L(a)边界加速度对推力的贡献沿着鱼体中线的分布图(a)Distribution diagram of the contribution of the boundaryacceleration to thrust along the center line of fishzonelzone2zone3zone4zone50.2583.25%0.200.150.1027.30%0.053.14%4.02%-3.39%0zonelzone2zone3zone4zone5(b)鱼体分区示意图以及不同分区对推力贡献的时均值(b)Schematic of fish body zoning and the time-averagedcontribution of different zone to thrust图4边界加速度对推力的贡献Fig.4 The contribution of the boundary acceleration是虚速度势.为了更好地展现流场信息，图5给出了4个典型时刻的流场涡量分布、Q值分布、Q对推力贡献分布，以及虚速度势的分布.从涡量场来看，鱼体上半部分边界层内为负涡量，下半部分为正涡量，尾迹呈现典型的反卡门涡街结构.从Q值图上可以看出，在鱼体上下部分的边界层处，Q值仅在尾部附近有明显的分布值，且两侧的值相反，意味着尾部拍动的过程中，“迎风面”（图中显示的鱼体尾部下表面）边界层内的流体应变速率强于旋转（Q0).在鱼体头部,存在很强的应变速率区（Q0).对于脱落到流场中的自由旋涡，无论是正涡还是负涡，中心都由流体旋转主导（Q0),且在旋涡周围存在较弱的应变速率区（Q0),故脱落出去的旋涡中心贡献了负推力，旋涡周围呈现一个较弱的正贡献.为了更好地定量展示流动结构对推力的贡献，我们基于流场中Q分布的特点，人为对整个流场进行了分区.如图6 所示，zoneo代表了距离鱼头部一定距离的流场，该区域的Q值都很小.zonel是鱼体头部附近的区域，其Q分布有正有负.zone2zone5分别代表身体波动的部分，且波动的振幅逐渐增强，这些区域分布的Q强度也逐渐增强.区域zone_wake则是尾迹区，该区域又被分为对应的5个子区，由近及远分别为wakelw a k e 5.图7(a)给出了各个分区Q分布对推力的贡献随Q-2p0Q-15一9-33815-200-120-4040120200-20-12-4412200.10-0.060.020.02 0.060.10图55t=0T,t=T/8,t=T/4,t=3T/8时刻的涡量场(第1 列)、Q分布图(第2 列)、Q对推力的贡献分布图(第3列)和虚速度势场(第4列)Fig.5 Vorticity field(first column),Q distribution diagram(second column),Q distribution diagram of contribution to thrust(third column)andvirtual velocity potential field (fourth column)at t=OT,t=T/8,t=T/4,t=3T/81445黄恺俊等：基于虚功率原理的科鱼类波状推进分析期第zoneozonelzone2zone3zone4zone5zone_wake0.50-0.5wakelwake2wake3wake4wake5-1.0-0.500.51.01.52.0X图6 流场分区示意图Fig.6Schematic of flow field zoning0.10-0.1CT.0.2oncoCT.0.2oclCT.0.20ne2-0.2CT.0.2one3CT.0.zonc4CT.0.zones-0.3Cr.0.zone_wakeCT.000.20.40.60.81.0t/T(a)各个分区分布对推力的贡献随着时间的变化(a)The time-dependent contribution ofQ distribution to thrust in each zone0.02-14.20%107.39%0.38%12.22%01.78%-16.88%-0.039.31%-0.05C-0.08-0.10-0.13ouoztouozSouozzone_(b)各个分区Q分布对推力贡献的周期均值，其中的百分占比由分区贡献比上总贡献求得(b)The time-averaged contribution of Q distribution to thrust ineach zone.The percentage is obtained from the contributionof each zone to the total contribution.图7 Q分布对推力的贡献Fig.7The contribution of Q distribution to thrust着时间变化的曲线，图7(b)给出了各个分区对推力贡献的周期均值，它们具有如下几个特点.（1）从时变曲线来看,zoneo 和zonel 的推力贡献值都在零附近小幅振荡.由于zoneo区的流场扰动小，Q值小，整体流动几乎没有推力贡献.zonel的头部区存在强正强负的贡献，但整个摆动周期中变化小，其贡献也较小.(2)鱼体波动的zone2z o n e 5这4个区因摆动振幅逐渐增大而对推力贡献值也波动增强.其中，zone2,zone3和zone4的波动曲线仍在零值附近振荡，从均值来看，zone2和zone3区域对推力有贡献，但zone4区呈现出阻力特性.zone5是鱼体尾尖附近的区域，振荡最强，其周期均值最大，具有大阻力特征.(3)尾迹区在摆动中大部分时间里都呈较弱的振荡，其周期均值为9.31%，体现的是阻力，远小于zone5的作用.由此可见，流场中的流动结构带来的阻力振荡也主要来自于尾部附近.为了更好地展示尾迹区的旋涡贡献，图8 给出了尾流区（wakelw a k e 5）的Q对推力贡献随时间的变化.结果表明，每一个分区对推力的贡献都是负的，越远离鱼尾的尾迹区域，曲线的涨落就越小，均值也越小.曲线的涨落主要是由于旋涡在区域内流进流出引起的，而远离鱼体的尾迹区的均值减小是由于虚速度势函数是随距离的衰减函数.整体来看，尾迹区对鱼体受力的影响是较弱的.本小节讨论了典型波状推进情况下的流体力分0.0200.02T.Q,wakel-0.04T.Q.wake2T.,wake3-0.06T.0,wake4T.Q.wakes00.20.40.60.81.0t/T图8 不同尾迹区的Q分布对推力的贡献Fig.8 The time-dependent contribution of Q distribution to thrust indifferentwake regions力14462023年第55卷报学学解.研究发现，鱼体波状摆动产生的推进力，其正贡献主要来自于边界加速度的瞬时响应，而流场中的流体旋转和应变速率相对大小对应的是阻力，在科类的摆尾推进中，两者都主要体现在尾部附近的流动.除此之外，边界上的黏性耗散效应主要体现在摩阻上.2.2推力与St数关系科鱼类通过改变鱼体波动的频率来实现推进的改变，这意味着推力系数将随St数改变.图9给出了不同雷诺数Re条件下按照式（9）计算的4部分时均值CT.a),CT.o),和随运动参数St的变化规律，其中频率从0.5增加到4.0.图中数据显示，CT,a)与St呈明显的二次关系，且与雷诺数无关；其余3项虽然受到雷诺数影响，但与St的关系也近似为二次方关系，CT.o)和(CT,)随着St增加而体现出阻力增长.由于 CT,v)和(CT,)与壁面摩阻有关,所以图9(c)和图9(d)中数据的绝对值随着雷诺数增大而减小.但在图9(b)中，CT.o)随雷诺数变化的规律有所不同.在鱼体高频波动条件下(即大St数下)，雷诺数越大，CT.o)的阻力系数减小,这与(CT.r)的规律一致.而在低频条件下,雷诺数越大，的阻力系数增大.1.00.Re=500Re=10000.80Re=2000+Re=50000.600.400.2000.1 0.20.30.40.50.60.7 0.8St(a)边界加速度对推力的贡献(a)The contribution of boundary acceleration to thrust-0.05Re=500-0.07Re=1000Re=2000-0.10Re=5000-0.12-0.15-0.17-0.20-0.230.10.20.30.40.50.60.70.8St(b)流场中Q分布对推力的贡献(b)The contribution of distribution to thrust in flow field0.10Re=500Re=10000.08Re=20000.06Re=50000.040.020-0.020.10.20.30.40.50.6 0.70.8St(c)壁面类摩阻分量(c)The wall friction-like component-0.05-0.10-0.15-0.20-0.25-0.30Re=500-0.35Re=1000Re=2000-0.40Re=50000.10.20.30.40.5 0.60.70.8St(d)壁面摩阻分量(d)The wall friction component图9不同雷诺数Re条件下，和(Ct,r)随运动参数St的变化规律Fig.9The time-averaged thrust coefficient(Cr.a),(Cr,o),(Cr.v),as a function of S t at different Reynolds number为了更好地显示4部分时均值对推力贡献的相对大小，图1 0 给出了Re=500和Re=5000情况下的各项贡献的堆叠柱状图.从图中可以看出，无论是低雷诺数还是高雷诺数情况，正推力的贡献都主要来自于边界变速运动.壁面切应力的类摩阻分量无论产生正贡献还是负贡献，相对于其他3项总是较小的.阻力分量主要来源于流场中Q的贡献和壁面摩阻的贡献.对于低雷诺数（Re=500）的情况，两者的阻力贡献相当，摩阻的阻力贡献略大于流场中Q的贡献.对于高雷诺数（Re=5000）的情况，流体中Q分布的阻力贡献强于摩阻的贡献.下面我们将详细探讨边界加速度项 CT,a）为什么与St成平方关系.根据式(1 4),边界加速度可以看作与振幅和频率的平方成正比，即y/atA.f，这只与波动有关，而鱼体上下表面的虚速度势差值d(x)与鱼体厚度分布和摆动幅值有关.在鱼体厚度分布不变的情况下，需考虑鱼尾摆动幅度改变对其值的影响.图1 1(a)给出了不同摆幅条件下（A=0.06,0.08,0.10,0.12,0.14)，在t=0T,t=T/8,t=T/4和t=3T/8时刻对应的d(x)沿着x的分布，当摆尾幅值A越大，对应的越大.图1 1(b)给出了p/A的分布1447黄恺俊等：基于虚功率原理的科鱼类波状推进分析第7 期1.000.750.500.250-0.25-0.500.20.40.60.81.01.21.41.6St(a)Re=5001.0（KC0.80.60.40.200.20.20.4 0.60.81.01.21.41.6St(b)Re=5000图1 0 不同雷诺数Re条件下，CT.a),(CT.),(CT,)和(CT,）的贡献堆叠柱状图Fig.10 Stacked histogram of the time-averaged thrust coefficient(CT.a),(CT,o),(Cr.v),(Cr.,)at different Reynolds number曲线几乎重合，这表明(x)与摆幅A成正比，即,d(x)A.因此，对于任意时刻，FT.a沿着鱼体中线的分布都与A?f?成正比，由此可得CT.aST2,CT.a)SP2(15)即由边界加速度得到的瞬时推力和推力均值都与St成正比，并且图4(b)给出的分区贡献占比不会随着运动参数St的变化而变化.0.200.150.100.500-0.05-0.100.06-0.150:0810810.-0.2001.00.10.80.20.60.30.4/70.40.2x/L0.50(a)g (x)0.150.100.500-0.05A1-0.06A0:080.100A8.12A1=8.141.00.1A0.80.20.60.30.4t/70.40.2X/L0.50(b)(x)/A图1 1 在t=0T,t=T/8,t=T/4和t=3T/8时刻以及不同摆幅条件下（A=0.0 6,0.0 8,0.1 0,0.1 2,0.1 4)，(x)以及d(x)/A沿着鱼体中线的分布Fig.11 Distribution of A(x)and Ao(x)/A along the center line of fishat different time(t=OT,t=T/8,t=T/4,t=3T/8)and under differentamplitude condition(A=0.06,0.08,0.10,0.12,0.14)进一步关于时均值 CT,a),CT,Q),(CT,v)和(CT,f)的分析可以结合已经发表的关于鱼体波状推进标度律1 5 进行分析.2.3波状推进的标度律分析Yu等1 5 基于涡动力学分析方法，给出了仿鱼类波状摆动推进器的推力系数与St和Re之间的标度关系，即)s2(1-)-Red2d3Cr)=(di-(16)VRe时均推力正比于St?(1-U/c)，比例系数有一部分和雷诺数无关，另一部分正比于雷诺数的-1/2 次方(1/V