首页/文章/ 详情

code_aster对汽轮机套环的疲劳损伤与裂纹扩展模拟

2年前浏览3514
图片

Code_Aster是法国电力集团(EDF)研发的一款开源有限元仿真软件,适用于固体力学、热学和声学等物理现象,具体细分为静力学、动力学、土-结构相互作用、流固耦合、传热学等方面的仿真分析,具备广阔的应用场景。Code_Aster通过核工业认证,满足法国核安全局和英国核监管部门要求,在各工业领域尤其是能源电力领域有大量的工程和研发应用案例。

01 研究背景


汽轮机是核电站常规岛中最关键的设备之一,负责将二回路水蒸气的热能转化为机械能。而套环是汽轮机的一种钢制旋转件,安装在转子轴的两端,用于固定感应器绕组的线圈头,正常运行工况下,不仅需要承担部分汽轮机的自重载荷,还需承受汽轮机叶片转动过程产生的高温和机械振动(1500转每分钟)。

图片

基于这种情况,细小的缺陷也很可能导致汽轮机发生故障,最终造成停堆。由于汽轮机运行工况的周期性特点,汽轮机套环易发生疲劳损伤现象,从而导致套环产生细微裂纹,或是导致已有的裂纹扩展。

法国电力公司为了评估套环的安全性和生命周期,使用Code_Aster的疲劳计算模块进行了静力加载的裂纹扩展分析和循环加载的疲劳裂纹扩展分析。

02 理论基础——断裂和疲劳扩展的原理


断裂的基本类型有三种:
  • 张开性裂纹(Ⅰ型)
  • 滑开型裂纹(Ⅱ型)
  • 撕开型裂纹(Ⅲ型)
材料和结构的断裂与否常由应力强度因子(Facteurs d‘intensité de contraintes)与材料的标准临界值Kc对比来判断,而Ⅰ型裂纹在实际工程应用是最容易断裂的裂纹类型。
图片
断裂的三种类型
应力强度因子计算公式:
图片
能量释放率公式
图片
 对于线弹性材料:
如果K < KIC,材料不发生断裂;
如果K > KIC,材料发生断裂,裂纹扩展。
疲劳问题的载荷通常以周期性的形式存在,此时的应力强度因子也随加载规律而呈现周期性变化,材料的裂纹扩展则需要用另一种方法来计算:
疲劳裂纹扩展速率  图片
即用裂纹长度的增量和交变应力的循环次数增量的比值来定量计算,表示交变应力每循环一次裂纹长度的平均增量。
根据Paris提出的公式图片
来推得每次循环的应力强度因子的峰值图片
和材料的疲劳断裂临界值图片对比
判断材料的裂纹是否会发生扩展(仅适用于短裂纹)。
图片
周期性加载下的应力强度因子变化曲线示意图
图片
Paris公式

03 模型建立


根据几何和边界条件的对称性,将汽轮机套环简化为1/8模型,通过Salome中的HOMARD模块将网格作密化处理,最终采用了1501970单元数量,包含286338节点数的TETRA4四面体网格。
图片
汽轮机套环的几何模型
图片
汽轮机套环的网格
汽轮机套环的裂纹萌生和扩展风险主要来源有两点:
  • 启堆和停堆期间的汽轮机轴速(离心力)变化,这将导致汽轮机套环的鼓胀和轴向裂纹发生;
  • 交变电流产生的热量引起线圈受热膨胀,由内而外地对套环进行挤压,这将导致汽轮机套环容易发生周向裂纹。
对比两种导致裂纹产生的原因,热膨胀引起的周向裂纹概率远大于轴速变化引起的轴向裂纹。
因此在本研究中仅考虑了周向裂纹,选取其常见的表征类型Ⅰ型(张开型)裂纹。
计算模型主要考虑:
  • 复杂的热力学边界条件;
  • 接触和摩擦问题(套环和线圈等);
  • 不同尺寸的椭圆形裂纹缺陷,使用扩展有限元方法(X-FEM,该方法无需画出裂纹形状的网格)。
图片
汽轮机套环的边界条件
图片
椭圆形裂纹缺陷

04 计算结果分析


图片
汽轮机套环的高斯应力云图 
code_aster的CALC_G模块提供了计算应力强度因子的功能,能够直接输出观测点上的应力强度因子。
在本研究中,我们对套环最容易产生裂纹的部位(槽)建立了X-FEM单元的椭圆张开型裂纹。下图为不同尺寸的裂纹和不同方向上的应力强度因子曲线:
图片
各方向上裂纹尖端的应力强度因子
从图中可观察到,左右两侧的裂纹存在应力集中现象,应力强度因子最大,即裂纹扩展速率最快;而随着初始裂纹尺寸的增大,应力强度因子的最大值也在增大。
针对交变电流产生的周期性热载荷所产生的疲劳问题,裂纹扩展分析主要分为以下步骤:
  • 应力强度因子的计算与周期性的热量输出峰值相关:,建立应力强度因子与热量的映射关系。

图片
应力强度因子峰峰值– 热量变化映射关系图
  • 根据已有的应力强度因子峰峰值计算材料寿命Nf

图片
  • 迭代计算随循环次数而变化裂纹尖端位置,得到下图(初始椭圆几何尺寸a=40mm,b=10mm)。
图片
不同尺寸的椭圆裂纹的裂纹尖端位置
图中黑色实线为初始的裂纹位置,红线和黑色圆形线为实验观测得到的结果,绿线和蓝线为计算得到的结果。
结果显示裂纹最终在套环边界处停止,和真实情况的裂纹相当吻合。在应力集中的左右两侧,真实情况的裂纹位置在计算的裂纹外侧,表明此处的真实裂纹扩展速率比计算的更大,计算结果相对保守。

05 结论


该案例主要对汽轮机套环进行了断裂-疲劳裂纹扩展分析,考虑了由交变电流引起的周期性热载荷,利用测得的交变电流的频率,可以推断裂纹箍在循环载荷下的疲劳寿命,并从而研究最佳的汽轮机运行控制策略,为汽轮机套环部件的寿命评估提供了重要分析依据。


更多信息关注【远算云学院】获取

仿真体系叶轮机械断裂疲劳code_aster
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2022-05-13
最近编辑:2年前
格物CAE
让创新发生
获赞 136粉丝 20文章 51课程 3
点赞
收藏
作者推荐

免费 5.0
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈