考虑结合面法向刚度的拉杆转子轴向振动特性

第２９卷第２４期
中国机械工程
V o l ．２９㊀N o ．２４２０１８年１２月
C H I N A M E C HA N I C A LE N G I N E E R I N G
p p
．２９５４Ｇ２９５８考虑结合面法向刚度的拉杆转子轴向振动特性
何㊀谦１,２
㊀王艾伦２㊀陈中祥１㊀杨㊀俊１
１．湖南师范大学工程与设计学院,长沙,４１００８１
２．中南大学高性能复杂制造国家重点实验室,长沙,４１００８３
摘要:针对接触刚度解析模型参数确定困难㊁精度难以保证等问题,提出了根据弹塑性粗糙表面微体单元受力变形的有限元分析结果确定轮盘结合面法向刚度的方法;为准确获取拉杆转子的轴向振动特性,建立了考虑轮盘结合面法向刚度的集中质量动力学模型;运用上述方法和模型计算了某型实验转子轴向振动的固有频率,并将结果与实测数据进行对比,误差低于５％,证明了该方法的有效性;改变拉杆预紧力,进一步研究预紧力对拉杆转子动力学行为的影响,结果表明:拉杆预紧力对转子的作用效果存在一个饱和区域,可为拉杆预紧力数值的确定提供重要的设计依据
.
关键词:拉杆转子;结合面法向刚度;轴向振动;拉杆预紧力中图分类号:T H １２４
D O I :１０．３９６９/j ．i s s n ．１００４ １３２X．２０１８．２４．００９开放科学(资源服务)标识码(O S I D )
:T h eA x i a lV i b r a t i o nC h a r a c t e r i s t i c s o fR o dF a s t e n i n g R
o t o rC o n s i d e r e d C o n t a c tN o r m a l S t i f f n e s s
H E Q i a n １,
２㊀WA N G A i l u n ２㊀C H E NZ h o n g x i a n g １㊀Y A N GJ
u n １
１．E n g i n e e r i n g &D e s i g nC o l l e g e ,H u n a nN o r m a lU n i v e r s i t y ,C h a n g
s h a ,４１００８１２．S t a t eK e y L a b o r a t o r y o fH i g hP e r f o r m a n c eC o m p l e x M a n u f a c t u r i n g
,C e n t r a l S o u t hU n i v e r s i t y ,C h a n g
s h a ,４１００８３A b s t r a c t :T h e e x p e r i m e n t a l a n d s t a t i s t i c a l p a r a m e t e r sw h i c h r e l i e d o n t h e a n a l y t i c a lm o d e l s i s h a r d t od e t e r m i n e ,t h e a c c u r a c y i s d i f f i c u l t t o g u a r a n t e e a n d o t h e r i s s u e s ．A m e t h o d o f f i n i t e e l e m e n t a n a l y
Ｇs i s o f t h eM i c r o Ｇe l e m e n t D e f o r m a t i o n o f t h e e l a s t i c Ｇp l a s t i c r o u g h s u r f a c ew a s p r o p o s e d h e r e i n t o d e t e r Ｇm i n e t h en o r m a l s t i f f n e s s o f t h e d i s ks u r f a c e ．I no r d e r t o a c c u r a t e l y o
b t a i n t h e a x i a l v i b r a t i o n
c h a r a c Ｇt e r i s t i c s o f t h e r o
d r o t o r ,a c o n c
e n t r a t e dm a s sd y n a m i cm o d e l c o n s i d e r i n g t
h en o r m a l s t i f f n e s so f t h e r o u l e t t e j o i n tw a s e s t a b l i s h e d ．T h en a t u r a l f r e q u e n c y o f t h e a x i a l v i b r a t i o no f a c e r t a i n t y p e o f e x p e r i Ｇm e n t a l r o t o rw a s c a l c u l a t e db y u s i n g t h ea b o v em e t h o da n d m o d e l ．T h e r e s u l t sw e r ec o m p
a r e dw i t h t h em e a s u r e dd a t a ,a n dt h ee r r o r i s l e s s t h a n５％,w h i c hd e m o n s t r a t e s t h ev a l i d i t y o
f t h i s m e t h o d ．T h e e f f e c to f t h e p r e Ｇt i
g
h t e n
i n g f o r c eo nt h ed y n a m i cb e h a v i o ro f t h er o t o rw a sf u r t h e rs t u d i e db y c h a n g i n g t h e v a l u e o f t h e p r e Ｇt i g h t e n i n g f o r c e ．T h e r e s u l t s s h o wt h a t t h e r e i s as a t u r a t i o nr e g
i o n f o r t h e e f f e c t o f t h e r o d p r e Ｇt i g h t e n i n g f o r c eo n t h e r o t o r ,w h i c h p r o v i d e s i m p o r t a n td e s i g
nb a s i s f o r t h e d e t e r m i n a t i o no f t h e r o d p r e Ｇt i g h t e n i n g f
o r c e ．K e y w
o r d s :r o df a s t e n i n g r o t o r ;n o r m a l s t i f f n e s so f j u n c t i o ns u r f a c e ;a x i a lv i b r a t i o n ;r o d p r e Ｇt i g h t e n i n g f
o r c e 收稿日期:２０１７１１２７基金项目:国家重点基础研究发展计划(９７３计划)
资助项目(２０１３C B ０３５７０６
)０㊀引言
拉杆转子具有质量小㊁刚度大㊁加工容易以及
便于冷却和维修等突出优势,被广泛应用于航空
发动机及大功率重型燃气轮机中[
１
].不同于整体转子,拉杆转子通过周向均匀分布的多根拉杆将轴头和轮盘预紧连接而成,因而在结构上存在多个非连续的接触界面,物理特性非常复杂.由于接触问题的本质是非线性的,零件的材质㊁表面形貌㊁载荷情况等因素都会对结合面的接触状态和
接触刚度造成影响[２]
,很难用一个统一的模型加
以描述,这给拉杆转子的深入研究带来了巨大的困难.现阶段对接触问题的研究主要有两种方法,即建立解析模型和开展数值计算.解析的方法通常依托G r e e n w o o d 等在１９６６年提出的GW 模型以及以其为基础的各种改良模型,如M B 分
形接触模型㊁GW 修正模型等.以上模型都基于某些假设,而且在参数的确定上需要依托实验和
统计的数据,操作起来比较困难[３
],结果的精度也
难以保证.随着计算机技术的飞速发展,以有限元为代表的数值计算方法在接触分析中得到了越
来越广泛的应用[４]
,其参数化的设计方法不仅能
４５９２ Copyright©博看网 . All Rights Reserved.
方便高效地获取实验样本,而且更容易发现或形成规律性的认识和结论.将其运用于拉杆转子的动力学特性分析,成效显著[５Ｇ７].
由于弯曲振动是拉杆转子动力学特性中最为重要的环节,故在业已开展的研究中,绝大部分都是针对这一问题进行的.事实上,拉杆转子的轴向振动问题也不容忽视.据报道,２００１年广东月亮湾燃机电厂燃气轮机发电机组投产一周后,就出现了机组转子因轴向振动剧烈而导致的故障[８].早在２０世纪５０年代,燃气轮机转子的轴向振动问题就引起人们注意[９],但由于检测困难等原因,研究未能深入.目前,国内外对转子轴向振动问题的研究相对较少,而关于拉杆转子轴向振动问题的研究更少.文献[１０]用传递矩阵法对多轮盘转子系统轴向振动的固有特性进行了研究,但其研究对象为整体转子.文献[１１]运用键合图的方法建立了拉杆转子轴向振动的动力学模型,但该方法界面接触刚度计算方面主要参照了文献[１２]的结论,给出的轮廓面积比ηc a在某一范围内,实际确定起来比较复杂.
本文在分析和总结现有研究成果的基础上,提出了一种基于有限元分析结果的轮盘结合面法向接触刚度的确定方法,将其与转子的集中质量模型结合,计算了某型实验转子轴向振动的固有频率,通过对比理论计算和实验实测的结果,验证该方法的有效性.
１㊀轮盘结合面法向刚度的确定
１．１㊀三维粗糙表面的生成
T HOMA S等[１３]已证明指数自相关函数能够很好地描述现实世界中的许多随机现象; WH I T E HO U S E等[１４]通过实验证明了许多工程表面轮廓具有指数自相关函数关系.一般给定指数形式的自相关函数为
R(τx,τy)＝σ２e－２．３㊀τxβx()２＋τyβy()２
τx＝１,２, ,N㊀τy＝１,２, ,N
式中,σ为表面粗糙度;βx㊁βy分别为x㊁y方向上的相关长度.
根据轮盘表面的粗糙度数值,按指数自相关值,以粗糙度R a＝３．２μm为例,在MA T
L A B中生成相应的微观表面,见图１.
１．２㊀基于微体单元的结合面法向刚度的确定由于零件的宏观尺寸均远远大于用以描述表面形貌的参数尺寸,故如果考虑表面的微观形貌直接对零件开展有限元分析,将会导致单元和节
图１㊀微体单元粗糙表面形貌
F i g．１㊀R o u g h s u r f a c e f e a t u r e s o f t h em i c r o b o d y
点的数目巨大,可能会因计算量过大而无法进行.当粗糙表面微型长方体的宏观尺寸达到一定数值时,所得的分析结果与继续扩大微体尺寸并无二致[１５],因此,借助具有合适几何尺寸的微体单元,所得的结果完全能反映实际零件的接触情况.据此将MA T L A B中所获得的点云数据导入C A D 软件中,生成一个左端面为粗糙表面的微型长方体,３个方向的尺寸分别为１mmˑ１mmˑ５mm.
考虑到两粗糙表面微元体接触和单个粗糙面与刚性平面接触的分析结果在保证微体单元足够长度后并无差别[１６],本文直接采用相对简单的微体单元与刚性平面接触的分析方法.将C A D中生成的模型导入A N S Y S,网格划分采用１０节点四面体单元,微体单元定义为弹塑性,各项力学指标依据轮盘材料给定.分析时,对长方体的底面加以约束,并限制粗糙表面上位于刚性平面内的节点的位移.在与粗糙表面相对的右端面上,施加一个由转子预紧力转换而成的分布力p.查看分析结果,获得微体单元在该方向的变形量,由此获得整个微体单元在该作用力下的等效法向刚度k e q:
k e q＝p A/Δl
式中,A为微体单元的横截面面积;Δl为微体单元在该方向的变形.
依据圣维南原理,粗糙表面的变形情况对距离其较远区域的影响很小.据此可按有限元软件的分析结果将微体划分为粗糙段和光滑段两部分,粗糙段为微体单元受力时发生塑性变形的区域,用l c表示;余下部分视作只发生弹性变形,用l s表示,见图２.
图２㊀微体单元等效模型
F i g．２㊀E q u i v a l e n tm o d e l o f t h em i c r o b o d y
５５９２
考虑结合面法向刚度的拉杆转子轴向振动特性 何㊀谦㊀王艾伦㊀陈中祥等
Copyright©博看网 . All Rights Reserved.
显然,微体单元在不同载荷作用下的粗糙段l c的长度也不相同.改变压力p,得到l c与p的变化关系曲线,见图３.
图３㊀粗糙段长度与作用力的关系曲线
F i g．３㊀R e l a t i o n c u r v e o f r o u g h s e c t i o n l e n g t h&p r e s s
由材料力学可知,弹性轴段l s的法向刚度
k s＝E A/l s
式中,E为轮盘材料的弹性模量;l s为弹性光滑段的长度.
用k c表示粗糙段的刚度.显然,微体单元各段的刚度k c㊁k s和等效法向刚度k e q三者之间存在以下关系:
１k e q＝１
k c＋１
k s
联立以上各式,即可求得粗糙段的刚度k c.
定义k j为轮盘结合面的法向刚度.考虑到实际情况为两粗糙表面接触,故有
k j＝k c/２
２㊀拉杆转子轴向振动集中参数模型
尽管不同厂家及不同型号的拉杆转子在几何尺寸上差异明显,但转子的组成和结构却基本相同,即均通过长螺栓将轴头和轮盘串连而成.为方便说明和计算,本文对拉杆转子的形态进行了适当简化,其物理结构见图４.
图４㊀拉杆转子结构示意图
F i g．４㊀S c h e m a t i c v i e wo f t h e r o d f a s t e n i n g r o t o r ２．１㊀拉杆转子质量的集中参数化
结合轴向振动分析的要求以及拉杆转子的结构特点,选取单个轮盘作为基本单元进行集中参数化:将轮盘的质量集中到两端,将其变形集中到中间.轮盘之间的接触刚度用一根无质量的弹簧表示.拉杆处理方式与轮盘类似.需要注意的是,由于位居拉杆两端的螺栓头和螺母与轴头没
有相对运动,故可将拉杆的集中质量看作与两轴头外侧的质量单元固连在一起.相应的力学模型见图５.
图５㊀拉杆转子轴向振动集中质量模型
F i g．５㊀L u m p e dＧm a s sm o d e l o f t h e r o d a s t e n i n g
r o t o r a x i a l v i b r a t i o n
图５中,k d㊁k r o d㊁k j分别为轮盘㊁拉杆以及轮盘接合面的法向刚度.k d㊁k r o d的确定不考虑塑性变形的影响,依据式(４)和零件的材料和几何尺寸获得,其中,k r o d为所有拉杆的刚度之和.２．２㊀数学模型的建立
由于转子系统一般都是弱阻尼系统,故在建模和计算过程中可以不计阻尼.依据图５中的集中参数模型,写成矩阵的形式,即
M X ＋K X＝F
X＝x１㊀x２ x２n
[]T㊀F＝f１㊀f２ f２n
[]T
M＝
m r o d
２＋
m D１
２ ０
m D１
２
⋮⋱⋮
m D
n
２
０ m r o d２＋m D n２é
ë
ê
ê
ê
ê
ê
ê
ê
ê
ê
ê
ê
ù
û
ú
ú
ú
ú
ú
ú
ú
ú
ú
ú
ú
K＝
k r o d＋k d－k d ０－k r o d －k d k j１＋k d k j１ ０
⋮⋱⋮
０ －k j n－１k j n－１＋k d－k d －k r o d０ －k d k r o d＋k d é
ë
ê
ê
ê
ê
ê
ê
ù
û
ú
ú
ú
ú
ú
ú据此可以求得拉杆转子轴向振动时的固有频率和振型[１７].
３㊀实验及数据分析
实验用拉杆转子由４轮盘８拉杆组成,每根拉杆的预紧力为４k N,轮盘的大圆直径与中心孔直径分别为１００mm和４０mm,轮盘厚度为３０mm,拉杆的直径为１０mm,轮盘中心与拉杆中心的距离为３５mm.为排除支承的影响,用两
６５９２
中国机械工程第２９卷第２４期２０１８年１２月下半月Copyright©博看网 . All Rights Reserved.
根绕在转子轴头的弹性绳将转子水平自由悬挂.用力锤沿轴线方向敲击转子的端面,用压电加速度计测量响应,并借助信号分析处理系统获得转子的固有频率和振型.
为验证２．２节转子模型的有效性,
将实验转子的相关参数代入该模型进行计算.作为对比,还建立了同形态整体转子的有限元模型.３种渠道所获得的该型转子前３阶轴向振动的固有频率见表１.
表１㊀轴向振动固有频率比较(预紧力为４k N )
T a b ．１㊀C o m p a r i s o n f o r n a t u r a l f r e q u e n c y o
f a x i a l v i b r a t i o n (p r e Ｇt i
g
h t e n
i n g f
o r c e ４k N )阶次拉杆转子模型理论值(H z )整体转子模型理论值(H z )实测值(H z )拉杆转子/整体转子模型误差(％)１６８５７
１３３８１６６４８
３．０,１０１．３２１１６９８２０７８４１１２２１４．
１,８５．２３
１３５３５２５８４６
１３０９２３．３,９７．４
㊀㊀为研究拉杆预紧力对拉杆转子动力学行为的
影响,采用１．２节的方法,从零开始逐步加大拉杆的预紧力,得到轮盘结合面法向刚度随作用力(压力)变化的曲线,见图６.将该结果代入集中质量模型进行计算,得到该型转子在不同预紧力作用下的前３阶固有频率,见图７.
图６㊀结合面法向刚度与压力的关系曲线F i g
．６㊀R e l a t i o n c u r v e o f t h e c o n t a c t s u r f a c e n o r m a l s t i f f n e s s&p r e s s
图７㊀固有频率与结合面压力的关系曲线F i g ．７㊀R e l a t i o n c u r v e o f t h e n a t u r a l f r e q u e n c y
&c o n t a c t s u r f a c e p r e s s
从表１中数据可以看出,拉杆转子与具有相
同外部形态的整体转子的固有频率差距明显,而实验实测值与用本文方法计算所得的结果非常接近,误差在５％范围内,远低于整体转子模型误差.这表明对拉杆转子开展动力学特性分析,必须考虑结合面法向刚度的影响.
由图６和图７可知,拉杆转子轮盘结合面的
法向刚度和轴向振动的固有频率主要取决于拉杆预紧力的大小,但它们之间并非简单的线性关系.特别需要注意的是,当预紧力达到一定数值(对应图６中约１００M P a 处)后,增大预紧力,结合面刚度和固有频率的变化很小,即进入所谓的预紧力
饱和阶段[１８]
.在这一阶段,拉杆转子的固有频率与整体转子趋于一致[
１１
].进一步加大拉杆的预紧力,当结合面间的压力超过２００M P a 时,无论结合面刚度还是固有频率都呈现下降的态势,其原因是压力值超出了轮盘材料的屈服极限,导致发生塑性变形的区域迅速扩大(见图３)
.此外,比较转子各阶次固有频率随结合面法向刚度变化的情况可以发现,低阶频率曲线从形态上与刚度曲线吻合程度更高,表明结合面刚度对低阶次固有频率的影响较大,而高阶次频率反映了转子整体刚度的变化.
４㊀结论
(１
)具有相同外部形态的拉杆转子和整体转子在动力学行为上存在较大的差别,分析拉杆转子的动力学特性必须考虑轮盘结合面法向刚度的影响.用本文方法得到的实验转子的轴向振动固有频率数值与实验实测值非常接近,误差低于
５％,
表明该方法切实可行,并且具有较高的精度.(２)拉杆的预紧力直接影响轮盘结合面的法向刚度和转子的固有频率,但并不一直随着拉杆预紧力的增大而增加,而是存在预紧力饱和阶段.
处于这一阶段的拉杆转子的固有频率也与整体转子趋于一致.
(３
)可借助结合面法向刚度的分析获取拉杆转子预紧力饱和的范围,此范围可以作为合理确定拉杆预紧力的数值以及保证拉杆转子结构完整性的重要设计依据.
参考文献:
[１]㊀林公舒,杨道刚．现代大功率发电用燃气轮机[M ]．
北京:机械工业出版社,２００７:１１５Ｇ１１６．
L I N G o n g s h u ,Y A N G D a o g a n g ．M o d e r n H i g h Ｇp o w e r G e n e r a t i o n G a s T u r b i n e [M ]．B e i j i n g
:M e Ｇ
７５９２ 考虑结合面法向刚度的拉杆转子轴向振动特性
何㊀谦㊀王艾伦㊀陈中祥等Copyright©博看网 . All Rights Reserved.
c h a n i c a l I n
d u s t r y P r
e s s,２００７:１１５Ｇ１１６．[２]㊀刘恒,刘意,王为民．接触界面法向刚度等效的新方法[J]．机械工程学报,２０１１,４７(１７):３７Ｇ４３．
L I U H e n g,L I U Y i,WA N G W e i m i n．N e wE q u i v aＧl e n tM e t h o df o rN o r m a lS t i f f n e s so fC o n t a c t I n t e rＧ
f a c e[J]．J o u r n a l o f M e c h a n i c a lE n
g i n e e r i n g,２０１１,
４７(１７):３７Ｇ４３．
[３]㊀高进,袁奇,李浦,等．燃气轮机拉杆转子考虑接触效应的扭转振动模态分析[J]．振动与冲击,２０１２,３１
(１２):９Ｇ１３．
G A OJ i n,Y U A N Q i,L I P u,e t a l．T o r s i o n a lV i b r aＧ
t i o n M o d a lA n a l y s i s f o r aR o dＧf a s t e n e dG a sT u r b i n e
R o t o rC o n s i d e r i n g C o n t a c tE f f e c t s[J]．J o u r n a lo f
V i b r a t i o na n dS h o c k,２０１２,３１(１２):９Ｇ１３．[４]㊀李辉光,刘恒,虞烈．考虑接触刚度的燃气轮机拉杆转子动力特性研究[J]．振动与冲击,２０１２,３１(７):４Ｇ８．
L IH u i g u a n g,L I U H e n g,Y U L i e．D y n a m i cC h a rＧ
a c t e r i s t i c s o f aR o dF a s t e n i n g R o t o r f o rG a sT u r
b i n e
C o n s i d e r i n g C o n t a c t S t i f f n e s s[J]．J o u r n a l o fV i b r aＧ
t i o na n dS h o c k,２０１２,３１(７):４Ｇ８．
[５]㊀高锐,袁奇,高进．燃气轮机拉杆转子有限元模型研究及临界转速计算[J]．热能动力学工程,２００９,２４
(３):３０５Ｇ３０８．
G A O R u i,Y U A N Q i,G A OJ i n．A S t u d y o f aF iＧ
n i t eE l e m e n tM o d e l f o r aG a sT u r b i n eT i eＧr o dR o t o r
a n d I t sC r i t i c a l S p e e dC a l c u l a t i o n[J]．J o u r n a l o fE nＧ
g i n e e r i n g f o rT h e r m a lE n e r g y a n dP o w e r,２００９,２４
(３):３０５Ｇ３０８．
[６]㊀A N T O N Y SJ,MO R E N OＧA T A N A S I O R,HA SＧS A N P O U R A．I n f l u e n c e o fC o n t a c tS t i f fN e s s e so n t h e M i c r o m e c h a n i c a lC h a r a c t e r i s t i c so fD e n s eP a rＧt i c u l a t eS y s t e m sS u b j e c t e dt oS h e a r i n g[J]．A p p l i e d
P h y s i c sL e t t e r s,２００６,８９(２１):１Ｇ３．[７]㊀周莹,程义悦．考虑接触效应的燃气轮机拉杆转子动力学模型研究[J]．装备机械,２０１４(３):５２Ｇ５７．
Z HO U Y i n g,C H E N GY i y u e．S t u d y o n t h eD y n a mＧ
i cM o d e l o fG a sT u r b i n ew i t hR o dF a s t e n i n g R o t o r
C o n s i d e r i n g t h eC o n t a c tE f f e c t[J]．T h e M a g a z i n e
o nE q u i p m e n tM a c h i n e r y,２０１４(３):５２Ｇ５７．[８]㊀巩桂亮,喻志强．M S６００１B燃气轮机发电机组轴向振动故障的分析处理[J]．华东电力,２００３,３０(５):４５Ｇ４７．
G O N G G u i l i a n g,Y U Z h i q i a n g．A x i a l V i b r a t i o n
F a u l t A n a l y s i s a n d P r o c e s s i n g o f M S６００１B
G a s
T u r b i n e[J]．E a s tC h i n a E l e c t r i cP o w e r,２００３,３０
(５):４５Ｇ４７．
[９]㊀邓哈陀JP．机械振动学[M]．北京:科学出版社,１９６１．
D E N H a r t o g JP．M e c h a n i c a l V i b r a t i o n[M]．B e iＧ
j i n g:S c i e n c eP r e s s,１９６１．
[１０]㊀唐贵基,王林,刘良玉,等．多轮盘转子系统轴向振动的固有特性研究[J]．汽轮机技术,２０１１,５３(２):５３Ｇ５５．
T A N G G u i j i,WA N G L i n,L I U L i a n g y u,e ta l．
T h eC h a r a c t e r i s t i co fL o n g i t u d i n a l V i b r a t i o n R eＧ
s e a r c ho fT u r b oＧg e n e r a t o rB a s e do nT r a n s f e r M aＧ
t r i x[J]．T u r b i n eT e c h n o l o g y,２０１１,５３(２):５３Ｇ５５．[１１]㊀王艾伦,骆舟．拉杆转子轴向振动的动力学模型[J]．中国机械工程,２００９,２０(１３):１５２４Ｇ１５２７．
WA N G A i l u n,L U O Z h o u．R e s e a r c ho nR o dF a sＧ
t e n i n g R o t o rD y n a m i c sＧa x i a lV i b r a t i o n[J]．C h i n a
M e c h a n i c a l E n g i n e e r i n g,２００９,２０(１３):１５２４Ｇ１５２７．[１２]㊀饶柱石,夏松波,汪光明．粗糙平面接触刚度的研究[J]．机械强度,１９９４,１６(２):７２Ｇ７５．
R A O Z h u s h i,X I A S o n g b o,WA N G G u a n g m i n g．
A S t u d y o fC o n t a c tS t i f f n e s so fF l a tR o u g eS u rＧ
f a c e s[J]．J o u r n a l o fM e c h a n i c a l S t r e n
g t h,１９９４,１６
(２):７２Ｇ７５．
[１３]㊀T HOMA STR．R o u g hS u r f a c e s[M]．２n d e d．L o nＧ
d o n:I m p
e r i a l C o l l e g eP r e s s,１９９９．
[１４]㊀WUJ i u n nＧJ o n g．S i m u l a t i o no fR o u g hS u r f a c e sw i t h
F F T[J]．T r i b o l o g y I n t e r n a t i o n a l,２０００,３３(１):
４７Ｇ５８．
[１５]㊀S E L L G R E N U,B J O R K L U N D S,A N D E R S S O N S．AF i n i t eE l e m e n t B a s e d M o d e l o fN o r m a l C o nＧ
t a c t b e t w e e nR o u g hS u r f a c e s[J]．W e a r,２００３,５４
(１１):１１８０Ｇ１１８８．
[１６]㊀李辉光．轴向拉杆转子结构强度及系统动力学研究
[D]．西安:西安交通大学,２０１１．
L IH u i g u a n g．S t u d y o nt h eS t r u c t u r eS t r e n g t ha n d
S y s t e m D y n a m i c so fC i r c u m f e r e n t i a lR o d F a s t e n i n g
R o t o r[D]．X i a n:X i a n J i a o t o n g U n i v e r s i t y,２０１１．[１７]㊀钟一谔,何衍宗,王正,等．转子动力学[M]．北京:清华大学出版社,１９８７．
Z HO N G Y i e,H EY a n z o n g,WA N GZ h e n g,e t a l．
R o t o rD y n a m i c s[M]．B e i j i n g:T s i n g h u aU n i v e r s i t y
P r e s s,１９８７．
[１８]㊀王为民．重型燃气轮机组合转子接触界面结构强度及系统动力学设计方法研究[D]．西安:西安交通
大学,２０１２．
W A N G W e i m i n．S t u d y o nD e s i g n M e t h o do fC o n t a c t
I n t e r f a c e S t r e n g t ha n dR o t o rＧb e a r i n g S y s t e m D y n a mＧ
i c s f o r H e a v yＧd u t y G a s T u r b i n e C o m b i n a t i o n R o t o r
[D]．X i a n:X i a n J i a o t o n g U n i v e r s i t y,２０１２．
(编辑㊀陈㊀勇)
作者简介:何㊀谦,男,１９７１年生,副教授㊁博士研究生.研究方向为机械系统动力学.发表论文１０余篇.EＧm a i l:h qＧ２００７＠q q．c o m.
８５９２
中国机械工程第２９卷第２４期２０１８年１２月下半月Copyright©博看网 . All Rights Reserved.。