翼型非对称拍动的水动力学研讨.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 翼型非对称拍动的水动力学研究木 胡文蓉 上海交通大学工程力学系，200240 E·mail：wrhu@sjtu．edu．cn 摘要： 通过求解不可压缩N—s方程模拟了翼型的非对称俯仰和沉浮运动．通过对升阻力 系数，流场涡结构和压强速度的分布情况的分析，探讨了各种非对称运动对流动结构和力学 性能的影响，并分析了柔性变形在非对称运动中的作用。 关键词：翼型；非对称拍动；非定常流动 l引言 鱼类常通过鳍的拍动来实现推进和机动运动。而机动运动则往往是靠鳍的非对称拍动来 产生较大的推进力和升力，从而实现快速机动。如果能弄清鱼鳍非对称拍动的水动力学特征， 则有助于更清楚地认识鱼类机动运动的机理，且能为仿生机器的研制提供理论依据。 关于翼型的数值和实验研究基本上是针对对称的简谐振荡而进行的，其中包括对振荡频 率瓯拍动幅度，最大攻角等参数对力学性能的影响的研究Il训。但目前关于翼型的非对称运 动的研究还很少见【5J，翼型的非对称运动的水动力学特性还并不清楚。 与对称运动相比较。通过对升阻力系数，流场涡结构和压强速度的分布情况的分析，探讨了 各种非对称运动对流动结构和力学性能的影响，-并分析了柔性变形在非对称运动中的作用。 2物理模型和数值方法 本研究从守恒型的非定常不可压缩流动的N．S方程出发，采用人T伪压缩方法求解。由 于同体边界的运动，物理平面内的求解域(包括儿何边界)及其网格都随时间改变。此时在 控制方程中也相应地引入网格速度，使之满足空间守恒定律，否则将可能引起人工质量的积累 或耗散。对于控制体n其无量纲方程为 ‘由国家自然科学基金(No和上海市重点学科建设项目(项目编号：B206)资助 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 矿￡(安)叭昙r￡四y+叮cM啪… (2) 其中／=(，+E，G+q)，月=(以，甩，)为外法向单位矢量，％为网格速度。 甜2+p Zn， 励 Q=[习口=[三]，= G恪卜去H睁去M ∥为伪压缩性系数；p为压强；“，’，为笛卡尔坐标系(x，y)中的速度分量；f为物理时间；f 为虚拟时间。 方程用有限体积法进行离散。对流项采用三阶精度的MUSCL格式：黏性项利用有限体积 法中高斯积分进行计算；在时间方向则采用欧拉隐式格式进行离散。对于变形边界，为了满 足空间守恒定律，在控制方程中引入了网格速度。并采用插值的方法生成适时动态网格。 本研究针对NACA0012翼型的非对称俯仰和沉浮运动进行了研究。仅当上挥和下拍运动 时间不对称时，翼型的非对称运动方式如下： f 厅：一竺sinf，丝]+竺f O≤f≤S丁 I 27【 Lsr／s丁 STstsT 卜挚(警詈)一尚¨。+滞 式中办为俯仰角度或沉浮运动y坐标，盯为上挥时间，△办为拍动幅度。当s=O．5时，则为 对称振荡运动。 3结果与讨论 3．1上挥下拍时间不对称 图l至图4分别给出了俯仰和沉浮运动中，上挥下拍时间不对称所引起的涡量压强分布 及其对力学性能的影响。 对于俯仰运动，当上挥下拍时间不对称时，能增人推力，不对称运动对升力也有着显著 的影响，上挥时间短，产生负升力，引起顺时针力矩，反之亦然。越不对称推力和升力都越 显著，当上挥下拍时间相差O．2丁时，能增加推力50％左右。上挥下拍时间不对称主要靠增强 附加惯性效应增加推力和升力，并且靠尾迹中旋涡的合并增强尾迹中压强，同时还增强了尾 迹射流，从而增加推力(图1和图2)。 对于沉浮运动，与对称沉浮运动相比，当上挥下拍时间不对称时，并非都能增加推力， 且并非是上挥时间短就一定产生负升力。当运动不对称性较弱时(如萨O．4或肛0．6)，推力都 略有下降，且上挥时间短产生正升力，下拍时间短产生负升力。而运动不对称性较强时，能