聶衛(wèi)健，鄧旺群，徐友良，楊　海，郭天才（中國(guó)航空動(dòng)力機(jī)械研究所航空發(fā)動(dòng)機(jī)振動(dòng)技術(shù)航空科技重點(diǎn)實(shí)驗(yàn)室，湖南株洲412002）

高速柔性轉(zhuǎn)子臨界轉(zhuǎn)速隨支承剛度和輪盤質(zhì)量的變化規(guī)律

聶衛(wèi)健，鄧旺群，徐友良，楊海，郭天才
（中國(guó)航空動(dòng)力機(jī)械研究所航空發(fā)動(dòng)機(jī)振動(dòng)技術(shù)航空科技重點(diǎn)實(shí)驗(yàn)室，湖南株洲412002）

摘要：以小型渦扇發(fā)動(dòng)機(jī)模擬低壓轉(zhuǎn)子為研究對(duì)象，采用梁?jiǎn)卧⑵鋭?dòng)力特性的有限元計(jì)算模型，并在不同的支承剛度及輪盤質(zhì)量下，運(yùn)用轉(zhuǎn)子動(dòng)力學(xué)分析軟件SAMCEF/ROTOR，對(duì)模擬低壓轉(zhuǎn)子的前三階臨界轉(zhuǎn)速進(jìn)行了系統(tǒng)計(jì)算及分析，揭示出了模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速隨各支承剛度、各輪盤質(zhì)量的變化規(guī)律，為模擬低壓轉(zhuǎn)子的臨界轉(zhuǎn)速設(shè)計(jì)和調(diào)整提供了理論依據(jù)，可為后續(xù)全尺寸模擬低壓轉(zhuǎn)子的動(dòng)力特性試驗(yàn)提供指導(dǎo)。

關(guān)鍵詞：渦扇發(fā)動(dòng)機(jī)；模擬低壓轉(zhuǎn)子；轉(zhuǎn)子動(dòng)力學(xué)；支承剛度；輪盤質(zhì)量；臨界轉(zhuǎn)速；變化規(guī)律

1　引言

轉(zhuǎn)子動(dòng)力特性分析是轉(zhuǎn)子動(dòng)力學(xué)研究的一項(xiàng)重要內(nèi)容，國(guó)內(nèi)學(xué)者在該領(lǐng)域開(kāi)展了大量的研究工作。黃太平［1］采用傳遞矩陣——擬模態(tài)綜合法，計(jì)算了渦扇發(fā)動(dòng)機(jī)雙轉(zhuǎn)子支承系統(tǒng)的臨界轉(zhuǎn)速。史亞杰等［2］對(duì)柔性轉(zhuǎn)子的動(dòng)力特性進(jìn)行了研究。鄧旺群等［3-6］對(duì)航空發(fā)動(dòng)機(jī)等高速柔性轉(zhuǎn)子的動(dòng)力特性及高速動(dòng)平衡技術(shù)，進(jìn)行了系統(tǒng)的分析和試驗(yàn)驗(yàn)證。白中祥、洪杰、盛步云等［7-9］研究了支承系統(tǒng)剛度或動(dòng)剛度對(duì)轉(zhuǎn)子動(dòng)力特性的影響。目前，常采用有限元法建立復(fù)雜轉(zhuǎn)子的動(dòng)力分析模型。繆輝等［10］運(yùn)用有限元法，對(duì)雙轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速進(jìn)行了分析。

臨界轉(zhuǎn)速設(shè)計(jì)是中小型航空發(fā)動(dòng)機(jī)轉(zhuǎn)子設(shè)計(jì)的核心內(nèi)容，在中小型航空發(fā)動(dòng)機(jī)研制中具有十分重要的地位。為滿足轉(zhuǎn)子臨界轉(zhuǎn)速設(shè)計(jì)的裕度要求，可通過(guò)選取合適支承剛度、改變質(zhì)量分布、優(yōu)化轉(zhuǎn)子結(jié)構(gòu)等方法來(lái)實(shí)現(xiàn)［11］。

本文以小型渦扇發(fā)動(dòng)機(jī)的模擬低壓轉(zhuǎn)子為對(duì)象，研究了其前三階臨界轉(zhuǎn)速隨各支承剛度、各輪盤質(zhì)量的變化規(guī)律，可為臨界轉(zhuǎn)速設(shè)計(jì)和調(diào)整提供理論依據(jù)，為后續(xù)全尺寸模擬低壓轉(zhuǎn)子的動(dòng)力特性試驗(yàn)提供指導(dǎo)。

2　模擬低壓轉(zhuǎn)子結(jié)構(gòu)

模擬低壓轉(zhuǎn)子結(jié)構(gòu)示意圖見(jiàn)圖1。整個(gè)轉(zhuǎn)子由兩級(jí)風(fēng)扇模擬盤、增壓級(jí)模擬盤、兩級(jí)低壓渦輪模擬盤、風(fēng)扇軸、拉桿及低壓渦輪軸等組成，其中風(fēng)扇模擬盤端懸臂。轉(zhuǎn)子采用4支點(diǎn)支承方案，軸承編號(hào)與發(fā)動(dòng)機(jī)中軸承編號(hào)一致。故該轉(zhuǎn)子是一個(gè)多盤、多支點(diǎn)、軸分段連接且一端懸臂的高速柔性轉(zhuǎn)子。

圖1　模擬低壓轉(zhuǎn)子結(jié)構(gòu)示意圖Fig.1 Structure sketch of the simulated LP rotor

3　有限元計(jì)算模型

針對(duì)模擬低壓轉(zhuǎn)子的復(fù)雜結(jié)構(gòu)，采用梁?jiǎn)卧⑵鋭?dòng)力特性的有限元計(jì)算模型（圖2）。建模時(shí)，簡(jiǎn)化轉(zhuǎn)子結(jié)構(gòu)，并忽略一些細(xì)小的局部結(jié)構(gòu)（如倒角、小孔等）；用集中質(zhì)量單元模擬部分兩級(jí)風(fēng)扇模擬盤和渦輪模擬盤，同時(shí)用軸承單元模擬支承。

圖2　模擬低壓轉(zhuǎn)子有限元計(jì)算模型Fig.2 Finite element calculation model of the simulated LP rotor

4　臨界轉(zhuǎn)速隨支承剛度的變化規(guī)律

通過(guò)動(dòng)力特性計(jì)算分析，初步確定了模擬低壓轉(zhuǎn)子1號(hào)、2號(hào)、5號(hào)和6號(hào)支承的剛度，見(jiàn)表1。在該組合支承剛度下，模擬低壓轉(zhuǎn)子的額定工作轉(zhuǎn)速高于前兩階臨界轉(zhuǎn)速，并滿足設(shè)計(jì)準(zhǔn)則要求。但轉(zhuǎn)子的實(shí)際臨界轉(zhuǎn)速是否滿足設(shè)計(jì)要求，還需全轉(zhuǎn)速范圍內(nèi)的試驗(yàn)驗(yàn)證。若不滿足設(shè)計(jì)要求，可通過(guò)調(diào)整轉(zhuǎn)子的支承剛度和/或質(zhì)量分布使其滿足。在不改變質(zhì)量分布和轉(zhuǎn)子結(jié)構(gòu)的前提下，以表1中的支承剛度為基準(zhǔn)，通過(guò)計(jì)算揭示模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速隨各支承剛度的變化規(guī)律，為基于支承剛度調(diào)整轉(zhuǎn)子的臨界轉(zhuǎn)速提供理論依據(jù)。

表1　模擬低壓轉(zhuǎn)子支承剛度Table 1 Supporting stiffness of the simulated LP rotor

4.1臨界轉(zhuǎn)速隨1號(hào)支承剛度的變化規(guī)律

在2號(hào)、5號(hào)及6號(hào)支承剛度保持不變的情況下，1號(hào)支承剛度在（0.1～50.0）×107N/m范圍內(nèi)變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的計(jì)算結(jié)果見(jiàn)表2。當(dāng)1號(hào)支承剛度在一定區(qū)段內(nèi)變化時(shí)，根據(jù)表2，可得到模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的變化率（表3）。

表2　改變1號(hào)支承剛度時(shí)的臨界轉(zhuǎn)速Table 2 Critical speeds of the rotor when changing No.1 supporting stiffness

表3　改變1號(hào)支承剛度時(shí)臨界轉(zhuǎn)速的變化率Table 3 Change rate of the critical speeds when changing No.1 supporting stiffness

從表2和表3可以看出：

（1）1號(hào)支承剛度從0.1×107N/m增大到5.0× 107N/m時(shí)，模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速顯著變化，提高了267.75%；隨著1號(hào)支承剛度的繼續(xù)增大，模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速幾乎不再發(fā)生變化。

（2）1號(hào)支承剛度從0.1×107N/m增大到5.0× 107N/m、從5.0×107N/m增大到15.0×107N/m、從15.0×107N/m增大到50.0×107N/m時(shí)，模擬低壓轉(zhuǎn)子第二階臨界轉(zhuǎn)速分別提高了5.84%、17.41%和6.91%，累計(jì)提高了30.16%。

（3）1號(hào)支承剛度從0.1×107N/m增大到5.0× 107N/m時(shí)，模擬低壓轉(zhuǎn)子第三階臨界轉(zhuǎn)速提高了10.28%；隨著1號(hào)支承剛度的繼續(xù)增大，模擬低壓轉(zhuǎn)子第三階臨界轉(zhuǎn)速幾乎不再發(fā)生變化。

4.2臨界轉(zhuǎn)速隨2號(hào)支承剛度的變化規(guī)律

在1號(hào)、5號(hào)及6號(hào)支承剛度保持不變的情況下，2號(hào)支承剛度在（0.1～50.0）×107N/m范圍內(nèi)變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的計(jì)算結(jié)果見(jiàn)表4。當(dāng)2號(hào)支承剛度在一定區(qū)段內(nèi)變化時(shí)，根據(jù)表4，可得到模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的變化率（表5）。

表4　改變2號(hào)支承剛度時(shí)的臨界轉(zhuǎn)速Table 4 Critical speeds of the rotor when changing No.2 supporting stiffness

表5　改變2號(hào)支承剛度時(shí)臨界轉(zhuǎn)速的變化率Table 5 Change rate of the critical speeds when changing No.2 supporting stiffness

從表4和表5可以看出：2號(hào)支承剛度從0.1× 107N/m增大到3.0×107N/m時(shí)，模擬低壓轉(zhuǎn)子第一、第二階臨界轉(zhuǎn)速分別增大了5.22%和18.83%，而第三階臨界轉(zhuǎn)速顯著提高了109.03%。隨著2號(hào)支承剛度的繼續(xù)增大，前三階臨界轉(zhuǎn)速不再有明顯變化。

4.3臨界轉(zhuǎn)速隨5號(hào)支承剛度的變化規(guī)律

在1號(hào)、2號(hào)及6號(hào)支承剛度保持不變的情況下，5號(hào)支承剛度在（0.1～50.0）×107N/m范圍內(nèi)變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的計(jì)算結(jié)果見(jiàn)表6。當(dāng)5號(hào)支承剛度在一定區(qū)段內(nèi)變化時(shí)，根據(jù)表6，可得到低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的變化率（表7）。

表6　改變5號(hào)支承剛度時(shí)的臨界轉(zhuǎn)速Table 6 Critical speeds of the rotor when changing No.5 supporting stiffness

表7　改變5號(hào)支承剛度時(shí)臨界轉(zhuǎn)速的變化率Table 7 Change rate of the critical speeds when changing No.5 supporting stiffness

從表6和表7可以看出：

（1）5號(hào)支承剛度從0.1×107N/m增大到50.0× 107N/m時(shí)，模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速?zèng)]有實(shí)質(zhì)性變化。

（2）5號(hào)支承剛度從0.1×107N/m增大到3.0× 107N/m時(shí)，模擬低壓轉(zhuǎn)子第二階臨界轉(zhuǎn)速顯著提高了186.62%；5號(hào)支承剛度從3.0×107N/m增大到8.0× 107N/m時(shí)，第二階臨界轉(zhuǎn)速提高了14.96%；隨著5號(hào)支承剛度的繼續(xù)增大，第二階臨界轉(zhuǎn)速不再明顯變化。

（3）5號(hào)支承剛度從0.1×107N/m增大到3.0× 107N/m時(shí)，模擬低壓轉(zhuǎn)子第三階臨界轉(zhuǎn)速?zèng)]有實(shí)質(zhì)性的變化；5號(hào)支承剛度從3.0×107N/m增大到8.0× 107N/m時(shí)，第三階臨界轉(zhuǎn)速提高了10.48%；隨著5號(hào)支承剛度的繼續(xù)增大，第三階臨界轉(zhuǎn)速不再有明顯變化。

4.4 臨界轉(zhuǎn)速隨6號(hào)支承剛度的變化規(guī)律

在1號(hào)、2號(hào)及5號(hào)支承剛度保持不變的情況下，6號(hào)支承剛度在（0.1～50.0）×107N/m范圍內(nèi)變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的計(jì)算結(jié)果見(jiàn)表8。當(dāng)6號(hào)支承剛度在一定區(qū)段內(nèi)變化時(shí)，根據(jù)表8，可得到低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的變化率（表9）。

表8　改變6號(hào)支承剛度時(shí)的臨界轉(zhuǎn)速Table 8 Critical speeds of the rotor when changing No.6 supporting stiffness

表9　改變6號(hào)支承剛度時(shí)臨界轉(zhuǎn)速的變化率Table 9 Change rate of the critical speeds when changing No.6 supporting stiffness

從表8和表9可以看出：

（1）6號(hào)支承剛度從0.1×107N/m增大到50.0× 107N/m，模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速?zèng)]有實(shí)質(zhì)性變化。

（2）6號(hào)支承剛度從0.1×107N/m增大到5.0× 107N/m，模擬低壓轉(zhuǎn)子第二、第三階臨界轉(zhuǎn)速分別提高了75.27%和41.47%。隨著6號(hào)支承剛度的繼續(xù)增大，第二、第三階臨界轉(zhuǎn)速都不再有明顯變化。

5　臨界轉(zhuǎn)速隨輪盤質(zhì)量的變化規(guī)律

在不改變支承剛度和輪盤結(jié)構(gòu)的前提下，通過(guò)計(jì)算揭示模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速隨各輪盤質(zhì)量（輪盤結(jié)構(gòu)不變的情況下，只取決于材料的密度）的變化規(guī)律，為模擬低壓轉(zhuǎn)子基于輪盤質(zhì)量的臨界轉(zhuǎn)速調(diào)整提供理論依據(jù)。計(jì)算使用的支承剛度見(jiàn)表1，選用的輪盤材料及其屬性見(jiàn)表10。

表10　選用的輪盤材料及其屬性Table 10 Properties of the disk materials

5.1臨界轉(zhuǎn)速隨風(fēng)扇一級(jí)模擬盤質(zhì)量的變化規(guī)律

風(fēng)扇一級(jí)模擬盤密度在1 780～8 680 kg/m3范圍變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速及其變化率的計(jì)算結(jié)果分別如表11、表12所示。可見(jiàn)：風(fēng)扇一級(jí)模擬盤密度從1 780 kg/m3增大到8 680 kg/m3時(shí)，模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速顯著變化，減小幅度達(dá)33.47%；而第二、第三階臨界轉(zhuǎn)速?zèng)]有實(shí)質(zhì)性變化，減小幅度均小于5%。

5.2臨界轉(zhuǎn)速隨風(fēng)扇二級(jí)模擬盤質(zhì)量的變化規(guī)律

風(fēng)扇二級(jí)模擬盤密度在1 780～8 680 kg/m3范圍變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速及其變化率的計(jì)算結(jié)果分別如表13、表14所示。可見(jiàn)：風(fēng)扇二級(jí)模擬盤的密度從1 780 kg/m3增加到8 680 kg/m3時(shí)，模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速減小了7.42%，第三階臨界轉(zhuǎn)速增大了14.33%，而第二階臨界轉(zhuǎn)速幾乎不變。

5.3臨界轉(zhuǎn)速隨增壓級(jí)模擬盤質(zhì)量的變化規(guī)律

增壓級(jí)模擬盤的密度在1 780～8 680 kg/m3范圍變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速及其變化率的計(jì)算結(jié)果分別如表15、表16所示?？梢?jiàn)：增壓級(jí)模擬盤密度從1 780 kg/m3增加到8 680 kg/m3時(shí)，模擬低壓轉(zhuǎn)子第一、第二階臨界轉(zhuǎn)速的變化均小于2%，第三階臨界轉(zhuǎn)速減小了7.78%。

5.4臨界轉(zhuǎn)速隨低壓渦輪一級(jí)模擬盤質(zhì)量的變化規(guī)律

低壓渦輪一級(jí)模擬盤的密度在1 780～8 680kg/m3范圍變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速及其變化率的計(jì)算結(jié)果分別如表17、表18所示?？梢?jiàn)：低壓渦輪一級(jí)模擬盤密度從1 780 kg/m3增大到8 680 kg/m3時(shí)，模擬低壓轉(zhuǎn)子第一、第三階臨界轉(zhuǎn)速的變化均小于1%，而第二階臨界轉(zhuǎn)速顯著變化，減小了27.89%。

表11　風(fēng)扇一級(jí)模擬盤密度變化時(shí)的臨界轉(zhuǎn)速Table 11 Critical speeds of the rotor when changing the density of the simulated first stage disk of fan

表12　風(fēng)扇一級(jí)模擬盤密度變化時(shí)臨界轉(zhuǎn)速的變化率Table 12 Change rate of the critical speeds when changing the density of the simulated first stage disk of fan

表13　風(fēng)扇二級(jí)模擬盤密度變化時(shí)的臨界轉(zhuǎn)速Table 13 Critical speeds of the rotor when changing the density of the simulated second stage disk of fan

表14　風(fēng)扇二級(jí)模擬盤密度變化時(shí)臨界轉(zhuǎn)速的變化率Table 14 Change rate of the critical speeds when changing the density of the simulated second stage disk of fan

表15　增壓級(jí)模擬盤密度變化時(shí)的臨界轉(zhuǎn)速Table 15 Critical speeds of the rotor when changing the density of the simulated booster stage disk

表16　增壓級(jí)模擬盤密度變化時(shí)臨界轉(zhuǎn)速的變化率Table 16 Change rate of the critical speeds when changing the density of simulated booster stage disk

表17　低壓渦輪一級(jí)模擬盤密度變化時(shí)的臨界轉(zhuǎn)速Table 17 Critical speeds of the rotor when changing the density of the simulated first stage disk of LP turbine

表18　低壓渦輪一級(jí)模擬盤密度變化時(shí)臨界轉(zhuǎn)速的變化率Table 18 Change rate of the critical speeds when changing the density of the simulated first stage disk of LP turbine

5.5臨界轉(zhuǎn)速隨低壓渦輪二級(jí)模擬盤質(zhì)量的變化規(guī)律

低壓渦輪二級(jí)模擬盤的密度在1 780～8 680 kg/m3范圍變化時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速及其變化率的計(jì)算結(jié)果分別如表19、表20所示。可見(jiàn)：低壓渦輪二級(jí)模擬盤密度從1 780 kg/m3增大到8 680 kg/m3時(shí)，模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速的變化均小于2.5%。

表19　低壓渦輪二級(jí)模擬盤密度變化時(shí)的臨界轉(zhuǎn)速Table 19 Critical speeds of the rotor when changing the density of the simulated second stage disk of LP turbine

表20　低壓渦輪二級(jí)模擬盤密度變化時(shí)臨界轉(zhuǎn)速的變化率Table 20 Change rate of the critical speeds when changing the density of the simulated second stage disk of LP turbine

6　結(jié)論

通過(guò)計(jì)算分析，揭示了某小型渦扇發(fā)動(dòng)機(jī)模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速隨支承剛度、輪盤質(zhì)量的變化規(guī)律，主要結(jié)論如下：

（1）1號(hào)支承剛度對(duì)模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速、2號(hào)支承剛度對(duì)模擬低壓轉(zhuǎn)子第三階臨界轉(zhuǎn)速、5號(hào)支承剛度對(duì)模擬低壓轉(zhuǎn)子第二階臨界轉(zhuǎn)速、6號(hào)支承剛度對(duì)模擬低壓轉(zhuǎn)子第二階和第三階臨界轉(zhuǎn)速，均有顯著影響。

（2）在一定支承剛度范圍內(nèi)，可通過(guò)改變一個(gè)或幾個(gè)彈性支承的支承剛度，來(lái)調(diào)整模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速。

（3）風(fēng)扇一級(jí)模擬盤質(zhì)量對(duì)模擬低壓轉(zhuǎn)子第一階臨界轉(zhuǎn)速有顯著影響，風(fēng)扇二級(jí)模擬盤質(zhì)量對(duì)第一、第三階臨界轉(zhuǎn)速有一定影響，增壓級(jí)模擬盤質(zhì)量對(duì)第三階臨界轉(zhuǎn)速有一定影響，低壓渦輪一級(jí)模擬盤質(zhì)量對(duì)第二階臨界轉(zhuǎn)速有顯著影響，低壓渦輪二級(jí)模擬盤質(zhì)量對(duì)前三階臨界轉(zhuǎn)速均無(wú)實(shí)質(zhì)性影響。

（4）如需通過(guò)改變轉(zhuǎn)子質(zhì)量分布來(lái)調(diào)整模擬低壓轉(zhuǎn)子前三階臨界轉(zhuǎn)速，在保證輪盤強(qiáng)度和壽命的前提下，可通過(guò)改變一個(gè)或幾個(gè)輪盤的材料來(lái)實(shí)現(xiàn)。

參考文獻(xiàn)：

［1］ 黃太平.渦扇發(fā)動(dòng)機(jī)轉(zhuǎn)子動(dòng)力特性計(jì)算與分析［C］//.福建省科協(xié)第四屆學(xué)術(shù)年會(huì)——提升福建制造業(yè)競(jìng)爭(zhēng)力的戰(zhàn)略思考專題學(xué)術(shù)年會(huì)論文集，2004.

［2］ 史亞杰，王孝利，洪杰，等.柔性轉(zhuǎn)子動(dòng)力特性研究［J］.航空發(fā)動(dòng)機(jī)，2005，31（1）：14—17.

［3］ 鄧旺群，郭飛躍，高德平.航空發(fā)動(dòng)機(jī)高速柔性轉(zhuǎn)子動(dòng)力特性計(jì)算［J］.振動(dòng)與沖擊，2006，25（5）：130—133.

［4］ 鄧旺群，郭飛躍，高德平.高速柔性轉(zhuǎn)子的主要零部件對(duì)其動(dòng)力特性的影響分析［J］.機(jī)械強(qiáng)度，2006，28（6）：813—819.

［5］ 鄧旺群，高德平.懸臂柔性轉(zhuǎn)子動(dòng)力特性及高速動(dòng)平衡試驗(yàn)［J］.航空動(dòng)力學(xué)報(bào)，2006，21（3）：556—562.

［6］ 鄧旺群，王楨，舒斯榮，等.渦軸發(fā)動(dòng)機(jī)細(xì)長(zhǎng)柔性轉(zhuǎn)子動(dòng)力特性及高速動(dòng)平衡技術(shù)研究［J］.振動(dòng)與沖擊，2012，31（7）：162—165.

［7］ 白中祥，吳偉亮，劉閎釗.支承系統(tǒng)的剛度對(duì)轉(zhuǎn)子動(dòng)力特性的影響［J］.機(jī)械設(shè)計(jì)與研究，2012，28（4）：18—21.

［8］ 洪杰，王華，肖大為，等.轉(zhuǎn)子支承動(dòng)剛度對(duì)轉(zhuǎn)子動(dòng)力特性的影響分析［J］.航空發(fā)動(dòng)機(jī)，2008，34（1）：23—27.

［9］ 盛步云，張濤，丁毓峰，等.支承剛度對(duì)汽輪機(jī)轉(zhuǎn)子動(dòng)力學(xué)特性的影響分析［J］.機(jī)械設(shè)計(jì)，2008，25（12）：38—40.

［10］繆輝，王克明，艾書民，等.雙轉(zhuǎn)子系統(tǒng)臨界轉(zhuǎn)速的有限元分析［J］.沈陽(yáng)航空航天大學(xué)學(xué)報(bào)，2011，28（5）：27—31.

［11］顧家柳，丁奎元，劉啟周，等.轉(zhuǎn)子動(dòng)力學(xué)［M］.北京：國(guó)防工業(yè)出版社，1985：102—117.

中圖分類號(hào)：V231.96

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1672-2620（2015）03-0019-06

收稿日期：2014-10-20

基金項(xiàng)目：航空科學(xué)基金（20112108001，2013ZB08001）

作者簡(jiǎn)介：聶衛(wèi)健（1991-），男，江西臨川人，助理工程師，碩士研究生，研究方向?yàn)楹娇瞻l(fā)動(dòng)機(jī)強(qiáng)度試驗(yàn)與轉(zhuǎn)子動(dòng)力學(xué)。

Analysis on the changes of high-speed flexible rotor critical speeds with supporting stiffness and disk mass

NIE Wei-jian，DENG Wang-qun，XU You-liang，YANG Hai，GUO Tian-cai
（Aviation Key Laboratory of Aero-engine Vibration Technology，China Aviation Powerplant Research Institute，Zhuzhou 412002，China）

Abstract：Taking a simulated low-pressure rotor of a small turbofan engine as the researching object，a finite element calculation model on the rotor characteristics was established.The first three stage critical speeds of the LP rotor were calculated and analyzed by SAMCEF/ROTOR software at the conditions of different supporting stiffness and different disk mass.The change laws of the first three stage critical speeds with supporting stiffness and disk mass were revealed.The research results provide theoretical bases for critical speed design and critical speed adjustment of the rotor.

Key words：turbofan engine；simulated low-pressure rotor；rotor dynamics；supporting stiffness；disk mass；critical speed；change laws