基于可插拔式接口的固定翼飞行器快速建模设计方法技术

本发明专利技术提出基于可插拔式接口的固定翼飞机气动参数编辑规范，设计可插拔式接口，将气动参数从动力学模型中剥离并抽象化，实现气动参数即插即用。使用VC2010开发工具，完成固定翼飞机的操纵输入、动力学方程和运动学方程模块开发。设计气动参数XML文件解析程序，将获取的气动参数同动力学模块相结合，实现飞行性能仿真模型的快速开发。用户通过在XML文件中替换不同飞行器的气动参数，在发动机类型与型号固定的前提下，对不同的飞行器建模时只需修改XML文件中的气动参数便可完成仿真模型设计。可快速建立针对不同固定翼飞机机型的空气动力学模型，缩短模型开发周期。减少在飞行器数学建模过程中的重复性劳动，提高了飞行器建模的灵活性、快速性，以及通用性。

【技术实现步骤摘要】

本专利技术具体涉及一种基于可插拔式接口的固定翼飞行器快速建模设计方法，属于计算机飞行器性能仿真领域。
技术介绍
XML为可扩展标记语言或可延伸，是一种置标语言。本专利技术的创新点就在于应用可扩展标记语言定于气动参数的XML文档，将飞行器的气动参数从飞行方程中抽离出来写入XML文档，在进行仿真的时候再将气动参数从文档中解析到飞行方程中就可以进行飞行器性能仿真。使用XML文档标记气动参数可以使得在对不同飞行器进行仿真时可共用飞行方程，只需要修改文档中的气动参数即可，如此可提高开发效率，实现代码的复用。XML语言具有简洁有效、高效扩充的特点，并且使用XML可以在不同的计算机系统间交换信息。XML支持复用文档片段，使用者可以专利技术和使用自己的标签，可伸缩性大，可以有效的进行XML文件的扩充，在飞行器平台不同的时候可以通过修改扩充XML文件快速的进行飞行器仿真建模。
技术实现思路
本专利技术提出了一种基于可插拔式接口的固定翼飞行器快速建模设计方法，提出基于XML的固定翼飞机气动参数编辑规范，设计可插拔式接口，将气动参数从动力学模型中剥离并抽象化，实现气动参数即插即用。使用VC2010开发工具，完成固定翼飞机的操纵输入、动力学方程和运动学方程模块开发。设计气动参数XML文件解析程序，将获取的气动参数同动力学模块相结合，实现飞行性能仿真模型的快速开发。用户通过在XML文件中替换不同飞行器的气动参数，在发动机类型与型号固定的前提下，对不同的飞行器建模时只需修改XML文件中的气动参数便可完成仿真模型设计。可快速建立针对不同固定翼飞机机型的空气动力学模型，缩短模型开发周期。减少在飞行器数学建模过程中的重复性劳动，提高了飞行器建模的灵活性、快速性，以及通用性。一种基于可插拔式接口的固定翼飞行器快速建模设计方法，主要包含以下内容：步骤一：定义XML文档，将飞行器的气动参数从飞行方程中剥离并写入文件中，实现飞行方程在初始化时从中读取气动参数；步骤二：根据飞行器的运动分解，将其运动分为横侧向运动与纵向运动，分别编写运动状态解算方程；步骤三：设计XML解析程序，将气动参数加载到空气动力学模型中，实现运动状态解算方程的初始化；步骤四：在发动机类型与型号固定的条件下，对于不同的飞行器可以复用同一飞行方程，只需修改XML文件中的气动参数再重新解析到动力学模型中即可；所述步骤一中定义XML文档的实现思路为：⑴.首先由XML编辑器创建XML文档，定义文档根节点；⑵.定义第一层节点，共8个节点，分别标记第一层节点数、升力解算相关气动参数、侧力解算气动参数、俯仰力矩解算气动参数、滚转力矩解算气动参数、偏航力矩解算气动参数、阻尼力解算气动参数、角加速度解算气动参数；⑶.为第一层的每个节点创建子节点，根据对应节点的气动参数个数创建子节点，涉及多少个气动参数就创建几个子节点，其中第一个子节点标记当前节点的子节点个数；⑷.一个子节点代表一个气动参数，由于存入文档的气动参数为一系列离散点，因此需要为子节点创建一系列里层节点。每个子节点扩展了21个里层节点，其中第一个里层节点标记气动参数离散点个数，其余里层节点标记气动参数离散值；(5).保存并关闭XML文档，实现气动参数的标记；所述步骤三中XML文档解析实现思路为：⑴.加载XML文档，开始进行气动参数解析；⑵.将文档指针定位到根节点，获取第一层节点的节点数；⑶.将文档指针依次定位到第一层的每个节点，开始遍历每个节点的子节点；⑷.遍历子节点扩展的里层节点，将里层节点中标记的气动参数值，读取到对应气动解算模块中对应的气动参数数组中；⑸.重复上述操作，直到所有节点的子节点中的里层节点全部读取为止；⑹.保存并关闭XML文档，将启动参数加载到计算机内存中；所述第一层节点代表对应的气动力、力矩解算模块；节点的子节点代表该气动解算模块涉及的气动参数；子节点的里层节点代表该气动参数标记的一维表数据个数；需要写入XML文件的气动参数表：为便理解，附如下机翼升力系数的XML文档标记格式。具体实现：1定义变量：变量模块包括飞行器的各气动参数对应的抽象变量、运动状态变量、力学状态变量以及控制输入。其中气动参数变量用于飞行器空气动力模型的解算；运动状态变量用于记录模型解算出的姿态信息、位置信息等；力学状态变量用于记录模型解算出的俯仰力矩、横滚力矩，偏航力矩、侧向力矩以及阻力、升力等状态信息；控制输入包括[δT δe δa δr]，对应油门推杆输入、升降舵偏角、副翼偏转角、方向舵偏角。2定义XML文档及XML解析程序：将飞行器空气动力学模型的气动参数写入XML文件中。从空气动力学模型中将气动参数剥离出来并抽象化，在各解算模块中以变量代替，通过XML解析程序将文件中的气动参数解析到对应的变量中，实现空气动力学模型的初始化。XML文档的定义是将表1中的气动参数按照解算方程以树状层次结构依次记录到文档中，这样飞行器的所有气动参数就构成了具有多层分叉结构的树，每一个气动参数对应为其中一个分支的一个元素。解析程序的作用就是分析树，通过进行树遍历就可以对每个元素进行操作，将其解析到飞行器飞行方程中，从而将飞行方程具体化。创建XML文档的方法见说明书附图1、图2，XML文档解析方法见说明书附图3。3状态解算模块：状态解算模块包括飞行器的力方程组、运动方程组、力矩方程组、导航方程组四部分，编写状态解算模块时将对应的方程组离散化转换成计算机程序。在确定了状态向量与控制输入[δT δe δa δr]之间的关系后，在已知飞行器相关的气动参数、特征参数，很据飞行高度h，马赫数Ma以及飞行状态，就可以确定力(Fx Fy Fz)和力矩应用状态解算模块便可以求解飞行器在任何时刻的运动状态。其中(u,v,w)为机体坐标系的三个速度分量；为姿态角，分别为滚转角，俯仰角，航向角；(p,q,r)为机体角速度分量，分别为机体滚转角速度、机体俯仰角速度以及机体偏航角速度；(xg,yg,zg)为地面坐标；δT,δe,δa,δr分别为油门开度输入、升降舵偏输入、副翼舵偏输入以及航向舵偏输入。1).根据控制输入[δT δe δa δr]解算力、力矩。力主要包括升力L、侧力Y、发动机推力T。具体解算如下：升力解算方程：其中有Q：动压；ρ：当前气压高度空气密度；V：飞机水平飞行速度；CLW：机翼升力系数；CLb：机身升力系数；Sb：机身横截面积；SW：机翼面积；CLt：平尾升力系数；St：平尾面积；侧力解算方程：ρ：当前气压高度空气密度；V：飞机水平飞行速度；SW：机翼面积；CYβ：侧滑角引起的侧力导数；β：侧滑角；：方向舵侧力导数；δr：方向舵偏角输入；：滚转角速度侧力导数；：以为量纲的滚转角速度，b为机翼展长；：偏航角速度侧力导数；：以为量纲的偏航角速度，b为机翼展长；俯仰力矩解算方程： M A = ( 1 2 ρV 2 ) ( C m , α = 本文档来自技高网...

【技术保护点】
一种基于可插拔式接口的固定翼飞行器快速建模设计方法，其特征在于，其包含以下步骤：步骤一，定义XML文档，将飞行器的气动参数从飞行方程中剥离并写入文件中，实现飞行方程在初始化时从中读取气动参数；步骤二，根据飞行器的运动分解，将其运动分为横侧向运动与纵向运动，分别编写运动状态解算方程；步骤三，设计XML解析程序，将气动参数加载到空气动力学模型中，实现运动状态解算方程的初始化；步骤四，在发动机类型与型号固定的条件下，对于不同的飞行器可以复用同一飞行方程，只需修改XML文件中的气动参数再重新解析到动力学模型中即可。

【技术特征摘要】
1.一种基于可插拔式接口的固定翼飞行器快速建模设计方法，其特征在于，其包含以下步骤：步骤一，定义XML文档，将飞行器的气动参数从飞行方程中剥离并写入文件中，实现飞行方程在初始化时从中读取气动参数；步骤二，根据飞行器的运动分解，将其运动分为横侧向运动与纵向运动，分别编写运动状态解算方程；步骤三，设计XML解析程序，将气动参数加载到空气动力学模型中，实现运动状态解算方程的初始化；步骤四，在发动机类型与型号固定的条件下，对于不同的飞行器可以复用同一飞行方程，只需修改XML文件中的气动参数再重新解析到动力学模型中即可。2.根据权利要求1所述的一种基于可插拔式接口的固定翼飞行器快速建模设计方法，其特征在于，所述步骤一中定义XML文档的具体步骤如下：⑴.首先由XML编辑器创建XML文档，定义文档根节点；⑵.定义第一层节点，共8个节点，分别标记第一层节点数、升力解算相关气动参数、侧力解算气动参数、俯仰力矩解算气动参数、滚转力矩解算气动参数、偏航力矩解算气动参数、阻尼力解算气动参数、角加速度解算气动参数；⑶.为第一层的每个节点创建子节点，根据对应节点的气动参数个数创建子节点，涉及多少个气动参数就创建几个子节点，其中第一个子节点标记当前节点的子节点个数；⑷.一个子节点代表一个气动参数，由于存入文档的气动参数为一系列离散点，因此需要为子节点创建一系列里层节点；每个子节点扩展了21个里层节点，其中第一个里层节点标记气动参数离散点个数，其余里层节点标记气动参数离散值；(5).保存并关闭XML文档，实现气动参数的标记。3.根据权利要求1所述的一种基于可插拔式接口的固定翼飞行器快速建模设计方法，其特征在于，所述步骤三中XML文档解析的具体步骤如下：⑴.加载XML文档，开始进行气动参数解析；⑵.将文档指针定位到根节点，获取第一层节点的节点数；⑶.将文档指针依次定位到第一层的每个节点，开始遍历每个节点的子节点；⑷.遍历子节点扩展的里层节点，将里层节点中标记的气动参数值，读取到对应气动解算模块中对应的气动参数数组中；⑸.重复上述操作，直到所有节点的子节点中的里层节点全部读取为止；⑹.保存并关闭XML文档，将启动参数加载到计算机内存中。4.根据权利要求1所述的一种基于可插拔式接口的固定翼飞行器快速建模设计方法，其特征在于，所述步骤三中飞行器运动状态解算方程，包括以下内容：所述状态解算方程包括飞行器的力方程组、运动方程组、力矩方程组、导航方程组四部分，编写状态解算模块时将对应的方程组离散化转换成计算机程序；在确定了状态向量与控制输入[δT δe δa δr]之间的关系，在已知飞行器相关的气动参数、特征参数，很据飞行高度h，马赫数Ma以及飞行状态，就可以确定力(Fx Fy Fz)和力矩应用状态解算模块便可以求解飞行器在任何时刻的运动状态；其中(u,v,w)为机体坐标系的三个速度分量；为姿态角，分别为滚转角，俯仰角，航向角；(p,q,r)为机体角速度分量，分别为机体滚转角速度、机体俯仰角速度以及机体偏航角速度；(xg,yg,zg)为地面坐标；δT,δe,δa,δr分别为油门开度输入、升降舵偏输入、副翼舵偏输入以及航向舵偏输入；1).根据控制输入[δT δe δa δr]解算力、力矩；力主要包括升力L、侧力Y、发动机推力T；具体解算如下：升力解算方程：其中有Q：动压；ρ：当前气压高度空气密度；V：飞机水平飞行速度；CLW：机翼升力系数；CLb：机身升力系数；Sb：机身横截面积；SW：机翼面积；CLt：平尾升力系数；St：平尾面积；侧力解算方程：ρ：当前气压高度空气密度；V：飞机水平飞行速度；SW：机翼面积；CYβ：侧滑角引起的侧力导数；β：侧滑角；方向舵侧力导数；δr：方向舵偏角输入；滚转角速度侧力导数；以为量纲的滚转角速度，b为机翼展长；偏航角速度侧力导数；以为量纲的偏航角速度，b为机翼展长；俯仰力矩解算方程：ρ：当前气压高度空气密度；V：飞机水平飞行速度；Cm,α＝0：静不定俯仰力矩；Cmα：迎角俯仰力矩导数；α：迎角；升降舵俯仰力矩导数；δe：升降舵偏角；平尾附加俯仰力矩系数；以为量纲的俯仰角速度，b为机翼展长；平尾下洗时差阻尼力矩导数；以为量纲的迎角速度升降舵偏角速率俯仰力矩导数；以为量纲的升降舵偏角速率；cA：机翼平均几何弦长；滚转力矩解算方程：Clβ：横滚静稳定导数；V：飞机水平飞行速度；b：b为机翼展长；ρ：当前气压高度空气密度；滚转操纵导数；δa：副翼偏转角；方向舵操纵交叉导数；横滚阻尼导数；以为量纲的滚转角速度，b为机翼展长；交叉动导数；以为量纲的偏航角速度，b为机翼展长；偏航力矩解算方程：ρ：当前气压高度空气密度；V：飞机水平飞行速度；Cnβ：偏航静稳定导数；副翼操纵交叉导数；δa：副翼偏转角；航向操纵导数；δr：航向舵偏角；交叉动导数；以为量纲的滚转角速度，b为机翼展长；航向阻尼导数；以为量纲的偏航角速度，b为机翼展长；阻尼力解算方程：Dk＝CDVk (6)CD：阻尼系数；由方程(1)～(6)可计算出飞行器沿机体坐标系力向量： F k x F k y F k z = T k 0 0 + 0 Y k - L k + - D k - m k g k s i n θ m k g k cos θ sin φ m k g k cos θ cos φ - - - ( 7 ) ]]>其中k为第k个解算周期；Dk为飞行阻力；mk为第k个解算周期飞机的质量；gk为第k个解算周期的重力加速度；由方程(1)～(6)可计算飞行器沿机体坐标系分解力矩向量： L k M k N k = 0 T k l z - T k l y + L ‾ k A M k A N k A - - - ( 8 ) ]]>其中MkA，NkA分别为第k个解算周期滚转操纵力矩、俯仰操纵力矩、偏航操纵力矩；Tk为第k个解算周期的发动机推力；ly,lz分别为y轴的转动惯量、z轴的转动惯量；Lk,Mk,Nk分别为滚转力矩、俯仰力矩、偏航力矩；2).由[Fkx Fky Fkz]T、[uk νk ωk]T、[φk θk ψk]T、[pk qk rk]T解算出当前状态加速度分量；加速度分量解算方程：由方程(9)解算出当前状态加速度分量3).由[Lk Mk Nk]T、[pk qk rk]T解算出当前状态角加速度分量；角加速度分量解算方程：其中：Ixz为惯性积； c 1 = ( I y - I z ) I z - I x z 2 I x I z - I x z 2 , c 2 = ( I x - I y + I z ) I x z I x I z - I x z 2 , c 3 = I z I x I z - I x z 2 ; ]]> c 4 = I x z I x I z - I x z 2 , c 5 = I z - I x I y , c 6 = I x z I y ; ]]> c 7 = 1 I y , c 8 = ( I x - I y ) I x + I x z 2 I x I z - I x z 2 , c 9 = I x I x I z - I x z 2 ; ]]>由方程(10)解算出当前状态角加速度分量4).由[uk νk ωk]T、解算出下一状态的速度分量[uk+1 νk+1 ωk+1]T下一状态的速度分量解算方程：其中Δτ为解算周期；5).由[pk qk rk]T、解算出下一状态的角速度分量[pk+1 qk+1 rk+1]T下一状态的角速度分量解算方程： p k + 1 q k + 1 r k + 1 = p k ...

【专利技术属性】
技术研发人员：孟光磊，罗元墙，梁宵，田丰，薛继佳，朱琳琳，喻勇涛，
申请(专利权)人：沈阳航空航天大学，
类型：发明
国别省市：辽宁;21