CN104401504A - 一种固定翼航测型无人机设计方法 - Google Patents

一种固定翼航测型无人机设计方法 Download PDF

Info

Publication number
CN104401504A
CN104401504A CN201410665303.5A CN201410665303A CN104401504A CN 104401504 A CN104401504 A CN 104401504A CN 201410665303 A CN201410665303 A CN 201410665303A CN 104401504 A CN104401504 A CN 104401504A
Authority
CN
China
Prior art keywords
wing
lift
zero
angle
coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410665303.5A
Other languages
English (en)
Other versions
CN104401504B (zh
Inventor
许伟
陈刚
吴鹏
赵茹玥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences
Original Assignee
China University of Geosciences
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences filed Critical China University of Geosciences
Priority to CN201410665303.5A priority Critical patent/CN104401504B/zh
Publication of CN104401504A publication Critical patent/CN104401504A/zh
Application granted granted Critical
Publication of CN104401504B publication Critical patent/CN104401504B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明一种固定翼航测型无人机设计方法,从机身结构和工作性能出发,设计出机翼翼型和机身,将设计出的参数通过Profili软件模型分析,在CAD中制做底图,将制作好的底图在SketchUp中进行3D建模,在得到载机3D模型后打印图纸进行手工制作,最后得到成本费用低、操作简单、稳定性强的XW-06航测型无人机载机。通过本发明提供的设计方法制作出的无人机,成本低、荷载能力强,能够针对不同的天气条件实现较强的抗干扰能力和稳定性,实现人性化的人机交互性。

Description

一种固定翼航测型无人机设计方法
技术领域
本发明属于航空摄影测量领域,具体涉及一种固定翼航测型无人机设计方法。 
背景技术
由于航测内业影像处理过程中,对影像质量(航向重叠度、旁向重叠度、旋偏角、俯仰角、横滚角、对比度、饱和度等)要求较高,这就使得作为作业平台的无人机载机系统具有抗风性强、飞行姿态稳定、巡航航线准确、操作简便等优良的性能。 
就现今社会无人机航测技术正在如火如荼的发展,其复杂的操控技术、高昂的成本与培训费用成为许多单位的绊脚石。若降低成本,设计出的无人机常常达不到规范要求,旋偏角过大、负载小,导致拍摄到的影响模糊不清晰。 
发明内容
为了解决现有技术的不足,本发明提供了一种固定翼航测型无人机设计方法,能够以最低成本费用制作出固定翼航测型无人机,并搭载航测设备获取优质的高分辨率数码影像数据。 
本发明为解决其技术问题所采用的技术方案是:提供了一种固定翼航测型无人机设计方法,具体包括以下步骤: 
(1)设计机翼参数; 
(1a)所述机翼采用S型翼型;根据以下公式确定雷诺系数Re: 
Re=ρVb/μ……(1) 
其中,ρ表示空气密度,V表示飞行速度,b表示机翼弦长,μ表示黏性系数;本定翼航测型无人机的飞行高度在1000米以下,空气密度、飞行速度和黏性系数均为已知量,机翼弦长b为设置值; 
(1b)根据步骤(1a)得到的雷诺系数Re确定以零升力迎角α为横轴、以机翼升力系数为纵轴的升力系数曲线Cl,所述升力系数曲线Cl用以下公式表示: 
C 1 = C 1 α α 1 + 18.2 C L αA = B 0 α 1 + 18.2 B 0 A . . . . . . ( 2 )
其中,B0为升力系数曲线的斜率,A为展弦比,通过以下公式计算: 
A=L2/S……(3) 
其中,L表示翼展,为设置值,S为主翼面积,为设置值; 
根据公式(2)得到当零升力迎角α为9°时,升力系数曲线Cl的斜率B0; 
(1c)根据步骤(1a)得到的雷诺系数Re确定以零升力迎角α为横轴、以翼型阻力系数为纵轴的翼型阻力系数曲线Cd0,再根据以下公式得到以零升力迎角α为横轴、以机翼阻力系数为纵轴的机翼阻力系数曲线Cd: 
Cd=Cd0+Cdi……(4) 
其中,Cdi表示诱导阻力系数,根据以下公式计算: 
C di = C 1 2 πA . . . . . . ( 5 )
则根据机翼阻力系数曲线得到当零升力迎角为9°时的机翼阻力系数Cd; 
(1d)机翼升力对于机翼上一点所产生的力矩不随零升力迎角改变而改变,以这一点作为支点,升力产生的力矩为常数,这一点为机翼焦点,升力对机翼焦点产生的力矩称为焦点力矩M0,通过以下公式计算: 
M 0 = 1 2 ρ V 2 SM z 0 . . . . . . ( 6 )
其中,Mz0为当零升力迎角α为9°时的焦点力矩系数; 
根据飞行速度V、翼弦长度L以及机翼面积S分别计算升力对重心产生的力矩M1、阻力对重心产生的力矩M2; 
根据以下公式计算重心距焦点的前后距离x: 
M1=1/2ρV2SClx……(7) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
根据以下公式计算重心距焦点的上下距离y: 
M2=1/2ρV2SCdy……(8) 
其中,Cl为零升力迎角α为9°时的机翼阻力系数; 
根据以下公式计算机翼对重心产生的力矩M: 
M=M1+M2+M0……(9) 
(1e)计算其他机翼设计参数: 
通过以下公式计算有力迎角α1: 
α1=MAX(Cl/Cd)……(10) 
通过以下公式计算升阻比R: 
R=Cl/Cd……(11) 
其中,Cl为零升力迎角α为9°时的机翼阻力系数,Cd为零升力迎角α为9°时的阻力系数; 
通过以下公式计算下洗角k: 
k=36.5Cl/A……(12) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算诱导迎角n: 
n=18.2Cl/A……(13) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算升力系数曲线的斜率B0: 
B0=Cl/(a+n)……(14) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算升力Y: 
Y=0.5vvpsCl……(15) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算理想阻力D1: 
D1=0.5vvpsCd……(16) 
其中,Cd为零升力迎角α为9°时的阻力系数; 
通过以下公式计算涡流诱导阻力D2: 
D2=0.5vvpsCdi……(17) 
其中,Cdi为零升力迎角α为9°时的诱导阻力系数; 
通过以下公式计算阻力D: 
D=D1+D2……(18) 
(2)机翼翼型建模: 
(2a)利用机翼翼型参数在CAD中进行制图,得到机翼翼型CAD数据; 
(2b)将机翼翼型CAD数据导入SketchUp进行3D建模,将翼展设置为2m,得到机翼3D模型; 
(3)设计机身参数; 
所述机身包括连接于机翼中部的机身前段、通过机身中段和机身尾段;根据机身前段的载重要求分别设计机身前段、机身中段和机身尾段的长度、宽度和高度; 
(4)机身建模: 
(4a)利用机身参数在CAD中进行制图,得到机身CAD数据; 
(4b)将机身CAD数据导入SketchUp进行3D建模,得到机身3D模型; 
(5)将机翼3D模型和机身3D模型打印为加工图纸,根据图纸进行加工制作和组装;机身采用epp与eps材质; 
(6)将搭载的电子设备安装在机身前段中; 
(7)对无人机进行飞行调试。 
进一步的改进是,在步骤(2)进行机翼翼型建模之前,将步骤(1)计算得到的机翼参数输入Profili软件进行模拟,得到不同迎角下的机翼受到的升力和阻力曲线分布图,验证设计参数是否正确;所述机翼参数包括零升力迎角α、飞行速度V、机翼弦长b、雷诺系数Re、翼弦长度L、前缘半径和最大弦厚位置。 
进一步的改进是,利用Profili软件进行模拟时,分别计算零升力迎角α为-3°、0°、3°、6°、9°和13°时,机翼受到的阻力和升力曲线分布图。 
进一步的改进是,该无人机的翼弦长度b=0.2m,翼展L=2m;计算雷诺系数Re时的空气密度p=1.226Kg/m3,飞行速度V=10m/s;计算机翼对重心产生的力矩M时的焦点力矩M0=0.4774N/m,升力产生的力矩M1=0.58848N/m,阻力产生的力矩M2=-0.0175N/m。 
进一步的改进是,步骤(6)所述的搭载的电子设备包括舵机、舵角、拉杆、接收机、遥控器、飞控、摄像头、Li-Po电池、GPS模块、图像传输发射机、图像传输接收机、数传电台发射机、数传电台接收机。 
本发明基于其技术方案所具有的有益效果在于: 
(1)本发明的设计方法根据空气动力学原理、模型空间结构力学原理、无线电通讯原理等理论知识,通过室内实验数据和多次外场试飞分析得到,采用S型翼型,并进一步确定展弦比、雷诺系数以及零升力迎角,翼型的此种组合避免了机翼交互在层流层与湍流层中过度,最大程度上保证了机翼在恶劣飞行环境下为飞机提供最大升力; 
(2)机身采用“卉”字型布局,机身前段的设置能够搭载所需的动力设备、测量设备等电子设备,2根中梁连接机身前段、机翼与机身尾端,保证了机身连接的牢固性;机身尾端包括两块尾梁和水平尾翼,能够控制飞行过程中飞机的飞行姿态;采用上述机身类型布局方式,相对于常规布局无人机在机身牢固性与稳定性方面大有提高; 
(3)本发明的整个机身材料在制作时可以采用高密度发泡的EPS材质,最大程度上保证了飞机高质量的同时使飞机自身重量降到最小; 
(4)利用本发明制成的无人机具有飞行风险低、成本费用低廉、操作灵活度高等优点,可以完全应用于开设与航测技术相关专业的高校及其学生的日常无人机教学课程和小面积范围航测作业。 
附图说明
图1是雷诺系数Re=137000时的升力系数曲线。 
图2是雷诺系数Re=137000时的机翼阻力系数曲线。 
图3是当零升力迎角α为-3°时,模拟得到的升力与阻力分布曲线图。 
图4是当零升力迎角α为0°时,模拟得到的升力与阻力分布曲线图。 
图5是当零升力迎角α为3°时,模拟得到的升力与阻力分布曲线图。 
图6是当零升力迎角α为6°时,模拟得到的升力与阻力分布曲线图。 
图7是当零升力迎角α为9°时,模拟得到的升力与阻力分布曲线图。 
图8是当零升力迎角α为13°时,模拟得到的升力与阻力分布曲线图。 
图9为机身结构示意图。 
图中:1-机身前段,2-机翼,3-中梁,4-尾梁,5-水平尾翼。 
具体实施方式
下面结合附图和实施例对本发明作进一步说明。 
本发明提供了一种固定翼航测型无人机设计方法,具体包括以下步骤: 
(1)设计机翼参数; 
(1a)所述机翼采用S型翼型; 
雷诺系数是空气密度、气流速度、翼弦长的乘积与空气黏性系数的比值,雷诺系数是判断机翼处于什么飞行状态的重要因素,雷诺系数越大流经机翼表面的边界层越早从层流边层过渡为湍流边界层,而湍流边界层不容易从翼表面分离,所以比较不容易失速。失速是指机翼提供的升力小于飞机的总阻力,导致飞机将要坠落的状态。雷诺系数小的机翼边界层尚未从层流边层过渡为湍流边界层时就先分离了,很容易造成机翼失速,所 以在设计该翼型时首要的要知道其最适合的雷诺系数。 
根据以下公式确定雷诺系数Re: 
Re=ρVb/μ……(1) 
其中,ρ表示空气密度,V表示飞行速度,b表示机翼弦长,μ表示黏性系数;本定翼航测型无人机的飞行高度在1000米以下,空气密度、飞行速度和黏性系数均为已知量,空气密度ρ=1.226Kg/m3,ρ/μ为定值68500,飞行速度V=10m/s,机翼弦长b设为0.2m;计算得出雷诺系数Re=137000。 
(1b)根据步骤(1a)得到的雷诺系数Re确定以零升力迎角α为横轴、以机翼升力系数为纵轴的升力系数曲线Cl,该曲线用以下公式表示: 
C 1 = C 1 α α 1 + 18.2 C L αA = B 0 α 1 + 18.2 B 0 A . . . . . . ( 2 )
其中,B0为升力系数曲线的斜率,A为展弦比,通过以下公式计算: 
A=L2/S……(3) 
其中,L表示翼展,S为主翼面积。由公式(3)可知,展弦比应该愈大愈好,但大展弦比机翼是很难制作得又轻又坚固的。对于无人机来说,考虑展弦比的时候还应该同时考虑到雷诺系数的影响。由于主翼的面积有一定的限制,所以用大展弦比就要求短翼弦,也就是小雷诺系数,小雷诺系数时机翼容易失速。雷诺系数愈大,机翼的性能便愈好,尤其是最大升力系数受雷诺系数的影响更大,所以在进行翼型设计时要综合考虑到二者的影响。L主要是参考了飞机的负载,因为L越大飞机的负载能力越大,但由此产生的阻力也会越大,基于本固定翼航测型无人机,其L取值为2m时,在制作出来后,配上任务荷载2.2kg进行测试飞行,完全可以符合荷载要求,飞机自重加任务荷载总共重4.6kg。综合考虑,翼展L设置为2m,主翼面积设置为0.4m2。 
如图1所示为雷诺系数Re=137000时的升力系数曲线,由该曲线可知,该翼型升力系数随着零升力迎角α的增大逐渐增大,在零升力迎角α达到9°时升力系数逐渐降低,这也就意味着机翼的升力逐渐开始减少,之后如果在大迎角飞行时将会损失大部分的升力;确定零升力迎角α后,根据公式(2)得到当零升力迎角α为9°时升力系数曲线Cl的斜率B0,B0=0.104。 
(1c)根据步骤(1a)得到的雷诺系数Re确定以零升力迎角α为横轴、以翼型阻 力系数为纵轴的翼型阻力系数曲线Cd0,再根据以下公式得到以零升力迎角α为横轴、以机翼阻力系数为纵轴的机翼阻力系数曲线Cd: 
Cd=Cd0+Cdi……(4) 
其中,Cdi表示诱导阻力系数,根据以下公式计算: 
C di = C 1 2 πA . . . . . . ( 5 )
所得的机翼阻力系数曲线如图2所示,由该曲线可知当机翼迎角低于零升力迎角时,阻力系数随着迎角的逐渐增大而减小,但当过了零升力迎角阻力系数随着迎角的增大逐渐增大,在迎角达到9°增大的斜率逐渐增大,如果飞机飞行迎角超过这个角度,那将会处于阻力大于升力的状态,飞机将处于失速状态会产生危险;则根据机翼阻力系数曲线得到当零升力迎角为9°时的机翼阻力系数Cd; 
(1d)机翼升力对于机翼上离前缘约1/4机翼弦长处的一点所产生的力矩不随零升力迎角改变而改变,以这一点作为支点,升力产生的力矩为常数,这一点为机翼焦点,升力对机翼焦点产生的力矩称为焦点力矩M0,通过以下公式计算: 
M 0 = 1 2 ρ V 2 SM z 0 . . . . . . ( 6 )
其中,Mz0为焦点力矩系数,当零升力迎角α为9°时焦点力矩系数Mz0为-0.13,负号表示力矩具有使模型飞机低头的趋势;计算得到焦点力矩M0=0.4774N/m。 
对重心产生的力矩一共有3个,分别为焦点力矩M0、升力对重心产生的力矩M1、阻力对重心产生的力矩M2。根据飞行速度V、翼弦长度L以及机翼面积S能够分别计算升力对重心产生的力矩M1=0.58848N/m、阻力对重心产生的力矩M2=-0.0175N/m。 
根据以下公式计算重心距焦点的前后距离x: 
M1=1/2ρV2SClx……(7) 
其中,Cl为零升力迎角α为9°时的机翼升力系数,为0.8;计算得到重心距焦点的前后距离x=0.03m。 
根据以下公式计算重心距焦点的上下距离y: 
M2=1/2ρV2SCdy……(8) 
其中,Cd为零升力迎角α为9°时的机翼阻力系数,为0.05;计算得到重心距焦点的上下距离y=0.01m。 
x与y的取值本身为不确定量,因为决定其数值的焦点的位置也是不确定的,在飞机飞行过程中,重心和焦点的位置在天空中是时刻在变的,具体值根据已有的机翼参数进行参考,然后在此基础上进行设计。x与y的取值是参考了T2机翼的参数,但是又在其上进行了改正,确定出重心和焦点位置,根据最大弦厚位置、下弧线、上弧线的位置计算出x,y的值,x,y值的表示在翼型设计图中可以看出相对位置关系。 
根据以下公式计算机翼对重心产生的力矩M: 
M=M1+M2+M0……(9) 
(1e)计算其他机翼设计参数: 
通过以下公式计算有力迎角α1: 
α1=MAX(Cl/Cd)……(10) 
通过以下公式计算升阻比R: 
R=Cl/Cd……(11) 
其中,Cl为零升力迎角α为9°时的机翼阻力系数,Cd为零升力迎角α为9°时的阻力系数; 
通过以下公式计算下洗角k: 
k=36.5Cl/A……(12) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算诱导迎角n: 
n=18.2Cl/A……(13) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算升力系数曲线的斜率B0: 
B0=Cl/(a+n)……(14) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算升力Y: 
Y=0.5vvpsCl……(15) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算理想阻力D1: 
D1=0.5vvpsCd……(16) 
其中,Cd为零升力迎角α为9°时的阻力系数; 
通过以下公式计算涡流诱导阻力D2: 
D2=0.5vvpsCdi……(17) 
其中,Cdi为零升力迎角α为9°时的诱导阻力系数; 
通过以下公式计算阻力D: 
D=D1+D2……(18) 
至此,得到如下表所示的机翼翼型参数: 
表1机翼翼型参数 
在步骤(2)进行机翼翼型建模之前,可以将步骤(1)计算得到的机翼参数输入Profili软件进行模拟,得到不同迎角下的机翼受到的升力和阻力曲线分布图,验证设计参数是否正确;所述机翼参数包括零升力迎角α、飞行速度V、机翼弦长b、雷诺系数Re、翼弦长度L、前缘半径和最大弦厚位置。 
利用Profili软件进行模拟时,分别计算零升力迎角α为-3°、0°、3°、6°、9°和13°时,机翼受到的阻力和升力曲线分布图。 
当零升力迎角α为-3°时,模拟得到的升力与阻力分布曲线图如图3所示。可以看出在零升力迎角α为负时,机翼产生向下的升力,此角度不利于飞机飞行。 
当零升力迎角α为0°时,模拟得到的升力与阻力分布曲线图如图4所示。在零升力迎角α逐渐增大的过程中机翼产生向上的升力并且逐渐增大,但还是有部分向下的负 升力。 
当零升力迎角α为3°时,模拟得到的升力与阻力分布曲线图如图5所示。此零升力迎角α下,升力逐渐增大,但同时阻力也随之增大。 
当零升力迎角α为6°时,模拟得到的升力与阻力分布曲线图如图6所示。在此零升力迎角α下,升力与阻力同时逐渐增大。 
当零升力迎角α为9°时,模拟得到的升力与阻力分布曲线图如图7所示。在此零升力迎角α下升力与阻力随着零升力迎角α的增大都增大,但是机翼上表面已经出现了升力不稳定的趋势,翼尖与翼根压差过大,气流将在这个阶段发生层流层与湍流层的过度,飞机将进入升力与阻力混乱的现象,也就是失速状态。 
当零升力迎角α为13°时,模拟得到的升力与阻力分布曲线图如图8所示。在此阶段流进机翼表面的气流已经相当紊乱,机翼侧端在这种情况下很容易使下表面的气流往上表面流动,形成翼尖涡流,从而使诱导阻力增大,飞机不应在此零升力迎角α下飞行。 
由模拟结果可以进一步验证设计机翼参数时的零升力迎角α取值应为9°。 
(2)机翼翼型建模: 
(2a)利用机翼翼型参数在CAD中进行制图,得到机翼翼型CAD数据; 
(2b)将机翼翼型CAD数据导入SketchUp进行3D建模,将翼展设置为2m,得到机翼3D模型; 
(3)设计机身参数;如图9所示: 
(3a)机身前段1:载重舱垂直于机翼2,其尾部与机翼中心处连接,考虑到搭载任务设备以及自身电子设备的能力,要保障在容纳下最多任务设备的同时使机身最小,将其长度设置为翼展的1/2; 
(3b)垂直于机翼的2根中梁3:分别位于机身前段的左边和右边,2根中梁的前部分别与机翼连接,其长度设置为翼展的2/5到3/5; 
(3c)机身后段:机身后段包括分别位于2根中梁尾部上方的尾梁4,以及安装在尾梁上方、平行于机翼的水平尾翼5;考虑到其与机翼及机身前段的链接紧密程度和质量强度,在此基础上设计两根尾梁也是出于飞行稳定性的考虑,一般设计而言机翼与机身的设计尺寸比例是1:0.6—1:1,本无人机的水平尾翼的长度设置为翼展的1/2; 
最后适当选择所搭载设备的尺寸,得到如下表所示的机身设计参数: 
表2机身设计参数 
(4)机身建模: 
(4a)利用机身参数在CAD中进行制图,得到机身CAD数据; 
(4b)将机身CAD数据导入SketchUp进行3D建模,得到机身3D模型; 
(5)将机翼3D模型和机身3D模型打印为加工图纸,根据图纸进行加工制作和组装;机身采用epp与eps材质; 
(6)将搭载的电子设备安装在机身前段中;所述的搭载的电子设备包括舵机、舵角、拉杆、接收机、遥控器、飞控、摄像头、Li-Po电池、GPS模块、图像传输发射机、图像传输接收机、数传电台发射机、数传电台接收机; 
(7)对无人机进行飞行调试。 

Claims (7)

1.一种固定翼航测型无人机设计方法,其特征在于具体包括以下步骤: 
(1)确定机翼翼型参数; 
(1a)所述机翼采用S型翼型;根据以下公式确定雷诺系数Re: 
Re=ρVb/μ……(1) 
其中,ρ表示空气密度,V表示飞行速度,b表示机翼弦长,μ表示黏性系数;本定翼航测型无人机的飞行高度在1000米以下,空气密度、飞行速度和黏性系数均为已知量,机翼弦长b为设置值; 
(1b)根据步骤(1a)得到的雷诺系数Re确定以零升力迎角α为横轴、以机翼升力系数为纵轴的升力系数曲线Cl,所述升力系数曲线Cl用以下公式表示: 
其中,B0为升力系数曲线的斜率,A为展弦比,通过以下公式计算: 
A=L2/S……(3) 
其中,L表示翼展,为设置值,S为主翼面积,为设置值; 
根据公式(2)得到当零升力迎角α为9°时,升力系数曲线Cl的斜率B0; 
(1c)根据步骤(1a)得到的雷诺系数Re确定以零升力迎角α为横轴、以翼型阻力系数为纵轴的翼型阻力系数曲线Cd0,再根据以下公式得到以零升力迎角α为横轴、以机翼阻力系数为纵轴的机翼阻力系数曲线Cd: 
Cd=Cd0+Cdi……(4) 
其中,Cdi表示诱导阻力系数,根据以下公式计算: 
则根据机翼阻力系数曲线得到当零升力迎角为9°时的机翼阻力系数Cd; 
(1d)机翼升力对于机翼上一点所产生的力矩不随零升力迎角改变而改变,以这一点作为支点,升力产生的力矩为常数,这一点为机翼焦点,升力对机翼焦点产生的力矩称为焦点力矩M0,通过以下公式计算: 
其中,Mz0为当零升力迎角α为9°时的焦点力矩系数; 
根据飞行速度V、翼弦长度L以及机翼面积S分别计算、阻力对重心产生的力矩M2; 
根据以下公式计算升力对重心产生的力矩M1: 
M1=1/2ρV2SClx……(7) 
其中,Cl为零升力迎角α为9°时的机翼升力系数,x为重心距焦点的前后距离; 
根据以下公式计算: 
M2=1/2ρV2SCdy……(8) 
其中,Cl为零升力迎角α为9°时的机翼阻力系数,y为重心距焦点的上下距离; 
根据以下公式计算机翼对重心产生的力矩M: 
M=M1+M2+M0……(9) 
(1e)确定其他机翼设计参数: 
通过以下公式计算有力迎角α1: 
α1=MAX(Cl/Cd)……(10) 
其中,Cl和Cd分别为机翼升力系数曲线和阻力系数曲线中的点的纵坐标,他们的取值使得Cl/Cd为最大值; 
通过以下公式计算升阻比R: 
R=Cl/Cd……(11) 
其中,Cl为零升力迎角α为9°时的机翼阻力系数,Cd为零升力迎角α为9°时的阻力系数; 
通过以下公式计算下洗角k: 
k=36.5Cl/A……(12) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算诱导迎角n: 
n=18.2Cl/A……(13) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算升力Y: 
Y=0.5V2ρSCl……(14) 
其中,Cl为零升力迎角α为9°时的机翼升力系数; 
通过以下公式计算理想阻力D1: 
D1=0.5V2ρSCd……(15) 
其中,Cd为零升力迎角α为9°时的阻力系数; 
通过以下公式计算涡流诱导阻力D2: 
D2=0.5V2ρSCdi……(16) 
其中,Cdi为零升力迎角α为9°时的诱导阻力系数; 
通过以下公式计算阻力D: 
D=D1+D2……(17) 
将前缘半径设置为翼弦长度的1/10; 
将最大弦厚位置设置为翼弦长度的1/3; 
(2)机翼翼型建模: 
(2a)利用步骤(1)得到的机翼翼型参数在CAD中进行制图,得到机翼翼型CAD数据; 
(2b)将机翼翼型CAD数据导入SketchUp进行3D建模,将翼展设置为2m,得到机翼3D模型; 
(3)确定机身参数; 
(3a)机身前段:载重舱垂直于机翼,其尾部与机翼中心处连接,其长度设置为翼展的2/5到3/5; 
(3b)垂直于机翼的2根中梁:分别位于机身前段的左边和右边,2根中梁的前部分别与机翼连接,其长度设置为翼展的2/5到3/5; 
(3c)机身后段:机身后段包括分别位于2根中梁尾部上方的尾梁,以及安装在尾梁上方、平行于机翼的水平尾翼,水平尾翼的长度设置为翼展的1/5到2/5; 
(4)机身建模: 
(4a)利用机身参数在CAD中进行制图,得到机身CAD数据; 
(4b)将机身CAD数据导入SketchUp进行3D建模,得到机身3D模型; 
(5)将机翼3D模型和机身3D模型打印为加工图纸,根据图纸进行加工制作和组 装; 
(6)将搭载的电子设备安装在机身前段中; 
(7)对无人机进行飞行调试。 
2.根据权利要求1所述的固定翼航测型无人机设计方法,其特征在于:在步骤(2)进行机翼翼型建模之前,将步骤(1)计算得到的机翼参数输入Profili软件进行模拟,得到不同迎角下的机翼受到的升力和阻力曲线分布图,验证设计参数是否正确;所述机翼参数包括零升力迎角α、飞行速度V、机翼弦长b、雷诺系数Re、翼弦长度L、前缘半径和最大弦厚位置。 
3.根据权利要求2所述的固定翼航测型无人机设计方法,其特征在于:利用Profili软件进行模拟时,分别计算零升力迎角α为-3°、0°、3°、6°、9°和13°时,机翼受到的阻力和升力曲线分布图。 
4.根据权利要求1所述的固定翼航测型无人机设计方法,其特征在于:该无人机的翼弦长度b=0.2m,翼展L=2m;计算雷诺系数Re时的空气密度p=1.226Kg/m3,飞行速度V=10m/s;计算机翼对重心产生的力矩M时的焦点力矩M0=0.4774N/m,升力产生的力矩M1=0.58848N/m,阻力产生的力矩M2=-0.0175N/m。 
5.根据权利要求1所述的固定翼航测型无人机设计方法,其特征在于:该无人机的机身前段长度设置为翼展的1/2,中梁的长度设置为翼展的1/5,水平尾翼的长度设置为翼展的1/2。 
6.根据权利要求1所述的固定翼航测型无人机设计方法,其特征在于:机身采用epp与eps材质。 
7.根据权利要求1所述的固定翼航测型无人机设计方法,其特征在于:步骤(6)所述的搭载的电子设备包括舵机、舵角、拉杆、接收机、遥控器、飞控、摄像头、Li-Po电池、GPS模块、图像传输发射机、图像传输接收机、数传电台发射机、数传电台接收机。 
CN201410665303.5A 2014-11-19 2014-11-19 一种固定翼航测型无人机设计方法 Expired - Fee Related CN104401504B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410665303.5A CN104401504B (zh) 2014-11-19 2014-11-19 一种固定翼航测型无人机设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410665303.5A CN104401504B (zh) 2014-11-19 2014-11-19 一种固定翼航测型无人机设计方法

Publications (2)

Publication Number Publication Date
CN104401504A true CN104401504A (zh) 2015-03-11
CN104401504B CN104401504B (zh) 2016-01-06

Family

ID=52639197

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410665303.5A Expired - Fee Related CN104401504B (zh) 2014-11-19 2014-11-19 一种固定翼航测型无人机设计方法

Country Status (1)

Country Link
CN (1) CN104401504B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106886625A (zh) * 2017-01-05 2017-06-23 北京航天自动控制研究所 一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法
CN107066636A (zh) * 2016-05-24 2017-08-18 成都前沿动力科技有限公司 一种多旋翼无人机増升增稳器
CN115783297A (zh) * 2022-12-09 2023-03-14 太原理工大学 一种人力低空飞行自行车设计方法及结构
CN118260873A (zh) * 2024-05-27 2024-06-28 北京舯迦科技有限公司 一种减小全机飞行阻力的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101458735A (zh) * 2008-12-31 2009-06-17 重庆大学 一种具有高升阻比的翼型
US20100123047A1 (en) * 2008-11-14 2010-05-20 Williams Aerospace, Inc. Blended Wing Body Unmanned Aerial Vehicle
CN101898644A (zh) * 2009-05-25 2010-12-01 丹麦技术大学 用来设计翼型的系统与方法
US20110202321A1 (en) * 2008-11-24 2011-08-18 Rolls-Royce Plc method for optimising the shape of an aerofoil

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100123047A1 (en) * 2008-11-14 2010-05-20 Williams Aerospace, Inc. Blended Wing Body Unmanned Aerial Vehicle
US20110202321A1 (en) * 2008-11-24 2011-08-18 Rolls-Royce Plc method for optimising the shape of an aerofoil
CN101458735A (zh) * 2008-12-31 2009-06-17 重庆大学 一种具有高升阻比的翼型
CN101898644A (zh) * 2009-05-25 2010-12-01 丹麦技术大学 用来设计翼型的系统与方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107066636A (zh) * 2016-05-24 2017-08-18 成都前沿动力科技有限公司 一种多旋翼无人机増升增稳器
CN106886625A (zh) * 2017-01-05 2017-06-23 北京航天自动控制研究所 一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法
CN106886625B (zh) * 2017-01-05 2020-04-14 北京航天自动控制研究所 一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法
CN115783297A (zh) * 2022-12-09 2023-03-14 太原理工大学 一种人力低空飞行自行车设计方法及结构
CN118260873A (zh) * 2024-05-27 2024-06-28 北京舯迦科技有限公司 一种减小全机飞行阻力的方法
CN118260873B (zh) * 2024-05-27 2024-08-02 北京舯迦科技有限公司 一种减小全机飞行阻力的方法

Also Published As

Publication number Publication date
CN104401504B (zh) 2016-01-06

Similar Documents

Publication Publication Date Title
Kontogiannis et al. Design, performance evaluation and optimization of a UAV
Selig Modeling full-envelope aerodynamics of small UAVs in realtime
Patel et al. CFD Analysis of an Aerofoil
Diehl Engineering aerodynamics
Selig Real-time flight simulation of highly maneuverable unmanned aerial vehicles
CN100541372C (zh) 一种无人机发动机意外停车下的自动归航控制方法
CN104401504B (zh) 一种固定翼航测型无人机设计方法
Saban et al. Simulation of wake vortex effects for UAVs in close formation flight
Stewart et al. Design of the air force research laboratory micro aerial vehicle research configuration
Stone The T-wing tail-sitter unmanned air vehicle: from design concept to research flight vehicle
Joels et al. Design, Analysis, and Testing of the Active Aeroelastic Aircraft Testbed (A3TB) Platform
Joels et al. Design, Analyses, and Flutter Testing of the Active Aeroelastic Aircraft Testbed (A3TB) Platform
Harasani et al. Initial conceptual design and wing aerodynamic analysis of a solar power-based UAV
Brinkworth On the aerodynamics of the Miles Libellula tandem-wing aircraft concept, 1941–1947
Koçak Effect of ground on flow structure of non-slender delta and reverse delta wings
Abdulrahim et al. Development of Mission capable Flexible-Wing Micro Air Vehicles
Pettersson Design of a drone system for maritime search and rescue missions
Koster Hyperion-three years of novel aircraft design
Moidel Development and Validation of Single-Engine General Aviation Aircraft Models within Merlin 521 Motion-Based Flight Simulators
Standridge et al. Design, Development, and Range Optimization of Flying Wing UAV
de Almeida Induced Drag Reduction with a Wing Tip Tail Configuration
Stefan et al. Skyhunter UAS Dynamic Model and Optimal Static Margin Analysis
Concepción Morales CFD-Based Method for Validation of DATCOMand Potential Flow Stability and Control Derivatives
Kanistras A comprehensive methodology for design of a circulation control small-scale unmanned aircraft
Morgan A new shape in the sky

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Chen Gang

Inventor after: Xu Wei

Inventor after: Wu Peng

Inventor after: Zhao Ruyue

Inventor before: Xu Wei

Inventor before: Chen Gang

Inventor before: Wu Peng

Inventor before: Zhao Ruyue

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: XU WEI CHEN GANG WU PENG ZHAO RUYUE TO: CHEN GANG XU WEI WU PENG ZHAO RUYUE

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160106

Termination date: 20171119

CF01 Termination of patent right due to non-payment of annual fee