一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法专利检索-升阻比空气动力学专利检索查询-专利查询网

一种基于固定翼鸭 舵的双旋稳定弹的 气动外形设计方法

阅读：845发布：2020-05-08

专利汇可以提供一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法专利检索，专利查询，专利分析的服务。并且本发明公开了一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法，通过分别计算固定翼鸭舵的舵面积SC、俯仰舵角设定值δZ和压心到弹体质心的距离XC后，根据其形成的约束条件求取升阻比 KL/D的最大值；再根据固定翼鸭舵内部与外部的摩擦力矩、电磁阻矩和气动力矩之间的力矩平衡关系，确定滚转舵角设定值δy和尾翼斜切角；在所述滚转舵角设定值δy和所述尾翼斜切角确定的情况下,调节所述升阻比KL/D，实现迫击炮的姿态调整。，下面是一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法专利的具体信息内容。

权利要求

1.一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法，其特征在于，分别计算固定翼鸭舵的舵面积SC、俯仰舵角设定值δZ和压心到弹体质心的距离XC后，根据其形成的约束条件求取升阻比KL/D的最大值；再根据固定翼鸭舵内部与外部的摩擦力矩、电磁阻矩和气动力矩之间的力矩平衡关系，确定滚转舵角设定值δy和尾翼斜切角；在所述滚转舵角设定值δy和所述尾翼斜切角确定的情况下,调节所述升阻比KL/D，实现迫击炮的姿态调整；
其中，根据固定翼鸭舵的舵面积SC、俯仰舵角设定值δZ和压心到弹体质心的距离XC和形成的约束条件求取升阻比KL/D的最大值的计算公式如下：
其中：
ScCNδ*(δz+α)xc≥-SCmαsinα；
ScCNδ*(δz+αb)xc≤-SCmαsinαb；
0.8≤δz/αb≤1.2；
xc(min)0≤α≤αb；0<δz；
式中：αb为平衡攻角；CLα表示气动升力系数；α表示攻角；CD表示气动阻力系数；Cmα表示气动静力矩系数；xc表示固定翼鸭舵压心到弹体质心的距离；Sc表示固定翼鸭舵的特征面积；CNδ表示固定翼鸭舵气动控制力系数；δz表示固定翼鸭舵的俯仰舵角；S表示弹体特征截面积。
2.如权利要求1所述的气动外形设计方法，其特征在于，计算所述固定翼鸭舵的舵面积SC的方法包括：
计算固定翼鸭舵的展弦比λC：
式中，dc为固定翼鸭舵直径；bav为平均翼弦；
计算固定翼鸭舵面积SC：
式中，df为控制舱直径。
3.如权利要求2所述的气动外形设计方法，其特征在于，所述俯仰舵角设定值δZ和压心到弹体质心的距离XC通过复形调优求目标函数的最大值计算方法确定。
4.如权利要求1所述的气动外形设计方法，其特征在于，计算所述滚转舵角设定值δy的方法如下：
固定翼鸭舵气动滚转力矩系数Cclα：
式中：为固定翼鸭舵气动滚转力矩向量；V为炮弹相对大气的速度；ρ为地面空气密度；
将所述固定翼鸭舵气动滚转力矩系数Cclα通过气动系数对应表，用插值法计算得出所述滚转舵角设定值δy。
5.如权利要求4所述的气动外形设计方法，其特征在于，计算所述尾翼斜切角的方法包括：
根据尾翼导转力矩模值大于外部极阻尼力矩模值与电机内部的电磁力矩Me之和的公式，计算尾翼导转系数Crα的对应值：
式中：Crα为尾翼导转系数；Clp为极阻尼力矩系数；pA为弹体转速；Me为电磁力矩；
将尾翼导转系数Crα的对应值通过气动系数对应表，用插值法计算得出所述尾翼斜切角。

说明书全文

一种基于固定翼鸭舵的双旋稳定弹的 气动外形设计方法

技术领域

[0001] 本发明涉及航天领域，尤其涉及一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法。

背景技术

[0002] 迫击炮具有弹道弯曲、射速快、结构简单、便于机动、易于操作等特点，适于城市巷战和山地作战，在世界各国得到了广泛的应用。对迫击炮弹的精确问题，现有技术中，其主要利用十字鸭舵实现炮弹的精确制导，此种方式成本较高，且现有技术中主要围绕十字鸭舵的气动外形参数设计，其并不适用于固定翼鸭舵的气动外形设计。

发明内容

[0003] 针对上述现有技术存在的缺陷，本发明提供一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法，合理利用炮弹外部空气动力，实现了对迫击炮的姿态调整。

[0004] 本发明提供的一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法，其改进之处在于，分别计算固定翼鸭舵的舵面积SC、俯仰舵角设定值δZ和压心到弹体质心的距离XC后，根据其形成的约束条件求取升阻比KL/D的最大值；再根据固定翼鸭舵内部与外部的摩擦力矩、电磁阻矩和气动力矩之间的力矩平衡关系，确定滚转舵角设定值δy和尾翼斜切角；在所述滚转舵角设定值δy和所述尾翼斜切角确定的情况下,调节所述升阻比KL/D，实现迫击炮的姿态调整。

[0005] 优选的，计算所述固定翼鸭舵的舵面积SC的方法包括：

[0006] 计算固定翼鸭舵的展弦比λC：

[0007]

[0008] 式中，dc为固定翼鸭舵直径；bav为平均翼弦；

[0009] 计算固定翼鸭舵面积SC：

[0010]

[0011] 式中，df为控制舱直径。

[0012] 较优选的，所述俯仰舵角设定值δZ和压心到弹体质心的距离XC通过复形调优求目标函数的最大值计算方法确定。

[0013] 较优选的，所述根据固定翼鸭舵的舵面积SC、俯仰舵角设定值δZ和压心到弹体质心的距离XC形成的约束条件求取升阻比KL/D的最大值，计算公式如下：

[0014]

[0015] 其中：

[0016] ScCNδ*(δz+α)xc≥-SCmαsinα；

[0017] ScCNδ*(δz+αb)xc≤-SCmαsinαb；

[0018] 0.8≤δz/αb≤1.2；

[0019] xc(min)

[0020] 0≤α≤αb；0<δz；

[0021] 式中：αb为平衡攻角；CLα表示气动升力系数；α表示攻角；CD表示气动阻力系数；Cmα表示气动静力矩系数；xc表示固定翼鸭舵压心到弹体质心的距离；Sc表示固定翼鸭舵的特征面积；CNδ表示固定翼鸭舵气动控制力系数；δz表示固定翼鸭舵的俯仰舵角；S表示弹体特征截面积。

[0022] 较优选的，计算所述滚转舵角设定值δy的方法如下：

[0023] 固定翼鸭舵气动滚转力矩系数Cclα：

[0024]

[0025] 式中：为固定翼鸭舵气动滚转力矩向量；V为炮弹相对大气的速度；ρ为地面空气密度；

[0026] 将所述固定翼鸭舵气动滚转力矩系数Cclα通过气动系数对应表，用插值法计算得出所述滚转舵角设定值δy。

[0027] 较优选的，计算所述尾翼斜切角的方法包括：

[0028] 根据尾翼导转力矩模值大于外部极阻尼力矩模值与电机内部的电磁力矩Me之和的公式，计算尾翼导转系数Crα的对应值：

[0029]

[0030] 式中：Crα为尾翼导转系数；Clp为极阻尼力矩系数；pA为弹体转速；Me为电磁力矩；

[0031] 将尾翼导转系数Crα的对应值通过气动系数对应表，用插值法计算得出所述尾翼斜切角。

[0032] 本发明通过设计固定翼鸭舵和斜切尾翼气动参数，能够合理利用炮弹外部空气动力，利用外部空气动力作为固定翼鸭舵的能量来源，类似于发电机的发电原理。同时由于采用固定舵角的固定翼鸭舵，如果固定舵角的角度较大，则制导迫击炮弹的静稳定性较差，如果固定舵角的角度较小，则固定翼鸭舵的弹道修正能力不足，本发明利用全弹升阻比建立相应气动外形参数的约束条件求取气动外形参数的最优值，此发明的计算方法即可保证迫击炮的稳定性又可实现迫击炮的姿态调整，对迫击炮的弹道进行了修正，从而提高了迫击炮的落地精准度。附图说明

[0033] 图1为本发明实施例的流程图；

[0034] 图2为本发明实施例的双旋稳定弹外形设计示意图；

[0035] 图3为本发明实施例的固定翼鸭舵的参数示意图。

具体实施方式

[0036] 为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举出优选实施例，对本发明进一步详细说明。然而，需要说明的是，说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解，即便没有这些特定的细节也可以实现本发明的这些方面。

[0037] 本实施例提供的一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法，双旋稳定弹外形的示意图如图2所示，其流程图如图1所示，分别计算固定翼鸭舵的舵面积SC、俯仰舵角设定值δZ和压心到弹体质心的距离XC后，根据其形成的约束条件求取升阻比KL/D的最大值；再根据固定翼鸭舵内部与外部的摩擦力矩、电磁阻矩和气动力矩之间的力矩平衡关系，确定滚转舵角设定值δy和尾翼斜切角；当滚转舵角设定值δy和尾翼斜切角确定时，即迫击炮的能量确定时，再通过升阻比和稳定性求取的相应气动外形参数，在保证炮弹静穩定性的基础下使气动升力值最大，从而通过气动升力调整弹体姿态，实现弹道修正。

[0038] 本实施例的各参数的求取方法如下：

[0039] 1.求取头部长度Ln；

[0040] 其中头部长细比λn的求取公式为：

[0041]

[0042] 式中：d为弹体特征直径。

[0043] 头部长细比对头部波阻影响较大，λn越大阻力越小，但是为了保证双旋弹二维修正控制舱的小型化标准，本实施例选取头部长细比λn为2，进行二维修正控制舱长度的设计。由头部长细比λn公式可得，当弹径d为120mm，头部长度Ln为240mm。

[0044] 2.固定翼鸭舵的示意图如图3所示，计算固定翼鸭舵的舵面积SC的方法包括：

[0045] 计算固定翼鸭舵的展弦比λC，主要由固定翼鸭舵的综合影响参考系数决定，其公式为：

[0046]

[0047] 式中，dc为固定翼鸭舵直径；bav为平均翼弦；

[0048] 本实施例基于120mm迫击炮弹改装的双旋弹在亚声速下旋转飞行，根据旋转弹翼式展弦比和亚声速飞行器展弦比的范围，选取固定翼鸭舵的展弦比的合理值，从而求取固定翼鸭舵的舵面积Sc。计算固定翼鸭舵面积SC公式为：

[0049]

[0050] 式中，df为控制舱直径。

[0051] 根据炮弹外形参考设计中提供的旋转弹翼式展弦比λC范围为2～4，亚声速飞行器的展弦比λC为4～6，由于基于120mm迫击炮弹改装的双旋弹在亚声速下旋转飞行，因此选取固定翼鸭舵的展弦比λC为4，本实施例根据迫击炮弹的出膛发射方式，为便于发射，设计翼展dc与弹径d相同，由展弦比公式λC可以得到平均翼弦bav为0.03m。

[0052] 根据120mm迫击炮弹头部控制舱安装螺纹的尺寸和控制舱内部控制机构尺寸，设计控制舱直径df为弹径d的1/2，由翼展dc和平均翼弦bav，设计舵片的形状为矩形可得固定翼鸭舵的固定翼鸭舵面积SC为0.9*10-3m2。

[0053] 3.本实施例俯仰舵角设定值δZ和压心到弹体质心的距离XC可通过复形调优求目标函数的最大值计算方法求取。

[0054] 4.根据固定翼鸭舵的舵面积SC、俯仰舵角设定值δZ和压心到弹体质心的距离XC形成的约束条件求取升阻比KL/D的最大值，计算公式如下：

[0055]

[0056] 其中，根据约束条件下n维极值的复形调优算法，列写约束条件：

[0057] 最大值：操纵性KL/D

[0058] 约束条件：机动性要求：ScCNδ*(δz+α)xc≥-SCmαsinα

[0059] 静稳定性要求：ScCNδ*(δz+αb)xc≤-SCmαsinαb

[0060] 操稳比：0.8≤δz/αb≤1.2

[0061] 固定翼鸭舵安装位置：xc(min)

[0062] 系数设定：0≤α≤αb；0<δz；

[0063] 式中：αb为平衡攻角；CLα表示气动升力系数；α表示攻角；CD表示气动阻力系数；Cmα表示气动静力矩系数；xc表示固定翼鸭舵压心到弹体质心的距离；Sc表示固定翼鸭舵的特征面积；CNδ表示固定翼鸭舵气动控制力系数；δz表示固定翼鸭舵的俯仰舵角；S表示弹体特征截面积。

[0064] 约束条件中求取升阻比KL/D的最大值，从而使气动力合成矢量方向与气动升力矢量之间的夹角最小，保证整弹的操纵性最优。为了保证固定翼鸭舵的机动性，需要保证固定翼鸭舵的控制力矩模值大于等于迫击炮弹的静力矩模值同时为了保证整弹的静稳定性，需要保证攻角α大于平衡攻角αb时，静力矩模值大于等于控制力矩模值根据固定翼鸭舵的操稳比δz/αb范围，确定固定翼鸭舵固定翼舵角δz与平衡攻角αb之间的约束范围。根据弹头距质心位置和控制舱安装螺纹距质心的位置，确定固定翼鸭舵压心位置距质心位置xc的安装位置范围。

[0065] 本实施例设定目标函数J＝-KL/D，采用复形调优法求目标函数J的极小值点的迭代过程如下：

[0066] 复形共有2n个顶点。假设给定初始复形中的第一个顶点坐标：

[0067] X(0)＝(x00,x10,…,xn-1,0)

[0068] 且此顶点坐标满足n个常数约束条件及m个函数约束条件。

[0069] 1)在n维变量空间中在确定出初始复形的其余2n-1个顶点。其方法如下：

[0070] 利用伪随机数按常量约束条件产生第j个顶点X(j)＝(x0j,x1j,…,xn-1,j)(j＝1,2,…,2n-1)中的各分量xij(i＝0,1,…,n-1)，即：

[0071] xij＝ai+r(bi-ai)

[0072] 式中：ai和bi为约束条件ai≤xi≤bi；r为[0,1]之间的一个伪随机数。

[0073] 显然，由上述方法产生的初始复形的各顶点满足常量约束条件。然后再检查它们是否符合函数约束条件，如果不符合，则需要作调整，直到全部顶点均符合函数约束条件及常量约束条件为止。调整的原则为：

[0074] 假设前j个顶点已满足所有的约束条件，而第j+1个顶点不满足约束条件，则作如下调整变换(j＝1，2，…，2n-1)：

[0075]

[0076] 式中这个过程一直运行到满足所有约束条件为止。

[0077] 初始复形的2n个顶点确定以后，计算各顶点处的目标函数值：

[0078] f(j)＝f(X(j)),j＝0,1,…,2n-1

[0079] 2)确定：

[0080]

[0081]

[0082] 其中X(R)称为最坏点。

[0083] 3)计算最坏点X(R)的对称点；

[0084] XT＝(1+a)XF-aX(R)

[0085] 其中

[0086]

[0087] a称为反射系数，一般取1.3左右。

[0088] 4)确定一个新的顶点替代最坏点X(R)以构成新的复形。其方法如下：

[0089] 如果f(XT)>f(G)为止，则用下式修改XT：

[0090] XT＝(XT+XF)/2

[0091] 直到f(XT)≤f(G)为止。

[0092] 然后检查XT是否满足所有约束条件。如果对于某个分量XT(j)不满足常量约束条件，即如果XT(j)bj，则令

[0093] XT(j)＝aj+δ或XT(j)＝bj-δ

[0094] 其中δ为很小的一个常数，一般取δ＝10-6。然后重复第四个步骤。直到f(XT)≤f(G)且XT满足所有约束条件为止。此时令：

[0095] X(R)＝XT,f(R)＝f(XT)

[0096] 重复2)～4)，直到复形中各顶点距离小于预先给定的精度要求为止。根据约束条件下n维极值的复形调优算法计算升阻比KL/D为最优解2.28时，对应的设计参数固定翼鸭舵压心到弹体质心的距离xc为0.3338m、平衡攻角αb为12°和固定翼舵偏角δz为10°。

[0097] 5.计算滚转舵角设定值δy的方法如下：

[0098] 固定翼鸭舵二维修正控制舱利用内部永磁电机的电磁力矩作为控制力矩，为了降低电子负载内部MOS管的功耗，设定内部永磁电机的电磁电枢电流的适用性范围，为使固定翼鸭舵相对弹体固定坐标系中的滚转角保持不变，需要内部电磁力矩最大值与外部固定翼鸭舵气动滚转力矩最大值保持一致，利用公式可以得到Cclα的对应值，从而确定固定翼鸭舵的滚转舵角δy设计值。

[0099]

[0100] 式中：为固定翼鸭舵气动滚转力矩向量；Cclα表示固定翼鸭舵气动滚转力矩系数；V为炮弹相对大气的速度；ρ为地面空气密度。

[0101] 将固定翼鸭舵气动滚转力矩系数Cclα通过气动系数对应表，用插值法计算得出所述滚转舵角设定值δy。

[0102] 本实施例设定内部永磁电机电磁电枢电流的适用性范围为0～1A，由于电磁力矩常数KM为0.1705Nm/A，因此电磁力矩最大值近似为0.1705Nm，当炮弹相对大气的速度V为300m/s，地面空气密度ρ值为1.2063kg/m3，为使固定翼鸭舵相对弹体固定坐标系中的滚转角保持不变，需要内部电磁力矩最大值与外部固定翼鸭舵气动滚转力矩最大值保持一致，因此滚转力矩模值最大值近似为0.1705Nm，利用固定翼鸭舵气动滚转力矩向量公式近似得到Cclα的对应值，从而确定固定翼鸭舵的滚转舵角δy设计值为5°。

[0103] 6.计算尾翼斜切角的方法包括：

[0104] 根据尾翼导转力矩模值需要大于外部极阻尼力矩模值与电机内部的电磁力矩Me之和的公式，计算尾翼导转系数Crα的对应值：

[0105]

[0106] 式中：Crα为尾翼导转系数；Clp为极阻尼力矩系数；pA为弹体转速；Me为电磁力矩；

[0107] 将尾翼导转系数Crα的对应值通过气动系数对应表，用插值法计算得出所述尾翼斜切角。本实施例电磁力矩Me的范围为0～0.1705Nm，弹体转速pA范围为0～1200rpm，根据弹体转动的动力学关系公式可以计算尾翼导转系数Crα在炮弹速度V为300m/s时的对应值，确定尾翼斜切角设计值为15°。

[0108] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

标题	发布/更新时间	阅读量
一种短航程低过载的再入轨迹设计方法	2020-05-14	61
运载器能量管理段FADS故障判断与控制方法、运载器及存储介质	2020-05-08	704
垂直轴风力机叶片、垂直轴风轮及垂直轴风力机	2020-05-08	808
强约束复杂任务条件下的三维剖面优化设计方法、系统及介质	2020-05-14	743
涡流发生器和叶片	2020-05-12	751
可回收探空仪、系统及回收方法	2020-05-13	63
一种分离式减摇鳍	2020-05-14	267
一种飞行器机翼	2020-05-15	267
一种双翼构型的新概念布局飞机	2020-05-11	88
一种螺旋桨飞机机翼及螺旋桨飞机	2020-05-14	166

一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法

一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：