CN105784220A - 一种提高压力敏感涂料测量精准度的方法 - Google Patents
一种提高压力敏感涂料测量精准度的方法 Download PDFInfo
- Publication number
- CN105784220A CN105784220A CN201610133744.XA CN201610133744A CN105784220A CN 105784220 A CN105784220 A CN 105784220A CN 201610133744 A CN201610133744 A CN 201610133744A CN 105784220 A CN105784220 A CN 105784220A
- Authority
- CN
- China
- Prior art keywords
- image
- wind
- pressure
- temperature
- blowing
- 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.)
- Granted
Links
- 239000011248 coating agent Substances 0.000 title claims abstract description 153
- 238000000576 coating method Methods 0.000 title claims abstract description 153
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000005259 measurement Methods 0.000 title abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 37
- 230000002277 temperature effect Effects 0.000 claims abstract description 12
- 238000007664 blowing Methods 0.000 claims description 89
- 230000005284 excitation Effects 0.000 claims description 27
- 238000001914 filtration Methods 0.000 claims description 18
- 238000003384 imaging method Methods 0.000 claims description 18
- 230000004044 response Effects 0.000 claims description 17
- 230000003068 static effect Effects 0.000 claims description 17
- 238000011088 calibration curve Methods 0.000 claims description 16
- 238000002372 labelling Methods 0.000 claims description 12
- 230000033228 biological regulation Effects 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000001105 regulatory effect Effects 0.000 claims description 9
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 8
- 238000005424 photoluminescence Methods 0.000 claims description 8
- 238000005507 spraying Methods 0.000 claims description 8
- YTAHJIFKAKIKAV-XNMGPUDCSA-N [(1R)-3-morpholin-4-yl-1-phenylpropyl] N-[(3S)-2-oxo-5-phenyl-1,3-dihydro-1,4-benzodiazepin-3-yl]carbamate Chemical compound O=C1[C@H](N=C(C2=C(N1)C=CC=C2)C1=CC=CC=C1)NC(O[C@H](CCN1CCOCC1)C1=CC=CC=C1)=O YTAHJIFKAKIKAV-XNMGPUDCSA-N 0.000 claims description 7
- 238000009434 installation Methods 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 6
- 238000003705 background correction Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000003825 pressing Methods 0.000 claims description 5
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 4
- 238000001035 drying Methods 0.000 claims description 4
- 238000005316 response function Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 6
- 239000003973 paint Substances 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- CSCPPACGZOOCGX-UHFFFAOYSA-N Acetone Chemical compound CC(C)=O CSCPPACGZOOCGX-UHFFFAOYSA-N 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 235000007926 Craterellus fallax Nutrition 0.000 description 1
- 240000007175 Datura inoxia Species 0.000 description 1
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000084 colloidal system Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 230000037452 priming Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L11/00—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00
- G01L11/02—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00 by optical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L19/00—Details of, or accessories for, apparatus for measuring steady or quasi-steady pressure of a fluent medium insofar as such details or accessories are not special to particular types of pressure gauges
- G01L19/04—Means for compensating for effects of changes of temperature, i.e. other than electric compensation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20132—Image cropping
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明提供了一种提高压力敏感涂料测量精准度的方法,该方法利用压力敏感涂料的温度敏感特性,基于同一次压力敏感涂料风洞试验采集到的图像数据,同时获得模型表面的温度分布与温度修正后的表面压力分布。本发明的提高压力敏感涂料测量精准度的方法可以获得更为丰富的结果数据,具有成本低,精度高的优点,同时考虑了压力敏感涂料的温度效应影响,提高了压力敏感涂料测量结果的精准度。
Description
技术领域
本发明属于航空航天工业空气动力学风洞试验技术领域,具体涉及一种提高压力敏感涂料测量精准度的方法。
背景技术
压力敏感涂料受到一定波长光源照射时,涂层中的探针分子会受到激发,进入电子激发态,并通过发射波长更长的光来耗散所吸收的能量,同时,激发态的探针分子与空气中渗透进入涂层并在其中扩散的氧分子相碰撞,转移所吸收的能量,回到基态。上述过程被称为Stern-Volmer过程,描述这一过程的关系式为Stern-Volmer关系式。Stern-Volmer关系式建立了表面压力与涂料发光强度之间的函数关系。工程应用中,Stern-Volmer关系式常用形式为:
(1)
其中,P和I分别为涂层表面压力及其光致发光强度,Ai(T)为Stern-Volmer系数,与温度有关,一般通过地面压力敏感涂料预先标定实验获得。“ref”表示参考状态,n为Stern-Volmer关系式阶数,一般情况下,n≤3。工程应用中,一般选择无风状态wind-off状态为参考状态,与风洞试验的有风状态wind-on时的发光强度做比来计算模型表面压力。
由于压力敏感涂料的发光强度本质上是温度相关的,很多压力敏感涂料温度灵敏度高于1%/K,对于跨超声速暂冲式风洞,吹风前后模型表面温度之差可达10oC以上,若不进行温度效应修正,压力敏感涂料风洞试验期间以及吹风前后的模型表面温度变化会引起较大的压力敏感涂料测量误差,必须予以修正。
目前,修正压力敏感涂料温度效应的方法主要有以下三种:(1)原位标定;利用模型表面测压孔的压力数据来确定Stern-Volmer系数,该方法可以有效地降低因温度效应和其他误差源导致的误差,经原位标定的测量数据和测压孔测量结果可以很好的吻合,但是,这种方法需要在模型不同区域埋设一定数量的测压孔,降低了压力敏感涂料试验的经济性,延长了试验准备周期,并且温度效应修正精度受测压孔开设的位置与数量影响较大。(2)通过温度敏感涂(Temperature
SensitivePaint,TSP)或红外相机获取模型表面温度分布;这种方法温度效应修正精度较高,可以获得更高精准度的压力敏感涂料试验结果,但是,极大地增加了试验成本与试验周期,修正方法较为复杂。(3)双组份压力敏感涂料;将具有相同温度敏感特性可发出不同波长荧光的光敏分子掺入涂料胶体之中,从而形成了双组份的压敏涂料,然而,双组份压敏涂料虽然可以抑制涂料的温度敏感性,解决光照不均匀或被测物体形变引起局部位置光照射强度发生变化的问题,但是各组份之间在同一涂层内容易产生相互作用,从而引起较为严重的光降解和较低的发光强度,同时,各个组份往往会存在一定程度的光谱重叠,也会降低涂料的压力灵敏度。
发明内容
本发明要解决的技术问题是提供一种提高压力敏感涂料测量精准度的方法。
本发明的提高压力敏感涂料测量精准度的方法,其特点是,包括以下步骤:
a.模型及校准样片表面喷涂和固化压力敏感涂料;
b.风洞中安装试验模型,调试压力敏感涂料测量系统;
c.采集背景图像,风洞试验前,关闭风洞光源和压力敏感涂料激发光源,按照风洞试验模型状态要求,每个模型状态采集20幅以上的背景图像;
d.风洞试验前,记录模型表面温度、当地大气压值;
e.采集吹风前的wind-off图像,打开压力敏感涂料激发光源系统,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风前的wind-off图像;
f.采集吹风时的wind-on图像,打开压力敏感涂料激发光源系统,开始风洞试验,在风洞流场和模型表面温度稳定后,采用定攻角、定马赫数的方式,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风时的wind-on图像;
g.采集吹风后的wind-off图像,风洞关车,排放风洞内余气,在风洞内气流稳定后,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风后的wind-off图像;
h.利用步骤a的校准样片进行压力敏感涂料静态校准,获得压力敏感涂料静态校准曲线;
i.利用步骤a的校准样片进行压力敏感涂料温度敏感响应特性校准,获得压力敏感涂料温度敏感响应校准曲线;
j.计算模型表面温度分布,根据步骤d记录的吹风前模型表面温度、背景图像、吹风前的wind-off图像、吹风后的wind-off图像和步骤i获得的压力敏感涂料温度响应特性校准曲线计算吹风后的模型表面温度分布;
k.计算模型表面压力分布,根据步骤d记录的当地大气压值、背景图像、吹风时的wind-on图像、吹风后的wind-off图像、步骤j获得的吹风后的模型表面温度分布和步骤h获得的压力敏感涂料静态校准曲线计算压力敏感涂料温度效应修正后的模型表面压力分布。
步骤a所述的压力敏感涂料厚度为40μm-60μm。
步骤b所述的压力敏感涂料测量系统包括激发光源系统、数字成像设备和测控系统;所述的激发光源系统提供波长为A的入射光,照射在压力敏感涂料表面,涂料发射波长为B的发射光,波长A<波长B,数字成像设备采集发射光的信号,信号传输至计算机进行图像处理;所述的测控系统用于完成激发光源系统、数字成像设备和风洞测控系统的通信与时序控制。
步骤h所述的压力敏感涂料静态校准系统包括校准箱、压力调节子系统、温度调节子系统、压力和温度控制子系统、激发光源系统以及数字成像设备。所述的校准箱提供足够的空间放置压力敏感涂料校准样片,压力调节子系统实现校准箱内压力的调节,温度调节子系统实现校准箱内校准样片表面的温度调节,压力和温度控制子系统实现压力、温度数据信息的反馈与控制,激发光源系统提供波长为A的入射光,照射在压力敏感涂料校准样片表面,涂料发射波长为B的发射光,波长A<波长B,数字成像设备采集波长为B的发射光的信号,传输至计算机进行图像处理。校准方法根据风洞试验条件确定压力变化范围和温度变化范围,采取等间隔取值方式将压力变化范围和温度变化范围离散化,分别采集温度值点与压力值点处校准样片的发光图像,以风洞试验时的当地大气压作为校准基准压力,按照Stern-Volmer关系式,采用最小二乘法拟合获得Stern-Volmer关系式系数Ai,完成压力敏感涂料静态校准,式中n为Stern-Volmer多项式的阶数,P为涂层表面压力,I为涂层表面光致发光强度,ref表示参考状态。
步骤i所述的压力敏感涂料温度敏感响应特性校准系统包括校准箱、压力调节子系统、温度调节子系统、压力和温度控制子系统、激发光源系统以及数字成像设备。所述的校准箱提供足够的空间放置压力敏感涂料校准样片,压力调节子系统实现校准箱内压力的调节,温度调节子系统实现校准箱内校准样片表面的温度调节,压力和温度控制子系统实现压力、温度数据信息的反馈与控制,激发光源系统提供波长为A的入射光,照射在压力敏感涂料校准样片表面,涂料发射波长为B的发射光,波长A<波长B,数字成像设备采集波长为B的发射光的信号,传输至计算机进行图像处理。校准方法以步骤d记录的模型表面温度为校准基准温度,根据风洞试验条件确定温度变化范围,在当地大气压环境下,采集不同温度值点处校准样片发光图像,用多项式函数描述压力敏感涂料温度敏感响应特性,采用最小二乘法拟合获得压力敏感涂料温度敏感响应函数,完成压力敏感涂料温度敏感响应特性校准,式中Ak为压力敏感涂料温度敏感响应函数系数,n为多项式函数的阶数,T为涂层表面压温度,I为涂层表面光致发光强度,ref表示参考状态。
步骤j包括以下步骤:
j1.背景图像减除,步骤e记录的吹风前的wind-off图像、步骤g记录的吹风后的wind-off图像扣除步骤c记录的背景图像;
j2.图像平均,将扣除背景图像后的步骤j1得到的吹风前后的wind-off图像进行平均运算,得到一幅吹风前的wind-off图像、一幅吹风后的wind-off图像;
j3.图像比运算,步骤j2获得的吹风后的wind-off图像与吹风前的wind-off图像进行比运算,获得发光强度比图像;
j4.图像滤波,利用中值滤波空间域滤波方法对步骤j3获得的发光强度比图像进行滤波,获得图像滤波后的发光强度比图像;
j5.发光强度与温度分布转换,利用步骤i获得的压力敏感涂料温度敏感响应校准曲线,即涂料发光强度比与表面温度之间的函数关系,将步骤j4获得的发光强度比图像转换为二维的模型表面温度图像。
步骤k包括以下步骤:
k1.背景图像减除,步骤f记录的吹风时的wind-on图像、步骤g记录的吹风后的wind-off图像扣除步骤c记录的背景图像;
k2.图像平均,将扣除背景图像后的步骤k1得到的吹风时的wind-on图像与吹风后的wind-off图像进行平均运算,得到一幅吹风时的wind-on图像与一幅吹风后的wind-off图像;
k3.标记点识别,确定步骤k2获得的一幅吹风时的wind-on图像与一幅吹风后的wind-off图像中标记点的序号与位置;
k4.图像配准,根据k3获得的吹风时的wind-on图像与吹风后的wind-off图像中标记点的序号与位置,将wind-on图像配准映射到wind-off图像;
k5.图像比运算,步骤k4获得的配准后的吹风后的吹风时的wind-on图像与吹风后的wind-off图像进行比运算,获得发光强度比图像;
k6.图像滤波,利用中值滤波空间域滤波方法对步骤k5获得的发光强度比图像进行滤波,获得图像滤波后的发光强度比图像;
k7.发光强度与压力分布转换,利用步骤h获得的压力敏感涂料静态校准曲线,即涂料发光强度比、表面温度与表面压力之间的函数关系,将步骤k6获得的发光强度比图像转换与步骤j得到的二维的模型表面温度图像转换为修正了涂料温度效应的二维表面压力分布;
k8.压力三维场重构,根据步骤k3获得的标记点二维坐标以及标记点在模型表面的三维坐标,获得二维图像到三维模型表面网格的映射关系,从而实现步骤k7获得的二维表面压力分布到三维模型表面网格的映射,最终获得三维空间压力场分布。
本发明的提高压力敏感涂料测量精准度的方法利用了压力敏感涂料的温度敏感特性,基于同一次压力敏感涂料风洞试验采集到的原始图像数据,同时获得了模型表面的温度分布与温度修正后的表面压力分布。本发明的提高压力敏感涂料测量精准度的方法可以获得更为丰富的结果数据,具有成本低,精度高的优点,考虑了压力敏感涂料的温度效应影响,提高了压力敏感涂料测量结果的精准度。
附图说明
图1为本发明的提高压力敏感涂料测量精准度的方法的试验流程图;
图2为本发明的提高压力敏感涂料测量精准度的方法的图像后处理流程图。
具体实施方式
下面结合附图和实施例具体说明本发明。
以下实施例仅用于说明本发明,而并非对本发明的限制。有关技术领域的人员在不脱离本发明的精神和范围的情况下,还可以做出各种变化、替换和变型,因此同等的技术方案也属于本发明的范畴。
实施例
图1为本发明的提高压力敏感涂料测量精准度的方法的试验流程图;图2为本发明的提高压力敏感涂料测量精准度的方法的图像后处理流程图。
具体工作步骤如下:
a.模型及校准样片表面喷涂和固化压力敏感涂料。
喷涂前,喷涂人员应仔细检查喷涂试件表面有无凹坑、毛刺或其他缺陷,采用800目以上的细砂纸打磨、抛光有缺陷区域,用腻子或原子灰封堵凹坑。处理完毕后,用丙酮将模型及试件表面清洗擦拭干净。之后喷涂美国ISSI公司SCR型底漆,厚度为20μm,喷涂完成后在90oC下烘烤5个小时。烘烤完成后,喷涂美国ISSI公司FIB型面漆,厚度为20μm,喷涂完成后在65oC下烘烤5个小时。
b.风洞中安装试验模型,调试压力敏感涂料测量系统。
在风洞的支撑机构上安装试验模型。将压力敏感涂料测量系统的设备安装在风洞中,设备安装完成后,进行测控系统与试验设施预运转,检查测控系统和试验设施的可靠性。同时检验光路设计、模型安装等相关项目。
压力敏感涂料测量系统包括激发光源系统、数字成像设备和测控系统;所述的激发光源系统提供波长为A的入射光,照射在压力敏感涂料表面,涂料发射波长为B的发射光,波长A<波长B,数字成像设备采集发射光的信号,信号传输至计算机进行图像处理,测控系统用于完成激发光源系统、数字成像设备和风洞测控系统的通信与时序控制。
c.采集背景图像。风洞试验前,关闭风洞环境光源以及压力敏感涂料激发光源,按照风洞试验模型状态要求,每个模型状态采集20幅以上的背景图像。
d.风洞试验前,记录模型表面温度、当地大气压值。
e.采集吹风前的wind-off图像。打开压力敏感涂料激发光源系统,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风前的wind-off图像。
f.采集吹风时的wind-on图像。打开压力敏感涂料激发光源系统,开始风洞试验,在风洞流场和模型表面温度稳定后,采用定攻角、定马赫数的方式,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风时的wind-on图像。
g.采集吹风后的wind-off图像。风洞关车,排放风洞内余气,在风洞内气流稳定后,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风后wind-off图像。
h.利用步骤a的校准样片进行压力敏感涂料静态校准,获得压力敏感涂料静态校准曲线。
压力敏感涂料静态校准系统主要由校准箱、压力调节子系统、温度调节子系统、压力和温度控制子系统、激发光源以及以及数字成像设备。安装调试压力敏感涂料静态校准系统后,根据试验条件预先确定模型表面的压力范围和温度范围,采取等间隔取值方式将压力范围和温度范围离散化,分别采集不同温度与压力值点处校准样片的发光图像,以风洞试验时当地大气压作为校准基准压力,完成涂料的温度与压力校准,按照Stern-Volmer关系式,采用最小二乘法拟合获得Stern-Volmer关系式系数,式中n为Stern-Volmer多项式的阶数,P为涂层表面压力,I为涂层表面光致发光强度,ref表示参考状态。本实施例采用二阶多项式函数拟合,即,式中P为涂层表面压力,I为涂层表面光致发光强度,ref表示参考状态,参考状态为步骤d测量的温度与大气压力状态。
i.利用步骤a的校准样片进行压力敏感涂料温度敏感响应特性校准,获得压力敏感涂料温度敏感响应校准曲线。
利用压力敏感涂料的温度敏感特性,将压力敏感涂料看作温度敏感涂料,按照常规温度敏感涂料静态校准流程进行压力敏感涂料温度敏感响应特性校准:以吹风前模型表面温度为校准基准温度,根据试验条件预先确定模型表面温度范围,在当地大气压环境下,采集不同温度下的校准样片发光强度,完成涂料校准。
本实施例中,温度敏感涂料发光强度与温度之间的函数关系由二阶多项式函数进行描述,即,式中,I为涂层表面光致发光强度,T为模型表面温度,ref表示参考状态,参考状态为步骤d测量的温度与大气压力状态。
j.计算吹风后模型表面温度分布,步骤如下:
j1.背景图像减除,步骤e记录的吹风前的wind-off图像、步骤g记录的吹风后的wind-off图像扣除步骤c记录的背景图像;
j2.图像平均,将扣除背景图像后的步骤j1得到的吹风前后的wind-off图像进行平均运算,得到一幅吹风前的wind-off图像、一幅吹风后的wind-off图像;
j3.图像比运算,步骤j2获得的吹风后的wind-off图像与吹风前的wind-off图像进行比运算,获得发光强度比图像;
j4.图像滤波,利用中值滤波空间域滤波方法对步骤j3获得的发光强度比图像进行滤波,获得图像滤波后的发光强度比图像;
j5.发光强度与温度分布转换,利用步骤i获得的压力敏感涂料温度敏感响应校准曲线,即涂料发光强度比与表面温度之间的函数关系,将步骤j4获得的发光强度比图像转换为二维的模型表面温度图像。
k.计算吹风后模型表面压力分布,步骤如下:
k1.背景图像减除,步骤f记录的吹风时的wind-on图像、步骤g记录的吹风后的wind-off图像扣除步骤c记录的背景图像;
k2.图像平均,将扣除背景图像后的步骤k1得到的吹风时的wind-on图像与吹风后的wind-off图像进行平均运算,得到一幅吹风时的wind-on图像与一幅吹风后的wind-off图像;
k3.标记点识别,确定步骤k2获得的一幅吹风时的wind-on图像与一幅吹风后的wind-off图像中标记点的序号与位置;
k4.图像配准,根据k3获得的吹风时的wind-on图像与吹风后的wind-off图像中标记点的序号与位置,将wind-on图像配准映射到wind-off图像;
k5.图像比运算,步骤k4获得的配准后的吹风后的吹风时的wind-on图像与吹风后的wind-off图像进行比运算,获得发光强度比图像;
k6.图像滤波,利用中值滤波空间域滤波方法对步骤k5获得的发光强度比图像进行滤波,获得图像滤波后的发光强度比图像;
k7.发光强度与压力分布转换,利用步骤h获得的压力敏感涂料静态校准曲线,即涂料发光强度比、表面温度与表面压力之间的函数关系,将步骤k6获得的发光强度比图像转换与步骤j得到的二维的模型表面温度图像转换为修正了涂料温度效应的二维表面压力分布;
k8.压力三维场重构,根据步骤k3获得的标记点二维坐标以及标记点在模型表面的三维坐标,获得二维图像到三维模型表面网格的映射关系,从而实现步骤k7获得的二维表面压力分布到三维模型表面网格的映射,最终获得三维空间压力场分布。
Claims (7)
1.一种提高压力敏感涂料测量精准度的方法,其特征在于,包括以下步骤:
a.模型及校准样片表面喷涂和固化压力敏感涂料;
b.风洞中安装试验模型,调试压力敏感涂料测量系统;
c.采集背景图像,风洞试验前,关闭风洞环境光源和压力敏感涂料激发光源,根据风洞试验模型状态要求,每个模型状态采集20幅以上的背景图像;
d.风洞试验前,记录模型表面温度、当地大气压值;
e.采集吹风前的wind-off图像,打开压力敏感涂料激发光源系统,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风前的wind-off图像;
f.采集吹风时的wind-on图像,打开压力敏感涂料激发光源系统,开始风洞试验,在风洞流场和模型表面温度稳定后,采用定攻角、定马赫数的方式,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风时的wind-on图像;
g.采集吹风后的wind-off图像,风洞关车,排放风洞内余气,在风洞内气流稳定后,按照步骤c所述的风洞试验模型状态要求,每个模型状态采集20幅以上的吹风后的wind-off图像;
h.利用步骤a的校准样片进行压力敏感涂料静态校准,获得压力敏感涂料静态校准曲线;
i.利用步骤a的校准样片进行压力敏感涂料温度敏感响应特性校准,获得压力敏感涂料温度敏感响应校准曲线;
j.计算吹风后模型表面温度分布,根据步骤c记录的背景图像、步骤d记录的风洞试验前模型表面温度、步骤e记录的吹风前的wind-off图像、步骤g记录的吹风后的wind-off图像和步骤i获得的压力敏感涂料温度响应特性校准曲线计算吹风后的模型表面温度分布;
k.计算吹风后模型表面压力分布,根据步骤c记录的背景图像、步骤d记录的风洞试验前当地大气压值、步骤f记录的吹风时的wind-on图像、步骤g记录的吹风后的wind-off图像、步骤h获得的压力敏感涂料静态校准曲线和步骤j获得的吹风后的模型表面温度分布计算进行了压力敏感涂料温度效应修正后的吹风后模型表面压力分布。
2.根据权利要求1所述的提高压力敏感涂料测量精准度的方法,其特征在于,步骤a所述的压力敏感涂料厚度为40μm-60μm。
3.根据权利要求1所述的提高压力敏感涂料测量精准度的方法,其特征在于,步骤b所述的压力敏感涂料测量系统包括激发光源系统、数字成像设备和测控系统;所述的激发光源系统提供波长为A的入射光,照射在压力敏感涂料表面,涂料发射波长为B的发射光,波长A<波长B,数字成像设备采集波长为B的发射光的信号,传输至计算机进行图像处理;所述的测控系统用于完成激发光源系统、数字成像设备和风洞测控系统的通信与时序控制。
4.根据权利要求1所述的提高压力敏感涂料测量精准度的方法,其特征在于,所述的步骤h的压力敏感涂料静态校准系统包括校准箱、压力调节子系统、温度调节子系统、压力和温度控制子系统、激发光源系统以及数字成像设备;所述的校准箱中放置压力敏感涂料校准样片,压力调节子系统实现校准箱内压力的调节,温度调节子系统实现校准箱内校准样片表面的温度调节,压力和温度控制子系统实现压力、温度数据信息的反馈与控制,激发光源系统提供波长为A的入射光,照射在压力敏感涂料校准样片表面,涂料发射波长为B的发射光,波长A<波长B,数字成像设备采集波长为B的发射光的信号,传输至计算机进行图像处理;
校准时根据风洞试验条件确定压力变化范围和温度变化范围,采取等间隔取值方式将压力变化范围和温度变化范围离散化,分别采集温度值点与压力值点处校准样片的发光图像,以风洞试验时的当地大气压作为校准基准压力,按照Stern-Volmer关系式,采用最小二乘法拟合获得Stern-Volmer关系式系数Ai,式中n为Stern-Volmer多项式的阶数,P为涂层表面压力,I为涂层表面光致发光强度,ref表示参考状态。
5.根据权利要求1所述的提高压力敏感涂料测量精准度的方法,其特征在于,步骤i所述的压力敏感涂料温度敏感响应特性校准系统包括校准箱、压力调节子系统、温度调节子系统、压力和温度控制子系统、激发光源系统以及数字成像设备;所述的校准箱放置压力敏感涂料校准样片,压力调节子系统实现校准箱内压力的调节,温度调节子系统实现校准箱内校准样片表面的温度调节,压力和温度控制子系统实现压力、温度数据信息的反馈与控制,激发光源系统提供波长为A的入射光,照射在压力敏感涂料校准样片表面,涂料发射波长为B的发射光,波长A<波长B,数字成像设备采集波长为B的发射光的信号,传输至计算机进行图像处理;
校准时以步骤d记录的模型表面温度为校准基准温度,根据风洞试验条件确定温度变化范围,在当地大气压环境下,采集温度值点处校准样片发光图像,用函数描述压力敏感涂料温度敏感响应特性,采用最小二乘法拟合获得压力敏感涂料温度敏感响应函数,式中Ak为压力敏感涂料温度敏感响应函数系数,n为多项式函数的阶数,T为涂层表面压温度,I为涂层表面光致发光强度,ref表示参考状态。
6.根据权利要求1所述的提高压力敏感涂料测量精准度的方法,其特征在于,步骤j包括以下步骤:
j1.背景图像减除,步骤e记录的吹风前的wind-off图像、步骤g记录的吹风后的wind-off图像扣除步骤c记录的背景图像;
j2.图像平均,将扣除背景图像后的步骤j1得到的吹风前后的wind-off图像进行平均运算,得到一幅吹风前的wind-off图像、一幅吹风后的wind-off图像;
j3.图像比运算,步骤j2获得的吹风后的wind-off图像与吹风前的wind-off图像进行比运算,获得发光强度比图像;
j4.图像滤波,利用中值滤波空间域滤波方法对步骤j3获得的发光强度比图像进行滤波,获得图像滤波后的发光强度比图像;
j5.发光强度与温度分布转换,利用步骤i获得的压力敏感涂料温度敏感响应校准曲线,即涂料发光强度比与表面温度之间的函数关系,将步骤j4获得的发光强度比图像转换为二维的模型表面温度图像。
7.根据权利要求1所述的提高压力敏感涂料测量精准度的方法,其特征在于,步骤k包括以下步骤:
k1.背景图像减除,步骤f记录的吹风时的wind-on图像、步骤g记录的吹风后的wind-off图像扣除步骤c记录的背景图像;
k2.图像平均,将扣除背景图像后的步骤k1得到的吹风时的wind-on图像与吹风后的wind-off图像进行平均运算,得到一幅吹风时的wind-on图像与一幅吹风后的wind-off图像;
k3.标记点识别,确定步骤k2获得的一幅吹风时的wind-on图像与一幅吹风后的wind-off图像中标记点的序号与位置;
k4.图像配准,根据k3获得的吹风时的wind-on图像与吹风后的wind-off图像中标记点的序号与位置,将wind-on图像配准映射到wind-off图像;
k5.图像比运算,步骤k4获得的配准后的吹风后的吹风时的wind-on图像与吹风后的wind-off图像进行比运算,获得发光强度比图像;
k6.图像滤波,利用中值滤波空间域滤波方法对步骤k5获得的发光强度比图像进行滤波,获得图像滤波后的发光强度比图像;
k7.发光强度与压力分布转换,利用步骤h获得的压力敏感涂料静态校准曲线,即涂料发光强度比、表面温度与表面压力之间的函数关系,将步骤k6获得的发光强度比图像转换与步骤j得到的二维的模型表面温度图像转换为修正了涂料温度效应的二维表面压力分布;
k8.压力三维场重构,根据步骤k3获得的标记点二维坐标以及标记点在模型表面的三维坐标,获得二维图像到三维模型表面网格的映射关系,从而实现步骤k7获得的二维表面压力分布到三维模型表面网格的映射,最终获得三维空间压力场分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610133744.XA CN105784220B (zh) | 2016-05-03 | 2016-05-03 | 一种提高压力敏感涂料测量精准度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610133744.XA CN105784220B (zh) | 2016-05-03 | 2016-05-03 | 一种提高压力敏感涂料测量精准度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105784220A true CN105784220A (zh) | 2016-07-20 |
CN105784220B CN105784220B (zh) | 2018-10-09 |
Family
ID=56386607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610133744.XA Expired - Fee Related CN105784220B (zh) | 2016-05-03 | 2016-05-03 | 一种提高压力敏感涂料测量精准度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105784220B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501227A (zh) * | 2016-10-25 | 2017-03-15 | 中国航空工业集团公司沈阳空气动力研究所 | 基于压力敏感涂料探针分子荧光寿命的测量方法 |
CN107702878A (zh) * | 2017-08-17 | 2018-02-16 | 上海交通大学 | 一种基于aao模板的柔性快响应psp测试装置、方法与应用 |
CN108318215A (zh) * | 2018-02-08 | 2018-07-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种压敏漆试验图像采集幅数的确定方法 |
CN109580080A (zh) * | 2018-12-29 | 2019-04-05 | 天津大学 | 一种基于飞秒激光成丝测量流场压力的装置及方法 |
CN111458101A (zh) * | 2020-03-12 | 2020-07-28 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞固定翼模型表面脉动压力测量方法 |
CN114323543A (zh) * | 2022-03-10 | 2022-04-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种提高压敏漆试验效率的方法 |
CN114323548A (zh) * | 2022-03-14 | 2022-04-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种适用于不同参考状态下的压敏漆涂料校准方法 |
CN114441090A (zh) * | 2022-04-11 | 2022-05-06 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种快速响应压敏漆温度效应修正方法 |
CN114459728A (zh) * | 2022-04-13 | 2022-05-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种低温温敏漆转捩测量试验方法 |
CN114739626A (zh) * | 2022-06-13 | 2022-07-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于快速响应压敏漆的旋转叶栅测压试验方法 |
CN115952826A (zh) * | 2023-03-09 | 2023-04-11 | 中国空气动力研究与发展中心低速空气动力研究所 | 基于人工神经网络的压力敏感涂料性能预测和测压方法 |
US11808058B2 (en) | 2013-07-12 | 2023-11-07 | Invue Security Products Inc. | Merchandise security devices for use with an electronic key |
CN117249967A (zh) * | 2023-11-20 | 2023-12-19 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于双组份压敏漆技术的压力与温度同步测量方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030111615A1 (en) * | 2001-12-17 | 2003-06-19 | Benne Michael E. | Method and apparatus to correct for the temperature sensitivity of pressure sensitive paint |
JP2005029767A (ja) * | 2003-07-11 | 2005-02-03 | National Aerospace Laboratory Of Japan | 感圧・感温複合機能塗料 |
CN101089575A (zh) * | 2007-05-11 | 2007-12-19 | 北京航空航天大学 | 用于压力敏感涂料测压系统的校准装置 |
CN101290259A (zh) * | 2008-06-13 | 2008-10-22 | 西北工业大学 | 光学压力敏感涂料测量压力的精度提高方法 |
CN102081001A (zh) * | 2009-11-26 | 2011-06-01 | 上海莫克电子技术有限公司 | 一种高精度智能压力变送器 |
CN102928192A (zh) * | 2012-10-10 | 2013-02-13 | 中国航空工业集团公司沈阳空气动力研究所 | 一种用于psp测量的压力分布修正方法 |
CN202939130U (zh) * | 2012-12-11 | 2013-05-15 | 中国航空工业集团公司沈阳空气动力研究所 | 激波法压敏涂料/温敏涂料压力/温度响应时间测量装置 |
CN103267622A (zh) * | 2013-05-13 | 2013-08-28 | 中国航空工业集团公司沈阳空气动力研究所 | 一种荧光油流摩擦力场测量装置及测量方法 |
-
2016
- 2016-05-03 CN CN201610133744.XA patent/CN105784220B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030111615A1 (en) * | 2001-12-17 | 2003-06-19 | Benne Michael E. | Method and apparatus to correct for the temperature sensitivity of pressure sensitive paint |
JP2005029767A (ja) * | 2003-07-11 | 2005-02-03 | National Aerospace Laboratory Of Japan | 感圧・感温複合機能塗料 |
CN101089575A (zh) * | 2007-05-11 | 2007-12-19 | 北京航空航天大学 | 用于压力敏感涂料测压系统的校准装置 |
CN101290259A (zh) * | 2008-06-13 | 2008-10-22 | 西北工业大学 | 光学压力敏感涂料测量压力的精度提高方法 |
CN102081001A (zh) * | 2009-11-26 | 2011-06-01 | 上海莫克电子技术有限公司 | 一种高精度智能压力变送器 |
CN102928192A (zh) * | 2012-10-10 | 2013-02-13 | 中国航空工业集团公司沈阳空气动力研究所 | 一种用于psp测量的压力分布修正方法 |
CN202939130U (zh) * | 2012-12-11 | 2013-05-15 | 中国航空工业集团公司沈阳空气动力研究所 | 激波法压敏涂料/温敏涂料压力/温度响应时间测量装置 |
CN103267622A (zh) * | 2013-05-13 | 2013-08-28 | 中国航空工业集团公司沈阳空气动力研究所 | 一种荧光油流摩擦力场测量装置及测量方法 |
Non-Patent Citations (1)
Title |
---|
李国帅: "压力敏感涂料特性及其校准技术研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11808058B2 (en) | 2013-07-12 | 2023-11-07 | Invue Security Products Inc. | Merchandise security devices for use with an electronic key |
CN106501227B (zh) * | 2016-10-25 | 2019-03-01 | 中国航空工业集团公司沈阳空气动力研究所 | 基于压力敏感涂料探针分子荧光寿命的测量方法 |
CN106501227A (zh) * | 2016-10-25 | 2017-03-15 | 中国航空工业集团公司沈阳空气动力研究所 | 基于压力敏感涂料探针分子荧光寿命的测量方法 |
CN107702878A (zh) * | 2017-08-17 | 2018-02-16 | 上海交通大学 | 一种基于aao模板的柔性快响应psp测试装置、方法与应用 |
CN108318215A (zh) * | 2018-02-08 | 2018-07-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种压敏漆试验图像采集幅数的确定方法 |
CN108318215B (zh) * | 2018-02-08 | 2019-09-13 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种压敏漆试验图像采集幅数的确定方法 |
CN109580080A (zh) * | 2018-12-29 | 2019-04-05 | 天津大学 | 一种基于飞秒激光成丝测量流场压力的装置及方法 |
CN111458101A (zh) * | 2020-03-12 | 2020-07-28 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞固定翼模型表面脉动压力测量方法 |
CN114323543A (zh) * | 2022-03-10 | 2022-04-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种提高压敏漆试验效率的方法 |
CN114323548A (zh) * | 2022-03-14 | 2022-04-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种适用于不同参考状态下的压敏漆涂料校准方法 |
CN114441090A (zh) * | 2022-04-11 | 2022-05-06 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种快速响应压敏漆温度效应修正方法 |
CN114459728A (zh) * | 2022-04-13 | 2022-05-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种低温温敏漆转捩测量试验方法 |
CN114459728B (zh) * | 2022-04-13 | 2022-06-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种低温温敏漆转捩测量试验方法 |
CN114739626A (zh) * | 2022-06-13 | 2022-07-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于快速响应压敏漆的旋转叶栅测压试验方法 |
CN115952826A (zh) * | 2023-03-09 | 2023-04-11 | 中国空气动力研究与发展中心低速空气动力研究所 | 基于人工神经网络的压力敏感涂料性能预测和测压方法 |
CN117249967A (zh) * | 2023-11-20 | 2023-12-19 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于双组份压敏漆技术的压力与温度同步测量方法 |
CN117249967B (zh) * | 2023-11-20 | 2024-02-06 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于双组份压敏漆技术的压力与温度同步测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105784220B (zh) | 2018-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105784220A (zh) | 一种提高压力敏感涂料测量精准度的方法 | |
CN112945501A (zh) | 层流机翼转捩位置测量试验方法 | |
CN108318215B (zh) | 一种压敏漆试验图像采集幅数的确定方法 | |
Fey et al. | Transition detection by temperature sensitive paint at cryogenic temperatures in the European Transonic Wind tunnel (ETW) | |
CN101290259B (zh) | 光学压力敏感涂料测量压力的精度提高方法 | |
Engler et al. | Application of PSP in low speed flows | |
CN110411699A (zh) | 用于激波风洞气动热试验的遮挡区域温敏热图试验装置 | |
Roozeboom et al. | Development of unsteady pressure-sensitive paint application on NASA space launch system | |
CN112945502B (zh) | 层流机翼转捩位置测量试验系统 | |
CN113155399A (zh) | 高速飞行器表面压力与变形三维连续分布同步测量方法 | |
Dong et al. | Simultaneous pressure and deformation field measurement on helicopter rotor blades using a grid-pattern pressure-sensitive paint system | |
CN113267132B (zh) | 一种基于双色成像的全局油膜厚度在线测量方法 | |
Yorita et al. | Time-resolved pressure-sensitive paint measurements for cryogenic wind tunnel tests | |
Wu et al. | Three‐dimensional deformation measurement of aero‐engine high‐speed rotating blade surface based on stroboscopic structure digital image correlation method | |
Bell | Pressure-sensitive paint measurements on the NASA common research model in the NASA 11-ft transonic wind tunnel | |
Davidson et al. | Time-resolved surface pressure and model deformation measurements in an industrial transonic wind tunnel | |
McLachlan et al. | Pressure-sensitive paint measurements on a supersonic high-sweep oblique wing model | |
Bencic | Rotating pressure and temperature measurements on scale-model fans using luminescent paints | |
Morris et al. | Application of pressure-and temperature-sensitive paints to high-speed flows | |
Watkins et al. | Using pressure-and temperature-sensitive paint on the aftbody of a capsule entry vehicle | |
Vardaki et al. | Pressure sensitive paint measurements at the ARA transonic wind tunnel | |
Strozik et al. | Measuring System for Pressure Sensitive Paint (PSP) Calibration in the Range of Low Pressure Gains | |
Dhanagopal et al. | High-Speed Pressure Sensitive Paint Measurements of the Initial Concept 3. X Vehicle at Mach 7 | |
Watkins et al. | The development and implementation of a cryogenic pressure sensitive paint system in the national transonic facility | |
Bencic | Application of pressure sensitive paint to ice accreted wind tunnel models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181009 Termination date: 20190503 |