基于有限元法和多島遺傳算法的飛輪結構參數優化設計

丁泉惠，王　森，黃修長，華宏星( .上海交通大學振動沖擊噪聲研究所機械系統與振動國家重點實驗室，上海0040；.上海航天控制技術研究所，上海0033）

?

基于有限元法和多島遺傳算法的飛輪結構參數優化設計

丁泉惠1，王森2，黃修長1，華宏星1
( 1.上海交通大學振動沖擊噪聲研究所機械系統與振動國家重點實驗室，上海200240；2.上海航天控制技術研究所，上海200233）

摘要：針對一種額定角動量為68 Nms的飛輪，以飛輪輪緣響應相對于基礎激勵輸入的加速度傳遞率為優化目標，以輪體轉動慣量、幾何尺寸、強度、剛度、質量等為約束條件，對截面形狀為工字型的輪輻和輪緣幾何尺寸參數的設計提出了優化方法。該優化方法基于有限元法建立飛輪輪體模型，同時利用多島遺傳算法對飛輪輪體進行了多學科設計優化研究。結果表明，在滿足設計要求的前提下，輪緣加速度傳遞率得到明顯降低。該方法對提高航天器飛輪輪體結構設計效率具有一定的現實意義。

關鍵詞：振動與波；飛輪；有限元法；多島遺傳算法；優化設計

近年來各國衛星數量增多，應用更加廣泛，表明衛星在各國經濟、社會發展和國防建設中發揮著越來越重要的作用。

衛星通過姿態控制系統控制自身的姿態，姿態控制方法可分為被動式和主動式兩類。衛星對姿態進行主動控制的系統由姿態敏感器、控制器和執行機構三部分組成。在衛星三個軸各設一個執行機構，就可以保持衛星主動三軸穩定，而飛輪是主動三軸穩定姿態控制系統的執行機構之一[1]，飛輪系統利用飛輪所儲存的角動量與飛行器進行角動量交換，來實現穩定飛行的三軸姿態，已成為衛星主動三軸穩定姿態控制系統中一種廣泛應用的機構執行方式[2]。

隨著中國航天事業不斷發展，對衛星姿態控制系統的可靠性、壽命、穩定性和精度等要求越來越高。在可靠性、壽命方面，在航天器發射階段振動、沖擊、加速度過載以及多種動態載荷作用[3]等嚴酷的動力學環境下，輪體結構產生較大的撓性形變，要保證輪體不會損傷，如果存在的振動慣性力不超過飛輪軸承組件的動負荷能力，軸承組件就不會損傷[4]。在穩定性、精度方面，飛輪是高精度衛星姿態控制的主要干擾源之一[5–7]，飛輪擾動力會引起航天器的微振動，而飛輪輪體結構撓性是飛輪擾動原因之一，可能與轉子不平衡引起的高頻振動激勵產生諧振，從而出現結構彈性振動[8–10]，如果由不平衡特性產生高頻激勵激發出飛輪的固有模態，振動將會更加劇烈，衛星的姿態控制和穩定將會受到更嚴重的影響。所以減小飛輪輪體結構撓性形變很有必要，這對衛星姿態控制系統的可靠性和壽命，甚至對其精度和穩定性具有重要意義。

目前飛輪輪體優化的主要工作是得到盡可能大的轉動慣量/質量比，文獻[11，12]提出飛輪在額定轉速達到一定角動量同時滿足結構剛度、強度、形狀等要求下，使轉子質量最小的優化方案。

文中提出一種優化方法，是在滿足輪體的轉動慣量、幾何約束、剛度約束、強度約束等多種約束條件下，以飛輪輪體輪緣響應為優化目標來減小輪體結構撓性形變。

1　飛輪輪體建模及簡化

飛輪是由飛輪（本體）和飛輪驅動與控制驅動組成，通常這兩部分結構上彼此獨立。而飛輪（本體）是由殼體組件、軸承組件、電機組件和輪體組件組成，如圖1所示。

圖1　飛輪組成

殼體組件是由底座、氣密外殼以及支撐上下外殼的芯軸和球冠組成，主要功能是保持真空密封環境，完成熱量傳導，提供機電接口和抵御外界機械沖擊等。

軸承組件是飛輪中關鍵的支撐部件，它包括主軸、軸承、預緊系統、潤滑供油系統、密封結構以及緩沖彈簧構件。

電機組件是電能轉換為機械能的機構，電機由永磁體的轉子和無鐵心電樞的定子構成，位置傳感器埋置在定子中，能可靠地反應轉子的位置。

輪體組件是飛輪系統的關鍵部件之一，常見的結構有兩種：盤式和輻條式結構。對于轉動慣量為小容量的飛輪，一般采取盤式結構的輪體。對于轉動慣量為大容量的飛輪，則采用輻條式結構輪體，以獲得比較好的慣量/質量比，降低飛輪總體質量。圖2為輻條式結構輪體剖視圖。輪體由輪轂、輪輻和輪緣三部分組成，輪輻一般采用斷面為工字型的薄壁結構，以減輕重量，并且具備足夠大的剛度和強度。同時為了提高軸向固有頻率，輪輻具有約9°左右的傾角。

圖2　輻條式結構輪體剖視圖

1.1建立模型

APDL是Ansys的參數化設計語言，可以實現有限元分析全過程的參數化，修改其中的參數可進行反復分析，只有創建參數化分析，才能對參數進行優化設計，運用APDL語言采用自頂向下建模方法建立有限元模型，建模時由關鍵點定義線、面。建立飛輪有限元模型時采用梁單元和殼單元，其中輪緣和輪輻采用梁單元，而輪轂采用殼單元建模。飛輪有限元模型如圖3所示，輪轂底部施加軸線Z向約束，輪轂施加徑向約束。

圖3　飛輪有限元模型

1.2有限元分析

1.2.1模態分析

根據約束邊界條件，進行模態分析，提取前4階模態。模態振型以及振型描述如圖4—圖7所示。

1.2.2諧響應分析

Ansys有限元軟件只能施加位移載荷和力載荷，無法在局部施加加速度載荷，而且施加非零位移時只能使用完全法計算諧響應，這樣不僅耗費計算時間，而且延長了分析時間。所以飛輪采用基于模態疊加法和大質量法計算諧響應，計算效率遠高于完全法，有利于飛輪輪體結構優化設計。一般取質量點的質量為飛輪結構質量的104～0.9×108倍。如圖8所示，底部模擬振動臺基礎，在底部質點上施加力載荷，使之產生所需大小的加速度載荷，從而實現加速度載荷的施加。

圖4　第1階振型（沿Z軸局部彎曲振型）

圖5　第2階振型（沿Z軸垂向振型）

圖6　第3階振型（沿Z軸局部扭轉振型）

圖7　第4階振型（沿XZ平面彎曲振型）

圖8　飛輪諧響應分析

定義飛輪輪緣響應相對于基礎激勵輸入的加速度傳遞率為飛輪輪緣響應加速度a1相對于基礎激勵輸入加速度a2的比值a1/a2，單位為無量綱；同時圖中加速度傳遞率均指飛輪輪緣響應相對于基礎激勵輸入的加速度傳遞率。圖9所示為諧響應分析結果。

在350 Hz左右頻率處，只有一個峰值，并且由模態分析知，此頻率下振型為第2階沿軸向“傘狀”振型。所以此階振型主要影響輪體結構撓性形變的整體振型，輪緣響應在這階頻率附近響應最大，故將其作為主要抑制的模態。

圖9　諧響應結果分析

2　優化方法及優化模型建立

2.1多島遺傳優化算法

遺傳算法是模擬生物在自然界中遺傳和進化過程而形成的一種自適應全局優化概率搜索算法，其進化過程如圖10所示。遺傳算法所有操作在每個島上進行，每個島上選定的個體定期地遷移到另外島上，然后繼續進行傳統遺傳算法操作（選擇、雜交、變異操作）。傳統的遺傳算法在優化過程中基因突變的概率較低，容易在進化幾代后出現早熟現象，導致優化結果收斂于局部最優解。

圖10　傳統遺傳算法

多島遺傳算法是在傳統遺傳算法基礎上建立的一種基于群體分組的并行性遺傳算法，其進化過程如圖11所示。多島遺傳算法將整個種群分解為多個子群（這些子群稱為“島”），并將各個子群互相隔絕于不同的“島嶼”上，對每個子群中的個體進行傳統遺傳算法操作（選擇、雜交、變異操作），各個子群獨立地進化，而不是全部種群采用相同的進化機制，并且各個“島嶼”間以一定的時間間隔進行“遷移”操作，使各個“島嶼”間進行信息交換．并且每隔幾代挑選子群中的個體進行交換（遷移操作），保證了進化過程中優化解的多樣性，從而有效抑制了早熟現象的發生，有利于找到全局最優解。

圖11　多島遺傳算法

從飛輪輪體結構優化設計角度來說，飛輪優化設計是包含有多種約束的非線性優化問題，其目標函數以及部分約束條件與結構參數之間的關系復雜，這里采用全局優化算法中的多島遺傳算法進行優化設計，邏輯結構簡單，直觀性強，易于程序化。

2.2優化模型建立

2.2.1優化目標

對飛輪進行結構優化設計時，要求飛輪在額定轉速達到一定角動量，同時在滿足幾何尺寸、剛度、強度、質量等要求下，以飛輪輪緣響應相對于基礎激勵輸入的加速度傳遞率為優化目標。

設計變量

式(1)表示設計變量x是由變量h2、r2、t2、b1、b2、b3、b4、b5、b6組成的集合，x為設計變量，如圖12所示，h2為輪緣高度，r2為輪緣半徑，t2為輪緣厚度，b1和b2為上截面和下截面寬度相對于b6的增量，b3為梁高度相比于b4+b5的增量，b4和b5為上截面和下截面的厚度；b6為支撐梁的厚度。

圖12　飛輪輪體輪緣和輪輻截面標注圖

目標優化函數為

式(2)為目標優化函數，即飛輪輪緣響應相對于基礎激勵輸入的加速度傳遞率取最小值的設計結果為最優；

約束條件為

約束條件可以用等式約束、不等式約束和邊界約束表示。式(3)為等式約束，式(4)為不等式約束，式(5)為邊界約束。

2.2.2約束條件

(1)飛輪輪體轉動慣量約束。飛輪在額定轉速65 00 r/min時需要提供68 Nms的角動量，要求轉動慣量I不大于0.098±0.002kg?m2。

(2)飛輪輪體幾何約束。為限制系統的總體體積和質量，給定的飛輪幾何約束如表1所示。

表1　飛輪幾何約束

(3)剛度約束。為了滿足飛輪輪體在工作轉速范圍內為剛性轉子的要求，需要其彈性1階共振頻率大于轉子最高工作轉速的1.4倍。對于最高轉速7 500 r/min的飛輪輪體，同時考慮控制系統的需要，則要求飛輪輪體的彈性1階共振頻率f1>175 Hz。

(4)強度約束。考慮飛輪旋轉實際工況，在安全系數2情況下，需要滿足飛輪在最高轉速下的最大Mises等效應力σmax,eq≤[σ]/2。

(5)強制約束。考慮飛輪與軸承組件、外殼等部件需要裝配的實際條件，軸承組件和外殼等尺寸已經確定，所以輪轂內徑和軸向方向尺寸為強制約束。

3　飛輪優化算例

飛輪輪體優化流程如圖13所示，首先利用有限元法建立飛輪輪體靜力學、動力學有限元模型，根據初始參數計算得到初步的分析結果，利用Isight基于多島遺傳優化算法，確定以輪緣響應相對于基礎輸入的加速度傳遞率為優化目標和以h2、r2、t2、b1、b2、b3、b4、b5、b6為優化設計變量，根據飛輪輪體的轉動慣量、幾何、剛度、強度等約束條件，對飛輪輪體結構進行優化設計。優化設計所使用的遺傳算法控制參數中種群數、交叉概率、變異概率最大代數分別為100、0.8、0.05和10。

算例中優化結果如表2所示。由表2可知，優化結果滿足設計要求，同時輪體材料為高強度鋼，強度σb≥1 000 Mpa，計算結果為219 Mpa，安全系數大于2。圖14所示為對飛輪優化前后輪緣響應相對于基礎激勵輸入的加速度傳遞率曲線的比較。優化后可見共振峰朝著低頻方向移動，從345 Hz變為301 Hz，共振峰的幅值也有所減小，從12 155.024變為1 458.980，優化后的質量為4.470 kg，優化前的質量為4.587 kg，優化后結果表明振動傳遞得到一定的控制，同時減小飛輪的質量，優化結果令人滿意。

圖13　飛輪輪體優化設計流程

圖14　加速度傳遞率-頻率曲線

4　結語

通過研究，可以得到以下結論：

(1)對輻條式輪體結構進行有限元分析，并基于有限元算法和多島遺傳算法，在飛輪額定轉速和同時滿足輪體的轉動慣量、幾何、剛度、強度等約束條件情況下，以工字梁和輪緣為優化參數，以輪緣響應相對于基礎激勵輸入的加速度傳遞率為優化目標，對輪體進行結構優化設計，振動得到了抑制，同時減小了飛輪的質量。

表2　飛輪優化設計結果

(2)利用有限元建立參數化模型對提高飛輪產品設計和研制的效率以及降低成本有很重要的意義。

(3)優化思路和方法可以應用于飛輪系列產品的設計和研制。

參考文獻：

[1]徐矛，康建勇.航天科技基礎[M].北京：國防工業出版社，2008.

[2]程顥.飛輪控制系統研究與設計[D].哈爾濱：哈爾濱工業大學，2009.

[3]王建偉，王剛，溫激鴻，等.航天器發射階段聲振環境載荷控制技術研究進展[J].噪聲與振動控制，2011，31(1)：1-8.

[4]于登云，楊建中.航天器機構技術[M].北京：中國科學技術出版社，2011.

[5] Bialke B. A compilation of reaction wheel induced spacecraft disturbances[C]//20 th Annual American Aeronautical Society Guidance and Control Conference. 1997:97-038.

[6] Neat G W, Melody J W, Lurie B J. Vibration attenuation approach for spaceborne optical interferometers[J]. Control Systems Technology, IEEE Transactions on, 1998, 6(6):689-700.

[7]李太平，姚煜中，李霖圣，等.不平衡飛輪的擾動特性[J].噪聲與振動控制，2014，34(3)：34-41.

[8]白爭鋒，趙陽，馬文來，等.飛輪擾動對航天器結構動態特性的影響分析[J].宇航學報，2009，30(5)：2073-2079.

[9] Masterson R A, Miller D W, Grogan R. Development of empirical and analytical reaction wheel disturbancemodels [J].AIAAPaper, 1999:99-1204.

[10]李連軍，戴金海.飛輪系統內干擾建模與仿真分析[J].系統仿真學報，2005，17(8)：1855-1858.

[11]韓邦成，虎剛，房建成，等. 50 Nms磁懸浮反作用飛輪轉子優化設計方法的研究[J].宇航學報，2006，3：536-540.

[12]趙麗濱，趙友選，張建宇，等.反作用飛輪結構的動態優化設計[J].中國機械工程，2007，(2)：139.

Optimization Design of Flywheel Structural Parameters Based on Finite Element Analysis and MIGAMethod

DING Quan-hui1, WANG Sen2, HUANG Xiu-chang1, HUA Hong-xing1
( 1. Instituteof Vibration, Shock and Noise, StateKey Laboratory of Mechanical Systemand Vibration, Shanghai Jiaotong University, Shanghai 200240, China; 2. Shanghai AerospaceControl Technology Research Institute, Shanghai 200233, China)

Abstract:An optimization design method for thegeometrical parametersof wheel spoke, consisting of I-shaped cross section and wheel flange, wasproposed for flywheels. A flywheel with 68Nmsnominal angular momentum wastaken asan example. In theoptimization, theconstraint conditionswerethemoment of inertiaof thewheel body, geometrical parameters and strength, stiffness and mass of the flywheel. And the goal of the optimization was to minimize the acceleration transmissibility (defined as the acceleration rate) of the rim output response to the input of base-excitation. The finite element method wasused to develop themodel and the Multi-Island Genetic Algorithm (MIGA) wasemployed to carry out the multi- disciplinary design optimization. It is demonstrated that by using the proposed method, the acceleration transmissibility is reduced significantly on the premise of meeting the design requirements. The proposed method is able to improvethedesignefficiency for spacecraft flywheels.

Key words:vibrationandwave; flywheel; finiteelement method; MIGA; optimizationdesign

通訊作者：華宏星，教授，博士生導師。

作者簡介：丁泉惠(1984- )，男，黑龍江人，碩士生，研究方向為振動分析與控制。E-mail:dingquanhui@163.com

收稿日期：2015-10-26

文章編號：1006-1355(2016)02-0056-05

中圖分類號：O328；O329

文獻標識碼：ADOI編碼：10.3969/j.issn.1006-1335.2016.02.012