關(guān)玉明，于 盼，李 朝，崔 佳，商 鵬

（河北工業(yè)大學(xué) 機(jī)械工程學(xué)院，天津 300130）

基于Workbench的濃密機(jī)中心軸吊點位置的優(yōu)化

關(guān)玉明，于 盼，李 朝，崔 佳，商 鵬

（河北工業(yè)大學(xué) 機(jī)械工程學(xué)院，天津 300130）

針對中心傳動式濃密機(jī)的中心軸整體變形、等效應(yīng)變和最大應(yīng)力較大的情況，介紹了一種基于Workbench優(yōu)化技術(shù)進(jìn)行中心軸吊點位置選取的方法.對中心軸進(jìn)行受力分析，定義優(yōu)化變量，將整體變形、等效應(yīng)變和最大應(yīng)力作為目標(biāo)函數(shù)，基于Workbench的不斷迭代計算得出幾種優(yōu)化方案和目標(biāo)函數(shù)曲線圖，經(jīng)比較選出合適的方案，并對選取的方案做了進(jìn)一步的驗證.結(jié)果表明，優(yōu)化后的整體變形、等效應(yīng)變和最大應(yīng)力有了顯著的改善.

Workbench優(yōu)化；中心軸；吊點位置；整體變形；等效應(yīng)變；最大應(yīng)力

濃密機(jī)是基于重力沉降作用的固液分離設(shè)備，通常為由混凝土、木材或金屬焊接板作為結(jié)構(gòu)材料建成帶錐底的圓筒形淺槽，可將含固量為10%～20%的礦漿通過重力沉降濃縮為含固量為45%～55%的底流礦漿.本文討論的濃密機(jī)為中心傳動式，由液壓馬達(dá)驅(qū)動，動力通過減速器后傳至中心軸使其轉(zhuǎn)動，并帶動安裝在中心軸上的耙架也隨之轉(zhuǎn)動，這樣便可攪動礦漿.在運(yùn)動過程中，漿料被攪動時所產(chǎn)生的扭矩直接作用在耙架上，造成中心軸應(yīng)發(fā)生整體變形，這對中心軸的設(shè)計帶來了一定的困擾.因此，在設(shè)計中心軸時，研究合理的耙架吊點位置使得中心軸吊點處的整體變形、等效應(yīng)變以及最大應(yīng)力達(dá)到最優(yōu)是十分必要的.

隨著計算機(jī)技術(shù)的飛速發(fā)展，結(jié)構(gòu)優(yōu)化算法也取得了很大的發(fā)展，目前結(jié)構(gòu)優(yōu)化經(jīng)歷了較低層次的尺寸優(yōu)化，較高層次的結(jié)構(gòu)形式優(yōu)化和如今更高層次的拓?fù)鋬?yōu)化.李純金等[1]利用遺傳算法，對隨車起重機(jī)折疊臂架鉸點位置進(jìn)行了優(yōu)化，不但降低了液壓缸在工作過程中所受到的最大值，而且使整個液壓系統(tǒng)壓力更趨于穩(wěn)定；王琪等[2]利用粒子群算法，對雙層剪式液壓升降臺液壓缸位置進(jìn)行仿真優(yōu)化分析，使升降臺升降更加平穩(wěn)，很好地改善了其工作性能，大大提高了其使用壽命.侯忠明等[3]在ADAMS中建立參數(shù)化伸縮臂叉車模型，利用設(shè)計研究及試驗設(shè)計功能對連桿變幅機(jī)構(gòu)進(jìn)行參數(shù)化分析，找到最優(yōu)鉸點位置，使得液壓缸最大受力減小，實現(xiàn)了優(yōu)化目標(biāo).曾剛[4]通過Pro/E軟件中的優(yōu)化設(shè)計功能對煤礦廠給料機(jī)的定位尺寸進(jìn)行優(yōu)化，使得更多的物料能以自溜的方式進(jìn)入給料機(jī)，從而減少了地面鏟車的數(shù)量及工作量，節(jié)約了成本.本文中主要針對中心軸與耙架的吊點位置建立數(shù)學(xué)模型，借助Workbench尋找最優(yōu)解，得出吊點的最優(yōu)位置，以達(dá)到優(yōu)化的目的.

1 靜力學(xué)分析

本文討論的新型濃密機(jī)為中心傳動式，中心軸與耙架的安裝關(guān)系如圖1所示，在濃密池的正中央安有一根中心軸1，中心軸的末端固定有一個十字形耙架，耙架有長短耙架之分，2組耙架對稱地安裝在中心軸上，耙架的下端裝有刮板，耙架與地面成9°角.

1.1 力學(xué)模型的建立

由于耙架對稱安裝，長短耙架的吊點位置位于同一水平高度，為了提高軟件的運(yùn)算效率，本文以一端耙架的吊點位置為主進(jìn)行分析計算.圖2為中心軸與耙架的機(jī)構(gòu)簡圖，以O(shè)點為坐標(biāo)原點，吊點處為A，在Y軸上，A點坐標(biāo)記為（0，d），短耙架吊點處為B，記為（xB，yB）=（5.685，1.641），OB與Y軸夾角為81°，AB與X軸的夾角記為θ.濃密機(jī)工作時，耙及耙架需要克服漿液對其產(chǎn)生的扭矩T1.

由圖2可以看出中心軸吊點處的拉力與短耙架吊點處的拉力互為一對作用力與反作用力.因此，可利用短耙架的求解來計算中心軸吊點處的受力.對中心軸和短耙架進(jìn)行靜力學(xué)分析得：

式中：T為驅(qū)動裝置對中心軸的扭矩，N·m；T1為漿液對短耙架整體的最大扭矩，N·m；F為中心軸吊點處的處所受到的拉力，N；F′為短耙架吊點處所受到的拉力，N.

1.2 短耙架靜力學(xué)分析

漿液對短耙架整體的最大扭矩T1可以認(rèn)為是由耙架所受的扭矩TZ與耙所受的扭矩Tl疊加作用產(chǎn)生[5]，即

1） 耙架扭矩分析

耙架和耙轉(zhuǎn)動范圍的漿液濃度較低，該區(qū)域漿液按流體分析.耙架與中心軸相連，其運(yùn)動軌跡是圓柱體，如圖3所示.耙架在運(yùn)動過程中要克服漿液對其的屈服應(yīng)力和周圍一定區(qū)域內(nèi)漿液對其的剪切應(yīng)力，因此，耙架所受到的扭矩TZ是剪切圓柱體側(cè)面的扭矩Ts和上下端面的扭矩Te的總和[5].將扭矩表示為剪力的函數(shù)為：

圖1 中心傳動式濃密機(jī)傳動件模型圖Fig.1 The model of the transmission part of the central drive thickener

圖2 中心軸與耙架的機(jī)構(gòu)簡圖Fig.2 Schematic diagram of the central shaft and harrow frame

圖3 耙架運(yùn)動軌跡及應(yīng)力分布Fig.3 Track and stress distribution of the harrow frame

式中：τs為耙架側(cè)面處的剪切應(yīng)力，Pa；τe為耙架上下端面處的屈服應(yīng)力，Pa；l為耙架的長度，m，l=5.917 m；r為剪切圓柱體的半徑，m，r=lcos81°；δ為剪切圓柱體的高度，m，δ=H=0.303 m.

由于耙架接觸部分流動的漿液濃度基本相同，而且耙架高度不大，因此認(rèn)為耙架上下端面以及側(cè)面所受剪切應(yīng)力大小相等，如圖3所示，即耙架所受剪切應(yīng)力為τ0=τs=τe=50 Pa[6-7].

2） 耙的扭矩分析

耙的運(yùn)動軌跡是圓環(huán)體，如圖4所示，因此會產(chǎn)生內(nèi)外圓面和上下端面4個剪切面.由于耙與耙架相連，在對耙架進(jìn)行分析時，已經(jīng)對圓環(huán)體的上端面進(jìn)行了分析，所以在這里省去這一部分的計算.耙運(yùn)動時的最大扭矩表示為屈服應(yīng)力的函數(shù)為

式中：τy為耙的屈服應(yīng)力，Pa，τy=50 Pa；h為耙的高度，m，h=0.6m；k為耙的個數(shù)，k=7；ri1和ri2分別為第i個剪切內(nèi)圓柱面與外圓柱面的半徑，m.

根據(jù)以上分析，可以得出漿液對短耙架整體的最大扭矩T1的應(yīng)力函數(shù)表達(dá)式

圖4 耙運(yùn)動軌跡及應(yīng)力分布Fig.4 Trajectory and stress distribution of harrow

短耙架受力分析如圖5所示.對其進(jìn)行靜力學(xué)分析，得出如下方程組：

式中：F′為鉸點處的作用力，N；G為短耙架與耙的重力，N，G=2 038 N；x0為耙架與耙的質(zhì)量中心，m，x0=3.098m；θ為拉桿與x軸的夾角，cosθ=0.999 1，sinθ=0.02.

聯(lián)立式（1） ～式（6），計算得出中心軸受到來自驅(qū)動的扭矩T=82 100 N·m，中心軸鉸點出的拉力F=15 472 N.

圖5 短耙架受力分析圖Fig.5 Force analysis of short frame

2 ANSYS Workbench優(yōu)化設(shè)計

ANSYS Workbench是利用APDL（參數(shù)化設(shè)計語言）對結(jié)構(gòu)設(shè)計的參數(shù)進(jìn)行調(diào)整，它要求優(yōu)化的模型必須為參數(shù)化模型.優(yōu)化的是通過算法不斷的循環(huán)“分析——評估——修正”[8]來實現(xiàn)的.ANSYS優(yōu)化模塊中常用的算法有零階法、一階法、隨機(jī)搜索法、最有梯度法[9]，其中，隨機(jī)搜索法完成指定次數(shù)的迭代分析，可以合理確定優(yōu)化空間[10-11]，是零階法和一階法的其前期處理步驟；零階法通過隨機(jī)搜索建立狀態(tài)變量，使其與目標(biāo)函數(shù)逐步逼近，廣泛適用于各類工程問題，但是其收斂速度較慢[12]；一階法的結(jié)果準(zhǔn)確，但是計算量大，容易陷入局部最優(yōu)解.

綜合考慮各種優(yōu)化算法的特點，本次優(yōu)化設(shè)計決定使用零階法優(yōu)化算法對模型進(jìn)行求解運(yùn)算.

圖6是中心軸的三維模型圖，為了提高運(yùn)行效率，在建模過程中對模型進(jìn)行了簡化，忽略了一些局部特征.中心軸的外形尺寸為：長度7 110 mm，直徑457 mm.

使用Workbench對原有中心軸進(jìn)行分析，得出如圖7所示整體變形云圖.中心軸在啟動瞬時所受到的扭矩最大，因此研究在該時刻中心軸的整體變形是有必要的.由圖7可以看出，中心軸最大變形量達(dá)到0.001 239 2 m，位于與減速器連接的部位，在結(jié)構(gòu)改進(jìn)時應(yīng)保證具有較小的整體變形量.

圖6 中心軸模型圖Fig.6 Central axis model

2.1 參數(shù)化建模

參數(shù)化建模基本步驟如下.

1） 定義材料屬性.中心軸的材料為Q235，材料屬性見表1.

2）施加約束條件.中心軸受到來自驅(qū)動的扭矩T=82 100 N·m和吊點處的拉力F=15 472 N，在Workbench施加力時將拉力分解為Y軸和X軸的力Fx=-15 458 N，F(xiàn)y=-309 N（說明：負(fù)號表示力的方向與坐標(biāo)軸正方向相反）.

3） 網(wǎng)格劃分.由于中心軸尺寸較大，結(jié)構(gòu)復(fù)雜，為了提高運(yùn)算效率，設(shè)置Sizing（單元尺寸）為Medium.劃分網(wǎng)格設(shè)置完成后，中心軸被劃分為73 509個節(jié)點，33 060個有限單元.

4） 參數(shù)化設(shè)置.參數(shù)化設(shè)置見表2.

表2 鉸點處參數(shù)化設(shè)置Tab.2 Parameterized setting of hinge points

2.2 結(jié)果分析

利用Workbench進(jìn)行結(jié)構(gòu)優(yōu)化，根據(jù)各變量取值范圍，優(yōu)化軟件進(jìn)行迭代計算得到12個設(shè)計點方案，計算得出整體變形對設(shè)計點的曲線圖8a）、最大應(yīng)力對設(shè)計點的關(guān)系曲線圖8b）和等效應(yīng)變對設(shè)計點的曲線圖8c） .

由圖8a）可以看出，在設(shè)計點2處，中心軸的整體變形達(dá)到最小，因此將設(shè)計點2作為一個候選點.由圖8b）和圖8c）可以看出，在設(shè)計點7處等效應(yīng)變和最大應(yīng)力均達(dá)到最小值，因此將設(shè)計點7作為一個候選點；另外，等效應(yīng)變和最大應(yīng)力對設(shè)計點的曲線最后都收斂于設(shè)計點12，表明Workbench將設(shè)計點12作為最優(yōu)方案，因此將設(shè)計點12也作為一個候選點.綜合分析圖8可以得出3個設(shè)計點，各自優(yōu)化后的數(shù)據(jù)如表3所示.

圖7 原有中心軸整體變形云圖Fig.7 Overall deformation of the original central axis

表1 Q235材料屬性表Tab.1 Q235 material properties table

圖8 目標(biāo)變量對設(shè)計點的曲線圖Fig.8 The curve of the target variable to the design point

根據(jù)表3數(shù)據(jù)可以看出，設(shè)計點12處的綜合效果更優(yōu)，因此選擇設(shè)計點12作為本次優(yōu)化的方案.

由圖9可以看出，優(yōu)化后的吊點拉力明顯降低，整體變形量的最大值減小了24%，等效應(yīng)變減小了31%，最大應(yīng)力減小了31%，這顯著提升了中心軸的工作穩(wěn)定性.

3 總結(jié)

使用ANSYS Workbench對濃密機(jī)中心軸吊點位置在滿足使用要求的情況下進(jìn)行了整體變形、等效應(yīng)變和最大應(yīng)力最小值的優(yōu)化設(shè)計，在得出的方案中根據(jù)實際要求選取了最佳方案，并進(jìn)一步證明了方案的可行性.同時，也證明了使用ANSYS Workbench對結(jié)構(gòu)進(jìn)行優(yōu)化的方法是可行的.這種優(yōu)化方法操作簡單，計算快捷，為結(jié)構(gòu)優(yōu)化提供了新的方法.

表3 優(yōu)化設(shè)計點Tab.3 Optimal design points

圖9 優(yōu)化前后數(shù)據(jù)對比Fig.9 Comparison of data before and after optimization

[1] 李純金，褚月橋，周海波.基于遺傳算法隨車起重機(jī)折疊臂架鉸點位置優(yōu)化[J].江蘇科技大學(xué)學(xué)報（自然科學(xué)版），2014，28（1）：30-36.

[2] 王琪，張任，劉建，等.基于粒子群算法的雙層剪式液壓升降臺液壓缸位置的多目標(biāo)優(yōu)化問題研究[J].機(jī)械科學(xué)與技術(shù)，2015，34（8）：1229-1234.

[3] 侯忠明，王勝軍，辛濤.基于ADAMS的伸縮臂叉車鉸點位置優(yōu)化設(shè)計[J].起重運(yùn)輸機(jī)械，2009，23（6）：42－44.

[4] 曾剛.基于Pro/E的堆煤場給料機(jī)位置優(yōu)化設(shè)計[J].煤礦技術(shù)，2014，33（10）：200-201.

[5] 吳愛祥，焦華喆，王洪江，等.深錐濃密機(jī)攪拌刮泥耙扭矩力學(xué)模型[J].中南大學(xué)學(xué)報（自然科學(xué)版），2012，43（4）：1469-1474.

[6] 李輝，王洪江，吳愛祥，等.基于尾砂沉降與流變特性的深錐濃密機(jī)壓耙分析[J].北京科技大學(xué)學(xué)報，2013，35（12）：1553-1558.

[7] 王勇，王洪江，吳愛祥.基于高徑比的深錐濃密機(jī)底流濃度數(shù)學(xué)模型[J].武漢理工大學(xué)學(xué)報，2011，33（8）：113-117.

[8] 徐釗，錢作勤，張平，等.基于ANSYS優(yōu)化技術(shù)的換熱器法蘭的結(jié)構(gòu)設(shè)計[J].石油和化工設(shè)備，2009（6）：15-18.

[9] 夏衛(wèi)明，駱桂林，嵇寬斌.ANSYS優(yōu)化算法的研究及其在液壓機(jī)優(yōu)化設(shè)計中的應(yīng)用[J].鍛壓裝備與制造技術(shù)，2010（1）：43-47.

[10]陸山，魯馮杰.基于ANSYS的整體葉盤結(jié)構(gòu)優(yōu)化設(shè)計[J].航空動力學(xué)報，2012，27（6）：1218-1224.

[11]李倫未，陸山.基于ANSYS的多輻板風(fēng)扇盤結(jié)構(gòu)優(yōu)化設(shè)計技術(shù)[J].航空動力學(xué)報，2011，26（10）：2245-2250.

[12]寇劍鋒，鄭錫濤，彭興國.基于零階優(yōu)化方法的無人機(jī)復(fù)合材料進(jìn)氣道結(jié)構(gòu)優(yōu)化與分析[J].機(jī)械強(qiáng)度，2008，30（6）：1018-1022.

Optimization of hanging point position of thickener center shaft based on Workbench

GUAN Yuming,YU Pan,LI Zhao,CUI Jia,SHANG Peng

（School of Mechanical Engineering,Hebei University of Technology,Tianjin 300130,China）

A method for the selection of the hanging point position of the center shaft based on the Workbench optimization technique is introduced on the matter of the center shaft of the central drive thickener,such as the total deformation,the equivalent strain and the equivalent stress maximum.Force analyses are conducted on the center shaft and the optimization variables are defined.Put the total deformation,the equivalent strain and the equivalent stress maximum as the output variables.Several optimization schemes and the diagrams of the hinge point position with the optimization goal are obtained through the continuous iterative calculation of Workbench.The optimal scheme is chosen in accordance with the chart.Finally,further validation of the selected scheme is made.The results show that there is a big improvement in the optimization of the total deformation,the equivalent strain and the equivalent stress maximum.

Workbench optimization;center shaft;hinge position;the total deformation;the equivalent strain;the equivalent stress maximum

TH133.2

A

1007-2373（2017） 05-0028-05

10.14081/j.cnki.hgdxb.2017.05.005

2017-03-28

天津市科技特派員項目（14JCTPJC00532）

關(guān)玉明（1957-），男，教授，gyuming@163.com.

[責(zé)任編輯 楊 屹]