申请/专利权人：西北工业大学

申请日：2024-05-09

公开（公告）日：2024-07-05

公开（公告）号：CN118296737A

主分类号：G06F30/15

分类号：G06F30/15;G06F30/20;G06F119/14

优先权：

专利状态码：在审-实质审查的生效

法律状态：2024.07.23#实质审查的生效;2024.07.05#公开

摘要：本发明公开了一种扑动翼的设计方法，属于涉及扑翼飞行器技术领域。扑动翼为仿生扑动翼，仿生扑动翼包括臂翼与手翼，臂翼为与鸟类翅膀的肱骨相对应的部分，手翼为与鸟类翅膀的尺骨和桡骨相对应的部分，方法包括以下步骤：a.确定扑动翼的初始设计参数组合；b.确定扑动翼的最优平面形状，令其为初始扑动翼平面形状；c.确定最大推力系数所需要的最佳变形；d.确定初始扑动翼的变形扭转角；e.确定扑动翼的最优刚度分布；f.确定扑动翼的最终扑动频率、巡航速度与迎角；g.通过全机风洞实验判断扑动翼设计是否符合设计指标，符合设计指标则设计流程结束；否则，进行迭代。它提高了扑动翼的设计迭代速度，加快了扑翼飞行器的研制进程。

主权项：1.一种扑动翼的设计方法，所述扑动翼为仿生扑动翼，所述仿生扑动翼包括臂翼与手翼，所述臂翼为与鸟类翅膀的肱骨相对应的部分，所述手翼为与鸟类翅膀的尺骨和桡骨相对应的部分，其特征在于：所述方法包括以下步骤：a.确定扑动翼的初始设计参数组合：a.依据扑翼飞行器的设计指标，根据尺度律公式，确定扑动翼的估算翼展b、翼面积S、扑动频率f，所述扑翼飞行器的设计指标为：起飞总重m和巡航速度Vcruise： 式中：m为起飞总重；b为翼展；S为翼面积；f为扑动频率；b.依据步骤a中所获取的扑动频率f及所述设计指标中的的巡航速度Vcruise，根据斯特劳哈尔数相似原则获取扑动幅度A： 式中：A为扑动幅度；St为斯特劳哈尔数，所述斯特劳哈尔数St与扑动翼的净推力正相关，扑翼飞行器在斯特劳哈尔数St为0.2～0.4时具有最高的推进效率，且St大于0.3时，平均功率将大幅增大，因此扑动幅度计算时，St取0.3以兼顾推力大小和推进效率；b.确定扑动翼的最优平面形状，令其为初始扑动翼平面形状：基于步骤a中得到的扑动幅度，以其作为初始设计点，进行扑动翼平面形状的优化设计，得到最优的扑动翼平面形状参数，所述最优平面形状参数包括翼面积S，展弦比AR，梢根比λ，臂翼占比1b以及控制手翼后缘形状的控制点坐标x，y，优化模型如下：minD＝fdesign_parameterss.t.L≥m0gs.t.b≤b0s.t.P≤P0design_parameters：S，AR，λ，lb，x，y式中：D为扑动翼所产生的阻力，其为设计变量的函数，函数关系由计算流体力学CFD或者风洞实验获得；L为扑动翼所产生的升力，其为设计变量的函数，函数关系由计算流体力学CFD或者风洞实验获得，约束为超过扑翼飞行器的起飞重量m0g；b为扑动翼的展长，约束为不超过扑翼飞行器的展长指标b0；P为扑动翼的气动功耗，约束为不超过扑翼飞行器的功耗指标P0；S为扑动翼的翼面积，AR为扑动翼的展弦比，λ为扑动翼的梢根比，臂翼占比1b，控制手翼后缘形状的控制点坐标x，y；c.确定最大推力系数所需要的最佳变形：扑动翼为骨架与膜的构型形式，其变形可等效为不同站位的二维翼型的扭转变形，因此，可通过二维翼型的沉浮扭转运动来模拟三维柔性翼的扑动；进而，通过CFD计算，确定最大推力系数CT与二维翼型的扭转变形角θ之间的关系，得到最大推力所需要的最佳变形θoptimal：h＝hmcos2πft θoptimal＝gf，hm，A，Vcruise式中：h为沉浮位移，用来模拟三维挥舞运动；hm为沉浮幅度，用来模拟不同展向站位以及扑动幅度的影响；θ为瞬时扭转角，用来模拟扑动翼的柔性；θn为扭转幅度，用来模拟不同的柔性变形；θoptimal为最大推力系数时所对应的最佳扭转角，它通过CFD计算上述沉浮扭转运动得出，是扑动频率f、沉浮幅度hm、扑动幅度A以及巡航速度Vcruise的函数；d.确定初始扑动翼的变形扭转角：通过数字图像相关技术DIC以及高速摄像机对初始扑动翼在扑动过程中的变形扭转角进行测量，在扑动翼表面粘贴标记点，对高速摄像机进行标定，之后对扑动过程进行拍摄以及标记点坐标x，y，z的输出；最终，对标记点坐标进行后处理得到初始扑动翼的变形θbase情况；a.通过DIC技术和高速摄像机测量扑动翼各标记点坐标：x0，y0，z0x1，y1，z1x2，y2，z2……xi，yi，zi式中：xi，yi，zi为扑动翼各标记点空间坐标，为扑动翼扑动角度的函数；b.根据扑动翼靠近翼根处三标记点确定扑动翼参考平面： aX+bY+cZ＝d式中：a，b，c，d分别为控制参考平面的系数，为扑动翼扑动角度的函数；c.获取扑动翼上其余各标记点距离参考平面的距离： d.获取各标记点扭转角，即为初始变形θbase：θbase＝arcsinLiclocal式中：Li为各标记点距离参考平面的距离，为扑动翼扑动角度的函数；clocal为标记点当地的弦长；e.确定扑动翼的最优刚度分布：通过对比初始变形θbase与最佳变形θoptimal之间的关系，对刚度分布进行设计，使得改进后的扑动翼变形θdesign逼近最佳变形θoptimal：θdesign＝fD，M式中：θdesign为设计后的扑动翼变形，其为刚度分布参数的函数，通过步骤d进行测量；D为前缘碳杆直径，刚度分布参数之一；M为翼膜材料，刚度分布参数之一；f.确定扑动翼的最终扑动频率、巡航速度与迎角：对经过上述步骤设计完成得到的扑动翼进行风洞实验，得到不同巡航速度下，扑动翼升推力以及功耗关于扑动频率和迎角的函数关系，求解升推力平衡方程，以功耗最低为目标，确定扑动翼的最终扑动频率、巡航速度与迎角：a.在不同的巡航速度下，风洞实验测量扑动翼升推力以及功耗关于扑动频率和迎角的函数关系：givenVcruiseiL＝gf，αT＝hf，αP＝sf，αb.在不同巡航速度下，求解升推力平衡方程，得到配平点迎角与扑动频率： 式中：ftrim为配平点扑动频率，αtrim为配平点迎角，Ptrim为配平点功耗；c.比较不同巡航速度下的配平点功耗Ptrim大小，得到功耗最低时的巡航速度以及对应的配平点频率ftrim与配平点迎角αtrim；g.通过全机风洞实验判断扑动翼设计是否符合设计指标，符合设计指标则设计流程结束；否则，重复步骤b～步骤f进行迭代，直至符合设计指标。

全文数据：

权利要求：

百度查询： 西北工业大学 一种扑动翼的设计方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD