宋志強(qiáng)，史青錄，彭萬(wàn)萬(wàn)，陳貫祥

（1.太原科技大學(xué) 機(jī)械工程學(xué)院，山西 太原 030024；2.中北大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院，山西 太原 030051）

0 引言

電機(jī)作為旋轉(zhuǎn)機(jī)械設(shè)備，在現(xiàn)代生產(chǎn)中扮演著不可或缺的角色，但電機(jī)振動(dòng)時(shí)有發(fā)生。振動(dòng)問(wèn)題往往是多種復(fù)雜因素綜合造成的，其中電機(jī)和底座之間發(fā)生共振就是其表現(xiàn)形式之一。當(dāng)轉(zhuǎn)子達(dá)到一定轉(zhuǎn)速時(shí)，轉(zhuǎn)軸頻率和底座固有頻率重合就會(huì)發(fā)生共振，共振將在電機(jī)上產(chǎn)生動(dòng)載荷，如果在該轉(zhuǎn)速下繼續(xù)運(yùn)行下去，必將加速轉(zhuǎn)子、軸承等部件的磨損，從而導(dǎo)致疲勞失效。通過(guò)模態(tài)分析，可以大致了解底座結(jié)構(gòu)固有頻率的分布，經(jīng)過(guò)設(shè)計(jì)改進(jìn)，使電機(jī)的工作轉(zhuǎn)速盡量避開(kāi)系統(tǒng)的共振頻率區(qū)域，以提高電機(jī)的使用壽命。

1 有限元前處理流程

建立有限元模型是整個(gè)分析過(guò)程的關(guān)鍵，合理的模型既能保證計(jì)算結(jié)果的精度，又不致使計(jì)算量太大和對(duì)計(jì)算機(jī)存儲(chǔ)容量的要求太高。有限元前處理流程如圖1所示。

2 建立底座有限元模型

2.1 單元類(lèi)型選擇

根據(jù)底座的形狀特征和精度要求，模型全部采用二階四面體單元，主要考慮以下兩點(diǎn)：一是高階單元的曲線或曲面邊界能更好地逼近結(jié)構(gòu)的邊界曲線或曲面；二是高階單元的高次位移函數(shù)能更好地逼近結(jié)構(gòu)復(fù)雜的位移分布［1］。

2.2 單元屬性定義

有限元中的每一個(gè)單元除了表現(xiàn)出一定的外部形狀外，還應(yīng)具備一組計(jì)算所需的內(nèi)部特征參數(shù)，這些參數(shù)用來(lái)定義結(jié)構(gòu)材料的性能、描述單元本身的物理特征和其他輔助幾何特征等。該模型屬性參數(shù)為：?jiǎn)卧?lèi)型Psolid，材料為Steel，彈性模量為210 GPa，泊松比為0.3，密度為7 900 kg／m3。

圖1 有限元前處理流程圖

2.3 網(wǎng)格質(zhì)量檢查

網(wǎng)格質(zhì)量是指網(wǎng)格幾何形狀的合理性。有些網(wǎng)格形狀是不允許的，它們會(huì)導(dǎo)致單元?jiǎng)偠染仃嚍榱慊蜇?fù)值，使有限元計(jì)算出現(xiàn)致命錯(cuò)誤，這種網(wǎng)格稱(chēng)為畸形網(wǎng)格。底座結(jié)構(gòu)有限元網(wǎng)格單元質(zhì)量控制為：縱橫比（aspect）的范圍為5～6；歪扭角（skew）的范圍為60～75；翹曲（warpage）的范圍為30～40；雅可比（jacobian）的范圍為0.5～0.6；三角形（tria）的范圍為20～120；坍塌（tet collapse）的范圍為0.15～0.5。

2.4 模型簡(jiǎn)化和施加約束

為了操作可行性，對(duì)模型進(jìn)行簡(jiǎn)化處理。由于電機(jī)部分彈性體模態(tài)不在考慮范圍內(nèi)，故用集中質(zhì)量單元（CONM2）代替，并用rbe3單元固定在底座支架上。用剛性單元rbe2模擬電機(jī)與底座以及底座與地面的螺栓連接，并約束其6個(gè)自由度。底座結(jié)構(gòu)有限元模型如圖2所示。

圖2 底座結(jié)構(gòu)有限元模型

3 底座結(jié)構(gòu)模態(tài)分析

3.1 模態(tài)分析理論及計(jì)算方法

模態(tài)分析是確定設(shè)計(jì)結(jié)構(gòu)或機(jī)械零部件的振動(dòng)特性并得到結(jié)構(gòu)固有頻率和振型的過(guò)程［2］。有阻尼的N自由度系統(tǒng)的強(qiáng)迫振動(dòng)運(yùn)動(dòng)微分方程可以表示為：

其中：［M］，［C］和［K］分別為系統(tǒng)質(zhì)量、阻尼和剛度矩陣；｛x″｝，｛x′｝，｛x｝和｛F｝分別為系統(tǒng)各點(diǎn)加速度列陣、速度列陣、位移列陣和激勵(lì)力列陣。模態(tài)分析時(shí)，要求解的是固有頻率和振型參數(shù)，與外載荷無(wú)關(guān)，即｛F｝＝｛0｝，則得到系統(tǒng)的自由振動(dòng)方程。在計(jì)算系統(tǒng)的固有頻率和振型時(shí)，結(jié)構(gòu)阻尼的影響很小，可以忽略不計(jì)，這時(shí)無(wú)阻尼自由振動(dòng)的微分方程為：

其對(duì)應(yīng)的特征方程為：

其中：ωi為系統(tǒng)的第i階固有頻率。N自由度的系統(tǒng)就有N 個(gè)固有頻率［3］。

3.2 模態(tài)分析結(jié)果

在Hypermesh中創(chuàng)建模態(tài)分析載荷步，利用求解器Radioss進(jìn)行求解。本模型為有約束的模態(tài)分析。模態(tài)的振型階次越高，與該模態(tài)頻率對(duì)應(yīng)的變形就越小，由于高階模態(tài)對(duì)結(jié)構(gòu)的動(dòng)力學(xué)影響很小，并考慮到外部激勵(lì)多處于低頻區(qū)域［4］，故底座結(jié)構(gòu)一般關(guān)心較低的幾階頻率。所求得的底座結(jié)構(gòu)前6階固有頻率見(jiàn)表1，各階振型如圖3所示。

在模態(tài)分析結(jié)果中，振型的大小只是一個(gè)相對(duì)的量值（位移相對(duì)值），它表征的是各點(diǎn)在某一階固有頻率上振動(dòng)量值的相對(duì)比值，反映該固有頻率上振動(dòng)的傳遞情況，并不反映實(shí)際振動(dòng)的數(shù)值［5］。在振型圖中，前兩階振型為沿某一坐標(biāo)軸的來(lái)回?cái)[動(dòng)；第3階為彎曲變形；第4和第5階為彎扭組合；第6階為局部模態(tài)，指模態(tài)振型在結(jié)構(gòu)的大部分自由度上接近零（模型大部分保持原狀，只在局部發(fā)生較大變形），僅在少部分自由度上存在明顯的幅值。局部模態(tài)對(duì)結(jié)構(gòu)響應(yīng)的貢獻(xiàn)幾乎為零，因此，對(duì)于結(jié)構(gòu)整體動(dòng)特性和響應(yīng)分析，可以忽略局部模態(tài)的影響。

表1 底座的前6階固有頻率和振型

圖3 底座的前6階模態(tài)振型圖

4 模型固有頻率優(yōu)化

提高結(jié)構(gòu)剛度是改善結(jié)構(gòu)振動(dòng)性能的有效措施，通過(guò)對(duì)結(jié)構(gòu)固有頻率進(jìn)行控制和調(diào)整，能有效降低結(jié)構(gòu)的動(dòng)力響應(yīng)。對(duì)于多自由度系統(tǒng)，避免共振的措施就是不使系統(tǒng)受到擾頻等于或接近于系統(tǒng)任一階固有頻率的外力作用［6］。為避免電機(jī)系統(tǒng)發(fā)生共振，提高其底座的固有頻率是行之有效的方法。究竟要提高哪一階模態(tài)的固有頻率，要視電機(jī)實(shí)際的工作轉(zhuǎn)速而定。本文以提高第1階固有頻率為例，闡述優(yōu)化底座結(jié)構(gòu)固有頻率的方法。

4.1 拓?fù)鋬?yōu)化簡(jiǎn)介

結(jié)構(gòu)拓?fù)鋬?yōu)化是一個(gè)在給定設(shè)計(jì)變量和約束條件下求解目標(biāo)函數(shù)最優(yōu)化的過(guò)程。目標(biāo)函數(shù)、設(shè)計(jì)變量、約束函數(shù)為優(yōu)化問(wèn)題的三要素，其主要思想是在給定的設(shè)計(jì)區(qū)域內(nèi)尋求材料的最優(yōu)分布問(wèn)題。經(jīng)過(guò)拓?fù)鋬?yōu)化，可以尋找到結(jié)構(gòu)合理的載荷傳遞路徑，將材料分布在最需要加強(qiáng)的區(qū)域，極大限度地提高材料的利用效率，使結(jié)構(gòu)的模態(tài)、剛度等性能大幅提升［7］。

4.2 確定設(shè)計(jì)變量、約束條件及目標(biāo)

這里采用變密度法，以1階模態(tài)作為對(duì)電機(jī)系統(tǒng)振動(dòng)影響最大的模態(tài)，因此把提高1階固有頻率設(shè)為優(yōu)化設(shè)計(jì)的目標(biāo)。設(shè)計(jì)變量（Design variables）為組成模型的所有四面體單元的密度值，根據(jù)經(jīng)驗(yàn)，約束響應(yīng)（Dconstraints）設(shè)置為模型體積分?jǐn)?shù)（Volfrac）上限的30%（即優(yōu)化后體積不超過(guò)原底座結(jié)構(gòu)體積30%作為約束），目標(biāo)（Objective）設(shè)置為 Max，模態(tài)階數(shù)（Mode number）為1。

4.3 提交計(jì)算

設(shè)置結(jié)果保存路徑后，運(yùn)行Optistruct進(jìn)行求解。經(jīng)過(guò)21次迭代后，得到最后一次迭代的單元密度云圖，見(jiàn)圖4。圖4中大部分區(qū)域的單元密度收斂于1或0，且最大密度值為1，說(shuō)明優(yōu)化迭代已進(jìn)行徹底。頻率與剛度及質(zhì)量相關(guān)，提高剛度成為提高頻率最直接的方法。紅色區(qū)域（A處）表示需要加強(qiáng)的區(qū)域，可以適當(dāng)增加厚度、加肋板或改變材質(zhì)。

圖4 材料最優(yōu)分布云圖

4.4 查看頻率迭代變化曲線和迭代數(shù)據(jù)

在結(jié)果文件中，打開(kāi)base＿h(yuǎn)ist.mvw文件，可以觀察頻率隨迭代次數(shù)變化的情況，從中可以計(jì)算出頻率增加的比例。圖5為1階頻率迭代變化曲線，表2為首、末兩次迭代數(shù)據(jù)。

圖5 1階頻率迭代變化曲線

由表2可以看出：經(jīng)過(guò)優(yōu)化處理，第1階固有頻率增加了79.82%（迭代數(shù)據(jù)表示頻率變化率，不同于模態(tài)振型圖中的固有頻率），由此可計(jì)算出優(yōu)化后的第1階固有頻率增加為380.72 Hz；第6階頻率在迭代過(guò)程中出現(xiàn)降低的情況，即在優(yōu)化過(guò)程中會(huì)優(yōu)先考慮某些因素，這是允許的。由此可知，只要處理得當(dāng)，可以明顯提高底座結(jié)構(gòu)的第1階固有頻率，使其避開(kāi)電機(jī)的工作頻率，從而有效控制共振的發(fā)生。

表2 首、末兩次迭代數(shù)據(jù)

5 結(jié)論

（1）通過(guò)模態(tài)分析，得到了電機(jī)底座低階振動(dòng)固有特性，并可以動(dòng)態(tài)、直觀地觀察每階振型的變化。

（2）優(yōu)化后第1階固有頻率提高了79.82%，其他幾階固有頻率也有不同程度的提高，避開(kāi)了可能發(fā)生共振的頻率。

（3）優(yōu)化結(jié)果顯示了滿(mǎn)足設(shè)計(jì)要求的材料最優(yōu)分布，為電機(jī)底座的最優(yōu)設(shè)計(jì)提供了參考，體現(xiàn)了拓?fù)鋬?yōu)化方法在機(jī)械零部件初始設(shè)計(jì)過(guò)程中的應(yīng)用價(jià)值。

［1］ 曾攀.有限元分析及應(yīng)用［M］.北京：清華大學(xué)出版社，2003.

［2］ 李瑞虎，王凱，謝國(guó)燦.基于ANSYS的平面凸輪有限元模態(tài)分析［J］.機(jī)械傳動(dòng)，2006，30（6）：61－63.

［3］ 翟文波.基于ANSYS技術(shù)的直線振動(dòng)電機(jī)模態(tài)分析［J］.微電機(jī)，2011，44（4）：29－31.

［4］ 宋桂霞.某大客車(chē)車(chē)架結(jié)構(gòu)模態(tài)分析［J］.計(jì)算機(jī)輔助工程，2011，20（2）：24－27.

［5］ 袁安富，陳俊.ANSYS在模態(tài)分析中的應(yīng)用［J］.設(shè)計(jì)與研究，2007（8）：79－82.

［6］ 侯之超，鄭兆昌.關(guān)于線性系統(tǒng)避免共振的條件［J］.振動(dòng)與沖擊，1997，16（3）：15－20.

［7］ 張鵬飛，董瑞強(qiáng).拓?fù)鋬?yōu)化在白車(chē)身概念設(shè)計(jì)中的應(yīng)用［J］.汽車(chē)技術(shù)，2010（7）：55－58.