CN112591132A - 一种飞机舵面转动惯量坐标变换方法 - Google Patents
一种飞机舵面转动惯量坐标变换方法 Download PDFInfo
- Publication number
- CN112591132A CN112591132A CN202011549356.2A CN202011549356A CN112591132A CN 112591132 A CN112591132 A CN 112591132A CN 202011549356 A CN202011549356 A CN 202011549356A CN 112591132 A CN112591132 A CN 112591132A
- Authority
- CN
- China
- Prior art keywords
- rotating shaft
- control surface
- rotational inertia
- airplane
- points
- 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 14
- 230000001131 transforming effect Effects 0.000 title claims abstract description 13
- 230000005484 gravity Effects 0.000 claims abstract description 19
- 238000013507 mapping Methods 0.000 abstract description 4
- 238000003672 processing method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Toys (AREA)
Abstract
本发明涉及一种飞机舵面转动惯量坐标变换的处理方法,属于飞机重量工程领域。一种飞机舵面转动惯量坐标变换方法,包括以下步骤:(1)简化转轴与原飞机坐标系的关系,在转轴上任意取两个点,基于两点三维坐标值,求得两点间直线段与各坐标系的投影长度,进而求得直线段与各投影的夹角;(2)构建重心点距转轴的关系图,计算舵面重心至转轴轴线的空间距离;(3)输入舵面在飞机坐标系下的重量、重心和转动惯量参数,以及所需转化的轴线上两个点坐标值,即可计算出绕转轴的转动惯量。本发明具有如下优点:用计算三角函数来代替作图法,解决了三维空间作图难以判读并容易出错的问题。
Description
技术领域
本发明涉及一种飞机舵面转动惯量坐标变换的处理方法,属于飞机重量工程领域。
背景技术
转动惯量是舵面作转动运动时其惯性的度量。转动惯量的大小取决于舵面的质量、质量分布和转动运动时所对应的转动轴位置。舵面转动惯量是飞机载荷计算,操稳特性和气动弹性分析不可缺少的原始数据。
常规的舵面转动惯量一般是平行于飞机自身坐标系的,但有时会需要取得各舵面的绕转轴的转动惯量,此时需对转动惯量进行坐标变换。《飞机设计手册第8册:重量平衡与控制》提供了转动惯量坐标变换的处理方法,如图1所示用作图法取得坐标轴旋转角度,再用公式进行坐标轴旋转变换,将转动惯量变换至与所需转轴方向一致的坐标系,计算公式如下:
IA=Ix cos2θa+Iy cos2γa+Iz cos2βa-2Ixy cosθa cosγa-2Ixz cosθa cosβa-2Iyzcosγa cosβa
IB=Ix cos2θb+Iy cos2γb+Iz cos2βb-2Ixy cosθb cosγb-2Ixz cosθb cosβb-2Iyzcosγb cosβb
IC=Ix cos2θc+Iy cos2γc+Iz cos2βc-2Ixy cosθc cosγc-2Ixz cosθc cosβc-2Iyzcosγc cosβc
IAB=Ix(-cosθa cosθb)+Iy(-cosγa cosγb)+Iz(-cosβa cosβb)+Ixy(cosθa cosγb+cosθb cosγa)+Iyz(cosγa cosβb+cosγb cosβa)+Ixz(cosβb cosθa+cosβa cosθb)
;最后进行如图2所示转动惯量平行移轴变换,得到绕转轴的转动惯量,计算公式如下:I′zm=Izm+M(Δx)2。但作图取得转轴与飞机坐标系的旋转夹角,以及重心位置与转轴的距离等参数,在三维空间结构中非常复杂,常常难以读取参数,并容易出错。
发明内容
本发明的目的是解决现有技术中存在的问题,提供一种飞机舵面转动惯量坐标变换方法,从而解决三维空间作图难以判读并容易出错的问题。
为了实现上述目的,本发明采用如下技术方案:一种飞机舵面转动惯量坐标变换方法,包括以下步骤:
(1)简化转轴与原飞机坐标系的关系,在转轴上任意取两个点,基于两点三维坐标值,求得两点间直线段与各坐标系的投影长度,进而求得直线段与各投影的夹角;
(2)构建重心点距转轴的关系图,计算舵面重心至转轴轴线的空间距离;
(3)输入舵面在飞机坐标系下的重量、重心和转动惯量参数,以及所需转化的轴线上两个点坐标值,即可计算出绕转轴的转动惯量。
优选的,利用反余弦公式求得直线段与各投影的夹角。
优选的,利用三角函数公式计算舵面重心至转轴轴线的空间距离。
与现有技术相比,本发明具有如下优点:用计算三角函数来代替作图法,解决了三维空间作图难以判读并容易出错的问题。
附图说明
图1是本发明背景技术飞机设计手册中的转动惯量坐标变换及计算公式;
图2是本发明背景技术飞机设计手册中的转动惯量平行移轴及计算公式;
图3是本发明实施例简化后的舵面转轴轴线与重心坐标系各夹角的关系图;
图4是本发明实施例舵面重心与转轴轴线的距离关系图。
具体实施方式
下面结合附图3-4对本发明做进一步详述:一种飞机舵面转动惯量坐标变换方法,包括以下步骤:
首先,采取如图3所示简化转轴与原飞机坐标系的关系,通过在转轴上任意取两个点,一个设置为“起点”,另一个为“终点”,然后基于两点三维坐标值,可求得该直线段与各坐标系的投影长度,再利用反余弦公式求得直线段与各投影的夹角,从而可直接计算转轴与飞机坐标系的各夹角;
然后,如图4所示构建重心点距转轴的关系图,利用三角函数公式计算舵面重心至转轴轴线的空间距离;
最后,通过输入舵面在飞机坐标系下的重量、重心和转动惯量参数,以及所需转化的轴线上两个点坐标值,即可计算出绕转轴的转动惯量。
以某舵面在飞机坐标系下的重量W,重心(x,y,z)和转动惯量(Ix,Iy,Iz,Ixy)值为例,对本发明一种飞机舵面转动惯量坐标变换方法的做进一步详细说明。
首先,采取如图3所示简化转轴与原飞机坐标系的关系,通过在转轴上任意取两个点,坐标值分别为(x1,y1,z1)和(x2,y2,z2),然后基于两点三维坐标值,可求得该直线段与各坐标系的投影长度L0,再利用反余弦公式求得直线段与各投影的夹角α、β、γ,即为转轴与飞机坐标系的各夹角。
然后,如图4所示构建重心点距转轴的关系图,利用三角函数公式计算舵面重心至转轴轴线的空间距离h,公式如下:
a=L0
最后,通过如下公式所示的将转动惯量换算至舵面转轴方向一致轴线上的转动惯量I1,再将I1进行平移至与舵面转轴重合,得到绕转轴的转动惯量I2即为所求。
I1=Ix·cos2α+Iy·cos2β+Iz·cos2γ-Ixy·cosα·cosβ
I2=I1+W·h2
上述实施例仅仅是本发明的优选实施方式,不构成对本发明的限制,本领域普通技术人员在不脱离本发明本质的基础上,所做的任何引申、变换、等同替换等均应包括在本发明的保护范围内。
Claims (3)
1.一种飞机舵面转动惯量坐标变换方法,其特征在于:包括以下步骤:
(1)简化转轴与原飞机坐标系的关系,在转轴上任意取两个点,基于两点三维坐标值,求得两点间直线段与各坐标系的投影长度,进而求得直线段与各投影的夹角;
(2)构建重心点距转轴的关系图,计算舵面重心至转轴轴线的空间距离;
(3)输入舵面在飞机坐标系下的重量、重心和转动惯量参数,以及所需转化的轴线上两个点坐标值,即可计算出绕转轴的转动惯量。
2.根据权利要求1所述的飞机舵面转动惯量坐标变换方法,其特征在于:在步骤(1)中,利用反余弦公式求得直线段与各投影的夹角。
3.根据权利要求1所述的飞机舵面转动惯量坐标变换方法,其特征在于:在步骤(2)中,利用三角函数公式计算舵面重心至转轴轴线的空间距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011549356.2A CN112591132A (zh) | 2020-12-24 | 2020-12-24 | 一种飞机舵面转动惯量坐标变换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011549356.2A CN112591132A (zh) | 2020-12-24 | 2020-12-24 | 一种飞机舵面转动惯量坐标变换方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112591132A true CN112591132A (zh) | 2021-04-02 |
Family
ID=75200981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011549356.2A Pending CN112591132A (zh) | 2020-12-24 | 2020-12-24 | 一种飞机舵面转动惯量坐标变换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112591132A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113295134A (zh) * | 2021-06-08 | 2021-08-24 | 北京普利永华科技发展有限公司 | 转轴无关性飞机舵面转角测量方法 |
CN114199457A (zh) * | 2022-02-17 | 2022-03-18 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种500公斤级飞机整机重量特性高精度测量方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846573A (zh) * | 2010-05-19 | 2010-09-29 | 中国科学院电工研究所 | 一种测试不完整球形超导转子转动惯量的方法 |
CN102354333A (zh) * | 2011-09-30 | 2012-02-15 | 中国海洋大学 | 一种钢悬链式立管出平面运动的分析方法 |
CN104229157A (zh) * | 2014-08-26 | 2014-12-24 | 中国直升机设计研究所 | 一种直升机极限重心配载方法 |
CN105527955A (zh) * | 2014-09-28 | 2016-04-27 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机质量特性建模方法 |
CN106777689A (zh) * | 2016-12-15 | 2017-05-31 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于有限元模型的飞机双铰链舵面偏转方法 |
US20170322101A1 (en) * | 2016-05-09 | 2017-11-09 | Fong Bong Enterprise Co., Ltd. | Calibrating Device for Measuring and Calibrating the Center of Gravity of a Remote Control Aircraft or an Airfoil Thereof |
CN111044220A (zh) * | 2019-11-26 | 2020-04-21 | 北京卫星制造厂有限公司 | 无人机质心惯量集成测试方法 |
-
2020
- 2020-12-24 CN CN202011549356.2A patent/CN112591132A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846573A (zh) * | 2010-05-19 | 2010-09-29 | 中国科学院电工研究所 | 一种测试不完整球形超导转子转动惯量的方法 |
CN102354333A (zh) * | 2011-09-30 | 2012-02-15 | 中国海洋大学 | 一种钢悬链式立管出平面运动的分析方法 |
CN104229157A (zh) * | 2014-08-26 | 2014-12-24 | 中国直升机设计研究所 | 一种直升机极限重心配载方法 |
CN105527955A (zh) * | 2014-09-28 | 2016-04-27 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机质量特性建模方法 |
US20170322101A1 (en) * | 2016-05-09 | 2017-11-09 | Fong Bong Enterprise Co., Ltd. | Calibrating Device for Measuring and Calibrating the Center of Gravity of a Remote Control Aircraft or an Airfoil Thereof |
CN106777689A (zh) * | 2016-12-15 | 2017-05-31 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于有限元模型的飞机双铰链舵面偏转方法 |
CN111044220A (zh) * | 2019-11-26 | 2020-04-21 | 北京卫星制造厂有限公司 | 无人机质心惯量集成测试方法 |
Non-Patent Citations (2)
Title |
---|
张家盛等: "飞机舵面液压伺服作动系统性能分析", 《航空科学技术》 * |
查佳韵等: "基于扭杆摆的汽车动力总成惯量参数测量实现", 《中国测试》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113295134A (zh) * | 2021-06-08 | 2021-08-24 | 北京普利永华科技发展有限公司 | 转轴无关性飞机舵面转角测量方法 |
CN114199457A (zh) * | 2022-02-17 | 2022-03-18 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种500公斤级飞机整机重量特性高精度测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112591132A (zh) | 一种飞机舵面转动惯量坐标变换方法 | |
CN103591949B (zh) | 三轴姿态测量系统非正交性误差的正交补偿方法 | |
CN103411589B (zh) | 一种基于四维实数矩阵的三维图像匹配导航方法 | |
CN109781059A (zh) | 星载点波束天线对地指向精度评估系统 | |
CN107038275B (zh) | 一种机械臂误差分析方法 | |
CN101944240A (zh) | 多机器人三维几何地图的融合方法 | |
CN111623768A (zh) | 一种基于克雷洛夫角奇异条件下的姿态角解算方法 | |
CN111123280A (zh) | 激光雷达的定位方法、装置、系统、电子设备及存储介质 | |
CN110134134A (zh) | 一种无人机悬停状态下的测风方法 | |
CN105043414A (zh) | 一种三轴惯性稳定平台系统的台体控制参数计算方法 | |
Xu et al. | A robust incremental-quaternion-based angle and axis estimation algorithm of a single-axis rotation using MARG sensors | |
CN106353033B (zh) | 一种航空发动机质心计算方法 | |
CN110888123A (zh) | 一种基于旋转矩阵的雷达坐标转换方法 | |
CN110954080A (zh) | 一种消除载体磁干扰的磁罗盘校准方法 | |
CN113084791A (zh) | 机械臂控制方法、机械臂控制装置及终端设备 | |
CN114329943B (zh) | 基于姿态旋转矩阵的控制性能边界设计方法、装置及介质 | |
CN114119684B (zh) | 基于四面体结构的标记点配准方法 | |
CN105043438A (zh) | 空间物体多维参数测量装置 | |
CN109737939B (zh) | 一种实现电子罗盘全方位安装的方法及装置 | |
Wang et al. | Kinematics of coupler curves for spherical four-bar linkages based on new spherical adjoint approach | |
CN111220114B (zh) | 一种单轴旋转载体转轴转角惯性测量系统及方法 | |
CN108253931B (zh) | 一种双目立体视觉测距方法及其测距装置 | |
Mityushov et al. | Quaternionic description of the Euler trajectories of a rigid body in various configuration spaces | |
Ge et al. | Computational geometry and motion approximation | |
CN103075994B (zh) | 一种基于螺旋线法确定回转体目标横滚角的测量方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210402 |
|
RJ01 | Rejection of invention patent application after publication |