陳 強(qiáng)，袁西冰

（1.貴州航天林泉電機(jī)有限公司，貴陽 550081；2.國家精密微特電機(jī)工程技術(shù)研究中心，貴陽 550081；3.貴州大學(xué)，貴陽 550000 ）

0 引言

飛輪儲能系統(tǒng)作為一種新型儲能技術(shù)，具有功率密度高、儲能密度大，性價(jià)比高，無污染，充電時(shí)間短且無過充電和過放電問題，壽命長且適合循環(huán)使用等優(yōu)點(diǎn)。在電網(wǎng)調(diào)峰、電動汽車、不間斷供電備用電源（UPS）、汽車供能、航空航天軍事等諸多方面都獲得了成功的應(yīng)用[1～7]。

國內(nèi)外學(xué)者對飛輪轉(zhuǎn)子的設(shè)計(jì)優(yōu)化進(jìn)行了大量研究。NASA Glenn中心在實(shí)驗(yàn)室環(huán)境研制成功磁懸浮復(fù)合材料飛輪儲能系統(tǒng)，轉(zhuǎn)速達(dá)到600000r/min[8]。UTCEM成功制作了復(fù)合材料飛輪轉(zhuǎn)子，最高轉(zhuǎn)速達(dá)到15000rpm，可儲存130kWh的能量[9]。HaSK基于平面應(yīng)力假設(shè)，計(jì)算了平面狀態(tài)下多環(huán)過盈裝配復(fù)合材料飛輪轉(zhuǎn)子的應(yīng)力分布，并以強(qiáng)度失效為準(zhǔn)則，最大儲能量為優(yōu)化目標(biāo)對多環(huán)過盈裝配復(fù)合材料飛輪進(jìn)行了優(yōu)化設(shè)計(jì)[10]。白越等驗(yàn)證了復(fù)合材料飛輪轉(zhuǎn)子中金屬輪轂與復(fù)合材料外環(huán)采用膠接的連接方式是合理的[11]。閆曉磊基于最優(yōu)控制理論，以最大儲能密度為優(yōu)化目標(biāo)對不同材料和結(jié)構(gòu)形狀的飛輪轉(zhuǎn)子進(jìn)行了結(jié)構(gòu)優(yōu)化設(shè)計(jì)[12]。目前，飛輪儲能系統(tǒng)正向著高速化和大功率方向發(fā)展，國際上飛輪轉(zhuǎn)速可達(dá)6×105r/min，線速度為8×102m/s以上。國內(nèi)飛輪轉(zhuǎn)速可達(dá)5.43×104r/min，線速度為796m/s以上[3]。在這樣的高速旋轉(zhuǎn)狀態(tài)下工作，飛輪轉(zhuǎn)子必然會產(chǎn)生強(qiáng)大的離心力，所引起的應(yīng)力若超過了飛輪轉(zhuǎn)子的強(qiáng)度極限，飛輪轉(zhuǎn)子被破壞，飛輪儲能系統(tǒng)出現(xiàn)故障。因此，對飛輪轉(zhuǎn)子進(jìn)行靜力學(xué)分析是非常重要的。

本文采用ANSYS Workbench對某公司所設(shè)計(jì)的飛輪儲能系統(tǒng)轉(zhuǎn)子進(jìn)行靜力學(xué)仿真分析，通過施加邊界條件，研究飛輪儲能系統(tǒng)轉(zhuǎn)子的應(yīng)力分布，對飛輪儲能系統(tǒng)的強(qiáng)度設(shè)計(jì)具有重要的意義。

1 飛輪儲能系統(tǒng)轉(zhuǎn)子有限元模型建立

1.1 飛輪儲能系統(tǒng)的結(jié)構(gòu)原理

飛輪儲能系統(tǒng)主要由飛輪轉(zhuǎn)子、電機(jī)、主軸、支撐系統(tǒng)和底座等部分組成，如圖1所示。其中，電機(jī)負(fù)責(zé)整個(gè)飛輪儲能系統(tǒng)的能量轉(zhuǎn)換，飛輪轉(zhuǎn)子實(shí)現(xiàn)能量存儲，軸承由角接觸球軸承和homopolar型永磁偏置徑向磁軸承共同作用。工作時(shí)首先由電機(jī)帶動飛輪加速旋轉(zhuǎn)，此時(shí)將驅(qū)動電機(jī)的電能轉(zhuǎn)化為機(jī)械能儲存在飛輪中，當(dāng)外界需要系統(tǒng)釋放能量時(shí)，飛輪利用慣性帶動電機(jī)進(jìn)行發(fā)電，將飛輪的機(jī)械能轉(zhuǎn)換成電能進(jìn)行輸出。當(dāng)飛輪空閑運(yùn)轉(zhuǎn)時(shí)，整個(gè)系統(tǒng)以最小的損耗運(yùn)行。

圖1 飛輪儲能系統(tǒng)結(jié)構(gòu)圖

1.2 有限元模型建立

1）材料物理參數(shù)。飛輪轉(zhuǎn)子及轉(zhuǎn)子兩端支撐均采用材料為高強(qiáng)度結(jié)構(gòu)鋼30CrMnSi，彈性模量E為200GPa，泊松比為0.3，密度為7750kg/m3。

2）有限元網(wǎng)格模型?？紤]到飛輪儲能系統(tǒng)結(jié)構(gòu)的軸對稱性，僅取周向四分之一模型進(jìn)行分析。通過UG建立飛輪儲能系統(tǒng)轉(zhuǎn)子的三維模型，并采用Hex Dominant劃分方法劃分網(wǎng)格，單元尺寸設(shè)置為10mm。飛輪轉(zhuǎn)子三維有限元網(wǎng)格模型如圖2所示，并對飛輪轉(zhuǎn)子兩端支撐采用同樣的方式進(jìn)行網(wǎng)格劃分。

圖2 飛輪轉(zhuǎn)子網(wǎng)格劃分模型

3）加載和邊界條件。飛輪儲能系統(tǒng)靜力學(xué)分析分為不考慮永磁體作用于轉(zhuǎn)子內(nèi)壁和考慮永磁體作用于轉(zhuǎn)子內(nèi)壁兩種情況。

（1）不考慮永磁體作用于轉(zhuǎn)子內(nèi)壁

不考慮永磁體作用于轉(zhuǎn)子內(nèi)壁，取對稱面約束，邊界條件如圖3所示。飛輪轉(zhuǎn)子是旋轉(zhuǎn)體，承受離心載荷。應(yīng)力分析時(shí)按轉(zhuǎn)速為6000rpm進(jìn)行分析，即：

（2）考慮永磁體作用于轉(zhuǎn)子內(nèi)壁

圖3 不考慮永磁體作用于內(nèi)壁的飛輪轉(zhuǎn)子模型

永磁體每片長度為80mm，安裝在轉(zhuǎn)子內(nèi)壁軸向長度為400mm的范圍，永磁體產(chǎn)生的離心力以外載荷形式作用于轉(zhuǎn)子內(nèi)壁。

每片永磁體體積為：

在轉(zhuǎn)速為6000rpm時(shí)，每片永磁體產(chǎn)生的離心力為：

每片永磁體的作用面積為：

作用于轉(zhuǎn)子內(nèi)壁的壓強(qiáng)為：

將壓強(qiáng)作用于轉(zhuǎn)子內(nèi)壁長度為400mm的區(qū)域上，載荷如圖4所示。

2 有限元仿真分析結(jié)果

1）不考慮永磁體作用于轉(zhuǎn)子內(nèi)壁的飛輪轉(zhuǎn)子靜力學(xué)分析

不考慮永磁體作用于轉(zhuǎn)子內(nèi)壁上時(shí)，飛輪儲能系統(tǒng)徑向位移如圖5所示，徑向應(yīng)力分布和環(huán)向應(yīng)力分布如圖6和圖7所示。由圖可以看出，飛輪變形最大位置在飛輪中部，最大變形量約0.3654mm。徑向應(yīng)力分布在-207.85MPa～440.41MPa，一部分為壓應(yīng)力，一部分為拉應(yīng)力。最大位置在內(nèi)壁倒角處。環(huán)向應(yīng)力分布均在330MPa左右，最大位置仍在螺栓孔處，此處為應(yīng)力集中。

圖4 考慮永磁體作用于內(nèi)壁的飛輪轉(zhuǎn)子模型

圖5 飛輪轉(zhuǎn)子徑向位移（不考慮永磁體作用于內(nèi)壁）

圖6 飛輪轉(zhuǎn)子徑向應(yīng)力分布（不考慮永磁體作用于內(nèi)壁）

圖7 飛輪轉(zhuǎn)子環(huán)向應(yīng)力分布（不考慮永磁體作用于內(nèi)壁）

2）考慮永磁體作用于轉(zhuǎn)子內(nèi)壁的飛輪轉(zhuǎn)子靜力學(xué)分析

考慮永磁體作用于轉(zhuǎn)子內(nèi)壁上時(shí)，飛輪儲能系統(tǒng)徑向位移如圖8所示，徑向應(yīng)力分布和環(huán)向應(yīng)力分布如圖9和圖10所示。由圖可知，飛輪本體在6000rpm時(shí)，其變形最大位置仍在飛輪中部，最大變形量上升為0.3877mm。徑向應(yīng)力分布在-206.43MPa～439.03MPa，略微下降，最大位置仍在內(nèi)壁倒角處。環(huán)向應(yīng)力分布均在343MPa左右，略有上升，最大位置仍在螺栓孔處，此處為應(yīng)力集中。

圖8 飛輪轉(zhuǎn)子徑向位移（考慮永磁體作用于內(nèi)壁）

圖9 飛輪轉(zhuǎn)子徑向應(yīng)力分布（考慮永磁體作用于內(nèi)壁）

圖10 飛輪轉(zhuǎn)子環(huán)向應(yīng)力分布（考慮永磁體作用于內(nèi)壁）

3 結(jié)論

飛輪儲能系統(tǒng)工作時(shí)處于高速旋轉(zhuǎn)狀態(tài)，因此飛輪所受離心應(yīng)力不容忽視。本文通過ANSYS靜力學(xué)分析模塊對6000rpm轉(zhuǎn)速下的飛輪轉(zhuǎn)子進(jìn)行應(yīng)力分析。通過對比不考慮永磁體作用于內(nèi)壁和考慮永磁體作用于內(nèi)壁兩種情況，位移變形最大量均集中在飛輪中部，最大變形量為0.3877mm；徑向應(yīng)力分布最大位置在內(nèi)壁倒角處；除去螺栓孔處應(yīng)力集中，環(huán)向應(yīng)力分布均在343MPa左右。在設(shè)計(jì)外轉(zhuǎn)子飛輪系統(tǒng)時(shí)，應(yīng)著重考慮轉(zhuǎn)子中部與內(nèi)壁倒角處的強(qiáng)度問題，選擇合適的轉(zhuǎn)子材料。因此，研究成果對于飛輪儲能系統(tǒng)的設(shè)計(jì)改進(jìn)與優(yōu)化具有一定的指導(dǎo)意義。

[1]張宇,俞國勤,施明融,等.電力儲能技術(shù)應(yīng)用前景分析[J].華東電力,2008,36(4):91-93.

[2]王黎,趙云麗,李衛(wèi)東.飛輪儲能的仿真系統(tǒng)研究[J].電網(wǎng)與清潔能源,2010,26(11):102-106.

[3]孫紅麗,吳士定,張攀峰,等.飛輪儲能系統(tǒng)轉(zhuǎn)子動力學(xué)分析[J]機(jī)械設(shè)計(jì),2016,33(6):72-76.

[4]張秀華,李光喜,王琬.高速儲能飛輪轉(zhuǎn)子的分析與設(shè)計(jì)[J].煤礦機(jī)械,2012,33(6):57-59.

[5]李保軍,王志新,吳定國.飛輪儲能系統(tǒng)充放電過程建模與仿真研究[J].工業(yè)控制計(jì)算機(jī),2011,24(12):107-111.

[6]楊鋒,于飛,張曉鋒,等.飛輪儲能系統(tǒng)建模與仿真研究[J].船電技術(shù),2011,31(4):3-7.

[7]蒙永民,李鐵才.基于飛輪儲能系統(tǒng)的控制方法仿真研究[J].電力電子技術(shù),2008,42(7):37-39.

[8]Mclallin K, Fause J.Advanced energy storage for NASA and USAF missions[R].AFRL/NASA Flywheel Program,2000,10:1- 18.

[9]Thelen RF, Herbst JD, Caprio MT. A 2MW flywheel foe hybrid locomotive power[A].IEEE Vehicular Technology Conference,2003 IEEE 58th Vehicular Technology Conference[C].VTC2003:3231-3235.

[10]Ha Sung K,Kim Seong J,Nasir Sana U et al.Design optimization and fabrication of a hybrid composite flywheel rotor[J].Composite Structures,2012,94(11):3290-3299.

[11]白越,黎海文,吳一輝,等.復(fù)合材料飛輪轉(zhuǎn)子設(shè)計(jì)[J].光學(xué)精密工程,2007,15(6):852-857.

[12]閆曉磊,鐘志華,張義,等.基于最優(yōu)控制理論的儲能飛輪轉(zhuǎn)子形狀優(yōu)化設(shè)計(jì)研究[J].機(jī)械工程學(xué)報(bào),2012,48(3):189-198.