干摩擦阻尼對失諧葉盤系統(tǒng)受迫振動的影響

劉雅琳,上官博,徐自力

(1.西安建筑科技大學(xué)環(huán)境與市政工程學(xué)院, 710055, 西安;2.西安熱工研究院有限公司, 710032, 西安;3.西安交通大學(xué)機(jī)械結(jié)構(gòu)強(qiáng)度與振動國家重點實驗室, 710049, 西安)

?

干摩擦阻尼對失諧葉盤系統(tǒng)受迫振動的影響

劉雅琳1,上官博2,徐自力3

(1.西安建筑科技大學(xué)環(huán)境與市政工程學(xué)院, 710055, 西安;2.西安熱工研究院有限公司, 710032, 西安;3.西安交通大學(xué)機(jī)械結(jié)構(gòu)強(qiáng)度與振動國家重點實驗室, 710049, 西安)

針對目前含非線性干摩擦阻尼結(jié)構(gòu)的失諧葉盤系統(tǒng)振動特性研究中,非線性干摩擦接觸模型的建立及高自由度非線性系統(tǒng)的求解問題,采用抗混疊時頻域融合算法和三維微滑移摩擦接觸模型,對某含有非線性干摩擦阻尼結(jié)構(gòu)的失諧葉盤系統(tǒng)的受迫振動響應(yīng)進(jìn)行計算分析,研究了干摩擦阻尼結(jié)構(gòu)對失諧葉盤系統(tǒng)受迫振動響應(yīng)的影響,并對含接觸面非線性干摩擦失諧和葉盤系統(tǒng)失諧耦合作用下的雙重失諧葉盤系統(tǒng)的振動特性進(jìn)行了研究。結(jié)果表明:非線性摩擦阻尼結(jié)構(gòu)對失諧葉盤系統(tǒng)的振動有抑制作用,在文中參數(shù)條件下振幅平均下降了9.53%;摩擦阻尼結(jié)構(gòu)對失諧葉盤中每支葉片的減振效果存在差異,每支葉片達(dá)到最佳減振效果時所對應(yīng)的最優(yōu)初始正壓力不同;接觸產(chǎn)生的摩擦失諧對失諧葉盤系統(tǒng)同樣有減振的作用;當(dāng)接觸面法向剛度失諧與葉盤系統(tǒng)的失諧強(qiáng)度和失諧量分布相同時,摩擦阻尼減振作用比摩擦力協(xié)調(diào)時減弱,在文中參數(shù)條件下振幅平均下降量減少了3.45%。

失諧葉盤系統(tǒng);干摩擦阻尼;抗混疊時頻域融合算法;三維摩擦接觸

透平機(jī)械葉片-輪盤(葉盤)系統(tǒng)的周期對稱性常常由于材質(zhì)不均勻、加工誤差、材料磨損、外物損傷等因素而遭到破壞,導(dǎo)致葉盤各子結(jié)構(gòu)之間的物理性質(zhì)存在一定差異,進(jìn)而造成葉盤的振動能量集中在少數(shù)的葉片上,產(chǎn)生葉片振動響應(yīng)局部化,加速葉盤系統(tǒng)的高周疲勞破壞。失諧問題的研究目前主要集中在對失諧葉盤振動特性的研究上,僅有少數(shù)是從抑制失諧所引起的振動局部化角度研究干摩擦阻尼結(jié)構(gòu)對失諧葉盤系統(tǒng)振動響應(yīng)的影響。實際上,開展這一方面的研究,可在葉盤系統(tǒng)的設(shè)計階段通過引入抑制失諧振動局部化的方法來減少葉盤系統(tǒng)對失諧的敏感性,從而達(dá)到減少高周疲勞破壞的目的,提高汽輪機(jī)、燃?xì)廨啓C(jī)等葉輪機(jī)械的安全可靠性。

現(xiàn)代透平機(jī)械中普遍采用拉金、凸肩、圍帶以及葉根阻尼器等結(jié)構(gòu)形式,通過人為地增加干摩擦阻尼達(dá)到減振的目的。國內(nèi)外學(xué)者曾對失諧葉盤系統(tǒng)中干摩擦阻尼結(jié)構(gòu)對振動響應(yīng)的影響進(jìn)行過研究[1-6],但是有的使用集中質(zhì)量模型對失諧葉盤進(jìn)行簡單的建模,有的對干摩擦接觸運動進(jìn)行簡化,并沒有得到較為統(tǒng)一的研究結(jié)論。此外,在實際運行的葉盤系統(tǒng)中,由于接觸干摩擦阻尼力為非線性,并且阻尼器的安裝、磨損等因素會導(dǎo)致不同的葉片干摩擦力產(chǎn)生隨機(jī)失諧(即各葉片的干摩擦力不再一致),因摩擦阻尼性質(zhì)不同而產(chǎn)生的摩擦失諧也是不可避免的,這也將與使用協(xié)調(diào)摩擦阻尼結(jié)構(gòu)進(jìn)行分析所得到的研究結(jié)果存在差異,給葉盤系統(tǒng)的振動特性研究和干摩擦阻尼器的結(jié)構(gòu)優(yōu)化設(shè)計帶來了新的困難。由于實際含干摩擦阻尼結(jié)構(gòu)的失諧葉盤系統(tǒng)非常復(fù)雜,現(xiàn)有研究對葉盤系統(tǒng)和干摩擦接觸模型進(jìn)行了大量的簡化,但是簡化得到的結(jié)果往往造成失諧葉盤系統(tǒng)許多重要動力學(xué)特性信息的丟失[7-11]。

本文將抗混疊時頻域融合算法和三維微滑移摩擦接觸模型應(yīng)用于失諧葉盤系統(tǒng)的研究中,計算分析了干摩擦阻尼結(jié)構(gòu)對失諧葉盤系統(tǒng)振動響應(yīng)的影響以及接觸面摩擦失諧情況下的雙重失諧葉盤系統(tǒng)的振動響應(yīng)。研究結(jié)論可為設(shè)計人員正確理解和設(shè)計含非線性干摩擦阻尼的葉盤系統(tǒng)、研制高性能汽輪機(jī)、燃?xì)廨啓C(jī)提供理論依據(jù)。

1 失諧葉盤系統(tǒng)模型及運動方程

不論是集中質(zhì)量參數(shù)模型還是分布質(zhì)量模型,都無法真正地模擬葉片復(fù)雜的彎扭結(jié)構(gòu),只有使用有限元方法進(jìn)行建模才能做到對復(fù)雜葉盤結(jié)構(gòu)的真實模擬。但是,由于失諧問題不同于協(xié)調(diào)問題可以采用某一個扇面模型進(jìn)行分析,必須采用整體有限元模型進(jìn)行分析,對隨機(jī)激勵和隨機(jī)失諧等問題而言,整個分析過程將非常繁復(fù)。為方便后續(xù)計算數(shù)據(jù)的比較和處理,本文在建立失諧葉盤系統(tǒng)的有限元模型時進(jìn)行了簡化,建立了一個包含18支葉片的模擬失諧葉盤系統(tǒng),如圖1所示。葉片圍帶采用平行圍帶結(jié)構(gòu),圍帶間的運動將產(chǎn)生干摩擦接觸。該模型包含72 943個節(jié)點和236 684個8節(jié)點6面體單元。材料的楊氏彈性模量Eb=206 GPa,泊松比ν=0.30,密度ρb=7 800 kg/m3。

圖1 失諧葉盤有限元模型

模型假定輪盤是周期對稱的,僅考慮葉片的失諧,并通過葉片的密度失諧來模擬葉片的質(zhì)量失諧。為不失一般性,從正態(tài)分布中隨機(jī)選取樣本作為密度失諧的偏差量,由于本文僅研究葉片質(zhì)量失諧一種形式,未涉及材質(zhì)不均勻、加工誤差、材料磨損、外物損傷等多種因素共同作用下的復(fù)合失諧,為在一定程度上體現(xiàn)出多種失諧因素共同作用對葉盤產(chǎn)生的影響,將失諧量正態(tài)分布的平均值設(shè)為0,標(biāo)準(zhǔn)差設(shè)為10%,每支葉片失諧量的具體數(shù)值詳見表1。在進(jìn)行模態(tài)分析時,將葉片和輪盤之間由于旋轉(zhuǎn)離心力產(chǎn)生的接觸應(yīng)力作為預(yù)應(yīng)力,即通過引入預(yù)應(yīng)力來考慮離心力載荷對葉盤振動的影響,并通過在葉根與輪緣處添加約束方程的形式考慮葉片和輪盤的耦合作用。

表1 失諧強(qiáng)度σρ=10%時各葉片的失諧量分布

葉盤系統(tǒng)的整體運動方程可表示為

(1)

式中:u(t)是位移向量;M、K、C分別表示葉盤系統(tǒng)的質(zhì)量矩陣、剛度矩陣和材料阻尼矩陣;fl(t)表示周期性激振力向量;fnl(t,u(t))表示作用在圍帶處的非線性干摩擦力向量。

2 系統(tǒng)求解方法及摩擦接觸模型

2.1 抗混疊時頻域融合算法

采用抗混疊時頻域融合算法[12]對式(1)進(jìn)行求解時,需同時對方程兩邊進(jìn)行傅里葉變換,并引入ωk以滿足離散傅里葉變換的需要。使用U(ω)表示位移向量u(t)的傅里葉變換所對應(yīng)的頻域形式,包含有Nk個諧波分量,Fl(ω)表示激振力fl(t)的傅里葉變換所對應(yīng)的頻域形式,Fnl(ω,U(ω))表示非線性摩擦力fnl(t,u(t))的傅里葉變換所對應(yīng)的頻域形式。這樣,得到如下非線性代數(shù)方程組

H(ω)U(ω)=Fl(ω)+Fnl(ω,U(ω))

(2)

H(ω)=-(kω)2M+jkωC+K

(3)

(4)

式中:H(ω)表示系統(tǒng)的動剛度矩陣,與頻率ω相關(guān);Δt和Nk分別表示離散傅里葉變換中所考慮的采樣時間及采樣點數(shù)。

圖2 抗混疊時頻域融合算法求解示意圖

2.2 三維微滑移摩擦接觸模型

在周期激勵作用下,干摩擦阻尼結(jié)構(gòu)摩擦接觸界面的運動形式往往很復(fù)雜,導(dǎo)致干摩擦接觸解析模型無法描述其界面約束力。當(dāng)葉片振動時,相鄰葉片的圍帶接觸面會發(fā)生三維的相對運動,沿界面切向的相對運動會產(chǎn)生切向摩擦力,消耗葉片振動能量,沿界面法向的相對運動會產(chǎn)生法向接觸正壓力會影響切向摩擦力的大小和分布,進(jìn)而影響圍帶阻尼結(jié)構(gòu)的減振特性。求解干摩擦阻尼結(jié)構(gòu)葉片系統(tǒng)非線性振動響應(yīng)的難點在于如何準(zhǔn)確地確定接觸界面摩擦約束力,這是因為界面摩擦約束力是葉片系統(tǒng)振動響應(yīng)的非線性函數(shù),它既取決于葉片系統(tǒng)的振動響應(yīng),又會影響葉片系統(tǒng)的振動響應(yīng)。

為了提高求解的準(zhǔn)確性,需要考慮圍帶接觸面的黏滯-滑移共存狀態(tài),以及接觸面法向正壓力分布不均勻的情況。本文在進(jìn)行接觸面運動的描述時,將圍帶接觸面離散成多個微小接觸區(qū)域,結(jié)合微滑移摩擦模型在接觸面之間建立多個接觸點對,如圖3所示。每個接觸點對都能夠描述三維接觸運動,且每個接觸點對的接觸運動狀態(tài)是單獨判斷的,整個圍帶接觸面的作用力為各個接觸點對的合力。

圖3 兩接觸面之間的三維微滑移摩擦接觸模型

以兩接觸面間任意一個接觸點對為例來說明求解摩擦約束力的方法。假設(shè)該摩擦接觸點對包含兩個摩擦節(jié)點A和B,以及用來模擬切向彈性接觸和法向彈性接觸的無質(zhì)量彈簧Kt和kn,如圖4所示。在面2上建立局部坐標(biāo)系x′y′z′,將x′軸和y′軸與面2的切向相重合,z′軸與面2的法向重合。

圖4 單個離散點上接觸對之間的相對運動

由于存在加工誤差、材料非均勻、非正常工況運行磨損等因素,將切向接觸剛度定義如下以表征上述原因所引起的各向異性

(5)

假設(shè)節(jié)點A始終與面1保持黏滯,節(jié)點B將在兩接觸面發(fā)生相對運動時沿著面2做黏滯-滑移運動,產(chǎn)生切向摩擦約束力。當(dāng)兩接觸面沒有振動時,節(jié)點A、B最初重合在一起,在總體坐標(biāo)系統(tǒng)中具有同樣的坐標(biāo)。當(dāng)兩接觸面發(fā)生相對運動時,由于假設(shè)節(jié)點A與面1始終保持黏滯,所以通過面1的運動軌跡就可得到節(jié)點A的運動軌跡,只需對節(jié)點B的運動狀態(tài)進(jìn)行計算,就可以得到接觸面之間的運動軌跡和摩擦約束力的分布。當(dāng)面1相對于面2的法向運動分量過大而導(dǎo)致點B跟面2發(fā)生分離時,此時的摩擦約束力用零表示。通過分析接觸點對之間的運動關(guān)系,不難發(fā)現(xiàn)它具有如下特點:接觸點對之間通過彈簧連接,在運動過程中,接觸點對之間的距離保持在某一變化范圍之內(nèi),無法直接寫出運動軌跡;此外,對接觸問題而言,求解滑動點的運動軌跡是為了最終求解摩擦約束力,而摩擦約束力的方向取決于兩接觸點之間的相對位置關(guān)系。摩擦力的具體求解方法詳見文獻(xiàn)[13]。

3 干摩擦對失諧葉盤振動響應(yīng)的影響

3.1 失諧葉盤非線性振動響應(yīng)

在協(xié)調(diào)葉盤系統(tǒng)中干摩擦阻尼結(jié)構(gòu)常被用來消耗葉片的振動能量,以達(dá)到抑制共振響應(yīng)幅度的目的。那么,在失諧葉盤系統(tǒng)中,干摩擦阻尼結(jié)構(gòu)還能否達(dá)到減振的作用,其作用效果將在本小節(jié)中進(jìn)行研究。文獻(xiàn)[14]曾在研究摩擦件對失諧葉片-輪盤系統(tǒng)的共振響應(yīng)的影響時指出,由于摩擦接觸的非線性特征,使得對具有摩擦阻尼結(jié)構(gòu)的葉片-輪盤系統(tǒng)振動響應(yīng)的分析變得較為困難。

3.1.1 參數(shù)設(shè)定 氣流通過進(jìn)氣口時因受非旋轉(zhuǎn)部件與旋轉(zhuǎn)部件的擾動,會在葉片上形成周期性的激振力載荷,通過傅里葉諧波分析可將周期性的激振力載荷表示為若干諧波分量的疊加。由于每個機(jī)組、每一級葉盤在實際運行過程中受到的激振力都不相同,因此本文通過在每支葉片頂部沿3個方向施加正弦激勵對葉片所受激振力進(jìn)行模擬。葉盤中第j支葉片所受的激勵力為

(6)

當(dāng)質(zhì)量失諧強(qiáng)度σρ=10%時,葉盤上18支葉片在一階固有頻率附近的幅頻響應(yīng)曲線如圖5所示,圖中的振動響應(yīng)是3個方向的合成響應(yīng)。

圖5 σρ=10%時的一階頻響曲線

從圖5可以看出,由于干摩擦阻尼結(jié)構(gòu)產(chǎn)生的摩擦約束作用,18支葉片的非線性響應(yīng)振幅與線性情況相比都產(chǎn)生了一定的變化。在本文參數(shù)條件下,葉片非線性響應(yīng)的最大共振振幅為0.35mm,最小共振振幅為0.21mm,比葉片線性響應(yīng)的最大共振振幅0.36mm和最小共振振幅0.27mm都有所下降。絕大多數(shù)非線性響應(yīng)的振幅都比線性情況有所下降,共振頻率增大。失諧葉盤上18支葉片的非線性響應(yīng)共振振幅比線性響應(yīng)的振幅最大降幅為21.50%,平均降幅為9.53%,減振效果顯著。總體而言,干摩擦阻尼結(jié)構(gòu)在失諧葉盤系統(tǒng)中仍能起到抑制共振幅值的作用。

(a)葉片1

(b)葉片2

(c)葉片3圖6 σρ=10%時不同初始正壓力下的頻響曲線

3.1.2 接觸面初始正壓力對振動響應(yīng)的影響 當(dāng)σρ=10%時,不同初始正壓力條件下葉盤上某3支葉片在一階固有頻率附近的幅頻響應(yīng)曲線如圖6所示,圖中的振動響應(yīng)是3個方向的合成響應(yīng)。葉片1、2、3分別對應(yīng)葉盤中共振響應(yīng)幅值最大、共振響應(yīng)幅值最小以及考慮非線性摩擦力之后與線性情況相比振幅降幅最大的葉片。

隨著接觸面初始正壓力的增大,3支葉片的共振頻率都在一直增加,而共振幅值的變化規(guī)律卻不盡相同。葉片2的共振幅值一直減小,葉片1和葉片3的共振幅值先減小后增大,都存在一個最優(yōu)的初始正壓力使其振動響應(yīng)最小,但二者所對應(yīng)的最優(yōu)初始正壓力并不相同:圖6a中的最優(yōu)初始正壓力n0=5N,圖6c中的最優(yōu)初始正壓力n0=4N。這就意味著,對失諧葉盤系統(tǒng)進(jìn)行優(yōu)化設(shè)計時,需要考慮每支葉片的振動狀況,完全統(tǒng)一的參數(shù)設(shè)計并不能使運行效果達(dá)到最佳。

3.2 摩擦失諧對失諧葉盤系統(tǒng)振動響應(yīng)的影響

在實際葉盤系統(tǒng)中,干摩擦阻尼結(jié)構(gòu)的性質(zhì)不會完全一致,也會出現(xiàn)與葉盤質(zhì)量失諧和剛度失諧形式類似的非線性摩擦阻尼失諧形式。文獻(xiàn)[2]對采用集中質(zhì)量模型建模的協(xié)調(diào)葉盤系統(tǒng)進(jìn)行了阻尼失諧條件下的振動響應(yīng)研究,結(jié)果顯示,阻尼失諧可導(dǎo)致與葉片剛度失諧程度相近似的振動幅值。文獻(xiàn)[8]結(jié)果顯示,接觸面摩擦失諧可顯著改變?nèi)~盤系統(tǒng)的受迫振動響應(yīng)特性。上述研究都是基于協(xié)調(diào)葉盤系統(tǒng)對摩擦阻尼失諧進(jìn)行的研究。實際上,由于失諧葉盤的振動局部化作用,葉片之間的摩擦接觸狀態(tài)并不相同,不可避免地產(chǎn)生摩擦阻尼結(jié)構(gòu)失諧現(xiàn)象。對于包含干摩擦阻尼結(jié)構(gòu)的葉片系統(tǒng)而言,葉盤失諧和摩擦阻尼結(jié)構(gòu)失諧應(yīng)該是同時存在的。無論是采用協(xié)調(diào)摩擦阻尼結(jié)構(gòu)還是采用摩擦組尼失諧對協(xié)調(diào)葉盤系統(tǒng)進(jìn)行減振研究都與實際情況存在差距。

3.2.1 接觸面摩擦失諧的形式 本節(jié)在考慮接觸面摩擦失諧時,選取接觸面法向正壓力為失諧項,分別研究了初始正壓力n0和法向接觸剛度kn失諧情況下,失諧葉盤系統(tǒng)的振動響應(yīng)。接觸面摩擦失諧情況下法向摩擦力的表達(dá)式為

k=1,2,…,Nk

(7)

式中:n0(1+σni)表示接觸面初始正壓力;σni為初始正壓力的失諧量;kn(1+σki)表示接觸面法向接觸剛度;σki為法向接觸剛度的失諧量;Δz′表示局部坐標(biāo)系下的法向位移分量。當(dāng)摩擦力協(xié)調(diào)時,σni=σki=0,接觸面初始正壓力n0(1+σni)取值相同,接觸面法向剛度kn(1+σki)取值亦相同。當(dāng)接觸面初始正壓力失諧時,σki=0,σni的值取隨機(jī)正態(tài)分布,失諧強(qiáng)度用失諧量的標(biāo)準(zhǔn)差表示。當(dāng)接觸面法向接觸剛度失諧時,σni=0,σki的值取隨機(jī)正態(tài)分布,失諧強(qiáng)度用失諧量的標(biāo)準(zhǔn)差表示。

3.2.2 雙重失諧情況下葉盤系統(tǒng)的振動特性 以σρ=10%的葉盤系統(tǒng)為例,對干摩擦阻尼結(jié)構(gòu)中含有初始正壓力失諧的雙重失諧系統(tǒng)進(jìn)行振動響應(yīng)分析。計算中,干摩擦阻尼結(jié)構(gòu)初始正壓力和接觸面法向剛度的失諧量均假設(shè)為10%,在初始正壓力和葉片質(zhì)量雙重失諧的情況下,葉盤中每支葉片的最大振幅與摩擦力協(xié)調(diào)時相比的變化情況如圖7所示。

圖7 初始正壓力和葉片質(zhì)量雙重失諧時的最大振幅

在考慮摩擦阻尼的前提下,接觸面初始正壓力失諧情況下葉片的最大振幅與摩擦力協(xié)調(diào)情況下葉片的最大振幅相比,有的增加,有的減少,還有的維持不變,沒有統(tǒng)一的變化規(guī)律。在本文參數(shù)條件下,接觸面初始正壓力失諧時,18支葉片的非線性響應(yīng)共振振幅與線性響應(yīng)的振幅相比,最大降幅為21.85%,平均降幅為9.20%。由此可知,接觸面初始正壓力失諧對失諧葉盤系統(tǒng)同樣有減振的作用,且與摩擦力協(xié)調(diào)情況的減振效果相當(dāng)。

圖8 法向接觸剛度和葉片質(zhì)量雙重失諧時的最大振幅

接觸面法向接觸剛度失諧對失諧葉盤的影響,見圖8。在都包含摩擦阻尼結(jié)構(gòu)的前提下,接觸面法向接觸剛度失諧情況下葉片最大振幅與摩擦力協(xié)調(diào)情況下葉片最大振幅相比,僅有1支葉片的振幅減小,其他17支葉片的振幅都有所增加。在本文參數(shù)條件下,接觸面法向接觸剛度失諧時,18支葉片的非線性響應(yīng)共振振幅與線性響應(yīng)的振幅相比,最大降幅為20.58%,平均降幅為6.08%。圖7、圖8相比,接觸面法向接觸剛度失諧時振幅變化更明顯,但是減振效果有所下降,與摩擦力協(xié)調(diào)時相比,振幅下降量減少了3.45%。

4 結(jié) 論

本文建立了失諧葉盤的有限元分析模型,采用葉片密度的正態(tài)隨機(jī)分布模擬葉盤的質(zhì)量失諧。基于抗混疊時頻域融合算法和三維微滑移摩擦接觸模型,研究了干摩擦阻尼結(jié)構(gòu)對失諧葉盤系統(tǒng)振動響應(yīng)的影響以及接觸面干摩擦失諧情況下失諧葉盤系統(tǒng)的振動特性,主要結(jié)論如下。

(1)考慮干摩擦阻尼結(jié)構(gòu)所產(chǎn)生的摩擦約束力后,葉片的共振幅值明顯降低,共振頻率增加。相鄰葉片間的摩擦約束力不同導(dǎo)致干摩擦阻尼結(jié)構(gòu)對每支葉片的減振效果存在差異,每支葉片達(dá)到最佳減振效果時所對應(yīng)的最優(yōu)初始正壓力不同,因此失諧葉盤的優(yōu)化設(shè)計需要考慮每支葉片的振動狀況,存在很多難點。對應(yīng)于相同的接觸面初始正壓力,阻尼結(jié)構(gòu)的減振效果對有的葉片比較好,對有的葉片比較差。從干摩擦減振控制的角度來講,完全統(tǒng)一的參數(shù)設(shè)計并不能使運行效果達(dá)到最佳,失諧葉盤系統(tǒng)的減振設(shè)計需要考慮摩擦控制參數(shù)與葉片失諧程度相互匹配的問題。

(2)接觸面初始正壓力失諧和接觸面法向剛度失諧所產(chǎn)生的摩擦失諧對失諧葉盤系統(tǒng)都有減振的作用。當(dāng)接觸面法向接觸剛度失諧與葉盤系統(tǒng)的失諧強(qiáng)度相同且失諧量分布相對應(yīng)時,法向接觸剛度失諧時干摩擦阻尼的減振效果與摩擦力協(xié)調(diào)時相比被減弱,在本文參數(shù)條件下,振幅平均下降量減少了3.45%。

[1] 白斌, 白廣忱, 童曉晨, 等. 整體葉盤結(jié)構(gòu)失諧振動的國內(nèi)外研究狀況 [J]. 航空動力學(xué)報, 2014, 29(1): 91-103. BAI Bin, BAI Guangchen, TONG Xiaochen, et al. Research on vibration problem of integral mistuned bladed disk assemblies at home and abroad [J]. Journal of Aerospace Power, 2014, 29(1): 91-103.

[2] LIN C C, MIGNOLET M P. An adaptive perturbation scheme for the analysis of mistuned bladed disks [J]. ASME Journal of Engineering for Gas Turbines and Power, 1997, 119: 153-160.

[3] GRIFFIN J H, SINHA A. The interaction between mistuning and friction in the forced response of bladed disk assemblies [J]. ASME Journal of Engineering for Gas Turbines and Power, 1985, 107: 205-211.

[4] PETROV E P. A method for forced response analysis of mistuned bladed disks with aerodynamic effects included [J]. ASME Journal of Engineering for Gas Turbines and Power, 2010, 132(6): 062502.

[5] 王紅建. 復(fù)雜耦合失諧葉片-輪盤系統(tǒng)振動局部化問題研究 [D]. 西安: 西北工業(yè)大學(xué), 2006.

[6] CHA D, SINHA A. Statistics of response of a mistuned and frictionally damped bladed disk assembly subjected to white noise and narrow band excitations [J]. Probabilistic Engineering Mechanics, 2006, 21: 384-396.

[7] PETROV E P, EWINS D J. Advanced modeling of under platform friction dampers for analysis of bladed disk vibration [J]. ASME Journal of Turbomachinery, 2007, 129: 143-150.

[8] PETROV E P, EWINS D J. Method for analysis of nonlinear multi harmonic vibrations of mistuned bladed disks with scatter of contact interface characteristics [J]. ASME Journal of Turbomachinary, 2005, 127: 128-136.

[9] JARELAND M H, CSABA G. Friction damper mistuning of a bladed disk and optimization with respect to work [C]∥ASME Turbo Expo 2000: Power for Sea, and Air. New York, USA: ASME, 2000: V004T03A009.

[10]王紅建, 賀爾銘, 余仕俠. 具有干摩擦散亂失諧的葉片-輪盤受迫響應(yīng)特性 [J]. 航空動力學(xué)報, 2006, 21(4): 711-715. WANG Hongjian, HE Erming, YU Shixia. Forced response characteristics of bladed disks with disordered dry friction [J]. Journal of Aerospace Power, 2006, 21(4): 711-715.

[11]王艾倫, 龍清. 具有非線性摩擦阻尼隨機(jī)失諧的葉盤系統(tǒng)響應(yīng)特性 [J]. 航空動力學(xué)報, 2011, 26(1): 178-184. WANG Ailun, LONG Qing. Forced response characteristics of bladed disks with random mistuned non-linear friction damping [J]. Journal of Aerospace Power, 2011, 26(1): 178-184.

[12]劉雅琳, 上官博, 徐自力. 抗混疊時頻域融合算法及在葉片響應(yīng)分析中的應(yīng)用 [J]. 航空動力學(xué)報, 2012, 27(6): 1238-1242. LIU Yalin, SHANGGUAN Bo, XU Zili. Anti-aliasing alternating frequency/time domain method and its application in response analysis of blade systems [J]. Journal of Aerospace Power, 2012, 27(6): 1238-1242.

[13]谷偉偉, 徐自力. 干摩擦阻尼葉片的界面約束力描述及振動響應(yīng)求解 [J]. 振動工程學(xué)報, 2012, 25(1): 64-67. GU Weiwei, XU Zili. Description of constraint force and prediction of vibration responses for dry friction damped blade [J]. Journal of Vibration Engineering, 2012, 25(1): 64-67.

[14]FEINER M D, GRIFFIN J H. Mistuning identification of bladed disks using a fundamental mistuning model: part 1 theory [J]. ASME Journal of Turbomachinary, 2004, 126(1): 150-158.

[本刊相關(guān)文獻(xiàn)鏈接]

孟慶虎,孟慶豐,朱永生,等.用于機(jī)械系統(tǒng)固有頻率及阻尼比計算的改進(jìn)頻域方法.2015,49(8):1-5.[doi:10.7652/xjtuxb201508001]

李應(yīng)剛,陳天寧,王小鵬,等.外部動態(tài)激勵作用下齒輪系統(tǒng)非線性動力學(xué)特性.2014,48(1):101-105.[doi:10.7652/xjtuxb201401017]

劉雅琳,徐自力,上官博,等.時頻域交互法在微滑移干摩擦阻尼葉片振動分析中的應(yīng)用.2009,43(11):22-26.[doi:10.7652/xjtuxb200911005]

徐自力,上官博,吳其林,等.松裝葉片燕尾型葉根干摩擦力模型及振動響應(yīng)分析.2009,43(7):1-5.[doi:10.7652/xjtuxb200907001]

尤晉閩,陳天寧.懸臂梁平面結(jié)合面參數(shù)的識別技術(shù)研究.2008,42(11):1323-1326.[doi:10.7652/xjtuxb200811002]

谷偉偉,徐自力,胡哺松,等.微滑移阻尼葉片最優(yōu)正壓力的影響因素分析.2008,42(7):828-832.[doi:10.7652/xjtuxb 200807010]

徐自力,谷偉偉,呂強(qiáng),等.基于微動滑移模型的葉片阻尼結(jié)構(gòu)參數(shù)研究.2007,41(5):507-511.[doi:10.7652/xjtuxb2007 05001]

(編輯 杜秀杰)

Effects of Dry Friction Damping on Forced Vibration Response of Mistuned Bladed Disk

LIU Yalin1,SHANGGUAN Bo2,XU Zili3

(1. School of Environmental and Municipal Engineering, Xi’an University of Architecture and Technology,Xi’an 710055, China; 2. Xi’an Thermal Power Research Institute Co., Ltd., Xi’an 710032, China;3. State Key Laboratory for Strength and Vibration of Mechanical Structures, Xi’an Jiaotong University, Xi’an 710049, China)

In the existing vibration characteristics research of mistuned bladed disk system with nonlinear dry friction damping structures, we have to encounter establishment of nonlinear friction contact model and solution of high degree of freedom nonlinear systems. Here the anti-aliasing hybrid frequency-time domain method and three-dimensional microslip friction contact model are chosen to evaluate the forced vibration response of a mistuned bladed disk with dry friction damping structure. The friction damping effect on forced vibration response of mistuned bladed disk, and the vibration characteristics of mistuned bladed disk under the condition of friction mistuned are investigated. It is found that the nonlinear friction vibration damping structure enables to inhibit vibration of mistuned bladed disk system with a 9.53% average amplitude decrease. The friction damping effect on each blade is different, and the corresponding optimal initial normal load of each blade reaching the optimal damping effect is not identical too. The normal contact stiffness mistuning and initial normal load mistuning also have the damping effect on the mistuned bladed disk. When the mistuned strength and mistuned distribution of the normal contact stiffness mistuning are the same as those of mistuned bladed disk, the friction damping effect attenuates in comparison with the situation of friction tuned, and the average decline of vibration amplitude reduces by 3.45%.

mistuned blade disc; dry friction damping; antialiasing hybrid frequency-time domain method; three-dimensional friction contact

2015-03-08。

劉雅琳(1983—),女,講師;徐自力(通信作者),男,教授,博士生導(dǎo)師。 基金項目:國家自然科學(xué)基金資助項目(51505362,51275385,51405384)。

時間:2015-11-11

10.7652/xjtuxb201602019

TK263.3

A

0253-987X(2016)02-0111-07

網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20151111.1181.006.html