張鵬飛,王志恒,席光

(西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安)

?

用于柔性轉(zhuǎn)子主動(dòng)控制的等幾何Timoshenko梁模型及其數(shù)值驗(yàn)證

張鵬飛,王志恒,席光

(西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安)

為了克服采用標(biāo)準(zhǔn)有限元方法建立的轉(zhuǎn)子模型自由度多不便直接用于控制器設(shè)計(jì)的問(wèn)題,結(jié)合等幾何分析方法的精度高、自由度少的特點(diǎn),提出了把低階等幾何Timoshenko梁模型運(yùn)用到柔性轉(zhuǎn)子的主動(dòng)控制中,并進(jìn)行了數(shù)值驗(yàn)證。首先,給出了半離散等幾何模型及在主動(dòng)控制中作為轉(zhuǎn)子控制輸入信號(hào)的各類邊界條件;其次,分別對(duì)比了高、低階等幾何模型的奇異值響應(yīng)以及簡(jiǎn)支條件下的高、低階模型的數(shù)值解與理論解;最后,在采用磁懸浮軸承支撐和給定分布不平衡力擾動(dòng)的條件下,對(duì)轉(zhuǎn)子進(jìn)行了分散比例微分(PD)仿真控制。結(jié)果表明:所采用的低階模型的奇異值響應(yīng)在前6階臨界轉(zhuǎn)速范圍內(nèi)與高階模型基本一致;高階模型的前10階模態(tài)頻率很好地吻合了理論解,低階模型前4階模態(tài)頻率誤差在0.2%以內(nèi);高、低階等幾何梁模型下的轉(zhuǎn)子不平衡振動(dòng)位移穩(wěn)態(tài)響應(yīng)的差別很小,該誤差可看成工作轉(zhuǎn)速下的同頻小擾動(dòng)。低階等幾何梁模型在低頻范圍的高精度驗(yàn)證了該方法所得低階模型直接用于控制器設(shè)計(jì)的可行性。

等幾何方法;柔性轉(zhuǎn)子模型;Timoshenko梁;振動(dòng)控制;磁懸浮軸承

在航空、電力、石化等領(lǐng)域,各類旋轉(zhuǎn)機(jī)械都向著高速、高功重比方向發(fā)展,為了限制轉(zhuǎn)子的重量和線速度,減小離心力,高速轉(zhuǎn)子往往采用細(xì)長(zhǎng)型的柔性轉(zhuǎn)子,工作轉(zhuǎn)速通常高于其一階、二階臨界轉(zhuǎn)速,甚至達(dá)到三、四階臨界轉(zhuǎn)速以上[1]。主動(dòng)磁懸浮軸承的高速、無(wú)磨損、無(wú)潤(rùn)滑及可控性等優(yōu)點(diǎn)使其成為高速轉(zhuǎn)子支承的理想方案[2]。主動(dòng)磁懸浮軸承支撐的轉(zhuǎn)子控制方法主要基于數(shù)學(xué)模型,雖然許多先進(jìn)算法(如魯棒控制、滑模控制等)可以降低對(duì)轉(zhuǎn)子建模精度的敏感性,但轉(zhuǎn)子的建模精度還是會(huì)直接影響其控制性能。

大量實(shí)踐表明,5自由度剛性轉(zhuǎn)子簡(jiǎn)單、有效,可直接用于主動(dòng)磁懸浮軸承的剛性轉(zhuǎn)子控制[3-5]。對(duì)于柔性轉(zhuǎn)子,其主動(dòng)控制要求建立更精確的轉(zhuǎn)子模型,采用標(biāo)準(zhǔn)有限元方法對(duì)柔性轉(zhuǎn)子進(jìn)行建模和分析已經(jīng)比較成熟,且有限元方法相比傳遞矩陣法更加精確[6]。但是,有限元方法所得模型自由度數(shù)目較大,不便于直接用于控制器的設(shè)計(jì),尤其當(dāng)控制器設(shè)計(jì)的復(fù)雜程度依賴于轉(zhuǎn)子的數(shù)學(xué)模型時(shí),例如魯棒控制。用于這類控制器的柔性轉(zhuǎn)子模型通常先采用標(biāo)準(zhǔn)有限元方法或集總參數(shù)法建立高階轉(zhuǎn)子模型,模型經(jīng)過(guò)修正后,利用子結(jié)構(gòu)技術(shù)(例如Guyan、Hurty模態(tài)截?cái)嗟确椒?或狀態(tài)空間技術(shù)對(duì)高階轉(zhuǎn)子進(jìn)行簡(jiǎn)化[7]。

標(biāo)準(zhǔn)有限元方法、子結(jié)構(gòu)技術(shù)和狀態(tài)空間方法已經(jīng)較為成熟,但由于有限元方法求解自由振動(dòng)問(wèn)題所得頻率譜的不連續(xù)性[12],基于有限元法的子結(jié)構(gòu)技術(shù)要求其所建模型具有較多的自由度數(shù)目,否則難以保證模型的精度。狀態(tài)空間方法仍然缺乏處理大規(guī)模系統(tǒng)的能力,需要采用子結(jié)構(gòu)技術(shù)先進(jìn)行一次降階[2]。

等幾何方法[8]是一種基于非均有理B樣條基函數(shù)的數(shù)值方法,主要實(shí)現(xiàn)了幾何模型與分析模型的一致表達(dá),已經(jīng)廣泛應(yīng)用于工程問(wèn)題的優(yōu)化,結(jié)構(gòu)的強(qiáng)度、振動(dòng)分析[9],流體分析[10-11]和電磁分析等領(lǐng)域。在轉(zhuǎn)子動(dòng)力學(xué)領(lǐng)域,文獻(xiàn)[12]研究了歐拉-伯努利梁的自由振動(dòng),與有限元方法對(duì)比,說(shuō)明等幾何方法精度高,且其頻率譜連續(xù);文獻(xiàn)[13]分析了Timoshenko梁的自由振動(dòng),成功避免了振動(dòng)分析中剪切自鎖現(xiàn)象;文獻(xiàn)[14]分析了各類邊界條件下方形Timoshenko梁的自由振動(dòng),發(fā)現(xiàn)網(wǎng)格k細(xì)分方法得到的結(jié)果最好,即在較少單元數(shù)目和基函數(shù)階次下可得到更高精度;文獻(xiàn)[15]提出了一種超收斂等幾何方法,并用于歐拉梁的振動(dòng)問(wèn)題;文獻(xiàn)[16]采用等幾何方法建立非線性幾何歐拉梁模型。

由于等幾何方法具有得到精度高且階數(shù)低的梁模型特性,本文提出采用等幾何方法獲得低階Timoshenko轉(zhuǎn)子模型以替代由有限元模型并進(jìn)行模型縮減后所得模型,由此可以把低階等幾何Timoshenko轉(zhuǎn)子模型直接用于柔性轉(zhuǎn)子的控制器設(shè)計(jì)。所得模型可以根據(jù)需要實(shí)時(shí)調(diào)整物理參數(shù),特別是像轉(zhuǎn)子轉(zhuǎn)速等實(shí)時(shí)變化的物理參數(shù)。數(shù)值分析中的邊界條件對(duì)應(yīng)于控制領(lǐng)域中的輸入信號(hào),為了能使其在主動(dòng)控制中運(yùn)用,本文將各類主動(dòng)控制中可能用到的邊界條件進(jìn)行了整理和推導(dǎo),并將數(shù)值解與理論解進(jìn)行對(duì)比,驗(yàn)證所建模型的正確性。最后,進(jìn)行了在高、低階等幾何Timoshenko梁模型下的奇異值響應(yīng)分析,并采用磁懸浮軸承支撐轉(zhuǎn)子,考慮分布質(zhì)量不平衡,分別進(jìn)行了分散PD控制仿真,通過(guò)對(duì)比仿真結(jié)果,初步驗(yàn)證該模型用于主動(dòng)控制的可行性。

1 等幾何Timoshenko梁建模

1.1 微分方程

假設(shè)梁微小變形且不考慮扭轉(zhuǎn)變形、軸向變形的前提下,對(duì)于變截面均質(zhì)對(duì)稱圓形截面Timoshenko梁-軸,不失一般性,令其中性軸位于z軸,且z∈(0,L),控制方程的微分形式[16]如下

(1)

式中:ux(z,t)、uy(z,t)為中性軸橫向位移分布;φx(z,t)、φy(z,t)為截面轉(zhuǎn)角分布;fx(z,t)、fy(z,t)為分布載荷力;mx(z,t)、my(z,t)為分布載荷力矩;A(z)為橫斷截面面積;Ix(z)、Iy(z)為截面慣性矩;Ip(z)為極慣性矩;ρ為材料的密度;E為材料的彈性模量;G為材料的剪切模量;μ為梁的剪切修正系數(shù);Ω為轉(zhuǎn)子工作轉(zhuǎn)速;L為轉(zhuǎn)子長(zhǎng)度。

1.2 半離散等幾何Timoshenko梁模型

與標(biāo)準(zhǔn)有限元相比,等幾何單元采用非均有理B樣條基函數(shù)作為插值基函數(shù),本文采用B樣條作為等幾何單元的插值基函數(shù)。式(1)的變分形式可直接采用伽遼金(Galerkin)方法推導(dǎo)得到。本文在采用等幾何方法時(shí),僅對(duì)空間進(jìn)行離散,保持時(shí)間域連續(xù),可得到半離散動(dòng)力學(xué)模型。由于篇幅所限,關(guān)于等幾何方法的詳細(xì)介紹可以參考文獻(xiàn)[9],這里僅列出最后的半離散模型

(2)

式中:M、G、K分別為4n×4n總質(zhì)量矩陣、總陀螺矩陣和總剛度矩陣;F為4n×1廣義力向量矩陣。

(3)

式中:W=diag(n(ξ),n(ξ),n(ξ),n(ξ)),n(ξ)為p次B樣條的n個(gè)基函數(shù)組成的假設(shè)模態(tài)列向量,可表示為

其中B樣條基函數(shù)Ni,p(ξ)的定義方法可參考文獻(xiàn)[18]。本文采用k型網(wǎng)格細(xì)分方法[14],等幾何單元數(shù)目ne、基函數(shù)數(shù)目n以及基函數(shù)階次p滿足以下關(guān)系式n=ne+p。

2 邊界條件處理

為了便于將等幾何梁模型用于柔性轉(zhuǎn)子的主動(dòng)振動(dòng)控制,本文整理和推導(dǎo)了在振動(dòng)控制、分析方面可能用到的邊界條件。邊界條件主要考慮了齊次邊界條件,包括簡(jiǎn)支、固持、導(dǎo)向、自由、拉伸彈簧、抗彎彈簧等6種邊界條件,還有轉(zhuǎn)子上的集中廣義力、集中廣義質(zhì)量單元、彈簧阻尼單元以及不平衡分布質(zhì)量引起的不平衡力非齊次邊界條件。

2.1 齊次邊界條件

齊次邊界條件是幾何、力邊界條件的組合,如表1所示。位移和轉(zhuǎn)角邊界條件對(duì)應(yīng)于有限元中的第一類邊界條件(Dirichlet邊界條件),力矩邊界條件對(duì)應(yīng)第二類邊界條件(Neumann邊界條件),力邊界條件為第三類邊界條件(Robin邊界條件),它們的處理方法與有限元方法相同。在表1中,k為拉伸彈簧剛度,κ為抗彎彈簧剛度,zp為邊界條件位置,特別地,在轉(zhuǎn)子兩端,zp=0,L。

表1 齊次邊界條件

2.2 集中的廣義力、質(zhì)量及彈簧阻尼單元

集中廣義力是轉(zhuǎn)子振動(dòng)控制的輸入,首先給出作用梁上在任意位置下徑向廣義力表達(dá)式。根據(jù)集中廣義力的位置,可以分為3種情形:①轉(zhuǎn)子兩端點(diǎn)處;②等幾何單元兩端點(diǎn)處;③等幾何單元兩端點(diǎn)間。經(jīng)過(guò)推導(dǎo),等幾何坐標(biāo)下任意位置的廣義力F都可以表示成

(4)

(5)

梁上的集中質(zhì)量、慣量可以做為外力引入到力學(xué)模型中去,結(jié)合式(5)有

(6)

式中

(7)

將式(6)、式(7)代入式(5),分別有

(8)

(9)

2.3 分布不平衡力

分布不平衡力由不平衡質(zhì)量引起。假設(shè)質(zhì)量偏心的幅值沿軸向分布為e(z),初始相位為0,則轉(zhuǎn)速恒定條件下分布不平衡力為

(10)

根據(jù)式(4),對(duì)式(10)積分,則有

(11)

則由分布不平衡質(zhì)量引起的不平衡力矩陣為

3 算例驗(yàn)證

以下從兩個(gè)方面來(lái)驗(yàn)證等幾何建模方法所得到的轉(zhuǎn)子模型用于柔性轉(zhuǎn)子控制的可行性:一是通過(guò)奇異值分析對(duì)比低階模型與高階模型的差異,并且在簡(jiǎn)支條件下進(jìn)行模型對(duì)比驗(yàn)證;二是假設(shè)Timoshenko轉(zhuǎn)子采用電磁軸承支撐,考慮轉(zhuǎn)子的分布質(zhì)量不平衡,將所得的低階、高階等幾何Timoshenko轉(zhuǎn)子模型與電磁軸承力學(xué)模型一起采用PD控制算法進(jìn)行控制仿真,并比較兩者差異。

3.1 狀態(tài)空間模型及轉(zhuǎn)子參數(shù)

令

結(jié)合式(3)、式(5),得到式(2)的狀態(tài)空間方程

(12)

式中

式(12)通過(guò)拉普拉斯變換得到系統(tǒng)的傳遞函數(shù)

G(s)=C(sI-A)-1B

(13)

本文所采用的Timoshenko轉(zhuǎn)子主要參數(shù)見(jiàn)表2。

3.2 奇異值分析

奇異值分析可以得到系統(tǒng)傳遞函數(shù)隨頻率變化的奇異值響應(yīng)。給定傳遞函數(shù)的奇異值具有唯一性,因此可用來(lái)刻畫(huà)系統(tǒng)降階前后輸入、輸出響應(yīng)的逼近程度。先將s=jω代入所得傳遞函數(shù)矩陣式(13),再進(jìn)行奇異值分解可得到奇異值響應(yīng)函數(shù)。

表2 Timoshenko轉(zhuǎn)子參數(shù)

注:帶有“#”的參數(shù)由不帶“#”的參數(shù)計(jì)算得到。

Timoshenko轉(zhuǎn)子兩端在x、y方向上都有拉伸彈簧以及拉伸阻尼器,且輸入廣義力的位置以及廣義位移觀測(cè)位置也都在轉(zhuǎn)子兩端。拉伸彈簧剛度系數(shù)為2×106N/m、拉伸阻尼器阻尼系數(shù)為1×103kg/s,其他參數(shù)見(jiàn)表2。使用的等幾何單元數(shù)量ne分別為2和32,B樣條基函數(shù)階次p=8。該轉(zhuǎn)子系統(tǒng)的奇異值響應(yīng)函數(shù)曲線如圖1所示。

圖1 高階、低階轉(zhuǎn)子系統(tǒng)傳遞函數(shù)的奇異值響應(yīng)

從圖1中可以看出,當(dāng)?shù)葞缀螁卧獢?shù)目取為32和2、基函數(shù)階次取為8時(shí),即式(2)的階數(shù)分別為160和40,系統(tǒng)傳遞函數(shù)的奇異值響應(yīng)在頻率小于3 000 rad/s范圍內(nèi)差異很小。采用等幾何方法對(duì)該模型進(jìn)行有阻尼模態(tài)分析(ne=32,p=8),可得到Campbell圖,從而得到轉(zhuǎn)子前6階臨界轉(zhuǎn)速,分別為123.76、418.40、736.40、1 110.2、1 748.7、2 722.5 rads-1。3 000 rad/s已超過(guò)轉(zhuǎn)子系統(tǒng)的第6階臨界轉(zhuǎn)速,所以研究前幾個(gè)低階模態(tài)頻率范圍內(nèi)的問(wèn)題時(shí),直接使用等幾何Timoshenko轉(zhuǎn)子模型作為控制模型是可行的。

3.3 模型驗(yàn)證

通過(guò)對(duì)比轉(zhuǎn)子模態(tài)頻率的數(shù)值解與理論解,以驗(yàn)證所建立等幾何Timoshenko梁模型的準(zhǔn)確性,采用模型參數(shù)如表1所示,邊界條件為兩端簡(jiǎn)支。

等截面Timoshenko梁在兩端簡(jiǎn)支條件下無(wú)量綱ω*模態(tài)頻率滿足下式[17]

(14)

式中

i∈N*;ω0=(π/L)2(EIy/(ρA))1/2

ω*=ω/ω0;Ω*=Ω/ω0

由式(2)、式(14)可分別得到無(wú)量綱模態(tài)頻率的數(shù)值解(ne=32,p=8)和理論解,其中前10階前向無(wú)量綱模態(tài)頻率如表3所示,兩者相對(duì)誤差絕對(duì)值最大為1.52×10-11。表4為高、低階等幾何Timoshenko梁模型的前6階前向無(wú)量綱模態(tài)頻率的數(shù)值解比較,前4階誤差在0.2%以內(nèi),但第5、6階相對(duì)誤差分別達(dá)到了6.42%與10.8%?？梢?jiàn),低階等幾何Timoshenko模型的準(zhǔn)確性在前4階彎曲模態(tài)頻率范圍內(nèi)的精度可以保證。若要求得更高的控制精度,則應(yīng)通過(guò)模態(tài)截?cái)嗟确椒ǐ@得更高精度的低階等幾何Timoshenko模型。

表3 前10階無(wú)量綱前向模態(tài)頻率的數(shù)值解與理論解

注:加粗?jǐn)?shù)值為2種解的不同部分。

3.4 理想PD控制仿真

Timoshenko轉(zhuǎn)子采用一對(duì)相同的主動(dòng)徑向磁懸浮軸承(AMB)支撐Timoshenko轉(zhuǎn)子,磁懸浮軸承采用分散重合PID控制,即每個(gè)磁懸浮軸承的x、

表4 高、低階模型無(wú)量綱前向模態(tài)頻率的

y方向的電磁力都用單通道PID控制,共4個(gè)通道。為方便考察高、低Timoshenko梁模型在相同控制下的差異,本文不考慮非并置問(wèn)題和重力效應(yīng),采用可行的最簡(jiǎn)單控制方法,作出如下設(shè)置:①傳感器軸向位置與磁懸浮位置相同,都在轉(zhuǎn)子兩端;②不考慮轉(zhuǎn)子重力影響;③PID控制采用最簡(jiǎn)單的理想PD控制器。磁懸浮系統(tǒng)的控制原理如圖2所示。磁懸浮軸承的電磁力采用線性化模型

Fx=-Ksx+Kiix;Fy=-Ksy+Kiiy

(15)

式中:位移剛度Ks=-2 MN/m;電流剛度Ki=250 N/A。

圖2 磁懸浮軸承-轉(zhuǎn)子系統(tǒng)反饋控制框圖

假設(shè)簡(jiǎn)化的功率放大器、傳感器傳遞函數(shù)為

P(s)=KPA;S(s)=Ksen

(16)

式中:功率放大器直流增益KPA=1 A·V-1;位移傳感器直流增益Ksen=7 800 V/m。

采用的理想PD控制器

C(s)=KP+KDs

(17)

式中:比例放大系數(shù)KP=2.051 3;微分系數(shù)KD=5.128 2×10-4。

比例放大系數(shù)KP使得磁軸承的等效剛度為K=2 MN/m,而微分系數(shù)KD使得磁軸承的等效阻尼為C=1 000 kg/s,此時(shí)磁懸浮軸承仍可看成拉伸彈簧-阻尼組合單元。

假設(shè)轉(zhuǎn)子有恒定工作轉(zhuǎn)速Ω,亞臨界轉(zhuǎn)速取第一階臨界轉(zhuǎn)速的一半,超臨界轉(zhuǎn)速分別取相鄰臨界轉(zhuǎn)速的算術(shù)平均值,如表5所示。

以分布不平衡質(zhì)量引起的不平衡力作為擾動(dòng)源,假設(shè)式(10)中的x、y方向上的偏心沿z軸分布ex(z)、ey(z)分別為

(18)

(a)Ω=61.880 rad/s

(b)Ω=271.080 rad/s

(c)Ω=577.265 rad/s圖3 高階轉(zhuǎn)子模型0.5 s內(nèi)不同轉(zhuǎn)速下轉(zhuǎn)子跨中橫截面形心軌跡

以式(17)的PD算法作為磁懸浮軸承的控制策略,假設(shè)在轉(zhuǎn)子初始位移、速度均為0的條件下,分別以高、低階等幾何Timoshenko梁模型作為被控對(duì)象,得到了0.5 s內(nèi)轉(zhuǎn)子跨中截面的形心軌跡以及在x、y軸方向上的位移隨時(shí)間的響應(yīng)偏差曲線,其中高階等幾何Timoshenko梁模型對(duì)應(yīng)ne=32、p=8,低階模型對(duì)應(yīng)于ne=2、p=8。以高階模型為例,3種不同工作轉(zhuǎn)速下0.5 s內(nèi)轉(zhuǎn)子跨中截面的形心軌跡如圖3所示,ux、uy分別為x、y方向上的位移。從圖中可以看出:截面形心軌跡在某一方向上有最大幅值,且隨工作轉(zhuǎn)速增加而增大;0.5 s內(nèi)亞臨界工作轉(zhuǎn)速下的幅值最大方向與超臨界轉(zhuǎn)速條件下明顯不同。

在高、低階等幾何Timoshenko梁模型條件下,轉(zhuǎn)子在3種工作轉(zhuǎn)速時(shí)x、y方向上位移響應(yīng)偏差d曲線見(jiàn)圖4。由于采用的模型階數(shù)不同,因此可以反映高、低階等幾何Timoshenko梁模型的偏離程度。如圖4所示,除在起始階段誤差曲線呈不規(guī)則振動(dòng)外,后續(xù)時(shí)刻曲線呈類似周期性波動(dòng),波動(dòng)周期隨工作轉(zhuǎn)速增加而減小,后續(xù)曲線的最大波動(dòng)幅值隨工作轉(zhuǎn)速增加而增大。6種工作轉(zhuǎn)速下x、y方向上平均波動(dòng)頻率(ω=2π/T,T為平均波動(dòng)周期)及最大位移波動(dòng)幅值xm、ym如表5所示,平均波動(dòng)周期由y方向誤差曲線的波峰得到。從表5中可以看出:在誤差允許條件下,平均波動(dòng)頻率與轉(zhuǎn)子工作轉(zhuǎn)速相一致,最大相對(duì)誤差為0.36%;x、y方向的位移波動(dòng)幅值隨工作轉(zhuǎn)速增加而增大,但幅值始終很小。因此,采用低階等幾何Timoshenko梁模型得到的不平衡振動(dòng)位移穩(wěn)態(tài)響應(yīng)與高階模型下的穩(wěn)態(tài)響應(yīng)偏差很小,在控制中可看成是工作轉(zhuǎn)速下的小擾動(dòng)。

表5 不同工作轉(zhuǎn)速下位移偏差的波動(dòng)頻率及幅值

由于不平衡質(zhì)量引起的振動(dòng)是一種同頻的強(qiáng)迫振動(dòng),所以穩(wěn)態(tài)位移偏差波動(dòng)頻率與工作轉(zhuǎn)速相近,轉(zhuǎn)子振動(dòng)的穩(wěn)態(tài)解的頻率與強(qiáng)迫振動(dòng)的頻率一致[18],而高、低階等幾何Timoshenko梁模型下的穩(wěn)態(tài)位移偏差相當(dāng)于同頻振動(dòng)的疊加,所以疊加后其穩(wěn)態(tài)振動(dòng)頻率不變。偏差幅值隨工作轉(zhuǎn)速增加而增大,這是因?yàn)楦?、低階模型偏差總體隨工作頻率增加而增大。在前6階臨界轉(zhuǎn)速范圍內(nèi)這種偏差很小(可以從圖1中反應(yīng)出來(lái)),這是高、低階模型下,兩者的穩(wěn)態(tài)位移偏差很小的原因。

(a)Ω=61.880 rad/s

(b)Ω=271.080 rad/s

(c)Ω=577.265 rad/s圖4 高、低階轉(zhuǎn)子模型下轉(zhuǎn)子跨中形心在x、y方向上的位移偏差

4 結(jié) 論

(1)在簡(jiǎn)支條件下,高階等幾何Timoshenko梁模型在前10階模態(tài)頻率范圍內(nèi)具有很高的精度,低階模型在前4階模態(tài)頻率范圍內(nèi)誤差較小,最大相對(duì)誤差小于0.2%。

(2)在前6階臨界轉(zhuǎn)速范圍內(nèi),低階等幾何Timoshenko梁模型與對(duì)應(yīng)的高階模型奇異值頻率響應(yīng)基本一致。

(3)在以分布質(zhì)量不平衡力作為擾動(dòng)、分散PD算法作為磁懸浮軸承仿真控制策略的條件下,高、低階等幾何Timoshenko梁模型下轉(zhuǎn)子的位移響應(yīng)頻率與工作頻率最大相對(duì)誤差為0.36%,位移波動(dòng)幅值最大誤差在1.9 nm以內(nèi),由低階等幾何Timoshenko梁模型引起的誤差可看成是工作轉(zhuǎn)速下的同頻小擾動(dòng)。

(4)低階等幾何Timoshenko梁模型在前4階模態(tài)頻率范圍內(nèi)能保證模型的精度,該方法所得的低階模型用于振動(dòng)控制的可行性得到初步驗(yàn)證,為模型在控制中應(yīng)用做好鋪墊。

[1] 孟光. 轉(zhuǎn)子動(dòng)力學(xué)研究的回顧與展望 [J]. 振動(dòng)工程學(xué)報(bào), 2002, 15(1): 1-2. MENG Guang. Review and prospect of the research on rotor dynamics [J]. Journal of Vibration Engineering, 2002, 15(1): 1-2.

[2] SCHWEITZER G, BLEULER H, MASLEN E H, et al. Magnetic bearings: theory, design, and application to rotating machinery [M]. Berlin, Germany: Springer, 2009: 1-11.

[3] 曹建榮, 虞烈, 謝友柏. 主動(dòng)磁懸浮軸承的解耦控制 [J]. 西安交通大學(xué)學(xué)報(bào), 1999, 33(12): 44-48. CAO Jianrong, YU Lie, XIE Youbai. Decoupling control for active magnetic bearing [J]. Journal of Xi’an Jiaotong University, 1999, 33(12): 44-48.

[4] 俞文伯, 欒勝, 房建成. CMG磁懸浮轉(zhuǎn)子系統(tǒng)的模型與控制規(guī)律 [J]. 航空學(xué)報(bào), 2003, 24(6): 541-545. YU Wenbo, LUAN Sheng, FANG Jiancheng. Model and control strategy of CMG magnetic bearing-rotor system [J]. Acta Aeronautica et Astronautica Sinica, 2003, 24(6): 541-545.

[5] 張剴, 趙雷, 趙鴻賓. 磁軸承飛輪控制系統(tǒng)設(shè)計(jì)中LQR方法的應(yīng)用研究 [J]. 機(jī)械工程學(xué)報(bào), 2004, 40(2): 127-131. ZHANG Kai, ZHAO Lei, ZHAO Hongbin. Study on the application of LQR method in the design of magnetic bearing flywheel control system [J]. Journal of Mechanical Engineering, 2004, 40(2): 127-131.

[6] 方鵬, 蔣科堅(jiān). 采用電磁軸承控制柔性轉(zhuǎn)子臨界轉(zhuǎn)速分布的建模和仿真分析 [J]. 浙江理工大學(xué)學(xué)報(bào)(自然科學(xué)), 2014, 31(3): 221-227. FANG Peng, JIANG Kejian. Modeling and simulation analysis of controlling distribution of critical speed of flexible rotor with electromagnetic bearing [J]. Journal of Zhejiang Sci-Tech University (Natural Sciences), 2014, 31(3): 221-227.

[7] MUSHI S E, LIN Zhongli, ALLAIRE P E. Design, construction, and modeling of a flexible rotor active magnetic bearing test rig [J]. IEEE/ASME Transactions on Mechatronics, 2012, 17(6): 1170-1182.

[8] HUGHES T J R, COTTRELL J A, BAZILEVS Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement [J]. Computer Methods in Applied Mechanics and Engineering, 2005, 194(39): 4135-4195.

[9] 尹碩輝, 余天堂, 劉鵬. 基于等幾何有限元法的功能梯度板自由振動(dòng)分析 [J]. 振動(dòng)與沖擊, 2013, 32(24): 180-186. YIN Shuohui, YU Tiantang, LIU Peng. Free vibration analysis of functionally graded plates using isogeometric finite element method [J]. Journal of Vibration and Shock, 2013, 32(24): 180-186.

[10]陳德祥, 徐自力, 劉石, 等. 求解Stokes方程的最小二乘等幾何分析方法 [J]. 西安交通大學(xué)學(xué)報(bào), 2013, 47(5): 51-55. CHEN Dexiang, XU Zili, LU Shi, et al. Least squares isogeometric analysis for Stokes equation [J]. Journal of Xi’an Jiaotong University, 2013, 47(5): 51-55.

[11]陳德祥, 徐自力. 多重網(wǎng)格在黏性流動(dòng)最小二乘等幾何模擬中的應(yīng)用 [J]. 西安交通大學(xué)學(xué)報(bào), 2014, 48(11): 122-127. CHEN Dexiang, XU Zili. An application of multigrid in viscous flow simulations with least squares isogeometric method [J]. Journal of Xi’an Jiaotong University, 2014, 48(11): 122-127.

[12]COTTRELL J A, HUGHES T J R, BAZILEVS Y. Isogeometric analysis: toward integration of CAD and FEA [M]. Hoboken, NJ, USA: John Wiley & Sons, 2009: 164-165.

[13]BEIRAO V L, LOVADINA C, REALI A. Avoiding shear locking for the Timoshenko beam problem via isogeometric collocation methods [J]. Computer Methods in Applied Mechanics and Engineering, 2012, 241: 38-51.

[14]LEE S J, PARK K S. Vibrations of Timoshenko beams with isogeometric approach [J]. Applied Mathematical Modelling, 2013, 37(22): 9174-9190.

[15]劉偉. 結(jié)構(gòu)振動(dòng)超收斂等幾何分析方法 [D]. 福建廈門(mén): 廈門(mén)大學(xué), 2014: 63-75.

[16]GENTA G. Dynamics of rotating systems [M]. Berlin, Germany: Springer Science & Business Media, 2007: 213-217.

[17]PIEGL L, TILLER W. The NURBS book [M]. Berlin, Germany: Springer-Verlag, 1997: 47-49.

[18]RAO S S. 機(jī)械振動(dòng) [M]. 4版. 李欣業(yè), 張明錄, 譯. 北京: 清華大學(xué)出版社, 2009: 357-358.

(編輯 杜秀杰)

Isogeometric Timoshenko Beam Model and Numerical Verification for Active Vibration Control of Flexible Rotor

ZHANG Pengfei,WANG Zhiheng,XI Guang

(School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China)

A Timoshenko beam model with fewer degrees of freedom (DOFs) by isogeometric analysis method is proposed to be applied to active control of the flexible rotor and is numerically verified. This strategy solves the difficulties in the model from the standard FEM analysis, which cannot be directly employed for controller design due to the more DOFs. The isogeometric beam model with the semi-discrete form and the various boundaries are considered as the input signals in the active control for the rotor, then the singular value responses of the isogeometric beam models with more and fewer DOFs are compared, and the numerical solutions of isogeometric beam models with more and fewer DOFs and theoretic solutions both with simple supports are also compared. The rotor supported by the magnetic bearing and under the disturbances by the given distributing unbalance forces is simulated by the decentralized PD method. The results show that the singular responses of the model with fewer DOFs almost agree with those with more DOFs within the range of the sixth critical speed; the numerical solutions with more DOFs agree well the theoretic solutions within the first ten modal frequencies and the relative error of the first four modal frequencies from the model with fewer DOFs is less than 0.2%; the deviation of the steady displacement responses from the unbalance vibration between the isogeometric models with more and fewer DOFs is small and their steady response deviation can be regarded as the small disturbances with the same frequency at the working speed. The feasibility of applying directly for controller design is verified by the higher accuracy of the fewer DOFs isogeometric beam model within the low frequency range.

isogeometric method; flexible rotor model; Timoshenko beam; vibration control; magnetic bearing

2016-01-15。

張鵬飛(1989—),男,博士生;席光(通信作者),男,教授,博士生導(dǎo)師。

國(guó)家自然科學(xué)基金資助項(xiàng)目(51236006);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(XJJ2013002)。

時(shí)間:2016-07-14

http:∥www.cnki.net/kcms/detail/61.1069.T.20160714.1721.014.html

10.7652/xjtuxb201610021

O347.6; TH133.3

A

0253-987X(2016)10-0139-08