全文HTML
--> --> -->着火理论研究主要集中在着火温度及着火延迟时间两个关键参数的求解. Frank-Kamenetskii[3]在谢苗诺夫着火热理论基础上建立了应用于非均相燃烧的着火理论; Merzhanov[4]提出了金属颗粒的着火理论模型, 并探讨了在不同气氛下临界参数的变化和外热源的变化; Elsayed[5]通过计算发现, 谢苗诺夫数的微小改变对临界着火温度的影响显著; Khaikin和Bloshenko[6]以及Aldushin等[7]建立了氧化层厚度与着火温度之间的临界关系模型, 计算发现随着颗粒半径的增加, 在8 μm左右会出现温度突变点; Chernenko等[8]在上述关系模型的基础上, 对钛颗粒的着火边界进行了理论计算, 同时结合Evans等提出的氧化机制得到了钛在空气环境下的氧化动力学方程; Shafirovich等[9]建立了钛颗粒的初始尺寸与燃烧时间的关系模型; Yuan等[10]针对纳米级钛颗粒建立了着火的数值模型; Breiter等[11]提出了扩散机理控制钛着火模型; Bolobov[12-14]实验研究了破裂条件下钛合金着火的微观机理, 并且基于非均相反应的热自燃理论建立了着火理论模型; Ouyang等[15]和弭光宝等[2,16,17]采用摩擦点燃方法开展了钛基材料的阻燃机理研究, 并初步分析了与等效压力的关系.
上述理论研究对于深入认识钛的摩擦着火现象产生了重要影响, 但受限于着火物化过程复杂性及其分析工具方法, 将摩擦作为外热源的着火理论研究鲜见报道, 尚没有成熟的理论模型给出满意的描述[2], 目前仍在不断探索. 因此, 本文在非均相着火理论基础上, 建立钛合金的摩擦着火模型, 推导出着火温度及着火延迟时间的理论计算公式, 并对摩擦系数、氧浓度、流速、接触半径以及阻燃层等因素的影响规律进行分析和讨论. 这对于发展钛合金燃烧基本理论及推动防钛火技术在先进航空发动机上应用具有重要意义.
2.1.物理模型
如图1所示, 在一定气流条件下, 转子与静子在环形区域接触后通过转子件的高速旋转产生摩擦热, 摩擦接触区(静子正面中心区域)局部温度急剧升高, 当反应区(靠近静子背面中心区域)达到临界着火温度时, 钛合金发生着火, 其物理模型假设如下:![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-1_mini.jpg)
Figure1. Schematic diagram of the model of friction between rotor and stator.
1)着火首先发生在钛合金反应面上, 摩擦面的热量通过反应区传导至反应面, 由于此区域厚度非常薄, 毕渥数趋近零, 因此忽略反应区内部热阻, 此时反应区温度与反应面、摩擦面近似相等;
2)气流流向与静子平行, 且流速不变, 反应区的散热通过与气流之间对流换热形式实现, 忽略热辐射的影响;
3)反应区表面满足单分子吸附模型;
4)放热反应为理想状态下钛与氧气充分燃烧产生二氧化钛的单步反应;
5)摩擦系数受温度影响, 随着温度的升高线性下降.
2
2.2.数学模型
32.2.1.着火温度计算公式推导
根据着火热理论, 通过求解生热平衡方程(生热速率与散热速率相等且对温度的变化率相等)得出临界着火温度, 能量方程由生热项QG与散热项QL组成, 数学表达如下:首先, 对生热项方程进行推导. 与一般非均相燃烧的生热平衡方程所不同的是, 生热项考虑了摩擦生热项QF, 由于在高温环境下, 摩擦系数是温度变化的函数, 因此QF对方程(1)和方程(2)都会产生影响.
反应区的QF可以通过库仑定律求出, 即
令
QH由阿伦尼乌斯方程获得
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M1.png)
单位时间内反应区的增重w根据如下公式[13,17]得到
令(11)式中包含吸附系数的常数项为
根据对流热交换系数计算公式[20]
强迫对流热交换系数按照经验公式计算:
最后, 推导整体方程组. 将(8)式、(13)式和(14)式代入(1)式和(2)式, 且令
联合(19)式和(20)式, 消去指数项得到
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M2.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M3.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M4.png)
为了进一步分析着火过程摩擦生热的影响, 根据上述推导公式绘制生热量QG与瞬时温度T的关系曲线, 如图2所示. 从图2可见, 从室温(288 K)开始, 生热量随着温度的升高而降低, 这是因为在较低温度条件下, 生热项主要受到摩擦热的影响, 而摩擦系数随着温度的提高而下降使得总体生热量减少; 反应热随着温度的升高呈指数增长, 在低温下对整体生热量影响较小, 在温度达到总生热量最低点之后占据生热量的主导地位.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-2_mini.jpg)
Figure2. Curves of heat generation in T-Q diagram.
因此, 定义Tg为dQG/dT = 0的临界生热温度, 其物理意义是生热过程中热量随温度的变化量为0. 当温度小于Tg时, 生热过程中热量随温度的变化量由摩擦热QF主导; 当温度大于Tg时, 生热过程中热量随温度的变化量由反应热QH主导, 如图2所示(摩擦系数采用YSZ材料, 正应力N为2650 kPa). 于是, 由(20)式得到Tg的计算公式, 即
3
2.2.2.着火延迟时间计算公式推导
根据能量方程[21]为了便于计算, 首先定义无因次温度[22,23]
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M5.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M6.png)
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-3_mini.jpg)
Figure3. Variations of θ within the 0?2 interval: (a) Comparison of numerical and approximate solutions; (b) errors between analytical and numerical solutions.
采用二次近似后, (28)式变为
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M7.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/21-20200304_M8.png)
1) 存在两个实根, 这种情况下生热量总是小于散热量, 因此不会发生着火现象, 此时
符号 | 物理意义 | 单位 | |
λ | 热传导系数 | W·m–1·K–1 | |
cp | 比热 | J·kg–1·K–1 | |
q | 单位反应热 | MJ·kg–1 | |
k | 指前因子 | kg·m–2·s–1 | |
E | 激活能 | kJ·mol–1 | |
f | 摩擦应力 | N·m–2 | |
τ | 着火延迟时间 | s | |
N | 接触应力 | N·m–2 | |
φ | 厚度 | m | |
d | 试验舱直径 | m | |
R2 | 接触面内径 | m | |
R1 | 接触面外径 | m | |
ω | 转子角速度 | r·min–1 | |
ci | 氧浓度 | 100% | |
v | 流速 | m·s | |
QG | 生热项 | J·s–1 | |
QH | 反应热项 | J·s–1 | |
QC | 对流散热项 | J·s–1 | |
QR | 辐射散热项 | J·s–1 | |
μ | 摩擦系数 | / | |
Sr | 反应区面积 | m2 | |
a | 吸附系数 | / | |
α | 总传热系数 | W·m–2·K–1 | |
Nu | 努塞尔数 | / | |
Pr | 普朗特数 | / | |
Re | 雷诺数 | / |
表A1术语表
TableA1.Nomenclature
2
2.3.边界条件
边界条件的设置主要包括转子与静子着火反应区的几何尺寸参数、环境条件以及材料的热物性参数等. 由于TC11钛合金(Ti-6.5 Al-3.5 Mo-1.5 Zr-0.3 Si, wt%)是我国目前航空发动机上用量最大的钛合金[24], 材料成熟度高, 其物理性能和机械性能等理论建模用基础数据积累丰富[25-27], 所以转子与静子材料选择接近使用状态的该类型钛合金作为研究对象, 对比材料为带ZrO2阻燃层的钛合金, 室温下合金组织主要含有α相和β相. 钛合金及其阻燃层的初始摩擦系数F、温度系数项B分别为1.11和0.97, –0.12 × 10–4和–0.45 × 10–4, 其他热性能参数[13,16,25-27]和环境参数如表1和表2所列. 其中, 带阻燃层钛合金的反应热、指前因子和激活能等边界条件的确定与不带阻燃层的钛合金相同, 主要原因在于二者的摩擦着火源均为静子中心反应区的钛合金, 即: 当摩擦生热达到临界条件时, 着火不会首先发生在转子与静子接触面的阻燃层上, 而是首先发生在静子背面中心反应区附近的新鲜钛合金表面上, 如图4所示.性能材料 | 密度 ρ/kg·m–3 | 比热容 cp/J·kg–1·K–1 | 反应热 q/MJ·kg–1 | 指前因子 k/kg·m–2·s–1 | 激活能 E/kJ·mol–1 | 导热系数 λ/W·m–1·K–1 | 吸附系数 a/MPa–0.5 |
钛合金 | 4500 | 493 | 33 | 4.2 | 44.5 | 17.8 | 0.52 |
阻燃层 | 5600 | 560 | — | — | — | 0.62 | — |
表1材料热物性参数[13,25]
Table1.Thermal property parameters of materials.
初始摩擦正应力N/kPa | 静子厚度φ/m | 特征长度d/m | 内半径R2/m | 外半径R1/m | 角速度ω/r·min–1 | 氧浓度ci/% | 流速v/m·s–1 |
265 | 0.002 | 0.016 | 0.002 | 0.004 | 5000 | 21 | 50 |
表2模型的初始边界条件[13,16]
Table2.Initial boundary conditions of the model.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-4_mini.jpg)
Figure4. In-situ observation of friction ignition process under critical condition (I, mirror image of ignition of reaction zone; II, spark of micro-bump near the centre hole of reaction zone): (a) 0.05 s, (b) 0.10 s, (c) 0.15 s, (d) 0.20 s, (e) 0.25 s, (f) 0.5 s after the rubbing ends.
2
3.1.摩擦系数的影响
图5是摩擦生热项中正应力N以及摩擦系数项中温度系数B对临界生热温度Tg影响的曲线. 可见, Tg随着B的升高而近似线性降低, 当静子为钛合金(B = –0.12 × 10–3)时, Tg为444 K, 当静子为带阻燃层的钛合金(B = –0.45 × 10–3)时, Tg为512 K, 二者的变化范围约为60 K, 如图5(a)所示; 当B增大时, 摩擦接触区的瞬态温度对摩擦系数的影响下降, 导致T-QG曲线中摩擦生热项随温度的变化减小, 温度到达Tg之前的斜率逐渐趋近于0, 如图5(b)所示.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-5_mini.jpg)
Figure5. (a) Effect of B on Tg; (b) effect of B on T-QG diagram.
图6为初始摩擦系数F及温度系数B与临界着火温度Tcr的关系, 选用26.5 kPa和2.65 kPa两个数量级的摩擦正应力N进行对比. 从图6(a)可见, 随着F的不断增大, Tcr呈线性下降趋势, N对F-Tcr曲线的初值及变化率均有较为明显的影响. 当静子为钛合金(F = 1.11)时, 对应N1和N2的Tcr分别为633 K和829 K; 当静子为带阻燃层的钛合金(F = 0.97)时, 对应N1和N2的Tcr分别为679 K和834 K, 后者比前者分别提高了46 K和5 K. 类似地, B对Tcr的影响如图6(b)所示, 可见Tcr同样随着B的不断增大而呈线性下降.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-6_mini.jpg)
Figure6. (a) Effect of F on Tcr; (b) effect of B on Tcr.
图7(a)为初始摩擦系数F与着火延迟时间τ的关系, 当静子为钛合金时, 对应N1和N2的τ分别为2.36 s和2.5 s; 当静子为带阻燃层的钛合金时, 对应N1和N2的τ分别为2.38 s和2.51 s, 变化范围在0.02 s以内; 类似地, B对τ的影响如图7(b)所示, 可见τ同样随着B的增大而呈线性规律下降, 但总体上即使将N提高至26.5 kPa, 其变化范围也在毫秒级别, 可以忽略不计.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-7_mini.jpg)
Figure7. (a) Effects of F on τ; (b) effects of B on τ.
综上可知, 摩擦系数F对Tcr的影响较大, 对τ的影响较小; 通过提高F和B可使摩擦生热明显增加, 从而使着火温度下降, 如当N = 26.5 kPa时, 摩擦系数每下降0.1, Tcr随之提高35 K. 因此, 在进行阻燃层设计时, 应尽量减小F和B的数值, 同时在实际应用中减少零件之间的摩擦接触应力, 从而降低钛火发生风险. 阻燃钛合金摩擦着火过程生成的948 K的低熔点相, 使摩擦系数大幅降低而显著提高阻燃性能的实验结果[17]进一步验证了摩擦系数对阻燃设计的重要影响.
2
3.2.氧浓度和流速的影响
根据摩擦着火的物理模型, 随着生热过程的持续发生, 摩擦接触面的瞬态温度迅速升高, 当反应区形成的微凸体突然出现在含氧环境(氧浓度ci)时, 氧的吸附起主导作用, 反应速率由ci的数值决定, 当ci达到某一临界值c0时, 反应速率急剧增大, 即发生着火, 于是根据前述计算模型可以得到ci与临界着火温度Tcr之间的关系曲线, 如图8所示.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-8_mini.jpg)
Figure8. Effects of oxygen concentrations on Tcr.
可见, Tcr随着氧浓度的提高而不断降低, 其下降的变化率在氧浓度增大至21%后趋于平缓, 此时氧浓度每提高0.01%时, Tcr降低约15 K. 在空气环境(ci = 21%), 钛合金及其阻燃层的Tcr分别为1109 K和1281 K, 阻燃层使钛合金的Tcr提高约172 K; 当ci提高至42%时, 二者的Tcr分别为896 K和1033 K, 阻燃层使Tcr提高约137 K; 当ci提高至100%时, 二者的Tcr分别为686 K和796 K, 阻燃层使Tcr提高约110 K. 如图9所示, 在静子接触区的背面(虚线红色圆圈附近)预先焊接连接两根K型热电偶, 然后对钛合金阻燃层进行摩擦着火实验(ci提高至42%), 得到Tcr的平均值为1126 K, 实验结果与理论计算值的相对误差为8.3%, 二者相一致, 从而验证了着火模型的可靠性. 因此, 将临界着火氧浓度的实验值作为模型输入参数, 可以对不同钛合金静子的Tcr进行预测评估.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-9_mini.jpg)
Figure9. In-situ observation of friction ignition reaction zone: (a) 0.15 s, (b) 0.25 s, (c) 0.5 s, (d) 1.0 s after the rubbing ends.
流速也是影响钛合金着火的重要环境因素, 比如航空发动机钛火主要发生在100—310 m/s, 当超过300 m/s以后, 钛火会发生闪燃, 在着火瞬间就发生熄灭[28]. 因而在理论计算过程中选取流速的数值在10—310 m/s范围. 图10是不同流速条件下临界着火温度的变化曲线. 从图10可见, 随着流速的不断提高, 着火温度呈上升趋势. 流速从10 m/s开始, 每提高1 m/s着火温度提升10 K左右, 而达到50 m/s之后, 着火温度的变化趋于稳定, 每提升1 m/s着火温度提升5 K左右, 当流速提高至310 m/s时, 钛合金及其阻燃层的着火温度分别达到1170 K和1339 K.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-10_mini.jpg)
Figure10. Effects of flow velocities on Tcr.
为了对氧浓度及流速这两个重要参数进行综合分析, 分别建立了钛合金及其阻燃层的氧浓度-流速着火边界图, 如图11和图12所示. 可见, 着火温度随着氧浓度的提升而下降, 在氧浓度小于21%范围内, 其着火温度变化显著, 钛合金着火温度从2500 K下降到约1000 K, 而带阻燃层的钛合金的着火温度从2500 K下降到1200 K左右; 同时着火温度随着流速的提升而上升, 在50 m/s以内变化最为剧烈, 钛合金的着火温度从700 K上升到900 K左右, 而带阻燃层的钛合金的着火温度从900 K上升到1100 K左右; 带阻燃层的钛合金的着火温度比钛合金平均高出200 K左右, 进一步说明摩擦系数及摩擦系数变化率对钛合金的着火温度具有重要影响.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-11_mini.jpg)
Figure11. Flow velocity-oxygen concentration critical condition of titanium alloy.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-12_mini.jpg)
Figure12. Flow velocity-oxygen concentration critical condition of titanium alloy with flame retardant layer.
2
3.3.摩擦接触半径的影响
摩擦实验中, 通过改变外径R1来改变接触面积, 接触面积对摩擦生热、反应产热、表面散热都会产生影响. 在谢苗诺夫非均相热自燃公式中, 由于反应产热与表面散热的面积相同, 所以在理论计算中相互抵消, 而本文建立的摩擦着火模型中, 考虑了单位摩擦应力的变化, 因此摩擦接触面积的变化会对着火温度产生明显影响, 如图13所示. 可见, 随着接触半径的增加, 着火温度呈现指数下降, 当接触半径增加至为0.007 m时, 钛合金及其阻燃层的着火温度分别为899 K和988 K, 而当继续提高接触半径后, 着火温度的变化趋于稳定, 每提高0.0001 m着火温度升高1 K左右, 当接触半径达到0.09 m时, 钛合金及其阻燃层的着火温度分别为891 K和959 K. 说明在开始阶段着火温度主要受反应产热的影响, 而随着半径的增加受到散热的影响逐渐增大.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-13_mini.jpg)
Figure13. Effect of friction contact radius on Tcr.
2
3.4.阻燃层的影响
为进一步揭示阻燃层对钛合金着火延迟时间的影响, 将上述理论模型计算出的钛合金及其阻燃层的着火温度(Tcr分别为1109 K和1281 K)作为边界条件设定参数, 采用摩擦着火热源模型[29]对带阻燃层的钛合金静子摩擦着火过程的温度场进行数值计算分析, 其流程如图14所示.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-14_mini.jpg)
Figure14. Flow chart of simulation
图15和图16分别为摩擦着火过程带阻燃层的钛合金静子背面温度场和接触面温度变化的计算结果(阻燃层厚度δ为0.4 mm、初始摩擦应力N为265 kPa). 由图15和图16可见, 随着摩擦接触时间的延长, 接触面的温度不断升高. 当摩擦接触时间小于1.7 s时, 由于接触面积小、摩擦应力高, 导致升温速率达到400 K/s; 在摩擦接触时间为1.7 s时, 接触面的温度在热累积的作用下上升至1100 K, 而由于阻燃层的存在, 此时静子背面的温度仅为433 K, 接触面与背面的温差达到667 K, 最高温度出现在中心孔反应区; 在摩擦接触时间为7 s时, 接触面的最高温度达到1483 K, 此时静子背面的温度为786 K; 而当摩擦接触时间为14.2 s时, 接触面的最高温度为1963 K, 此时静子背面中心反应区的温度已经达到了Tcr, 即着火延迟时间约为14.2 s.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-15_mini.jpg)
Figure15. Temperature field of back surface of titanium alloy stator with flame retardant layer.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-16_mini.jpg)
Figure16. Temperature history of back surface and contact surface of titanium alloy stator with flame retardant layer.
同样地, 可以得到摩擦着火过程不带阻燃层的钛合金静子背面及摩擦接触面的温度变化情况, 如图17所示. 由图17可见, 由于摩擦升温过程传热比较均匀, 静子背面及摩擦接触面的升温近似呈线性增加, 升温速率约为10 K/s, 摩擦接触面与背面的温差保持在55 K左右. 比如, 当摩擦接触时间为11.3 s时, 摩擦接触面的最高温度升至1164 K, 背面的最高温度为1110 K, 二者相差54 K, 而此时静子中心反应区的温度已经达到了Tcr, 即着火延迟时间约为11.3 s.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-17_mini.jpg)
Figure17. Temperature history of back surface and contact surface of titanium alloy stator without flame retardant layer
通过比较静子摩擦接触面与背面温差可知, 带阻燃层的钛合金静子的温差平均值为637 K, 而不带阻燃层的钛合金静子的温差约为55 K, 前者比后者高12倍, 如图18所示. 说明带阻燃层的钛合金静子摩擦生热量主要累积在阻燃层, 而传输到静子背面的热量很少, 即大幅度降低了反应区的瞬态温升, 从而在着火延迟时间内降低了钛合金着火的概率, 起到了阻燃作用. 需要说明的是, 摩擦过程应避免出现开裂和剥落等问题, 阻燃层应具有比较致密的结构, 且与钛合金基体之间具有良好的界面结合.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/21/PIC/21-20200304-18_mini.jpg)
Figure18. Difference in temperature between friction contact surface and back surface of stator.
感谢圣彼得堡国立矿业大学机械工程系Болобов В教授的讨论.
附录A术语表