CN108733858B - 应用于高空飞行器***的建模方法和装置 - Google Patents
应用于高空飞行器***的建模方法和装置 Download PDFInfo
- Publication number
- CN108733858B CN108733858B CN201710267450.0A CN201710267450A CN108733858B CN 108733858 B CN108733858 B CN 108733858B CN 201710267450 A CN201710267450 A CN 201710267450A CN 108733858 B CN108733858 B CN 108733858B
- Authority
- CN
- China
- Prior art keywords
- altitude aircraft
- aircraft system
- rigid body
- generalized
- matrix
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims description 97
- 238000005381 potential energy Methods 0.000 claims description 24
- 230000009977 dual effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 23
- 230000005484 gravity Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000005437 stratosphere Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- 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
-
- 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
-
- 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)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种应用于高空飞行器***的建模方法和装置。其中,该方法包括:将高空飞行器***的物理模型简化为双刚体***,其中,浮空器和吊舱均被简化为一个刚体;确定高空飞行器***的广义坐标,其中,广义坐标至少包括高空飞行器***中每个刚体的姿态角信息;求解高空飞行器***的拉格朗日函数并建立拉格朗日方程;依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息。本发明解决了相关技术中对高空飞行器***进行建模仿真的结果相对于实际情况的误差较大的技术问题。
Description
技术领域
本发明涉及仿真建模分析领域,具体而言,涉及一种应用于高空飞行器***的建模方法和装置。
背景技术
高空飞行器***包括浮空器和与浮空器柔性连接的吊舱,浮空器可以是气球、飞艇等,吊舱通常为载荷舱,浮空器与吊舱之间通常通过旋转连接器或系绳等球铰连接方式连接。在对高空飞行器***的设计放飞过程中,高空飞行器***的动力学模型是控制***设计的重点与基础,动力学性能决定了高空飞行器***所能够达到的性能指标。目前在对高空飞行器***进行建模时,通常将浮空器与吊舱当作一个整体,使用经典的牛顿力学方法,求取其动量和其对质心的动量矩,然后对时间求导,进而建立飞行器在合外力F作用下的平移运动和在合外力矩M作用下的转动运动方程,只有高空飞行器的空间位置与姿态角6个自由度,甚至在某些时刻简化为位置3个自由度,而不考虑飞行器姿态,导致对高空飞行器***进行建模仿真的结果相对于实际情况的误差较大。而且在动力***安装在载荷舱或者载荷舱载人飞行时,通常需要确切知道载荷舱的姿态信息,载荷舱与气球通过旋转连接器连接,实现姿态运动分离。现有建模方法无法体现气球姿态与载荷舱姿态分离的特点,无法得到载荷舱的姿态角信息,以便进行旅游观光,也无法有效进行载荷舱与气球之间的动力学耦合分析。
针对上述的技术问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种应用于高空飞行器***的建模方法和装置,以至少解决相关技术中对高空飞行器***进行建模仿真的结果相对于实际情况的误差较大的技术问题。
根据本发明实施例的一个方面,提供了一种应用于高空飞行器***的建模方法,该方法包括:将高空飞行器***的物理模型简化为双刚体***,其中,浮空器和吊舱均被简化为一个刚体;确定高空飞行器***的广义坐标,其中,广义坐标至少包括高空飞行器***中每个刚体的姿态角信息;求解高空飞行器***的拉格朗日函数并建立拉格朗日方程;依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息。
进一步地,确定双刚体***的广义坐标包括:建立高空飞行器***中每个刚体的位置坐标系和姿态角坐标系;确定广义坐标,其中,广义坐标包括每个刚体的姿态角和高空飞行器***的位置。
进一步地,求解高空飞行器***的拉格朗日函数包括:求解高空飞行器***中每个刚体的动能矩阵和势能矩阵;根据高空飞行器***中每个刚体的动能矩阵和势能矩阵确定高空飞行器的拉格朗日函数矩阵。
进一步地,求解高空飞行器***中每个刚体的动能矩阵和势能矩阵包括:基于广义坐标确定高空飞行器***中每个刚体的广义速度矩阵,其中,广义速度矩阵为通过广义坐标的求导式表示的刚体的速度和角速度的矩阵;通过广义速度矩阵表示高空飞行器***中每个刚体的动能矩阵。
进一步地,建立拉格朗日方程包括:获取高空飞行器的广义受力向量,其中,广义受力向量为高空飞行器***所受到的力的集合通过广义坐标表示的向量;根据拉格朗日函数矩阵和广义受力向量建立拉格朗日方程。
进一步地,依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息包括:利用矩阵乘法运算规则执行拉格朗日方程中的求导运算。
根据本发明实施例的另一方面,还提供了一种应用于高空飞行器***的建模装置,高空飞行器***包括柔性连接的浮空器和吊舱,该装置包括:逻辑单元,用于将高空飞行器***的物理模型简化为双刚体***,其中,浮空器和吊舱均被简化为一个刚体;确定单元,用于确定高空飞行器***的广义坐标,其中,广义坐标至少包括高空飞行器***中每个刚体的姿态角信息;建模单元,用于求解高空飞行器***的拉格朗日函数并建立拉格朗日方程;仿真单元,用于依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息。
进一步地,确定单元包括:第一建立模块,用于建立高空飞行器***中每个刚体的位置坐标系和姿态角坐标系;第一确定模块,用于确定广义坐标,其中,广义坐标包括每个刚体的姿态角和高空飞行器***的位置。
进一步地,建模单元包括:求解模块,用于求解高空飞行器***中每个刚体的动能矩阵和势能矩阵;第二确定模块,用于根据高空飞行器***中每个刚体的动能矩阵和势能矩阵确定高空飞行器的拉格朗日函数矩阵。
进一步地,求解模块包括:确定子模块,用于基于广义坐标确定高空飞行器***中每个刚体的广义速度矩阵,其中,广义速度矩阵为通过广义坐标的求导式表示的刚体的速度和角速度的矩阵;逻辑子模块,用于通过广义速度矩阵表示高空飞行器***中每个刚体的动能矩阵。
进一步地,建模单元包括:获取模块,用于获取高空飞行器的广义受力向量,其中,广义受力向量为高空飞行器***所受到的力的集合通过广义坐标表示的向量;第二建立模块,用于根据拉格朗日函数矩阵和广义受力向量建立拉格朗日方程。
进一步地,仿真单元包括:运算模块,用于利用矩阵乘法运算规则执行拉格朗日方程中的求导运算。
根据本发明实施例的另一方面,还提供了一种存储介质,该存储介质包括存储的程序,其中,在程序运行时控制存储介质所在设备执行本发明的应用于高空飞行器***的建模方法。
根据本发明实施例的另一方面,还提供了一种处理器,该处理器用于运行程序,其中,程序运行时执行本发明的应用于高空飞行器***的建模方法。
在本发明实施例中,通过将高空飞行器***的物理模型简化为双刚体***,其中,浮空器和吊舱均被简化为一个刚体;确定高空飞行器***的广义坐标,其中,广义坐标至少包括高空飞行器***中每个刚体的姿态角信息;求解高空飞行器***的拉格朗日函数并建立拉格朗日方程;依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息,解决了相关技术中对高空飞行器***进行建模仿真的结果相对于实际情况的误差较大的技术问题,进而实现了能够减小对高空飞行器***进行建模仿真的误差的技术效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的一种可选的应用于高空飞行器***的建模方法的流程图;
图2是根据本发明实施例的一种可选的高空飞行器***的坐标系示意图;
图3是根据本发明实施例的一种可选的高空飞行器***的建模方法的仿真结果的示意图;
图4是根据本发明实施例的另一种可选的高空飞行器***的建模方法的仿真结果的示意图;
图5是根据本发明实施例的另一种可选的高空飞行器***的建模方法的仿真结果的示意图;
图6是根据本发明实施例的一种可选的应用于高空飞行器***的建模装置的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
本申请提供了一种应用于高空飞行器***的建模方法的实施例。需要说明的是,该实施例提供的方法应用的高空飞行器***包括柔性连接的浮空器和吊舱。浮空器可以是气球、飞艇等,吊舱为载荷舱,浮空器与吊舱之间可以通过旋转连接器或系绳等球铰连接方式连接。
图1是根据本发明实施例的一种可选的应用于高空飞行器***的建模方法的流程图,如图1所示,该方法包括如下步骤:
步骤S101,将高空飞行器***的物理模型简化为双刚体***,其中,浮空器和吊舱均被简化为一个刚体。浮空器和吊舱之间的柔性连接方式如果是系绳连接或旋转连接器连接,则可以将柔性连接简化为球铰连接。
步骤S102,确定高空飞行器***的广义坐标。
广义坐标要至少包括高空飞行器***中每个刚体的姿态角信息。确定广义坐标要先建立坐标系。具体的,可以在建立高空飞行器***中每个刚体的位置坐标系和姿态角坐标系之后确定广义坐标,其中,广义坐标包括每个刚体的姿态角和高空飞行器***的位置。
如图2所示为对一种可选的高空飞行器***建立的坐标系的示意图。该实施例的方法应用在图2中所示的高空气球***中,高压气球***中的浮空器(飞行器)为气球,图2中高压气球***的气球为超压气球101,吊舱(载荷舱)为桁架103,超压气球101与桁架103之间通过旋转连接器102连接,桁架103还设置有助推器,具体为图2中的两个螺旋桨:螺旋桨104和螺旋桨105,这两个螺旋桨可以调节吊舱的前进方向。
对图2中所示的高空飞行器***定义如下坐标系:
a)地面坐标系Oexeyeze:原点Oe为地面某一固定点,xe轴的正向沿正东方向,ze轴的正向垂直地面指向地心,ye轴与xe轴、ze轴构成右手系。
b)气球本体坐标系Oxyz:该坐标系与超压气球101固连。原点O位于超压气球101的浮心,x轴的正向为垂直于两个螺旋桨安装中心连线向前,x轴的正向表示桁架103前进的方向,z轴的正向为在过轴的对称面内垂直轴向下,y轴与x轴、z轴构成右手系。
c)载荷舱本体坐标系Oaxayaza:该坐标系与桁架103固连。原点Oa位于桁架103的质心,xa轴的正向垂直于两个螺旋桨安装中心连线且向前,xa轴的正向表示桁架103前进的方向,za轴的正向在过轴的对称面内垂直轴向下,ya轴与xa轴、za轴构成右手系。
在建立坐标系之后,选取广义坐标:
选取气球质心位置作为高空飞行器***的位置,广义坐标包括气球质心位置、气球姿态角、载荷舱姿态角共9个广义坐标。气球质心位置在地面惯性坐标系中表示为{x,y,z},姿态角θ、ψ分别为气球的滚转角、俯仰角、偏航角,载荷舱姿态角/>θ1、ψ1分别为载荷舱的滚转角、俯仰角、偏航角,广义坐标的表达式如下:
步骤S103,求解高空飞行器***的拉格朗日函数并建立拉格朗日方程。
求解高空飞行器***的拉格朗日函数可以通过如下方式求解:求解高空飞行器***中每个刚体的动能矩阵和势能矩阵,根据高空飞行器***中每个刚体的动能矩阵和势能矩阵确定高空飞行器的拉格朗日函数矩阵。其中,求解高空飞行器***中每个刚体的动能矩阵和势能矩阵可以包括:基于广义坐标确定高空飞行器***中每个刚体的广义速度矩阵,其中,广义速度矩阵为通过广义坐标的求导式表示的刚体的速度和角速度的矩阵;通过广义速度矩阵表示高空飞行器***中每个刚体的动能矩阵。
具体而言,在图2所示的高空气球***中,求解高空飞行器***的拉格朗日函数的方法如下:
(1)确定气球和载荷舱的广义速度矩阵:
假设气球的角速度在气球本体坐标系下表示为ω1,载荷舱的角速度在载荷舱本体坐标系下表示为ω2,由角速度定义可以得到
根据气球的速度v1、载荷舱的速度v2、气球的角速度ω1和载荷舱的角速度ω2的关系得到图2所示的高空气球***的广义速度矩阵如下:
其中,E是单位矩阵,A1为地面坐标系Oexeyeze到气球本体坐标系Oxyz的坐标变换矩阵,A2为地面坐标系Oexeyeze到载荷舱本体坐标系Oaxayaza的坐标变换矩阵,S(r)表示向量r的叉乘矩阵,l1为在气球本体系Oxyz下表示的气球质心到球铰位置的矢量,l2为在吊舱本体坐标系表示的吊舱质心到球铰位置矢量。
(2)基于广义速度矩阵表示高空飞行器***中每个刚体的动能矩阵:
1)气球的动能T1:
其中,m1为气球质量,I1为气球关于其质心的惯性张量矩阵在气球本体坐标系下的矩阵表达式。
2)吊舱的动能T2:
其中,m2为吊舱质量,I2为吊舱关于其质心的惯性张量矩阵在吊舱本体坐标系下的矩阵表达式。
(3)确定高空飞行器***中每个刚体的势能矩阵:
1)气球的势能
以海平面为零势能参考面,气球的势能V1表示为:
V1=m1r1 Tg1=m1g1z
其中,g1=[0 0 g1]T为气球所在高度的重力加速度在地面坐标系下的矩阵表达式。
2)吊舱的势能
以海平面为零势能参考面,吊舱的势能V2表示为:
其中,g2=[0 0 g2]T为吊舱所在高度的重力加速度在地面坐标系下的矩阵表达式。
(4)根据高空飞行器***中每个刚体的动能矩阵和势能矩阵确定高空飞行器***的拉格朗日函数矩阵:
高空飞行器***的拉格朗日函数为:
L=T1+T2-V1-V2
在求解高空飞行器***的拉格朗日函数之后,可以建立拉格朗日方程。建立拉格朗日方程的步骤可以是:获取高空飞行器的广义受力向量,其中,广义受力向量为高空飞行器***所受到的力的集合通过广义坐标表示的向量,然后根据拉格朗日函数矩阵和广义受力向量建立拉格朗日方程。
具体而言,在图2所示的高空气球***中,建立拉格朗日方程的步骤如下:
(1)获取广义受力向量QN
假设高空气球***受到非保守力系{FN1 … FNm}的作用(受到的非保守力系中包括m个力,例如,浮力、重力、附加惯性力、气动力和主动控制力等),其中,力FNi(i=1,…,m)的作用点矢径pi基于高空气球***的广义坐标q可以表示为pi(q)(i=1,…,m),则高空气球***所受的对应于q的非保守的广义受力向量QN为
(2)根据拉格朗日函数矩阵L和广义受力向量QN建立拉格朗日方程
步骤S104,依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息。
依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息可以包括:利用矩阵乘法运算规则执行拉格朗日方程中的求导运算。
广义力中包括主动控制力和非主动控制力,其中,浮力、重力、附加惯性力以及气动力都是高空气球***中一般不会主动加以控制的,是非主动控制力,对于图2所示的高空气球***,主动控制力为设置在桁架103双侧的两个螺旋桨的推力,应用建模分析时通常考虑双侧推进螺旋桨的不同工作方式加以仿真分析。
采用上述实施例提供的高空飞行器***的建模方法对图2所示的高空气球***在双侧推进螺旋桨的不同工作方式下的载荷舱姿态角进行仿真分析,得到如下结果:
(1)双侧推进螺旋桨不工作,高空气球***悬浮在平流层,浮力和重力平衡,考虑风速为0,仿真得到载荷舱的姿态角信息,如图3所示。根据仿真结果,在无风情况下,浮力重力平衡,载荷舱的姿态角都为0,不会变化。
(2)双侧推进螺旋桨提供相等的前进推力,使得高空气球***前进,同样考虑浮力重力平衡,风速为0,得到的载荷舱姿态角,如图4所示。在无风情况下,螺旋桨推力向前作用时,会使得载荷舱有一定的俯仰角,通过旋转连接器102连接气球向前运动,俯仰角呈现周期性变化,而滚转角与偏航角不会发生变化。
(3)双侧推进螺旋桨提供不同的前向推力,以实现载荷舱的航向改变,浮力重力平衡,考虑风速为0,得到载荷舱的姿态角如图5所示。在无风情况下,双侧推进螺旋桨实现偏航时,载荷舱偏航角会变化,同时由于有一定的前向推力,载荷舱俯仰角也会有一定变化。
由上面三个具体的算例可以看出,该实施例提供的建模方法可以有效实现载荷舱姿态与气球姿态的分离,分析载荷舱姿态、载荷舱与气球之间的动力学耦合关系。
该实施例提供的建模方法将高空气球***简化为气球与载荷舱通过球铰连接的双刚体***,通过定义气球位置和姿态角以及载荷舱姿态角为广义坐标,然后求解气球与载荷舱的动能与势能,采用拉格朗日方法进行动力学建模,在需要分析载荷舱姿态角信息时,可以有效地将气球与载荷舱的姿态角分离分析,不同于通常的单刚体动力学模型,
该实施例提供的建模方法选取的广义坐标由于包括了气球的姿态角与载荷舱的姿态角信息,使得建立的动力学模型可以模拟仿真分析载荷舱的姿态信息,以及气球与载荷舱的动力学耦合关系,在载人旅行与航向控制时得到的仿真结果更为准确。
在依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息时,可以利用矩阵乘法运算规则执行拉格朗日方程中的求导运算。由于现有技术中的建模方法在获得高空飞行器***的动能之后,需要对动能求偏导之后再对时间求导得到广义速度的导数,如果采用这种计算方式逐项展开求导,求解算法复杂,计算量大。本发明实施例采用矩阵表达形式表示气球与载荷舱的速度和角速度,写成广义速度矩阵的表达形式,并将动能表示为矩阵形式,使得在求导时可以利用矩阵乘法进行求导运算,分化运算量,对动能进行求导运算(在求解动能偏导数之后求解对于时间的导数)更简便。
需要说明的是,该实施例提供的建模方法不仅适用于对图2所示的高空飞行器***的动力学建模,也适用于其它的浮空器与载荷舱通过旋转连接器连接的高空飞行器***,例如,飞艇携带载荷舱,高空气球携带载荷舱等高空飞行器***。此外,载荷舱的推进***也可以是不同数量与不同布局,该实施例提供的建模方法同样可以适用。
需要说明的是,在附图的流程图虽然示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
本申请还提供了一种存储介质的实施例。该存储介质包括存储的程序,其中,在程序运行时控制存储介质所在设备执行本发明实施例的应用于高空飞行器***的建模方法。
本申请还提供了一种处理器的实施例。该处理器用于运行程序,其中,程序运行时执行本发明实施例的应用于高空飞行器***的建模方法。
本申请还提供了一种应用于高空飞行器***的建模装置的实施例。
图2是根据本发明实施例的一种可选的应用于高空飞行器***的建模装置的示意图,如图2所示,该装置包括逻辑单元10,确定单元20,建模单元30和仿真单元40,其中:
逻辑单元,用于将高空飞行器***的物理模型简化为双刚体***,其中,浮空器和吊舱均被简化为一个刚体;确定单元,用于确定高空飞行器***的广义坐标,其中,广义坐标至少包括高空飞行器***中每个刚体的姿态角信息;建模单元,用于求解高空飞行器***的拉格朗日函数并建立拉格朗日方程;仿真单元,用于依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息。
该实施例通过逻辑单元将高空飞行器***的物理模型简化为双刚体***,确定单元确定高空飞行器***的广义坐标,建模单元求解高空飞行器***的拉格朗日函数并建立拉格朗日方程,仿真单元依据建立的拉格朗日方程确定高空飞行器***中每个刚体的姿态角信息,解决了相关技术中对高空飞行器***进行建模仿真的结果相对于实际情况的误差较大的技术问题,进而实现了能够减小对高空飞行器***进行建模仿真的误差的技术效果。
可选地,确定单元包括:第一建立模块,用于建立高空飞行器***中每个刚体的位置坐标系和姿态角坐标系;第一确定模块,用于确定广义坐标,其中,广义坐标包括每个刚体的姿态角和高空飞行器***的位置。
可选地,建模单元包括:求解模块,用于求解高空飞行器***中每个刚体的动能矩阵和势能矩阵;第二确定模块,用于根据高空飞行器***中每个刚体的动能矩阵和势能矩阵确定高空飞行器的拉格朗日函数矩阵。
可选地,求解模块包括:确定子模块,用于基于广义坐标确定高空飞行器***中每个刚体的广义速度矩阵,其中,广义速度矩阵为通过广义坐标的求导式表示的刚体的速度和角速度的矩阵;逻辑子模块,用于通过广义速度矩阵表示高空飞行器***中每个刚体的动能矩阵。
可选地,建模单元包括:获取模块,用于获取高空飞行器的广义受力向量,其中,广义受力向量为高空飞行器***所受到的力的集合通过广义坐标表示的向量;第二建立模块,用于根据拉格朗日函数矩阵和广义受力向量建立拉格朗日方程。
可选地,仿真单元包括:运算模块,用于利用矩阵乘法运算规则执行拉格朗日方程中的求导运算。
上述的装置可以包括处理器和存储器,上述单元均可以作为程序单元存储在存储器中,由处理器执行存储在存储器中的上述程序单元来实现相应的功能。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flashRAM),存储器包括至少一个存储芯片。
上述本申请实施例的顺序不代表实施例的优劣。
在本申请的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。在本申请所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。
其中,以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个***,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,单元或模块的间接耦合或通信连接,可以是电性或其它的形式。
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本申请的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本申请原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本申请的保护范围。
Claims (12)
1.一种应用于高空飞行器***的建模方法,其特征在于,所述高空飞行器***包括柔性连接的浮空器和吊舱,所述吊舱上设置有助推器,所述方法包括:
将所述高空飞行器***的物理模型简化为双刚体***,其中,所述浮空器和所述吊舱均被简化为一个刚体;
确定所述高空飞行器***的广义坐标,其中,所述广义坐标至少包括所述高空飞行器***中每个刚体的姿态角信息;
求解所述高空飞行器***的拉格朗日函数并建立拉格朗日方程;
依据建立的所述拉格朗日方程确定所述高空飞行器***中每个刚体的姿态角信息;
其中,所述建立拉格朗日方程,通过以下步骤进行:
获取所述高空飞行器的广义受力向量,其中,所述广义受力向量为所述高空飞行器***所受到的力的集合通过所述广义坐标表示的向量,所述高空飞行器***所受到的力至少包括:主动控制力和非主动控制力;所述主动控制力为所述高空飞行器***上的助推器产生的推力;
根据拉格朗日函数矩阵和所述广义受力向量建立所述拉格朗日方程。
2.根据权利要求1所述的方法,其特征在于,确定所述高空飞行器***的广义坐标包括:
建立所述高空飞行器***中每个刚体的位置坐标系和姿态角坐标系;
确定所述广义坐标,其中,所述广义坐标包括每个刚体的姿态角和所述高空飞行器***的位置。
3.根据权利要求2述的方法,其特征在于,求解所述高空飞行器***的拉格朗日函数包括:
求解所述高空飞行器***中每个刚体的动能矩阵和势能矩阵;
根据所述高空飞行器***中每个刚体的动能矩阵和势能矩阵确定所述高空飞行器的拉格朗日函数矩阵。
4.根据权利要求3所述的方法,其特征在于,求解所述高空飞行器***中每个刚体的动能矩阵和势能矩阵包括:
基于所述广义坐标确定所述高空飞行器***中每个刚体的广义速度矩阵,其中,所述广义速度矩阵为通过所述广义坐标的求导式表示的刚体的速度和角速度的矩阵;
通过所述广义速度矩阵表示所述高空飞行器***中每个刚体的动能矩阵。
5.根据权利要求1所述的方法,其特征在于,依据建立的所述拉格朗日方程确定所述高空飞行器***中每个刚体的姿态角信息包括:
利用矩阵乘法运算规则执行所述拉格朗日方程中的求导运算。
6.一种应用于高空飞行器***的建模装置,其特征在于,所述高空飞行器***包括柔性连接的浮空器和吊舱,所述吊舱上设置有助推器,所述装置包括:
逻辑单元,用于将所述高空飞行器***的物理模型简化为双刚体***,其中,所述浮空器和所述吊舱均被简化为一个刚体;
确定单元,用于确定所述高空飞行器***的广义坐标,其中,所述广义坐标至少包括所述高空飞行器***中每个刚体的姿态角信息;
建模单元,用于求解所述高空飞行器***的拉格朗日函数并建立拉格朗日方程;
仿真单元,用于依据建立的所述拉格朗日方程确定所述高空飞行器***中每个刚体的姿态角信息;
所述建模单元包括:
获取模块,用于获取所述高空飞行器的广义受力向量,其中,所述广义受力向量为所述高空飞行器***所受到的力的集合通过所述广义坐标表示的向量;所述高空飞行器***所受到的力至少包括:主动控制力和非主动控制力;所述主动控制力为所述高空飞行器***上的助推器产生的推力;
第二建立模块,用于根据拉格朗日函数矩阵和所述广义受力向量建立所述拉格朗日方程。
7.根据权利要求6所述的装置,其特征在于,所述确定单元包括:
第一建立模块,用于建立所述高空飞行器***中每个刚体的位置坐标系和姿态角坐标系;
第一确定模块,用于确定所述广义坐标,其中,所述广义坐标包括每个刚体的姿态角和所述高空飞行器***的位置。
8.根据权利要求7述的装置,其特征在于,所述建模单元包括:
求解模块,用于求解所述高空飞行器***中每个刚体的动能矩阵和势能矩阵;
第二确定模块,用于根据所述高空飞行器***中每个刚体的动能矩阵和势能矩阵确定所述高空飞行器的拉格朗日函数矩阵。
9.根据权利要求8所述的装置,其特征在于,所述求解模块包括:
确定子模块,用于基于所述广义坐标确定所述高空飞行器***中每个刚体的广义速度矩阵,其中,所述广义速度矩阵为通过所述广义坐标的求导式表示的刚体的速度和角速度的矩阵;
逻辑子模块,用于通过所述广义速度矩阵表示所述高空飞行器***中每个刚体的动能矩阵。
10.根据权利要求6所述的装置,其特征在于,所述仿真单元包括:
运算模块,用于利用矩阵乘法运算规则执行所述拉格朗日方程中的求导运算。
11.一种存储介质,其特征在于,所述存储介质包括存储的程序,其中,在所述程序运行时控制所述存储介质所在设备执行权利要求1至5中任意一项所述的应用于高空飞行器***的建模方法。
12.一种处理器,其特征在于,所述处理器用于运行程序,其中,所述程序运行时执行权利要求1至5中任意一项所述的应用于高空飞行器***的建模方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710267450.0A CN108733858B (zh) | 2017-04-21 | 2017-04-21 | 应用于高空飞行器***的建模方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710267450.0A CN108733858B (zh) | 2017-04-21 | 2017-04-21 | 应用于高空飞行器***的建模方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108733858A CN108733858A (zh) | 2018-11-02 |
CN108733858B true CN108733858B (zh) | 2023-08-11 |
Family
ID=63933638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710267450.0A Active CN108733858B (zh) | 2017-04-21 | 2017-04-21 | 应用于高空飞行器***的建模方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108733858B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111413995B (zh) * | 2020-03-24 | 2021-03-05 | 北京科技大学 | 双刚体特征点间相对位置跟踪与姿态同步控制方法及*** |
CN112052516B (zh) * | 2020-08-13 | 2021-11-23 | 中国人民解放军军事科学院国防科技创新研究院 | 基于序列摆放的组件布局随机采样方法 |
CN117588358A (zh) * | 2023-12-26 | 2024-02-23 | 中国电力工程顾问集团有限公司 | 基于混合状态误差的高空风力发电装置防碰撞方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4795113A (en) * | 1984-02-06 | 1989-01-03 | Minovitch Michael Andrew | Electromagnetic transportation system for manned space travel |
RU2009072C1 (ru) * | 1989-12-14 | 1994-03-15 | Анатолий Александрович Снегов | Дирижабль |
US9146557B1 (en) * | 2014-04-23 | 2015-09-29 | King Fahd University Of Petroleum And Minerals | Adaptive control method for unmanned vehicle with slung load |
CN105700542A (zh) * | 2016-03-30 | 2016-06-22 | 北京航空航天大学 | 一种基于矢量场制导和最小二乘法的平流层飞艇控制分配方法 |
CN106569502A (zh) * | 2016-05-20 | 2017-04-19 | 上海铸天智能科技有限公司 | 多旋翼飞行器抓捕目标后的复合体姿态自适应控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10114911B2 (en) * | 2010-05-24 | 2018-10-30 | Fujitsu Limited | Fluid structure interaction simulation method and apparatus, and computer-readable storage medium |
-
2017
- 2017-04-21 CN CN201710267450.0A patent/CN108733858B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4795113A (en) * | 1984-02-06 | 1989-01-03 | Minovitch Michael Andrew | Electromagnetic transportation system for manned space travel |
RU2009072C1 (ru) * | 1989-12-14 | 1994-03-15 | Анатолий Александрович Снегов | Дирижабль |
US9146557B1 (en) * | 2014-04-23 | 2015-09-29 | King Fahd University Of Petroleum And Minerals | Adaptive control method for unmanned vehicle with slung load |
CN105700542A (zh) * | 2016-03-30 | 2016-06-22 | 北京航空航天大学 | 一种基于矢量场制导和最小二乘法的平流层飞艇控制分配方法 |
CN106569502A (zh) * | 2016-05-20 | 2017-04-19 | 上海铸天智能科技有限公司 | 多旋翼飞行器抓捕目标后的复合体姿态自适应控制方法 |
Non-Patent Citations (1)
Title |
---|
平流层飞艇运动分析与仿真;胡国昌;吴美平;;哈尔滨工程大学学报(第11期);第1501-1508页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108733858A (zh) | 2018-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Koehl et al. | Aerodynamic modelling and experimental identification of a coaxial-rotor UAV | |
Zhang et al. | Modeling and flight control simulation of a quadrotor tailsitter vtol uav | |
Peloquin et al. | Design of a passive vertical takeoff and landing aquatic UAV | |
Lu et al. | Real-time simulation system for UAV based on Matlab/Simulink | |
CN108733858B (zh) | 应用于高空飞行器***的建模方法和装置 | |
Pshikhopov et al. | Mathematical model of robot on base of airship | |
CN111976974A (zh) | 一种飞行控制方法、无人机及存储介质 | |
Jeong et al. | Dynamic modeling and analysis of a single tilt-wing unmanned aerial vehicle | |
Caetano et al. | Rigid-body kinematics versus flapping kinematics of a flapping wing micro air vehicle | |
CN116360497B (zh) | 一种基于强化学***流层飞艇轨迹跟踪方法 | |
Ahranjani et al. | Applied flight dynamics modeling and stability analysis of a nonlinear time-periodic mono-wing aerial vehicle | |
Kemper et al. | Impact of center of gravity in quadrotor helicopter controller design | |
Husseyin et al. | Design considerations for a stopped-rotor cyclocopter for Venus exploration | |
Montalvo et al. | Avoiding lockout instability for towed parafoil systems | |
Gavrilovic et al. | Performance improvement of small UAVs through energy-harvesting within atmospheric gusts | |
Yoder et al. | Aerodynamic Tethered Sails for Scientific Balloon Trajectory Control: Small-Scale Experimental Demonstration | |
Battipede et al. | Flight simulator for the control law design of an innovative remotely-piloted airship | |
Zhao et al. | Design and dynamic modelling of a coaxial-rotor system | |
Merz et al. | Feasibility study of a circularly towed cable-body system for uav applications | |
Hang et al. | Dynamic Modeling and Longitudinal Stability Analysis of a Stratospheric Dual-Hull Airship | |
Pan et al. | Design and simulation of high altitude air-launched automatic underwater vehicles | |
Moller | Automated landing of a quadrotor unmanned aerial vehicle on a translating platform | |
Zhang et al. | Numerical simulation of postlaunching behaviors for a “balloon-borne UAV system” | |
Razzak et al. | Computational multiphysics platform for virtual controlled flight of quadrotor unmanned aerial vehicles | |
Coulter et al. | Evaluation and Comparison of Sailing Flight Optimization Algorithms for a Stratospheric Dual Aircraft Platform Concept |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231211 Address after: 518000 room 301-b077, building 2, No.1, Mawu Road, Baoan community, Yuanshan street, Longgang District, Shenzhen City, Guangdong Province Patentee after: SHENZHEN KUANG-CHI SPACE TECH. Co.,Ltd. Address before: 570100 building a, business incubation center, Haikou national high tech Zone, no.6, Yaogu 1st Road, Nanhai Avenue, Haikou City, Hainan Province Patentee before: HAIKOU INSTITUTE OF FUTURE TECHNOLOGY |