一种非奇异多约束弹道快速优化方法技术

本发明专利技术公开了一种非奇异多约束弹道快速优化方法，首先根据已有的发射系下的超远程弹道动力学方程，得到考虑地球引力摄动，附加哥式惯性力、离心惯性力、高空风等长时作用力的影响；其次姿态矩阵将体坐标系和发射系变量转换到速度系下，得到速度系下垂直发射非奇异弹道方程；再次，设置控制变量输出为攻角和侧滑角的导数，设置状态量为发射系速度、弹道倾角、弹道偏角、攻角和侧滑角；最后，这些状态量交由伪谱法进行快速弹道优化，并对优化结果进行三次样条曲线拟合，得到适用于内环跟踪的标称弹道。该方法适用于超远程火箭完整建模的弹道优化，并适合内环弹道响应，可广泛应用于各类要求垂直发射和长时间飞行的弹道导弹和巡航导弹中。弹中。弹中。

【技术实现步骤摘要】
一种非奇异多约束弹道快速优化方法


[0001]本专利技术属于导弹
，具体涉及一种弹道快速优化方法。

技术介绍

[0002]弹道优化技术经过多年发展逐步趋于成熟，相比于经验化的弹道设计，弹道优化技术寻找到满足各种约束下且能实现给定指标最优的弹道，正在迈向工程化应用。已有的弹道优化技术通常采用速度坐标系进行微分约束的建模，但为了适应惯性导航装置的定位模型，超远程导弹的六自由度模型通常建立在发射坐标系，气动参数差值表通常建立在弹体坐标系，不能直接提供速度坐标下的变量计算；其次，传统的速度坐标系方程不考虑摄动、哥式惯性力、离心惯性力、高空风等长时作用；最后，传统的弹道优化结果由于优化变量的离散化会导致优化变量稀疏且容易出现跳变，只作为导弹能力验证，优化结果不能直接用于导弹内环响应进行跟踪。故如何开发一种兼顾导航装置模型，使其利于设置约束，对垂直发射不奇异，且能直接用于跟踪的弹道优化技术非常关键。

技术实现思路

[0003]为了克服现有技术的不足，本专利技术提供了一种非奇异多约束弹道快速优化方法，首先根据已有的发射系下的超远程弹道动力学方程，得到考虑地球引力摄动，附加哥式惯性力、离心惯性力、高空风等长时作用力的影响；其次利用化简以后的姿态矩阵将体坐标系和发射系变量转换到速度系下，得到速度系下垂直发射非奇异弹道方程；再次，为了得到平滑的有利于内环控制响应的曲线，设置控制变量输出为攻角和侧滑角的导数，设置状态量为发射系速度、弹道倾角、弹道偏角、攻角和侧滑角；最后，这些状态量交由伪谱法进行快速弹道优化，并对优化结果进行三次样条曲线拟合，得到适用于内环跟踪的标称弹道。该方法适用于超远程火箭完整建模的弹道优化，并适合内环弹道响应，可广泛应用于各类要求垂直发射和长时间飞行的弹道导弹和巡航导弹中。
[0004]本专利技术解决其技术问题所采用的技术方案包括如下步骤：
[0005]步骤1：建立发射系动力学方程：
[0006][0007]其中，表示导弹质量，表示导弹速度，、、分别为发射系下的
坐标；为气动力，P为推力；为地球引力；为离心惯性力；为哥式惯性力；
[0008]气动力和推力需要将体坐标系的分量转化为发射系，即：
[0009][0010]其中为体坐标系和发射系的转换矩阵，表示动压，表示参考面积，、、分别表示轴向力、法向力、侧向力气动系数；
[0011]地球引力项包含径向引力和切向引力：
[0012][0013]其中径向引力，切向引力，为地心距，为地理纬度，、、分别表示发射点在地心固连坐标系下的坐标，、、分别表示地球旋转角速度地心固连坐标系下的坐标；为地球引力常数，第二摄动常数，为地球长半轴；
[0014]离心惯性力为：
[0015][0016]其中为地球旋转角速度；
[0017]哥式惯性力：
[0018][0019]体坐标系到发射系的转换矩阵需要获得姿态角信息，滚转角、偏航角、俯仰角三个欧拉姿态角分别为：
[0020][0021]其中，表示侧滑角，表示弹道倾角，表示攻角，表示弹道倾角；
[0022]因此，体坐标系到发射系的转换矩阵为
[0023][0024]海拔高度采用地心高度差进行计算：
[0025][0026]其中弹道任一点的地心距，为发射点地球半径分量，为地球半径；
[0027]步骤2：风攻角、风侧滑、风合攻角计算；
[0028]风攻角、风侧滑、风合攻角由以下公式计算：
[0029][0030]其中，分别为考虑高空风影响下的弹体系在xyz三个方向的速度分量，通过发射系下速度分量转换而成：
[0031][0032]其中为发射系下受高空风影响的速度分量，通过风速、风向和发射系下速度合成：
[0033][0034]其中为发射系北向角，、分别为风速和风向；
[0035]步骤3：建立速度系下的微分约束；
[0036]设优化状态变量，控制变量为，得到优化问题微分约束为：
[0037][0038]步骤4：弹道优化；
[0039]优化约束包含高度限幅、热流密度限幅、速度限幅、弹道倾角限幅、弹道偏角限幅、攻角限幅、侧滑角限幅和速度限幅；优化指标包括最大落速、最大落角、最大射程、最小总加热量和最小攻角；设置伪谱法弹道优化微分约束和状态约束，采用伪谱法工具箱进行优化求解；
[0040]步骤5：对优化求解结果进行三次样条插值，对不同的状态优化变量输出结果进行拟合。
[0041]本专利技术的有益效果如下：
[0042]本专利技术方法通过深入分析远程导弹弹道特性，创新性的提出将发射系下的方程转化为速度坐标系下的坐标进行弹道优化，解决了垂直发射速度系弹道方程奇异的问题，一并解决了发射系下的约束过于复杂的问题。本专利技术方法设计的优化弹道在力热载荷非常严酷的情况下能够实现垂发导弹的优化弹道的直接跟踪响应，具有非常广阔的应用前景。
附图说明
[0043]图1为本专利技术实施例射程和弹道高度曲线。
[0044]图2为本专利技术实施例时间对应飞行速度曲线。
[0045]图3为本专利技术实施例时间对应法向过载曲线。
[0046]图4为本专利技术实施例时间对应攻角曲线。
[0047]图5为本专利技术实施例时间对应驻点热流密度曲线。
[0048]图6为本专利技术实施例时间对应倾角曲线。
具体实施方式
[0049]下面结合附图和实施例对本专利技术进一步说明。
[0050]为解决现有弹道优化技术在远程弹道上精度低，垂直发射奇异，不能直接用于内环跟踪等一系列问题，需研究一种适用性更广的优化策略。本专利技术提供了一种适应各种角度发射、能考虑远程中的长时作用、能直接用于内环跟踪的弹道优化技术。
[0051]具体步骤如下：
[0052]步骤一，建立发射系动力学方程
[0053][0054]其中，表示导弹质量，表示导弹速度，、、分别为发射系下的坐标；为气动力，P为推力；为地球引力；为离心惯性力；为哥式惯性力；
[0055]气动力和推力需要将体坐标系的分量转化为发射系，即：
[0056][0057]其中为体坐标系和发射系的转换矩阵。
[0058]引力项包含径向和切相引力：
[0059][0060]其中径向，切向，为地理纬度。
[0061]离心惯性力为：
[0062][0063]其中为地球旋转角速度。
[0064]哥式惯性力：
[0065][0066]体坐标系到发射系的转换矩阵首先需要获得姿态角信息，根据由于弹道优化过程忽略内环响应，因此需采用“瞬时平衡假设”。三个欧拉姿态角分别为：
[0067][0068]因此，体坐标系到发射系的转换矩阵为
[0069][0070]高度采用地心高度差进行计算：
[0071][0072]其中弹道任一点的地心距，为发射点地球半径分量，为地球半径。
[0073]步骤二，根据风攻角、风侧滑、风合攻角计算体坐标系下的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种非奇异多约束弹道快速优化方法，其特征在于，包括以下步骤：步骤1：建立发射系动力学方程：；其中，表示导弹质量，表示导弹速度，、、分别为发射系下的坐标；为气动力，P为推力；为地球引力；为离心惯性力；为哥式惯性力；气动力和推力需要将体坐标系的分量转化为发射系，即：；其中为体坐标系和发射系的转换矩阵，表示动压，表示参考面积，、、分别表示轴向力、法向力、侧向力气动系数；地球引力项包含径向引力和切向引力：；其中径向引力，切向引力，为地心距，为地理纬度，、、分别表示发射点在地心固连坐标系下的坐标，、、分别表示地球旋转角速度地心固连坐标系下的坐标；为地球引力常数，第二摄动常数，为地球长半轴；离心惯性力为：；其中为地球旋转角速度；哥式惯性力：
；体坐标系到发射系的转换矩阵需要获得姿态角信息，滚转角、偏航角、俯仰角三个欧拉姿态角分别为：；其中，表示侧滑角，表示弹道倾角，表示攻角，表示弹道倾角；因此，体坐标系到发射系的转换矩阵...

【专利技术属性】
技术研发人员：刘钧圣，苗昊春，高登巍，栗金平，邓海鹏，李琪，常江，潘瑞，范中平，张一杰，张哲明，李雅君，张梦典，张建松，李昊远，杜运理，
申请(专利权)人：西安现代控制技术研究所，
类型：发明
国别省市：