CN101929856B - 一种测量热处理过程水轮机叶片翘曲变形的方法 - Google Patents

一种测量热处理过程水轮机叶片翘曲变形的方法 Download PDF

Info

Publication number
CN101929856B
CN101929856B CN2010102458124A CN201010245812A CN101929856B CN 101929856 B CN101929856 B CN 101929856B CN 2010102458124 A CN2010102458124 A CN 2010102458124A CN 201010245812 A CN201010245812 A CN 201010245812A CN 101929856 B CN101929856 B CN 101929856B
Authority
CN
China
Prior art keywords
blade
point
blades
coordinate
buckling deformation
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.)
Active
Application number
CN2010102458124A
Other languages
English (en)
Other versions
CN101929856A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN2010102458124A priority Critical patent/CN101929856B/zh
Publication of CN101929856A publication Critical patent/CN101929856A/zh
Application granted granted Critical
Publication of CN101929856B publication Critical patent/CN101929856B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

一种测量热处理过程水轮机叶片翘曲变形的方法,叶片入炉前,在叶片的根部区域中,选择三个点A、B、C点为入炉前标记点,建立任意的三维坐标系OXYZ,分别测量出点A、B、C的三维坐标,然后在叶片上选取其它点进行测量,所测位置至少覆盖叶片的出水边表面位置,以所测各点的数据建立叶片t0时刻的表面形状;对叶片加热出炉后在炉外冷却过程中不同时间段叶片形状进行测量,获取各时间段叶片上各点的三维坐标,以所测各点的三维坐标数据建立叶片tj时刻的表面形状;建立新的坐标系O′X′Y′Z′,将任意一点坐标转化成在新坐标系O′X′Y′Z′下的坐标;最后统一所有测量数据的坐标系统,建立叶片各个时间段的形状,对比分析各时间段叶片形状变化规律,从而获取叶片翘曲变形情况。

Description

一种测量热处理过程水轮机叶片翘曲变形的方法
技术领域
本发明属于大型水轮机制造技术领域,涉及热处理过程叶片翘曲变形控制,具体涉及一种测量热处理过程水轮机叶片翘曲变形的方法。
背景技术
随着社会对清洁能源需求的增加,发展水电是其中的一个重要方向。水轮机是水力发电的基本设备,而叶片是水轮机的重要组成部件,它的质量关系到发电的效率、功率与水轮机的运行寿命。在叶片生产过程中,热处理工艺极为关键,不仅决定叶片最终的组织分布,也决定机加工前的几何形状。在热处理中,由于叶片厚薄分布差别较大,温度冷却不均匀,导致叶片内部产生热应力和相变应力,在它们的综合作用下,叶片容易发生翘曲变形。如果叶片翘曲变形控制不好,其可能需要增加补焊或校形等额外工序,增加生产成本。因此,掌握叶片的翘曲变形规律,能够为叶片生产提供合理的科学依据,对企业具有重要的意义。然而,至今仍没有相关的技术用于测量水轮机叶片三维整体的翘曲变形行为。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种测量热处理过程水轮机叶片翘曲变形的方法,该方法可准确测量热处理过程叶片的翘曲变形情况,为水轮机叶片生产中的变形控制提供指导,实现叶片形状的精确控制,避免补焊、校形等补救工艺,减少机加工余量。
为了解决上述技术问题,本发明采用的技术方案是:一种测量热处理过程水轮机叶片翘曲变形的方法,包括以下步骤:
第一步,叶片入炉前,在叶片的根部区域1中,选择三个点A点、B点和C点为入炉前标记点;
第二步,建立任意的三维坐标系OXYZ,分别测量出A点、B点和C点的三维坐标
Figure BSA00000217587900021
然后在叶片上选取其它点进行测量,所测位置至少覆盖叶片的出水边表面位置,所测坐标依次为
Figure BSA00000217587900022
以所测各点的数据
Figure BSA00000217587900023
(i=1,2,3……N),建立t0时刻叶片的表面形状;
第三步,对叶片加热出炉后在炉外冷却过程中不同时间段叶片形状进行测量,获取风冷过程各时间段叶片上各点的三维坐标
Figure BSA00000217587900024
Figure BSA00000217587900025
以所测各点的三维坐标数据
Figure BSA00000217587900027
建立叶片tj时刻的表面形状;
第四步,建立新的坐标系O′X′Y′Z′,设定A点位于原点O′,以向量
Figure BSA00000217587900028
的方向为X′轴的方向,向量
Figure BSA00000217587900029
与向量
Figure BSA000002175879000210
做叉乘得到向量
Figure BSA000002175879000211
以该向量的方向作为Z′轴的方向,由已得到的X′轴与Z′轴,按照右手坐标系法则,建立直角坐标系O′X′Y′Z′;
第五步,将任意一点坐标
Figure BSA000002175879000212
转化成在新坐标系O′X′Y′Z′下的坐标
Figure BSA000002175879000213
第六步,统一所有测量数据的坐标系统,建立叶片各个时间段的形状,对比分析各时间段叶片形状变化规律,从而获取叶片翘曲变形情况。
所述第二步中对叶片其他位置进行测量,测量的点数N在100个以上。
所述第四步中获取风冷过程各时间段叶片上各点的三维坐标
Figure BSA000002175879000214
所选取的时间段在3个以上。
所述测量三维坐标的方式为激光检测。
本发明与现有技术相比具有以下优点:
1)可通过多个点的测量实现面的测量,可在热处理过程中以一定时间间隔测量,从而准确反映整个叶片在热处理炉外冷却过程中的动态翘曲变形;还可通过比较叶片入炉前和出炉后的形状,得到叶片在热处理炉内的变形;
2)过程科学,非接触式,成本较低。
附图说明
图1为本发明叶片表面相对位置的三维测量。
图2为本发明叶片翘曲变形测量;其中:2为覆盖了叶片出水边表面位置的区域。
图3为本发明坐标转换示意图。
图4为本发明的效果示意图;其中:3为第一时间叶片测量结果;4为第二时间叶片测量结果。
具体实施方式
下面结合附图和公式详细说明本发明的实施方式。
一种测量热处理过程水轮机叶片翘曲变形的方法,其特征在于,包括以下步骤:
第一步,叶片入炉前,在叶片的根部区域1中,选择三个点A点、B点和C点为入炉前标记点;
第二步,建立任意的三维坐标系OXYZ,分别测量出A点、B点和C点的三维坐标然后在叶片上选取其它点进行测量,所测位置至少覆盖叶片的出水边表面位置,所测坐标依次为以所测各点的数据
Figure BSA00000217587900033
建立叶片t0时刻的表面形状;
第三步,对叶片加热出炉后在炉外冷却过程中不同时间段叶片形状进行测量,获取风冷过程各时间段叶片上各点的三维坐标
Figure BSA00000217587900034
Figure BSA00000217587900041
Figure BSA00000217587900042
以所测各点的三维坐标数据
Figure BSA00000217587900043
建立叶片tj时刻的表面形状;
第四步,建立新的坐标系O′X′Y′Z′,设定A点位于原点O′,以向量的方向为X′轴的方向,向量
Figure BSA00000217587900045
与向量
Figure BSA00000217587900046
做叉乘得到向量
Figure BSA00000217587900047
以该向量的方向作为Z′轴的方向,由已得到的X′轴与Z′轴,按照右手坐标系法则,建立直角坐标系O′X′Y′Z′;
第五步,将任意一点坐标
Figure BSA00000217587900048
转化成在新坐标系O′X′Y′Z′下的坐标
第六步,统一所有测量数据的坐标系统,建立叶片各个时间段的形状,对比分析各时间段叶片形状变化规律,从而获取叶片翘曲变形情况。
所述第二步中对叶片其他位置进行测量,测量的点数在100个以上。
所述第四步中获取风冷过程各时间段叶片上各点的三维坐标所选取的时间段在3个以上。
所述测量三维坐标的方式为激光检测。
其中所述第五步中将任意一点坐标
Figure BSA000002175879000411
转化成在新坐标系O′X′Y′Z′下的坐标
Figure BSA000002175879000412
其转换方法如下:
X ′ i , t j Y ′ i , t j Z ′ i , t j = a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 X i , t j - X 1 , t j Y i , t j - Y 1 , t j Z i , t j - Z 1 , t j
其中
a 11 = X 2 , t j - X 1 , t j l 12 , t j ,
a 12 = Y 2 , t j - Y 1 , t j l 12 , t j
a 13 = Z 2 , t j - Z 1 , t j l 12 , t j
a21=a32a13-a33a12
a22=a32a13-a33a12
a23=a32a13-a33a12
a 31 = ( Y 2 , t j - Y 1 , t j ) ( Z 3 , t j - Z 1 , t j ) - ( Z 2 , t j - Z 1 , t j ) ( Y 3 , t j - Y 1 , t j ) w
a 32 = ( Z 2 , t j - Z 1 , t j ) ( X 3 , t j - X 1 , t j ) - ( X 2 , t j - X 1 , t j ) ( Z 3 , t j - Z 1 , t j ) w
a 32 = ( X 2 , t j - X 1 , t j ) ( Y 3 , t j - Y 1 , t j ) - ( Y 2 , t j - Y 1 , t j ) ( X 3 , t j - X 1 , t j ) w
其中,
l 12 , t j = ( X 2 , t j - X 1 , t j ) 2 + ( Y 2 , t j - Y 1 , t j ) 2 + ( Z 2 , t j - Z 1 , t j ) 2
w = ( ( ( Y 2 , t j - Y 1 , t j ) ( Z 3 , t j - Z 1 , t j ) - ( Z 2 , t j - Z 1 , t j ) ( Y 3 , t j - Y 1 , t j ) ) 2
+ ( ( Z 2 , t j - Z 1 , t j ) ( X 3 , t j - X 1 , t j ) - ( X 2 , t j - X 1 , t j ) ( Z 3 , t j - Z 1 , t j ) ) 2
+ ( ( X 2 , t j - X 1 , t j ) ( Y 3 , t j - Y 1 , t j ) - ( Y 2 , t j - Y 1 , t j ) ( X 3 , t j - X 1 , t j ) ) 2 ) 1 / 2
本发明的主要原理是以某点为相对坐标原点建立三维坐标系,通过红外检测准确测量叶片表面各点的三维坐标(i=1,2,……N;j=0,1,2,……M),其中i为特征点顺序号,为了保证测量精确度,N≥100,tj为测量时间的标记,M≥3。在测量过程中,各点测量误差均±1mm。叶片表面相对位置三维测量如图1所示。
图2为本发明叶片翘曲变形测量,根部区域1在叶片的根部位置,叶片在该位置处相对厚大,翘曲变形非常小,为叶片不容易翘曲的区域,在根部区域1内选择三个点A、B、C作为入炉前标记点,用来统一坐标系,它们被测相对坐标分别为
Figure BSA000002175879000511
然后,在叶片上选取其它点进行测量,获取其坐标
Figure BSA000002175879000512
所取的点的数量至少为100个并保证包含叶片的出水边表面位置。以所测的各点的数据
Figure BSA00000217587900061
建立叶片tj时刻的表面形状。在此基础上,以
Figure BSA00000217587900062
Figure BSA00000217587900063
为空间坐标,对叶片形状进行转换,建立一个统一的坐标系统,如图3所示。
通过测量不同时间段叶片的形状,将它们进行比较,对比的实际效果如图4所示,可以明显地反映出在热处理过程中的翘曲变形行为。

Claims (4)

1.一种测量热处理过程水轮机叶片翘曲变形的方法,其特征在于,包括以下步骤:
第一步,叶片入炉前,在叶片的根部区域(1)中,选择三个点A点、B点和C点为入炉前标记点;
第二步,建立任意的三维坐标系OXYZ,分别测量出A点、B点和C点的三维坐标 ( X 1 , t 0 , Y 1 , t 0 , Z 1 , t 0 ) , ( X 2 , t 0 , Y 2 , t 0 , Z 2 , t 0 ) , ( Z 3 , t 0 , Y 3 , t 0 , Z 3 , t 0 ) , 然后在叶片上选取其它点进行测量,所测位置至少覆盖叶片的出水边表面位置,所测坐标依次为 ( X 4 , t 0 , Y 4 , t 0 , Z 4 , t 0 ) , ( X 5 , t 0 , Y 5 , t 0 , Z 5 , t 0 ) , ( Z 6 , t 0 , Y 6 , t 0 , Z 6 , t 0 ) . . . . . . ( X N , t 0 , Y N , t 0 , Z N , t 0 ) , 以所测各点的数据
Figure FSB00000736921000013
(i=1,2,3……N),建立叶片t0时刻的表面形状;
第三步,对叶片加热出炉后在炉外冷却过程中不同时间段叶片形状进行测量,获取风冷过程各时间段叶片上各点的三维坐标 ( X 1 , t 1 , Y 1 , t 1 , Z 1 , t 1 ) , ( X 2 , t 1 , Y 2 , t 1 , Z 2 , t 1 ) , ( Z 3 , t 1 , Y 3 , t 1 , Z 3 , t 1 ) . . . . . . ( Z N , t 1 , Y N , t 1 , Z N , t 1 ) , ( Z 1 , t 2 , Y 1 , t 2 , Z 1 , t 2 ) , ( Z 2 , t 2 , Y 2 , t 2 , Z 2 , t 2 ) , ( Z 3 , t 2 , Y 3 , t 2 , Z 3 , t 2 ) . . . . . . ( X N , t 2 , Y N , t 2 , Z N , t 2 ) , . . . . . . , ( Z 1 , t j , Y 1 , t j , Z 1 , t j ) , ( Z 2 , t j , Y 2 , t j , Z 2 , t j ) , ( X 3 , t j , Y 3 , t j , Z 3 , t j ) . . . . . . ( X N , t j , Y N , t j , Z N , t j ) , 以所测各点的三维坐标数据建立叶片tj时刻的表面形状;
第四步,建立新的坐标系O′X′Y′Z′,设定A点位于原点O′,以向量的方向为X′轴的方向,向量与向量做叉乘得到向量以该向量的方向作为Z′轴的方向,由已得到的X′轴与Z′轴,按照右手坐标系法则,建立直角坐标系O′X′Y′Z′;
第五步,将任意一点坐标
Figure FSB000007369210000112
转化成在新坐标系O′X′Y′Z′下的坐标 ( X ′ i , t j , Y ′ i , t j , Z ′ i , t j ) ;
第六步,统一所有测量数据的坐标系统,建立叶片各个时间段的形状,对比分析各时间段叶片形状变化规律,从而获取叶片翘曲变形情况。
2.根据权利要求1所述的一种测量热处理过程水轮机叶片翘曲变形的方法,其特征在于,所述第二步中对叶片其他位置进行测量,测量的点数N在100个以上。
3.根据权利要求1所述的一种测量热处理过程水轮机叶片翘曲变形的方法,其特征在于,所述第三步中获取风冷过程各时间段叶片上各点的三维坐标
Figure FSB00000736921000021
所选取的时间段在3个以上。
4.根据权利要求1所述的一种测量热处理过程水轮机叶片翘曲变形的方法,其特征在于,所述测量三维坐标的方式为激光检测。
CN2010102458124A 2010-08-04 2010-08-04 一种测量热处理过程水轮机叶片翘曲变形的方法 Active CN101929856B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102458124A CN101929856B (zh) 2010-08-04 2010-08-04 一种测量热处理过程水轮机叶片翘曲变形的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102458124A CN101929856B (zh) 2010-08-04 2010-08-04 一种测量热处理过程水轮机叶片翘曲变形的方法

Publications (2)

Publication Number Publication Date
CN101929856A CN101929856A (zh) 2010-12-29
CN101929856B true CN101929856B (zh) 2012-07-25

Family

ID=43369167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102458124A Active CN101929856B (zh) 2010-08-04 2010-08-04 一种测量热处理过程水轮机叶片翘曲变形的方法

Country Status (1)

Country Link
CN (1) CN101929856B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2679778B1 (en) * 2012-06-27 2019-08-21 Ansaldo Energia IP UK Limited A method for measuring geometry deformations of a turbine component
US10345179B2 (en) * 2017-02-14 2019-07-09 General Electric Company Passive strain indicator
CN113155023B (zh) * 2021-04-02 2023-02-07 甘肃旭盛显示科技有限公司 液晶基板玻璃翘曲度测量方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101315270A (zh) * 2007-05-31 2008-12-03 通用电气公司 用于测量涡轮中叶片变形的方法和系统
RU2379626C1 (ru) * 2008-06-16 2010-01-20 Институт проблем управления сложными системами Российской академии наук Способ оценки деформации статора и параметров биения ротора газотурбинного двигателя
CN101696875A (zh) * 2009-11-04 2010-04-21 清华大学 大型铸件在热处理炉外冷却过程中的变形动态检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55131711A (en) * 1979-04-02 1980-10-13 Hitachi Ltd Solid object deformation measuring method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101315270A (zh) * 2007-05-31 2008-12-03 通用电气公司 用于测量涡轮中叶片变形的方法和系统
RU2379626C1 (ru) * 2008-06-16 2010-01-20 Институт проблем управления сложными системами Российской академии наук Способ оценки деформации статора и параметров биения ротора газотурбинного двигателя
CN101696875A (zh) * 2009-11-04 2010-04-21 清华大学 大型铸件在热处理炉外冷却过程中的变形动态检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JP昭55-131711A 1980.10.13
王立民 等.403Nb叶片钢热压缩变形研究.《钢铁》.2010,第45卷(第6期),63-66. *
程建国 等.从凝固过程的位移数值模拟结果中提取铸件的变形.《铸造技术》.2008,第29卷(第10期),第1322-1326页. *

Also Published As

Publication number Publication date
CN101929856A (zh) 2010-12-29

Similar Documents

Publication Publication Date Title
CN101672637B (zh) 一种复杂曲面的数字化检测方法
CN102735204B (zh) 一种基于弦线的航空薄壁叶片加工扭曲度误差测量方法
CN102778403B (zh) 一种焊缝材料参数识别方法
CN103246774B (zh) 数值模拟p92钢管焊接热影响区宽度的方法
CN104268343A (zh) 一种用于端铣切削的切削力预测及温度预测的方法
CN106354094A (zh) 基于空间标准球的机床随动激光扫描坐标标定方法
CN102288123B (zh) 一种用于焊接变形三维非接触式测量的精确定位方法
CN105760641B (zh) 面向非标准型风场的风电机组疲劳寿命高效评估方法
CN102554705B (zh) 一种光学自由曲面补偿加工方法
CN101929856B (zh) 一种测量热处理过程水轮机叶片翘曲变形的方法
CN102728932A (zh) 一种自动优化电弧焊工艺参数的方法及其装置
CN104776827B (zh) Gps高程异常数据的粗差探测方法
CN103617369B (zh) 超大型曲面零部件加工精度评估方法和设备
CN105404350A (zh) 一种基于二分梯度法的热电发电mppt方法
CN109612441B (zh) 一种基于最优分类和矢量平均的波浪方向浮标主波向计算方法
Janjua et al. Blade profile optimization of kaplan turbine using CFD analysis
CN105334873A (zh) 一种太阳能定日镜转动精度的检测方法
CN201514207U (zh) 一种复杂曲面的数字化检测系统
Auger et al. Investigating the quality inspection process of offshore wind turbine blades using B-spline surfaces
CN105512441A (zh) 一种物体内热源分布重构系统及方法
CN103399273B (zh) 波浪能装置实海况测试方法
CN103278118B (zh) 确定闭式整体叶盘三坐标测量机测头角度的方法
CN105242546B (zh) 一种基于材料特性的高速铣削力建模方法
CN104657548A (zh) 一种平板裂缝阵天线辐射阵面误差的建模方法
CN101498572A (zh) 凸轮二维曲线测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant