安曉衛(wèi)，穆 珊，徐文彬，李秀艷，李素妍

(1.沈陽理工大學(xué) 機械工程學(xué)院，遼寧 沈陽 110159；2.鞍山重型礦山機械股份有限公司，遼寧 鞍山 114051)

振動篩是篩分機械中的重要設(shè)備，大處理量、高效率是其發(fā)展趨勢。振動篩在惡劣的工作條件下常常由于激勵引起的共振使側(cè)板斷裂、篩框扭曲變形等故障，對振動篩的結(jié)構(gòu)布局進(jìn)行優(yōu)化設(shè)計，提高其結(jié)構(gòu)剛度，不僅能延長振動篩的使用壽命，還可以帶來巨大的經(jīng)濟效益。

文獻(xiàn)[1]用Solidworks軟件建立篩體的三維實體模型，對圓振動篩進(jìn)行模態(tài)分析，得到篩箱的固有頻率和振型；文獻(xiàn)[2]對大型直線振動篩進(jìn)行了動態(tài)仿真和仿真測試，并根據(jù)計算結(jié)果對大型篩體做了改進(jìn)設(shè)計；文獻(xiàn)[3]以ZK450直線振動篩為研究對象，運用有限元分析軟件ANSYS進(jìn)行動態(tài)特性分析，計算出固有頻率和模態(tài)振型及在額定載荷下的應(yīng)力分布情況，為以后振動篩的設(shè)計提供了參考和借鑒；文獻(xiàn)[4]對振動篩在ANSYS中進(jìn)行強度仿真，得到振動篩的強度規(guī)律，從試驗結(jié)果出發(fā)，對振動篩進(jìn)行壽命估計；文獻(xiàn)[5]嵌入解析靈敏度法將板厚和截面寬度作為設(shè)計變量，以減小側(cè)板質(zhì)量為目的，對振動篩進(jìn)行優(yōu)化設(shè)計；文獻(xiàn)[6]以篩箱各板厚為設(shè)計變量，以總質(zhì)量為目標(biāo)函數(shù)對篩箱結(jié)構(gòu)進(jìn)行優(yōu)化計算；文獻(xiàn)[7]對大型圓振動篩進(jìn)行動力學(xué)分析和動態(tài)優(yōu)化設(shè)計。但是文獻(xiàn)[5-7]只將某些板厚作為設(shè)計變量對振動篩進(jìn)行優(yōu)化，并沒有考慮各構(gòu)件位置布局對其固有特性的影響；文獻(xiàn)[8]中為使固有頻率遠(yuǎn)離激勵頻率，在篩框底部增加了兩條縱向矩形加強梁以提高整體結(jié)構(gòu)的縱向剛度。文獻(xiàn)[1-8]對振動篩的修改只是憑設(shè)計者的經(jīng)驗，缺少一定的理論依據(jù)；結(jié)構(gòu)優(yōu)化設(shè)計大多停留在尺寸優(yōu)化層次上，對于振動篩這樣大型復(fù)雜的工程機械進(jìn)行整機優(yōu)化并涉及靈敏度這方面的研究資料尚少[9]。而關(guān)于振動篩結(jié)構(gòu)形狀優(yōu)化方面的研究幾乎是空白。

振動篩結(jié)構(gòu)復(fù)雜，影響結(jié)構(gòu)動態(tài)特性的因素多。為避免結(jié)構(gòu)修改的盲目性，提高優(yōu)化設(shè)計的計算效率，減少計算成本，本文通過靈敏度分析，尋找出結(jié)構(gòu)的薄弱環(huán)節(jié)，以對結(jié)構(gòu)特性影響大的參數(shù)作為設(shè)計變量，對某振動篩的結(jié)構(gòu)布局進(jìn)行形狀優(yōu)化設(shè)計。

1 振動篩模型處理及模態(tài)分析

1.1 有限元模型處理

振動篩的外形尺寸長、寬、高分別為9.117m、4.9m和3.49m。其結(jié)構(gòu)主要是側(cè)板、上下橫梁、主梁、前后彈簧座耳、U型加強板、后擋板和支撐角鋼等組成，其材料為Q235,側(cè)板是一塊整鋼板，上、下橫梁、主梁、U型加強板和后擋板是由鋼板焊接而成，彈簧座耳采用鉚焊結(jié)構(gòu)，用兩塊端板與側(cè)板鉚釘連接到一起。振動篩體積大，連接件多，結(jié)構(gòu)復(fù)雜。振動篩結(jié)構(gòu)上有一些工藝小孔，由于其僅對很小的局部區(qū)域的剛度有影響，對整個振動篩的計算精度影響很小，因此，綜合考慮其結(jié)構(gòu)形狀、承載方式、邊界條件、計算成本與可行性，在滿足工程精度的前提下，建模時將這些小孔忽略不計。采用板單元模擬側(cè)板、上下橫梁、主梁等結(jié)構(gòu)件；側(cè)板與橫梁、主梁等各結(jié)構(gòu)件之間的連接鉚釘，采用考慮剪切變形影響的Timoshenko梁單元模擬；篩箱與地面之間的連接金屬螺旋彈簧用彈簧單元模擬。離散化后，共生成251997個節(jié)點及158532個單元。有限元模型如圖1所示。

1.2 模態(tài)分析

選取Block Lanczos法對振動篩進(jìn)行模態(tài)分析，計算得振動篩的第一階固有頻率為36.108Hz，其第一階振型如圖2所示。

圖1 振動篩有限元模型

圖2 第一階振型圖

由圖2可見第一階振型主要是側(cè)板發(fā)生的變形，是局部振動。振動篩工作時，側(cè)板是結(jié)構(gòu)主要承載件，分布著很大的慣性力，是易損壞部件。因此，提高振動篩的固有頻率，使其固有頻率遠(yuǎn)離激勵頻率，可減小側(cè)板的動應(yīng)力，延長其使用壽命。

2 靈敏度分析

在結(jié)構(gòu)設(shè)計中，常常需要對結(jié)構(gòu)某些部位做出適當(dāng)修改，以提高結(jié)構(gòu)的動態(tài)特性。對于振動篩這樣復(fù)雜結(jié)構(gòu)，影響其剛度和強度的因素有很多，并且整機優(yōu)化的計算量非常大。如果僅憑設(shè)計人員的設(shè)計經(jīng)驗進(jìn)行優(yōu)化設(shè)計，由于缺乏理論依據(jù)，自然避免不了盲目性，也不能做到以最小的改動達(dá)到最佳的效果。

為解決這一問題，對振動篩進(jìn)行靈敏度分析，確定敏感變量，即找出對系統(tǒng)第一階固有頻率影響最大的因素作為設(shè)計變量，提高整機優(yōu)化的效率。

2.1 基本原理

靈敏度研究動態(tài)特性對結(jié)構(gòu)參數(shù)的敏感程度，即系統(tǒng)目標(biāo)函數(shù)對設(shè)計變量的偏微分。

系統(tǒng)無阻尼自由振動微分方程為

(1)

[K]{u}=ω2[M]{u}

(2)

式中，ω為固有頻率,{μ}為固有振型向量。

根據(jù)式(2)對設(shè)計變量xi求偏導(dǎo)，得

(3)

將式(3)左乘{(lán)μ}T，由于[K]為對稱矩陣，整理得

(4)

由方程(2)[K]{u}=ω2[M]{u}，可將式(4)簡化為

(5)

將振型向量對質(zhì)量矩陣做歸一正則化處理，并簡化式(5)，且ω=2πf，得到系統(tǒng)的固有頻率對設(shè)計變量的靈敏度關(guān)系式為

(6)

公式(6)反映了系統(tǒng)固有頻率f隨設(shè)計參數(shù)變化的關(guān)系。靈敏度值為正數(shù)時，說明增加相應(yīng)變量的值會提高第一階固有頻率，反之則應(yīng)減小變量值來提高頻率。靈敏度的絕對值越大，對結(jié)構(gòu)的影響就越顯著。

2.2 設(shè)計變量的選擇

對振動篩這類復(fù)雜結(jié)構(gòu)，修改方案有很多種，可供調(diào)整的結(jié)構(gòu)參數(shù)也有很多，為確定何種方案最為有效，需分析結(jié)構(gòu)參數(shù)或設(shè)計變量對結(jié)構(gòu)特性的敏感程度。

振型幅值較大的區(qū)域的變動對該階固有頻率的影響較明顯；若某部分結(jié)構(gòu)處于振型節(jié)點附近，則這部分變動不會對該階頻率有明顯的影響[10]。振動篩第一階的變形主要是側(cè)板的局部變形，為降低優(yōu)化問題的難度，減少設(shè)計變量個數(shù)的同時確定有效的設(shè)計變量，盡量多選取側(cè)板上變形區(qū)域附近的變量進(jìn)行靈敏度分析。根據(jù)上述原則，選取如圖3所示與側(cè)板相關(guān)的12個變量進(jìn)行模態(tài)靈敏度分析。

圖3 變量名及含義

注：X1、Y1—上橫梁1的X和Y方向的位置改變量；X2、Y2—上橫梁2的X和Y方向的位置改變量；X3、Y3—上橫梁3的X和Y方向的位置改變量；X5、Y5—下橫梁5的X和Y方向的位置改變量；X6、Y6—下橫梁6的X和Y方向的位置改變量；JGWZ1—斜角鋼12在Y方向的位置改變量；JGJD—斜角鋼12繞頂點A沿逆時針方向的角度變化。

根據(jù)公式(6)，編寫出相應(yīng)的振動篩靈敏度計算程序。經(jīng)過計算，得到各變量對第一階固有頻率的靈敏度，如圖4所示。

圖4 一階靈敏度的柱形圖

由圖4可以看出，X1、X2、Y1、Y2、JGWZ1、JGJD這6個變量的靈敏度數(shù)值較大，說明在變形區(qū)域這些零件的位置變化對第一階固有頻率敏感；而變量X3、Y3、X5、Y5、X6、Y6靈敏度值較小，說明其改變對結(jié)構(gòu)第一階固有頻率影響非常小。變量JGJD的值為負(fù)，說明當(dāng)JGJD減小時，頻率會有所增加。是在確保部件之間不相互干涉的條件下，變量JGJD可取的變化范圍非常小，所以當(dāng)JGJD改變時，很容易與兩側(cè)的U型加強板及彈簧座耳干涉；并且變量JGJD對第一階固有頻率的貢獻(xiàn)并不非常明顯，故在優(yōu)化時沒有把變量JGJD列為設(shè)計變量。變量X2的靈敏度值相對于其他五個變量的靈敏度值小，但X2、Y2是影響上橫梁位置的一對變量，由于結(jié)構(gòu)裝配工藝要求，Y2位置改變時，X2也必需隨之改變。

由靈敏度的計算，綜合考慮各個變量對結(jié)構(gòu)的影響程度，選擇X1、Y1、Y2、JGWZ1這4個變量進(jìn)行優(yōu)化設(shè)計。計算得到的靈敏度系數(shù)均為正值，說明這些變量值增加時振動篩的第一階固有頻率也隨之提高，這樣確定了設(shè)計變量的變化趨勢，縮小了設(shè)計變量的取值范圍。

3 振動篩的優(yōu)化設(shè)計

根據(jù)結(jié)構(gòu)裝配工藝要求，位于側(cè)板內(nèi)外側(cè)的斜角鋼12和上橫梁2是用同一組鉚釘相連接，二者位置始終相互對應(yīng)。結(jié)構(gòu)上要求鉚釘不能隨其中某一部件位置的改變而移動，所以他們必須是同時變動。故優(yōu)化設(shè)計時將變量JGWZ1和Y2作為一個設(shè)計變量，即令Y2=JGWZ1。同樣設(shè)計變量X2與Y2也存在相應(yīng)的幾何約束關(guān)系。因為斜角鋼與水平(x軸)的夾角為33.5度，故X2與Y2之間的比例關(guān)系為tan33.5=Y2/X2,所以X2=Y2/tan33.5≈1.5×Y2。因此用設(shè)計變量Y2代替JGWZ1和X2，則設(shè)計變量由5個變成3個：X1、Y1和Y2。這樣處理，不僅避免了由于設(shè)計變量過多而收斂到局部最優(yōu)點的危險，同時也可提高計算效率。

在保證結(jié)構(gòu)件之間不相互干涉的條件下，確定設(shè)計變量的變化范圍。各變量的上限值均為各部件之間不發(fā)生干涉情況下的最大變動量。振動篩優(yōu)化的數(shù)學(xué)模型為

Minf(X)=min(-freq_1)

61-Y2≥0

47-X1≥0

s.t. 42-Y1≥0

Y2≥0

X1≥0

Y1≥0

設(shè)計變量和目標(biāo)函數(shù)freq_1的迭代曲線分別如圖5和圖6所示。根據(jù)計算結(jié)果并通過數(shù)據(jù)圓整，得到各設(shè)計變量值分別為X1=42mm，Y1=36mm，Y2=59mm，則JGWZ1=Y2=59mm，X2=1.5×JGWZ1=1.5×59=88.5mm。各部件的位置改變在圖7中標(biāo)出。

圖5 X1、Y1、Y2迭代曲線

圖6 目標(biāo)函數(shù)f(X)迭代曲線

圖7 優(yōu)化后的振動篩模型

圖8 優(yōu)化后的振動篩第一階振型圖

按照優(yōu)化結(jié)果修改結(jié)構(gòu)，再對振動篩進(jìn)行模態(tài)分析，其優(yōu)化后的固有頻率列于表1，對應(yīng)的第一階振型如圖8所示。

表1 修改前后固有頻率

通過對結(jié)構(gòu)進(jìn)行靈敏度分析可以找出最佳的修正位置，在此基礎(chǔ)上進(jìn)行結(jié)構(gòu)優(yōu)化設(shè)計，可提高優(yōu)化效率。結(jié)果表明，在靈敏度較大部位進(jìn)行結(jié)構(gòu)修改的方法是行之有效的。

4 動力學(xué)分析

對振動篩優(yōu)化設(shè)計后，需進(jìn)行動力學(xué)分析，確定所受的應(yīng)力以及位移是否滿足工作要求。

振動篩所受的簡諧激振力在水平方向與鉛垂方向的分力分別為349865 sin(ωt)N和-349865 sin(ωt)N，對優(yōu)化后的振動篩進(jìn)行動力學(xué)分析后，得到振動篩關(guān)鍵部件的應(yīng)力以及位移，如表2所示。

由表2可以看出，下橫梁所受的應(yīng)力最大，為40.093MPa。圖9為下橫梁應(yīng)力隨時間變化曲線，由圖9可以看出，最大應(yīng)力出現(xiàn)在t=0.00165s時，此時對應(yīng)的應(yīng)力云圖如圖10所示。

圖9 下橫梁隨時間t變化曲線

圖10 t=0.00165s時下橫梁應(yīng)力云圖

5 結(jié)論

建立香蕉型直線振動篩的有限元模型，對其進(jìn)行模態(tài)分析，得到固有頻率以及相應(yīng)的振型。通過靈敏度分析篩選出對第一階固有頻率影響較大的因素作為設(shè)計變量，進(jìn)行優(yōu)化設(shè)計，避免了優(yōu)化設(shè)計的盲目性，提高計算效率。優(yōu)化圓整后，第一階固有頻率提高了5.58%。對優(yōu)化之后的振動篩進(jìn)行強度校核，驗證其滿足強度要求。

[1] 楊獻(xiàn)輝，趙俊利，薛雷澤.基于ANSYS的SDLl645-1圓周振動篩模態(tài)特性分析[J].煤礦機械，2010,31(6):74-76.

[2] 王永巖，邰英樓，孫奇涵.大型直線振動篩動態(tài)仿真的研究[J].計算力學(xué)學(xué)報，2001,18(1):94-98.

[3] 吳炳勝，蘇茹茹，李文選.基于ANSYS的ZK450振動篩動態(tài)特性的分析[J].煤礦機械，2012,33(2):97-99.

[4] 孟彩茹，李磊，馮忠緒.振動篩的強度分析和壽命估算[J].機械設(shè)計，2012,29(7):20-24.

[5] 賀孝梅,劉初升,張成勇.基于多頻約束和解析靈敏度法的大型振動篩優(yōu)化設(shè)計[J].中南大學(xué)學(xué)報,2011,42(3):664-670.

[6] 方治華，趙爽.HFZS1640型振動篩強度的優(yōu)化設(shè)計[J].煤礦機械，2007,28(11):23-24.

[7] 喬崇全，王克杰.大型振動篩動力學(xué)分析及動態(tài)優(yōu)化設(shè)計[J].煤礦機械，2011,32(11):3-5.

[8] 周磊，劉建慧.振動篩側(cè)板及篩框動態(tài)設(shè)計[J].機械設(shè)計與研究，2010,26(4):108-110.

[9] 趙海軍，鄧兆祥，胡玉梅.基于靈敏度分析的微型貨車車身結(jié)構(gòu)動力修改[J].汽車工程，2008,30(6):534-538.

[10]張軍挪，王瑞林，李永建.基于靈敏度分析的機槍結(jié)構(gòu)動力修改研究[J].中北大學(xué)學(xué)報，2007,28(2):112-116.

[11]成大先.機械設(shè)計手冊[M].(第一卷).北京：化學(xué)工業(yè)出版社,2008:107-109.