关键词:
陶瓷基复合材料
各向异性导热系数
概率性热分析
随机有限元
温度场波动
涡轮叶片
摘要:
陶瓷基复合材料(Ceramic Matrix Composite,CMC)作为提升未来航空发动机高温部件耐温能力的重要解决途径,相关CMC高温部件的热分析及冷却结构设计需求得到越来越多研究工作的关注。然而,CMC材料由于非均质、非均匀的内部结构特征,其热物性参数呈现各向异性及分散性,需要建立适用于CMC材料高温部件的热分析方法。本文针对2.5维编织结构陶瓷基复合材料涡轮叶片,重点考虑CMC材料内部结构参数和各向异性导热系数的分散性,统计分析了细观结构几何特征参数的概率分布,建立了各向异性导热系数参数化全尺寸预估模型,分析了多种几何特征的分散性对各向异性导热系数的影响规律;开发了CMC涡轮叶片的随机有限元温度场计算程序,探究了各向异性导热系数分散性对CMC涡轮叶片温度场概率分布特征的影响;最终,通过开展CMC材料平板温度场测试试验及CMC涡轮叶片综合冷效试验,验证了本文所建立的概率性热分析方法精度。
研究中首先针对2.5维编织结构CMC材料,应用X-CT/SEM扫描技术进行细观结构识别,批量获取样本照片,通过样本统计的方法建立涵盖几何特征分散性的数据库,并计算概率分布函数;针对所提取的几何特征,建立了反应细观结构分散性的参数化全尺寸模型,应用Monte-Carlo模拟的方法,探究了编织CMC材料细观结构几何特征参数的随机性对材料各向异性导热系数概率分布的影响规律,并从导热机理的角度分析其影响因素,通过正交试验验证了规律的合理性。研究结果表明,基体和纤维束的体积分数发生改变是影响等效导热系数波动的主要原因;研究中考虑了六种几何特征的影响,其中编织角的分散性对等效导热系数波动的影响最大,为k~N(2.912,0.0312);
在此基础上,考虑CMC涡轮叶片表面曲率的影响,引入主坐标和计算坐标的空间转换关系,编写了考虑各向异性导热系数概率分布的CMC涡轮叶片Neumann随机有限元温度场计算程序,分析了材料结构分散性导致的导热系数概率分布对CMC涡轮叶片温度场波动的影响规律,并与传统Monte-Carlo模拟的温度场计算结果进行对比,探究了Neumann随机有限元温度场计算的优势。研究结果表明,在导热系数波动下,叶片前缘点温度波动较大,在导热系数k11=9W/(m·K),变异系数为σ/μ=0.1时,前缘点温度波动标准差为9.06K,波动超过均值18.12K的概率约为8%;随着导热系数的增大,叶片高温敏感区域也逐渐变大,当变异系数σ/μ=0.01时,叶片高温敏感区域由ΔShot=0.41%增大为ΔShot=21.61%;Neumann随机有限元的计算时间为1295.78s,相比Monte-Carlo法计算效率提升7.7%;
随后开展了复合材料涡轮叶片综合冷效试验及平板验证试验,获得了不同冷气/燃气流量比下,CMC气冷涡轮叶片温度场分布规律及综合冷却效果,同时探究了不同主流燃气温度下,CMC平板的温度场分布规律,并基于此试验工况进行仿真计算。研究结果表明,随着冷气/燃气流量比的增大,叶片表面综合冷却效率较高区域面积逐渐增大,表现为沿气膜孔逐渐向外扩散的趋势,试验相对误差小于6.5%,验证了温度场预估精度;通过CMC平板验证试验与计算结果的比对,三种工况下的相对误差分别为δ1=5.031%、δ2=5.141%、δ3=3.812%,均小于10%,验证本文所建立概率性热分析方法的合理性。