CN104636549A - 一种冲击式水轮机水斗动态应力分析方法 - Google Patents

一种冲击式水轮机水斗动态应力分析方法 Download PDF

Info

Publication number
CN104636549A
CN104636549A CN201510056000.8A CN201510056000A CN104636549A CN 104636549 A CN104636549 A CN 104636549A CN 201510056000 A CN201510056000 A CN 201510056000A CN 104636549 A CN104636549 A CN 104636549A
Authority
CN
China
Prior art keywords
bucket
dynamic stress
stress
dynamic
pressure
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510056000.8A
Other languages
English (en)
Other versions
CN104636549B (zh
Inventor
刘晶石
姜铁良
贾伟
田超
王珏
吕桂萍
钟苏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Electric Machinery Co Ltd
Original Assignee
Harbin Electric Machinery Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Electric Machinery Co Ltd filed Critical Harbin Electric Machinery Co Ltd
Priority to CN201510056000.8A priority Critical patent/CN104636549B/zh
Publication of CN104636549A publication Critical patent/CN104636549A/zh
Application granted granted Critical
Publication of CN104636549B publication Critical patent/CN104636549B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及一种冲击式水轮机水斗动态应力分析方法。利用流体计算软件分析水斗周围流场的压力分布,提取水斗过流表面动态压力场数据,编写程序进行数据转换,生成水斗结构分析压力载荷文件并加以修正,循环加载求解,绘制水斗动态应力变化曲线。该方法的优点是:以压力场数据作为结构分析载荷,并考虑了压力场载荷的修正,计算精度较高。利用该方法可以得到冲击式水轮机运行过程中,水斗的应力随时间的变化情况,对掌握水斗工作受力状态以及评估水斗的动应力和静应力水平具有重要意义。

Description

一种冲击式水轮机水斗动态应力分析方法
技术领域:
本发明涉及冲击式水轮机水斗动态应力分析方法。
背景技术:
水斗是冲击式水轮机的核心部件,水斗结构复杂且曲面变化剧烈,在机组运行过程中承受巨大的脉冲载荷,极易发生断斗事故,必须严格控制水斗的动应力和静应力水平。
目前,在冲击式水轮机水斗应力分析中,通常考虑水轮机额定出力工况和离心力工况,只采用简化的射流力载荷进行计算,作用在转轮上的射流力按照一定比例分别施加在相邻的几个水斗上,同时考虑机组转速引起的离心力载荷以及转轮与主轴连接螺栓的固定约束。利用该方法只能得到水斗运行过程中的最高应力和最低应力。
这种分析方法的缺点是:只采用简化的射流力载荷作为水斗应力分析有限元模型载荷,计算精度较差;分析方法简单,不能精确描述水斗在机组运行过程中的应力水平和应力变化情况;随着机组容量和设计水头的日益增加,冲击式水轮机水斗的设计难度越来越大,这种分析方法已经不能满足当前设计需求。
发明内容:
本发明提供了一种冲击式水轮机水斗动态应力分析方法。利用流体计算软件分析水斗周围流场的压力分布,提取水斗过流表面动态压力场数据,编写程序进行数据转换,生成水斗结构分析压力载荷文件并加以修正,循环加载求解,绘制水斗动态应力变化曲线。该方法的优点是:以压力场数据作为结构分析载荷,并考虑了压力场载荷的修正,计算精度较高。利用该方法可以得到冲击式水轮机运行过程中,水斗的应力随时间的变化情况,对掌握水斗工作受力状态以及评估水斗的动应力和静应力水平具有重要意义。本技术方案的技术方案为:
1)运用商业软件建立冲击式水轮机水斗周围流体的数学模型,设置计算域的计算方程和边界条件等计算参数,求解并提取水斗过流表面的动态压力场数据;
2)运用商业软件建立冲击式水轮机水斗结构的有限元模型,在水斗连接螺栓分布圆施加固定约束,对水斗转动中心施加转速载荷;
3)编写数据转换程序,利用流体压力场数据获得水斗动态应力分析压力载荷文件;
4)在水斗过流表面循环加载不同时刻单个水斗表面压力载荷文件并求解,提取结果文件,绘制单个水斗力矩曲线;
5)对水斗压力载荷进行修正,修正系数计算公式如下:
Kp=Nr/Nc
Nc=N0×Zb×Zp
N0=M0×ω
ω=2π×nr/60
M0=A0/360
式中Kp为压力载荷修正系数;
Nr为额定出力;
Nc为计算出力;
N0为单个水斗出力;
Zb为水斗数;
Zp为喷管数;
ω为角速度;
nr为额定转速;
M0为单个水斗平均力矩;
A0为单个水斗力矩曲线积分。
6)在水斗过流表面循环加载不同时刻全部水斗表面压力载荷文件并求解,提取结果文件,绘制水斗应力随时间变化的曲线,上述步骤即为冲击式水轮机水斗动态应力分析方法。
技术效果:
与现有技术相比,本发明的有益效果是:
现有的冲击式水轮机水斗应力分析方法,只采用简化的射流力载荷作为水斗应力分析有限元模型载荷,计算精度较差;分析方法简单,不能精确描述水斗在机组运行过程中的应力水平和应力变化情况;随着机组容量和设计水头的日益增加,冲击式水轮机水斗的设计难度越来越大,这种分析方法已经不能满足当前设计需求。
本发明利用流体计算软件分析水斗周围流场的压力分布,提取水斗过流表面动态压力场数据,编写程序进行数据转换,生成水斗结构分析压力载荷文件并加以修正,循环加载求解,绘制水斗动态应力变化曲线。该方法的优点是:以压力场数据作为结构分析载荷,并考虑了压力场载荷的修正,计算精度较高。利用该方法可以得到冲击式水轮机运行过程中,水斗的应力随时间的变化情况,对掌握水斗工作受力状态以及评估水斗的动应力和静应力水平具有重要意义。
附图说明:
图1为冲击式水轮机水斗动态应力分析方法流程图;
图2为冲击式水轮机水斗动态应力分析程序流程图;
图3为冲击式水轮机转轮轴向视图;
图4为冲击式水轮机转轮水斗主视图;
图5为单个水斗力矩曲线;
图6为水斗动态应力变化曲线。
具体实施方式:
如图1所示,冲击式水轮机水斗动态应力分析方法,操作步骤如下:
1)运用商业软件建立冲击式水轮机水斗周围流体的数学模型,设置计算域的计算方程和边界条件等计算参数,求解并提取水斗过流表面的动态压力场数据;
2)运用商业软件建立冲击式水轮机水斗结构的有限元模型,在水斗连接螺栓分布圆施加固定约束,对水斗转动中心施加转速载荷;
3)编写数据转换程序,利用流体压力场数据获得水斗动态应力分析压力载荷文件;
4)在水斗过流表面循环加载不同时刻单个水斗表面压力载荷文件并求解,提取结果文件,绘制单个水斗力矩曲线;
5)对水斗压力载荷进行修正,修正系数计算公式如下:
Kp=Nr/Nc
Nc=N0×Zb×Zp
N0=M0×ω
ω=2π×nr/60
M0=A0/360
式中Kp为压力载荷修正系数;
Nr为额定出力;
Nc为计算出力;
N0为单个水斗出力;
Zb为水斗数;
Zp为喷管数;
ω为角速度;
nr为额定转速;
M0为单个水斗平均力矩;
A0为单个水斗力矩曲线积分。
6)在水斗过流表面循环加载不同时刻全部水斗表面压力载荷文件并求解,提取结果文件,绘制水斗应力随时间变化的曲线,上述步骤即为冲击式水轮机水斗动态应力分析方法。
如图2所示,冲击式水轮机水斗动态应力分析程序流程图,编写数据转换程序,利用流体压力场数据获得水斗动态应力分析压力载荷文件。
如图3所示,冲击式水轮机转轮轴向视图,在水斗连接螺栓分布圆施加固定约束,对水斗转动中心施加转速载荷,在水斗过流表面施加压力场载荷。
如图4所示,冲击式水轮机转轮水斗主视图,在水斗过流表面施加压力场载荷。
如图5所示,单个水斗力矩曲线,通过对曲线积分,并利用机组参数计算压力场载荷修正系数。
如图6所示,水斗动态应力变化曲线,利用曲线数据评估机组运行过程中水斗的动应力和静应力水平。

Claims (1)

1.一种冲击式水轮机水斗动态应力分析方法,其特征是:该操
作步骤如下:
1)运用商业软件建立冲击式水轮机水斗周围流体的数学模型,设置计算域的计算方程和边界条件等计算参数,求解并提取水斗过流表面(3)的动态压力场数据;
2)运用商业软件建立冲击式水轮机水斗结构的有限元模型,在水斗连接螺栓分布圆(1)施加固定约束,对水斗转动中心(2)施加转速载荷;
3)编写数据转换程序,利用流体压力场数据获得水斗动态应力分析压力载荷文件;
4)在水斗过流表面(3)循环加载不同时刻单个水斗表面压力载荷文件并求解,提取结果文件,绘制单个水斗力矩曲线;
5)对水斗压力载荷进行修正,修正系数计算公式如下:
Kp=Nr/Nc
Nc=N0×Zb×Zp
N0=M0×ω
ω=2π×nr/60
M0=A0/360
式中Kp为压力载荷修正系数;
Nr为额定出力;
Nc为计算出力;
N0为单个水斗出力;
Zb为水斗数;
Zp为喷管数;
ω为角速度;
nr为额定转速;
M0为单个水斗平均力矩;
A0为单个水斗力矩曲线积分;
6)在水斗过流表面(3)循环加载不同时刻全部水斗表面压力载荷文件并求解,提取结果文件,绘制水斗应力随时间变化的曲线,上述步骤即为冲击式水轮机水斗动态应力分析方法。
CN201510056000.8A 2015-02-03 2015-02-03 冲击式水轮机水斗动态应力分析方法 Active CN104636549B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510056000.8A CN104636549B (zh) 2015-02-03 2015-02-03 冲击式水轮机水斗动态应力分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510056000.8A CN104636549B (zh) 2015-02-03 2015-02-03 冲击式水轮机水斗动态应力分析方法

Publications (2)

Publication Number Publication Date
CN104636549A true CN104636549A (zh) 2015-05-20
CN104636549B CN104636549B (zh) 2017-10-10

Family

ID=53215290

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510056000.8A Active CN104636549B (zh) 2015-02-03 2015-02-03 冲击式水轮机水斗动态应力分析方法

Country Status (1)

Country Link
CN (1) CN104636549B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108021760A (zh) * 2017-12-08 2018-05-11 中国电建集团中南勘测设计研究院有限公司 一种基于数值模拟的流体中曲面受力分布计算方法
CN111120181A (zh) * 2019-12-17 2020-05-08 河海大学 一种测量冲击式水轮机表面动应力试验系统及方法
CN112711882A (zh) * 2020-12-29 2021-04-27 重庆建筑科技职业学院 一种冲击式水轮机转轮疲劳破坏模型构建方法
CN113868796A (zh) * 2021-09-28 2021-12-31 哈尔滨电机厂有限责任公司 一种冲击式水轮机水斗根部过渡自适应方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101404043A (zh) * 2008-09-10 2009-04-08 陕西电力科学研究院 水轮机及其有压过水系统仿真方法
US20110091325A1 (en) * 2009-10-18 2011-04-21 Faris Sadeg M Tool and Method for Rapid Design and Reduction of Rotor Mass
CN103853884A (zh) * 2014-02-24 2014-06-11 昆明理工大学 一种水轮机活动导叶振动特性预测方法
CN104298806A (zh) * 2014-08-13 2015-01-21 浙江工业大学 一种水电站动力特性计算机辅助测试方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101404043A (zh) * 2008-09-10 2009-04-08 陕西电力科学研究院 水轮机及其有压过水系统仿真方法
US20110091325A1 (en) * 2009-10-18 2011-04-21 Faris Sadeg M Tool and Method for Rapid Design and Reduction of Rotor Mass
CN103853884A (zh) * 2014-02-24 2014-06-11 昆明理工大学 一种水轮机活动导叶振动特性预测方法
CN104298806A (zh) * 2014-08-13 2015-01-21 浙江工业大学 一种水电站动力特性计算机辅助测试方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
E.帕金森等: ""冲击式转轮的现代设计技术"", 《水利水电快报》 *
李庆刚等: ""水斗式水轮机喷嘴及水斗应力和振形研究"", 《农业机械学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108021760A (zh) * 2017-12-08 2018-05-11 中国电建集团中南勘测设计研究院有限公司 一种基于数值模拟的流体中曲面受力分布计算方法
CN111120181A (zh) * 2019-12-17 2020-05-08 河海大学 一种测量冲击式水轮机表面动应力试验系统及方法
CN112711882A (zh) * 2020-12-29 2021-04-27 重庆建筑科技职业学院 一种冲击式水轮机转轮疲劳破坏模型构建方法
CN112711882B (zh) * 2020-12-29 2024-02-27 重庆建筑科技职业学院 一种冲击式水轮机转轮疲劳破坏模型构建方法
CN113868796A (zh) * 2021-09-28 2021-12-31 哈尔滨电机厂有限责任公司 一种冲击式水轮机水斗根部过渡自适应方法
CN113868796B (zh) * 2021-09-28 2024-04-09 哈尔滨电机厂有限责任公司 一种冲击式水轮机水斗根部过渡自适应方法

Also Published As

Publication number Publication date
CN104636549B (zh) 2017-10-10

Similar Documents

Publication Publication Date Title
Quaranta et al. Output power and power losses estimation for an overshot water wheel
Yelmule et al. CFD predictions of NREL phase VI rotor experiments in NASA/AMES wind tunnel
CN104636549A (zh) 一种冲击式水轮机水斗动态应力分析方法
Campobasso et al. Effects of flow instabilities on the linear analysis of turbomachinery aeroelasticity
Zhou et al. Investigation of rotating stall for a centrifugal pump impeller using various SGS models
CN106529048A (zh) 空调压缩机配管振动仿真方法
CN104298869A (zh) 一种弹性水翼的流固耦合特性数值预测方法
Kan et al. Dynamic stress of impeller blade of shaft extension tubular pump device based on bidirectional fluid-structure interaction
CN105260580A (zh) 一种确定混流式水轮机尾水管压力脉动幅值的方法
CN111159941A (zh) 一种用于汽车液力变矩器内流场瞬态数值模拟的方法
CN111523182A (zh) 一种基于流固耦合的叶片振动响应分析方法
Zhang et al. Application of the hybrid RANS/LES method on the hydraulic dynamic performance of centrifugal pumps
Li et al. Stability optimization and analysis of a bidirectional shaft extension pump
CN105653767A (zh) 水轮机导叶水力矩的测控方法
Sharma et al. Numerical simulation for pressure distribution in pelton turbine nozzle for the different shapes of spear
CN112283051A (zh) 一种基于升力线模型的振动信号特征优化方法及系统
Huang et al. A new setting criterion of tailrace surge chambers for pumped-storage power plants
Bahod et al. Performance Improvement of Pump as Turbine by Back Cavity Filling through Simulational Approach
Amano et al. Multi-Disciplinary Design Optimization of a Horizontal Micro Kaplan Hydro Turbine
CN202420564U (zh) 高压流量自控仪反推圆弧叶轮
Saeed et al. Complete Francis turbine flow simulation at Derbendikan power station
CN104422820B (zh) 能量回馈抽油机真实电能输入功率的测试计算方法
Hu et al. Computation of stress distribution in a mixed flow pump based on fluid-structure interaction analysis
He et al. Numerical Simulation of Fluid-Structure Interaction of Anisotropic Composite Propeller
Shi et al. Design optimization of axial-flow pump blades based on iSIGHT

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant