一种热-离心载荷作用下多层薄膜耦合应力仿真计算方法技术

本发明专利技术提出一种热

【技术实现步骤摘要】
一种热
‑
离心载荷作用下多层薄膜耦合应力仿真计算方法


[0001]本专利技术涉及一种应用于涡轮叶片的多层薄膜热偶传感器的耦合应力的仿真计算方法，具体适用于多层薄膜结构在涡轮叶片工况下热载荷与离心载荷作用下的耦合应力仿真计算。

技术介绍

[0002]发动机涡轮叶片表面温度将对发动机的性能和寿命产生重要影响，准确测量其工作状态下的表面的温度至关重要。薄膜热电偶直接沉积在涡轮叶片表面并进行实时温度测量，其厚度在微纳米量级。典型的薄膜热偶传感器由多层膜构成，按功能从内到外划分为过渡层、绝缘层、功能层和保护层四部分。在涡轮叶片上适用的薄膜的制备方法为多源离子束共溅，可在大尺寸、大角度曲面合金基底工件制备连续、均匀、致密、粘附力强的薄膜。
[0003]涡轮发动机运行时的工作温度可达到1300℃，由于薄膜传感器各层材料热膨胀系数失配，外界温度载荷将在多层薄膜结构中引入较高的热应力。除了热失配产生的热应力外，基底变形还将引入薄膜内的弹性应力。发动机叶片工作时，离心载荷的作用将对叶片产生径向的拉应力作用，进而产生径向拉应变。薄膜本身质量极小，其内部产生的离心应力可忽略不计，但薄膜传感器直接镀于叶片上，叶片表面的应变将直接传到薄膜之上，在薄膜内部产生应力作用。由于各层薄膜材料弹性模量不同，基底表面应变也将在膜层结构中产生应力失配作用。
[0004]薄膜内应力引起结构发生翘曲、剥离、龟裂等失效模式，开裂后敏感材料有可能发生氧化、断裂，进而使薄膜传感器失效。

技术实现思路

[0005]本专利技术的目的是提出了一种适用于涡轮叶片多层薄膜热偶传感器在热载荷与离心载荷作用下，由于各层材料的热膨胀系数与弹性模量失配而产生的内应力的理论计算与仿真计算方法。具体涉及到通过仿真获得的涡轮叶片表面温度场与离心应变场，理论计算镀于叶片表面的多层薄膜传感器各层薄膜的耦合应力大小。并提出了叶片基底
‑
薄膜传感器模型的耦合应力仿真方法。
[0006]本专利技术由以下步骤实现：
[0007]步骤一：建立涡轮叶片基底模型：
[0008]S11、获得涡轮叶片材料参数：通过文献调研查询材料性能参数表或者采用脉冲激振法、声频共振法、静态法等方法进行测量得到涡轮叶片的弹性模量E
f
，泊松比v
f
，热膨胀系数α
f
，密度G
f
；
[0009]S12、利用ANSYS软件创建涡轮叶片结构模型，建立的模型如图1a和图1b所示。完成几何建模并定义材料参数，具体需定义的参数类型为S11中获取的弹性模量E
f
，泊松比v
f
，热膨胀系数α
f
，密度G
f
，对模型进行网格划分，采用“Element Size”划分类型，Element Size参数设定为5e
‑
3，划分后的结果如图2。
[0010]S13、施加位移约束条件，将榫头施加轴向与周向位移约束；轮盘圆柱面内侧施加圆柱约束，如图3。
[0011]步骤二：仿真得到热载荷下涡轮叶片表面温度分布：
[0012]选取测量点，测量得到涡轮叶片工况下不同位置的表面温度值；对叶片模型对应的测量点位置施加温度值，通过ANSYS中的求解器求解获得模型表面温度场T，如图4a和图4b所示。
[0013]步骤三：仿真得到离心载荷下涡轮叶片表面应变分布：
[0014]通过查询涡轮叶片发动机工况，得到涡轮叶片工作时的转速(6000rad/s)；对轮盘圆柱面施加转速，利用ANSYS workbench的求解器仿真求解涡轮叶片表面离心应变场ε，如图5a和图5b所示。
[0015]步骤四：理论计算各层薄膜内应力：
[0016]S41、获得基底及薄膜各层材料参数：通过S11中的方式测量或查表，获得各层薄膜材料的弹性模量，泊松比，热膨胀系数，密度，具体如下表；
[0017][0018]S42、根据步骤二获得的表面温度场T与步骤三获得的表面离心应变场ε，代入下式获得各层薄膜内应力σ
x,i
。
[0019][0020]式中，E
i
为第i层材料的弹性模量，α
s
为基底的热膨胀系数，α
i
为第i层材料的热膨胀系数，v
i
为第i层材料的泊松比。T
x
为位置x处叶片表面温度，T0为叶片即传感器初始温度，一般取室温，ε
x
为位置x处叶片表面应变。
[0021]步骤五：构建叶片基底
‑
多层薄膜结构模型仿真模型：
[0022]S51、利用ANSYS软件创建叶片基底
‑
多层薄膜结构模型，涡轮叶片主要结构由叶身、缘板与榫头组成，内部包含冷却孔与冷却腔结构。薄膜传感器镀于叶身之上，由叶根延伸到叶尖。对涡轮叶片模型进行简化建模，将缘板与榫头简化为底部为圆柱面的块体，记为叶片底部；将叶身简化为长方体板。
[0023]薄膜传感器为基底
‑
过渡层
‑
绝缘层
‑
功能层
‑
保护层的结构，底部过渡层为1um的NiCrAlY和2um的Al2O3薄膜，绝缘层为2um的Ta2O5和2um的SiO2薄膜，功能层采用0.5um的NiCr/NiSi薄膜，保护层采用4um的SiO2薄膜。
[0024]在简化叶片模型中，涡轮叶片叶身简化为尺寸为10cm*2cm的长方体结构，厚度为1cm，薄膜传感器尺寸为8cm*1cm；自基底从下往上各层薄膜的厚度和参数具体见S41中表，
完成几何建模并定义表中材料参数。对模型进行网格划分，并对薄膜部分网格进行细化，对6层薄膜部分采用“Element Size”划分类型，Element Size参数设定为5e
‑
4，建模及网格划分结果如图6所示；设置基底与薄膜之间、基底与薄膜之间为绑定接触。
[0025]薄膜传感器结构示意图如图7。
[0026]S52、施加位移约束条件，对简化叶片底面施加全约束(即完全固定)；对叶片底部侧面施加轴向与周向位移约束；对底部圆柱面内侧施加圆柱约束。
[0027]步骤六：仿真计算多层薄膜结构内应力大小：
[0028]S61、根据叶片表面温度场，对薄膜
‑
基底模型关键点施加温度值，通过ANSYS workbench的求解器求解获得模型温度场，如图8。
[0029]S62、对底部圆柱面施加转速为6000rad/s，利用ANSYS workbench的求解器仿真得到温度载荷与离心载荷耦合作用下的各层薄膜内应力σ
x,i
，如图9。
[0030]本专利技术所述的应用于涡轮叶片的多层薄膜热偶传感器的热
‑
离心耦合应力的理论计算方法，首先采用ANSYS软件，建立涡轮叶片结构的有限元模型，通过仿真分析得到其表面的温度场本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种热
‑
离心载荷作用下多层薄膜耦合应力仿真计算方法，其特征在于，具体步骤如下：步骤一：建立涡轮叶片基底模型：S11、获得涡轮叶片材料参数：S12、利用ANSYS软件创建涡轮叶片结构模型，S13、施加位移约束条件，将榫头施加轴向与周向位移约束；轮盘圆柱面内侧施加圆柱约束；步骤二：仿真得到热载荷下涡轮叶片表面温度分布：选取测量点，测量得到涡轮叶片工况下不同位置的表面温度值；对叶片模型对应的测量点位置施加温度值，通过ANSYS中的求解器求解获得模型表面温度场T；步骤三：仿真得到离心载荷下涡轮叶片表面应变分布：通过查询涡轮叶片发动机工况，得到涡轮叶片工作时的转速；对轮盘圆柱面施加转速，利用ANSYS workbench的求解器仿真求解涡轮叶片表面离心应变场ε；步骤四：理论计算各层薄膜内应力：S41、获得基底及薄膜各层材料参数：S42、根据步骤二获得的表面温度场T与步骤三获得的表面离心应变场ε，代入下式获得各层薄膜内应力σ
x,i
；式中，E
i
为第i层材料的弹性模量，α
s
为基底的热膨胀系数，α
i
为第i层材料的热膨胀系数，v
i
为第i层材料的泊松比；T
x
为位置x处叶片表面温度，T0为叶片即传感器初始温度，一般取室温，ε
x
为位置x处叶片表面应变；步骤五：构建叶片基底
‑
多层薄膜结构模型仿真模型：S51、利用ANSYS软件创建叶片基底
‑
多层薄膜结构模型；S52、施加位移约束条件，对简化叶片底面施加全约束；对叶片底部侧面施加轴向与周向位移约束；对底部圆柱面内侧施加圆柱约束；步骤六：仿真计算多层薄膜结构内应力大小：S61、根据叶片表面温度场，对薄膜
‑
基底模型关键点施加温度值，通过ANSYS workbench的求解器求解获得模型温度场；S62、对底部圆柱面施加转速，利用ANSYS workbench的求解器仿真得到温度载荷与离心载荷耦合作用下的各层薄膜内应力σ
x,i
。2.根据权利要求1所述的一种热
‑
离心载荷作用下多层薄膜耦合应力仿真计算方法，其特征在于：在步骤S11中，通过文献调研查询材料性能参数表或者采用脉冲激振法、声频共振法、静态法进行测量得到涡轮叶片的弹性模量E
f
，泊松比v
f
，热膨胀系数α
f
，密度G
f
。3.根据权利要求1或2所述的一种热
‑
离心载荷作用下多层薄膜...

【专利技术属性】
技术研发人员：孙宇锋，王佳飞，胡听春，薛雨晴，王一丹，赵广燕，胡薇薇，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：