楊旭娟，李宇龍，李兆軍，顏 瑞，張 振

（廣西制造系統(tǒng)與先進(jìn)制造技術(shù)重點(diǎn)實(shí)驗(yàn)室（廣西大學(xué)機(jī)械工程學(xué)院），廣西 南寧530004）

上懸式離心機(jī)是一種間歇式離心機(jī)，它是制糖工藝流程中分蜜工段的重要生產(chǎn)設(shè)備，其工作原理是利用離心力把白砂糖從母液糖膏中分離出來(lái)[1]。上懸式離心機(jī)分離蔗糖晶體和糖膏的能力隨著其轉(zhuǎn)速的提高而增加，而隨著轉(zhuǎn)速的提高，離心機(jī)的振動(dòng)問題就越發(fā)嚴(yán)重，甚至?xí)斐刹槐匾娜藛T傷亡和企業(yè)財(cái)產(chǎn)損失[2]。因此，如何減少離心機(jī)在正常運(yùn)轉(zhuǎn)時(shí)的振動(dòng)已成為一個(gè)亟待解決的難題。目前，有關(guān)離心機(jī)優(yōu)化問題的研究主要集中在離心機(jī)轉(zhuǎn)鼓的優(yōu)化分析、主軸的結(jié)構(gòu)優(yōu)化方面。如韋堯兵，姜永濤等人利用ANSYS軟件對(duì)轉(zhuǎn)鼓進(jìn)行優(yōu)化設(shè)計(jì)，在滿足材料強(qiáng)度和剛度的前提下，減輕了轉(zhuǎn)鼓的質(zhì)量，提高了轉(zhuǎn)鼓的分離性能[3-5]；J.PARK，C.R.WASSGREN以轉(zhuǎn)鼓驅(qū)動(dòng)功率為優(yōu)化目標(biāo)，通過動(dòng)力學(xué)仿真提高了轉(zhuǎn)鼓的節(jié)能性能[6]；李巖舟等人利用ANSYS軟件對(duì)上懸式離心機(jī)轉(zhuǎn)鼓的壁厚進(jìn)行優(yōu)化，降低了轉(zhuǎn)鼓的壁厚，達(dá)到了節(jié)能減排的效果[7]；W.Wilsmann對(duì)離心機(jī)轉(zhuǎn)鼓進(jìn)行應(yīng)力分析，得出轉(zhuǎn)鼓采用薄底板可以減少邊緣應(yīng)力的結(jié)論，對(duì)轉(zhuǎn)鼓的優(yōu)化設(shè)計(jì)具有一定的指導(dǎo)意義[8]；韓玉坤等人利用ANSYS軟件對(duì)離心機(jī)主軸進(jìn)行優(yōu)化分析，提高了主軸的力學(xué)性能，為設(shè)備的升級(jí)和改造提供了參考[9]。以上研究都是關(guān)于離心機(jī)單個(gè)零部件的結(jié)構(gòu)優(yōu)化，而關(guān)于上懸式離心機(jī)主軸系統(tǒng)的整體優(yōu)化分析目前還鮮有報(bào)道，因而有必要對(duì)其進(jìn)行整體優(yōu)化研究。

本文以上懸式離心機(jī)主軸系統(tǒng)為研究對(duì)象，在主軸系統(tǒng)機(jī)電液全局耦合動(dòng)力學(xué)模型的基礎(chǔ)上，以主軸系統(tǒng)輸出端振幅最小為優(yōu)化目標(biāo)，采用遺傳算法對(duì)上懸式離心機(jī)主軸系統(tǒng)的結(jié)構(gòu)參數(shù)進(jìn)行合理的優(yōu)化，以期降低離心機(jī)工作時(shí)的振動(dòng)，提高離心機(jī)正常運(yùn)轉(zhuǎn)時(shí)的穩(wěn)定性。

1 系統(tǒng)機(jī)電液耦合動(dòng)態(tài)方程的建立

研究的上懸式離心機(jī)主要由電機(jī)、主軸、轉(zhuǎn)轂和轉(zhuǎn)鼓等組成，其結(jié)構(gòu)簡(jiǎn)圖如圖1所示。

圖1 上懸式離心機(jī)主軸系統(tǒng)示意圖

不考慮系統(tǒng)的非線性因素，利用有限單元法得到上懸式離心機(jī)主軸系統(tǒng)的機(jī)電液全局耦合動(dòng)力學(xué)方程[10]，該方程反映了系統(tǒng)的結(jié)構(gòu)參數(shù)、電磁參數(shù)和流體參數(shù)與其動(dòng)態(tài)性能之間的關(guān)系。

式中，

U、U觶和U咬分別為系統(tǒng)的廣義位移、廣義速度和廣義加速度列向量；

M為系統(tǒng)質(zhì)量矩陣；

Ma為考慮轉(zhuǎn)鼓與糖膏之間流固耦合作用而引入的質(zhì)量矩陣；

C為系統(tǒng)阻尼矩陣；

K為系統(tǒng)剛度矩陣；

K01考慮系統(tǒng)的機(jī)電耦合關(guān)系而引入的剛度矩陣；

F為系統(tǒng)受到的外加載荷的廣義力列向量；

U咬r為系統(tǒng)的剛體加速度列向量。

2 上懸式離心機(jī)主軸系統(tǒng)的優(yōu)化

2.1 優(yōu)化算法的選擇

在實(shí)際工程運(yùn)用中，許多最優(yōu)化問題由于其系統(tǒng)的復(fù)雜性，很難用傳統(tǒng)的優(yōu)化設(shè)計(jì)方法對(duì)其進(jìn)行求解。傳統(tǒng)優(yōu)化算法通常是采用梯度信息的方法來(lái)尋找最優(yōu)點(diǎn)，容易陷入局部最優(yōu)，難以尋求出真正的最優(yōu)解。遺傳算法是一種全局遍歷的隨機(jī)尋優(yōu)方法[11]，與傳統(tǒng)的優(yōu)化算法相比，其本質(zhì)特征在于群體搜索策略和簡(jiǎn)單的遺傳算子[12]，在計(jì)算的過程中模擬生物在自然環(huán)境中的遺傳和進(jìn)化機(jī)理，反復(fù)將選擇算子、交叉算子、變異算子作用于新的群體，經(jīng)過如此的往復(fù)迭代，最終可以求得問題的最優(yōu)解或者近似最優(yōu)解[13]。該優(yōu)化算法的計(jì)算過程簡(jiǎn)單明了，對(duì)所要求解的優(yōu)化問題沒有太多的數(shù)學(xué)要求，尤其適合用來(lái)處理傳統(tǒng)搜索方法難以解決的高度復(fù)雜的非線性問題[11]，具有較高的使用價(jià)值。

2.2 變量的設(shè)定

將與離心機(jī)主軸系統(tǒng)有關(guān)的結(jié)構(gòu)參數(shù)作為設(shè)計(jì)變量，通過將相關(guān)參數(shù)在適度范圍內(nèi)變動(dòng)取值，來(lái)實(shí)現(xiàn)整個(gè)系統(tǒng)在運(yùn)轉(zhuǎn)時(shí)輸出端振幅值的減小。結(jié)合圖1，與離心機(jī)旋轉(zhuǎn)主軸各段長(zhǎng)度相關(guān)的參數(shù)分別為L(zhǎng)1，L2，L3，其中L1表示聯(lián)軸器中心到軸承1所在軸端面的距離，L2表示軸承1所在軸端面到軸承2所在軸端面的距離，L3表示軸承2所在軸端面到主軸輸出端與轉(zhuǎn)轂交界處的距離；與離心機(jī)旋轉(zhuǎn)主軸各段直徑相關(guān)的參數(shù)為D1，D2，D3，其中D1表示與L1段相對(duì)應(yīng)的主軸外徑，D2表示與L2段相對(duì)應(yīng)的主軸外徑，D3表示與L3段相對(duì)應(yīng)的主軸外徑；與轉(zhuǎn)鼓有關(guān)的參數(shù)為 r，Dlow，Dup，其中 r表示轉(zhuǎn)鼓的內(nèi)徑，Dlow表示轉(zhuǎn)鼓下底板的開口直徑，Dup表示轉(zhuǎn)鼓上頂板開口直徑。則離心機(jī)主軸系統(tǒng)的設(shè)計(jì)變量可以表示為：

2.3 目標(biāo)函數(shù)的設(shè)定

在對(duì)上懸式離心機(jī)主軸系統(tǒng)進(jìn)行優(yōu)化時(shí)，建立一個(gè)合理的數(shù)學(xué)優(yōu)化模型對(duì)于問題的求解十分必要。考慮到離心機(jī)主軸系統(tǒng)的彈性位移列陣u，以及由于離心機(jī)主軸系統(tǒng)在轉(zhuǎn)速發(fā)生變換時(shí)，糖膏往往會(huì)不可避免的對(duì)軸系產(chǎn)生的沖擊，選擇對(duì)離心機(jī)主軸系統(tǒng)速度發(fā)生變換的時(shí)刻進(jìn)行分析，以離心機(jī)主軸系統(tǒng)動(dòng)態(tài)變換后的輸出端振幅的絕對(duì)值最小作為本次優(yōu)化設(shè)計(jì)的目標(biāo)，則離心機(jī)主軸系統(tǒng)的目標(biāo)優(yōu)化函數(shù)可以表示為：

式中，u28和u30皆為主軸系統(tǒng)輸出端節(jié)點(diǎn)6的橫向位移。

2.4 約束條件的設(shè)定

以某型號(hào)的上懸式離心機(jī)主軸系統(tǒng)為研究對(duì)象，已知該型號(hào)的離心機(jī)額定工作轉(zhuǎn)速為1 060 r/min，最大工作轉(zhuǎn)速為1 300 r/min，糖膏的初始密度為1 500 kg/m3，離心機(jī)在工作過程中因布料不均而造成的最大偏心質(zhì)量為1 600 kg，主軸因偏心力造成的變形控制在0.02m以內(nèi)，偏心引起的交變載荷作用在轉(zhuǎn)鼓的質(zhì)心處。與該型號(hào)上懸式離心機(jī)主軸系統(tǒng)有關(guān)的設(shè)計(jì)變量初始值分別為：L1=260 mm，L2=282 mm，L3=1484 mm，D1=130 mm，D2=145 mm，D3=160 mm，r=800 mm，Dlow=530 mm，Dup=570 mm.在此將上懸式離心機(jī)主軸系統(tǒng)參數(shù)變量的約束條件設(shè)定為：

式中，xU、xL分別為x的上限和下限。

最后求出使目標(biāo)函數(shù)g（x）最小時(shí)的相應(yīng)參數(shù)的值。

運(yùn)用matlab遺傳算法工具箱對(duì)上懸式離心機(jī)主軸系統(tǒng)進(jìn)行優(yōu)化，種群類型設(shè)定為為雙精度，種群大小設(shè)定為20，交叉概率設(shè)定為0.8，變異概率設(shè)定為0.01，設(shè)定遺傳算法的終止進(jìn)化代數(shù)100，停滯代數(shù)為100，并給定計(jì)算精度為10-6.則主軸系統(tǒng)的優(yōu)化模型可以表示為：

與離心機(jī)主軸系統(tǒng)相關(guān)的其他參數(shù)的值保持不變。

2.5 優(yōu)化結(jié)果的分析

運(yùn)用matlab遺傳算法工具箱進(jìn)行優(yōu)化，可得主軸系統(tǒng)運(yùn)動(dòng)狀態(tài)發(fā)生變換時(shí)輸出端彈性位移的最小值及其所對(duì)應(yīng)的設(shè)計(jì)變量的值，如表1所示。

表1 初始設(shè)計(jì)與優(yōu)化設(shè)計(jì)后比較

在表1中，變化百分比的正或負(fù)表示在優(yōu)化設(shè)計(jì)后設(shè)計(jì)參數(shù)是增加還是減少。由表1可知，優(yōu)化后，在系統(tǒng)其他零件參數(shù)變化不大的情況下，能實(shí)現(xiàn)上懸式離心機(jī)主軸系統(tǒng)軸端輸出點(diǎn)的最大彈性位移減少8.33%，提高了系統(tǒng)運(yùn)轉(zhuǎn)時(shí)的穩(wěn)定性；優(yōu)化后的系統(tǒng)旋轉(zhuǎn)主軸所受到的糖膏對(duì)其的沖擊力大大減小，進(jìn)而延長(zhǎng)了轉(zhuǎn)軸的使用壽命；優(yōu)化后上懸式離心機(jī)的總重量減少2.1%，節(jié)約了制造成本，提高了經(jīng)濟(jì)效益。

圖2 初始設(shè)計(jì)與優(yōu)化設(shè)計(jì)后主軸系統(tǒng)工作過程輸出點(diǎn)振幅值比較

離心機(jī)主軸系統(tǒng)經(jīng)過優(yōu)化后，可以通過改變?cè)休斎爰?lì)的方法來(lái)改變主軸系統(tǒng)的運(yùn)動(dòng)軌跡以滿足實(shí)際工程任務(wù)的需要，只要離心機(jī)旋轉(zhuǎn)主軸輸入端的轉(zhuǎn)速不超過初始設(shè)計(jì)時(shí)所規(guī)定的輸入端的最大轉(zhuǎn)速，優(yōu)化后的主軸系統(tǒng)在運(yùn)動(dòng)狀態(tài)發(fā)生變換時(shí)輸出端的振幅值不會(huì)有太大變化。

3 結(jié)束語(yǔ)

本文以主軸系統(tǒng)輸出端振幅最小為優(yōu)化目標(biāo)，采用遺傳算法建立離心機(jī)主軸系統(tǒng)的優(yōu)化設(shè)計(jì)數(shù)學(xué)模型，對(duì)系統(tǒng)的結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化。優(yōu)化后的離心機(jī)主軸系統(tǒng)輸出端振幅得到降低，系統(tǒng)的動(dòng)態(tài)性能得到提高，離心機(jī)的整機(jī)重量得以降低，進(jìn)而節(jié)約了其制造成本，提高了企業(yè)的競(jìng)爭(zhēng)力。該優(yōu)化對(duì)上懸式離心機(jī)的設(shè)計(jì)再制造具有一定的借鑒意義。

[1]孫啟才，金鼎五.離心機(jī)原理結(jié)構(gòu)與設(shè)計(jì)計(jì)算[M].北京:機(jī)械工業(yè)出版社，1987.

[2]湯新文，陳海輝.離心機(jī)轉(zhuǎn)鼓爆裂事故分析[J].機(jī)電工程技術(shù)，2005，34（3）:104-105.

[3]韋堯兵，姜永濤，姜明星，等.上懸式離心機(jī)轉(zhuǎn)鼓的應(yīng)力分析與結(jié)構(gòu)優(yōu)化[J].機(jī)械設(shè)計(jì)，2009，26（10）:55-58.

[4]韋堯兵，姜永濤，劉 軍.上懸式離心機(jī)轉(zhuǎn)鼓的動(dòng)態(tài)性能分析與優(yōu)化[J].機(jī)械設(shè)計(jì)，2010，27（08）:69-72.

[5]韋堯兵，姜永濤，剡昌鋒，等.GD—1300上懸式離心機(jī)轉(zhuǎn)鼓的應(yīng)力分析與優(yōu)化[J].機(jī)械設(shè)計(jì)與制造，2009，（06）:220-222.

[6]J.PARK，C.R.WASSGREN.Modeling the Dynamics of Fabric in a Rotating Horizontal Drum Using the Discrete Element Method[J].Particulate Science and technology，2003，（21）:157-175.

[7]李巖舟，陳 淵，王 偉，等.上懸式制糖離心機(jī)轉(zhuǎn)鼓壁厚優(yōu)化[J].機(jī)械設(shè)計(jì)與制造，2013，（09）:178-180.

[8]馮 卉.臥式螺旋沉降卸料離心機(jī)動(dòng)力學(xué)仿真與結(jié)構(gòu)優(yōu)化[D].南京:東南大學(xué)，2011.

[9]韓玉坤，段 非，姜永濤.XZ-1400上懸式離心機(jī)主軸的應(yīng)力分析與結(jié)構(gòu)優(yōu)化[J].食品與機(jī)械，2011，27（04）:98-100.

[10]徐光恒.上懸式離心機(jī)主軸系統(tǒng)動(dòng)力學(xué)特性研究[D].南寧:廣西大學(xué)，2015.

[11]玄光男，程潤(rùn)偉.遺傳算法與工程優(yōu)化[M].北京:清華大學(xué)出版社，2004.

[12]馬永杰，云文霞.遺傳算法研究進(jìn)展[J].計(jì)算機(jī)應(yīng)用研究，2012，（04）:1201-1206.

[13]陳輝慶.新型可控變胞式碼垛機(jī)器人機(jī)構(gòu)動(dòng)態(tài)穩(wěn)定性研究[D].南寧:廣西大學(xué)，2015.