无人飞机总体设计---设计过程及算例

无人机总体设计算例任务要求飞行高度30200M，飞行速度4090KM/H，巡航速度18M/S，最大飞行速度28M/S，爬升率4M/S，续航时间1H，最大过载17，任务载荷重量05KG，背包式运输，发射方式手抛式，回收方式机腹着陆设计过程1布局形式及布局初步设计无尾布局【方法参考已有同类无人机】确定布局形式主要是机翼、垂尾、动力、起落架等。（1）机翼根据经验或同类飞机确定展弦比556，尖削比0405，后掠角28，下反角15，安装角2展弦比【展弦比增大，升致阻力减小，升阻比增大】【展弦比增大，弦长减小，雷诺数降低，气动效率降低】【展弦比增大，弦长减小，翼型厚度减小，机翼结构重量上升】尖削比【尖削比影响升力展向分布，当展向升力分布接近椭圆时，升致阻力最小，低速机翼一般取0405】后掠角【后掠角增加，横向稳定性增大，配下反角】【后掠角增加，尾翼舵效增加】【后掠角增加，纵向阻尼增强，纵向动稳定性增强】下反角【上反角增加，横向稳定性增加，下反角相反】安装角【巡航阻力最小对应机翼的迎角，通用航空飞机和自制飞机的安装角大约为2，运输机大约为1，军用飞机大约为0，在以后的设计阶段，可通过气动计算来检查设计状态所需要的机翼实际的安装角。】机翼外型草图（2）垂尾垂尾形式翼尖垂尾尾空系数CVT004/2002【双重尾】（3）动力系统形式电动无人机推进系统安装位置主要有机头拉进式、机尾推进式、单发机翼前缘拉进式、双发形式、单发机翼后缘推进式。下面研究各种布置形式对布局设计的影响。动力形式优点缺点实例机头拉进式螺旋桨前方进气稳定未被干扰；容易实现重心位置设计；手抛发射不会对发射员造成危害；排气被机身和机翼阻止，影响动力系统的效率；回收降落时，电动机和螺旋桨容易触地损坏机尾推进式机头可以安装任务设备；螺旋桨也不容易在着陆时触地损坏；对螺旋桨的干扰较小；重心配置在设计重心点非常困难；单发翼前缘拉进式电动机不在占用机头位置；以便在机头安装任务设备；机身的阻力会产生一个较大的低头力矩；过高的机身也增大的结构重量，浸润面积也比较大双发翼前缘拉进式机头安装摄像设备布置需要两台电动机，增加了系统的复杂性单发机翼后缘推进式机头安装摄像设备螺旋桨的滑流直接吹在尾翼上，造成无人机的稳定性变化本方案为机尾推进式2无人机升阻特性（极曲线）估算前面确定了机翼的基本参数，要确定无人机的具体机翼参数，还需要知道“起飞重量”、“翼载荷”，然后进行布局缩放。确定起飞重量，关键是电池重量，电池重量由飞机需要的能量决定，能量由飞机升阻特性决定。升阻特性由飞机布局形式决定，可参考同类飞机，进行初步估算。飞机的极曲线20,0DDILCCK（1）零升阻力系数，一般可取为2X（一张纸打比方）DFESC浸湿参考【参考面积统一为机翼面积】对于机身34/2SS浸湿侧俯对于机翼、尾翼，一般以翼型最大相对厚度为基础计算。/050319752/TCSTCS浸湿外露浸湿外露也可以直接根据各类飞机的统计值，选取参考值。这里假设机翼，则/01TC209S浸湿机身取，则5S侧俯341072浸湿垂尾，则01外露09S浸湿052972132DFESC浸湿参考（2）升致阻力因子1KAE对于后掠翼飞机0680154615COS32718LEEA173583058E至此，可以估算得到飞机的极曲线210DLCC0020406081121416000100200300400500600700800901011012013014015016017018CDCL（3）飞机极曲线20137DLCC升阻比最大时，2/45064LD最大升阻比MAX/16LD3功重比与翼载荷的确定如果飞机重量知道，获得了升阻特性，根据速度可以得到功率需求，根据航时要求可以得到能量要求，即起飞重量决定功率能量但是起飞重量主要包括机体结构、任务设备、动力装置、电池。而电池重量又决定它包含的能量的多少。即功率能量决定起飞重量确定其中一个需要依靠对方，从而提出功重比的概念。起飞重量决定机翼大小，机翼大小又决定起飞重量，从而提出翼载荷的概念。根据功率需求，可推出飞机功重比与翼载荷的约束分析方程220GSWQNKCSQVGWPTDTYT一般情况下，可先根据经验值确定翼载，然后在无人机巡航、爬升、盘旋、最大飞行速度等多个工况下，由翼载计算功重比。表41无尾布局小型电动无人机参数统计名称翼展M机长M机翼面积M2重量KG翼载荷KG/M2DRAGONEYE114090352777DUIGAN3096572P150351510605252946552876UAVZALA42108080410251768从统计值可知，翼载可取7KG/M2代入上式，可得到巡航状态V18M/S功重比为1119W/KG爬升状态手抛速度V10M/S22MAX11134LLCWVS起飞V05（人手抛速度巡航速度）12M/S，VY4M/S功重比为484W/KG巡航盘旋状态V18M/S，N173功重比为201W/KG最大平飞速度状态V28M/S功重比为339W/KG工况功重比巡航状态1119W/KG爬升状态484W/KG巡航盘旋状态201W/KG最大平飞速度状态339W/KG由上得出最大功重比为484W/KG，巡航功重比为1119W/KG0501001502002500123456789101112131415WT/SP/WT12345实际上，各种工况下，翼载与功重比之间关系图可以画出来，然后根据一些限制条件（起飞距离。），找范围，确定相应满足条件的翼载和功重比若干组。4起飞重量确定1234TW其中，是结构重量，是动力装置重量，是电池重量，是航1W24W空电子与任务设备。其中，在重量设计中是不变的，是任务要求4中给定的。（1）飞机结构重量1TWF其中，为结构重量系数。一般起飞重量在几公斤范围内的小型无1F人机结构重量系数在025035范围内，作为初步分析，可取为03。常规飞机种类结构重量系数飞机种类1F轻型030032中型028030亚音速干线客机重型025027。（2）动力装置重量动力装置包括电机、减速器、螺旋桨等。电动飞机起飞重量不随飞行发生变化。2TWF推导过程MAXMAX2/TTDJDJP其中，为电机的最大输出功率，为飞机最大功重比，MAXPAX/TW为动力装置的比功率（功率/动力装置重量）。这一参数可以取统DJ计值。【分析最大功重比为484W/KG，小型手抛电动无人机重量不大于5KG，因此，最大需求的功率250W】注通常手抛电动无人机300W的电机重量约为100G，电调约为50G，电机与螺旋桨连接器为30G。从而有，动力装置的重量约为205WKG（3）电池重量电池重量能量/能量密度3/WEE其中，为飞行中电池提供的能量，为电池实际比能量（能量密度）E。/4EPT其中，为飞行中电池提供的平均功率，为飞行时间。/4P由于飞机在爬升段需要较高功率，在飞行高度不高（相对地面最大平飞功率/最大效率。（如果多桨，则可以根据功率情况进行选择，以”可用功率需用功率且两者接近”为准则，排除一部分）【分析由前可知，最大需用功率为434W，可在最大飞行速度下，螺旋桨功率满足大于且接近的要求。最小功率需求是在12M/S下为13W，在12M/S下，螺桨最小转速功率为29W，较为接近。】第五步一旦选定螺旋桨，则根据巡航速度V下的效率最大化确定巡航最佳转速。（这就为电机选择提出了要求）A巡航状态昌敏做法巡航速度18M/S，推力09605N0203040506070809105055060650707508J功功D0810/ND12/ND14/ND16/ND18/ND20/ND22/ND24/ND26/ND28/NDJ02953035430413404724053150590606496070870767708268051059069071075075074072072069NV18M/S29997522500217214278418751631666654149987613636501249932115387110713921000017999150011285711251100000899908182075000692306428T232096156272114113765945644637164248941621710483058730608112141618X1040404505055060650707508功功功功0608112141618X1040510152025功功T巡航转速约6500转，效率约072从功率角度也可以，避免了以下的迭代再由算出功率，计算出扭矩。242535TQPCCNDNND为电机选择作输入，选取效率最高的电机。电机最大工作电压169V下，计算不同速度下的可用推力或功率，进而确定最大最小飞行速度，即速度范围。需要迭代计算，迭代出合适的转速。效率就不考虑了。巡航速度18M/S满足效率最高，螺旋桨功率17923405PROLJPW186/MIN05240VNSRJD功率校核（实际上不用校推力，只要功率即可）2234W35350621560248103PCNDW0995N24792TN以N9600R/MIN为巡航转速，效率最高，但螺旋桨功率过高，不匹配，不合适。【降低转速，损失一点效率，换取功率】取N7500R/MIN，J07087，效率为072，螺旋桨功率1792380PROLJPW功率校核2388W3535048912524086PCND取N6000R/MIN，J08858，效率为066螺旋桨功率72605PROLJPW功率校核2388W3535041257/6240876PCND反复迭代，约6900转为最佳转速，可以满足效率与功率兼得。巡航最佳转速N6800R/MIN，J07816，070螺旋桨功率1792450PROLJPW功率校核3535081256/02483PCND巡航扭矩40/PROQNM【分析如果想进一步提高效率，则需换桨，因此要准备尽可能多的螺桨作为备选桨。如果选择了效率最高的桨，仍想再提高效率，则需要改进飞机升阻特性。也就是说，一方面改桨，一方面改飞机升升阻特性】【总结为什么不能用需用功率、推力反推转速，因为这是一个隐式关系，无法事前确定CT,CP】B最大飞行速度状态飞机需用功率434W，飞行速度28M/S由前面的功率转速速度表可得出，取转速N1625R/S，9750R/MIN，前进比为J08480，效率065，螺旋桨功率436805PROLJPW功率校核611W353506125240819PCND取N167R/S（10000R/MIN），前进比为J08274，效率068螺旋桨功率43680PROLJPW功率校核1443W35350712240834PCND功率不接近。取转速N1875R/S（11250R/MIN），J03307，0542146/VMS螺旋桨功率7361054PROLJPW功率校核1443W353507182724089PCND取转速N175R/S（10500R/MIN），J03543，0592146/VMS螺旋桨功率736124059PROLJPW功率校核1247W3535071272408163PCND取转速N1625R/S（9750R/MIN），J03816，0642146/VMS螺旋桨功率7361504PROLJPW功率校核1150W353507162502481304PCND取转速N150R/S（9000R/MIN），J04134，0692146/VMS螺旋桨功率736109PROLJPW功率校核1067W3535071250248102PCND反复迭代，爬升最佳转速约为9000R/MIN爬升扭矩为106713229/PROPQNMN状态螺桨效率螺桨功率最佳转速扭矩爬升0691067900001132巡航070246680000340最大速0706201020000585（3）电机的选择电机的主要性能参数有，内阻，空载电流VKMR0I电机的效率PRODJPUI选择电机的要求巡航效率高，电机的输出扭矩0TQKI电动机的电压MVNUR电机扭矩常数与KV值的关系3095TVK根据上面公式QIU备选电动机的性能参数型号KV空载电流I0（A）内阻RM）重量（G）HACKERA2034S1500090147425HACKERA2022L924080109567HIMAXHC281206506500360285602HPZ300726124018006765转速单位R/MIN飞机巡航状态下电机的电流、电压、功率、效率型号电流电压功率效率HACKERA2034S626845454834192907194HACKERA2022L410697807032062607672HIMAXHC281206502686311227130159408157HPZ300726623795858136542206732飞机爬升状态下电机的电流、电压、功率、效率型号电流电压功率效率HACKERA2034S18773787597164452006488HACKERA2022L118102110276130238208193HIMAXHC2812065081053161562130950808148HPZ30072616575682526136791807800飞机最大飞行速度状态下电机的电流、电压、功率、效率型号电流电压功率效率HACKERA2034S1013688290184035107378HACKERA2022L6489911746476233008133HIMAXHC281206504362616935673883208392HPZ300726943588792082959607474【分析由上面可以看出，电机效率最高为HIMAXHC28120650，其工作电压最大，工作电流最小。】综上所述，本方案螺旋桨采用TAIPAN86，电机采用HIMAXHC28120650，巡航状态桨的效率070，电机效率08157，巡航状态电机电压113V，电机最大工作电压169V，采用5节聚合物锂电池串联，电压为375185，电池重量031KG。近似等于与之前估计值。【注如果严重大于前面估计值，还得重新走一遍前面的设计工作。】6飞机布局几何参数确定（1）机翼几何参数根据翼载可得机翼面积2/1536/70195TSWM翼展029BA几何平均弦长21B平均气动弦长205805823703BZDZDZMSS根弦长1961B尖弦长0133机翼视图（2）翼型的选择本方案设计的无尾布局电动无人机尺寸小，飞行速度低，雷诺数很小。翼型厚度相对较小【不能太小，重量过大】。起飞和着陆段可能需要人工遥控飞行，飞机必须具有很好的自然飞行稳定性，由于飞翼布局无平尾，这要求机翼具有正的零升俯仰力矩。总体对翼型的要求S型翼型，较高的升阻比，低雷诺数下的翼型效率较高，在整个飞行速度范围内力矩线性变化。现有的小型无尾式无人机和飞翼模型的翼型有EMX07、MH62、E186、S5010、HS510。备选翼型翼型相对厚度最大厚度位置相对弯度最大弯度位置EMX0799297253206MH62929269159366E18610232913029S50109822721927HS5108792721927A设计点设计升力系数【设计升力系数是指飞机常用的升力系数，通常指巡航飞行时的升力系数。】21179803457058LXHWGCSQ设计雷诺数【采用几何平均气动弦长S/B】5512819/RE21070VBB翼型气动性能分析从CMALPHA曲线上可以看出只有EMX07、E186零升俯仰力矩系数为正，其它的均为负，纵向配平较难。E186零升俯仰力矩系数大，但从升阻比曲线上可以看出，EMX07最大升阻比大，从CLALPHA曲线上看出，EMX07失速迎角大。从CLCD曲线上可以看出，在设计升力03457附近，阻力基本不变。而且在不同雷诺数下，EMX07的零升俯仰力矩系数变化不大。综上分析，本方案选用翼型为EMX07最大厚度零升迎角零升力矩系数最大升力系数最大升阻比失速迎角990800151217511翼型升力线斜率线性迎角范围6959（3）垂尾设计尾翼详细参数计算采用典型飞机的尾翼容量系数法，本方案尾容系数初步为。02VTC【尾容系数尾翼升力系数尾翼产生的力矩系数】从机翼俯视图上看，可得/2TAN1/2TAN80924VTLBM由可得VTWSLCB垂尾面积2020394VTWSBM展弦比20；垂尾后掠角45翼展2038M垂尾平均气动弦长03758VTSM根梢比05；根弦长005，梢弦长0025重尾视图（4）舵面设计小型无尾布局电动无人机大多采用升降副翼混合控制实现俯仰和偏航控制，一般在机翼后缘布置舵面，利用控制系统实现副翼和升降舵的功能。舵面设计在前期阶段不重要，要根据后期操纵性能来进一步修改。对速度不高的飞机，舵面相对面积约取为0304。副翼面积相对机翼面积一般57；副翼相对弦长约为2025；一般副翼偏角，不超过25。本方案无人机的升降副翼布置在翼尖。弦长取机翼平均弦长的12，为0025M升降副翼面积为02250011M2展长044M后缘上下偏角257重心位置确定由于本方案飞机起飞着陆时需人工操纵，所以需要有较好的静定性。初步确定纵向静稳定裕度为。即,08LMC,LCMCCXB其中，为全机重心位置与全机焦点位置间的距离与平均气MCX动弦长之比。重心位置由内部装载布置确定，焦点则由气动布局确定。利用AAA飞机设计软件计算无人机的焦点位于机翼根弦前缘点后距离。（使用软件来确定飞焦点）对于本方案的飞翼布局，机翼焦点可近似为全机的焦点，具体确定后掠翼焦点的方法如下05805802101327TAN04935829067CXXZBDSZDZZ9316CXM从而可以求得0978023198XM即重心位于机翼根弦前缘点后01982M，重心位置确定。对于小型电动无人机，其重心位置可以根据操稳特性计算后，通过移动电池位置来调整。8三维模型建立及内部装载布置（1）三维模型本方案三维数学模型的建立使用CATIA完成。三面图前视图俯视图侧视图效果图（2）内部装载布置电动无人机机身内部装载有电池、自动驾驶仪、数据传输设备、图象传输设备、侦察设备。在机翼中段的分置见图所示。内部装载布置9无人机气动特性分析可以工程解析法计算，也可以涡格法ALV软件计算。估计AAA也能计算。气动特性包括飞机的升力特性、阻力特性和力矩特性。工程估算分析结果将作为性能计算的输入，用于飞行性能的分析。（1）全机升力特性分析A全机升力线斜率对于低速常规的直线边梯形机翼，机翼升力线斜率可用下式,LAWC估算,21/2222TAN4354081/509769/8LAWACKRAD其中，为翼型升力线斜率，为展弦比。/LACLA2MA由于全机没有平尾，因此，机翼的升力线斜率就是全机的升力线斜率。,4081/LAWCRADB全机零升力系数亚音速时，对于具有等翼型、线性扭转角分布的机翼，其零升迎角可用下式估计0,00,8018COMLWWINO其中，翼型零升迎角；每度扭转角引起的零升迎角0W0增量；压缩性修正因子。INCOM,0【注速度低空气压缩性不考虑，由于飞翼布局忽略机翼扭转】【注机翼的零升迎角不是全机的零升迎角，因为存在安装角。】零升力系数【零升力系数零升迎角升力线斜率】通常机翼的零升力系数为机翼零升力系数与平尾升力系数之和。0,0,LWLHC本方案无平尾。则0,0,284802151LWLWCI其中，为机翼安装角。WIC全机零升迎角有了全机的零升力系数以及升力线斜率，可以求得零升迎角02152848LCRADO实际上，有了机翼的零升迎角，由于无尾翼，但机翼存在安装角，可知，全机的零升迎角为28度。D升力系数线性的迎角范围初步估计可采用下式927WIOE最大升力系数及失速迎角在雷诺数差不多的情况下，干净机翼的最大升力系数通常取由二维翼型数据确定的翼型最大升力系数的90左右。机翼后掠使最大升力系数减小，由无后掠机翼的最大升力值乘以1/4弦长处的后掠角的余弦得到下式MAX,AX,1409COS092COS809615LWLWC对于本方案巡航状态全机的雷诺数为300000，翼型在此雷诺数下的最大升力系数由翼型选择可知。机翼最大升力对应的迎角MAXAX,0,MAX09615801274LLWCLCO其中，分离引起的迎角增量为【通过查表，查什么表】AXLCO全机的最大升力系数MAXAX,09615LWC全机的最大升力系数对应的失速迎角MAXAX,27LLCWIOF全机升力曲线（有点问题，没考虑安装角）（2）全机阻力特性分析阻力分为零升阻力和升致阻力，对于低速电动机，零升阻力主要为压差阻力和摩擦阻力。A全机零升阻力系数部件构成法是用平板摩擦阻力系数以及形状因子FF来估算飞FC机每一部件的亚音速零升阻力。然后用因子Q来考虑部件阻力的相互干扰，总的部件阻力等于浸湿面积、FF和Q的乘积。【可WETSF以用于计算机翼、平尾垂尾等的零升阻力系数】采用部件构成法，亚音速飞机零升阻力估算公式为0FWETDCFQS其中，部件表面摩擦系数；部件形状因子；QFC为干扰因子；部件的浸湿面积；参考面积WETSS对于大部分飞机，流过部件的气流可认为是紊流，但对于低雷诺数飞行器，气流大部分可能是层流。一般地，当雷诺数在50万时，气流流过平板会从层流变为紊流，转捩点位置为555104672100488EVRXV机翼机翼雷诺数为551023E/71046VLV机翼处于层流层和紊流层的摩擦阻力系数为层流5MIN1328/RE1328/2FLARC紊流25820651058206545/LGE473198FTURBLENTOMA从而，机翼的平板摩擦阻力系数为MIN10458205948045FFFCXLARXCTURBLENT机翼的形状因子40180281028061103COS/9959COS270MMTTFMAXCC【注这里近似将】025M机翼浸湿面积与参考面积比【从CATIA三维设计图/27WETS中测量，S为三视图外露平面面积】干扰因子【由于干扰较小】1Q机翼零升阻力系数0,045921702FWETDWCFS垂尾MIN10458205948045FFFCXLARXTURBLENT418002806103COS/5152945079MTFMAXCC【注意参考面积需统一】03/28159WETHETS干扰因子【由于干扰较小】Q垂尾零升阻力系数50,2045829108210FWETDHCFS总的废阻力还包括飞机特殊部件的杂项阻力，如襟翼、固定式起落架、上翘的后机身及底部面积，并且把估计的漏泄及鼓包阻力一起加到总阻力中。杂项阻力可以使用大量的经验图表及公式分别确定，然后把结果加到上面已确定的零升阻力中去。本机的杂项阻力取为总废,DMISCC阻力的4。【飞翼布局全机零升阻力系数可用机翼零升阻力系数近似，这里不计两个垂尾的零升阻力系数】全机零升阻力系数0124016DCB全机升致阻力系数可以采用涡格法求诱导阻力因子。也可以用解析法对于后掠翼飞机0680154615COS32718LEEA17353058KAEC全机极曲线2016DLCC升阻比最大时，/739803LD最大升阻比MAX/12LD0020406081121416000100200300400500600700800901011012013014015016017018CDCL（3）全机俯仰力矩特性分析全机俯仰力矩由机翼和尾翼俯仰力矩组成，但飞翼布局没有平尾，则机翼俯仰力矩则为全机俯仰力矩。可以用涡格法计算。A全机零升俯仰力矩系数机翼计算机翼俯仰力矩系数的参考面积为机翼面积，参考长度取为平均气动弦长。力矩参考点取为设计重心处。机翼的零升俯仰力矩系数一部分来源于翼型，另一部分来0,MWC源于机翼的扭转和平面形状。当整个机翼处于零升力迎角时，局部剖面的升力并不都为零，会引起附加的零升力矩。当展弦比大于25，后掠角小于线性扭转时，可用下式估算机翼低速时的450,MWC20,1/4002COS851895S2MRTMMCACCB机翼俯仰力矩系数对升力导数机翼的俯仰力矩系数为,0MWMACLCXC因此，机翼俯仰力矩系数对升力导数为【如何查表获得】08ACMRREFACREFLADCXX垂尾没有垂尾全机零升俯仰力矩系数为,0890MWMHWLCCC10无人机操纵性与稳定性分析飞机稳定性是飞机设计的一项重要指标。在评价飞机稳定性过程中主要是通过飞机的气动导数来判断。【要求理解各导数物理意义，掌握一种计算方法】（1）横航向静导数计算横航向静导数是指飞机因侧滑而引起的横向力、滚转力矩和偏航力矩等系数对侧滑角的导数。A横向力对侧滑角的导数YC翼身组合体主要包括机翼、机身、翼身干扰、机翼上反角等。小迎角YC时，机翼贡献是小量，机身贡献包括干扰，则,YWFYYFC对于亚音速飞机，,057305731086YWC,12289FOYFIACKS其中，为翼身干扰因子、为机身位于横向力作用点的横IFOA截面积。S为机翼面积【如何得来】,08651067YWFYYFC垂尾本机机翼两端设置垂尾，此处按照双垂尾计算,0322981754692YVWFHVYVYEFEFCS式中机翼机身平尾对双垂尾侧滑导数的干扰因子，VEFYF,双垂尾横向力系数导数的有效值，VEFYC,单垂尾面积，S为机翼面积；从而，全机的横向力对侧滑角的导数为,067490136YWFYVYWYFYVCCB滚转力矩对侧滑角的导数L滚转力矩主要来源为翼身组合体和垂尾，其中，机翼占主体。翼身组合体亚音速小迎角范围，翼身组合体的滚转力矩导数由下式给出,141412LLLLLLWFLFMFMLZWLWACCCKKCTGT【不知如何计算得来】垂尾重心位于机翼根弦前缘点后01982M。【没明白是迎角可能因为机体坐标系与气流坐标系不一致的原因。】,COSIN0524SI4690181VVLVYVZLCB全机滚转力矩对侧滑角的导数,04LLWFLVCC偏航力矩系数对侧滑角的导数N偏航力矩系数对侧滑角导数主要由机翼、机身、垂尾等部件产生。翼身组合体小迎角下，可以忽略机翼产生的偏航力矩。,0315480757297NWFNFWFSRCLKB其中，为机身侧面积、为机身长度。FSFL垂尾,COSIN024COS05IN2691031VVNVYVLZCB全机偏航力矩对侧滑角的导数（1/RAD）,0570346NWFNVC（2）纵向动导数计算包括俯仰角速度引起的升力系数动导数和迎角角度率引起的升力系数动导数。A俯仰角速度引起的升力系数动导数【俯仰阻尼导数】一般飞机包括平尾和机翼两部分。飞翼布局只有机翼，机身可忽略。机翼对的贡献，考虑机身存在时可按下式估算LQC12052025761840840253391/AWFLAWCGFLQFXCXRD焦点到翼根前缘距离0923016CXM机翼平均气动弦长前缘到翼根前缘距离014M机翼平均气动弦长前缘到焦点距离0216101400761MFX重心位于机翼根弦前缘点后01982M机翼平均气动弦长前缘到重心距离0198201400582MCG全机俯仰角速度引起的升力系数动导数为291/LQWFCRADB俯仰角速度引起的引起俯仰力矩系数动导数一般飞机包括平尾和机翼两部分，而且平尾是主要的，机翼贡献为平尾的10，对于飞翼布局，无平尾，忽略不计。0MQCC迎角变化率引起的引起升力系数动导数机翼对的贡献，考虑机身存在时可按下式估算LAC32312WFLALALAKCM【不知道符号是什么意义，怎么算】LAWFCD迎角变化率引起的引起俯仰力矩系数动导数又称为时差导数，主要来自于平尾，因此0LMC（3）横向动导数计算横向动导数计算包括滚转角速度引起的横向力系数动导数、滚转角速度引起的滚转力矩系数动导数、滚转角速度引起的偏航力矩系数动导数。A横向力系数动导数此动导数在动态分析中不重要，忽略。B滚转角速度引起的滚转力矩系数动导数【滚转阻尼导数】对于一般飞机，主要由机翼、平尾、垂尾产生。机翼DRAGLPLPWCLAWCLLPLPWK,0,00,SIN134231SINLPLPWWDLWCDILPDRAGL0,25【不知道怎么算，参数哪里来】垂尾垂尾贡献也是一个小量。估算表达式为205469031LPVYVCZB平尾无平尾。全机滚转角速度引起的滚转力矩系数动导数LPLVLPWCC滚转角速度引起的偏航力矩系数动导数【交叉导数】主要垂尾和机翼。机翼机翼的贡献主要是左右翼之间有效迎角引起的阻力差而引起的。公式太复杂了，不适合用解析法。垂尾VYVVVVNPVCAZLALZBCSINCOSINCO2全机8LLPLVLPWCD偏航角速度引起的横向力系数动导数不太重要，忽略。E偏航角速度引起的滚转力矩系数导数【交叉导数】（4）航向动导数计算航向动导数计算包括偏航角速度引起的横向力系数动导数、偏航角速度引起的滚转力矩系数动导数、偏航角速度引起的偏航力矩系数动导数。建议采用软件设计。（5）操纵导数计算A副翼操纵导数副翼偏转引起的横向力系数导数【忽略】副翼偏转引起的滚转力矩系数导数副翼偏转引起的偏航力矩系数导数B升降舵操纵导数升降舵偏转引起的升力系数导数升降舵偏转引起的俯仰力矩系数导数C方向舵操纵导数方向舵偏转引起的横向力系数导数方向舵偏转引起的滚转力矩系数导数方向舵偏转引起的偏航力矩系数导数（6）全机稳定性分析全机稳定性导数LCDMCLMCLQMQC纵向迎角迎角角速率俯仰角速率YLNLPNPLRNR横航向侧滑角滚转角速率偏航角速率A纵向静稳定性与动稳定性静稳定性要求对于小型电动无人机，静裕度取值一般在1015。动稳定性要求按飞行品质规范GJB18586对长短周期评价。这里采用CATIA结构装载建模，定义材料属性和设备重量属性，由CATIA测量41得到飞机的质量惯性矩。飞机六个转动惯量飞机纵向动稳定性由扰动运动的典型模态表示，通常，它由两个快慢相差较大的振荡运动组成，按振荡周期短、长分别命名为短周期模态和长周期模态。在扰动运动的最初阶段，主要特征是以迎角和角速度变化为代表的短周期运动，飞机速度基本保持不变；而在扰动运动的后一阶段，主要特征是以速度和俯仰角变化为代表的长周期运动，飞机迎角基本保持不变。利用飞机纵向小扰动运动方程得到的特征方程有两对共轭复根，一对大的和一对小的共轭复根（或两个小实根）。大根描述快速运动，对应于飞机短周期运动，小根对应于长周期运动（沉浮运动）。飞机纵向小扰动运动方程2220SPSPPPSS短周期模态对应短周期模态的特征方程220SPSPS特征方程的一对大的共轭复根为21,21SPSSPI无阻尼振荡频率MLYASPKCISCQ式中自由流动压；。QQCMKL,AEQSC4阻尼比24LAEQMQYSPSPCVSI振荡周期21SPSPT倍增时半衰期1/2,0693SPSP短周期模态对应长周期模态的特征方程220PPS特征方程的一对小的共轭复根为23,41PPI无阻尼振荡频率EQPVG2阻尼比振荡1/2PDLC周期2/1PPT半衰期1/2,0693PPB横航向静稳定性与动稳定性静稳定性要求具有航向静稳定性的飞机其，具有横向0NC静稳定的飞机；0LC动稳定性要求按飞行品质规范GJB18586进行评价。飞机的横航向动态特性用三个典型扰动运动模态来表征，一般以大的实根（一般为负）表示快速滚转运动模态；以小实根（可正可负）表示缓慢螺旋运动模态；以一对共轭复根表示荷兰滚运动模态。其中滚转模态和螺旋模态为两个非周期模态，另一个是振荡模态即荷兰滚模态。利用飞机横航向小扰动运动方程组，可得到特征方程220RSDDSS滚转模态滚转模态特性参数通常用滚转模态时间常数来表示，与特RTR征方程中滚转模态大的负实根成反比。当横向和航向耦合不显著R时，【如果不是大实根，则不行】LPEQXRCSBVIT2螺旋模态当0时，螺旋模态稳定。RRLN式中；2EQLXSVBCILRXEQRCIBSVL42；。NZEQIN2NRZEQRIN2荷兰滚模态荷兰滚模态是一种具有偏航、滚转和侧滑3个运动同时存在的短周期振荡的二阶模态。对应荷兰滚模态的特征方程为220DDS特征方程的根21,21DDDI式中荷兰滚模态无阻尼振荡频率，；D/2DN荷兰滚模态阻尼比，；D/REQDDYV荷兰滚模态总阻尼，；D/2REQDN；YEQCMSVY2。YREQRB4【分析根的性质不能变】（7）全机操纵性分析全机操纵导数LECME纵向升降舵YRLRNRCYLCN横航向方向舵副翼【操纵性如何评价找相关文献。】说明升降舵有足够的操纵能力1MEC说明副翼有足够的操纵能力1LC说明方向舵有足够的操纵能力NR11无人机飞行性能分析（1）平飞需用推力速度取832M/S，间隔2M/S。速度平飞所需升力系数阻力系数升阻比推力（阻力）800001750102368739092037510000011201010481068941408712000007778005741355901110614000005715003701542840976016000004375002721610070935318000003457002191576800955020000002800001891479701017722000002314001711352581113324000001945001601218371236026000001657001521089781381828000001429001479725415484300000012450014386844173405101520253008112141618222平飞需用推力曲线（2）最大平飞速度最大平飞速度可通过可用需用推力曲线来确定，也可以通过可用需用功率曲线来确定。电池最大电压185V电机参数型号KV空载电流I0（A）内阻RM）重量（G）HIMAXHC281206506500360285602空速2634M/S,间隔2M/S。由电机的输出扭矩和电动机的电压0TQKI可得到MVNUIRK00959518/637943562VVTMMNNUUKKQIIRRNN由上式可画出电压为185V时的QN曲线080911112131415RPM/MINQNM最大电压下电机转速扭矩曲线电机最大转速空载转速MAX09435/MIN9/7NRPRS计算收敛条件【与之前不同，之前是飞机功率与螺MPROQ旋桨功率相近，这里不要求功率】【从最大转速依次减小，试算】空速26M/S取N180R/S10800R/MIN,J26/18002540807108电机输出扭矩07918043509MQNM螺旋桨需用扭矩2525250481504801PPROQCNDNM，说明电机在此转速下扭矩不够，带不动螺旋桨转动。MR因此，减小转速。取N175R/S10500R/MIN,J26/17502540807312电机输出扭矩0791504301MQNM螺旋需用扭矩25252541570480931PPROQCND，可以适当减小转速。MR取N1775R/S10650R/MIN,J26/177502540807209电机输出扭矩079165043102MQNM螺旋需用扭矩2525254157048097PPROQCND。MPRO螺旋桨输出功率以下是等价的35352072681027580244802OUTPROTPQNWCD空速28M/S取N180R/S10800R/MIN,J28/18002540807655电机输出扭矩07918043509MQNM螺旋需用扭矩25252541580480PPROQCNDMPROQ螺旋桨输出功率2072689031749OUTPROPNW空速30M/S取N185R/S11100R/MIN,J30/1850254080777电机输出扭矩0791043506MQNM螺旋需用扭矩252525381504807PPROQCND，说明电机在此转速下扭矩不够，带不动螺旋桨转动。MR因此，减小转速。取N1825R/S10950R/MIN，J30/182502540808090电机输出扭矩079150430785MQNM螺旋需用扭矩252525