一种发动机尾焰温度测试装置及温度场重构方法与流程

1.本发明属于高温高速火焰温度探测技术领域，具体涉及一种发动机尾焰温 度测试装置及温度场重构方法。


背景技术：

2.研究发动机工作时喷出的高温、高速燃气射流特性，对火箭红外探测、热 防护、发动机燃烧动力性能检测等方面具有重要意义。发动机喷射火焰涉及高 温高速流动、复燃化学反应等一系列复杂的物理化学效应，火焰温度范围在 800-3000k。目前通过实验研究喷射火特性是喷射火焰研究的重要手段，但是由 于高速流场、复杂冲击波和边界层之间的相互作用使得喷射火焰温度测量较为 困难。
3.目前火焰温度测试主要有热电偶接触式测温和红外光学非接触式测温两种 方式。热电偶测温技术已经比较成熟，测试方法简单，结果直接稳定。但是目 前热电偶测温范围较窄，一般不能承受2300k以上的高温；此外，裸露在测试 环境中的热电偶丝在发动机喷射火测试中，极其容易被高速气流冲刷变形导致 短路甚至冲断，得不到有效温度数据。目前红外光学非接触测温常用的方法有 单色测温法和比色测温法，其中单色测温法受距离和测试环境影响较大，导致 测试重复性差、误差大；比色测温法采用双通道探测，精度较高，测温上限高， 但在低温下信噪比低，不适用于探测低于1300k的温度测试；此外，光学法测 温只能探测光路上的最高温度，无法得到精确位置的温度值。中国专利公开文 本cn112539939a公开了一种相变发动机尾焰温度测试装置及控制方法，但是 其测量温度的范围有一定的限制。
4.因此，对于温度上限高，流速高的发动机尾焰温度测试，目前尚无有效的 测试装置和温度场重构方法，无法准确得到发动机尾焰的温度分布。


技术实现要素：

5.有鉴于此，本发明的目的是提供一种发动机尾焰温度测试装置及温度场重 构方法，可以较为精确的得到发动机喷管外喷射火及高温燃烧产物的温度分布。
6.通过如下技术手段实现：
7.一种发动机尾焰温度测试装置，包括待测发动机、双波温度变送仪、热电 偶单点测温传感器、双波温度变送仪适配器，热电偶适配器、高频数据采集仪、 测控计算机。
8.设置5～7枚所述双色温度变送仪(优选6枚)于待测发动机紧靠喷口处的 中轴线上，用于探测靠近喷口处中轴线上的测试点温度，设置4～6枚所述热电 偶单点测温传感器(优选5枚)于待测发动机远离喷口处的中轴线上，用于测 量中轴线远端的测试点温度，再设置4～6枚所述热电偶单点测温传感器(优选 5枚)分别弧形布置在待测发动机的火焰区域，用于根据火焰形状呈弧形布置 在火焰区域而探测火焰区域的各点温度，布置时双色温度变送仪的镜头瞄准待 测点，热电偶单点测温传感器的热电偶探头与待测发动机的尾焰流速方向共线， 所述双色温度变送仪通过信号线连接到所述双色温度变送仪适配器，所有
的所 述热电偶单点测温传感器通过信号线连接到所述热电偶适配器，双色温度变送 仪适配器和热电偶适配器得到的信号通过所述高频数据采集仪连接至所述测控 计算机。
9.所述双色温度变送仪采用si和ge半导体红外光子探测器，探测中心波长 分别为0.95～0.96μm和1.470～1.478μm，测温范围为1300k-3300k，高温物体发 射出的辐射光能经凸透镜聚焦后进入所述双色温度变送仪的双色合成传感器， 产生电信号，再由放大电路转换放大并被数据采集系统接收，最后以电压响应 曲线的形式显示于计算机界面上。
10.所述热电偶单点测温传感器的测温范围为300k-2000k。
11.所述高频数据采集器设置为与发动机的点火装置同步触发。
12.进一步的，所述双色温度变送仪采用si和ge半导体红外光子探测器，探 测中心波长分别为0.957μm和1.474μm；所述热电偶单点测温传感器采用 0.02mm钨-铼热电偶丝。
13.所述双色温度变送仪(3)采用以下方法获得其电信号-温度响应关系：si 和ge两探测通道的光电转换系数为k1和k2，标定得到的具体工作波段为(λ
1a
， λ
1b
)μm与(λ
2a
，λ
2b
)μm，测量温度为t的物体时输出的电流信号为i1和i2， 则比值与温度的对应关系满足如下公式：
[0014][0015][0016][0017]
其中，c1为第一辐射常数，取0.595521
×
10-16w·
m2；c2为第一辐射常数， 取1.438775
×
10-16w·
m2；λ为波长；i
10
和i
20
分别为实验测定的si和ge两探测 通道的零点电信号。该部分以λa、λb和k为求解对象，为获得准确的唯一解则 需建立三元方程组，需要三组对应的标准温度和输出电信号；以光电转换系数 k和波段上下极限波长λa、λb为对象的三点标定法可以实现快速准确的标定求 解，以梯度速降法求解非线性方程组，以gauss积分法计算定积分，得到三个 参数的具体数值和电信号-温度响应关系。
[0018]
进一步的，探测高速火焰流场温度，需要对热电偶采取有效防护措施，防 止热电偶丝短路或损坏；所述热电偶单点测温传感器的热电偶探头包括热电偶 丝、信号输出线、热塑管、不锈钢保护壳和填充型陶瓷化硅橡胶材料，所述热 电偶丝与所述信号输出线连接，且热电偶丝和信号输出线的外部套设有一层所 述热塑管，所述不锈钢保护壳套设在最外部，所述热塑管从不锈钢保护壳内部 穿出，所述热塑管与不锈钢保护壳的间隙采用所述填充型陶瓷化硅橡胶材料进 行密封，所述填充型陶瓷化硅橡胶材料用于固定热电偶丝并保护信号输出线， 防止短路或损坏。
[0019]
作为优选，热电偶探头的不锈钢保护壳端部为缩口结构，使得高温高速流 体进入保护壳后膨胀减速，降低对热电偶丝的冲击。
[0020]
作为优选，热电偶探头的不锈钢保护壳侧面对称开有四个圆孔型出气口， 用于气
体流出，同时增加气流湍流度，使其与热电偶充分接触，得到精确的动 态温度值。
[0021]
作为优选，热电偶节点四周端面涂覆一层硅橡胶缓冲层，用于缓冲高速气 流对热电偶探头的振动和冲击。
[0022]
一种发动机尾焰温度测试装置的温度场重构方法，其操作步骤如下：
[0023]
(1)采用上述发动机尾焰温度测试装置测试得到多个测试点的单点温度 值，比色测温结果和热电偶测试结果首先经过小波分解去除噪声；
[0024]
(2)将步骤(1)得到的监测点温度值的单点温度值采用多维极值非线性 拟合法进行插值拟合，进行火焰温度场重构，得到火焰温度场分布图和温度等 值线图；
[0025]
(3)根据推进剂燃烧火焰面特征温度，在不同时刻的火焰温度场分布图中 得到火焰面位置。
[0026]
进一步的，步骤(3)所述的多维极值非线性拟合法的估计值q(x,y)满足下 式公式：
[0027]
q(x,y)＝∑a
idi2
logdi+a+bx+cy
[0028]
式中，x和y为由插值得到的点的坐标，di为(x，y)和(xi，yi)两点距离， xi、yi分别为控制点i的x、y坐标，ai、a、b、c为待拟合系数。多维极值非线 性拟合法包括两部分：a+bx+cy表示局部趋势函数，它与线性或一阶趋势面具 有相同的形状，∑a
idi2
logdi为基函数，可获得最小曲率面。
[0029]
进一步的，多维极值非线性拟合法包括两部分：a+bx+cy表示局部趋势函 数，它与线性或一阶趋势面具有相同的形状，∑a
idi2
logdi为基函数，可获得最 小曲率面；系数由以下线性方程组确定：
[0030][0031][0032][0033][0034]
式中，n为控制点数目，fi为已知控制点坐标，系数计算由n+3个方程联立 求解。
[0035]
进一步的，步骤(1)中，对同一型号发动机进行5～15次试验，将每次试 验得到的数据整体作为多个测试点的单点温度值。进一步的，测试点增多得到 插值拟合结果越接近真实工况，但为了尽可能减小测试仪器对火焰流场的影响， 每次测试不宜架设过多传感器，应结合高速摄影图像和火焰形状预测结果进行 热电偶传感器布置，对同一型号发动机进行多次试验，将各发实验数据进行统 一处理，得到均一化的温度场分布。
[0036]
温度场重构结果的精确度用由sse(和方差、误差平方和)和r-square(确定 系数)两个指标来评价，较佳的，当sse小于10-20
，r-square大于0.98时，拟合 结果可以反映真实工况。
[0037]
本发明具有如下有益效果：
[0038]
本发明提供的发动机尾焰温度测试装置，通过对各个部件进行具体设置以 及设
置测温过程中各部件的位置关系以及测温部件的具体结构设置，使得其测 温范围可以涵盖273-3300k，能够实现各类火箭尾焰各点温度的实时准确测试， 测试装置稳定可靠，重复性强。
[0039]
本发明提供的与测试装置对应的温度场重构方法，通过使用本发明特定设 置的发动机尾焰温度测试装置，使得温度场重构方法能够更好的实现，通过对 各个步骤进行具体控制和设置，使得关键测试点能够重构火箭尾焰温度场分布。 同时通过本发明特定的温度场重构方法的各步骤设置，能够为发动机的火焰导 流槽结构和热防护设计提供有效数据支持。
附图说明
[0040]
图1是本发明一种实施方式的发动机尾焰温度测试装置示意图。
[0041]
图2是本发明的发动机尾焰温度测试装置的热电偶结构示意图。
[0042]
图3是本发明实施例提供的热电偶安装架设示意图。
[0043]
图4是本发明实施例提供的温度测点布置图。
[0044]
图5是本发明实施例提供的热电偶测试数据。
[0045]
图6是本发明实施例提供的三维极值非线性差值拟合结果。
[0046]
图7是本发明实施例提供的火焰温度场重构等值线图。
[0047]
其中：1-发动机，2-热电偶单点测温传感器，3-双波温度变送仪，4-双波 温度变送仪适配器，5-热电偶适配器，6-高频数据采集仪，7-测控计算机，8
‑ꢀ
喷射火焰，9-热电偶丝，10-出气口，11-硅橡胶缓冲层，12-填充型陶瓷化硅橡 胶材料，13-不锈钢保护壳，14-信号线。
具体实施方式
[0048]
下面结合附图并举实施例，对本发明进行详细描述。
[0049]
本发明主要解决的技术问题是提供一种发动机尾焰温度测试装置及温度场 重构方法，可以较为精确的得到发动机喷管外喷射火及高温燃烧产物的温度分 布。
[0050]
本实施例提供的发动机尾焰温度测试装置如图1所示，包括待测发动机 1，六枚双波温度变送仪3，用于探测靠近喷口处中轴线上的六个测试点温度， 五枚热电偶单点测温传感器2，测量中轴线远端的五个测试点温度，另外五枚 热电偶单点测温传感器2，如图1所示的根据火焰形状呈弧形布置在火焰区 域，探测各点温度。热电偶安装架设示意图如图3所示，布置时热电偶探头与 尾焰流速方向共线。六枚双波温度变送仪3通过信号线连接到双波温度变送仪 适配器4，十枚热电偶单点测温传感器2通过信号线连接到热电偶适配器5；双 波温度变送仪适配器4和热电偶适配器5得到的信号传输至高频数据采集仪6， 再连接至测控计算机7。
[0051]
本实施例采用的双波温度变送仪3测温范围为1300-3300k，用于探测靠近 推进器出口处的高温火焰，双波温度变送仪是利用两个很窄的相近波段测量同 一物体，选取短波段信号与较长波段信号的比值，该比值随着温度升高而增大。 设热力学温度为t的非黑体物质在同一点的波长为λ1，λ2下的单色辐射出射度 分别为m(λ1,t),m(λ2,t)，则两波长处辐射功率比值r(t)为：
[0052][0053]
选取的非常接近两个波段，认为被测物体在该波段内发射率无变化，即 ε(λ1,t)≈ε(λ2,t)时，则其发射率和气体吸收对两个波段信号的衰减相同，其比值 不变：
[0054][0055]
本实施例采用si和ge双通道红外光子探测器，对应探测中心波长分别为 0.957μm和1.474μm，根据上述关系，读出两通道信号比值就可以计算得到待测 物体的真实温度。
[0056]
本实施例采用的热电偶单点测温传感器2采用0.02mm钨铼热电偶丝，测 温范围300k-2000k，用于探测中轴线远端温度和外围火焰温度。
[0057]
本实施例的热电偶单点测温传感器如图2所示。探测高速火焰流场温度， 需要对热电偶采取有效防护措施，防止热电偶丝短路或损坏。热电偶丝9和信 号输出线14外部覆盖一层热塑管15，置于不锈钢保护壳13内部，热塑管与不 锈钢保护壳连接间隙采用填充型陶瓷化硅橡胶材料12进行密封，用于固定热 电偶丝9和保护后面的信号线14；保护壳13端部为缩口结构，使得高温高速 流体进入保护壳后膨胀减速，减小对热电偶的冲击力；保护壳13侧面对称开 有四个圆孔状出气口10，用于气体流出，同时增加气流湍流度，使其与热电偶 充分接触，得到更为精确的动态温度值，热电偶四周端面涂覆一层硅橡胶缓冲 层11，用于缓冲高速气流对热电偶探头的振动和冲击。
[0058]
由于采用高频数据采集器6，而测控计算机7内存及读写速度有限，需设 置同步触发装置，与发动机点火装置同步触发。
[0059]
为了尽可能减小测试仪器对火焰流场的影响，每次测试不宜架设过多传感 器，应结合高速摄影图像和火焰形状理论预测结果进行热电偶传感器布置。本 实施例对同一型号发动机进行多次试验，将各发实验数据进行统一处理，得到 均一化的温度场分布，本实施例的四次实验测点布置情况如图4所示。
[0060]
测试得到多个测试点的单点温度值，比色测温结果和热电偶测试结果首先 经过小波分解去除噪声，再转换为对应温度值，由于火焰存在明显波动，读取 数据平稳段的中值作为监测点温度值，实施例中一个热电偶测得原始数据曲线 和小波降噪转换后的温度曲线分别如图5、图6所示。得到的众多单点温度值 采用多维极值非线性拟合法进行插值拟合，进行火焰温度场重构。多维极值非 线性拟合法是一种差值拟合算法，建立一个通过控制点的面，并使所有点的坡 度变化最小，即多维极值非线性拟合法以最小曲率面拟合控制点。多维极值非 线性拟合法的估计值由下式计算：
[0061]
q(x,y)＝∑a
idi2
logdi+a+bx+cy
[0062]
式中，x和y为由插值得到的点的坐标，di为(x，y)和(xi，yi)两点距离， xi、yi分别为控制点i的x、y坐标。
[0063]
多维极值非线性拟合法包括两部分：a+bx+cy表示局部趋势函数，它与线 性或一阶趋势面具有相同的形状，∑a
idi2
logdi为基函数，可获得最小曲率面。 有关系数可以由以下线性方程组确定：
[0064]
[0065][0066][0067][0068]
式中，n为控制点数目，fi为已知控制点坐标，系数计算需要n+3个方程联 立求解。
[0069]
多维极值非线性拟合法差值拟合的结果由sse(和方差、误差平方和)和 r-square(确定系数)两个指标来评价。sse为拟合数据与原始数据对应点的误差 的平方和，计算公式如下：
[0070][0071]
sse越接近于0，表征拟合结果更好，数据预测结果越准确。
[0072]
r-square是通过数据的变化来表征拟合的好坏，由ssr和sst两个参数计 算得到，其中ssr为预测数据与原始数据均值之差的平方和，sst为原始数据 与其均值之差的平方和，计算公式如下：
[0073]
r-square＝ssr/sst
[0074][0075][0076]
由上式可知，r-square取值范围为[0,1]，越接近于1，表征方程的变量对y 的解释能力越强，数据拟合结果也较好。
[0077]
本实施例得到的三维极值非线性差值拟合结果和火焰温度场重构等值线图 如图6和图7所示，其中图7中横坐标z表示距离发动机尾部端面的轴向距离， 纵坐标r表示距离发动机尾部端面中心的径向距离。拟合结果的sse值为 1.2648e-24，r-square值为0.9999，表明数据拟合结果良好，该方法可以准确进 行火焰温度场重构。
[0078]
综上所述，以上仅为本发明的较佳实施例，并非用于限定本发明的保护范 围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均 应包含在本发明的保护范围之内。