等离子体喷涂锆酸钆涂层残余应力的仿真计算.pdfVIP

第十四届国际热喷涂研讨会论文 (ITSS’2011) 第十五届全国热喷涂年会论文 (CNTSC’2011) 2. 实验过程 2.1 涂层制备 GZ 的合成原料为市售 Gd O (99.5%)和 ZrO (99.89%)粉末，按化学配比，经高温合成并球磨破 2 3 2 碎、喷雾造粒后，获得适合喷涂的粉体。采用 Sulzer Metco 公司的 PT-A2000 等离子体喷涂系统和 F4-MB 喷枪制备涂层，喷涂工艺参数见文献[6]。热震试验的试样采用尺寸为 20mm ×30mm ×1.5mm 的镍基高温合金为基材，涂层厚度为200~250μm 。热震试验是将涂层试样置于 1100℃恒温的管式炉 中~10min ，迅速取出并投入室温去离子水中，1min 后将试样从水中取出，用压缩空气吹干，这样重 复一次为一个循环。 2.2 有限元模型建立 2.2.1 几何模型和网格划分 本文采用的模型如图 1 所示，为基材（Sub ），结合层（BC ）和陶瓷层（TC ）构成的三层结构， 其厚度分别为 1.5mm，20μm 和 100μm 。采用理想的正弦波形（波长λ，波幅A ）模拟实际制备GZ 涂层 TC/BC 界面的起伏不平，模型中波长 λ取值 30μm 。为减少模拟工作量，采用周期性边界条件， 取其中一个完整正弦波周期作为模拟计算对象。在模型的上端施加 x 方向位移约束条件，用于防止 模型的刚性位移；模型的下端施加周期性边界条件，使得该边上所有节点在 x 方向上位移相等。为 保证模拟精确性，在 TC /BC 界面结合处采用密集单元作为过渡，每个周期波内共有 16900 个四边 形网格单元。 图 1 有限元模型与网格划分 Fig. 1 Geometry of the model used and the detailed view of the finite element mesh 2.2.2 材料属性 表 1 给出了涂层内各层材料随温度变化的物性参数：包括热膨胀系数(α) 、热导率(λ) 、比热(c)、 密度(ρ) 、弹性模量(E)及泊松比(ν) 。考虑到高温下材料的塑性变形对应力再分布的影响，模型中采 用双线性随动强化本构关系描述耐热合金和结合层的非弹性属性。 - 16 - 第十四届国际热喷涂研讨会论文 (ITSS’2011) 第十五届全国热喷涂年会论文 (CNTSC’2011) 表 1 模型中输入的材料物性参数 Table 1 Material properties used in the model Material T E ρ α ν κ C σy (℃) (GPa) (kg/m3) (10-6/K) (W/m K) (J/kg K) (MPa) GZ 25 140 [1] 6019 7.4 0.3 … (a) 415 138 400