关键词:褶皱;柔性;蜻蜓前翅;流固耦合;滑翔 Abstract Dragonflies are capable of carrying out dramatic flight manoeuvres, gliding flight is a common mode of flight for dragonfly, and dragonfly wings are the source of dragonflies dramatic flight manoeuvres. Unlike typical engineered airfoil, dragonfly wings are not smooth, wing cross-section are highly corrugated. It has been shown that corrugations could enhance the spanwise stiffness in the wings and influence aerodynamic performance of the dragonfly wings. Flexibility is another characteristic of the dragonfly's wings, which is mainly manifested as the flexible deformation of the wings during the flight. To explore corrugations and flexibility effect on aerodynamic performance of the dragonfly forewings in gliding flight, a computational fluid dynamics (CFD) model and a computational structural mechanics (CSD) model of the corrugated dragonfly forewing are established based on current research, and the modal analysis verified that the model has sufficient accuracy. The corrugated rigid and flexible dragonfly forewing are acquired by using CFD method and CFD/CSD coupling method respectively. The simulation indicated that flexible and corrugated forewings is subjected to aerodynamic load, which only produces bending deformation without torsion deformation in gliding flight, and the aerodynamic response time is short. Compared with the aerodynamic performance of the rigid and corrugated forewings, the result showed that veins and cuticular membrane of flexible forewing are deformed which caused the lift coefficient and drag coefficient decrease, the leading edge vortex of flexible forewing is much higher than rigid forewings at large angle of attack because of deformed vein. The aerodynamic performance of rigid forewings is better below 10 degree angle of attack,and the aerodynamic performance of flexible forewings is better at large angle of attack as result.
Keywords:corrugation;flexibility;dragonfly forewing;fuid-structure interaction;gliding flight -->0 PDF (19491KB)元数据多维度评价相关文章收藏文章 本文引用格式导出EndNoteRisBibtex收藏本文--> 刘惠祥, 何国毅, 王琦. 蜻蜓滑翔时柔性褶皱前翅气动特性分析1)[J]. 力学学报, 2019, 51(1): 94-102 https://doi.org/10.6052/0459-1879-18-157 LiuHuixiang, HeGuoyi, WangQi. NUMERICAL STUDY ON THE AERODYNAMIC PERFORMANCE OF THEFLEXIBLE AND CORRUGATED FOREWING OF DRAGONFLY IN GILDINGFLIGHT1)[J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(1): 94-102 https://doi.org/10.6052/0459-1879-18-157 微型飞行器低雷诺数飞行时,由于机翼表面的层流附面层难以维持,气流不易稳定,导致其存在失速迎角小和抗阵风性能差等缺点[1-6],而昆虫飞行时能克服低雷诺数($Re$)的气动局限性,实现高稳定性和低能耗的飞行.研究表明,昆虫柔性翅膀的变形是其用来克服低雷诺数飞行时气动局限性的主要方式之一[7-9].有昆虫界战斗机之称的蜻蜓,作为自然界最优秀的飞行者之一[10],其翅膀上没有神经和肌肉,所以蜻蜓对翅膀的主动控制相对较少,在飞行过程中的蜻蜓翅膀的变形是被动的,变形的大小取决于翅膀本身的结构和材料的属性以及所受气动力的分布情况[11],研究蜻蜓柔性翅膀的变形对其气动特性的影响,对微型扑翼飞行器的设计具有一定的指导 意义. 蜻蜓飞行过程中柔性翅膀的变形是一个流固耦合问题[12],目前国内外对蜻蜓柔性翅膀进行了相关的研究,Hamamoto等[13-14]基于有限元分析方法处理蜻蜓悬停时的流固耦合问题,结果表明柔性翅膀的变形对其气动特性具有较小的影响.孟令兵等[15]基于计算流体力学/计算结构力学(computational fluiddynamics,CFD/computational structuralmechanics,CSD)方法对蜻蜓扑翼飞行时柔性翅膀的变形进行了数值模拟,结果显示柔性翅膀的变形降低了阻力,并增大推力,提高了蜻蜓的飞行性能,但其建立的模型实际上为具有蜻蜓翅膀外轮廓的平板.郝淑文[16]利用Matlab软件建立了适合于求解柔性膜扑翼气动特性的流固耦合整体法求解程序,分析了刚性和柔性蜻蜓翅膀模型前飞时动态变形的瞬时、时均气动响应,但此翅膀模型忽略了翅脉、翅膜以及褶皱等重要结构.以上研究只考虑了柔性翅膀的变形,而褶皱是蜻蜓翅膀上最为显著的结构之一,文献[11,17-18]对蜻蜓褶皱翼的研究表明,与平板相比,褶皱翼有比平板具有更高的升阻比;而Rees[19]、Meng等[20]对蜻蜓褶皱翼的研究则表明,褶皱翼的升阻比低于平板;Zhang等[21]的研究结果显示,蜻蜓褶皱翼的升阻比低于平板,通过改变褶皱的起皱幅度可显著改善其气动效能,使褶皱翼的升阻比高于平板.在蜻蜓实际飞行中,皱褶与柔性翅膀的变形同时存在,两者共同作用下对蜻蜓翅膀气动特性的影响如何?本文基于真实的蜻蜓前翅,根据前人对蜻蜓前翅的实际测量数据和研究成果,经过适当的简化,构建出有褶皱结构的蜻蜓三维前翅CFD模型和CSD模型,基于CFD/CSD双向耦合的计算方法,对蜻蜓三维柔性褶皱前翅的变形问题展开研究,以期揭示蜻蜓滑翔时柔性对褶皱前翅气动特性的影响.
利用三维建模软件CATIA,通过逆向工程的方法,将高清蜻蜓前翅数码图片导入CATIA,利用样条曲线拟合纵横交错的翅脉,建立蜻蜓前翅的轮廓,然后填充轮廓形成翅膜曲面,根据实际测量皱褶结构深度和夹角的数据,偏转与之相对应翅膜面,翅脉和翅膜的空间分布形成前翅上的皱褶结构,赋予翅膜厚度和翅脉的形状与尺寸,翅膜镶嵌在圆管状翅脉上,最终建立蜻蜓三维褶皱前翅的CFD模型和CSD模型分别如图3和图4所示.有限元分析时,翅脉和翅膜采用均质的实体单元,因此前翅CSD模型的翅脉是由厚度为0.035mm的空心圆管构成,而CFD模型的翅脉为实心圆管,其他尺寸参数以及空间构型与前翅的CFD模型一致. 显示原图|下载原图ZIP|生成PPT 图3蜻蜓褶皱前翅CFD模型... -->Fig. 3CFD model of the corrugated dragonfly forewing... -->
显示原图|下载原图ZIP|生成PPT 图4蜻蜓褶皱前翅CSD模型... -->Fig. 4CSD model of the corrugated dragonfly forewing... -->
1.3 模型的验证
为了验证蜻蜓三维褶皱前翅模型的准确性,对此模型进行模态分析,与弯艳玲[38]通过三维激光扫描技术建立的蜻蜓三维褶皱前翅模型进行对比.真实前翅翼根处的翅脉和翅膜是固定在蜻蜓躯体上的,因此在模态分析时,将前翅翼根处的翅脉和翅膜设置为固定约束,如图5所示,所有在翼根处的网格节点被设置为固定约束,约束此网格节点的所有自由度.蜻蜓前翅CSD模型有限元网格如图5所示,翅脉和翅膜都采用均质的实体单元. 显示原图|下载原图ZIP|生成PPT 图5有限元分析边界条件及网格... -->Fig. 5Boundary conditions and grids of finite element analysis... -->
模态分析的结果如表1所示,结果表明蜻蜓前翅的主要模态(一阶模态,二阶模态和三阶模态)与弯艳玲的结果差别不大,因弯艳玲建立的蜻蜓前翅是通过激光扫描系统建立的精确模型,对比本文建立的蜻蜓三维褶皱前翅模型前三阶模态频率,其最大相对误差在10%以下,说明建立的此模型能保证足够的精度. Table 1 表1 表1蜻蜓褶皱前翅频率对比 Table 1Frequency comparison of the corrugated dragonfly forewing
建立外流域,前翅的翼根处所在的平面为对称面,前翅的边界条件为壁面,湍流模型采用的是$k$-$\omega$ SST低雷诺数模型.采用非结构网格离散计算域,对于非结构网格,目前常用的网格形式有四面体网格、八面体网格以及笛卡尔网格,在相同条件下,笛卡尔的网格数量最少,其次为八面体网格,但由于笛卡尔网格无法适用于网格变形技术--动网格技术,因此用于柔性前翅的CFD网格采用适用动网格技术且高效的八面体网格对计算域进行离散.为使计算的结果不依赖于网格选取,对蜻蜓前翅CFD模型分别设置了3套不同大小的计算网格:第一套网格的蜻蜓前翅表面网格尺寸为0.04 mm;第二套网格为 0.03mm;第三套网格为0.02 mm.计算结果表明,使用第一套和第二套网格计算得到蜻蜓前翅的升力系数相对误差为14.9%,第二套和第三套网格的相对误差为4.7%.在权衡计算的准确性和计算时间后,最终选用第二套网格,如图7所示,即蜻蜓前翅尺寸为0.03mm. 固体求解时有限元网格如图5所示,蜻蜓前翅的表面为流固耦合交界面. 显示原图|下载原图ZIP|生成PPT 图7前翅周围网格局部图... -->Fig. 7Grids of the corrugateddragonfly forewing... -->
1.6 算例验证
选用第二套网格,计算了平板在$Re$=10 000时攻角($\alpha)$介于0$^\circ$~10$^\circ$(间隔2$^\circ$)下的升力系数($C_{l}$),与Kesel的实验结果对比如图8所示.从图8中可以看出,平板的升力系数数值计算值与实验值基本相近,表明此网格能够保证足够的精度. 显示原图|下载原图ZIP|生成PPT 图8数值计算平板的升力系数与Kesel[17]的实验结果对比图... -->Fig. 8Comparison between the numerical and the experimental lift coefficient of Kesel[17] ... -->
对比刚性和柔性蜻蜓褶皱前翅在滑翔时不同攻角下的气动特性差异,气动特性的差异主要根据升力系数($C_{ l}$)、阻力系数($C_{ d}$)和升阻比($C_{ l}/C_{d}$)来判定,刚性和柔性前翅的气动特性比较如表2所示,其中负值表示柔性前翅比刚性前翅的气动特性差,正值表示柔性前翅的的气动特性优于刚性前翅.在0$^\circ\!\!\sim$25$^\circ$攻角范围内,与刚性前翅相比,柔性减小了前翅的升力系数,同时也减小了阻力系数. Table 2 表2 表2柔性前翅和刚性前翅气动特性对比 Table 2Aerodynamic performance comparison of the rigid and flexible forewing
a
0。
5。
10。
15。
20。
25。
Cl
-2.98%
-5.21%
-1.83%
-0.48%
-0.20%
-0.19%
Cd
-0.58%
-1.86%
-2.21%
-3.78%
-3.53%
-7.15%
Cl/Cd
-4.34%
-3.53%
0.41%
3.32%
3.34%
3.15%
新窗口打开 与刚性前翅相比,随着攻角的增大,柔性前翅升力系数的减少量逐渐减小,而阻力系数则相反,柔性前翅阻力系数的减少量随攻角的增加逐渐增大.在10$^\circ$攻角以下时,刚性前翅的升阻比高于柔性前翅,而当攻角继续增大时,则柔性前翅的升阻比高于刚性前翅.因此攻角的变化决定了柔性能否改善蜻蜓前翅的气动特性,只有在大攻角(大于10$^\circ$)时,蜻蜓的柔性翅膀有助于改善其气动特性. 如图9所示,蜻蜓前翅褶皱被许多涡填充,褶皱上表面凹槽内的驻留涡沿顺时针旋转,下表面凹槽内的驻留涡沿逆时针旋转,这与Rees[19]的实验结果一致,且涡主要分布在蜻蜓前翅前缘区域起伏较大的褶皱处,因此褶皱结构能够滞留涡. 显示原图|下载原图ZIP|生成PPT 图9攻角为0$^\circ$时柔性前翅的涡量图及流线图... -->Fig. 9The vortex and streamline of flexible forewing at 0$^\circ$angles of attack... -->
如图10所示,对比刚性和柔性褶皱前翅的压力云图,柔性前翅的翅膜由于气动力的作用,其表面变得凹凸不平,上表面的低压区出现相对更高的压力分布,下表面高压区出现相对更低的压力分布,对比图10(d)和图10(h),攻角为25$^\circ$时这种差别更为明显.翅膜的上下表面压力差相较于刚性前翅的翅膜有所降低,更低的压力差意味着翅膀所受的合力更小,合力的分量之一--阻力也更小,因此柔性前翅的阻力系数比刚性前翅的低,攻角为25$^\circ$时阻力系数的减小量比攻角为15$^\circ$时的更大. 柔性前翅翅膜的上下表面压力差相较于刚性前翅的翅膜有所降低,更低的压力差意味着翅膀所受的合力更小,合力的分量之一--升力也更小,数值模拟的结果也显示在0$^\circ\sim$25$^\circ$攻角时,柔性前翅的升力系数比刚性前翅小,但从表2中攻角为5$^\circ$时柔性前翅比刚性前翅的升力系数减小5.21%,到攻角为25$^\circ$时柔性前翅比刚性前翅的升力系数减小0.19%,这是因为与攻角为5$^\circ$时蜻蜓柔性前翅的前缘涡相比,攻角为25$^\circ$时,柔性前翅的前缘翅脉处存在比刚性前翅更强的前缘涡,如图11所示.刚性前翅的前缘脉为圆管状,而柔性前翅的前缘脉产生变形,成类圆管状,如图11(d)所示,变形后的前缘脉诱导出更强的前缘涡.攻角为25$^\circ$时,前缘涡的存在使蜻蜓柔性前翅能够保持与刚性前翅相当的升力系数,但由于柔性前翅具有更低的阻力系数,相较于刚性蜻蜓前翅,柔性使得前翅气动特性得到改善,升阻比提高3.15%. 显示原图|下载原图ZIP|生成PPT 图10不同攻角时不同属性的褶皱前翅翅尖处压力云图... -->Fig. 10The pressure of corrugated forewing on wingtipwith different properties at different angles of attack... -->
显示原图|下载原图ZIP|生成PPT 图11不同攻角时不同属性60%展长处横截面的涡量图对比... -->Fig. 11The vortex of the cross sections with differentproperties at different angles of attack shown below at 0.6$l_{rel}$, where $l_{ rel}$ is the relative span length... -->
2.2 蜻蜓前翅的变形与应力
基于CFD/CSD的双向流固耦合计算方法,数值求解$Re$=2500时0$^\circ\!\!\sim$25$^\circ$攻角(间隔5$^\circ$)蜻蜓滑翔时前翅在流场中受气动载荷后的变形和应力.蜻蜓柔性褶皱前翅在0$^\circ\!\!\!\sim$25$^\circ$攻角(间隔5$^\circ$)时最大变形值依次为:0.1694,1.265,2.169,2.602,2.721,2.910,单位为mm;柔性褶皱前翅最大应力值依次为:1.602,9.537, 16.24,19.76,20.82,22.34,单位为MPa.变形和应力总体上随着攻角的增大而增大,攻角10$^\circ$以后,变形和应力增长率减小,特别是攻角在15$^\circ$和25$^\circ$之间,变形和应力增大,但增幅不大,此攻角范围内褶皱前翅在流场中表现为升力系数变化也不明显,与此时蜻蜓前翅出现大范围的流动分离有关. 如图12所示,在攻角为0$^\circ$时蜻蜓柔性褶皱前翅的最大变形位于靠近后缘的翅尖处,而攻角为5$^\circ\sim$25$^\circ$时前翅的最大变形位于翅尖处,柔性只是增加了垂直于蜻蜓前翅平面方向上的变形,但并没有使前翅发生扭转,即柔性对蜻蜓褶皱前翅的攻角的影响很小.如图10和图11所示,翅脉和翅膜因空气作用力产生变形,翅膜受气动载荷后表面成凹凸不平状,从而改变了蜻蜓柔性前翅的上下表面压力差.翅脉受气动载荷后由圆管状变成类圆管状,变形后的前缘翅脉,改变了柔性前翅的前缘涡强度. 显示原图|下载原图ZIP|生成PPT 图12不同攻角($\alpha )$下柔性前翅的变形(单位:m) ... -->Fig. 12The deformation(unit: m) of flexible dragonfly forewingat different angles of attack ($\alpha )$ ... -->
蜻蜓柔性褶皱前翅最大应力位置在翼根处,如图13所示,攻角为0$^\circ$时最大应力位置位于第五根横脉上表面,而攻角为5$^\circ$~25$^\circ$时蜻蜓前翅的最大应力位于第三根横脉上表面.随着攻角的逐渐增大,应力的分布发生变化,在第三根横脉和第五根横脉的根部出现了最大应力,而第三根横脉和第五根横脉都处于蜻蜓前翅褶皱的的上凸部分,应力主要集中于此区域,因此褶皱上凸部分承受蜻蜓滑翔时前翅的主要载荷. 显示原图|下载原图ZIP|生成PPT 图13不同攻角($\alpha )$下柔性前翅翼根部的应力(单位:Pa)云图... -->Fig. 13The stress(unit: Pa) of flexible forewing at differentanglesof attack ($\alpha )$ ... -->
2.3 蜻蜓前翅的气动响应
首先数值模拟蜻蜓刚性褶皱前翅的流场,在刚性前翅周围流场稳定后,改变前翅的属性,赋予前翅材料,基于CFD/CSD双向流固耦合数值计算开始,此时前翅将产生变形,刚性褶皱前翅的数值模拟相当于为流固耦合加载了一个稳定的初始流场.如图14(a)所示,在蜻蜓前翅周围流场稳定后,褶皱前翅的升力系数趋于稳定,其升力系数为0.1247,此时前翅属性为刚性.当将前翅属性由刚性改变为柔性时,在受到气动载荷后,蜻蜓柔性前翅出现振动,如图14(b)所示,在变形响应曲线中的表现为变形出现震荡,但震 荡的幅值逐渐减小最终稳定,前翅的升力 显示原图|下载原图ZIP|生成PPT 图14攻角0$^\circ$蜻蜓前翅气动响应图... -->Fig. 14The pneumatic response of corrugated dragonfly forewingat0$^\circ$ angles of attack... -->
蜻蜓在滑翔飞行时,三维褶皱刚性前翅和柔性前翅相比,在小攻角(小于10$^\circ$)时刚性前翅的气动特性优于柔性前翅,继续增大攻角,柔性前翅的气动特性则优于刚性前翅.蜻蜓扑翼飞行时,翅膀的攻角变化幅度大且频率高,而蜻蜓柔性翅膀的气动特性同时还受攻角变化的影响,因此蜻蜓扑翼飞行时,柔性也是其不能忽略的特性. 蜻蜓滑翔过程中,褶皱结构能够滞留涡,柔性前翅在翅尖处产生了垂直于翅膀所在平面方向上的最大变形,而没有发生扭转,翼根处受到应力最大,褶皱上凸部分承受蜻蜓滑翔时前翅的主要载荷.蜻蜓柔性前翅受气动载荷后产生变形,表现为翅脉和翅膜的外形改变以及整个前翅弯曲变形.翅膜的外形改变,其上下表面压力差相较于刚性前翅减小了,从而柔性前翅的升力和阻力也减小了.但在大攻角时,柔性前翅由于外形改变的前缘脉,诱导出比刚性前翅更强的前缘涡,维持了柔性前翅的升力,使柔性前翅的升力减少量变小,因此大攻角时柔性对蜻蜓褶皱前翅的气动特性的改善约3%.蜻蜓褶皱前翅在受气动载荷后,气动响应时间短,表现出优异的力学性能.本研究有助于解释蜻蜓翅膀优秀的功能表现,对仿生扑翼飞行器的柔性机翼设计有一定的参考意义. The authors have declared that no competing interests exist.
(LiuQiang, LiuZhou, BaiPeng, et al.Numerical study about aerodynamic characteristics and flow field structures for a skin of airfoil with active oscillation at low Reynolds number . Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(2): 269-277 (in Chinese))URL [本文引用: 1]
(MengLingbing, AngHaisong, XiaoTianhang.Numerical simulation of fluid-structure interaction for flexible wing MAV . Journal of Nanjing University of Aeronautics & Astronautics, 2013, 45(5): 621-627 (in Chinese))
[3]
LianY, WeiS, ViieruD, et al.Membrane wing aerodynamics for micro air vehicles . Progress in Aerospace Sciences, 2003, 39(6): 425-465
[4]
AIAA. Computational modeling of highly flexible membrane wings in micro air vehicles . AIAA Paper2006-1661, 2006
[5]
JaimanRK, JiaoX, GeubellePH, et al.Conservative load transfer along curved fluid-solid interface with non-matching meshes . Journal of Computational Physics, 2006, 218(1): 372-397
[6]
LianY, Wei S. Three-dimensional fluid-structure interactions of a membrane wing for micro air vehicle applications . AIAA Paper2003-1726, 2003 [本文引用: 1]
[7]
MuellerTJ.Fixed and flapping wing dynamics for MAV applications . AIAA Progress in Astron and Aeron, Massachusetts: AIAA Press, 2001: 195 [本文引用: 1]
[8]
LiuH, KawachiK.A numerical study of insect flight . Journal of Computational Physics, 1998, 146(1): 124-156
[9]
RaynerJMV.A new approach to animal flight mechanics . Journal of Experimental Biology, 1979, 80(1): 17-54 [本文引用: 1]
[10]
HuZ, DengXY.Aerodynamic interaction between forewing and hindwing of a hovering dragonfly . Acta Mechanica Sinica, 2014, 30(6): 787-799URL [本文引用: 1]
[11]
KwokM, MittalR.Experimental investigation of the aerodynamics of a modeled dragonfly wing section . Bulletin of Faculty of Engineering Toyama University, 2005, 5(1): 50-53 [本文引用: 2]
[12]
NodaR, NakataT, LiuH.Effects of wing deformation on aerodynamic performance of a revolving insect wing . Acta Mechanica Sinica, 2014, 30(6): 819-827URL [本文引用: 1]
[13]
HamamotoM, OhtaY, HaraK, et al.Application of fluid--structure interaction analysis to flapping flight of insects with deformable wings . Advanced Robotics, 2007, 21(1-2): 1-21 [本文引用: 1]
[14]
HamamotoM, OhtaY, HaraK, et al.Free-flight analysis of dragonfly hovering by fluid--structure interaction analysis based on an arbitrary Lagrangian--Eulerian method . Advanced Robotics, 2013, 27(9): 657-666 [本文引用: 1]
(MengLingbing, AngHaisong, XiaoTianhang.Analysis of aerodynamic characteristics of a flexible wing of dragonfly based on CFD/CSD method . Journal of Aerospace Power, 2014, 29(9): 2063-2069 (in Chinese)) [本文引用: 1]
(HaoShuwen.Fluid-structure interaction of unsteady aerodynamic of membrane flapping wings at low reynolds number. [PhD Thesis] . Beijing: Beijing Institute of Technology, 2015 (in Chinese)) [本文引用: 1]
[17]
KeselAB.Aerodynamic characteristics of dragonfly wing sections compared with technical aerofoils . Journal of Experimental Biology, 2000, 203: 3125-3135 [本文引用: 3]
[18]
VargasA, MittalR, DongH.A computational study of the aerodynamic performance of a dragonfly wing section in gliding flight . Bioinspiration & Biomimetics, 2008, 3(2): 026004 [本文引用: 1]
[19]
ReesCJC.Form and function in corrugated insect wings . Nature, 1975, 256(5514): 200-203 [本文引用: 2]
[20]
MengXG, SunM.Aerodynamic effects of wing corrugation at gliding flight at low Reynolds numbers . Physics of Fluids, 2013, 25(7): 071905 [本文引用: 1]
[21]
ZhangZ, YinY, ZhongZ, et al.Aerodynamic performance of dragonfly wing with well-designed corrugated section in gliding flight . Computer Modeling in Engineering & Sciences, 2015, 109(3): 285-302 [本文引用: 1]
[22]
DomprehKA, EghanMJ, KotsediL, et al.Comprehensive optical study of the dragonfly Aeshna cyanea, transparent wing . Optics Communications, 2013, 297(12): 176-181 [本文引用: 1]
(ZhaoHongxiao, ZhongZheng.Research advance in mechanical research of dragonfly wings . Chinese Quarterly of Mechanics, 2009, 30(3): 398-404 (in Chinese)) [本文引用: 1]
[24]
OkamotoM, YasudaK, AzumaA.Aerodynamic characteristics of the wings and body of a dragonfly . Journal of Experimental Biology, 1996, 199(2): 281-294 [本文引用: 1]
(SunMao.Aerodynamics of inset flight . Chinese Journal of Theoretical and Applied Mechanics, 2015, 45(2): 1-28 (in Chinese)) [本文引用: 1]
[26]
Hertel H.Structure, Form, Movement. New York: Reinhold Publishing Corporation, 1966
[27]
WoottonRJ, NewmanADJS.An approach to the mechanics of pleating in dragonfly wings . Journal of Experimental Biology, 1986, 125(1): 361-372
[28]
NewmanDJS.The functional wing morphology of some. [PhD Thesis] . Odonata: University of Exeter, 1983
[29]
AntoniaBK, UteP, WernerN.Biomechanical aspects of the insect wing: an analysis using the finite element method. Computers in Biology and Medicine}.1998, 28: 423-437
[30]
陈应龙. 蜻蜓翅膀微结构与力学行为的仿生分析研究. [博士论文] . 北京:清华大学, 2012
(ChenYinglong.Bionicresearch on the relation between the multi-scaled structure and microstructure and mechanical behavior of dragonfly wing. [PhD Thesis] . Beijing: Tsinghua University, 2012 (in Chinese))
[31]
DongHB, ChristopherK, LiangZX, et al.An integrated analysis of a dragonfly in free flight//Proceedings of 28th AIAA Applied Aerodynamics Conference, 2010: 4390-4400 [本文引用: 1]
[32]
TongJ, ZhaoYR, SunJY, et al.Nanomechanical properties of the stigmaof dragonfly Anax parthenope julius Brauer . Journal of Materials Science, 2007, 42(8):2894-2898URL [本文引用: 1]
[33]
RenHH, WangXS, LiXD, et al.Effects of dragonfly wing structure on the dynamic per formances . Journal of Bionic Engineering, 2013, 10(1): 28-38URL
[34]
JongeriusSR, LentinkD.Structural analysis of a dragonfly wing . Experimental Mechanics, 2010, 50(9): 1323-1334URL [本文引用: 1]
[35]
SudoS, TsuyukiK, IkohagiT, et al.A study on the wing structure and flapping behavior of a dragonfly. JSME International Journal Series C Mechanical Systems, Machine Elements and Manufacturing, 1999, 42(3): 721-729 [本文引用: 2]
(LiXiujuan.Bionic investigation on mechanical mechanism of dragonfly wings functional characteristics. [PhD Thesis] . Jilin: Jilin University, 2013 (in Chinese)) [本文引用: 2]
[37]
ShahnaziR, ParizN, KamyadAV.Adaptive fuzzy output feedback control for a class of uncertain nonlinear systems with unknown backlash-like hysteresis . Communications in Nonlinear Science & Numerical Simulation, 2010, 15(8): 2206-2221URL [本文引用: 1]
(WanYanling.Dynamics and fatigue life of three-dimensional structure of dragonfly wings. [PhD Thesis] . Jilin: Jilin University, 2010 (in Chinese)) [本文引用: 2]