CN114646663A - High-efficiency measurement system and method for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hood - Google Patents
High-efficiency measurement system and method for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hood Download PDFInfo
- Publication number
- CN114646663A CN114646663A CN202210324793.7A CN202210324793A CN114646663A CN 114646663 A CN114646663 A CN 114646663A CN 202210324793 A CN202210324793 A CN 202210324793A CN 114646663 A CN114646663 A CN 114646663A
- Authority
- CN
- China
- Prior art keywords
- infrared
- radiation
- temperature
- head cover
- hood
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Radiation Pyrometers (AREA)
Abstract
高温红外头罩不同厚度材料热辐射特性高效测量系统及方法,属于热辐射测量技术领域。为了解决现有的采用数值模拟方法存在因为受到待测材料厚度以及测量误差的影响较大导致的重建结果精度交叉的问题。本发明沿辐射传递方向将红外头罩试件等分,根据辐射传输原理和能量守恒关系依次得到透过头罩材料第1层、第1~2层、……、1~n层的总辐射;随后基于能量法推演出厚度为Δ且温度为T的红外头罩材料的透过率和自身辐射与整个红外探测头罩的透过率和自身辐射之间所满足的代数关系。通过实验测量出厚度x的红外头罩的透射率和自身辐射等热辐射特性,基于代数关系计算单位厚度材料的的热辐射特性数据。主要用于获得红外头罩的热辐射特性。
The invention discloses a high-efficiency measuring system and method for thermal radiation characteristics of materials with different thicknesses of a high-temperature infrared hood, belonging to the technical field of thermal radiation measurement. In order to solve the problem that the existing numerical simulation method is greatly influenced by the thickness of the material to be measured and the measurement error, the accuracy of the reconstruction results is crossed. The invention divides the infrared head cover test piece into equal parts along the radiation transmission direction, and sequentially obtains the total radiation passing through the first layer, the first to the second layer, ..., the 1 to n layers of the head cover material according to the radiation transmission principle and the energy conservation relationship; Then based on the energy method, the algebraic relationship satisfied between the transmittance and self-radiation of the infrared head cover material with a thickness of Δ and a temperature of T and the transmittance and self-radiation of the entire infrared detection head cover was deduced. The thermal radiation characteristics such as the transmittance and self-radiation of the infrared head cover of thickness x are measured experimentally, and the thermal radiation characteristics data of the material per unit thickness are calculated based on the algebraic relationship. It is mainly used to obtain the thermal radiation characteristics of the infrared hood.
Description
技术领域technical field
本发明属于热辐射测量技术领域,具体涉及一种高温红外头罩不同厚度材料热辐射特性测量方法。The invention belongs to the technical field of thermal radiation measurement, in particular to a method for measuring thermal radiation characteristics of materials with different thicknesses of a high-temperature infrared head cover.
背景技术Background technique
飞行器在大气层内超高声速飞行时,高温红外头罩迅速成为红外探测系统气动热辐射效应的主要因素。高温使红外头罩自身热辐射增加,红外头罩表面的辐射主要集中于红外波段,高温激波气体和红外头罩产生强烈的气动热辐射效应,对探测器形成辐射干扰,对目标的红外信号形成干扰,使探测器背景亮度增加,降低目标的探测信噪比,降低了系统对目标的检测和跟踪能力,甚至导致红外探测器饱和而不能准确分辨来自目标的信号,引起目标探测、跟踪与识别能力的减弱,使红外探测系统成像质量下降。When the aircraft is flying at hypersonic speed in the atmosphere, the high temperature infrared hood quickly becomes the main factor of the aerodynamic thermal radiation effect of the infrared detection system. The high temperature increases the thermal radiation of the infrared hood itself, and the radiation on the surface of the infrared hood is mainly concentrated in the infrared band. Form interference, increase the background brightness of the detector, reduce the detection signal-to-noise ratio of the target, reduce the detection and tracking ability of the system to the target, and even cause the infrared detector to be saturated and unable to accurately distinguish the signal from the target, causing target detection, tracking and The weakening of the recognition ability reduces the imaging quality of the infrared detection system.
目前对红外头罩材料的测量研究多集中在强度、硬度﹑熔点、折射率、热导率和耐腐蚀等物理化学性能,而透过率、衰减系数等热辐射传输特性的研究较少,高温状态的数据更少,影响甚至制约了高温红外光学窗口的气动热辐射效应研究的发展,阻碍了红外探测系统在高超声速飞行器领域的应用。At present, the measurement research on infrared head cover materials mostly focuses on the physical and chemical properties such as strength, hardness, melting point, refractive index, thermal conductivity and corrosion resistance, while the research on thermal radiation transmission characteristics such as transmittance and attenuation coefficient is less. The data of the state is less, which affects or even restricts the development of the research on the aerothermal radiation effect of the high-temperature infrared optical window, and hinders the application of the infrared detection system in the field of hypersonic vehicles.
传统的红外头罩辐射特性测试方法针对同一种材料,材料厚度改变时需要开展一次乃至多次重复性实验,每次实验的温度等环境工况不能保证完全一致,红外头罩材料辐射特性测量值与真实值误差较大且造成额外的实验成本浪费。传统的思路为采用数值模拟的方法包括离散坐标法、有限体积法等对实际情况进行模拟,结合不同角度的表观辐射强度测量结果,反演重建介质内部的辐射特性场,如吸收系数场、折射率场等。重建过程受到测量误差等影响较大,且重建的精度随着介质的厚度增加而降低。因此急需一种在高温状态下实现不同厚度红外头罩材料热辐射特性一次性测量的方法,更加高效地建立高温红外头罩传输特性量化模型。The traditional infrared hood radiation characteristic test method is aimed at the same material, and it is necessary to carry out one or even multiple repeated experiments when the material thickness changes. The environmental conditions such as temperature of each experiment cannot be guaranteed to be completely consistent. The measured value of the infrared hood material radiation characteristics The error with the real value is large and the extra experiment cost is wasted. The traditional idea is to use numerical simulation methods including discrete coordinate method, finite volume method, etc. to simulate the actual situation, and combine the apparent radiation intensity measurement results from different angles to invert and reconstruct the radiation characteristic field inside the medium, such as the absorption coefficient field, Refractive index field, etc. The reconstruction process is greatly affected by measurement errors, and the reconstruction accuracy decreases with the increase of the thickness of the medium. Therefore, there is an urgent need for a method to measure the thermal radiation characteristics of infrared head cover materials with different thicknesses at one time under high temperature conditions, and to establish a more efficient quantitative model of the transmission characteristics of high temperature infrared head cover materials.
发明内容SUMMARY OF THE INVENTION
本发明为了解决现有的采用数值模拟方法存在因为受到待测材料厚度以及测量误差的影响较大导致的重建结果精度交叉的问题。The present invention solves the problem that the accuracy of the reconstruction result is crossed due to the large influence of the thickness of the material to be measured and the measurement error in the existing numerical simulation method.
一种高温红外头罩不同厚度材料热辐射特性高效测量系统,包括傅里叶变换红外光谱仪、加热炉、黑体炉、温控巡检仪和数据采集及处理系统;An efficient measurement system for thermal radiation characteristics of materials with different thicknesses of a high-temperature infrared hood, comprising a Fourier transform infrared spectrometer, a heating furnace, a blackbody furnace, a temperature control inspection instrument, and a data acquisition and processing system;
测量时,傅里叶红外光谱仪的探测镜头中心、加热炉中心以及黑体炉腔体中心设置在同一条水平线上;During measurement, the center of the detection lens of the Fourier transform infrared spectrometer, the center of the heating furnace and the center of the blackbody furnace cavity are set on the same horizontal line;
黑体炉,用于发射黑体红外辐射;在测量工作过程中,调整黑体炉改变黑体温度,以发射不同黑体温度下的红外辐射;The black body furnace is used to emit black body infrared radiation; during the measurement process, the black body furnace is adjusted to change the black body temperature to emit infrared radiation at different black body temperatures;
加热炉,用于红外头罩样品加热;在测量工作过程中,调整加热炉的温度,为红外头罩样品材料提供不同温度;Heating furnace, used for infrared hood sample heating; during the measurement process, adjust the temperature of the heating furnace to provide different temperatures for the infrared hood sample material;
温控巡检仪,用于检测和控制加热炉内的温度;Temperature control inspection instrument, used to detect and control the temperature in the heating furnace;
傅里叶红外光谱仪,用于获取透过红外头罩的黑体红外辐射;Fourier infrared spectrometer, used to obtain black body infrared radiation through the infrared hood;
数据采集及处理系统,用于采集傅里叶红外光谱仪的数据和温控巡检仪,并利用傅里叶红外光谱仪得到的信号计算得到材料在温控巡检仪所显示的温度下的法相光谱表观辐射强度。The data acquisition and processing system is used to collect the data of the Fourier transform infrared spectrometer and the temperature control inspection instrument, and use the signal obtained by the Fourier infrared spectrometer to calculate the normal phase spectrum of the material at the temperature displayed by the temperature control inspection instrument. apparent radiation intensity.
一种高温红外头罩不同厚度材料热辐射特性高效测量方法,包括以下步骤:An efficient method for measuring thermal radiation characteristics of materials with different thicknesses of a high-temperature infrared hood, comprising the following steps:
步骤一、搭建权利要求1所述的高温红外头罩不同厚度材料热辐射特性高效测量系统;
步骤二、初始阶段,不启动加热炉,加热炉内不放置样品,启动黑体炉,设定黑体炉温度为Tb,用傅里叶红外光谱仪获取黑体的红外辐射Lobj;
步骤三、将红外头罩样品放置于高温加热炉内加热,待红外头罩样品材料温度达到指定温度Twin且分布均匀后,用红外探测器获取透过红外头罩样品材料的红外辐射Ltot;Step 3: Place the infrared hood sample in a high-temperature heating furnace for heating. After the temperature of the infrared hood sample material reaches the specified temperature T win and is uniformly distributed, use an infrared detector to obtain the infrared radiation L tot passing through the infrared hood sample material. ;
步骤四、控制样品材料温度,使其温度保持Twin不变,改变黑体温度Tb状态下,重复步骤二与步骤三,获取多组黑体温度Tb,i状态下的红外辐射Lobj,i和Ltot,i;其中下标i表示第i次测量;Step 4: Control the temperature of the sample material to keep the temperature T win unchanged, and repeat
步骤五、当材料的温度不变时,其辐射特性参数为固定值,通过对试验结果的统计,利用最小二乘法拟合多组黑体温度Tb,i状态下的红外辐射Lobj,i和Ltot,i,进而得出温度Twin均匀分布的红外头罩样品材料的透过率τT,win和自身辐射LT,win;
步骤六、沿厚度方向将红外头罩样品等分成n层,根据能量守恒关系获得单位厚度Δ的红外头罩样品材料的表观光谱透过率自身辐射并基于获得表观法向光谱发射率 Step 6: Divide the infrared head cover sample into n layers equally along the thickness direction, and obtain the apparent spectral transmittance of the infrared head cover sample material with unit thickness Δ according to the energy conservation relationship self-radiation and based on Obtain the apparent normal spectral emissivity
步骤七、根据辐射传输逆问题求解算法,假设红外头罩样品材料的折射率为吸收系数为通过求解辐射传输方程计算得到该红外头罩样品材料的出射界面上任意角度的表观光谱辐射强度表观法向光谱发射率估计值以及表观光谱透过率估计值 Step 7. According to the algorithm for solving the inverse problem of radiation transmission, it is assumed that the refractive index of the sample material of the infrared head cover is The absorption coefficient is The apparent spectral radiation intensity at any angle on the exit interface of the infrared mask sample material is calculated by solving the radiation transfer equation Apparent Normal Spectral Emissivity Estimates and an estimate of apparent spectral transmittance
步骤八、将步骤六得到的红外头罩样品材料的表观法向光谱发射率与表观光谱透过率和步骤七得到的红外头罩样品材料的表观法向发射率估计值与表观光谱透过率估计值代入如下目标函数计算公式,计算得到目标函数值Fobj;Step 8. Calculate the apparent normal spectral emissivity of the infrared hood sample material obtained in Step 6 and apparent spectral transmittance and the estimated apparent normal emissivity of the infrared hood sample material obtained in step 7 and apparent spectral transmittance estimates Substitute the following objective function calculation formula, and calculate the objective function value F obj ;
步骤九、判断步骤八中的目标函数值Fobj是否小于设定阈值ξ,Step 9: Determine whether the objective function value F obj in step 8 is less than the set threshold ξ,
若是,则步骤八中所假设的红外头罩样品材料的折射率吸收系数即为该的红外头罩样品材料的真实折射率、吸收系数;If yes, then the refractive index of the infrared mask sample material assumed in step 8 absorption coefficient That is, the real refractive index and absorption coefficient of the infrared head cover sample material;
若不是,则返回步骤七,根据逆问题算法更新的红外头罩样品材料的折射率吸收系数重新设定红外头罩样品材料的折射率和吸收系数重新计算,直至步骤八中的目标函数值Fobj小于设定阈值ξ,得到该的红外头罩样品材料的真实折射率吸收系数 If not, go back to step 7 and update the refractive index of the infrared hood sample material according to the inverse problem algorithm absorption coefficient Re-set the refractive index and absorption coefficient of the infrared mask sample material and recalculate until the objective function value F obj in step 8 is less than the set threshold ξ, and obtain the real refractive index of the infrared mask sample material absorption coefficient
结合步骤六,目前获得了温度为Twin的样品材料的自身辐射折射率吸收系数 Combined with step 6, the self-radiation of the sample material with temperature T win is now obtained refractive index absorption coefficient
步骤十、改变样品材料的温度Twin,重复步骤二至步骤九,以获得不同温度Twin,j下,单位厚度Δ红外头罩样品材料的折射率吸收系数与自身辐射其中下标j表示第j组测量;Step 10. Change the temperature T win of the sample material, and repeat
通过计算得到单位厚度Δ红外头罩样品材料的不同方向的自身辐射 The self-radiation in different directions of the sample material of the infrared head cover per unit thickness Δ is obtained by calculation
通过不同温度Twin,j下的样品材料的折射率吸收系数与自身辐射建立不同温度单位厚度Δ红外头罩材料的辐射物性数据库;By the refractive index of the sample material at different temperatures T win,j absorption coefficient with self-radiation Establish a database of radiation properties of infrared head cover materials with different temperature unit thickness Δ;
步骤十一、利用物理离散的思想,将待测红外头罩均匀分成m层厚度为Δ的薄层;使用红外热像仪测量,待测工况下,各个红外头罩薄层的温度,记为Twin,k,其中下标k=1,2,…,m,表示第k个薄层;根据测量结果建立红外头罩的温度场,查询步骤十中建立的不同温度单位厚度Δ红外头罩材料的辐射物性数据库,得到各个红外头罩薄层的折射率、吸收系数、自身辐射分布场以及不同方向辐射分布场;进而得到红外头罩内不同位置处的折射率以及吸收系数,通过沿厚度方向叠加获取待测红外头罩的方向辐射强度以及方向发射率。Step 11. Using the idea of physical dispersion, divide the infrared head cover to be measured into m layers with a thickness of Δ evenly; use an infrared thermal imager to measure the temperature of each infrared head cover thin layer under the working condition to be measured, record is T win,k , where the subscript k=1,2,...,m represents the kth thin layer; according to the measurement results, the temperature field of the infrared head cover is established, and the infrared head with different temperature unit thickness Δ established in step 10 is queried. The radiation physical property database of the cover material is used to obtain the refractive index, absorption coefficient, self-radiation distribution field and radiation distribution field in different directions of each infrared head cover thin layer; and then the refractive index and absorption coefficient at different positions in the infrared head cover are obtained. The thickness direction superposition obtains the directional radiation intensity and directional emissivity of the infrared head cover to be tested.
有益效果:Beneficial effects:
传统的计算方法在模型层面对待测红外光学窗口进行离散,但是传统的计算方法受到待测材料厚度以及测量误差的影响。因此本发明另辟蹊径,采用在物理层面对材料进行离散,将待测材料分成数个厚度较小的薄层,对薄层的透过率以及自身辐射进行测量,并引入辐射传输反问题计算方法,对材料的辐射特性,诸如吸收系数以及折射率进行反演重建。此时,由于反演重建的模型厚度较小,因此能保证重建结果具有较高精度。The traditional calculation method discretizes the infrared optical window to be measured at the model level, but the traditional calculation method is affected by the thickness of the material to be measured and the measurement error. Therefore, the present invention takes a different approach, using discrete materials at the physical level, dividing the material to be measured into several thin layers with smaller thickness, measuring the transmittance and self-radiation of the thin layers, and introducing a calculation method for the inverse problem of radiation transmission, Inversion reconstruction of the material's radiative properties, such as absorption coefficient and refractive index. At this time, since the thickness of the model reconstructed by inversion is small, the reconstruction result can be guaranteed to have high accuracy.
同时本发明还根据红外头罩的目标特性传输机理提供了一种无需多次重复实验即可针对不同材料厚度红外头罩的透过率和自身辐射等热辐射特性测量方法,即在前期建立红外头罩传输特性量化模型时只要基于不同温度下的样品材料的折射率、吸收系数与自身辐射建立单位厚度红外头罩材料的辐射物性数据库之后就以不用重复实验即可针对不同材料厚度获得红外头罩的透过率和自身辐射等热辐射特性。从而为评估弹载条件对红外探测图像以及信息处理的影响提供关键数据支持,形成更加逼真的红外探测仿真图像,促进高超声速飞行红外探测系统的设计与优化。At the same time, the present invention also provides a method for measuring the thermal radiation characteristics such as transmittance and self-radiation of the infrared head cover with different material thicknesses without repeated experiments, according to the transmission mechanism of the target characteristics of the infrared head cover. In the quantitative model of the transmission characteristics of the head cover, as long as the radiation physical database of the infrared head cover material per unit thickness is established based on the refractive index, absorption coefficient and self-radiation of the sample material at different temperatures, the infrared head can be obtained for different material thicknesses without repeating the experiment. The transmittance of the cover and the thermal radiation characteristics such as self-radiation. Therefore, it provides key data support for evaluating the impact of missile loading conditions on infrared detection images and information processing, forms more realistic infrared detection simulation images, and promotes the design and optimization of hypersonic flight infrared detection systems.
附图说明Description of drawings
图1为高温红外头罩不同厚度材料热辐射特性高效测量系统示意图。Figure 1 is a schematic diagram of a high-efficiency measurement system for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hoods.
具体实施方式Detailed ways
下面将结合附图对本发明的工作原理进行进一步详细说明。The working principle of the present invention will be further described in detail below with reference to the accompanying drawings.
具体实施方式一:结合图1说明本实施方式,Embodiment 1: This embodiment is described with reference to FIG. 1 ,
本实施方式为一种高温红外头罩不同厚度材料热辐射特性高效测量系统,包括傅里叶变换红外光谱仪1、加热炉2、黑体炉3、温控巡检仪4和数据采集及处理系统5;This embodiment is an efficient measurement system for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hoods, including a Fourier transform
工作时,傅里叶红外光谱仪的探测镜头中心、加热炉中心以及黑体炉腔体中心设置在同一条水平线上;When working, the center of the detection lens, the center of the heating furnace and the center of the blackbody furnace cavity of the Fourier transform infrared spectrometer are set on the same horizontal line;
黑体炉,用于发射黑体红外辐射;在测量工作过程中,调整黑体炉改变黑体温度,以发射不同黑体温度下的红外辐射。The black body furnace is used to emit black body infrared radiation; during the measurement process, the black body furnace is adjusted to change the black body temperature to emit infrared radiation at different black body temperatures.
加热炉,用于红外头罩样品加热;在测量工作过程中,调整加热炉的温度,为红外头罩样品材料提供不同温度;Heating furnace, used for infrared hood sample heating; during the measurement process, adjust the temperature of the heating furnace to provide different temperatures for the infrared hood sample material;
温控巡检仪,用于检测和控制加热炉内的温度;Temperature control inspection instrument, used to detect and control the temperature in the heating furnace;
傅里叶红外光谱仪,用于获取透过红外头罩的黑体红外辐射;Fourier infrared spectrometer, used to obtain black body infrared radiation through the infrared hood;
数据采集及处理系统,用于采集傅里叶红外光谱仪的数据和温控巡检仪,并利用傅里叶红外光谱仪得到的信号计算得到材料在温控巡检仪所显示的温度下的法相光谱表观辐射强度。The data acquisition and processing system is used to collect the data of the Fourier transform infrared spectrometer and the temperature control inspection instrument, and use the signal obtained by the Fourier infrared spectrometer to calculate the normal phase spectrum of the material at the temperature displayed by the temperature control inspection instrument. apparent radiation intensity.
本实施方式中,In this embodiment,
傅里叶变换红外光谱仪为FTIR-6100型傅里叶变换红外光谱仪,主要指标有:(1)最高光谱分辨率为0.022nm;(2)扫描光谱范围为1.25-25;(3)扫描频率为20Hz;(4)信噪比为50000/1。The Fourier transform infrared spectrometer is FTIR-6100 Fourier transform infrared spectrometer. The main indicators are: (1) the highest spectral resolution is 0.022nm; (2) the scanning spectral range is 1.25-25; (3) the scanning frequency is 20Hz; (4) The signal-to-noise ratio is 50000/1.
加热炉为SGM.M6/14AE型加热炉,其为双开门加热炉,炉膛深、宽、高尺寸(mm)为230×180×150。加热功率:4KW,加热元件为硅钼棒。最高温度:1700℃,温度稳定性为±1℃,温度均匀性为±6℃,从室温升到1400度需要45-50分钟。The heating furnace is a SGM.M6/14AE heating furnace, which is a double-door heating furnace, and the depth, width and height (mm) of the furnace chamber are 230×180×150. Heating power: 4KW, the heating element is silicon molybdenum rod. Maximum temperature: 1700°C, temperature stability is ±1°C, temperature uniformity is ±6°C, and it takes 45-50 minutes to rise from room temperature to 1400°C.
黑体炉为RT1500型黑体炉,主要指标参数为:最高温度1450℃,有效发射率0.99,辐射孔径φ50,温度范围:0—1450℃范围内任意设定,温度控制精度±0.5℃,稳定度1℃/3min,升温时间:室温至1450℃不大于1小时。The black body furnace is RT1500 type black body furnace. The main index parameters are: maximum temperature 1450℃, effective emissivity 0.99, radiation aperture φ50, temperature range: arbitrarily set within the range of 0-1450℃, temperature control accuracy ±0.5℃,
具体实施方式二:Specific implementation two:
本实施方式为一种高温红外头罩不同厚度材料热辐射特性高效测量方法,包括以下步骤:The present embodiment is an efficient method for measuring the thermal radiation characteristics of materials with different thicknesses of a high-temperature infrared head cover, including the following steps:
步骤一、搭建高温红外头罩不同厚度材料热辐射特性高效测量系统。Step 1: Build an efficient measurement system for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hoods.
步骤二、初始阶段,不启动加热炉,加热炉内不放置样品,启动黑体炉,设定黑体炉温度为Tb,用傅里叶红外光谱仪获取黑体的红外辐射Lobj。
步骤三、将红外头罩样品放置于高温加热炉内加热,待红外头罩样品材料温度达到指定温度Twin且分布均匀后,用红外探测器获取透过红外头罩样品材料的红外辐射Ltot;Step 3: Place the infrared hood sample in a high-temperature heating furnace for heating. After the temperature of the infrared hood sample material reaches the specified temperature T win and is uniformly distributed, use an infrared detector to obtain the infrared radiation L tot passing through the infrared hood sample material. ;
红外头罩材料即红外光学窗口材料;Infrared head cover material is infrared optical window material;
步骤四、控制样品材料温度,使其温度保持Twin不变,改变黑体温度Tb状态下,重复步骤二与步骤三,获取多组黑体温度Tb,i状态下的红外辐射Lobj,i和Ltot,i;其中下标i表示第i次测量。Step 4: Control the temperature of the sample material to keep the temperature T win unchanged, and repeat
步骤五、当材料的温度不变时,其辐射特性参数为固定值,因此可以通过对试验结果的统计,利用最小二乘法拟合多组黑体温度Tb,i状态下的红外辐射Lobj,i和Ltot,i,进而得出温度Twin均匀分布的红外头罩样品材料的透过率τT,win和自身辐射LT,win。
步骤六、沿厚度方向将红外头罩样品等分成n层,根据能量守恒关系获得单位厚度Δ的红外头罩样品材料的表观光谱透过率自身辐射并基于获得表观法向光谱发射率 Step 6: Divide the infrared head cover sample into n layers equally along the thickness direction, and obtain the apparent spectral transmittance of the infrared head cover sample material with unit thickness Δ according to the energy conservation relationship self-radiation and based on Obtain the apparent normal spectral emissivity
表观法向光谱发射率可以通过下式计算得到The apparent normal spectral emissivity can be calculated by
式中,LT,b表示温度与单位厚度的样品材料的温度相同,即温度为Twin的黑体辐射强度。In the formula, L T,b indicates that the temperature is the same as the temperature of the sample material per unit thickness, that is, the black body radiation intensity with temperature T win .
步骤七、根据辐射传输逆问题求解算法,假设红外头罩样品材料的折射率为吸收系数为通过求解辐射传输方程计算得到该红外头罩样品材料的出射界面上任意角度的表观光谱辐射强度表观法向光谱发射率估计值以及表观光谱透过率估计值 Step 7. According to the algorithm for solving the inverse problem of radiation transmission, it is assumed that the refractive index of the sample material of the infrared head cover is The absorption coefficient is The apparent spectral radiation intensity at any angle on the exit interface of the infrared mask sample material is calculated by solving the radiation transfer equation Apparent Normal Spectral Emissivity Estimates and an estimate of apparent spectral transmittance
步骤八、将步骤六得到的红外头罩样品材料的表观法向光谱发射率与表观光谱透过率和步骤七得到的红外头罩样品材料的表观法向发射率估计值与表观光谱透过率估计值代入如下目标函数计算公式,计算得到目标函数值Fobj;Step 8. Calculate the apparent normal spectral emissivity of the infrared hood sample material obtained in Step 6 and apparent spectral transmittance and the estimated apparent normal emissivity of the infrared hood sample material obtained in step 7 and apparent spectral transmittance estimates Substitute the following objective function calculation formula, and calculate the objective function value F obj ;
步骤九、判断步骤八中的目标函数值Fobj是否小于设定阈值ξ,Step 9: Determine whether the objective function value F obj in step 8 is less than the set threshold ξ,
若是,则步骤八中所假设的红外头罩样品材料的折射率吸收系数即为该的红外头罩样品材料的真实折射率、吸收系数;If yes, then the refractive index of the infrared mask sample material assumed in step 8 absorption coefficient That is, the real refractive index and absorption coefficient of the infrared head cover sample material;
若不是,则返回步骤七,根据逆问题算法更新的红外头罩样品材料的折射率吸收系数重新设定红外头罩样品材料的折射率和吸收系数值重新计算,直至步骤八中的目标函数值Fobj小于设定阈值ξ,得到该的红外头罩样品材料的真实折射率吸收系数结合步骤六,目前获得了温度为Twin的样品材料的自身辐射折射率吸收系数 If not, go back to step 7 and update the refractive index of the infrared hood sample material according to the inverse problem algorithm absorption coefficient Re-set the refractive index and absorption coefficient values of the infrared hood sample material and recalculate until the objective function value F obj in step 8 is less than the set threshold ξ, and obtain the true refractive index of the infrared hood sample material absorption coefficient Combined with step 6, the self-radiation of the sample material with temperature T win is now obtained refractive index absorption coefficient
步骤十、改变样品材料的温度Twin,重复步骤二至步骤九,以获得不同温度Twin,j下,单位厚度Δ红外头罩样品材料的折射率吸收系数与自身辐射其中下标j表示第j组测量。Step 10. Change the temperature T win of the sample material, and repeat
通过计算可以得到单位厚度Δ红外头罩样品材料的不同方向的自身辐射(这里引入的不同方向的辐射强度是通过折射率和吸收系数计算得到的,不是测量得到的)。通过不同温度Twin,j下的样品材料的折射率吸收系数与自身辐射建立不同温度单位厚度Δ红外头罩材料的辐射物性数据库。The self-radiation in different directions of the sample material of the infrared hood per unit thickness Δ can be obtained by calculation (The radiation intensity in different directions introduced here is calculated from the refractive index and absorption coefficient, not measured). By the refractive index of the sample material at different temperatures T win,j absorption coefficient with self-radiation A database of radiation properties of infrared head cover materials with different temperature and unit thickness Δ is established.
步骤十一、利用物理离散的思想,将待测红外头罩均匀分成m层厚度为Δ的薄层。使用红外热像仪测量,待测工况下,各个红外头罩薄层的温度,记为Twin,k,其中下标k=1,2,…,m,表示第k个薄层。根据测量结果建立红外头罩的温度场,查询步骤十中建立的不同温度单位厚度Δ红外头罩材料的辐射物性数据库,可以得到各个红外头罩薄层的折射率、吸收系数、自身辐射分布场以及不同方向辐射分布场。进而可以得到红外头罩内不同位置处的折射率以及吸收系数,通过沿厚度方向叠加可以获取待测红外头罩的方向辐射强度以及方向发射率。Step 11: Using the idea of physical dispersion, divide the infrared head cover to be measured into m layers with a thickness of Δ. Using an infrared thermal imager to measure, under the working conditions to be measured, the temperature of each infrared hood thin layer is denoted as T win,k , where the subscript k=1,2,...,m represents the kth thin layer. According to the measurement results, the temperature field of the infrared head cover is established, and the radiation properties database of the infrared head cover material with different temperature unit thickness Δ established in step 10 is inquired, and the refractive index, absorption coefficient, and self-radiation distribution field of each infrared head cover thin layer can be obtained. And the radiation distribution field in different directions. Then, the refractive index and absorption coefficient at different positions in the infrared head cover can be obtained, and the directional radiation intensity and directional emissivity of the infrared head cover to be measured can be obtained by superimposing along the thickness direction.
具体实施方式三:Specific implementation three:
本实施方式为一种高温红外头罩不同厚度材料热辐射特性高效测量方法,步骤六获得单位厚度Δ的红外头罩样品材料的表观光谱透过率自身辐射的过程包括以下步骤:This embodiment is an efficient method for measuring the thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hoods. Step 6 obtains the apparent spectral transmittance of the infrared hood sample material with unit thickness Δ self-radiation The process includes the following steps:
红外探测系统(傅里叶变换红外光谱仪中的系统)所测得的红外头罩材料红外辐射Ltot是材料内表面的辐射亮度Lλ(s),由到达红外头罩试样外表面的目标红外辐射Lobj以及红外头罩自身辐射Lwin共同作用所产生,即The infrared radiation L tot of the infrared head cover material measured by the infrared detection system (the system in the Fourier transform infrared spectrometer) is the radiance L λ (s) of the inner surface of the material, which is determined by the target reaching the outer surface of the infrared head cover sample. The combined action of the infrared radiation L obj and the radiation L win of the infrared hood itself, namely
Ltot=Lobjτwin+Lwin (3)L tot =L obj τ win +L win (3)
式中,τwin为红外光学窗口的透过率。where τ win is the transmittance of the infrared optical window.
利用辐射传输方程描述目标辐射能量透过红外头罩试件的传递过程,沿辐射传输方向能量守恒,将红外头罩材料等分为n层,推演出单位厚度红外头罩材料的透过率和自身辐射与整个红外头罩的透过率和自身辐射之间存在固定的代数关系;根据所得代数关系可实现直接计算出其他不同厚度红外头罩材料在相同温度等工况下的透过率、衰减系数等热辐射特性数据。The radiation transfer equation is used to describe the transfer process of the target radiant energy through the infrared hood specimen. The energy is conserved along the radiation transmission direction. The infrared hood material is divided into n equal layers, and the transmittance and the unit thickness of the infrared hood material are deduced. There is a fixed algebraic relationship between self-radiation and the transmittance of the entire infrared head cover and self-radiation; according to the obtained algebraic relationship, the transmittance, Thermal radiation characteristic data such as attenuation coefficient.
为了实现该目的,本发明采用的技术方案是:沿厚度方向将红外头罩样品等分成n层,假定红外光学窗口的温度T分布均匀,且每一层的表观光谱透过率和自身辐射各向同性,根据辐射传输原理和能量守恒关系可以得到透过红外头罩材料第1层的总辐射为In order to achieve this purpose, the technical solution adopted in the present invention is: divide the infrared head cover sample into n layers equally along the thickness direction, assuming that the temperature T distribution of the infrared optical window is uniform, and the apparent spectral transmittance of each layer is and self-radiation Isotropic, according to the principle of radiation transfer and the relationship of energy conservation, the total radiation passing through the first layer of the infrared head cover material can be obtained as
透过第1~2层的红外辐射为The infrared radiation transmitted through the 1st to 2nd layers is
同理可得,透过第1~n层的总辐射即透过整个红外头罩的总辐射为In the same way, the total radiation passing through the first to n layers, that is, the total radiation passing through the entire infrared head cover is:
进而可以通过由能量法推演出厚度为Δ且温度为T的红外头罩材料的表观光谱透过率和自身辐射与整个红外探测红外头罩的透过率τT,win和自身辐射LT,win之间所满足的代数关系式为Furthermore, the apparent spectral transmittance of the infrared head cover material with a thickness of Δ and a temperature of T can be deduced by the energy method. and self-radiation The algebraic relationship satisfied between the transmittance τ T,win of the entire infrared detection infrared hood and the self-radiation L T,win is:
通过测量出厚度为x的红外头罩的透射率和红外辐射特性,进而获得单位厚度Δ的红外光学窗口材料的表观光谱透过率和自身辐射等热辐射特性数据。By measuring the transmittance and infrared radiation characteristics of the infrared head cover with a thickness of x, the apparent spectral transmittance and self-radiation characteristics of the infrared optical window material with a unit thickness Δ are obtained.
由此无需再多次重复实验即可根据所得代数关系直接计算出相同温度等环境工况下其他不同厚度的红外光学窗口材料热辐射特性数据,建立红外头罩传输特性量化模型。Therefore, without repeating the experiment many times, the thermal radiation characteristic data of other infrared optical window materials with different thicknesses under the same temperature and other environmental conditions can be directly calculated according to the obtained algebraic relationship, and the quantitative model of the transmission characteristics of the infrared head cover is established.
具体实施方式四:Specific implementation four:
本实施方式为一种高温红外头罩不同厚度材料热辐射特性高效测量方法,步骤七中,通过求解辐射传输方程计算得到该红外头罩样品材料的出射界面上任意角度的表观光谱辐射强度表观法向光谱发射率估计值以及表观光谱透过率估计值的过程包括以下步骤:This embodiment is an efficient method for measuring the thermal radiation characteristics of materials with different thicknesses of a high-temperature infrared head cover. In step 7, the apparent spectral radiation intensity at any angle on the exit interface of the infrared head cover sample material is calculated by solving the radiation transfer equation. Apparent Normal Spectral Emissivity Estimates and an estimate of apparent spectral transmittance The process includes the following steps:
假设各向同性介质温度为折射率为吸收系数为反射率为介质厚度为Δ;在一维条件下,当介质处于稳态条件下且不考虑介质散射时,辐射传递方程可化简为:Assuming that the temperature of the isotropic medium is The refractive index is The absorption coefficient is reflectivity The thickness of the medium is Δ; under one-dimensional conditions, when the medium is in a steady state and the scattering of the medium is not considered, the radiative transfer equation can be simplified to:
对辐射传递方程进行求解可以得到:Solving the radiation transfer equation yields:
其中,表示前向辐射,表示后向辐射;θ为前向辐射、后向辐射分别与面法线的夹角;in, represents forward radiation, Represents the backward radiation; θ is the angle between the forward radiation and the backward radiation respectively and the surface normal;
在内部边界x=0处,沿正方向传播的辐射强度包括将入射在界面上的辐射强度反射的部分,因此有:At the inner boundary x=0, the radiation intensity propagating in the positive direction includes the part that reflects the radiation intensity incident on the interface, thus:
同理,在内部边界x=L处,沿负方向传播的辐射强度也包括将入射在界面上的辐射强度反射的部分,即:Similarly, at the inner boundary x=L, the radiation intensity propagating in the negative direction also includes the part that reflects the radiation intensity incident on the interface, namely:
结合介质内沿正向、负向的辐射传递方程,以及两个边界条件,可以化简得到:Combining the radiation transfer equations along the positive and negative directions in the medium and the two boundary conditions, it can be simplified to get:
其中,in,
进而出射界面上任意角度的表观光谱辐射强度为:Then the apparent spectral radiation intensity at any angle on the exit interface is:
表观法向光谱发射率为:The apparent normal spectral emissivity is:
式中,LT,b是对应温度为T、波长为λ的黑体辐射强度;where L T,b is the black body radiation intensity corresponding to temperature T and wavelength λ;
表观光谱透过率由贝尔定律可以得到:Apparent spectral transmittance From Bell's law we can get:
由于各向同性介质温度折射率吸收系数反射率为为假设值,因此对应得到的表观法向光谱发射率表观光谱透过率即为表观光谱发射率估计值表观光谱透过率估计值 Due to isotropic medium temperature refractive index absorption coefficient reflectivity is an assumed value and therefore corresponds to the apparent normal spectral emissivity obtained Apparent spectral transmittance is the apparent spectral emissivity estimate Apparent spectral transmittance estimate
在本发明中,通过测量一种厚度的红外头罩材料热辐射传输特性,计算出单位厚度红外头罩材料的相应辐射特性,进而计算得到相同温度工况下其他不同厚度红外头罩材料的透过率等热辐射特性数据,实现了红外头罩传输量化模型的建立。In the present invention, by measuring the thermal radiation transmission characteristics of an infrared head cover material of one thickness, the corresponding radiation characteristics of the infrared head cover material per unit thickness are calculated, and then the transmission characteristics of other infrared head cover materials with different thicknesses under the same temperature condition are calculated. The thermal radiation characteristic data such as over-rate has realized the establishment of the quantitative model of infrared head cover transmission.
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。The present invention can also have other various embodiments. Without departing from the spirit and essence of the present invention, those skilled in the art can make various corresponding changes and deformations according to the present invention, but these corresponding changes and deformations are all It should belong to the protection scope of the appended claims of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210324793.7A CN114646663A (en) | 2022-03-30 | 2022-03-30 | High-efficiency measurement system and method for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hood |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210324793.7A CN114646663A (en) | 2022-03-30 | 2022-03-30 | High-efficiency measurement system and method for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hood |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114646663A true CN114646663A (en) | 2022-06-21 |
Family
ID=81994576
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210324793.7A Pending CN114646663A (en) | 2022-03-30 | 2022-03-30 | High-efficiency measurement system and method for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hood |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114646663A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115825145A (en) * | 2022-11-18 | 2023-03-21 | 哈尔滨工业大学 | Joint measurement device and inversion method for heat radiation and heat conduction parameters of high-temperature liquid molten salt |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5358333A (en) * | 1991-05-06 | 1994-10-25 | Hoechst Aktiengesellschaft | Method of and measuring arrangement for contactless on-line measurement |
CN101000264A (en) * | 2007-01-15 | 2007-07-18 | 哈尔滨工业大学 | Spectral extreme temp. measuring method |
JP2015225034A (en) * | 2014-05-29 | 2015-12-14 | 株式会社超高温材料研究センター | Measurement method of thermal diffusivity of translucent material |
CN105203494A (en) * | 2015-09-21 | 2015-12-30 | 哈尔滨工业大学 | Inversion measurement device and method for refractive index and absorption coefficient of high-temperature translucent material |
CN105223230A (en) * | 2015-09-29 | 2016-01-06 | 北京航天自动控制研究所 | A kind of infrared electromagnetic wave transparent material radiation measurement of transmission characterist method |
RU2676385C1 (en) * | 2017-08-31 | 2018-12-28 | Акционерное общество "Обнинское научно-производственное предприятие "Технология" им. А.Г. Ромашина" | Method for thermal control during heat tests of antenna domes of missiles |
-
2022
- 2022-03-30 CN CN202210324793.7A patent/CN114646663A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5358333A (en) * | 1991-05-06 | 1994-10-25 | Hoechst Aktiengesellschaft | Method of and measuring arrangement for contactless on-line measurement |
CN101000264A (en) * | 2007-01-15 | 2007-07-18 | 哈尔滨工业大学 | Spectral extreme temp. measuring method |
JP2015225034A (en) * | 2014-05-29 | 2015-12-14 | 株式会社超高温材料研究センター | Measurement method of thermal diffusivity of translucent material |
CN105203494A (en) * | 2015-09-21 | 2015-12-30 | 哈尔滨工业大学 | Inversion measurement device and method for refractive index and absorption coefficient of high-temperature translucent material |
CN105223230A (en) * | 2015-09-29 | 2016-01-06 | 北京航天自动控制研究所 | A kind of infrared electromagnetic wave transparent material radiation measurement of transmission characterist method |
RU2676385C1 (en) * | 2017-08-31 | 2018-12-28 | Акционерное общество "Обнинское научно-производственное предприятие "Технология" им. А.Г. Ромашина" | Method for thermal control during heat tests of antenna domes of missiles |
Non-Patent Citations (3)
Title |
---|
RAINE K.W.DENG: "Instrument for measuring the thermal coefficient of refractive index in the infrared", 《JOURNAL OF PHYSICS E: SCIENTIFIC INSTRUMENTS》, 31 December 1985 (1985-12-31), pages 593 - 597 * |
孙建平: "复杂环境背景下典型飞行器表面红外反射特性研究", 《激光与红外》, vol. 51, no. 4, 30 April 2021 (2021-04-30), pages 493 - 502 * |
宋敏敏;王碧云;王爽;吕弢;: "高速飞行器头罩气动热辐射传输效应计算", 上海航天, no. 04, 25 August 2016 (2016-08-25), pages 53 - 58 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115825145A (en) * | 2022-11-18 | 2023-03-21 | 哈尔滨工业大学 | Joint measurement device and inversion method for heat radiation and heat conduction parameters of high-temperature liquid molten salt |
CN115825145B (en) * | 2022-11-18 | 2023-07-11 | 哈尔滨工业大学 | High-temperature liquid molten salt heat radiation and heat conduction parameter combined measurement device and inversion method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | A method for reducing the influence of measuring distance on infrared thermal imager temperature measurement accuracy | |
CN103558243B (en) | A kind of high-speed aircraft hot surface full field deformation measure device based on optical means | |
CN103674888B (en) | High temperature semitransparent material spectrum direction apparent emissivity backstepping measurement mechanism and method | |
Ochs et al. | A novel calibration method for an infrared thermography system applied to heat transfer experiments | |
CN105203494B (en) | High temperature semitransparent material refractive index and absorption coefficient inversion measuring device and method | |
CN109655161B (en) | Target infrared integral radiation intensity testing method and device based on thermal infrared imager | |
CN105547469B (en) | The general yardstick optical field detection method of flame temperature based on microlens array and pulse laser | |
Cai et al. | Acquisition of kHz-frequency two-dimensional surface temperature field using phosphor thermometry and proper orthogonal decomposition assisted long short-term memory neural networks | |
CN110567591A (en) | A Temperature/Emissivity Inversion Method Applicable to Surface Thermal Infrared Data | |
CN114646663A (en) | High-efficiency measurement system and method for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hood | |
Wen et al. | Real-time estimation of time-dependent imposed heat flux in graded index media by KF-RLSE algorithm | |
Che et al. | Development of ReFaST pyrometer for measuring surface temperature with unknown emissivity: methodology, implementation, and validation | |
CN111581720B (en) | Method for evaluating uncertainty of temperature data of all-dimensional program control thermal test of aircraft | |
JP3568271B2 (en) | Method and apparatus for measuring thermal constant using laser flash method | |
CN113390602B (en) | Full-field pneumatic heat measurement system and method based on embedded temperature-sensitive paint | |
CN105571741B (en) | The general scale optical field detection method of flame temperature based on microlens array and continuous laser | |
CN113063819B (en) | A system and method for studying radiation characteristics of engine environmental thermal resistance coating | |
CN215953424U (en) | A Research System for Radiation Characteristics of Engine Environmental Thermal Resistance Coatings | |
CN113065275B (en) | An online prediction method for stagnant heat flow during aircraft flight | |
CN113218872B (en) | A method for simultaneous identification of multiple parameters of optical properties of high-temperature translucent materials | |
CN116183036A (en) | Method for correcting background radiation response of short wave infrared band of polarized remote sensor | |
CN117109752A (en) | Medium-length wave infrared sea surface emissivity measuring method based on measured data | |
CN113008389A (en) | Method for measuring temperature of divertor target plate of tokamak device based on single-wave filtering | |
CN211855419U (en) | A system for measuring flight attitude and surface temperature of ultra-high-speed aircraft | |
CN110132874B (en) | Multi-angle measurement-based detection device and method for optical parameter field of dispersion medium |
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 |