CN107688705A - 基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法 - Google Patents
基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法 Download PDFInfo
- Publication number
- CN107688705A CN107688705A CN201710744273.0A CN201710744273A CN107688705A CN 107688705 A CN107688705 A CN 107688705A CN 201710744273 A CN201710744273 A CN 201710744273A CN 107688705 A CN107688705 A CN 107688705A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msup
- mfrac
- msqrt
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Nitrogen And Oxygen Or Sulfur-Condensed Heterocyclic Ring Systems (AREA)
Abstract
基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,本发明涉及旋翼系统诱导流场的轴向诱导速度计算方法。本发明的目的是为了解决现有计算方法难以同时兼顾快速计算和快速收敛的要求的问题。过程为:一、获取旋翼平面上方空间内测试点集在椭圆坐标系下的坐标;二、设定最大谐波参数,根据最大谐波参数求解相应矩阵参数,并根据 及相应矩阵参数,求解三、求解四、求解旋翼平面上方测试点处诱导流场的轴向诱导速度;五、根据重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度;六、根据优化后重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度。本发明用于旋翼系统诱导速度计算领域。
Description
技术领域
本发明涉及旋翼系统诱导流场的轴向诱导速度计算方法。
背景技术
伴随着直升机,大型风机发电设备的不断发展,为满足相应的系统的设计指标和要求,深入了解系统动态特性,准确的模型已经显得越来越重要。而在这类系统中,无疑旋翼系统的模型和动态将起到重要的影响,而准确计算旋翼的动态特性的关键之一是精确表达旋翼桨盘及其周围的空气诱导速度,它与旋翼存在着作用与反作用的关系。当旋翼发生扰动时,旋翼气动载荷的变化必然引起诱导速度的变化,而这一变化又会反过来影响旋翼的气动载荷,进而影响旋翼的运动特性。旋翼系统坐标系的建立如图1所示。
准确的计算空气诱导速度场是旋翼系统准确建模的关键问题,国内外很多旋翼系统的建模过程中,采用动量定理或者Pitt-Peters等较低阶次的动态入流模型,导致建模误差较大,无法准确地反映旋翼系统相关特性。虽然也有高阶次的动态入流模型提出和应用,如Peters-He模型和Peters-Morillo模型,但这两个模型均有其固有的不足和缺陷,制约了他们的应用。而运用大型流体计算软件计算诱导流场,往往边界条件设置复杂,计算量大,耗时耗力且无法保证实时性。与此同时,现代对直升机,大型风力发电机的要求越来越高,疲劳特性和机动特性等应该在设计阶段就得到预先分析、计算和检验,必然需要较高精度的模型予以验证;以上问题导致现有的计算方法难以同时兼顾快速计算和快速收敛的要求。
发明内容
本发明的目的是为了解决现有计算方法难以同时兼顾快速计算和快速收敛的要求的问题,而提出基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法。
基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法具体过程为:
步骤一、获取旋翼平面上方空间内测试点集在椭圆坐标系ν,η,下的坐标;
步骤二、设定最大谐波参数N(设定范围一般可选为8-15),根据最大谐波参数求解相应矩阵参数,并根据正弦分量的动态压力系数余弦分量的动态压力系数及相应矩阵参数,求解
相应矩阵参数为[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls];
其中[Mc]、[Ms]代表余弦和正弦分量的质量矩阵,[Dc]、[Ds]代表余弦和正弦分量的阻尼矩阵,[Lc]为余弦分量的关联系数矩阵、[Ls]为正弦分量的关联系数矩阵, 为Peters-Morillo模型展开系数;
步骤三、根据(包含奇偶数量)求解(只包含奇数量);
其中为混合模型中分别对应余弦和正弦奇数量的系数;
步骤四、根据和步骤一获取的测试点集在椭圆坐标系ν,η,下的坐标构建混合模型,基于混合模型求解旋翼平面上方测试点处诱导流场的轴向诱导速度
步骤五、根据混合模型计算得到的旋翼平面上方测试点处的轴向诱导速度和Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度建立重构模型,根据重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度vz;
步骤六、利用基因算法、蚂蚁算法或粒子算法对重构模型中的系数进行优化,使重构模型误差达到最小,最终求出优化后重构模型,根据优化后重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度v′z。
本发明的有益效果为:
本发明方法通过建立旋翼平面上方测试点集在椭圆坐标系ν,η,下的坐标,整合了现有的Peters-He模型与Peters-Morillo模型,求解相应矩阵参数,构建混合模型,基于混合模型求解旋翼平面上方测试点处诱导流场的轴向诱导速度并根据旋翼平面上方测试点处诱导流场的轴向诱导速度和Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度建立重构模型,根据重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度vz,对重构模型中的a、c进行优化,使误差e达到最小,并求出优化后重构模型,根据优化后重构模型最终得到旋翼平面上方测试点处诱导流场的轴向诱导速度v′z,完成本发明基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,达到了同时兼顾快速计算和快速收敛的要求。
结合实施例一选取a=19.8,c=0.3,作为一组假定优化后的参数,比较改进后模型的优越性。其中图4-图9为旋翼系统z=0平面上,不同角度下,在静态压力(ω=0)和动态压力(ω=4)下的轴向诱导速度对比图。可以看出在-1<x<1区间上,重构模型比Peters-Morillo模型接近准确理论解,尤其在大角度时(χ=89°),改进效果更明显。在边沿处,可以有效计算轴向诱导速度。在-1<x<1区间外与Peters-Morillo模型收敛速度基本一致。图10-图13为旋翼系统z=-0.5平面上不同角度下,在静态压力(ω=0)和动态压力(ω=4)下的轴向诱导速度对比图。可以看出在旋翼平面上方,重构模型和Peters-Morillo模型具有相似的收敛速度和准确性。以上结果也充分说明,本发明的计算方法在大角度、旋翼边缘有更好的收敛特性。
附图说明
图1为建立的旋翼系统坐标系示意图;
图2为重构模型计算算法流程图;
图3为优化过程流程图;
图4为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=0及时轴向诱导速度对比结果图(χ=0°,静态压力ω=0);x,y,z为笛卡尔坐标系下的坐标,为对应模态的动态压力系数,ω为动态压力的频率,Peter-Morrilo模型为(11)式,χ指尾流倾角;
图5为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=0及时轴向诱导速度对比结果图(χ=45°,静态压力ω=0);
图6为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=0及时轴向诱导速度对比结果图(χ=89°,静态压力ω=0);
图7为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=0及时轴向诱导速度对比结果图(χ=0°,动态压力ω=4);
图8为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=0及时轴向诱导速度对比结果图(χ=45°,动态压力ω=4);
图9为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=0及时轴向诱导速度对比结果图(χ=89°,动态压力ω=4);
图10为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=-0.5及时轴向诱导速度对比结果图(χ=0°,静态压力ω=0);
图11为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=-0.5及时轴向诱导速度对比结果图(χ=89°,静态压力ω=0);
图12为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=-0.5及时轴向诱导速度对比结果图(χ=0°,动态压力ω=4);
图13为Peter-Morrilo模型方法、本发明重构模型方法、准确理论解方法在y=0,z=-0.5及时轴向诱导速度对比结果图(χ=89°,动态压力ω=4)。
具体实施方式
具体实施方式一:本实施方式为基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法具体过程为:
步骤一、获取旋翼平面上方空间内(是旋翼平面上方的半球形空间区域内,半径可以无穷大,旋翼平面把空间分为上半球和下半球空间)测试点集在椭圆坐标系ν,η,下的坐标;
步骤二、根据精度要求和实际情况,设定最大谐波参数N(设定范围一般可选为8-15),根据最大谐波参数求解相应矩阵参数,并根据正弦分量的动态压力系数余弦分量的动态压力系数及相应矩阵参数,结合方程(4),(5)求解
相应矩阵参数为[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls];
其中[Mc]、[Ms]代表余弦和正弦分量的质量矩阵,[Dc]、[Ds]代表余弦和正弦分量的阻尼矩阵,[Lc]为余弦分量的关联系数矩阵、[Ls]为正弦分量的关联系数矩阵, 为Peters-Morillo模型展开系数;
步骤三、根据(包含奇偶数量)求解(只包含奇数量);
其中为混合模型中分别对应余弦和正弦奇数量的系数;
步骤四、根据和步骤一获取的测试点集在椭圆坐标系ν,η,下的坐标构建混合模型,基于混合模型求解旋翼平面上方测试点处诱导流场的轴向诱导速度
步骤五、根据混合模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度和Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度建立重构模型,根据重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度vz;如图2所示;
步骤六、利用基因算法、蚂蚁算法或粒子算法对重构模型中的系数进行优化,使重构模型误差达到最小,最终求出优化后重构模型,根据优化后重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度v′z。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中获取旋翼平面上方测试点集在椭圆坐标系ν,η,下的坐标;具体过程为:
以旋翼中心为原点,建立如图1所示直角坐标系,χ代表尾流倾角,利用如下公式,将笛卡尔坐标系下的坐标转换为椭圆坐标系下的坐标
其中r′2=x2+y2+z2,r′为测试点与原点间的距离;x,y,z代表测试点集在笛卡尔坐标系下的坐标,ν,η,代表测试点集在椭圆坐标系下的坐标;其中η=0代表了旋翼所在的圆形区域。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤二中根据精度要求和实际情况,设定最大谐波参数N(设定范围为8-15),根据最大谐波参数求解中相应矩阵参数,并根据正弦分量的动态压力系数余弦分量的动态压力系数及相应矩阵参数,结合方程(4),(5)求解
相应矩阵参数为[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls];
其中[Mc]、[Ms]代表余弦和正弦分量的质量矩阵,[Dc]、[Ds]代表余弦和正弦分量的阻尼矩阵,[Lc]为余弦分量的关联系数矩阵、[Ls]为正弦分量的关联系数矩阵, 为Peters-Morillo模型展开系数;具体过程为:
在此,我们通过如下公式,分别计算Peters-Morillo模型的余弦和正弦分量:
其中代表系数的一阶导数,[]-1表示对关联系数矩阵求逆运算;且[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls]矩阵的构造均遵从如下的关系:
其中m、n、j、r为矩阵元素所在位置的索引信息,且满足r≤j≤N,m≤n≤N,odd代表加和为奇数,even代表加和为偶数;m、n、j、r取值为正整数;
[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls]矩阵中各元素计算公式如下:
[Mc]和[Ms]矩阵中各元素的计算公式为([Mc]和[Ms]矩阵元素的计算是相同的):
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
r=m;j+r=odd;n+m=even;
为中间变量;
(A)!!代表数A的双阶乘,有如下定义
(A)!!=A(A-2)(A-4)…2,A=even
(A)!!=A(A-2)(A-4)…1,A=odd
0!!=1;(-1)!!=1;(-2)!!=∞;(-3)!!=-1;
(A)!!表示(n+m-1)!!、(n-m-1)!!、(n+m)!!、(n-m)!!、(j+m-1)!!、(j-m-1)!!、(j+m)!!、(j-m)!!;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置),
当r=m;j=n±1;j+r=odd;n+m=even;
当r=m;j=n±1;j+r=even;n+m=odd;
为中间变量;
(A)!!代表数A的双阶乘,有如下定义
(A)!!=A(A-2)(A-4)…2,A=even
(A)!!=A(A-2)(A-4)…1,A=odd
0!!=1;(-1)!!=1;(-2)!!=∞;(-3)!!=-1;
(A)!!表示(j+r-1)!!、(j-r-1)!!、(j+r)!!、(j-r)!!;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置),
当r=m;j+r=even;n+m=even;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r≠m;
其中[Mc]和[Ms]矩阵的差别在于[Ms]矩阵中不含有r=0的所有行和列。
[Dc]和[Ds]矩阵中各元素的计算公式为([Dc]和[Ds]矩阵元素的计算是相同的):
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r=m;j+r=odd;n+m=odd;
当r=m;j+r=even;n+m=even;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r=m;j+r=odd;n+m=even;
当r=m;j+r=even;n+m=odd;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r≠m;
其中为中间变量,
其中[Dc]和[Ds]矩阵的差别在于[Ds]矩阵中不含有r=0的所有行和列;
[Lc]矩阵和[Ls]矩阵中各元素的计算公式为:
且
式中Xm、X|m-r|、X|m+r|、为中间变量;
其中min(r,m)代表比较r、m,取两者中较小的整数;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r+m=odd;j+r=odd;n+m=odd;
当r+m=odd;j+r=even;n+m=even;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r+m=even;j+r=odd;n+m=odd;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r+m=even;j+r=even;n+m=even;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r+m=odd;j+r=odd;n+m=even;
当r+m=odd;j+r=even;n+m=odd;
(每个表达式下面的m、n、j、r的约束条件不同,决定了元素的不同位置)
当r+m=even;j+r=odd;n+m=even;
当r+m=even;j+r=even;n+m=odd;
其中为中间变量。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述最大谐波参数N取值范围一般为8-15。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中根据(包含奇偶数量)求解(只包含奇数量);具体过程为:
在中取奇数量的Peters-Morillo模型展开系数和将 转换为此处将由下式计算获得
其中,f≤N,为混合模型中的奇数项系数,为由Peters-Morillo模型展开系数中元素索引的奇数项组成;
其中中的元素由以下公式计算获得
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤四中根据和步骤一获取的旋翼平面上方测试点集在椭圆坐标系ν,η,下的坐标构建混合模型,基于混合模型求解旋翼平面上方测试点处诱导流场的轴向诱导速度具体过程为:
利用如公式(9)的混合模型计算旋翼平面上方测试点处诱导流场的轴向速度:
其中表示基于混合模型计算得到的旋翼平面上方测试点处诱导流场的轴向速度,表示第一类归一化勒让德函数,表示第二类归一化勒让德函数,i为复数的虚部;
该模型混合了Peter-He模型的基底,保留了该模型在大角度入流时更为准确的优点,同时克服了其在旋翼平面边缘(即桨叶外缘)无法有效计算的问题。混合模型(9)更好的解决了大角度时,模型在旋翼平面内的快速收敛问题。但由于该模型中混合了Peters-He模型的基底,故该模型仅能用于准确计算旋翼所在平面内轴向诱导速度的缺陷也将影响到新模型。因此在非旋翼平面内,其收敛是病态的,不能快速有效的计算旋翼上方的诱导流场。下面将利用重构模型,完成旋翼平面上方测试点处诱导流场的计算。
其它步骤及参数与具体实施方式一至四之一相同。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤五中根据混合模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度和Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度建立重构模型(10),根据重构模型(10)计算旋翼平面上方测试点处诱导流场的轴向诱导速度vz;具体过程为:
其中vz代表旋翼平面上方测试点处诱导流场的轴向诱导速度,表示由Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度,a,c为与重构模型收敛速度相关的系数;
其它步骤及参数与具体实施方式一至五之一相同。
具体实施方式七:本实施方式与具体实施方式一至六之一不同的是:所述步骤六中利用基因算法、蚂蚁算法或粒子算法对重构模型(10)中的系数(a、c)进行优化,使重构模型误差(e)达到最小,最终求出优化后重构模型(10),根据优化后重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度v′z;如图3所示。具体过程为:
上述模型中存在两组待定参数a,c,其决定了重构模型中和的权重,显然在旋翼平面内更为精确,而在平面外,计算结果将更依赖于为了进一步减少重构模型计算误差,可对参数a,c进行优化。首先,根据需要解决的问题和实际需求,选取兴趣区域。在区域内,随机均匀选取P个参照点,用以更全面反应区域内的特性。
利用卷积公式计算旋翼平面上方测试点中参照点(优化的时候选取几个参照点就行)处在给定压力系数时的轴向诱导速度的准确理论值作为与重构模型计算轴向诱导速度的参照;
利用重构模型(10)计算旋翼平面上方测试点中参照点处(优化的时候选取几个参照点就行)诱导流场的轴向诱导速度,并利用如下公式计算重构模型误差:
其中||e||2代表旋翼平面上方参照点的误差二范数的平方之和,为由重构模型计算得到的旋翼平面上方参照点k处的轴向诱导速度,为由卷积公式计算得到的旋翼平面上方参照点k处的轴向诱导速度的准确理论值;
利用基因算法、蚂蚁算法或粒子算法等优化算法,对下式进行参数优化
min||e||2subject to a>0,0<c≤1, (13)
subject to a>0,0<c≤1表示使a、c服从a>0,0<c≤1;
使得a,c的选取在旋翼平面上方(兴趣区内)最优,得到优化后重构模型(10);
利用优化后重构模型(10)计算旋翼平面上方其他测试点处的诱导流场的轴向诱导速度v′z。
其它步骤及参数与具体实施方式一至六之一相同。
采用以下实施例验证本发明的有益效果:
实施例一:
本实施例基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法具体是按照以下步骤制备的:
由于选取兴趣区域,参照点,参考压力系数,容许误差等选取的不同,优化参数a,c会存在一定的差别。在此仅选取a=19.8,c=0.3,作为一组假定优化后的参数,比较改进后模型的优越性。这里仅考虑无量纲化后的情形,且仅针对压力系数进行说明,实际中正余弦动态压力系数可以由不同频率分量的动态压力系数的线性组合计算得到。
图4-图9为旋翼系统z=0平面上,不同角度下,在静态压力(ω=0)和动态压力(ω=4)下的轴向诱导速度对比图。可以看出在-1<x<1区间上,重构模型比Peters-Morillo模型接近准确理论解,尤其在大角度时(χ=89°),改进效果更明显。在边沿处,可以有效计算轴向诱导速度。在区间外与Peters-Morillo模型收敛速度基本一致。
图10-图13为旋翼系统z=-0.5平面上不同角度下,在静态压力(ω=0)和动态压力(ω=4)下的轴向诱导速度对比图。可以看出在旋翼平面上方,重构模型和Peters-Morillo模型具有相似的收敛速度和准确性。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (7)
1.基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,其特征在于:所述方法具体过程为:
步骤一、获取旋翼平面上方空间内测试点集在椭圆坐标系ν,η,下的坐标;
步骤二、设定最大谐波参数N,根据最大谐波参数求解相应矩阵参数,并根据正弦分量的动态压力系数余弦分量的动态压力系数及相应矩阵参数,求解
相应矩阵参数为[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls];
其中[Mc]、[Ms]代表余弦和正弦分量的质量矩阵,[Dc]、[Ds]代表余弦和正弦分量的阻尼矩阵,[Lc]为余弦分量的关联系数矩阵、[Ls]为正弦分量的关联系数矩阵, 为Peters-Morillo模型展开系数;
步骤三、根据求解
其中为混合模型中分别对应余弦和正弦奇数量的系数;
步骤四、根据和步骤一获取的测试点集在椭圆坐标系ν,η,下的坐标构建混合模型,基于混合模型求解旋翼平面上方测试点处诱导流场的轴向诱导速度
步骤五、根据混合模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度和Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度建立重构模型,根据重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度vz;
步骤六、利用基因算法、蚂蚁算法或粒子算法对重构模型中的系数进行优化,使重构模型误差达到最小,最终求出优化后重构模型,根据优化后重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度v′z。
2.根据权利要求1所述基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,其特征在于:所述步骤二中设定最大谐波参数N,根据最大谐波参数求解相应矩阵参数,并根据正弦分量的动态压力系数余弦分量的动态压力系数及相应矩阵参数,求解具体过程为:
分别计算Peters-Morillo模型的余弦和正弦分量:
<mrow>
<mo>&lsqb;</mo>
<msup>
<mi>M</mi>
<mi>c</mi>
</msup>
<mo>&rsqb;</mo>
<mo>{</mo>
<mover>
<msubsup>
<mi>a</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>&CenterDot;</mo>
</mover>
<mo>}</mo>
<mo>+</mo>
<mo>&lsqb;</mo>
<msup>
<mi>D</mi>
<mi>c</mi>
</msup>
<mo>&rsqb;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mi>L</mi>
<mi>c</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&lsqb;</mo>
<msup>
<mi>M</mi>
<mi>c</mi>
</msup>
<mo>&rsqb;</mo>
<mo>{</mo>
<msubsup>
<mi>a</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>}</mo>
<mo>=</mo>
<mo>&lsqb;</mo>
<msup>
<mi>D</mi>
<mi>c</mi>
</msup>
<mo>&rsqb;</mo>
<mo>{</mo>
<msubsup>
<mi>&tau;</mi>
<mi>n</mi>
<mrow>
<mi>m</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>}</mo>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mi>M</mi>
<mi>s</mi>
</msup>
<mo>&rsqb;</mo>
<mo>{</mo>
<mover>
<msubsup>
<mi>b</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>&CenterDot;</mo>
</mover>
<mo>}</mo>
<mo>+</mo>
<mo>&lsqb;</mo>
<msup>
<mi>D</mi>
<mi>s</mi>
</msup>
<mo>&rsqb;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mi>L</mi>
<mi>s</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&lsqb;</mo>
<msup>
<mi>M</mi>
<mi>s</mi>
</msup>
<mo>&rsqb;</mo>
<mo>{</mo>
<msubsup>
<mi>b</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>}</mo>
<mo>=</mo>
<mo>&lsqb;</mo>
<msup>
<mi>D</mi>
<mi>s</mi>
</msup>
<mo>&rsqb;</mo>
<mo>{</mo>
<msubsup>
<mi>&tau;</mi>
<mi>n</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>}</mo>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中代表系数的一阶导数,[]-1表示对关联系数矩阵求逆运算;且[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls]矩阵的构造均遵从如下的关系:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中m、n、j、r为矩阵元素所在位置的索引信息,且满足r≤j≤N,m≤n≤N,odd代表加和为奇数,even代表加和为偶数;
[Mc]、[Ms]、[Dc]、[Ds]、[Lc]、[Ls]矩阵中各元素计算公式如下:
[Mc]和[Ms]矩阵中各元素的计算公式为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>M</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>m</mi>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>-</mo>
<mn>2</mn>
<mi>m</mi>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>r</mi>
<mo>=</mo>
<mi>m</mi>
<mo>;</mo>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
<mo>;</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
为中间变量;
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>m</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
</mrow>
</mfrac>
</mrow>
(A)!!代表数A的双阶乘,有如下定义
(A)!!=A(A-2)(A-4)…2,A=even
(A)!!=A(A-2)(A-4)…1,A=odd
0!!=1;(-1)!!=1;(-2)!!=∞;(-3)!!=-1;
(A)!!表示(n+m-1)!!、(n-m-1)!!、(n+m)!!、(n-m)!!、(j+m-1)!!、(j-m-1)!!、(j+m)!!、(j-m)!!;
当r=m;j=n±1;j+r=odd;n+m=even;
<mrow>
<msubsup>
<mi>M</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
当r=m;j=n±1;j+r=even;n+m=odd;
<mrow>
<msubsup>
<mi>M</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
为中间变量;
<mrow>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mi>r</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mi>r</mi>
<mo>)</mo>
<mo>!</mo>
<mo>!</mo>
</mrow>
</mfrac>
</mrow>
(A)!!代表数A的双阶乘,有如下定义
(A)!!=A(A-2)(A-4)…2,A=even
(A)!!=A(A-2)(A-4)…1,A=odd
0!!=1;(-1)!!=1;(-2)!!=∞;(-3)!!=-1;
(A)!!表示(j+r-1)!!、(j-r-1)!!、(j+r)!!、(j-r)!!;
当r=m;j+r=even;n+m=even;
<mrow>
<msubsup>
<mi>M</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>8</mn>
<mrow>
<msup>
<mi>&pi;</mi>
<mn>2</mn>
</msup>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>m</mi>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>-</mo>
<mn>2</mn>
<mi>m</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
当r≠m;
<mrow>
<msubsup>
<mi>M</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mn>0</mn>
</mrow>
其中[Mc]和[Ms]矩阵的差别在于[Ms]矩阵中不含有r=0的所有行和列;[Dc]和[Ds]矩阵中各元素的计算公式为:
当r=m;j+r=odd;n+m=odd;
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>K</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
</mfrac>
</mrow>
当r=m;j+r=even;n+m=even;
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>K</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
</mfrac>
</mrow>
当r=m;j+r=odd;n+m=even;
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mi>&pi;</mi>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>m</mi>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mfrac>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>3</mn>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
当r=m;j+r=even;n+m=odd;
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mi>&pi;</mi>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>m</mi>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mfrac>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>3</mn>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
当r≠m;
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mn>0</mn>
</mrow>
其中为中间变量,
<mrow>
<msubsup>
<mi>K</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
</mrow>
</msup>
</msup>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>,</mo>
</mrow>
其中[Dc]和[Ds]矩阵的差别在于[Ds]矩阵中不含有r=0的所有行和列;[Lc]矩阵和[Ls]矩阵中各元素的计算公式为:
且
<mrow>
<mo>(</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>)</mo>
<mo>=</mo>
<mo>(</mo>
<msup>
<mi>X</mi>
<mrow>
<mo>|</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mi>r</mi>
</mrow>
<mo>|</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msup>
<mi>X</mi>
<mrow>
<mo>|</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>r</mi>
</mrow>
<mo>|</mo>
</mrow>
</msup>
<mo>)</mo>
<mo>(</mo>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>)</mo>
<mo>=</mo>
<mo>(</mo>
<msup>
<mi>X</mi>
<mrow>
<mo>|</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mi>r</mi>
</mrow>
<mo>|</mo>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msup>
<mi>X</mi>
<mrow>
<mo>|</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>r</mi>
</mrow>
<mo>|</mo>
</mrow>
</msup>
<mo>)</mo>
<mo>(</mo>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
式中Xm、X|m-r|、X|m+r|、为中间变量;
其中min(r,m)代表比较r、m,取两者中较小的整数;
当r+m=odd;j+r=odd;n+m=odd;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mi>K</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>K</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
当r+m=odd;j+r=even;n+m=even;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mi>K</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>K</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
<mrow>
<msubsup>
<mi>K</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mi>r</mi>
</mrow>
</msup>
</msup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
当r+m=even;j+r=odd;n+m=odd;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>-</mo>
<mn>2</mn>
<mi>r</mi>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
当r+m=even;j+r=even;n+m=even;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>8</mn>
<mrow>
<msup>
<mi>&pi;</mi>
<mn>2</mn>
</msup>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>-</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
当r+m=odd;j+r=odd;n+m=even;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>4</mn>
<mrow>
<mi>&pi;</mi>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mn>3</mn>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mi>r</mi>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
当r+m=odd;j+r=even;n+m=odd;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>4</mn>
<mrow>
<mi>&pi;</mi>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mn>3</mn>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mi>r</mi>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
当r+m=even;j+r=odd;n+m=even;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
当r+m=even;j+r=even;n+m=odd;
<mrow>
<msubsup>
<mi>&Gamma;</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>r</mi>
</msubsup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
其中为中间变量。
3.根据权利要求2所述基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,其特征在于:所述最大谐波参数N取值范围为8-15。
4.根据权利要求3所述基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,其特征在于:所述步骤三中根据求解具体过程为:
在中取奇数量的Peters-Morillo模型展开系数和将 转换为此处将由下式计算获得
<mrow>
<msub>
<mrow>
<mo>{</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>}</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msubsup>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>f</mi>
</mrow>
<mi>m</mi>
</msubsup>
<mo>&rsqb;</mo>
<msub>
<mrow>
<mo>{</mo>
<msubsup>
<mi>a</mi>
<mi>f</mi>
<mi>m</mi>
</msubsup>
<mo>}</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mrow>
<mo>{</mo>
<msubsup>
<mi>&beta;</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>}</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msubsup>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>j</mi>
</mrow>
<mi>m</mi>
</msubsup>
<mo>&rsqb;</mo>
<msub>
<mrow>
<mo>{</mo>
<msubsup>
<mi>b</mi>
<mi>f</mi>
<mi>m</mi>
</msubsup>
<mo>}</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,f≤N,为混合模型中的奇数项系数,为由Peters-Morillo模型展开系数中元素索引的奇数项组成;
其中中的元素由以下公式计算获得
<mrow>
<msubsup>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>f</mi>
</mrow>
<mi>m</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<msqrt>
<mrow>
<msubsup>
<mi>H</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>f</mi>
<mi>m</mi>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>f</mi>
<mo>-</mo>
<mi>m</mi>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<msqrt>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>f</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>f</mi>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
5.根据权利要求4所述基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,其特征在于:所述步骤四中根据和步骤一获取的旋翼平面上方测试点集在椭圆坐标系ν,η,下的坐标构建混合模型,基于混合模型求解旋翼平面上方测试点处诱导流场的轴向诱导速度具体过程为:
利用如公式(9)的混合模型计算旋翼平面上方测试点处诱导流场的轴向速度:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</msub>
<mo>=</mo>
<munderover>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</munder>
<mrow>
<mi>m</mi>
<mo>&le;</mo>
<mi>n</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<msubsup>
<mi>&alpha;</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mfrac>
<mrow>
<msubsup>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>v</mi>
</mfrac>
<msubsup>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mover>
<mi>&psi;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</munder>
<mrow>
<mi>m</mi>
<mo>&le;</mo>
<mi>n</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<msubsup>
<mi>a</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mover>
<mi>&psi;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<munderover>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>o</mi>
<mi>d</mi>
<mi>d</mi>
</mrow>
</munder>
<mrow>
<mi>m</mi>
<mo>&le;</mo>
<mi>n</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<msubsup>
<mi>&beta;</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mfrac>
<mrow>
<msubsup>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>v</mi>
</mfrac>
<msubsup>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mover>
<mi>&psi;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
<mo>=</mo>
<mi>e</mi>
<mi>v</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</munder>
<mrow>
<mi>m</mi>
<mo>&le;</mo>
<mi>n</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<msubsup>
<mi>b</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mover>
<mi>&psi;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中表示基于混合模型计算得到的旋翼平面上方测试点处诱导流场的轴向速度,表示第一类归一化勒让德函数,表示第二类归一化勒让德函数,i为复数的虚部。
6.根据权利要求5所述基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,其特征在于:所述步骤五中根据混合模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度和Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度建立重构模型,根据重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度vz;具体过程为:
<mrow>
<msub>
<mi>v</mi>
<mi>z</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>a&eta;</mi>
<mi>c</mi>
</msup>
</mrow>
</mfrac>
<msub>
<mi>v</mi>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>a&eta;</mi>
<mi>c</mi>
</msup>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>a&eta;</mi>
<mi>c</mi>
</msup>
</mrow>
</mfrac>
<msub>
<mi>v</mi>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</msub>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中vz代表旋翼平面上方测试点处诱导流场的轴向诱导速度,表示由Peters-Morillo模型计算得到的旋翼平面上方测试点处诱导流场的轴向诱导速度,a,c为与重构模型收敛速度相关的系数;
<mrow>
<msub>
<mi>v</mi>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mi>m</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<msubsup>
<mi>a</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mover>
<mi>&psi;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mi>m</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<msubsup>
<mi>b</mi>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mover>
<mi>&psi;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
7.根据权利要求6所述基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法,其特征在于:所述步骤六中利用基因算法、蚂蚁算法或粒子算法对重构模型中的系数进行优化,使重构模型误差达到最小,最终求出优化后重构模型,根据优化后重构模型计算旋翼平面上方测试点处诱导流场的轴向诱导速度vz′;具体过程为:
利用卷积公式计算旋翼平面上方测试点中参照点处在给定压力系数时的轴向诱导速度的准确理论值
利用重构模型计算旋翼平面上方测试点中参照点处诱导流场的轴向诱导速度,并利用如下公式计算重构模型误差:
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>e</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>P</mi>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>v</mi>
<mi>z</mi>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mi>z</mi>
<mi>k</mi>
</msubsup>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中||e||2代表旋翼平面上方参照点的误差二范数的平方之和,为由重构模型计算得到的旋翼平面上方参照点k处的轴向诱导速度,为由卷积公式计算得到的旋翼平面上方参照点k处的轴向诱导速度的准确理论值;
利用基因算法、蚂蚁算法或粒子算法,对下式进行参数优化
min||e||2subject to a>0,0<c≤1, (13)
subject to a>0,0<c≤1表示使a、c服从a>0,0<c≤1;
使得a,c的选取在旋翼平面上方最优,得到优化后重构模型(10);
利用优化后重构模型(10)计算旋翼平面上方其他测试点处的诱导流场的轴向诱导速度v′z。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710744273.0A CN107688705A (zh) | 2017-08-25 | 2017-08-25 | 基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710744273.0A CN107688705A (zh) | 2017-08-25 | 2017-08-25 | 基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107688705A true CN107688705A (zh) | 2018-02-13 |
Family
ID=61155497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710744273.0A Pending CN107688705A (zh) | 2017-08-25 | 2017-08-25 | 基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107688705A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110162933A (zh) * | 2019-06-26 | 2019-08-23 | 上海交通大学 | 一种共轴多旋翼仿真方法及系统 |
CN111079235A (zh) * | 2019-12-11 | 2020-04-28 | 内蒙动力机械研究所 | 一种固体火箭发动机内流场仿真快速收敛方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130197871A1 (en) * | 2012-01-27 | 2013-08-01 | Fujitsu Limited | Predicted value calculation method and design support apparatus |
CN105117517A (zh) * | 2015-07-28 | 2015-12-02 | 中国电力科学研究院 | 一种基于改进粒子群算法的配网重构方法 |
US20160070833A1 (en) * | 2014-09-04 | 2016-03-10 | Synopsys, Inc. | Exponentially fitted approximation for anisotropic semiconductor equations |
CN106777739A (zh) * | 2016-12-28 | 2017-05-31 | 南京航空航天大学 | 一种倾转旋翼机倾转过渡过程的求解方法 |
-
2017
- 2017-08-25 CN CN201710744273.0A patent/CN107688705A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130197871A1 (en) * | 2012-01-27 | 2013-08-01 | Fujitsu Limited | Predicted value calculation method and design support apparatus |
US20160070833A1 (en) * | 2014-09-04 | 2016-03-10 | Synopsys, Inc. | Exponentially fitted approximation for anisotropic semiconductor equations |
CN105117517A (zh) * | 2015-07-28 | 2015-12-02 | 中国电力科学研究院 | 一种基于改进粒子群算法的配网重构方法 |
CN106777739A (zh) * | 2016-12-28 | 2017-05-31 | 南京航空航天大学 | 一种倾转旋翼机倾转过渡过程的求解方法 |
Non-Patent Citations (3)
Title |
---|
JIANZHE HUANG 等: "Converged Velocity Field for Rotors by a Blended Potential Flow Method", 《70TH AMERICAN HELICOPTER SOCIETY INTERNATIONAL ANNUAL FORUM 2014》 * |
JORGE A.MORILLO 等: "Velocity field above a rotor disk by a new dynamic inflow model", 《JOURNAL OF AIRCRAFT》 * |
吴文祥: "多自由度串联机器人关节摩擦分析与低速高精度运动控制", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110162933A (zh) * | 2019-06-26 | 2019-08-23 | 上海交通大学 | 一种共轴多旋翼仿真方法及系统 |
CN111079235A (zh) * | 2019-12-11 | 2020-04-28 | 内蒙动力机械研究所 | 一种固体火箭发动机内流场仿真快速收敛方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103530482B (zh) | 一种非均匀入流中螺旋桨噪声数值预报方法 | |
CN106640546A (zh) | 一种对风力发电设备的塔筒进行监测的系统和方法 | |
Wang et al. | Typhoon-induced vibration response and the working mechanism of large wind turbine considering multi-stage effects | |
CN106777783B (zh) | 一种航空发动机叶片裂纹预测方法 | |
CN106286151B (zh) | 一种风电机组低速轴扭转载荷监测方法及载荷分析方法 | |
Firdaus et al. | Numerical and experimental studies of a small vertical-axis wind turbine with variable-pitch straight blades | |
Ghadirian et al. | Considering induction factor using BEM method in wind farm layout optimization | |
CN107688705A (zh) | 基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法 | |
Anagnostopoulos et al. | Optimal design and experimental validation of a Turgo model Hydro turbine | |
Miller et al. | Solidity effects on the performance of vertical-axis wind turbines | |
Bottasso et al. | Wind tunnel testing of wind turbines and farms | |
Peng et al. | High-solidity straight-bladed vertical axis wind turbine: Numerical simulation and validation | |
Rogowski et al. | Numerical analysis of a small-size vertical-axis wind turbine performance and averaged flow parameters around the rotor | |
CN205779882U (zh) | 一种叶片泵及用于叶片泵的叶轮叶片 | |
Kelly et al. | Predicting blade vortex interaction, airloads and acoustics using the vorticity transport model | |
CN109143893B (zh) | 一种基于真机实测的水轮机特性获取方法 | |
CN106503375B (zh) | 一种基于cn群理论确定汽轮机转子临界转速的方法及系统 | |
CN115577493A (zh) | 一种水平轴潮流能水轮机叶片非定常特性预测方法 | |
CN110046420B (zh) | 一种用于确定斜轴泵在不同叶片角度下飞逸转速的方法 | |
Ivanell et al. | Validation of methods using EllipSys3D | |
Rocchio et al. | Development of a BEM-CFD tool for Vertical Axis Wind Turbines based on the Actuator Disk model | |
He et al. | Analysis of a propeller wake flow field using viscous fluid mechanics | |
CN106050730A (zh) | 一种叶片泵及用于叶片泵的叶轮叶片 | |
CN204116020U (zh) | 非对称双臂离心加载装置 | |
Lee et al. | Performance of a direct drive hydro turbine for wave power generation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180213 |