91学术服务平台

您好,欢迎来到91学术官网!站长邮箱:91xszz@sina.com

发布论文

论文咨询

小破口失水背景下屏蔽泵轴系动力学研究

  2020-08-28    254  上传者:管理员

摘要:以屏蔽式核主泵的主要设计参数为依据,在对屏蔽式核主泵合理简化的基础上进行水力设计与三维模型建立。采用计算流体力学软件(CFX)对正常工况下、进口段小破口失水工况下和出口段小破口失水工况下的屏蔽式核主泵流场进行数值分析,计算叶轮受力与变形。将叶轮受力作为屏蔽式核主泵谐响应分析的输入载荷,分析屏蔽式核主泵轴系的振动情况。计算结果表明:由于进口段破口导致冷却剂流失,叶轮受力、核主泵轴系振幅大幅减小;出口段破口导致核主泵流量短时间内增加,致使叶轮受力略有增加,核主泵轴系振幅最大值仅有小幅增大。

  • 关键词:
  • 叶轮受力
  • 失水事故
  • 屏蔽式核主泵
  • 核工程
  • 谐响应
  • 加入收藏

核主泵是核反应堆冷却剂循环泵的简称,主要用于推动冷却剂在反应堆冷却剂系统中循环运动,将核反应堆产生的热量带到蒸汽发生器中产生蒸汽推动汽轮机旋转发电[1,2,3,4]。失水事故(LOCA)是由回路压力边界出现破口导致的。失水事故会导致冷却剂泄漏、堆芯无法及时冷却、放射性物质进入核电站外部世界[5]。因此,研究核主泵在小破口失水工况下的工作状态有利于采取相应安全措施防止发生核事故。

国内对失水工况下核主泵工作状态已有一些研究,于健[6]研究了核主泵在小破口失水事故中的振动工况、惰转工况及气液两相工况;王学吉[7]通过数值模拟和试验验证研究了CAP1400核主泵在发生下破口失水事故下核主泵的性能和含气率的关系;王海彬[8]通过试验研究和数值模拟的方法对CAP1400核主泵进口段发生失水事故时的泵内流场的压力脉动特性进行了研究。对于小破口工况下核主泵轴系响应尚未有相关报道。

本文以屏蔽式核主泵为研究对象,对其进行有限元建模、流固耦合数值分析,计算其在小破口失水工况下的轴系谐响应,为小破口失水工况下核主泵轴系动力学响应特征提供数据基础。


1、数值计算


图1屏蔽式核主泵的结构图

1.1 模型简化

以屏蔽式核主泵为研究对象,其设计参数如表1[9]所示。屏蔽式核主泵的结构如图1所示,其主要部件包括泵壳,叶轮,上、下飞轮,上、下径向轴承,双向推力轴承,定子外壳,定子铁芯,转子铁芯。

表1屏蔽式核主泵的主要设计参数

屏蔽式核主泵的结构比较复杂,对其进行适当简化再建模。简化内容包括去除与计算无关的部件及尺寸小、影响小的结构,如冷却循环系统、屏蔽电动机及屏蔽套,将上、下飞轮等效为圆盘等[10]。且叶轮、下飞轮、上飞轮有各自对应的质量不平衡量,如表2所示。屏蔽式核主泵的简化三维模型图如图2所示。

表2屏蔽式核主泵的质量不平衡量

图2屏蔽式核主泵的三维模型

1.2 边界条件

在所有需要进行模拟计算区域的外表面上,均需设置边界条件。流体介质的材料选择液态水。入口处采用质量流量边界条件;出口处采用压力边界条件。叶轮处采用移动壁面,边界条件为无滑移壁面;泵壳壁面采用静止壁面,边界条件为无滑移壁面。边界条件参数如表3所示。

1.3 数值分析方法

表3边界条件参数

选用Workbench中的CFX模块进行流体数值分析。将叶轮流体计算域设置为旋转域,旋转轴为Z轴,设置旋转速度为1800r/min;其他流体计算域设置为静止域。湍流模型选为k-Epsilon,k-Epsilon是湍流模式理论中的一种,简称k-ε模型。因为标准k-ε模型可用于复杂几何外部流动问题,又可用于可压缩和不可压缩流体,且精度合理,所以选用标准k-ε湍流模型。数值分析的流程图如图3所示。

图3数值分析流程图


2、计算结果及分析


采用CFD方法计算出失水事故工况下核反应堆冷却剂对核主泵叶轮的作用力,将这个作用力作为核主泵谐响应分析的输入载荷,研究失水事故工况对核主泵轴系振动的影响。按破口所在位置可将破口失水事故分为进口段破口失水事故和出口段破口失水事故;按破口面积可将失水事故分成小破口、中破口及大破口失水事故[11]。破口面积在0.046m2以下的破口被称为小破口。进口段破口位置如图4(a)所示,出口段破口位置如图4(b)所示,本文研究方法适用所有破口面积,选择临界面积0.0457m2进行研究。

图4破口所在位置

2.1 正常工况下核主泵的谐响应分析

图5为正常工况下流体压力分布图,压力最大值为10.2648MPa。图6(a)为正常工况下叶轮的变形图。从图中可以看出,叶轮的变形是沿着径向由内向外增大,最大变形量为0.081785mm。图6(b)为正常工况下叶轮受力图,叶片外缘受力最大,最大受力值为57.438MPa。

图5正常工况下流体压力分布图

从图7可以看出,在7Hz时核主泵轴系的振幅达到最大值2.7253mm,随后振幅一直减小,在30Hz时振幅为0.0138mm。

2.2 进口段破口工况下核主泵谐响应分析

图8为进口段破口工况下流体压力分布图,压力最大值为2.43515MPa,比正常工况下压力减小了76.277%。图9(a)为进口段破口工况下叶轮变形图。从图中可以看出,叶轮外缘变形最大,为0.013671mm,比正常工况下叶轮变形减小了83.284%。图9(b)为进口段破口工况下叶轮受力,受力最大值为10.051MPa,比正常工况下叶轮最大受力减小了82.501%。

图6正常工况下叶轮变形与受力

图7正常工况下谐响应频谱图

图8进口段小破口工况下流体压力分布图

从图10可以看出,在7Hz时核主泵轴系的振幅达到最大值1.891mm,比正常工况下核主泵轴系最大振幅下降了30.613%,随后振幅一直减小,在30Hz时振幅为0.00987mm,比正常工况下振幅下降了28.478%。

图9进口段小破口工况下叶轮变形与受力

图10进口段小破口工况下谐响应频谱图

2.3 出口段小破口工况下核主泵谐响应分析

图11为出口段破口工况下流体压力分布图,压力最大值为10.4741MPa,比正常工况下压力增大了2.039%。图12(a)为出口段破口工况下叶轮变形图,叶轮外缘变形最大,变形值为0.08311mm,比正常工况下叶轮变形增大1.62%。图12(b)为出口段破口工况下叶轮受力图,叶片受力最大,最大受力为60.977MPa,比正常工况下叶轮受力增加了6.161%。从图13可以看出,在7Hz时核主泵轴系的振幅达到最大值(2.7504mm),比正常工况下核主泵振幅增加了0.921%,随后振幅一直减小,在30Hz时振幅为0.0138mm,与正常工况下核主泵振幅一致。图14为3种工况下核主泵轴系谐响应图,从图中可以看出进口段破口工况下核主泵轴系振幅最小,出口段破口工况下核主泵振幅略大于正常工况下核主泵振幅。

图11出口段破口工况下流体压力分布图

2.4 计算结果分析

1)进口段破口工况下冷却剂从破口流失无法到达叶轮处,导致叶轮最大受力比正常工况减小了82.501%,最大变形减小了83.284%,最大振幅下降了28.478%。

图12出口段破口工况下叶轮变形与受力

2)出口段出现破口导致出口面积增大,冷却剂流速在短时间内增大,所以叶轮受力增大。与正常工况相比叶轮最大受力增加了6.161%,最大变形增大1.62%,轴系最大振幅增加了0.921%。

图13出口段破口工况下谐响应频谱图

图143种工况下谐响应频谱图

3)通过监测7Hz时核主泵轴系振幅变化可判断核主泵工况。当核主泵轴系在7Hz下的振幅大幅下降时进口段出现小破口;当核主泵轴系在7Hz下的振幅小幅增加时出口段出现破口。


3、结论


通过对核主泵在正常工况下和在小破口工况下的谐响应分析后,得到了以下结论:1)与正常工况下核主泵工作状态相比,进口段小破口工况下叶轮受力、变形及振幅大幅减小;2)与正常工况下核主泵工作状态相比,出口段小破口工况下,叶轮受力、变形略有增大,核主泵轴系振幅仅有小幅变化。


参考文献:

[1]朱荣生,刘永,王秀礼,等.失水事故初期核主泵气液两相流动特性[J].哈尔滨工程大学学报,2017,38(9):1399-1405.

[2]付强,习毅,朱荣生,等.含气率对AP1000核主泵影响的非定常分析[J].振动与冲击,2015,34(6):132-136.

[3]朱荣生,郑宝义,袁寿其,等.1000MW核主泵失水事故工况下气液两相流分析[J].原子能科学技术,2012,46(10):1202-1206.

[4]王俊,张永超,王达,等.混流式核主泵内部复杂流动结构分析[J].流体机械,2017,45(10):58-63.

[5]博金海,王飞.小破口失水事故研究综述[J].核科学与工程,1998,18(2):172-179.

[6]于健.小破口事故工况下核主泵性能研究[D].大连:大连理工大学,2012.

[7]王学吉.小破口失水工况下CAP1400核主泵气液两相流模拟[D].镇江:江苏大学,2017.

[8]王海滨.进口破口事故下CAP1400核主泵内压力脉动特性研究[D].镇江:江苏大学,2018.

[9]孙汉虹,程平东,缪鸿兴,等.第三代核电技术AP1000[M].北京:中国电力出版社,2010.

[10]孙丹.大型屏蔽电动机泵轴向力分析[D].沈阳:沈阳工业大学,2018.

[11]丁书华,钱立波,吴丹.反应堆冷却机泵水力特性对大破口失水事故的影响研究[J].核动力工程,2013,34(1):192-195.


王硕,赵晶,孙丹,夏海明.小破口失水工况下屏蔽泵轴系动力学分析[J].机械工程师,2020(08):57-59+62.

基金:辽宁省重点研发项目(2017106016).

分享:

91学术论文范文

相关论文

推荐期刊

网友评论

加载更多

我要评论

核安全

期刊名称:核安全

期刊人气:2730

期刊详情

主管单位:生态环境部

主办单位:国家环境保护总局核安全中心

出版地方:北京

专业分类:科技

国际刊号:1672-5360

国内刊号:11-5145/TL

邮发代号:无

创刊时间:2003年

发行周期:双月刊

期刊开本:大16开

见刊时间:4-6个月

论文导航

查看更多

相关期刊

热门论文

【91学术】(www.91xueshu.com)属于综合性学术交流平台,信息来自源互联网共享,如有版权协议请告知删除,ICP备案:冀ICP备19018493号

微信咨询

返回顶部

发布论文

上传文件

发布论文

上传文件

发布论文

您的论文已提交,我们会尽快联系您,请耐心等待!

知 道 了

登录

点击换一张
点击换一张
已经有账号?立即登录
已经有账号?立即登录

找回密码

找回密码

你的密码已发送到您的邮箱,请查看!

确 定