发明名称 一种无人机久航飞行参数确定方法
摘要 本发明公开一种无人机久航飞行参数确定方法,从实际使用角度出发,以寻找飞机最小平飞需用功率对应的发动机最小耗油量为目标,在飞机、螺旋桨以及发动机的工作点互相匹配的前提下使耗油量最小,达到了增加无人机的航时的目的。本发明原理简单,便于工程应用。
申请公布号 CN101634612A 申请公布日期 2010.01.27
申请号 CN200910091068.4 申请日期 2009.08.21
申请人 北京航空航天大学 发明人 李仪;向锦武;马东立;林海英;郭伟
分类号 G01M19/00(2006.01)I 主分类号 G01M19/00(2006.01)I
代理机构 北京永创新实专利事务所 代理人 周长琪
主权项 1、一种无人机久航飞行参数确定方法通过以下四步实现:步骤一:确定飞机最小平飞需用功率Pxy·min;飞机的极曲线一般为:Cx=Cx0+Cx1Cy+Cx2Cy2 (1)其中Cy为升力系数,Cx为阻力系数,Cx0、Cx1、Cx2为方程系数,可由吹风试验得到;假设飞机重量为G,飞行高度为H,根据简单推力法:<maths id="math0001" num="0001" ><math><![CDATA[ <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>G</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>&rho;V</mi> <mn>2</mn> </msup> <msub> <mi>SC</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mi>T</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>&rho;V</mi> <mn>2</mn> </msup> <msub> <mi>SC</mi> <mi>x</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>]]></math></maths>式中,S为机翼面积、,ρ空气密度;由公式(2)得到飞机平飞所需用的功率Pxy:<maths id="math0002" num="0002" ><math><![CDATA[ <mrow> <msub> <mi>P</mi> <mi>xy</mi> </msub> <mo>=</mo> <mi>T</mi> <mo>&CenterDot;</mo> <mi>V</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>&rho;V</mi> <mn>3</mn> </msup> <msub> <mi>SC</mi> <mi>x</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>]]></math></maths>对公式(3)求导,并令一阶导数为零,即:<maths id="math0003" num="0003" ><math><![CDATA[ <mrow> <mfrac> <msub> <mi>dP</mi> <mi>xy</mi> </msub> <mrow> <mi>d</mi> <msub> <mi>C</mi> <mi>y</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>dP</mi> <mi>xy</mi> </msub> <mi>dV</mi> </mfrac> <mfrac> <mi>dV</mi> <msub> <mi>dC</mi> <mi>y</mi> </msub> </mfrac> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>]]></math></maths>即可得到给定高度H下飞机最小平飞需用功率P*对应的Cy*;根据公式(2)可以求出对应的速度V*;处于飞行安全考虑,一般对各个飞行高度的飞行速度V范围都进行了限定,假设此速度范围为[Vmin,Vmax];如果V*∈[Vmin,Vmax],则取飞行速度V=V*;如果V*≤Vmin,则取飞行速度V=Vmin;如果V*≥Vmax,则取飞行速度V=Vmax;通过(2)、(3)式求出V对应的功率P,则Pxygmin=P;步骤二:确定发动机功率P;选取一组发动机节风门位置TP=[TP1L TPiL TPr]和一组转速值n=[n1L njL ns],根据发动机特性曲线P=f2(H,TP,n)、q=f3(H,TP,n),求出每一个节风门TPi(i=1,2,…,r)对应的所有转速下的发动机功率P和耗油量q:<maths id="math0004" num="0004" ><math><![CDATA[ <mrow> <mi>P</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>11</mn> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mn>1</mn> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>ij</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>is</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>rj</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>rs</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>]]></math></maths><maths id="math0005" num="0005" ><math><![CDATA[ <mrow> <mi>q</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>q</mi> <mn>11</mn> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>q</mi> <mi>ij</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>q</mi> <mi>is</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>q</mi> <mi>rj</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>q</mi> <mi>rs</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>]]></math></maths>步骤三:确定螺旋桨效率η与桨距角θ;根据步骤一中确定的飞行速度V和步骤二选取的转速n以及发动机输出功率P,可求出所有节风门TPi对应的每一个转速nj下的前进比λ和功率系数Cp,由螺旋桨特性曲线η=f1(Cp,λ,θ),可插值求出对应的螺旋桨桨距角θ和螺旋桨效率η:<maths id="math0006" num="0006" ><math><![CDATA[ <mrow> <mi>&theta;</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>&theta;</mi> <mn>11</mn> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&theta;</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&theta;</mi> <mrow> <mn>1</mn> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&theta;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&theta;</mi> <mi>is</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&theta;</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&eta;</mi> <mi>rj</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&theta;</mi> <mi>rs</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mi></mi> <mo>)</mo> </mrow> </mrow>]]></math></maths><maths id="math0007" num="0007" ><math><![CDATA[ <mrow> <mi>&eta;</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>&eta;</mi> <mn>11</mn> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&eta;</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&eta;</mi> <mrow> <mn>1</mn> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&eta;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&eta;</mi> <mi>ij</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&eta;</mi> <mi>is</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&eta;</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&eta;</mi> <mi>rj</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>&eta;</mi> <mi>rs</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>]]></math></maths>步骤四:确定最小耗油量qxh·min;根据步骤三求出的螺旋桨效率η,计算所有节风门TPi对应的每一个转速下的发动机可用功率<maths id="math0008" num="0008" ><math><![CDATA[ <mrow> <msub> <mi>P</mi> <mi>ky</mi> </msub> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mn>11</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mn>1</mn> <mi>j</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mn>1</mn> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mi>ij</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mi>is</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mi>rj</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mi>ky</mi> <mo>&CenterDot;</mo> <mi>rs</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>11</mn> </msub> <msub> <mi>&eta;</mi> <mn>11</mn> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> </msub> <msub> <mi>&eta;</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mrow> <mn>1</mn> <mi>s</mi> </mrow> </msub> <msub> <mi>&eta;</mi> <mrow> <mn>1</mn> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&eta;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>ij</mi> </msub> <msub> <mi>&eta;</mi> <mi>ij</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>is</mi> </msub> <msub> <mi>&eta;</mi> <mi>is</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&eta;</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>rj</mi> </msub> <msub> <mi>&eta;</mi> <mi>rj</mi> </msub> <mi>L</mi> </mtd> <mtd> <msub> <mi>P</mi> <mi>rs</mi> </msub> <msub> <mi>&eta;</mi> <mi>rs</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>]]></math></maths>根据飞机最小平飞需用功率Pxy·min插值求出Pxy·min在每一个转速nj时,发动机发出Pxy·min对应的节风门[TP1L TPjL TPs],再由[TP1L TPjL TPs]以及η、θ、q插值求出螺旋桨效率[η1L ηjL ηs]、螺旋桨桨距[θ1L θjL θs]、耗油量[q1L qjL qs],取qxh·min=min[q1L qjL qs],则对应的飞行参数TP、θ、n即是最省油状态下的飞行参数;假设经过时间Δt,飞机重量变为G′=G-ΔG=G-q·Δt,令G=G′,重复步骤一到步骤四,就可以确定出下一重量G和高度H时的飞行参数,即可确定整个飞行过程中的飞行参数。
地址 100083北京市海淀区学院路37号