一种船尾外形旋转弹马格努斯力滚转导数获取方法，主要步骤包括：根据旋转弹几何特征参数和来流条件计算得到等效长度雷诺数、马赫数和几何外形影响因子，最后代入解析估算公式得到全弹的马格努斯滚转导数。本发明专利技术针对船尾外形的旋转弹，采用等效长度思想，利用解析估算公式快速高效的得到全弹的马格努斯力滚转导数，相较于非定常计算方法，本发明专利技术大大节省了计算资源和计算机时。
A method to obtain the Magnus force roll derivative of the stern shape rotating bomb
【技术实现步骤摘要】
一种船尾外形旋转弹马格努斯力滚转导数获取方法
本专利技术涉及一种船尾外形旋转弹马格努斯力滚转导数获取方法，属于旋转弹气动设计
，用于榴弹或火箭炮发射的炮弹的气动外形设计和气动性能预测。
技术介绍
炮弹、火箭弹、导弹等战术武器为提高其飞行稳定性，通常采用绕体轴旋转的的方式飞行。但是当自转和攻角同时存在时，由于流场边界层的畸变和离心力因素，弹体附近流场相对于攻角平面不对称，因此会产生垂直于攻角平面的力，通常称为马格努斯力、也叫面外力。马格努斯力会影响弹丸的航向动稳定性，降低打靶精度。因此，准确预测旋转弹的马格努斯特性成为弹道计算、弹箭设计和稳定性研究的必然要求。马格努斯特性是由于弹体纵-横向耦合运动而产生的一类复杂的气动现象，其诱发的因素较多，如非对称湍流剪切层、非对称旋涡、非对称转捩等。针对这类复杂的气动问题，国内外采用的研究手段主要有两种，一是风洞试验，二是数值计算。风洞试验具有代表性的工作包括：①Leroy在ArnoldEngineeringandDevelopmentCenter(AEDC)研究了马赫数0.2～2.5，攻角0～90°的舵身组合外形旋转气动特性，并发现在攻角大于20°、马赫数小于1.5时滚转阻尼和马格努斯特性随旋转速率呈现非线性(LeroyMJ.ExperimentalRoll-Damping,Magnus,andStaticStabilityCharacteristicsofTwoSlenderMissileConfigurationsatHighAnglesofAttack(0to90Deg)andMachNumber0.2Through2.5[R].AEDC-TR-76-58,1976.)，②雷娟棉利用风洞试验研究了马赫数0.4～4.5，攻角-2～12°条件下正、反装卷弧形尾翼及平板型尾翼的面外力特性，认为卷弧形尾翼反装且反向旋转是抑制旋转火箭弹锥形失稳的有效措施(雷娟棉,吴甲生.尾翼稳定大长径比无控旋转火箭弹的锥形运动与抑制[J].空气动力学学报,2005,23(4):455-458)。对于数值计算，专利技术人所在团队做了大量的工作：采用完全时间相关的非定常Navier-Strokes方程，研究了马赫数2.5，攻角-5～40°翼身组合弹箭马格努斯特性及产生机理，发现切应力产生的马格努斯力只占压力产生马格努斯力的1％；研究了不同湍流模型模拟的旋转弹近壁湍流粘性系数差异，发现一方程湍流模型和两方程湍流模型的物面摩阻系数最大相差35％(①石磊,杨云军,周伟江.两种湍流模型在高速旋转翼身组合弹箭中的对比研究[J],力学学报,2017,49(1):84-92.②石磊,刘周,杨云军,周伟江.弯曲和旋转修正的湍流模型在旋转翼身组合弹箭中的应用研究[J].空气动力学学报,2018,36(4):620-625)。针对旋转弹开展的试验，主要为了获得指定状态下的全弹马格努斯特性，进而研究其随马赫数、转速和攻角的变化规律；数值计算集中在典型状态下的流场分析，进而研究马格努斯效应产生的机理。不管是试验还是数值计算，都需要花费大量的时间和费用。经过课题组大量的调研和研究，认为建立一种适用于旋转弹马格努斯特性的快速计算方法，能够快速高效的预估旋转弹的马格努斯特性，可为旋转弹设计和旋转弹气动特性评估提供有效支撑，同时也是风洞试验和数值计算的有力补充。而关于旋转弹气动特性方面的专利，只搜索到2013年闫晓勇利用试验数据进行气动参数辨识的专利技术(专利名称:一种旋转弹气动参数获取方法,专利号:CN201310143748.2,专利技术人:闫晓勇,钟凌伟)，与本专利技术相关的专利还未能查到。
技术实现思路
本专利技术的技术解决问题：克服现有技术的不足，提出一种旋转弹马格努斯力滚转导数快速计算方法，通过解析的经验公式预测旋转弹在飞行状态或风洞状态下的马格努斯特性。本专利技术的技术解决方案：一种船尾外形旋转弹马格努斯力滚转导数获取方法，1)获取旋转弹的几何尺寸，包括弹长L、表面积S、弹径D、弹头长度Lh、船尾长度LB、船尾锥角θ；2)获取飞行状态参数，包括单位长度雷诺数Re和马赫数Ma；3)计算得到等效长度雷诺数；4)使用马赫数按照线性公式计算得到马赫数影响因子ηM；5)根据步骤1)得到的旋转弹的几何尺寸，计算得到几何外形影响因子ηgeom；6)根据有效长度雷诺数Lval、马赫数影响因子ηM、几何外形影响因子ηgeom，计算得到马格努斯力滚转导数所述步骤3)计算得到等效长度雷诺数的具体过程为：用旋转弹表面积除以弹径与圆周率的乘积得到旋转弹的等效长度Lval；将单位长度雷诺数Re乘以等效长度Lval，得到等效长度雷诺数Reval。所述步骤4)马赫数影响因子ηM的计算公式如下：ηM＝1-0.0927Ma。所述步骤6)几何外形影响因子ηgeom的计算公式如下：所述步骤6)马格努斯力滚转导数计算的计算公式如下：本专利技术与现有技术相比的优点在于：(1)本专利技术获得解析估算公式可以体现马格努斯力滚转导数的非线性特征，与传统的线性估算公式相比，准确性更好；(2)本专利技术无需进行风洞试验和非定常计算，可以快速高效的估算出马格努斯力滚转导数，大大的节省费用和计算时间。附图说明：图1为本专利技术流程图。图2为船尾外形旋转弹几何外形图，图中L表示弹长、D表示弹径、Lh表示弹头长度、LB表示船尾长度、θ表示船尾锥角。具体实施方式：具体实施方式如下：本示例选取的外形、风洞条件和实验结果均来自于文献(LeroyM.ExperimentalMagnusCharacteristicsofBallisticProjectileswithanawithoutAnti-magusVanesatMachNumbers1.5through2.5[R].AD771807,1973)(1)获取旋转弹的几何尺寸，包括弹长L、表面积S、弹径D、弹头长度Lh、船尾长度LB、船尾锥角θ，表1给出三种外形的几何尺寸。表1几何尺寸表外形L(mm)D(mm)S(mm2)Lh(mm)LB(mm)θ(deg)1728139.62.4369E+5384.869.852686.8132.082.1153E+5363.2132.087.53853.44142.242.5928E+5504.85320.043.72(2)获取飞行(风洞)状态参数，单位长度雷诺数Re和马赫数Ma。表2给出3个飞行(风洞)状态的单位雷诺数和马赫数。表2飞行(风洞)状态参数本文档来自技高网...
【技术保护点】
1.一种船尾外形旋转弹马格努斯力滚转导数获取方法，其特征在于：/n1)获取旋转弹的几何尺寸，包括弹长L、表面积S、弹径D、弹头长度L
【技术特征摘要】
1.一种船尾外形旋转弹马格努斯力滚转导数获取方法，其特征在于：
1)获取旋转弹的几何尺寸，包括弹长L、表面积S、弹径D、弹头长度Lh、船尾长度LB、船尾锥角θ；
2)获取飞行状态参数，包括单位长度雷诺数Re和马赫数Ma；
3)计算得到等效长度雷诺数；
4)使用马赫数按照线性公式计算得到马赫数影响因子ηM；
5)根据步骤1)得到的旋转弹的几何尺寸，计算得到几何外形影响因子ηgeom；
6)根据有效长度雷诺数Lval、马赫数影响因子ηM、几何外形影响因子ηgeom，计算得到马格努斯力滚转导数
2.根据权利要求1所述的一种船尾外形旋转弹马格努斯力滚转导数获取方法，其特征在于：所述步骤3)计算得到等效长度雷诺数的具体...
【专利技术属性】
技术研发人员：石磊，豆国辉，王利，
申请(专利权)人：中国航天空气动力技术研究院，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。