CN117708985A - Machine learning-based high-speed double-spin unsteady aerodynamic prediction method - Google Patents
Machine learning-based high-speed double-spin unsteady aerodynamic prediction method Download PDFInfo
- Publication number
- CN117708985A CN117708985A CN202311745542.7A CN202311745542A CN117708985A CN 117708985 A CN117708985 A CN 117708985A CN 202311745542 A CN202311745542 A CN 202311745542A CN 117708985 A CN117708985 A CN 117708985A
- Authority
- CN
- China
- Prior art keywords
- aerodynamic force
- rotating
- aerodynamic
- double
- speed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000010801 machine learning Methods 0.000 title claims abstract description 21
- 238000012360 testing method Methods 0.000 claims abstract description 19
- 238000004088 simulation Methods 0.000 claims abstract description 13
- 238000012549 training Methods 0.000 claims abstract description 12
- 238000013528 artificial neural network Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000010168 coupling process Methods 0.000 claims description 12
- 238000009987 spinning Methods 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 11
- 238000005859 coupling reaction Methods 0.000 claims description 11
- 241000272517 Anseriformes Species 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000013473 artificial intelligence Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 229920001971 elastomer Polymers 0.000 description 1
- 239000000806 elastomer Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Software Systems (AREA)
- Fluid Mechanics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
提供一种基于机器学习的高速双旋弹非定常气动力预测方法,属于飞行力学和人工智能技术领域。包括生成高速双旋弹气动力建模样本数据集;将数据集划分为训练集和测试集,确定模型输入输出;用神经网络建立从高速双旋弹状态量到所受气动力和力矩的预测模型;将模型用于高速双旋弹气动力和力矩的快速预测。本发明基于机器学习技术,采用非线性自回归神经网络,可很好的学习到双旋弹飞行状态参数到其所受气动力和力矩的映射关系,实现对高速双旋弹非定常气动力和力矩的快速精确预测,减小大规模耦合仿真耗时。本发明建模时考虑了当前时刻和历史时刻的飞行状态对当前时刻气动力和力矩的影响,并考虑了历史时刻气动力和力矩的影响,有效提高预测精度。
Provides a machine learning-based unsteady aerodynamic prediction method for high-speed dual-rotation missiles, which belongs to the fields of flight mechanics and artificial intelligence technology. It includes generating a sample data set for high-speed dual-rotation missile aerodynamic modeling; dividing the data set into a training set and a test set to determine the model input and output; and using a neural network to establish a prediction model from the high-speed dual-rotation missile state quantity to the aerodynamic force and torque received. ;Use the model for rapid prediction of high-speed dual-rotation elastomeric aerodynamic forces and moments. This invention is based on machine learning technology and uses a nonlinear autoregressive neural network, which can well learn the mapping relationship between the flight state parameters of the double-rotating projectile and the aerodynamic force and torque it receives, and realize the unsteady aerodynamic force and torque of the high-speed double-rotating projectile. Fast and accurate prediction, reducing the time-consuming of large-scale coupled simulation. When modeling, the present invention considers the influence of the flight status at the current moment and the historical moment on the aerodynamic force and torque at the current moment, and also considers the influence of the aerodynamic force and moment at the historical moment, thereby effectively improving the prediction accuracy.
Description
技术领域Technical field
本发明属于飞行力学和人工智能技术领域,具体涉及一种基于机器学习的高速双旋弹非定常气动力预测方法。The invention belongs to the technical fields of flight mechanics and artificial intelligence, and specifically relates to a method for predicting the unsteady aerodynamic force of a high-speed dual-rotation missile based on machine learning.
背景技术Background technique
在高速双旋弹设计过程中,其性能参数计算对掌握双旋弹的气动性能和飞行性能都至关重要。CFD/RBD耦合方法可以精确计算和模拟飞行器的气动特性、弹体姿态和运动轨迹,但由于在CFD/RBD耦合计算中需要反复调用CFD分析,使得这种方法效率极低,对计算条件和资源要求高。建立一个高效且精确的非定常气动力模型,并用其代替耦合仿真中的CFD模块,进而对飞行器飞行性能快速、准确地预测。实际上,飞行状态量(如速度、攻角等)和其受到的空气动力和力矩之间存在一定的映射关系,使用机器学习方法可以学习并建立该映射关系并实现双旋弹非定常气动力的快速预测。因此,发展一种高速双旋弹非定常气动力预测方法,是必要并具有广泛需求的。In the design process of high-speed twin-spinning missiles, the calculation of its performance parameters is crucial to mastering the aerodynamic performance and flight performance of the twin-spinning missiles. The CFD/RBD coupling method can accurately calculate and simulate the aerodynamic characteristics, missile body attitude and motion trajectory of the aircraft. However, due to the need to repeatedly call CFD analysis in the CFD/RBD coupling calculation, this method is extremely inefficient and has high requirements on computing conditions and resources. High standard. Establish an efficient and accurate unsteady aerodynamic model and use it to replace the CFD module in coupled simulation to quickly and accurately predict aircraft flight performance. In fact, there is a certain mapping relationship between the flight state variables (such as speed, angle of attack, etc.) and the aerodynamic force and torque it receives. Using machine learning methods, you can learn and establish this mapping relationship and realize the unsteady aerodynamic force of the dual-rotation missile. rapid prediction. Therefore, it is necessary and widely needed to develop a method for predicting unsteady aerodynamic forces of high-speed dual-rotating missiles.
发明内容Contents of the invention
本发明解决的技术问题:提供一种基于机器学习的高速双旋弹非定常气动力预测方法,本发明目的在于为了减小在高速双旋弹设计和性能评估中CFD/RBD耦合计算耗时问题,而通过该方法建立双旋弹飞行状态参数和其受到的气动力/力矩的映射模型,实现对高速双旋弹气动力和力矩的快速准确预测,提升耦合计算效率。The technical problem solved by this invention is to provide a method for predicting the unsteady aerodynamic force of high-speed dual-rotating missiles based on machine learning. The purpose of the invention is to reduce the time-consuming problem of CFD/RBD coupling calculation in the design and performance evaluation of high-speed dual-rotating missiles. , and through this method, a mapping model is established between the flight state parameters of the dual-gyro missile and the aerodynamic force/torque it receives, to achieve rapid and accurate prediction of the aerodynamic force and torque of the high-speed dual-gyro missile, and to improve the efficiency of coupling calculations.
为达到上述目的,本发明采用的技术方案:In order to achieve the above object, the technical solution adopted by the present invention is:
一种基于机器学习的高速双旋弹非定常气动力预测方法,包括以下步骤:A method for predicting unsteady aerodynamic force of high-speed dual-rotating projectiles based on machine learning, including the following steps:
步骤1):生成高速双旋弹非定常气动力建模数据集;Step 1): Generate a high-speed dual-rotation projectile unsteady aerodynamic modeling data set;
步骤2):确定模型的输入输出参数,由于双旋弹的气动力存在显著的非定常特性,因此输入参数包括当前时刻和前m个时刻的飞行状态量,前n个时刻双旋弹所受的气动力/力矩,输出为当前时刻双旋弹所受的气动力/力矩;Step 2): Determine the input and output parameters of the model. Since the aerodynamic force of the twin-spinning projectile has significant unsteady characteristics, the input parameters include the flight state quantities at the current moment and the previous m moments. The aerodynamic force/torque, the output is the aerodynamic force/torque exerted by the double spinning bomb at the current moment;
步骤3):将步骤1)中的数据集划分为训练集和测试集,其中训练集用以训练模型,测试集用于检验模型预测精度;Step 3): Divide the data set in step 1) into a training set and a test set, where the training set is used to train the model, and the test set is used to test the prediction accuracy of the model;
步骤4):利用NARX神经网络,建立从步骤2)中的输入参数到输出参数的预测模型,完成用于高速双旋弹非定常气动力快速预测的气动力模型构建;Step 4): Use the NARX neural network to establish a prediction model from the input parameters in step 2) to the output parameters, and complete the construction of an aerodynamic model for rapid prediction of unsteady aerodynamic forces of high-speed dual-rotation missiles;
步骤5):将步骤4)中的气动力模型用于高速双旋弹非定常气动力和力矩快速预测。Step 5): Use the aerodynamic model in step 4) to quickly predict the unsteady aerodynamic force and torque of the high-speed double-rotating projectile.
上述步骤1)中,具体步骤如下:In the above step 1), the specific steps are as follows:
步骤1-1):选择某一双旋弹构型,生成双旋弹的气动计算网格;Step 1-1): Select a certain double-rotating bomb configuration and generate the aerodynamic calculation grid of the double-rotating bomb;
步骤1-2):确定双旋弹的初始状态参数,包括初始速度、角速度、攻角、侧滑角和姿态角;Step 1-2): Determine the initial state parameters of the double-rotating projectile, including initial velocity, angular velocity, angle of attack, sideslip angle and attitude angle;
步骤1-3):采用步骤1-1)中的气动计算网格,对步骤1-2)中的初始状态进行计算流体力学CFD/刚体动力学RBD耦合仿真,得到双旋弹非定常气动力建模数据集。Step 1-3): Use the aerodynamic calculation grid in step 1-1) to perform computational fluid dynamics CFD/rigid body dynamics RBD coupling simulation on the initial state in step 1-2) to obtain the unsteady aerodynamic force of the double-rotating projectile Modeling dataset.
上述步骤1-1)中,所述双旋弹构型为155mm固定翼鸭式布局双旋弹;步骤1-1)中,所述气动计算网格为非结构混合网格,整套网格共包含60000个表面网格和6613683个体网格;附面层网格使用结构网格,第一层高度为4e-7m,共50层,增长率为1.17;其他区域为非结构网格。In the above step 1-1), the configuration of the double-rotating bomb is a 155mm fixed-wing canard layout double-rotating missile; in the step 1-1), the aerodynamic calculation grid is a non-structural hybrid grid, and the entire set of grids has a total of Contains 60,000 surface grids and 6,613,683 individual grids; the boundary layer grid uses structural grids, the first layer height is 4e-7m, a total of 50 layers, and the growth rate is 1.17; other areas are unstructured grids.
上述步骤1-2)中,双旋弹的初始状态参数如下:马赫数为1.08、前体滚转角速度为0、后体滚转角速度为1570.7rad/s、俯仰角速度和偏航均为0、俯仰角为31.98°、偏航角为0、前体和后体的滚转角均为0、攻角为0°、侧滑角为0°。In the above step 1-2), the initial state parameters of the dual-rotation bomb are as follows: Mach number is 1.08, front body roll angular velocity is 0, rear body roll angular velocity is 1570.7rad/s, pitch angular velocity and yaw are both 0, The pitch angle is 31.98°, the yaw angle is 0, the roll angles of the front and rear bodies are both 0, the angle of attack is 0°, and the sideslip angle is 0°.
上述步骤1-3)中,所述CFD计算时空间离散采用Roe格式,湍流模型为S-A模型;七自由度RBD方程求解采用四阶变步长Adams预估校正方法;对高速双旋弹飞行的0-7秒内进行CFD/RBD耦合仿真,仿真时间步长为0.1ms,共得到70000个样本。In the above steps 1-3), the Roe format is used for spatial discretization in the CFD calculation, and the turbulence model is the S-A model; the seven-degree-of-freedom RBD equation is solved using the fourth-order variable step size Adams prediction correction method; for high-speed dual-rotation missile flight CFD/RBD coupling simulation is performed within 0-7 seconds, the simulation time step is 0.1ms, and a total of 70,000 samples are obtained.
上述步骤2)中,输入参数包括:当前时刻和前4个时刻的飞行状态量,前4个时刻双旋弹所受的气动力/力矩,输出为当前时刻双旋弹所受的气动力/力矩。In the above step 2), the input parameters include: the flight status quantity at the current moment and the previous 4 moments, the aerodynamic force/torque exerted by the dual-rotating bomb at the previous 4 moments, and the output is the aerodynamic force/torque exerted by the dual-rotating bomb at the current moment. moment.
其中,所述飞行状态量包括:前后体各自的三轴速度、三轴角速度、姿态角、攻角和侧滑角;双旋弹所受的气动力/力矩包括前后体各自所受的轴向力、侧向力、法向力、滚转力矩、偏航力矩和俯仰力矩。Among them, the flight state quantities include: the three-axis speed, three-axis angular velocity, attitude angle, angle of attack and sideslip angle of the front and rear bodies respectively; the aerodynamic force/moment suffered by the double-rotating bomb includes the axial direction of the front and rear bodies. force, lateral force, normal force, roll moment, yaw moment and pitch moment.
上述步骤3)中,将步骤1)中的样本划分为训练集和测试集,其中训练集占75%,测试集占25%。In the above step 3), the samples in step 1) are divided into a training set and a test set, of which the training set accounts for 75% and the test set accounts for 25%.
上述步骤4)中,利用NARX神经网络,按步骤2)所述,将当前时刻和前4个时刻的飞行状态量,前4个时刻双旋弹所受的气动力和力矩作为网络的输入,当前时刻双旋弹所受的气动力/力矩作为NARX网络的输出,建立高速双旋弹非定常气动力快速预测的气动力模型并进行训练。In the above step 4), the NARX neural network is used, as described in step 2), and the flight state quantities at the current moment and the previous four moments, and the aerodynamic force and torque experienced by the double-rotating bomb at the previous four moments are used as inputs to the network. The aerodynamic force/torque experienced by the double-rotating projectile at the current moment is used as the output of the NARX network. An aerodynamic model for rapid prediction of the unsteady aerodynamic force of the high-speed double-rotating projectile is established and trained.
上述步骤5)中,具体步骤如下:In the above step 5), the specific steps are as follows:
步骤5-1):从步骤3)中的测试集数据中选择输入参数;Step 5-1): Select input parameters from the test set data in step 3);
步骤5-2):将步骤5-1)中的输入参数作为步骤4)中训练好的气动力模型的输入,预测得到当前时刻双旋弹所受到的气动力和力矩。Step 5-2): Use the input parameters in step 5-1) as the input of the aerodynamic model trained in step 4), and predict the aerodynamic force and torque experienced by the dual-rotation projectile at the current moment.
本发明与现有技术相比的优点:Advantages of the present invention compared with existing technology:
1、本方案基于机器学习技术,采用基于带外源输入的非线性自回归神经网络,可以建立双旋弹飞行状态参数和其受到的气动力/力矩的映射模型,从而能很好的学习到双旋弹飞行状态参数到其所受的气动力和力矩的映射关系,实现对高速双旋弹非定常气动力和力矩的快速、精确预测,提升耦合计算效率,减小在高速双旋弹设计和性能评估中CFD/RBD耦合仿真计算耗时问题;1. This solution is based on machine learning technology and uses a nonlinear autoregressive neural network with external source input. It can establish a mapping model of the flight state parameters of the dual-rotary missile and the aerodynamic force/torque it receives, so that it can learn well The mapping relationship between the flight state parameters of the twin-spinning projectile and the aerodynamic force and torque it receives enables fast and accurate prediction of the unsteady aerodynamic force and torque of the high-speed twin-spinning projectile, improves the efficiency of coupling calculations, and reduces the cost of high-speed twin-spinning projectile design. And the time-consuming problem of CFD/RBD coupling simulation calculation in performance evaluation;
2、本方案所提出的双旋弹非定常气动力预测方法在建立气动力预测模型时,不仅考虑了当前时刻和历史时刻的飞行状态对当前时刻气动力和力矩的影响,还考虑了历史时刻气动力和力矩的影响,有效提高了预测精度;2. When establishing the aerodynamic prediction model for the unsteady aerodynamic prediction method of the double-rotating missile proposed in this plan, it not only considers the influence of the flight status at the current moment and the historical moment on the aerodynamic force and torque at the current moment, but also considers the historical moment. The influence of aerodynamic force and torque effectively improves the prediction accuracy;
3、本方案中建立了双旋弹状态量和其所受的气动力和力矩的映射关系,因此所建立的预测模型可以替代CFD模块与RBD模块进行耦合仿真。3. In this plan, the mapping relationship between the state quantity of the double spin bomb and the aerodynamic force and torque it receives is established. Therefore, the established prediction model can replace the CFD module and the RBD module for coupled simulation.
附图说明Description of the drawings
图1为本发明中基于机器学习的高速双旋弹非定常气动力预测方法的建模及预测的流程框图;Figure 1 is a flow chart of the modeling and prediction of the high-speed dual-rotation projectile unsteady aerodynamic force prediction method based on machine learning in the present invention;
图2为本发明中气动力预测模型的输入输出示意图;Figure 2 is a schematic diagram of the input and output of the aerodynamic prediction model in the present invention;
图3为本发明中采用的双旋弹几何外形示意图;Figure 3 is a schematic diagram of the geometric shape of the double-rotating bomb used in the present invention;
图4为本发明中所采用的双旋弹气动计算网格示意图,(a)为面网格,(b)为滑移网格整体示意图;Figure 4 is a schematic diagram of the double-rotation aerodynamic calculation grid used in the present invention, (a) is a surface grid, (b) is an overall schematic diagram of the sliding grid;
图5为本发明中测试集上高速双旋弹气动力和力矩的部分预测结果,(a)为后体法向力的真实和预测结果对比,(b)为后体侧向力的真实和预测结果对比,(c)为后体俯仰力矩的真实和预测结果对比,(d)为后体偏航力矩的真实和预测结果对比。Figure 5 is a partial prediction result of the aerodynamic force and torque of a high-speed dual-rotation projectile on the test set of the present invention. (a) is a comparison of the real and predicted results of the rear body normal force, and (b) is the real sum of the rear body lateral force. Comparison of predicted results, (c) is the comparison between the actual and predicted results of the rear body pitching moment, (d) is the comparison of the actual and predicted results of the rear body yaw moment.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of the present invention.
请参阅图1-5,详述本发明的实施例。Please refer to Figures 1-5 for a detailed description of the embodiments of the present invention.
一种基于机器学习的高速双旋弹非定常气动力预测方法,参阅图1所示,包括以下步骤:A method for predicting unsteady aerodynamic force of high-speed dual-rotating projectiles based on machine learning, as shown in Figure 1, including the following steps:
步骤1):生成高速双旋弹非定常气动力建模数据集;具体步骤如下:Step 1): Generate a high-speed dual-rotation projectile unsteady aerodynamic modeling data set; the specific steps are as follows:
步骤1-1):选择某一双旋弹构型,生成双旋弹的气动计算网格;Step 1-1): Select a certain double-rotating bomb configuration and generate the aerodynamic calculation grid of the double-rotating bomb;
本实施例中选择的双旋弹构型为155mm固定翼鸭式布局双旋弹,几何外形示意图如图3所示;所述气动计算网格为非结构混合网格,如图4所示,整套网格共包含60000个表面网格和6613683个体网格;附面层网格使用结构网格,第一层高度为4e-7m,共50层,增长率为1.17;其他区域为非结构网格。The configuration of the double-rotating missile selected in this embodiment is a 155mm fixed-wing canard layout double-rotating missile. The geometric shape diagram is shown in Figure 3; the aerodynamic calculation grid is a non-structural hybrid grid, as shown in Figure 4. The entire set of grids contains a total of 60,000 surface grids and 6,613,683 individual grids; the boundary layer grid uses a structural grid, with a first layer height of 4e-7m, a total of 50 layers, and a growth rate of 1.17; other areas are non-structural grids. grid.
步骤1-2):确定双旋弹的初始状态参数,包括初始速度、角速度、攻角、侧滑角和姿态角;Step 1-2): Determine the initial state parameters of the double-rotating projectile, including initial velocity, angular velocity, angle of attack, sideslip angle and attitude angle;
双旋弹的初始状态参数如下:马赫数为1.08、前体滚转角速度为0、后体滚转角速度为1570.7rad/s、俯仰角速度和偏航均为0、俯仰角为31.98°、偏航角为0、前体和后体的滚转角均为0、攻角为0°、侧滑角为0°。The initial state parameters of the dual-rotation bomb are as follows: Mach number is 1.08, front body roll angular velocity is 0, rear body roll angular velocity is 1570.7rad/s, pitch angular velocity and yaw are both 0, pitch angle is 31.98°, yaw The angle is 0, the roll angles of the front and rear bodies are both 0, the angle of attack is 0°, and the sideslip angle is 0°.
步骤1-3):采用步骤1-1)中的气动计算网格,对步骤1-2)中的初始状态进行计算流体力学CFD/刚体动力学RBD耦合仿真,得到双旋弹非定常气动力建模数据集。Step 1-3): Use the aerodynamic calculation grid in step 1-1) to perform computational fluid dynamics CFD/rigid body dynamics RBD coupling simulation on the initial state in step 1-2) to obtain the unsteady aerodynamic force of the double-rotating projectile Modeling dataset.
具体的,以步骤1-2)中的初始状态参数对高速双旋弹飞行的0-7秒内进行CFD/RBD耦合仿真,仿真时间步长为0.1ms,共得到70000个样本;所述CFD计算时空间离散采用Roe格式,湍流模型为S-A模型;七自由度RBD方程求解采用四阶变步长Adams预估校正方法。Specifically, the CFD/RBD coupling simulation was performed on the high-speed dual-rotary missile flight within 0-7 seconds using the initial state parameters in step 1-2). The simulation time step was 0.1ms, and a total of 70,000 samples were obtained; the CFD The Roe format is used for spatial discretization during calculation, and the turbulence model is the S-A model; the fourth-order variable step size Adams prediction correction method is used to solve the seven-degree-of-freedom RBD equation.
步骤2):确定模型的输入输出参数,由于双旋弹的气动力存在显著的非定常特性,因此输入参数包括当前时刻和前m个时刻的飞行状态量,前n个时刻双旋弹所受的气动力/力矩,输出为当前时刻双旋弹所受的气动力/力矩,参阅图2所示;Step 2): Determine the input and output parameters of the model. Since the aerodynamic force of the twin-spinning projectile has significant unsteady characteristics, the input parameters include the flight state quantities at the current moment and the previous m moments. The aerodynamic force/torque of , the output is the aerodynamic force/torque exerted by the double-rotating bomb at the current moment, as shown in Figure 2;
输入参数包括:当前时刻和前4个时刻的飞行状态量,前4个时刻双旋弹所受的气动力/力矩,输出为当前时刻双旋弹所受的气动力/力矩。The input parameters include: the flight state quantity at the current moment and the previous 4 moments, the aerodynamic force/torque experienced by the dual-rotating bomb at the previous 4 moments, and the output is the aerodynamic force/torque experienced by the dual-rotating bomb at the current moment.
其中,飞行状态量包括:前后体各自的三轴速度、三轴角速度、姿态角、攻角和侧滑角;双旋弹所受的气动力/力矩包括前后体各自所受的轴向力、侧向力、法向力、滚转力矩、偏航力矩和俯仰力矩。Among them, the flight state quantities include: the three-axis velocity, three-axis angular velocity, attitude angle, angle of attack and sideslip angle of the front and rear bodies respectively; the aerodynamic force/moment experienced by the double-rotating bomb includes the axial force, lateral force, normal force, roll moment, yaw moment and pitch moment.
步骤3):将步骤1-3)中的数据集划分为训练集和测试集,训练集用以训练模型,测试集用于检验模型预测精度;其中训练集占75%,测试集占25%。Step 3): Divide the data set in step 1-3) into a training set and a test set. The training set is used to train the model, and the test set is used to test the model prediction accuracy; the training set accounts for 75% and the test set accounts for 25%. .
步骤4):利用NARX神经网络,建立从步骤2)中的输入参数到输出参数的预测模型,将当前时刻和前4个时刻的飞行状态量,前4个时刻双旋弹所受的气动力和力矩作为网络的输入,当前时刻双旋弹所受的气动力/力矩作为NARX网络的输出,建立高速双旋弹非定常气动力快速预测的气动力模型并进行训练,完成用于高速双旋弹气动力/力矩快速预测的非定常气动力模型构建。Step 4): Use the NARX neural network to establish a prediction model from the input parameters in step 2) to the output parameters, and combine the flight state quantities at the current moment and the previous 4 moments, and the aerodynamic force experienced by the dual-spinning missile at the previous 4 moments. and torque are used as the input of the network, and the aerodynamic force/torque exerted by the double-rotating projectile at the current moment is used as the output of the NARX network. An aerodynamic model for rapid prediction of the unsteady aerodynamic force of the high-speed double-rotating projectile is established and trained, and the aerodynamic force model for the high-speed double-rotating projectile is completed. Construction of unsteady aerodynamic model for rapid prediction of elastomer force/torque.
步骤5):将步骤4)中的气动力模型用于高速双旋弹非定常气动力和力矩快速预测;具体步骤如下:Step 5): Use the aerodynamic model in step 4) to quickly predict the unsteady aerodynamic force and torque of the high-speed double-rotating projectile; the specific steps are as follows:
步骤5-1):从步骤3)中的测试集数据中选择输入参数;Step 5-1): Select input parameters from the test set data in step 3);
步骤5-2):将步骤5-1)中的输入参数作为步骤4)中训练好的气动力模型的输入,预测得到当前时刻双旋弹所受到的气动力和力矩。参阅图5所示,给出了测试集上高速双旋弹气动力和力矩的部分预测结果(以双旋弹后体为例):(a)为后体法向力的真实和预测结果对比,(b)为后体侧向力的真实和预测结果对比,(c)为后体俯仰力矩的真实和预测结果对比,(d)为后体偏航力矩的真实和预测结果对比。Step 5-2): Use the input parameters in step 5-1) as the input of the aerodynamic model trained in step 4), and predict the aerodynamic force and torque experienced by the dual-rotation projectile at the current moment. Referring to Figure 5, some prediction results of the aerodynamic force and torque of the high-speed double-rotating projectile on the test set are given (taking the rear body of the double-rotating projectile as an example): (a) Comparison of the real and predicted results of the normal force of the rear body , (b) is the comparison of the real and predicted results of the rear body lateral force, (c) is the comparison of the real and predicted results of the rear body pitching moment, (d) is the comparison of the real and predicted results of the rear body yaw moment.
本技术方案基于机器学习技术,采用基于带外源输入的非线性自回归神经网络,可以很好的学习到双旋弹飞行状态参数到其所受的气动力和力矩的映射关系,实现对高速双旋弹非定常气动力和力矩的快速、精确预测。本发明在建立双旋弹气动力预测模型时,不仅考虑了当前时刻和历史时刻的飞行状态对当前时刻气动力和力矩的影响,还考虑了历史时刻气动力和力矩的影响,有效提高了预测精度。This technical solution is based on machine learning technology and uses a nonlinear autoregressive neural network with external source input. It can well learn the mapping relationship between the flight state parameters of the dual-rotary missile and the aerodynamic force and torque it receives, and realize high-speed Fast and accurate prediction of unsteady aerodynamic forces and moments of dual-rotating projectiles. When establishing the dual-gyro aerodynamic prediction model, the present invention not only considers the influence of the flight status at the current moment and the historical moment on the aerodynamic force and torque at the current moment, but also considers the influence of the aerodynamic force and moment at the historical moment, effectively improving the prediction Accuracy.
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其它的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。It is obvious to those skilled in the art that the present invention is not limited to the details of the above-described exemplary embodiments, and that the present invention can be implemented in other specific forms without departing from the spirit or essential characteristics of the invention. Therefore, the embodiments should be regarded as illustrative and non-restrictive from any point of view, and the scope of the present invention is defined by the appended claims rather than the above description, and it is therefore intended that all claims falling within the claims All changes within the meaning and scope of equivalent elements are included in the present invention. Any reference signs in the claims shall not be construed as limiting the claim in question.
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其它实施方式。In addition, it should be understood that although this specification is described in terms of implementations, not each implementation only contains an independent technical solution. This description of the specification is only for the sake of clarity, and those skilled in the art should take the specification as a whole. , the technical solutions in each embodiment can also be appropriately combined to form other implementations that can be understood by those skilled in the art.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311745542.7A CN117708985A (en) | 2023-12-18 | 2023-12-18 | Machine learning-based high-speed double-spin unsteady aerodynamic prediction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311745542.7A CN117708985A (en) | 2023-12-18 | 2023-12-18 | Machine learning-based high-speed double-spin unsteady aerodynamic prediction method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117708985A true CN117708985A (en) | 2024-03-15 |
Family
ID=90155007
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311745542.7A Pending CN117708985A (en) | 2023-12-18 | 2023-12-18 | Machine learning-based high-speed double-spin unsteady aerodynamic prediction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117708985A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117952037A (en) * | 2024-03-22 | 2024-04-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | High-speed aircraft aerodynamic engineering estimation correction method based on deep learning |
CN118849000A (en) * | 2024-09-27 | 2024-10-29 | 摩多利智能传动(江苏)有限公司 | Torque control method of planetary gear reduction device for robot |
-
2023
- 2023-12-18 CN CN202311745542.7A patent/CN117708985A/en active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117952037A (en) * | 2024-03-22 | 2024-04-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | High-speed aircraft aerodynamic engineering estimation correction method based on deep learning |
CN117952037B (en) * | 2024-03-22 | 2024-05-28 | 中国空气动力研究与发展中心超高速空气动力研究所 | High-speed aircraft aerodynamic engineering estimation correction method based on deep learning |
CN118849000A (en) * | 2024-09-27 | 2024-10-29 | 摩多利智能传动(江苏)有限公司 | Torque control method of planetary gear reduction device for robot |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117708985A (en) | Machine learning-based high-speed double-spin unsteady aerodynamic prediction method | |
CN104298805B (en) | Hypersonic aircraft CFD aerodynamic modeling methods | |
CN104881510B (en) | A kind of lifting airscrew/tail-rotor aerodynamic interference numerical value emulation method | |
CN108182328A (en) | A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter | |
CN108090302B (en) | Helicopter flight mechanics simulation method and system | |
CN114355777B (en) | A dynamic gliding method and system based on distributed pressure sensors and segmented attitude control | |
CN109760046A (en) | Motion planning method for capturing rolling target of space robot based on reinforcement learning | |
CN114004023B (en) | Method for identifying aerodynamic parameters of aircraft based on cyclic neural network | |
CN105426955A (en) | Disturbance-based elite reverse learning particle swarm optimization implementation method | |
CN116560384A (en) | Variant aircraft robust control method based on deep reinforcement learning | |
Wang et al. | Glide performance analysis of underwater glider with sweep wings inspired by swift | |
CN117452957A (en) | A method for identifying aircraft aerodynamic parameters based on physical information neural network | |
Wang et al. | Deep learning based missile trajectory prediction | |
CN102566446A (en) | Method for establishing full-envelope mathematical model of unmanned helicopter based on linear model group | |
CN104914736A (en) | Method for building general simulation model for hypersonic flight vehicles | |
Li et al. | Deep learning-based surrogate model for flight load analysis | |
Baayen | Vortexje-an open-source panel method for co-simulation | |
CN111506959A (en) | Universal flight control system and method for flight simulation training | |
Cheng et al. | Study on Small UAVs' Deep Stall Landing Procedure | |
CN107944201B (en) | A fast modeling method for fixed-wing aircraft | |
CN109376482B (en) | Projectile motion envelope refined prediction method based on digital prototype | |
Han et al. | Investigation of the unsteady aerodynamic characteristics of an unmanned aerial vehicle with variable-sweep morphing | |
CN118485005B (en) | A tiltrotor thrust prediction method and system based on multi-precision neural network | |
Katz | Classical and potential flow based aerodynamics-do we need them? | |
CN113761653B (en) | Shrapnel separation safety evaluation method and system and shrapnel separation simulation platform |
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 |