某型发动机涡轮叶片的蠕变寿命分析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
89. 77 4 3 【 44 3 l 03 .
() 3 在实 际叶片模 型上施加 的约束 。 如下 3 : 有 项
一
是在各个榫头上齿接触面施加沿径向约束 ; 二是 由于涡轮叶片工作状态下相互抵紧 ,在座
台与叶冠施加周 向约束 ; 三是 为 了防 止 叶片 在离 心 载 荷 作 用下 的翘 曲变
动机 ,涡轮转子及叶片作为燃气轮机关键零部件 , 在 工作中承受着很大的热负荷、 离心力 、 气流力。 由于其 工作 条 件恶 劣 , 得 涡轮 盘 叶片成 为故 障多 发 的零部 使 有限元分析过程流程如图 1 所示。
件, 其寿命直接决定 了燃气轮机的使 用寿命 , 涡轮叶 片在高温燃气和循环载荷下的工作可靠性 , 将直接影
22 材 料 属性 .
涡轮 叶片所使用 的金属材料 为在 R n ’0 ee8 基础 上 , 当提高合金元素铝和钛 的含量 , 适 并添加合金元 素铌发展而成的以 N —c 固溶体为基体 、有较高的 i r 热强性 、 良好 的综合性 能、 足够 的抗氧化性能和满意 的组 织稳 定 性 , 于耐 高温 的强 韧化材 料 。 属 在A SS N Y 中输人材料参量。密度 D n , es泊松 比 P X 弹性模量 E 导热系数 K x 如 图 3 R Y, X, x ( 所示 ) , 热膨胀 系数 A P ( 图 4 LX如 所示 )比热容进行有 限元 ,
24 离 心载 荷和 温度场 综 合分析 .
由图 7叶片 变形 情 况 、 图 8叶 身 等 效应 力 云 图
以及 图 9叶片等效应ຫໍສະໝຸດ Baidu云图综合分析 ,叶身部分应 力主要集中于叶身与叶冠接触部分 、叶身与榫头接
温度( ℃)
触部分 ,即应力危险 区域主要分 布于涡轮 叶片上 3
个 部位 :
在图 1 l温 度 一屈 服 / 限 应 力 曲 线 中 , 极 将 到 8 OI是屈 服 极 限为 6 87MP ; 90~90℃之 25 不 同工 况下 的应 力分析 比较 7c = 1. a 将 0 5 . 间 的屈服 极 限视 为 线性 变 化 , 则可 得 到 9 0C 2  ̄是屈 服
温 度( ℃)
图 1 1
温度 一屈 服 / 限应 力曲线 极
叶身部分最大应力出现于叶尖截 面 ,应力值为
4 67 a如 图 1 所 示 。 5 . MP , 9 2
() 1 叶身与榫头接触部位。 此处为叶身部分受到 离 心 载 荷最 大 的部 位 ,由于 约束 作 用会 产 生应 力 集 中还有弯曲应力 、 扭转应力存在 。而且叶身部分与榫 头部分温度分布相差较大 , 进而引起很大的热应力 。 () 2 叶身与叶冠接触部位。此处叶身较薄, 与叶
摘 要 : 对某型航 空发 动机 的涡轮 叶片 , 针 开展 了应力计算 , 进行 给定转速和 温度场条件 下的弹 塑性有限元分析 , 到应 得
力场; 对叶片上的危险位置和最大应力值进行分析 , 计算提高工况时的涡轮叶片的应力场, 分析温度场提高和转速提高
对叶身应力分布的影响和最 大应 力值 的变化程度 ; 对温度场 和离心栽荷分析计算的结果 , 针 对于涡轮 叶片结构设计 、 故
() 3 叶身温度最高部位。 此处产生离心应力水平 中等 , 由于 高 温会 使 得材 料 性 能 的改 变 , 易 产生 但 容 较大的危险。 由图 l O叶片温度分布情况可 以看出 , 叶身与榫 头 接 触 部 位 以 及 叶 身 与 叶 冠 接 触 部 位 的 温 度 在 80 8 0℃之间 , 3~ 7 为提高安全性 , 8OC 以 7 q标准进行
,
00, . 0 共计划分单元数 4 2 个 , 8 537 节点数 8 8 7 4 个。 5
将 H prm s ye eh的 网格 模 型 导 入 A ss软 件 中 ny
23 载 荷施 加 .
() 1 离心 载荷 。在 叶片模 型上施 加 绕 轴 的
时, 进行涡轮叶片的应力分析。由于使得导人的模型
响飞行安全l以及发动机的使用寿命[] 生 2。 - 5 因此 , 对于涡轮叶片蠕变寿命 的研究 , 将对提高 涡轮 叶片设计 以及航空发 动机整体可靠 性提升 , 都 有十分显著的作用 。
图 1 有 限 元分 析流 程
1 研 究 方 法
2 涡轮叶片的有 限元分析
通 过采用 比较研 究法 以及数 值模 拟 的研 究方 2 1 有 限元模 型 . 法, 对涡轮叶片进行模拟计算 , 利用 已选择 的蠕变寿 首先利用 U G软件 建立某型航空发动机涡轮叶 命预测方法 , 对所 给的叶片模型进行计算模拟 ; 首先 片部件的几何模型。根据分析问题的需要 , 叶片模 对 对 所 给 的模 型进 行合 理 的简 化 ,根 据 叶片 的结构 特 型采用 一定 的简化措 施 。研究 过 程 中 , 点考 察 叶身 重 点 和 载荷 分 布 划 分 网格 , 导人 A S S进 行 模拟 计 部 分 负 载情 况 ,对 榫 头 和轮 盘 榫槽 间 的摩 擦 力 不 做 再 NY 算 。对 叶 片进行 热 分 析 以及 应 力分 析 。 关 注 , 榫槽 部 位 做简 化 处理 , 略一 些 非关 键 部 位 对 忽 () 1 热分析 。选取最高温度作为主要计算状态。 小 的销 钉孔 以及 小 的倒 角 。 由涡轮叶片的热源和结构 ,在有限元模型 中构成热 将 叶片模型通过 Pr o d a sl 格式从 U a i G导入 H — y 边界 条 件进 行模 拟 计算 。 pr e 软件 中, em s h 使用 四面体网格和六面体 网格 的混 ( ) 力 分 析 。 由资料 中可 以得 到 叶片 的载 荷 , 合 网格划 分 , 2应 在保 证精 度 的 同时大 大减 少 了工作 量 。 并结合飞行状态计算 ,可知各个状态下的温度场和 混合 网格模式下 , 利用三维 2 节点六面体单元 0 应力载荷 , 进而进行有限元分析。计划选取温度最高 与 1 0节点 四面体单元综合分网 , 两种单元过渡处为 部位 、叶根叶背叶尖 、叶身应力最大点进行考核计 金字塔单元 ;单元大小为 1 0 ,最小单元大小 为 .0 0
图 4 R n ’0材 料 的温 度 一热 膨 胀 系数 曲线 e e8 6
《 装备制造技术)02 2 1 年第 6 期
图 7 叶片变形情况
50 0 60 0 70 0 80 0 9O O l0 0 0 I 10 0
; I 一 一
图8 叶身等效应力云图 图9 叶片等效应力云图
极 限为 4 5 1 a 9 .4MP 。
4 67 5 .9
应力危险区域最大应力
4 I1 0 .l 4 31 3 .9
温度最高处 榫头连 接处 叶冠连接处
安全系数
I 3 .4 2
I 2 .7 4
1 5 .3 3
因此实 际叶片工作状态下叶身部分应力水平处 80~90o之 间的屈 服极 限视为 线 性变 化 ,则 可得 于 较为 安全 的范 围 。 5 0 C
障预防 维修 以及航 空发动机整 ̄ - 靠性有重要意 义。 q . -
关键词 : 涡轮 叶片 ; 应力分析 ; 有限元 ; 数值模拟
中图 分 类 号 : 2 1 V 3 文 献 标识 码 : A 文 章 编 号 :6 2 5 5 ( 0 2 0 - 0 5 0 1 7 — 4 X 2 1 )6 0 0 - 4
收 稿 日期 :0 2 0 — 9 2 1 — 3 1 作者简介 : 柯铭 亮( 97 ) 男 , 18 一 , 福建福州人 , 助理工程师 , 本科学历 , 研究方向为机械 制造 与设计 。
5
E u p n Ma u a t n e h o o y N . 2 1 q i me t n f cr g T c n l g o 6 0 2 i
冠部分由于约束作用 , 会产生应力集 中, 同时 由于温
图 1 2
叶尖截面等效应 力云 图
但在叶片最高温部位( 图 1 所示 )高温使得 如 3 , 度分布不同 , 引起热应力。工作载荷会使此处产生较 大塑性变形 , 而使应力重新分布 , 对叶身寿命产生较 材料性能恶化 , 虽然此处应力并非最大值 , 但相对叶 根、 叶尖部分 , 此处安全系数低 , 危险程度大。 大 的 影 响。
《 装备制造技术) 02 2 1 年第 6 期
某型发 动机涡轮叶片的蠕变寿命分析
柯 铭 亮 杨伟锋 , 丁玉亭 薛伦生 ’ , , , /
(. 1 空军工程大学 导弹学院, 陕西 三原 7 30 ;. 180 2 解放军 974 队装备处 , 30 部 北京 通州 114 ) 0 19
形 , 座 台与 叶冠施 加轴 向约 束 ( 图 6 。 在 如 )
3 0 4 0 6 o 7 0 9 O 10 0 0 5 0 5 O 5
温度( ℃)
图 3 R n ’0材 料的温度 一弹性模量 曲线与导热系数 曲线 e e8
图 6 叶 片 约 束 施 加 情 况
,
L:= := 壤蘸再豫 触 f == = =
\
I
{
5O O
60 0
70 0
8t 0}
90 0
l0 0 0
I10 0
分析的准备。
表 1 材 料 参 量
温度 , ℃
图 5 温 度 场 加 载 对 比情 况
密 (/ ) 比 容C l 泊 度Pkm 1 热 松比 g3
Ⅲ 一
位置
叶根截 面
轴心距( m)
O 4 . O3 2
温度( ) ℃
85 7
叶高 1 处 / 4 叶高 1 处 / 2 叶高 3 处 / 4 叶尖截 面
05 . 69 2 0 7 . 34 2 0 9O .O 2 0366 . O
98 l 9O 2 9o 0 88 3
u( 、 y 航空发动机涡轮叶片构件 , 长期在高温和循环载 算 。主要使用的软件包括 : G 建立实体模型 )H — em s  ̄模型进行 网格划分 ) h 以及 A s S 进行仿 NY ( 荷共同作用下工作 ,疲劳 一蠕变交互作 用导致的断 pr e ( 。 裂, 是此 类 热锻构 件 失效 的主要 原 因【 而对 于航 空 发 真模拟计算 ) l 】 。
1 2 . d 旋转角速度, 3 1a/ 3 6r s 使得叶片产生离心拉应力。 在叶背位置出现明显缺块和网格划分单元品质不高 () 2 温度载荷。选取 叶身部位的关键截面温度 , 等原因 ,选择 s i17 o d8 以及 cm il 单元对实 际模 如表 2 l ob 4 n 。 型 网格划分进行修正 ,共计 2 8 3 6 单元 , l 17 2 s i 8 单 0d 表 2 关键截面温度数据 元数 2 , 4个 ,o b l 单元数 2 个 , 32 4 cm i 4 n 4 整个有限元 模型共计 4 7 2 个节点 , 3 如图 2 所示。
=
蕾
…
将关键截面温度分布与实际工作温度场进行 比
{ l
图2 涡 轮 叶片 有 限元 模 型
对发现 ,利用关键截面的温度数据所施加的温度场 与实际工作温度场 的差别属于允许范 围内 ( 图 5 如 所示 )可 以利用这组数据进行温度载荷的施加。 ,
; 一实际工 作温 度场 ; ~
分 析计 算 。
图 1 叶身最 高温度截面应 力云 图 3
;
图 1 叶 片 温 度分 布 情 况 0
在实际 叶片工作载荷条件下 ,叶身部分较为安 全, 各部分最大应力都没有超过屈服极限 , 安全系数 较为合理 , 如表 3 。
表3
应 力 ( a MP )
~ 一~ .
() 3 在实 际叶片模 型上施加 的约束 。 如下 3 : 有 项
一
是在各个榫头上齿接触面施加沿径向约束 ; 二是 由于涡轮叶片工作状态下相互抵紧 ,在座
台与叶冠施加周 向约束 ; 三是 为 了防 止 叶片 在离 心 载 荷 作 用下 的翘 曲变
动机 ,涡轮转子及叶片作为燃气轮机关键零部件 , 在 工作中承受着很大的热负荷、 离心力 、 气流力。 由于其 工作 条 件恶 劣 , 得 涡轮 盘 叶片成 为故 障多 发 的零部 使 有限元分析过程流程如图 1 所示。
件, 其寿命直接决定 了燃气轮机的使 用寿命 , 涡轮叶 片在高温燃气和循环载荷下的工作可靠性 , 将直接影
22 材 料 属性 .
涡轮 叶片所使用 的金属材料 为在 R n ’0 ee8 基础 上 , 当提高合金元素铝和钛 的含量 , 适 并添加合金元 素铌发展而成的以 N —c 固溶体为基体 、有较高的 i r 热强性 、 良好 的综合性 能、 足够 的抗氧化性能和满意 的组 织稳 定 性 , 于耐 高温 的强 韧化材 料 。 属 在A SS N Y 中输人材料参量。密度 D n , es泊松 比 P X 弹性模量 E 导热系数 K x 如 图 3 R Y, X, x ( 所示 ) , 热膨胀 系数 A P ( 图 4 LX如 所示 )比热容进行有 限元 ,
24 离 心载 荷和 温度场 综 合分析 .
由图 7叶片 变形 情 况 、 图 8叶 身 等 效应 力 云 图
以及 图 9叶片等效应ຫໍສະໝຸດ Baidu云图综合分析 ,叶身部分应 力主要集中于叶身与叶冠接触部分 、叶身与榫头接
温度( ℃)
触部分 ,即应力危险 区域主要分 布于涡轮 叶片上 3
个 部位 :
在图 1 l温 度 一屈 服 / 限 应 力 曲 线 中 , 极 将 到 8 OI是屈 服 极 限为 6 87MP ; 90~90℃之 25 不 同工 况下 的应 力分析 比较 7c = 1. a 将 0 5 . 间 的屈服 极 限视 为 线性 变 化 , 则可 得 到 9 0C 2  ̄是屈 服
温 度( ℃)
图 1 1
温度 一屈 服 / 限应 力曲线 极
叶身部分最大应力出现于叶尖截 面 ,应力值为
4 67 a如 图 1 所 示 。 5 . MP , 9 2
() 1 叶身与榫头接触部位。 此处为叶身部分受到 离 心 载 荷最 大 的部 位 ,由于 约束 作 用会 产 生应 力 集 中还有弯曲应力 、 扭转应力存在 。而且叶身部分与榫 头部分温度分布相差较大 , 进而引起很大的热应力 。 () 2 叶身与叶冠接触部位。此处叶身较薄, 与叶
摘 要 : 对某型航 空发 动机 的涡轮 叶片 , 针 开展 了应力计算 , 进行 给定转速和 温度场条件 下的弹 塑性有限元分析 , 到应 得
力场; 对叶片上的危险位置和最大应力值进行分析 , 计算提高工况时的涡轮叶片的应力场, 分析温度场提高和转速提高
对叶身应力分布的影响和最 大应 力值 的变化程度 ; 对温度场 和离心栽荷分析计算的结果 , 针 对于涡轮 叶片结构设计 、 故
() 3 叶身温度最高部位。 此处产生离心应力水平 中等 , 由于 高 温会 使 得材 料 性 能 的改 变 , 易 产生 但 容 较大的危险。 由图 l O叶片温度分布情况可 以看出 , 叶身与榫 头 接 触 部 位 以 及 叶 身 与 叶 冠 接 触 部 位 的 温 度 在 80 8 0℃之间 , 3~ 7 为提高安全性 , 8OC 以 7 q标准进行
,
00, . 0 共计划分单元数 4 2 个 , 8 537 节点数 8 8 7 4 个。 5
将 H prm s ye eh的 网格 模 型 导 入 A ss软 件 中 ny
23 载 荷施 加 .
() 1 离心 载荷 。在 叶片模 型上施 加 绕 轴 的
时, 进行涡轮叶片的应力分析。由于使得导人的模型
响飞行安全l以及发动机的使用寿命[] 生 2。 - 5 因此 , 对于涡轮叶片蠕变寿命 的研究 , 将对提高 涡轮 叶片设计 以及航空发 动机整体可靠 性提升 , 都 有十分显著的作用 。
图 1 有 限 元分 析流 程
1 研 究 方 法
2 涡轮叶片的有 限元分析
通 过采用 比较研 究法 以及数 值模 拟 的研 究方 2 1 有 限元模 型 . 法, 对涡轮叶片进行模拟计算 , 利用 已选择 的蠕变寿 首先利用 U G软件 建立某型航空发动机涡轮叶 命预测方法 , 对所 给的叶片模型进行计算模拟 ; 首先 片部件的几何模型。根据分析问题的需要 , 叶片模 对 对 所 给 的模 型进 行合 理 的简 化 ,根 据 叶片 的结构 特 型采用 一定 的简化措 施 。研究 过 程 中 , 点考 察 叶身 重 点 和 载荷 分 布 划 分 网格 , 导人 A S S进 行 模拟 计 部 分 负 载情 况 ,对 榫 头 和轮 盘 榫槽 间 的摩 擦 力 不 做 再 NY 算 。对 叶 片进行 热 分 析 以及 应 力分 析 。 关 注 , 榫槽 部 位 做简 化 处理 , 略一 些 非关 键 部 位 对 忽 () 1 热分析 。选取最高温度作为主要计算状态。 小 的销 钉孔 以及 小 的倒 角 。 由涡轮叶片的热源和结构 ,在有限元模型 中构成热 将 叶片模型通过 Pr o d a sl 格式从 U a i G导入 H — y 边界 条 件进 行模 拟 计算 。 pr e 软件 中, em s h 使用 四面体网格和六面体 网格 的混 ( ) 力 分 析 。 由资料 中可 以得 到 叶片 的载 荷 , 合 网格划 分 , 2应 在保 证精 度 的 同时大 大减 少 了工作 量 。 并结合飞行状态计算 ,可知各个状态下的温度场和 混合 网格模式下 , 利用三维 2 节点六面体单元 0 应力载荷 , 进而进行有限元分析。计划选取温度最高 与 1 0节点 四面体单元综合分网 , 两种单元过渡处为 部位 、叶根叶背叶尖 、叶身应力最大点进行考核计 金字塔单元 ;单元大小为 1 0 ,最小单元大小 为 .0 0
图 4 R n ’0材 料 的温 度 一热 膨 胀 系数 曲线 e e8 6
《 装备制造技术)02 2 1 年第 6 期
图 7 叶片变形情况
50 0 60 0 70 0 80 0 9O O l0 0 0 I 10 0
; I 一 一
图8 叶身等效应力云图 图9 叶片等效应力云图
极 限为 4 5 1 a 9 .4MP 。
4 67 5 .9
应力危险区域最大应力
4 I1 0 .l 4 31 3 .9
温度最高处 榫头连 接处 叶冠连接处
安全系数
I 3 .4 2
I 2 .7 4
1 5 .3 3
因此实 际叶片工作状态下叶身部分应力水平处 80~90o之 间的屈 服极 限视为 线 性变 化 ,则 可得 于 较为 安全 的范 围 。 5 0 C
障预防 维修 以及航 空发动机整 ̄ - 靠性有重要意 义。 q . -
关键词 : 涡轮 叶片 ; 应力分析 ; 有限元 ; 数值模拟
中图 分 类 号 : 2 1 V 3 文 献 标识 码 : A 文 章 编 号 :6 2 5 5 ( 0 2 0 - 0 5 0 1 7 — 4 X 2 1 )6 0 0 - 4
收 稿 日期 :0 2 0 — 9 2 1 — 3 1 作者简介 : 柯铭 亮( 97 ) 男 , 18 一 , 福建福州人 , 助理工程师 , 本科学历 , 研究方向为机械 制造 与设计 。
5
E u p n Ma u a t n e h o o y N . 2 1 q i me t n f cr g T c n l g o 6 0 2 i
冠部分由于约束作用 , 会产生应力集 中, 同时 由于温
图 1 2
叶尖截面等效应 力云 图
但在叶片最高温部位( 图 1 所示 )高温使得 如 3 , 度分布不同 , 引起热应力。工作载荷会使此处产生较 大塑性变形 , 而使应力重新分布 , 对叶身寿命产生较 材料性能恶化 , 虽然此处应力并非最大值 , 但相对叶 根、 叶尖部分 , 此处安全系数低 , 危险程度大。 大 的 影 响。
《 装备制造技术) 02 2 1 年第 6 期
某型发 动机涡轮叶片的蠕变寿命分析
柯 铭 亮 杨伟锋 , 丁玉亭 薛伦生 ’ , , , /
(. 1 空军工程大学 导弹学院, 陕西 三原 7 30 ;. 180 2 解放军 974 队装备处 , 30 部 北京 通州 114 ) 0 19
形 , 座 台与 叶冠施 加轴 向约 束 ( 图 6 。 在 如 )
3 0 4 0 6 o 7 0 9 O 10 0 0 5 0 5 O 5
温度( ℃)
图 3 R n ’0材 料的温度 一弹性模量 曲线与导热系数 曲线 e e8
图 6 叶 片 约 束 施 加 情 况
,
L:= := 壤蘸再豫 触 f == = =
\
I
{
5O O
60 0
70 0
8t 0}
90 0
l0 0 0
I10 0
分析的准备。
表 1 材 料 参 量
温度 , ℃
图 5 温 度 场 加 载 对 比情 况
密 (/ ) 比 容C l 泊 度Pkm 1 热 松比 g3
Ⅲ 一
位置
叶根截 面
轴心距( m)
O 4 . O3 2
温度( ) ℃
85 7
叶高 1 处 / 4 叶高 1 处 / 2 叶高 3 处 / 4 叶尖截 面
05 . 69 2 0 7 . 34 2 0 9O .O 2 0366 . O
98 l 9O 2 9o 0 88 3
u( 、 y 航空发动机涡轮叶片构件 , 长期在高温和循环载 算 。主要使用的软件包括 : G 建立实体模型 )H — em s  ̄模型进行 网格划分 ) h 以及 A s S 进行仿 NY ( 荷共同作用下工作 ,疲劳 一蠕变交互作 用导致的断 pr e ( 。 裂, 是此 类 热锻构 件 失效 的主要 原 因【 而对 于航 空 发 真模拟计算 ) l 】 。
1 2 . d 旋转角速度, 3 1a/ 3 6r s 使得叶片产生离心拉应力。 在叶背位置出现明显缺块和网格划分单元品质不高 () 2 温度载荷。选取 叶身部位的关键截面温度 , 等原因 ,选择 s i17 o d8 以及 cm il 单元对实 际模 如表 2 l ob 4 n 。 型 网格划分进行修正 ,共计 2 8 3 6 单元 , l 17 2 s i 8 单 0d 表 2 关键截面温度数据 元数 2 , 4个 ,o b l 单元数 2 个 , 32 4 cm i 4 n 4 整个有限元 模型共计 4 7 2 个节点 , 3 如图 2 所示。
=
蕾
…
将关键截面温度分布与实际工作温度场进行 比
{ l
图2 涡 轮 叶片 有 限元 模 型
对发现 ,利用关键截面的温度数据所施加的温度场 与实际工作温度场 的差别属于允许范 围内 ( 图 5 如 所示 )可 以利用这组数据进行温度载荷的施加。 ,
; 一实际工 作温 度场 ; ~
分 析计 算 。
图 1 叶身最 高温度截面应 力云 图 3
;
图 1 叶 片 温 度分 布 情 况 0
在实际 叶片工作载荷条件下 ,叶身部分较为安 全, 各部分最大应力都没有超过屈服极限 , 安全系数 较为合理 , 如表 3 。
表3
应 力 ( a MP )
~ 一~ .