基于最优输运无网格方法的热障涂层氧化层生长模拟  

Simulation of oxide growth in thermal barrier coating based on optimal transport meshless method

作　　者：许洪斌 樊江[1] 荆甫雷 廖祜明 黎波 樊宗岳 XU Hongbin;FAN Jiang;JING Fulei;LIAO Huming;LI Bo;FAN Zongyue(School of Energy and Power Engineering,Beihang University,Beijing 102206,China;Aero Engine Academy of China,Aero Engine Corporation of China,Beijing 101304,China;College of Engineering,Peking University,Beijing 100871,China;Yunyi Supercomputing(Beijing)Software Technology Company,Beijing 100027,China)

机构地区：[1]北京航空航天大学能源与动力工程学院,北京102206 [2]中国航发中国航空发动机研究院,北京101304 [3]北京大学工学院,北京100871 [4]云翼超算(北京)软件科技有限公司,北京100027

出　　处：《航空动力学报》2022年第10期2104-2111,共8页Journal of Aerospace Power

基　　金：国家自然科学基金(51905510)。

摘　　要：在自研的最优输运无网格(OTM)框架下,针对热障涂层氧化层(TGO)生长力学等效模型开发了厚度增长算法,模拟TGO各向异性氧化生长过程。利用该方法以TGO层典型的转折段为对象,研究了热循环载荷作用下的TGO层生长过程中应力、位移随周期载荷的变化规律。仿真结果与热障涂层试验结果吻合。结果表明:该方法能很好模拟模拟界面生长过程中出现的褶皱现象,与有限元相比单元变形均匀,适合障涂层TGO层生长过程的数值模拟。热障涂层应力最大位置主要出现在凸起位置,凸起处的横向变化具有加剧TGO层界面大变形的趋势。Under the self-developed optimal transport meshless(OTM) framework,a thickness growth algorithm was developed and the anisotropic oxidation growth process of thermal oxide layer(TGO) was simulated.Using this method,the typical transition section of TGO layer was taken as an object to study the change of stress and displacement under thermal cycling load.The simulation results coincided with the test.The results showed that this method can well simulate the wrinkle phenomenon in the process of interface growth.Compared with the finite element method,the deformation of the element was uniform,making it suitable for numerical simulation of the growth process of TGO.The maximum stress of the thermal barrier coating mainly occurred at the convex position,and the lateral change of the convex position had a tendency to aggravate the large deformation of the oxide layer interface.

关 键 词：热障涂层 无网格方法 形函数 热循环 大变形 

分 类 号：V231[航空宇航科学与技术—航空宇航推进理论与工程]