張大斌，余朝靜，張紀利，羅建欽，韋建玉*，雷　焱

(1.貴州大學 機械工程學院，貴州 貴陽 550025；2.廣西中煙工業(yè)有限責任公司，廣西 南寧 530001)

?

山地烤煙苗期揭膜機機架結構優(yōu)化

張大斌1，余朝靜1，張紀利2，羅建欽2，韋建玉2*，雷焱1

(1.貴州大學 機械工程學院，貴州 貴陽 550025；2.廣西中煙工業(yè)有限責任公司，廣西 南寧 530001)

為降低山地烤煙苗期揭膜機機架質量，本文通過DesignModeler建立了山地烤煙苗期揭膜機機架的幾何模型，對揭膜機勻速作業(yè)和調頭兩種工況下的機架進行了應力分析，采用多目標遺傳算法對機架結構進行了優(yōu)化。結果表明：機架最大應力出現在揭膜機掉頭工況下橫梁方鋼管與立柱連接處，應力值為123.89 MPa；優(yōu)化后機架總質量為5.9564 kg，比優(yōu)化前減輕了19.84%；優(yōu)化后機架最大應力為129.84 MPa，雖比優(yōu)化前略重，但遠小于機架所用材料Q235的最大許用應力值235 MPa。

山地烤煙揭膜機；苗期揭膜機；揭膜機機架；有限元分析

國內已有不少研究機構對揭膜機進行研究，新疆建設兵團123團修造廠生產了2JMSD-4.5型棉花揭膜機[1]。新疆兵團農機推廣站和134團共同研發(fā)了CSM-130B型齒鏈式懸掛收膜機[2]。國外的研究要早于國內，在羅伯森和索依爾研制的地膜回收機中，其工作原理也是通過卷膜輥將地膜直接卷起完成收膜[3]。

然而，國內外現有的揭膜機重量均較大，需要大功率動力機構進行牽引，不適合山區(qū)復雜地形作業(yè)。針對山區(qū)地形研究山地煙草苗期揭膜機將幫助山區(qū)煙農快速高效地回收地膜，減輕煙農的作業(yè)強度，降低地膜的殘留率，減輕地膜對土壤的污染。揭膜機機架作為各部件的承載者，其結構的合理設計對于整機的性能好壞至關重要。面對復雜的山地作業(yè)工況，機架的強度要求較高，同時山區(qū)地形的不規(guī)則性又對整機質量和體積有一定的要求，為了便于機器的轉移，機器質量不宜過重。

機架結構作為揭膜機的主要構件，其結構設計是否合理關系到整機的質量和性能[4]。因此，對機架結構進行優(yōu)化設計對山地烤煙揭膜機輕量化至關重要。本文將對機架方鋼管壁厚和鋼板厚度進行優(yōu)化設計，在機架最大應力處于安全范圍的前提下，求出最輕機架結構。

1　山地揭膜機機架結構靜態(tài)有限元分析

1.1山地揭膜機機架結構有限元模型建立

對揭膜機機架結構三維建模前，在能保證充分反映實際機架結構的前提下，為節(jié)省計算時間，需對機架機構做一定的簡化。主要包括忽略一些非承載件，比如螺栓、螺母等，這些非承載件對結構強度、剛度并沒有多大影響[5]。

ANSYS Workbench結構優(yōu)化設計必須進行參數化建模，而由其他大部分通用三維軟件建模導入的模型不能進行參數化，所以為了利用ANSYS Workbench對機架進行結構優(yōu)化最佳選擇是在ANSYS DesignModeler中建立三維模型，并對優(yōu)化對象進行參數化。

1仿形輪固定板；2松土鏟臂；3底盤方鋼管；4立柱；5橫梁方鋼管；6調節(jié)臂；7上肩；8切膜刀臂圖1　山地煙草苗期揭膜機機架模型

機架結構由薄壁方鋼管構成，會產生彎曲、拉壓和扭轉等變形。因此，根據其結構特點，在對機架結構強度、剛度進行分析時，單元的類型選擇以六面體為主，四面體與六面體結合的混合單元。采用這種混合單元所得到的網格質量較好，對提高計算精度、節(jié)省計算時間均有很大的幫助。同時，混合單元還能很好的解決局部過渡所產生的病態(tài)單元。機架結構方鋼管之間的焊縫采用共用節(jié)點來模擬[5]，所得有限元模型含95430個網格單元和455906個節(jié)點。

1.2機架結構典型工況下應力分析

對于機架結構質量，只需在ANSYS Workbench中輸入材料密度，程序便會根據所輸入的參數進行計算。機架材料為Q235，最大許用應力為235 MPa。進行機架結構的有限元分析前，需添加載荷和定義邊界條件。因機架與各零部件多以面接觸為主，所以在各接觸面施加均布載荷，對機架與仿形輪連接處和機架與田園管理機聯接處施加固定約束。選擇勻速作業(yè)、調頭兩種典型工況進行模擬計算，得出機架在勻速作業(yè)工況下的應力分布和最大應力位置如圖2、圖3所示，機架在調頭工況下的應力分布和最大應力位置如圖4、圖5所示。

圖2　勻速作業(yè)下應力分布

圖3　勻速作業(yè)下最大應力位置

圖4　調頭時應力分布

圖5　調頭時最大應力位置

綜合兩種分析工況可以看出，機架的應力分布不均勻。勻速作業(yè)時，機架底盤方鋼管存在應力集中現象，最大應力分布于底盤方鋼管開孔處，最大值為26.368 MPa。機器調頭工況下機架立柱與底盤方鋼和橫梁聯接處存在應力集中現象，最大應力出現在立柱與橫梁聯接處，最大值為123.89 MPa。Q235鋼材的屈服應力為235 MPa，兩種工況下機架都保持在安全范圍內。這只是機架結構的初步設計，該方案還存在著一定的弊端，為了最大限度地減小機架質量，同時減少應力集中現象，需要對整個機架各部分方鋼管壁厚進行結構優(yōu)化。

2　山地揭膜機機架結構優(yōu)化設計

2.1機架結構優(yōu)化模型

優(yōu)化設計以數學中最優(yōu)化理論為基礎，以計算機為手段，從多種方案中選擇最佳方案的設計方法。即在所選方案滿足了全部的設計要求的情況下，同時所需的支出(如質量、體積、應力應變等)最小[6]。ANSYS Workbench平臺提供多種優(yōu)化方法，本文中采用Multi-Objective Genetic Algorithm--多目標遺傳算法中的基于精英策略的多目標非支配排序遺傳算法[6]。它是在多目標遺傳算法的基礎上進行了改進，通過定義擁擠距離估計某個點周圍的解密度取代適應值共享來進一步提高計算效率和算法的魯棒性[7]。

設計變量、狀態(tài)變量和目標函數總稱為優(yōu)化分析內容[8]。在ANSYS Workbench優(yōu)化中，由用戶定義的參數來指定設計變量、狀態(tài)變量和目標函數[9]，其數學模型如下：

Vmin≤V≤VmaxV=[V1，V2…，Vn]

Wmin≤W≤WmaxW=[W1，W2…Wn]

minf(V)

式中：v為設計變量；w為狀態(tài)變量；f(v)為目標函數。

設計變量作為自變量，通過改變它的數值來獲得優(yōu)化結果。優(yōu)化設計過程中，為改善機架結構的應力分布和質量分布情況，可通過改變機架結構的各部位尺寸來實現，文中選取方鋼管的壁厚為設計變量。狀態(tài)變量是設計變量的函數。本文以減少機架結構應力集中為目的，應盡可能的減少機架結構的最大應力，因此，以機架結構方鋼管的最大應力為狀態(tài)變量[9]。

目標函數作為設計變量的函數，應盡量減小其數值[10]。本文優(yōu)化設計的目標主要是為了實現機架結構的輕量化，同時能滿足強度、剛度要求。因此，文中把機架結構的總質量作為目標函數[8]。

2.2機架結構優(yōu)化設計

由靜態(tài)結構分析可知，機架在調頭工況下最大應力明顯高于勻速作業(yè)工況，所以本文將在調頭工況對機架進行結構優(yōu)化設計。通過在ANSYS Workbench中設置優(yōu)化參數來對機架結構進行優(yōu)化設計，便可得到設計變量的優(yōu)化結果，結果如表1。

以車架底盤方鋼為例，其厚度在優(yōu)化過程中的變化情況如圖6所示。

表1　設計變量參數選擇及優(yōu)化結果

圖6　底盤方鋼管壁厚變化曲線

從圖6可知機架底盤方鋼壁厚在1.8～2.5 mm之間變化。在優(yōu)化的過程中，ANSYS Workbench不斷的對優(yōu)化參數進行調整，并隨著迭代次數的增加，方鋼管壁厚值變化的范圍會遞減，最終達到一個穩(wěn)定值。通過求解計算得到最優(yōu)結果為1.8344 mm，迭代次數為220次。

文中機架結構的最大應力值作為狀態(tài)變量，在機架底盤方鋼與立柱交接處初始最大應力為123.89 MPa。在優(yōu)化過程中，所得到的應力最大值為130 MPa。在優(yōu)化后，應力的最大值值變?yōu)?29.84 MPa。目標函數為機架總質量，初始質量為7.4308 kg，優(yōu)化后的質量為5.9564 kg，目標函數的變化情況如圖7所示。

圖7　質量變化曲線

從優(yōu)化結果可以看出，機架結構的總質量減少了19.84%，其效果明顯。雖然在實際的生產中，會受到方鋼管材料規(guī)格的限制和為滿足機架結構生產工藝要求，很難達到最優(yōu)的理想狀態(tài)，跟理論結果之間會存在一定的誤差[6]。但這為機架結構輕量化、合理化設計提供了理論基礎和參考。

2.3優(yōu)化設計前后對比分析

將優(yōu)化結果中的最佳候選點輸入參數復制到當前結構中，更新機架有限元模型，分別對兩種工況下機架靜態(tài)結構分析模型進行求解計算，對比前后兩種工況下的最大應力值，其結果如表2。

表2　改進前后兩種工況最大應力

從表2的比較結果可以看出，勻速作業(yè)工況下機架結構改進后最大應力有所增加，但機架的許用應力遠小于材料的許用應力，兩種工況下最大應力值變化很小，且都小于材料的最大許用應力值235 MPa。由兩種工況下的應力云圖可以看出，改進后應力分布更加的均勻。

3　結語

本文基于ANSYS Workbench平臺，對機架結構進行有限元模型的建立，選擇兩種典型工況進行有限元分析，得到機架結構的應力分布情況，機架

所受最大應力為123.89 MPa，且有明顯的應力集中現象。采用多目標遺優(yōu)化傳算法對機架進行優(yōu)化設計，得到機架最大應力值為129.48 MPa，小于材料的最大許用應力值235 MPa，機架質量減少19.84%，優(yōu)化效果明顯。現實生產中受材料規(guī)格等因素的制約，實際情形與理論結果還存在一定的差距，但本文結論可為機架的設計提供重要參考。

[1] 湯愛民.2JMSD-4.5型揭膜機、回收滴灌帶機[J].新疆農機化，2008(5)：12、25.

[2] 王能勇. CSM-130B型齒鏈式懸掛收膜機[J]. 新疆農機化，2000(3)：38-39.

[3] R.L.Parish.An automated machine for remove of plastic mulch.T-rans activns of ASAE，1998，42(1):49-51.

[4] 雷焱，張大斌，曹陽，等. 山地烤煙苗期地膜回收機卷膜輥結構設計[J]. 現代機械，2015，6:41-43.

[5] 許佩霞，蔡炳芳. 基于ANSYS的全地形車車架結構優(yōu)化設計[J]. 機械設計，2008，12:56-58.

[6] 寧波，呂志軍，婁文斌. 基于ANSYS Workbench的堆垛機結構分析與優(yōu)化[J]. 機械設計與制造，2012，6:71-73.

[7] 呂志軍.基于ANSYS Workbench的堆垛機結構優(yōu)化設計[C].見：全國先進制造技術高層論壇暨第十屆制造業(yè)自動化與信息化技術研討會，2011.8，山西太原，機械工業(yè)自動化分會或制造技術專業(yè)委員會，2011，8.

[8] 王燈，喻道遠，舒彪. 基于ANSYS的退火爐臺結構優(yōu)化設計[J]. 機械設計與制造，2012，4:41-43.

[9] 劉艷珍. 應用相對剛度相似法優(yōu)化大型塔內桁架梁結構[D].天津大學，2013.

[10]秦斌. 拱橋斜拉扣掛施工中的扣索索力優(yōu)化研究[D].重慶交通大學，2011.

(責任編輯：江龍)

Optimization Design of Mountain Flue-cured Tobacco Seedling Plastic MulchRemoving Machine Chassis' Structure Based on ANSYS Workbench

ZHANG Dabin1，YU Chaojing1，ZHANG Jili2，LUO Jianqin ，WEI Jianyu2*，LEI Yan1

(1. College of Mechanical Engineering，Guizhou University，Guiyang 550025，China; 2. Guangxi ZhongYan Co.， LTD，Nanning 530001，China)

To reduce the mass of mountain flue-cured tobacco seedling plastic mulch removing machine chassis， the geometric model of mountain flue-cured tobacco seedling plastic mulch removing machine chassis was established by Design Modeler. Then the stress on the chassis in constant speed working condition and turn around condition was analyzed and the maximum stress of the chassis at the connection between steel beam and column in the turn around condition was obtained. The maximum stress was 123.89MPa which was within the safty margin. The multi-objective genetic algorithms optimization method was used by putting the chassis'mass as objective function to optimize the structure of chassis in the turn around condition to obtain the optimal thickness of the key square steel tube's wall. The optimal mass of the chassis was 5.9564kg which was decreased 19.84%. The maximum stress of chassis was 129.84MPa which was less than 235MPa--the maximum allowable stress of Q235 which chassis used.

workbench; multi-objective genetic algorithms; chassis; optimal design

A

1000-5269(2016)03-0031-04

10.15958/j.cnki.gdxbzrb.2016.03.09

2016-03-11

廣西中煙工業(yè)有限責任公司科技項目“山地烤煙地膜、煙桿機械化回收關鍵技術及多功能機的研制”(合同編號：1212014015)貴州省科技計劃《貴州山地殘膜撿拾機具研發(fā)及應用研究》(黔科合GZ字(2015)3007)

張大斌(1976-)，男，副教授，博士，研究方向：農業(yè)機械、先進裝備及制造等，Email：332397268@qq.com.

韋建玉，Email：jtx_wjy@163.com.

TH122