1、文章编号:0258-1825(2023)09-0030-08基于物理信息神经网络的飞机气动参数辨识方法付军泉1,2,钟伯文1,2,*,钟运琴3,刘赟1,2(1.南昌航空大学飞行器工程学院,南昌330063;2.江西省飞行器设计与气动仿真重点实验室,南昌330063;3.中国科学院大学中国科学院大数据挖掘与知识管理重点实验室,北京100190)摘要:在飞机设计与研制过程中,通过气动参数辨识建立可靠的飞行动力学模型非常重要。传统的气动参数辨识工程算法,诸如极大似然法,需要给出合理的飞行动力学模型以及待辨识参数的初值。基于传统神经网络的气动参数辨识可以避免飞行动力学建模过程,这种方法需要通过增量法、
2、导数法间接地从神经网络提取气动参数。本文提出了一种基于物理信息神经网络的飞机气动参数辨识方法,可将含待辨识参数的飞行动力学模型作为正则项加入损失函数,直接辨识得到气动参数。该方法可以显著减少建模数据需求,也能提高建模精度。飞行仿真数据验证结果表明,该方法的无噪声、含 2%噪声仿真数据,纵向飞行状态空间模型辨识最大相对误差分别为 1.80%、4.64%,表明了基于物理信息神经网络的飞机气动参数辨识方法具有可行性,并对含噪声的飞行数据具有泛化性。关键词:深度学习;基于物理信息的神经网络;飞行动力学;状态空间;气动参数辨识中图分类号:TP183;V212.1文献标识码:Adoi:10.7638/kq
3、dlxxb-2022.0201A physics informed neural network based method for aircraftaerodynamic parameter identificationFUJunquan1,2,ZHONGBowen1,2,*,ZHONGYunqin3,LIUYun1,2(1.School of Aircraft Engineering,Nanchang Hangkong University,Nanchang330063,China;2.Jiangxi Key Laboratory for Aircraft Design and Aerody
4、namic Simulation,Nanchang330063,China;3.Key Laboratory of Big Data Mining and Knowledge Management,University of Chinese Academy of Sciences,Beijing100190,China)Abstract:Itisofgreatimportancetoestablishareliableflightdynamicmodelbasedonaerodynamicparameter identification in the process of aircraft d
5、esign.Traditional engineering methods for aerodynamicparameteridentification,suchasthemaximumlikelihoodmethodetc.,requireareasonabledynamicmodelandinitial values for the parameters to be identified.The identification of aerodynamic parameters based ontraditionalneuralnetworkshasnorequirementforthefl
6、ightdynamicmodelandparameterinitialization.Inthisstudy,anaerodynamicparameteridentificationmethodbasedonphysicsinformedneuralnetwork(PINN)isproposed,whichcandirectlyidentifyaerodynamicparametersbyincorporatingtheflightdynamicmodelasaregularizationtermintothelossfunction.Thismethodcansignificantlyred
7、ucetherequiredtrainingdataandimprove the modeling accuracy.Verification results for the flight simulation data show that the maximumrelativeerrorsofthePINNbasedmethodonthelongitudinalflightstatespacemodelidentificationwithnonoiseand 2%noise are 1.80%and 4.64%,respectively.It suggests that the PINN b
8、ased aerodynamic parameteridentificationmethodisfeasibleandhasthegeneralizationcapacitytonoisyflightdata.Keywords:deep learning;physics informed neural network;flight dynamics;state space;parameteridentification收稿日期:2022-12-28;修订日期:2023-02-27;录用日期:2023-03-13;网络出版时间:2023-04-07基金项目:国家自然科学基金(12262023);
9、江西省重点研发计划项目(20223BBE51001)作者简介:付军泉(1993),男,江西抚州人,讲师,研究方向:空气动力学.E-mail:通信作者:钟伯文*,研究员,研究方向:空气动力学.E-mail:引用格式:付军泉,钟伯文,钟运琴,等.基于物理信息神经网络的飞机气动参数辨识方法J.空气动力学学报,2023,41(9):3037.FU J Q,ZHONG B W,ZHONG Y Q,et al.A physics informed neural network based method for aircraft aerodynamic parameteridentificationJ.Ac
10、taAerodynamicaSinica,2023,41(9):3037(inChinese).doi:10.7638/kqdlxxb-2022.0201第41卷第9期空气动力学学报Vol.41,No.92023年9月ACTA AERODYNAMICA SINICASept.,2023 0 引言在飞机设计与研制过程中,建立精确的飞行动力学模型是飞行控制系统设计、飞行品质评估、实物及半实物飞行仿真的基础与前提。飞机气动参数一般可通过计算流体力学1、风洞试验2、参数辨识3等方法获得。计算流体力学方法由于存在 N-S 方程离散误差、计算截断误差、湍流模型误差等因素,导致精确的飞行流场很难求解。风洞试
11、验中由于存在洞壁效应、相似性效应、传感器测量误差等因素,也使得精确的气动参数获取存在困难。随着传感器技术与飞行试验技术迅速发展,通过飞机气动参数辨识方法、基于飞行试验数据来获取更高精度的气动参数成为一种可行的、必要的手段。飞机气动参数辨识是一个系统工程,主要包括飞行试验设计、飞行动力学模型结构假设、飞机气动参数辨识、模型辨识结果验证四个步骤。1920 年,Waener 和 Norton 对飞行器的气动参数辨识方法进行了早期研究4。1951 年,Harry5提出了一种以飞机小扰动线性化模型为基础,基于最小二乘法辨识飞机稳定导数与操纵导数的方法。目前,应用最广泛的飞机气动参数辨识方法是极大似然法与
12、牛顿拉夫逊迭代的组合方法。基于极大似然估计的模型参数辨识问题本质上是参数优化问题,但该方法是一种收敛缓慢的算法。Lin 等6通过改进极大似然估计方法的参数化梯度的递推中间计算过程,提高了极大似然方法的收敛速度。Kalman 滤波算法能将参数估计问题转化为状态估计问题,依据飞行器的实际飞行情况,采用连续估计模型,并通过离散时间测量与离散滤波算法达到实时估计,其中应用最广泛的非线性气动参数辨识方法是扩展 Kalman 滤波(extendedKalmanfilter,EKF)方法。此外,无迹卡尔曼滤波(unscentedKalmanfilter,UKF)方法能降低模型线性化所带来的误差,因此在辨识非
13、线性气动模型时有较明显的优势7。极大似然法是目前应用于飞机气动参数辨识的主要算法。极大似然法辨识需要给出合理的飞行动力学模型结构以及待辨识参数初始值,不准确的动力学模型会导致气动数据拟合结果不理想,不合理的待辨识参数初始值则容易导致优化陷入局部最小值,甚至导致辨识过程发散8。近年来,人工智能技术的日益成熟为飞机气动参数辨识提供了一种新方法。神经网络可以逼近任何非线性输入输出系统,是一种系统辨识建模的有效工具9。基于神经网络的系统辨识可以避免飞行动力学建模过程,并且在气动参数辨识过程中不需要给出待辨识参数的初始值估计。21 世纪初期开始,国内外学者均开展了传统神经网络在飞机气动参数辨识的应用研究
14、,这类方法用神经网络模拟飞行动力学系统的输入输出关系,并利用增量法、导数法从神经网络间接地提取气动参数,是一种数据驱动的辨识建模方法。Ghosh 等10基于真实飞行数据建立舵面控制输入与气动力、力矩之间的前馈神经网络模型,并研究舵面控制输入叠加微小增量后的神经网络输出变化,通过增量法提取气动参数。Kuttieri 等11研究了基于前馈神经网络的弹性体飞机的气动参数辨识,在建立升降舵、飞行状态输入与升力、俯仰力矩输出之间神经网络模型的基础上,利用深度学习框架自带的微分求导功能,直接求解升力、俯仰力矩输出相对于升降舵、飞行状态输入的导数,得到气动参数。2020 年,Verma 等12基 于 极 限
15、 学 习 机 器(extreme learningmachine,ELM),一种改进的前馈神经网络,并利用导数法辨识基于真实飞行数据的飞机气动参数,结果显示,该方法能显著加速气动力、力矩的神经网络模型的训练过程。随着深度学习理论与技术的持续发展,2017 年美国布朗大学 Raissi 等13综合并改进了前人在神经网络求解微分方程方面的技术方法,正式提出了基于物理信息的神经网络(physicsinformedneuralnetwork,PINN)算法框架,并首次提出 PINN 方法不仅可以求解微分方程正问题(即建模问题),也可用于求解微分方程逆问题(即参数辨识问题)。PINN 是一种数据与物理双
16、驱动的建模框架,可以将物理约束作为正则项加入到传统神经网络建模过程中,从而减少对建模数据的需求。同时,PINN 相对于传统神经网络,也具有更高的辨识精度14。PINN 被正式提出后,引发了流体力学、热力学、生物医学、材料力学等领域大量的相关应用研究15。比较典型的网络结构有应用于不可压 N-S 方程的分数阶 fPINNs 网络16,应用于心房颤动诊断分布型 DPINNs 网络17,应用于新材料设计的贝叶斯型 B-PINNs 网络18等。基于 PINN 网络的系统辨识方法应用研究目前尚不多见。本文研究的重点是根据待辨识动力学系统的特征,设置 PINN网络的结构和损失函数。文章首先给出通用 PIN
17、N 网络的基本结构以及工作原理,然后给出基于 PINN 网络的飞机纵向状态空间模型的气动参数辨识建模方法,最后利用基于JSBSim 的飞行仿真数据验证上述建模方法的正确性第9期付军泉等:基于物理信息神经网络的飞机气动参数辨识方法31与有效性。1 PINN设含参数的偏微分方程的通用形式为:f(x,t,u,xu,tu,;)=0,x ,t 0,T(1)方程的初始条件和边界条件定义为:u(x,t0)=g(x),x (2)u(x,t)=h(t),x ,t 0,T(3)x RdtRdfxu tu =1,2,.fu(x,t)fg(x)fh(x)f其中:空间变量;为时间变量;为欧式空间的子集;为的边界;微分方
18、程 包含微分算子、;为微分方程 的参数;为偏微分方程 的理论解;为 的初始条件;为 的狄利克雷边界条件、黎曼边界条件或混合边界条件。u(x,t)fu(x,t)LPINN 用一个神经网络逼近偏微分方程的解。不妨定义 PINN 网络逼近微分方程 的解,该 PINN 网络的损失函数 定义为:L=LPDE+LData+LIC+LBC(4)LPDELIC其中:为损失函数中模型驱动的部分;为损失LBCLData函数中初始条件驱动的部分;为损失函数中边界条件驱动的部分;为损失函数中其他数据驱动的部分。损失函数各部分定义如下:LPDE=f(u,t(k)u,x(k)u;)(5)LData=u|u|Data(6)
19、LIC=u|,t0g|,t0(7)LBC=(n u|nh|)+(u|h|)(8)tx u(x,t)u(x,t)在损失函数(4)的约束下,在尽可能满足偏微分方程的条件下,使 PINN 网络的训练数据的误差尽量小。含参数的偏微分方程(1)对应的 PINN 网络结构如图 1 所示,在输入时间 和空间 数据后,首先用全连接神经网络逼近函数,然后利用神经网络自动微分技术19,求出偏微分方程残差、初值、边值残差,并将其作为正则项加入损失函数中,最后利用梯度下降法更新神经网络权重参数和偏微分方程物理参数。需要说明的是,PINN 广义上是一种基于物理信息神经网络的框架,主要用于求解偏微分方程,但也适用于常微分
20、方程的求解。txIt(k)xn(k)u(k)(k)LPDE=f(u,t u,x u;)LData=u|u|DataLIC=u|,t0g|,t0LBC=(nu|nh|)+(u|h|)图 1 PINN 网络结构示意图Fig.1 Schematic of PINN 2 PINN 气动参数辨识方法以飞机纵向运动为例,研究基于 PINN 网络的飞机气动参数辨识方法。假设飞机在定常平飞状态只受到升降舵激励输入扰动,飞机的纵向运动与横航向运动在该小扰动作用下可以完全解耦,则典型的飞机纵向运动小扰动线性化状态方程可表示为20:u w q=XuXw0gcos0ZuZwu0gsin0MuMwMq00010uwq+
21、XeZeMe0e(9)uwqe其中:为流向速度;为垂向速度;为俯仰角速率;为俯仰角;为升降舵偏转角。如无特殊说明,文中u w q e所有与飞机飞行状态相关的参数,如、等,均为扰动增量。进一步假设在上述飞机纵向激励输入下,观测变量与状态变量完全一致,则此处观测方程可略去。y=f(t,y)本文飞机气动参数辨识方法研究采用离散型PINN 网络。参考 Raissi 等14提出的 PINN 求解离散型偏微分方程逆问题的方法,构造离散型 PINN 系统辨识网络一般需要采用高阶 Runge-Kutta 算法离散偏微分方程。本文采用三阶 Runge-Kutta 公式离散飞机纵向小扰动线性化状态方程模型。记一阶
22、微分方程通用形式为,其三阶 Runge-Kutta 公式为:yk+1=yk+(K1+4K2+K3)dt/6(10)K1=f(tk,yk)(11)32空气动力学学报第41卷K2=f(tk+dt/2,yk+K1dt/2)(12)K3=f(tk+dt,ykK1dt+2K2dt)(13)uwqfufwfq为后文表述方便,将飞机纵向运动小扰动线性化状态方程(9)中与速度、速度、俯仰角速率 相关的微分方程分别记为、:fu=Xuu+Xwwgcos0+Xee(14)fw=Zuu+Zwwgsin0+Zee(15)fq=Muu+Mww+Mqq+Zee(16)fufwfqu w q按式(10)三阶 Runge-Ku
23、tta 公式将微分方程、依次离散化,得到关于、的递推公式:uk+1=uk+(K1,u+4K2,u+K3,u)dt/6(17)wk+1=wk+(K1,w+4K2,w+K3,w)dt/6(18)qk+1=qk+(K1,q+4K2,q+K3,q)dt/6(19)ukwkqkkekuk+1wk+1qk+1由飞机纵向小扰动线性化状态方程的离散化公式(1719),需要构造 PINN 网络,其输入为、,输出为、。基于 PINN 的飞机气动参数辨识全连接神经网络结构参数见表 1,该神经网络共有 9 层,激活函数均采用 Relu。表 1 PINN 气动参数辨识网络结构参数Table 1 Structure pa
24、rameters of PINN for aerodynamicparameter identification层数输入节点输出节点激活函数15128Relu2128256Relu3256512Relu4512512Relu5512512Relu6512256Relu7256128Relu812864Relu9643-ukwkqkkek uk+1 wk+1 qk+1LDataLODEuLODEwLODEqPINN 气动参数辨识网络结构如图 2 所示。输入、后,首先用 PINN 全连接神经网络预测输出、,预测输出与真实输出的残差构成数据驱动的损失函数,然后将飞机纵向小扰动线性化状态方程的离散化公
25、式(1719)中计算输出与预测输出的残差、作为模型驱动的部分加入损失函数中,最后利用梯度下降法更新 PINN 全连接网络权重参数和状态方程中的待辨识气动参数。wkukuk+1LDataq=(uk+1uk+1)+(wk+1wk+1)+(qk+1qk+1)wk+1qk+1ekLODEu=uk+1uk(K1,u+4K2,u+K3,u)dt/6LODEw=wk+1wk(K1,w+4K2,w+K3,w)dt/6LODEq=qk+1qk(K1,q+4K2,q+K3,q)dt/6图 2 PINN 气动参数辨识网络结构示意图Fig.2 Schematic of PINN for aerodynamic par
26、ameter identification 3 飞行仿真研究采用的飞行仿真数据基于 JSBSim 求解器产生。JSBSim 是一个可在多种操作系统下编译和运行的开源飞行动力学模型库。JSBSim 既可以作为一个独立的程序运行,从 xml 格式的脚本文件和配置文件中获取输入,也可以集成到一个支持交互输入的飞行模拟平台中21,如 FlightGear 和 OpenEaagles。飞机气动参数辨识需要在飞机配平状态下施加合理的舵面控制输入,以充分激励出飞行动力学响应过程。目前基于时域的气动参数辨识常用的激励输入信号为倍脉冲信号与 3211 信号22。图 3 所示为采用的 PINN 网络纵向气动参数辨
27、识升降舵含噪声倍脉冲激励信号与 3211 激励信号。倍脉冲激励信号可近似为单位周期内的正弦信号,通过对倍脉冲激励信号宽度的选择,可以将激励所覆盖的频段移至飞行仿真所需要激发的模态频段。同理,在 3211 激励信号中,通常将脉冲“2”的宽度设为飞机固有模态频率一半所对应的时间,脉冲“3”与“1”分别设置在飞机固有模态频率的两侧,以形成频带更宽的激励输入。在默认设置下,JSBSim 输出的飞行动力学响应数据是无噪声的,因为 JSBSim 会将激励输入信号送入一个低通滤波器滤波,再迭代求解飞行动力学方程。本文主要研究基于无噪声飞行仿真数据的 PINN第9期付军泉等:基于物理信息神经网络的飞机气动参数
28、辨识方法33网络纵向气动参数方法,同时为了研究该 PINN 网络辨识方法的泛化性,也研究了该方法在含噪声的飞行仿真数据上的辨识效果。含噪声的飞行仿真数据u(t)定义如下:u(t)=u(t)+rand(t)(20)u(t)(0,1)u(t)rand(t)N(0,1)其中:为原始信号;表示噪声水平;为原始信号的标准差;标准随机噪声。MuXe=0.14ZeMe飞行仿真所采用机型为“DELTA”23,与纵向运动相关的定常平飞配平状态为:u0=75m/s,w0=0m/s,q0=0rad/s,0=0.047rad。飞机纵向小扰动状态方程(9)中,“DELTA”飞机的值相对于其他气动参数为小量,在飞行仿真中
29、将该值设为 0。其他待辨识纵向气动参数的真实值为:Xu=0.02,Xw=0.1,Zu=0.23,Zw=0.634,=2.9,Mw=0.005,Mq=0.61,=0.44。在 JSBSim 飞机气动力配置文件中设置上述飞行动力学气动参数,在JSBSim 初始化配置文件中设置上述定常平飞配平状态为仿真初始状态,在 JSBSim 脚本设置文件中设置飞行状态参数的输出频率为 64Hz。分别以 10 个幅值与周期不同的升降舵含噪声倍脉冲信号和 10 个幅值与周期不同的升降舵含噪声 3211 信号激励作为输入,通过 JSBSim 求解器得到无噪声飞行仿真数据。含噪声飞行仿真数据通过在无噪声数据中添加噪声水
30、平 2%的随机噪声获得。图 4 为升降舵含噪声 3211激励输入下的纵向飞行参数随时间变化图。5.02.5000.51.01.52.02.5时间/s(a)倍脉冲激励2.5升降舵偏转角/()5.05.02.500213457689时间/s(b)3211激励2.5升降舵偏转角/()5.0图 3 升降舵倍脉冲信号与 3211 激励信号Fig.3 Doublet and 3211 activation signal of the elevator5.02.50024(a)升降舵偏转角随时间变化图(b)速度u随时间变化图(c)俯仰角速率随时间变化图(d)俯仰角随时间变化图68时间/s2.5升降舵偏转角/
31、()5.04200速度u速度w2468时间/s2速度/(ms1)64202468时间/s0俯仰角/()2002468时间/s2俯仰角速率/(s1)图 4 升降舵 3211 激励下纵向飞行参数随时间变化图Fig.4 Time history of longitudinal flight parameters under 3211 activation of the elevator34空气动力学学报第41卷 4 模型验证气动参数辨识算法模型基于 PyTorch 开源机器学习框架实现。随机从 JSBSim 飞行仿真数据中选择5000 组数据作为 PINN 网络气动参数辨识方法研究的数据集,并将该数
32、据集以 80%和 20%的比例随机划分为训练集与测试集。PINN 网络训练过程中主要超参数设置为:训练数据批量大小为 512,学习率为0.001,所有训练集数据迭代训练次数为 10000。基于无噪声飞行仿真数据的 PINN 网络飞机纵向气动参数辨识结果及对应误差如表 2 所示,基于含噪声飞行仿真数据的气动参数辨识结果及对应误差Xe如表 3 所示。无噪声时,9 个辨识参数中辨识误差最大的是 Mw,其相对真值的误差为 1.80%,这主要是因为 Mw相对于其他待辨识气动参数小一个数量级,在PINN 网络误差反向传播、梯度下降及参数更新过程中相对更容易产生累积误差。含噪声时,待辨识参数中辨识结果误差最
33、大的是,其相对真值的误差为4.64%。同时,含噪声时的待辨识参数的相对误差比无噪声时的结果整体上略微偏大。为进一步验证基于 PINN 网络的飞机气动参数辨识方法的效果,分别选取一组不在训练数据集中的完整升降舵激励信号,测试无噪声、含噪声 PINN 网络的预测输出速度 u、速度 w、俯仰角速率 q。图 5、表 2 无噪声时 PINN 气动参数辨识结果及对应误差Table 2 Results and errors of aerodynamic parameteridentification by PINN without noise参数真实值预测值相对误差/%Xu0.02000.020030.15
34、Xw0.10000.10000Xe0.14000.14030.21Zu0.23000.23010.04Zw0.63400.63400Ze2.90002.90400.14Mw0.00500.004911.80Mq0.61000.61360.60Me0.44000.44150.34表 3 含噪声时 PINN 气动参数辨识结果及对应误差Table 3 Results and errors of aerodynamic parameteridentification by PINN with noise参数真实值预测值相对误差/%Xu0.02000.019950.25Xw0.10000.09871.3
35、0Xe0.14000.14654.64Zu0.23000.22681.39Zw0.63400.63190.33Ze2.90002.93641.26Mw0.00500.005132.60Mq0.61000.60391.00Me0.44000.43940.145.02.500123时间/s2.5升降舵偏转角/()5.00.01250.01000.00750.5无噪声真实值无噪声预测值0.70.80.90.61.0时间/s0.00500.0025速度u/(ms1)0.20.40.5无噪声真实值无噪声预测值0.70.80.90.61.0时间/s0.60.8速度w/(ms1)0.50.5无噪声真实值无噪
36、声预测值0.70.80.90.61.0时间/s1.01.5俯仰角速率/(s1)(a)升降舵偏转角随时间变化图(b)速度u随时间变化图(c)速度w随时间变化图(d)俯仰角速率随时间变化图图 5 无噪声时 PINN 网络的预测输出与真实值对比图Fig.5 Comparison between the PINN prediction and the real value without noise第9期付军泉等:基于物理信息神经网络的飞机气动参数辨识方法35图 6 分别为无噪声时一组升降舵倍脉冲激励、含噪声时一组升降舵 3211 激励输入下,PINN 网络 u、w、q 的预测输出与真实值对比图。为了
37、更清晰地展示对比情况,无噪声时选取 t=0.51.0s 时间范围,含噪声时选取 t=2.02.5s 时间范围。由图可见,无论纵向动力学飞行仿真数据中是否含有噪声,相对速度 u、速度 w、俯仰角速率 q 的预测误差总是偏大,这也与表 2、表 3 中气动参数辨识相对误差的分布表现一致。无噪声时,与 u、w、q 相关的气动参数辨识的最大相对误差分别为 0.15%、0.14%、1.80%;含噪声时,与u、w、q 相关的气动参数辨识的最大相对误差分别为4.64%、1.39%、2.60%;无论是否含有噪声,与 w 相关的气动参数的最大相对误差总是小于与 u、q 相关的气动参数的最大相对误差。5002468
38、时间/s升降舵偏转角/()50.1750.1500.1252.0含噪声真实值含噪声预测值2.22.32.42.12.5时间/s0.1000.075速度u/(ms1)3.02.0含噪声真实值含噪声预测值2.22.32.42.12.5时间/s3.5速度w/(ms1)2.92.0含噪声真实值含噪声预测值2.22.32.42.12.5时间/s3.13.23.0俯仰角速率/(s1)(a)升降舵偏转角随时间变化图(b)速度u随时间变化图(c)速度w随时间变化图(d)俯仰角速率随时间变化图图 6 含噪声时 PINN 网络的预测输出与真实值对比图Fig.6 Comparison between the PIN
39、N prediction and the real value with noise由于本文采用梯度下降法训练 PINN 网络,其中与 u、w、q 相关的物理方程正则项约束在加入损失函数时权重均设为 1,而与 w 相关的物理方程正则项以更大的梯度更新相关气动参数,所以导致与 w 相关的气动参数的相对误差较小。Wang 等24指出,在PINN 网络训练过程中,可以通过调整损失函数中不同组成成分之间的权重,使得其梯度值大小趋于一致,来改善 PINN 网络的收敛情况。一般来说,在飞机气动参数辨识工程应用领域,气动参数辨识结果的相对误差不大于 10%就可采纳25。综上所述,基于 PINN 网络的飞机气
40、动参数辨识方法可行且有效,并对含噪声飞行数据的气动参数辨识具有一定的泛化能力。5 结论本文在传统神经网络气动参数辨识方法基础上,提出了一种基于 PINN 网络的飞机气动参数辨识方法。该方法可以将含待辨识参数的飞机飞行动力学方程作为正则项加入到损失函数,通过 PINN 网络梯度下降与参数更新,直接辨识得到气动参数。与极大似然法等工程辨识算法相比,该方法不需要给出待辨识参数的初始估计。与基于传统神经网络气动参数辨识方法相比,基于 PINN 网络的辨识建模方法可以显著减少建模数据需求,也能提高建模精度。同时,针对纵向飞行状态空间模型辨识,通过 JSBSim 飞行仿真验证,获得了在无噪声和含 2%噪声
41、仿真数据时9 个待辨识参数中辨识结果的最大相对误差(分别为1.80%、4.64%)。研究结果表明,基于 PINN 网络的飞机气动参数辨识方法具有可行性,并对含噪声飞行数据的气动参数辨识具有一定的泛化能力。后期,针对基于 PINN 网络的非线性飞行动力学系统(特别是含有时间延迟、过程噪声等复杂场景)36空气动力学学报第41卷的气动参数辨识,需要结合真实飞行数据进行进一步验证。参 考 文 献:张淼,刘铁军,马涂亮,等.基于CFD方法的大型客机高速气动设计J.航空学报,2016,37(1):244254.ZHANGM,LIUTJ,MATL,etal.Highspeedaerodynamicdesig
42、noflarge civil transporter based on CFD methodJ.Acta Aeronautica etAstronauticaSinica,2016,37(1):244254(inChinese).1刘斌,刘沛清,王亮.某小型边条翼无人飞机的气动参数估算及风洞试验分析J.飞机设计,2010,30(3):15,22.LIU B,LIU P Q,WANG L.A small trim-wing unmanned aircraftaerodynamicparameterestimationandanalysisofwindtunneltestJ.AircraftDes
43、ign,2010,30(3):15,22(inChinese).doi:10.3969/j.issn.1673-4599.2010.03.0012蔡金狮.飞行器气动参数辨识进展J.力学进展,1987,17(4):467478.CAI J S.Advances in identification of aircraft aerodynamicparametersJ.Advances in Mechanics,1987,17(4):467478(inChinese).3WARNER E,NORTON F.Preliminary report on free flight testsR.NACA-T
44、R-70,1920.https:/ntrs.nasa.gov/api/citations/19930091132/downloads/19930091132.pdf4HARRYG.AsurveyofmethodsfordeterminingstabilityparametersofanairplanefromdynamicflightmeasurementsR.NACA-TN-2340,1951.https:/ntrs.nasa.gov/api/citations/19930082979/downloads/19930082979.pdf5LIN D M,MAKIS V.On-line par
45、ameter estimation for a partiallyobservable system subject to random failureJ.Naval ResearchLogistics,2006,53(5):477483.doi:10.1002/nav.201566杨枝山.基于组合导航的助推滑翔导弹气动在线辨识与弹道规划方法研究D.长沙:国防科学技术大学,2017.YANGZS.Researchonaerodynamiconlineidentificationandtrajectoryoptimizationbasedonintegratednavigationforboos
46、t-glidemissileD.Changsha:National University of Defense Technology,2017(inChinese).7崔平远.飞行器参数辨识方法与应用研究J.飞行力学,1996,14(3):18.CUI P Y.Research on the parameter identification method andapplication for vehicleJ.Flight Dynamics,1996,14(3):18(inChinese).doi:10.13645/ki.f.d.1996.03.0018杨广慧,杜立夫,李辉,等.基于BP神经网
47、络的飞行器参数辨识与自适应控制J.航天控制,2021,39(5):37.YANGGH,DULF,LIH,etal.ParameteridentificationandadaptivecontrolofaircraftbasedonBPneuralnetworkJ.AerospaceControl,2021,39(5):37(inChinese).9GHOSH A K,KHAN A S,RAISINGHANI S C.Estimation oflongitudinalparametersfromreal-flightdatausingneuralnetworksJ.JournaloftheAer
48、onauticalSocietyofIndia,2002,54(1):1622.10KUTTIERIRA,SINHAM.Neuralpartialdifferentiationforparameterestimation of flexible aircraft dynamicsJ.Journal of AerospaceSciencesandTechnologies,2016,68(4):276286.11VERMAHO,PEYADANK.AircraftparameterestimationusingELMnetworkJ.Aircraft Engineering and Aerospac
49、e Technology,2020,92(6):895907.12doi:10.1108/aeat-01-2019-0003RAISSIM,PERDIKARISP,KARNIADAKISGE.Physicsinformeddeep learning(part I):data-driven solutions of nonlinear partialdifferentialequationsR/OL.arXiv:1711.10561,2017.https:/arxiv.org/pdf/1711.10561.pdf13RAISSIM,PERDIKARISP,KARNIADAKISGE.Physic
50、sinformeddeep learning(part II):Data-driven discovery of nonlinear partialdifferentialequationsR/OL.arXiv:1711.10566,2017.https:/arxiv.org/pdf/1711.10566.pdf14李野,陈松灿.基于物理信息的神经网络:最新进展与展望J.计算机科学,2022,49(4):254262.LIY,CHENSC.Physics-informedneuralnetworks:recentadvancesand prospectsJ.Computer Science,2