CN108507754B - 抗眩光的风洞试验模型弹性变形视频测量方法 - Google Patents
抗眩光的风洞试验模型弹性变形视频测量方法 Download PDFInfo
- Publication number
- CN108507754B CN108507754B CN201810315108.8A CN201810315108A CN108507754B CN 108507754 B CN108507754 B CN 108507754B CN 201810315108 A CN201810315108 A CN 201810315108A CN 108507754 B CN108507754 B CN 108507754B
- Authority
- CN
- China
- Prior art keywords
- test model
- wind tunnel
- elastic deformation
- tunnel test
- camera
- 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
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/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
-
- 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/08—Aerodynamic models
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种抗眩光的风洞试验模型弹性变形视频测量方法,包括如下步骤:在试验模型上喷涂形成圆形荧光标记点,采用紫光照明,激发荧光标记点发出绿光;用装有绿光带通滤光片的相机采集试验模型运动至指定迎角αi时对应荧光标记点图像G’αi,计算G’αi中各个圆形荧光标记点在风洞坐标系oxyz下的三维坐标,解算各截面的弯曲与扭转变形。与现有技术相比,本发明的积极效果是:既确保试验模型在其整个姿态变化过程中,标记点不受眩光与镜面成像干扰,又确保相机标定光照环境与试验时光照环境的一致性,提升风洞试验模型弹性变形视频测量的精准度和可靠性,具有巨大的工程应用前景。
Description
技术领域
本发明属于基于机器视觉和摄影测量的风洞试验技术领域,具体涉及一种抗眩光的风洞试验模型弹性变形视频测量方法。
背景技术
风洞试验中模型及其支撑系统在气动载荷作用下会发生弹性变形。如2.4米跨声速风洞试验时模型承受的气动载荷高达数吨,即使是高强度钢制的机翼也会发生明显弹性变形,而大量研究表明:转捩、分离以及激波/边界层干扰等复杂的流动现象对形状变化非常敏感,模型形状细微的变化可能导致气动特性产生较大变化。
因此,准确测量试验模型的扭转和弯曲变形,掌握实测气动数据与其试验模型气动外形间的对应关系,是风洞试验数据实现模型弹性影响修正的前提,也是基于试验数据验证CFD数值模拟结果的必然要求。
鉴于视频测量(Videogrammetric Measurement,VM)技术对试验模型设计无特殊要求,仅需在试验模型上粘印标记点,即可利用共线方程求解标记点的三维坐标,进而获得试验模型的扭转和弯曲变形数据,因此,受到国内外风洞试验机构的青睐。
但按照风洞模型设计规范(GJB569-88和GJB180-86),国外常规高速风洞试验模型表面粗糙度一般取Ra=0.4~0.8微米。如图1所示,飞行器试验模型总装后其光泽度一般达到镜面反光水准,导致模型在如图2所示的光照环境下必然出现眩光现象和镜面成像干扰,而且,眩光位置和镜面成像干扰位置随试验模型的姿态改变而变化,当标记点处出现眩光现象,轻则降低标记点三维坐标的视频测量精准度,重则使图像算法无法准确检测出标记点,导致测量失败。
虽然Optotrak光学商用测量系统采用主动发光方式,可解决眩光现象和镜面成像干扰问题,但需要在试验模型上平齐嵌装主动发光MARKER点,既破坏模型的外形,又必须在试验模型上开孔布线为MARKER点供电,导致模型设计与制造困难,且造价高。另一方面,因Optotrak外形尺寸长达1米,3个线阵CCD间距达0.45米,测量物体必须置于离Optotrak1.5米到6米以内,导致现有的高速风洞的观察窗很难具备测量条件。
发明内容
为了克服现有技术的上述缺点,本发明提出了一种抗眩光的风洞试验模型弹性变形视频测量方法,通过在模型表面喷涂圆形荧光材料标记点,利用紫外光激发荧光材料标记点发出绿光,再通过相机镜头前的(绿光带通)滤光片消除紫外光产生的眩光现象,确保试验模型在其整个姿态变化过程中,标记点不受眩光现象影响,提升风洞试验模型弹性变形视频测量的精准度,具有巨大的工程应用前景。
本发明解决其技术问题所采用的技术方案是:一种抗眩光的风洞试验模型弹性变形视频测量方法,包括如下步骤:
步骤一、采用绿色荧光喷漆在试验模型上喷涂形成圆形荧光标记点;
步骤二、在地面实验室,采用紫光照明,并利用在镜头前安装有绿光带通滤光片的相机采集圆形荧光标记点的图像,得到各个圆形荧光标记点在风洞坐标系oxyz下的三维坐标;
步骤三、将试验模型安装在风洞中,再打开风洞试验模型弹性变形视频测量系统的紫光光源,并在风洞试验模型弹性变形视频测量系统的相机镜头前安装绿光带通滤光片;
步骤四、对风洞试验模型弹性变形视频测量系统的相机内外参数进行标定;
步骤五、风洞试验时,依次将试验模型运动至指定迎角αi并利用步骤三的相机采集αi对应荧光标记点图像G’αi;
步骤六、利用步骤四得到的相机内外参数、给定αi时对应的荧光标记点图像G’αi,计算G’αi中标记点的三维坐标;
步骤七、利用步骤六计算得到的荧光标记点的三维坐标,计算各截面的弯曲与扭转变形。
与现有技术相比,本发明的积极效果是:
本发明提出抗眩光的风洞试验模型弹性变形视频测量方法,消除了传统方法因试验模型镜面反射所致的眩光和镜面成像干扰,可确保试验模型在其整个姿态变化过程中,标记点不受眩光和镜面成像干扰,提升风洞试验模型弹性变形视频测量的精准度,具有巨大的工程应用前景。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是风洞试验模型的眩光与镜面成像干扰图;
图2是风洞试验模型弹性变形视频测量示意图;
图3是喷涂的圆形荧光标记点示意图;
图4为相机的结构示意图。
具体实施方式
一种抗眩光的风洞试验模型弹性变形视频测量方法,包括如下内容:
1、如图3所示,采用绿色的荧光喷漆,在试验模型上喷涂圆形荧光标记点,喷涂厚度为10微米,要求所有荧光标记点在图1所示的相机中成像图像,至少包含30个像素;要求每列标记点位于的截面(即用与xoz坐标轴平面平行的平面,沿翼展方向切割机翼而成的截面)与X轴平行;其中标记点位于中性轴对应的翼表面上,用于测量弯曲变形;机翼的前缘的标记点和对应后缘标记点 用于测量扭转变形。
2、在地面实验室,采用紫光照明(一般用365纳米波长的紫光灯),并按照图3方式在相机镜头前安装(绿光带通)滤光片,采集圆形荧光标记点图像;利用成熟的机器视觉和摄影测量方法或商业软件(如德国GOM的摄影测量软件),得到图1和图2所示各个荧光标记点在风洞坐标系oxyz下的三维坐标。
3、在风洞中,将图2所示的紫光光源打开,按照图4方式,在图2的相机镜头前安装(绿光带通)滤光片;将试验模型安装在风洞中。
4、将试验的姿态运动范围按升序排序得到{α1,α2,…,αn},考虑到试验时气动力导致的弹性变形,将姿态运动范围扩展为{α0,α1,α2,…,αn,αn+1},其中α0=α1-2度,αn+1=αn+3度。
5、风洞试验现场的相机内外参数标定。试验前,先启动迎角机构,按照{α0,α1,α2,…,αi,…,αn,αn+1}顺序;依次将试验模型运动至指定迎角αi并控制图1的相机采集αi对应荧光标记点图像Gαi,得到{Gα0,Gα1,Gα2,…,Gαi,…,Gαn,Gαn+1};利用第二步测得的标记点在风洞坐标系oxyz下的三维坐标,计算图2的相机内参数;通过试验模型迎角为零时的荧光标记点图像G0,计算图2的相机的外参数。
6、试验时,依次将试验模型运动至指定迎角αi并控制图1的相机采集αi对应荧光标记点图像G’αi,得到{G’α1,G’α2,…,G’αi,…,G’αn}。
7、利用第五步得到的相机内外参数、给定αi时对应的荧光标记点图像G’αi,利用成熟的机器视觉和摄影测量方法或商业软件(如德国GOM的摄影测量软件),可以得到G’αi中标记点的三维坐标。
8、利用荧光标记点的三维坐标,计算各截面的弯曲与扭转变形。计算第j个截面标记点弯曲与扭转变形的方法如下:
由于机身刚度很大,与机身连接的翼根部的第1个截面基本不发生弯曲与扭转变形,故将其作为计算第j个截面弯曲与扭转变形的基准。设G’αi中第j个截面(即第j列)中性轴上的标记点为第1个截面中性轴上的标记点为则第j个截面弯曲变形
式中表示点的Z坐标,表示点的Z坐标,表示点的Z坐标,表示点的Z坐标。
第j个截面扭转变形为
式中
其中,与为G’αi中第j列荧光标记点与在xoz平面上的投影点,与为G0中第j列荧光标记点与在xoz平面上的投影点。
本发明的工作原理为:与现有的风洞试验模型弹性变形光学测量方法不同,本发明通过在模型表面喷涂10微米厚的圆形荧光材料标记点,利用紫外光激发荧光材料,使标记点发出绿光,再经装有(绿光带通)滤光片的相机成像,消除传统方法因试验模型镜面反射所致的眩光与镜面成像干扰;同时,利用风洞迎角机构变换荧光标记点与相机的视角关系,获得测量空间内的多视角荧光标记点图像,在风洞试验现场完成相机内外参数的标定,既确保试验模型在其整个姿态变化过程中,标记点不受眩光与镜面成像干扰,又确保相机标定光照环境与试验时光照环境的一致性,以提升风洞试验模型弹性变形视频测量的精准度和可靠性。
Claims (5)
1.一种抗眩光的风洞试验模型弹性变形视频测量方法,其特征在于:包括如下步骤:
步骤一、采用绿色荧光喷漆在试验模型上喷涂形成圆形荧光标记点;
步骤二、在地面实验室,采用紫光照明,并利用在镜头前安装有绿光带通滤光片的相机采集圆形荧光标记点的图像,得到各个圆形荧光标记点在风洞坐标系oxyz下的三维坐标;
步骤三、将试验模型安装在风洞中,再将风洞试验模型弹性变形视频测量系统的光源换成紫光光源,并在风洞试验模型弹性变形视频测量系统的相机镜头前安装绿光带通滤光片;
步骤四、对风洞试验模型弹性变形视频测量系统的相机内外参数进行标定;
步骤五、风洞试验时,依次将试验模型运动至指定迎角αi并利用步骤三的相机采集αi对应荧光标记点图像G’αi,得到{G’α1,G’α2,…,G’αi,…,G’αn};
步骤六、利用步骤四得到的相机内外参数、给定αi时对应的荧光标记点图像G’αi,计算G’αi中标记点的三维坐标;
步骤七、利用步骤六计算得到的荧光标记点的三维坐标,计算各截面的弯曲与扭转变形。
2.根据权利要求1所述的抗眩光的风洞试验模型弹性变形视频测量方法,其特征在于:所述喷涂厚度为10微米,所有荧光标记点在风洞试验模型弹性变形视频测量系统的相机中成像图像至少包含30个像素;每列标记点位于的截面与X轴平行。
3.根据权利要求2所述的抗眩光的风洞试验模型弹性变形视频测量方法,其特征在于:步骤二所述紫光照明采用365纳米波长的紫光灯进行照明。
4.根据权利要求3所述的抗眩光的风洞试验模型弹性变形视频测量方法,其特征在于:步骤四所述对风洞试验模型弹性变形视频测量系统的相机内外参数进行标定的方法为:
设位于中性轴对应的翼表面上的标记点为位于机翼的前缘的标记点为位于机翼的后缘的标记点为设试验的姿态运动范围按升序排序为{α1,α2,…,αn},
(1)将试验的姿态运动范围扩展为{α0,α1,α2,…,αn,αn+1},其中α0=α1-2度,αn+1=αn+3度;
(2)试验前,先启动迎角机构,按照{α0,α1,α2,…,αi,…,αn,αn+1}的顺序,依次将试验模型运动至指定迎角αi,并利用步骤三的相机采集αi对应荧光标记点图像Gαi,得到{Gα0,Gα1,Gα2,…,Gαi,…,Gαn,Gαn+1};
(3)利用步骤二得到的各个圆形荧光标记点在风洞坐标系oxyz下的三维坐标和{Gα0,Gα1,Gα2,…,Gαi,…,Gαn,Gαn+1},标定风洞试验模型弹性变形视频测量系统的相机内参数;
(4)通过试验模型迎角为零时的荧光标记点图像G0,计算风洞试验模型弹性变形视频测量系统的相机的外参数。
5.根据权利要求4所述的抗眩光的风洞试验模型弹性变形视频测量方法,其特征在于:步骤七、利用荧光标记点的三维坐标,计算第j个截面标记点弯曲与扭转变形的方法如下:
将与机身连接的翼根部的第1个截面作为基准,设G’αi中第j个截面中性轴上的标记点为第1个截面中性轴上的标记点为则第j个截面弯曲变形
式中表示点的Z坐标,表示点的Z坐标,表示点的Z坐标,表示点的Z坐标;
第j个截面扭转变形为
式中
其中,与为G’αi中第j列荧光标记点与在xoz平面上的投影点,与为G0中第j列荧光标记点与在xoz平面上的投影点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810315108.8A CN108507754B (zh) | 2018-04-10 | 2018-04-10 | 抗眩光的风洞试验模型弹性变形视频测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810315108.8A CN108507754B (zh) | 2018-04-10 | 2018-04-10 | 抗眩光的风洞试验模型弹性变形视频测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108507754A CN108507754A (zh) | 2018-09-07 |
CN108507754B true CN108507754B (zh) | 2019-12-10 |
Family
ID=63381272
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810315108.8A Active CN108507754B (zh) | 2018-04-10 | 2018-04-10 | 抗眩光的风洞试验模型弹性变形视频测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108507754B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109100112B (zh) * | 2018-09-29 | 2020-01-17 | 浙江大学 | 一种基于计算机视觉的气弹模型风洞试验平面风振响应测试方法 |
CN111272380B (zh) * | 2020-02-24 | 2022-02-11 | 中国空气动力研究与发展中心高速空气动力研究所 | 风洞试验模型位姿视频测量的风轴系自标定方法 |
CN111413064A (zh) * | 2020-03-27 | 2020-07-14 | 智方达(天津)科技有限公司 | 一种对风洞中飞行器模型的响应测量方法 |
CN112082749A (zh) * | 2020-09-23 | 2020-12-15 | 吉林德蕴电气集团股份有限公司 | 一种隔离开关监测装置 |
CN112333890B (zh) * | 2020-11-10 | 2022-08-19 | 四川九洲光电科技股份有限公司 | 一种防眩光调节控制方法、装置、存储介质及电子设备 |
CN113405762A (zh) * | 2021-06-04 | 2021-09-17 | 中国航空工业集团公司沈阳空气动力研究所 | 一种风洞模型形变姿态视觉测量方法 |
CN113237628B (zh) * | 2021-07-08 | 2021-09-21 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种低速风洞水平自由飞模型姿态测量方法 |
CN114926071B (zh) * | 2022-06-02 | 2023-05-09 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种风洞群真空系统的动态调度方法及存储介质 |
CN116593121B (zh) * | 2023-07-12 | 2023-10-24 | 中国航空工业集团公司沈阳空气动力研究所 | 一种基于监控相机的飞行器模型振动测量方法 |
CN117232395B (zh) * | 2023-11-14 | 2024-01-23 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种压敏漆图像激波位置自动识别方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104776861B (zh) * | 2015-04-15 | 2017-05-24 | 大连理工大学 | 用于风洞视觉测量系统标定的柔性标定装置 |
CN105444982B (zh) * | 2015-11-24 | 2017-12-19 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种外挂物分离轨迹风洞试验的单目视频测量方法 |
CN107452054B (zh) * | 2017-08-08 | 2020-05-22 | 西华大学 | 一种外挂物风洞模型环绕式编码标记点的编解码方法 |
-
2018
- 2018-04-10 CN CN201810315108.8A patent/CN108507754B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108507754A (zh) | 2018-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108507754B (zh) | 抗眩光的风洞试验模型弹性变形视频测量方法 | |
CN108759714B (zh) | 一种多线激光轮廓传感器坐标系融合及转轴标定方法 | |
CN103453849B (zh) | 多光学传感器协同的复杂曲面零件三维测量方法与系统 | |
ES2399883T3 (es) | Procedimiento y sistema para la visualización de errores de superficie | |
CN101915563A (zh) | 一种飞机舵偏转角测量方法 | |
CN113155399B (zh) | 高速飞行器表面压力与变形三维连续分布同步测量方法 | |
Schneider | Analysis of SPR measurements from HART II | |
CN104567679A (zh) | 一种涡轮叶片视觉检测的系统 | |
CN104748693A (zh) | 一种基于双目立体视觉的桨叶剖面扭转刚度测量系统 | |
CN105300339A (zh) | 用于将区段对准的方法和设备 | |
CN109072881A (zh) | 用于测量风能设备的转子叶片的表面的测量系统 | |
CN105447910A (zh) | 一种航空发动机压气机叶片叶尖缺损部位三维重建方法 | |
CN102735204A (zh) | 一种基于弦线的航空薄壁叶片加工扭曲度误差测量方法 | |
CN113375594B (zh) | 飞机蒙皮型面数字化检测方法 | |
CN113405762A (zh) | 一种风洞模型形变姿态视觉测量方法 | |
CN207675408U (zh) | 应用于高超声速低密度风洞的光纤气动力测量天平 | |
CN112197713A (zh) | 直升机旋翼桨叶变形及表面压力同步测量装置及方法 | |
CN103292760A (zh) | 一种薄壁叶片误差分析方法 | |
Abrego et al. | Blade displacement measurement technique applied to a full-scale rotor test | |
CN106323587B (zh) | 机翼风洞试验模型弹性变形的单目视频高精度测量方法 | |
CN205373634U (zh) | 风力发电机组叶片表面坐标测量定位系统 | |
Davidson et al. | Time-resolved surface pressure and model deformation measurements in an industrial transonic wind tunnel | |
CN112665504B (zh) | 一种基于外型面扫描的水平测量方法 | |
Abrego et al. | Summary of full-scale blade displacement measurements of the UH-60A airloads rotor | |
CN111380476B (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 |