CN108680336B - 一种低速风洞试验中直升机旋翼振动位移测量方法 - Google Patents
一种低速风洞试验中直升机旋翼振动位移测量方法 Download PDFInfo
- Publication number
- CN108680336B CN108680336B CN201810728394.0A CN201810728394A CN108680336B CN 108680336 B CN108680336 B CN 108680336B CN 201810728394 A CN201810728394 A CN 201810728394A CN 108680336 B CN108680336 B CN 108680336B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- coding mark
- rotor
- vibration displacement
- mark 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.)
- Active
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 26
- 238000012360 testing method Methods 0.000 title claims abstract description 16
- 238000000691 measurement method Methods 0.000 title claims abstract description 8
- 230000001052 transient effect Effects 0.000 claims abstract description 14
- 230000003068 static effect Effects 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 6
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 10
- 230000010355 oscillation Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 239000003550 marker Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/022—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness by means of tv-camera scanning
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开一种低速风洞试验中直升机旋翼振动位移测量方法,属于低速风洞试验测量技术领域。包括以下步骤:在直升机旋翼桨毂上布置一定数量的编码标记点;获取直升机旋翼静止状态下的图像数据;计算图像数据中编码标记点的三维坐标,并构建旋翼坐标系;获取直升机旋翼高速旋转状态下的瞬态图像数据;计算桨毂中心编码标记点在旋翼坐标系中的三维坐标;计算得到直升机旋翼高速旋转状态下的瞬态振动位移向量。本发明所提供的方法可有效测量直升机旋翼高速旋转造成的振动位移偏差,具有操作简单、非接触式、精度高、安全性好等优点。
Description
技术领域
本发明涉及低速风洞试验测量技术领域,具体指一种低速风洞试验中直升机旋翼振动位移测量方法。
背景技术
低速风洞试验是测量直升机旋翼性能参数,研究旋翼空气动力特性的重要手段。风洞试验中,通过采用各种手段测量旋翼的各项性能参数,如挥舞角、摆振角、扭转变形量等,可以确定旋翼的运动范围、变形状况以及桨尖运动轨迹,从而为旋翼结构设计、桨毂上下限动块设计、桨叶载荷设计、旋翼共锥度调整等提供重要依据。
然而,由于动力源和结构的特点,旋翼在高速旋转过程中会不可避免地发生振动,目前采用的振动控制方法也只能将振动减小到最小水平,无法实现完全消除。并且,在风洞试验中,旋翼高速旋转还会引起整个旋翼试验台发生振动。因此,在对旋翼性能参数进行测量的同时,有必要同步测量旋翼高速旋转造成的振动位移偏差,以对测量得到的参数数据进行修正。但到目前为止,还没有一种方法能够对直升机旋翼高速旋转造成的振动位移进行有效测量,因此,如何测量直升机旋翼振动位移成为了一个亟需解决的问题。
发明内容
本发明提供了一种低速风洞试验中直升机旋翼振动位移测量方法,旨在克服现有技术的不足,其能有效测量直升机旋翼高速旋转造成的振动位移偏差,从而能够对旋翼的各项性能参数数据进行有效修正。
为解决上述技术问题,本发明提出的技术方案是:
一种低速风洞试验中直升机旋翼振动位移测量方法,其特征在于,包括以下步骤:
步骤1:在直升机旋翼的桨毂上布置一定数量的编码标记点;
所述步骤1中桨毂上编码标记点布置于桨毂中心及其周围,编码标记点数量N≥4;
步骤2:通过同步控制器同步触发两台高速CCD相机拍摄获取直升机旋翼静止状态下的图像数据;
步骤3:对步骤2中获取的图像数据进行编码标记点识别与匹配,计算编码标记点在摄像机坐标系中的三维坐标;
步骤4:基于步骤3求得的静止状态下编码标记点三维坐标,构建旋翼坐标系;
所述步骤4中旋翼坐标系构建方法具体如下:
步骤4-1:对静止状态下编码标记点的三维坐标进行最小二乘平面拟合,得到拟合平面S;
步骤4-2:以静止状态下桨毂中心编码标记点为原点、以拟合平面S的法向量为z轴方向,构建旋翼坐标系,计算得到从摄像机坐标系转换到旋翼坐标系的平移向量T和旋转矩阵R;
所述步骤4-2中平移向量T计算如下:
其中,表示静止状态下桨毂中心编码标记点在摄像机坐标系中的三维坐标;
所述步骤4-2中旋转矩阵R计算如下:
其中,(a,b,c)表示拟合平面S的法向量,表示旋转角度;
步骤5:保持两台高速CCD相机不动,采用高频激光器在纳秒级时间内照亮测量区域,通过同步控制器同步触发两台高速CCD相机拍摄获取直升机旋翼高速旋转状态下的瞬态图像数据;
步骤6:对步骤5中获取的瞬态图像数据进行编码标记点识别与匹配,计算桨毂中心编码标记点在摄像机坐标系中的三维坐标,并转换到旋翼坐标系中;
步骤7:计算得到直升机旋翼高速旋转状态下的瞬态振动位移向量d;
所述步骤7中瞬态振动位移向量d计算如下:
其中,表示高速旋转状态下桨毂中心编码标记点在旋翼坐标系中的三维坐标。
本发明有益效果:
本发明所提供的方法能够测量得到直升机旋翼高速旋转造成的振动位移偏差,具有操作简单、非接触式、精度高、安全性好等优点。测量获得的振动位移向量被转换到了旋翼坐标系中,去除了与成像系统的关联,具有数据独立性,因而能够对采用其他手段同步测量的旋翼性能参数进行修正,从而实现旋翼性能参数的精确测量,为直升机旋翼的空气动力特性分析提供精准的参数依据。
附图说明
图1是本发明方法的整体流程图;
图2是直升机旋翼桨毂编码标记点布局示意图;
图3是直升机旋翼静止状态下图像数据获取示意图;
图4是直升机旋翼高速旋转状态下瞬态图像数据获取示意图;
图5是直升机旋翼高速旋转状态下瞬态振动位移示意图。
其中,1为桨毂,2为高速CCD相机,3为同步控制器,4为计算服务器,5为高频激光器。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本实施例提供了一种低速风洞试验中直升机旋翼振动位移测量方法。
步骤1:如图2所示,在直升机旋翼的桨毂1上布置5个编码标记点,分别位于桨毂中心以及与桨毂中心等距的上、下、左、右位置。
步骤2:如图3所示,将两台高速CCD相机2布置在直升机旋翼的上方,在普通光照条件下,通过计算服务器4控制同步控制器3同步触发两台高速CCD相机2拍摄获取直升机旋翼静止状态下的图像数据。
步骤3:对获取的图像数据进行编码标记点识别与匹配,计算编码标记点在摄像机坐标系中的三维坐标如下:
其中,表示静止状态下第i个编码标记点在左图像数据中的图像坐标,表示静止状态下第i个编码标记点在右图像数据中的图像横坐标,Q表示双目立体视觉系统的重投影矩阵,和表示求得的静止状态下第i个编码标记点在摄像机坐标系中的三维齐次坐标和三维坐标。
步骤4:基于静止状态下编码标记点的三维坐标,构建旋翼坐标系,具体方法如下:
步骤4-1:对静止状态下编码标记点三维坐标进行最小二乘平面拟合,得到拟合平面S:
ax+by+cz+d=0 (2)
其中,a,b,c,d表示拟合得到的平面参数。
步骤4-2:以静止状态下桨毂中心编码标记点为原点、以拟合平面S的法向量为z轴方向,构建旋翼坐标系,计算得到从摄像机坐标系转换到旋翼坐标系的平移向量T和旋转矩阵R:
其中,表示静止状态下桨毂中心编码标记点在摄像机坐标系中的三维坐标,(a,b,c)表示拟合平面S的法向量,表示旋转角度。
步骤5:如图4所示,保持两台高速CCD相机2不动,将高频激光器5布置在直升机旋翼的上方,通过计算服务器4控制高频激光器5在纳秒级时间内照亮测量区域,并同步触发两台高速CCD相机2拍摄获取直升机旋翼高速旋转状态下的瞬态图像数据。
步骤6:对获取的图像数据进行编码标记点识别与匹配,计算桨毂中心编码标记点在摄像机坐标系中的三维坐标如下:
其中,表示高速旋转状态下桨毂中心编码标记点在左图像数据中的图像坐标,表示高速旋转状态下桨毂中心编码标记点在右图像数据中的图像横坐标,Q表示双目立体视觉系统的重投影矩阵,和表示求得的高速旋转状态下桨毂中心编码标记点在摄像机坐标系中的三维齐次坐标和三维坐标。
将桨毂中心编码标记点三维坐标从摄像机坐标系转换到旋翼坐标系中:
其中,和表示高速旋转状态下桨毂中心编码标记点在摄像机坐标系和旋翼坐标系中的三维坐标,T表示平移向量,R表示旋转矩阵。
步骤7:如图5所示,当直升机旋翼在高速旋转状态下发生振动位移时,桨毂中心编码标记点必然会随之发生同样的位移变化,因此,可以通过测量桨毂中心编码标记点的位移偏差来得到旋翼的位移偏差。桨毂中心编码标记点的瞬态振动位移向量d计算如下:
其中,表示高速旋转状态下桨毂中心编码标记点在旋翼坐标系中的三维坐标。
以上所揭露的仅为本发明一种较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (3)
1.一种低速风洞试验中直升机旋翼振动位移测量方法,其特征在于,包括以下步骤:
步骤1:在直升机旋翼的桨毂上布置一定数量的编码标记点;
步骤2:通过同步控制器同步触发两台高速CCD相机拍摄获取直升机旋翼静止状态下的图像数据;
步骤3:对步骤2中获取的图像数据进行编码标记点识别与匹配,计算编码标记点在摄像机坐标系中的三维坐标;
步骤4:基于步骤3求得的静止状态下编码标记点的三维坐标,构建旋翼坐标系;
所述步骤4中旋翼坐标系构建方法具体如下:
步骤4-1:对静止状态下编码标记点的三维坐标进行最小二乘平面拟合,得到拟合平面S;
步骤4-2:以静止状态下桨毂中心编码标记点为原点、以拟合平面S的法向量为z轴方向,构建旋翼坐标系,计算得到从摄像机坐标系转换到旋翼坐标系的平移向量T和旋转矩阵R;所述步骤4-2中平移向量T计算如下:
其中,表示静止状态下桨毂中心编码标记点在摄像机坐标系中的三维坐标;
所述步骤4-2中旋转矩阵R计算如下:
其中,(a,b,c)表示拟合平面S的法向量,表示旋转角度;
步骤5:保持两台高速CCD相机不动,采用高频激光器在纳秒级时间内照亮测量区域,通过同步控制器同步触发两台高速CCD相机拍摄获取直升机旋翼高速旋转状态下的瞬态图像数据;
步骤6:对步骤5中获取的瞬态图像数据进行编码标记点识别与匹配,计算桨毂中心编码标记点在摄像机坐标系中的三维坐标,并转换到旋翼坐标系中;
步骤7:计算得到直升机旋翼高速旋转状态下的瞬态振动位移向量d。
2.根据权利要求1所述的低速风洞试验中直升机旋翼振动位移测量方法,其特征在于:所述步骤1中桨毂上编码标记点布置于桨毂中心及其周围,编码标记点数量N≥4。
3.根据权利要求1所述的低速风洞试验中直升机旋翼振动位移测量方法,其特征在于:所述步骤7中瞬态振动位移向量d计算如下:
其中,表示高速旋转状态下桨毂中心编码标记点在旋翼坐标系中的三维坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810728394.0A CN108680336B (zh) | 2018-07-05 | 2018-07-05 | 一种低速风洞试验中直升机旋翼振动位移测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810728394.0A CN108680336B (zh) | 2018-07-05 | 2018-07-05 | 一种低速风洞试验中直升机旋翼振动位移测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108680336A CN108680336A (zh) | 2018-10-19 |
CN108680336B true CN108680336B (zh) | 2019-12-03 |
Family
ID=63813308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810728394.0A Active CN108680336B (zh) | 2018-07-05 | 2018-07-05 | 一种低速风洞试验中直升机旋翼振动位移测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108680336B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109612681B (zh) * | 2018-11-16 | 2020-07-07 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种直升机共轴刚性旋翼模型的气动干扰测量方法 |
CN109657989A (zh) * | 2018-12-20 | 2019-04-19 | 南京航空航天大学 | 直升机高速重载输入级健康状态评估方法 |
CN110095089B (zh) * | 2019-03-11 | 2020-10-13 | 新拓三维技术(深圳)有限公司 | 一种飞行器旋转角度的测量方法及系统 |
CN111292375B (zh) * | 2020-02-11 | 2023-05-09 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种基于位置约束的直升机桨叶标记点识别匹配方法 |
CN113534093B (zh) * | 2021-08-13 | 2023-06-27 | 北京环境特性研究所 | 飞机目标的螺旋桨叶片数量反演方法及目标识别方法 |
CN113701644B (zh) * | 2021-09-15 | 2023-02-03 | 南京航空航天大学 | 基于风洞试验高速摄像技术位移测量系统及位移测量方法 |
CN113834625B (zh) * | 2021-10-13 | 2022-02-25 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种飞行器模型表面压力测量方法及系统 |
CN114295290A (zh) * | 2022-01-04 | 2022-04-08 | 北京航空航天大学 | 一种平流层浮空器螺旋桨在线动平衡调节装置 |
CN114397905B (zh) * | 2022-03-25 | 2022-06-03 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种倾转旋翼机倾转过渡风洞飞行实验方法与系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5137353A (en) * | 1991-10-10 | 1992-08-11 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Angular displacement measuring device |
US6644355B1 (en) * | 2002-12-19 | 2003-11-11 | Daimlerchrysler Corporation | Diffusing corner for fluid flow |
CN105784313B (zh) * | 2016-03-01 | 2019-03-26 | 中国航天空气动力技术研究院 | 一种基于压电技术的风洞实验支杆振动抑制系统及方法 |
CN205642788U (zh) * | 2016-05-11 | 2016-10-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 高速风洞大攻角俯仰动态失速试验装置 |
-
2018
- 2018-07-05 CN CN201810728394.0A patent/CN108680336B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108680336A (zh) | 2018-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108680336B (zh) | 一种低速风洞试验中直升机旋翼振动位移测量方法 | |
CN108844524B (zh) | 一种基于振动修正的直升机旋翼桨叶运动参数测量方法 | |
CN103134477B (zh) | 一种基于双目立体视觉的直升机旋翼桨叶运动参数测量方法 | |
CN111028220B (zh) | 一种点云铆钉齐平度自动检测方法 | |
CN204992418U (zh) | 一种无人机输电线路缺陷自动巡检装置 | |
CN104748693B (zh) | 一种基于双目立体视觉的桨叶剖面扭转刚度测量系统 | |
CN106092057B (zh) | 一种基于四目立体视觉的直升机旋翼桨叶动态轨迹测量方法 | |
CN109684709B (zh) | 一种航空发动机叶片修复的自适应定位方法 | |
CN107966112A (zh) | 一种大尺寸旋翼运动参数测量方法 | |
CN105067011A (zh) | 一种基于视觉标定及坐标转换的测量系统整体校准方法 | |
CN102818544B (zh) | 汽车轮毂螺栓孔节圆中心和中孔中心偏心距在线测量方法 | |
CN107869960B (zh) | 一种飞机机翼变形测量装置、安装方法及测试方法 | |
CN108318009B (zh) | 一种基于无人机视频的通讯塔垂直度检测方法 | |
CN108204879B (zh) | 一种转动惯量的光学测量方法及系统 | |
JP7188832B1 (ja) | 無人飛行体、および、風力発電設備の点検システム、並びに、風力発電設備の点検方法 | |
CN110849331B (zh) | 基于三维点云数据库模型的单目视觉测量与地面试验方法 | |
CN107390704A (zh) | 一种基于imu姿态补偿的多旋翼无人机光流悬停方法 | |
CN111977025B (zh) | 一种飞行试验外挂式直升机桨叶运动状态测量装置及方法 | |
CN114993608B (zh) | 一种风洞模型三维姿态角测量方法 | |
CN106774365A (zh) | 一种基于外部视觉定位的无人机调试系统 | |
CN108638497A (zh) | 一种3d打印机打印模型外表面的全方位检测系统和方法 | |
CN113031462A (zh) | 无人机用港机巡检航线规划系统及方法 | |
CN112792814A (zh) | 基于视觉标志的机械臂零点标定方法 | |
CN114062265B (zh) | 一种用于视觉系统的支撑结构稳定性的评估方法 | |
CN114639025B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |