，，，

(福建工程學(xué)院 福建省汽車電子與電驅(qū)動(dòng)重點(diǎn)實(shí)驗(yàn)室，福建 福州 350118)

林業(yè)機(jī)械不斷向輕量化、高效率方向發(fā)展，摩擦卷筒作為林業(yè)絞盤機(jī)的重要?jiǎng)恿Σ考苯佑绊懼謽I(yè)絞盤機(jī)的工作效率及性能，因此，對(duì)摩擦卷筒進(jìn)行進(jìn)一步優(yōu)化很有必要[1]。國(guó)內(nèi)外學(xué)者對(duì)此已開展了一系列研究: Usabiaga等[2]在分析卷筒直徑對(duì)鋼絲繩影響的基礎(chǔ)上，對(duì)卷筒直徑進(jìn)行優(yōu)化設(shè)計(jì)；羅健康等[3]基于有限元方法，針對(duì)卷筒壁厚進(jìn)行了優(yōu)化設(shè)計(jì)；聶金寧等[4]采用微積分與有限元仿真分析相結(jié)合的方法，研究卷筒在不同載荷條件下的應(yīng)力狀態(tài)，并對(duì)卷筒進(jìn)行了優(yōu)化設(shè)計(jì)；吳傳宇等[5]借助有限元分析軟件對(duì)摩擦卷筒進(jìn)行了分析與改進(jìn)；張嘉璐等[6]基于有限元方法，對(duì)拋物線形和雙錐形摩擦卷筒結(jié)構(gòu)形式進(jìn)行了應(yīng)力及形變分析。以上研究?jī)H對(duì)摩擦卷筒進(jìn)行了結(jié)構(gòu)分析或尺寸參數(shù)優(yōu)化設(shè)計(jì)。為了得到力學(xué)性能更加良好的結(jié)構(gòu)，拓?fù)鋬?yōu)化成為了最常用的設(shè)計(jì)方法[7]。張小珍等[8]基于Optistruct軟件，以摩擦卷筒體積最小為優(yōu)化目標(biāo)，對(duì)摩擦卷筒進(jìn)行了拓?fù)浣Y(jié)構(gòu)優(yōu)化，但得到的優(yōu)化結(jié)果邊界不清晰。本文嘗試?yán)靡环N改進(jìn)的拓?fù)鋬?yōu)化方法對(duì)摩擦卷筒進(jìn)行優(yōu)化設(shè)計(jì)。

常見的拓?fù)鋬?yōu)化方法有水平集方法、均勻化方法、SIMP方法[9-10]、漸近結(jié)構(gòu)優(yōu)化方法(BESO)、獨(dú)立連續(xù)映射法(ICM)、移動(dòng)可變形組件法(MMC)等。BESO方法具有邊界清晰、易于實(shí)現(xiàn)的特點(diǎn)，但實(shí)際應(yīng)用中仍存在靈敏度推導(dǎo)繁瑣復(fù)雜[11-12]、刪除單元導(dǎo)致的數(shù)值奇異[13-14]、基于準(zhǔn)則法優(yōu)化時(shí)性能難以評(píng)價(jià)等一些問題[15]。針對(duì)這些問題，本文嘗試改進(jìn)基于von Mises應(yīng)力的雙向漸進(jìn)結(jié)構(gòu)拓?fù)鋬?yōu)化方法，并將其應(yīng)用于林業(yè)絞盤機(jī)摩擦卷筒的輕量化設(shè)計(jì)。

1 基于米塞斯應(yīng)力的BESO算法改進(jìn)

1.1 材料插值模型

BESO方法的設(shè)計(jì)變量是離散的(0或1)，結(jié)構(gòu)有限元網(wǎng)格劃分后，每個(gè)單元對(duì)應(yīng)一個(gè)設(shè)計(jì)變量，當(dāng)某個(gè)單元對(duì)應(yīng)的設(shè)計(jì)變量為0時(shí)，就意味著將該單元的材料刪掉，這在原理上行得通，但在實(shí)際的有限元數(shù)值計(jì)算時(shí)，直接將單元材料賦值為0可能會(huì)造成結(jié)構(gòu)剛度矩陣奇異[16-18]。為避免這種數(shù)值奇異現(xiàn)象，本文采用修改的材料插值模型：

(1)

其中，xi為i個(gè)單元的設(shè)計(jì)變量；E(xi)、ρ(xi)分別代表單元的彈性模量和密度；E0、ρ0分別代表xi=1時(shí)的彈性模量和密度；Emin、ρmin分別代表xi=0時(shí)的彈性模量和密度，通常取Emin=10-6E0、ρmin=10-6ρ0；p為懲罰系數(shù)(通常取3)。此時(shí)，單元“刪除”是通過賦予單元一種很“軟”的材料實(shí)現(xiàn)的，從而可避免結(jié)構(gòu)剛度數(shù)值奇異的問題。

1.2 應(yīng)力靈敏度分析

應(yīng)力水平是結(jié)構(gòu)設(shè)計(jì)的一個(gè)重要判斷標(biāo)準(zhǔn)，一個(gè)理想的結(jié)構(gòu)，內(nèi)部應(yīng)力應(yīng)處于同一應(yīng)力水平。以應(yīng)力水平為優(yōu)化準(zhǔn)則的BESO方法，原理簡(jiǎn)單、計(jì)算效率高，更適合于實(shí)際工程應(yīng)用。結(jié)構(gòu)中某一點(diǎn)處的Von Mises應(yīng)力為：

(2)

其中，σ為某一點(diǎn)處的應(yīng)力向量，T為系數(shù)矩陣。

根據(jù)有限元理論，第i個(gè)單元內(nèi)的應(yīng)力狀態(tài)可表示為：

σi=DBui

(3)

其中，D為單元的彈性矩陣；B為應(yīng)變矩陣；ui為單元節(jié)點(diǎn)位移。

將公式(1)代入公式(3)整理可得：

(4)

其中，D0表示單元為實(shí)體材料時(shí)的彈性矩陣。

再將公式(4)代入公式(2)，得到第i個(gè)單元內(nèi)等效的米塞斯應(yīng)力：

(5)

將公式(5)兩邊對(duì)xi求導(dǎo)得：

(6)

因此，基于von Mises應(yīng)力的單元靈敏度數(shù)可定義為：

(7)

可見，基于應(yīng)力準(zhǔn)則的拓?fù)鋬?yōu)化，其靈敏度計(jì)算很容易。在實(shí)際應(yīng)用中，單元的靈敏度數(shù)可通過單元各節(jié)點(diǎn)處的靈敏度數(shù)求平均得到：

(8)

式中，M為第i個(gè)單元的節(jié)點(diǎn)數(shù)，對(duì)于四邊形單元M=4，對(duì)于六面體單元M=8。

1.3 性能評(píng)價(jià)指標(biāo)

基于米塞斯應(yīng)力的BESO算法的基本原理就是逐步刪除結(jié)構(gòu)中應(yīng)力低的材料，從而使結(jié)構(gòu)的應(yīng)力分布逐漸均勻。在統(tǒng)計(jì)學(xué)中，均方差是用來衡量數(shù)據(jù)離散程度的，均方差越小說明數(shù)據(jù)越集中，反之，則越離散。結(jié)構(gòu)應(yīng)力是否均勻，同樣可以用應(yīng)力均方差來衡量[19]。因此，本文將應(yīng)力均方差作為基于應(yīng)力的BESO算法的性能評(píng)價(jià)指標(biāo)PI，即

(9)

1.4 改進(jìn)的BESO算法步驟及軟件實(shí)現(xiàn)

改進(jìn)后的BESO算法包括以下計(jì)算步驟：

1) 初始化設(shè)計(jì)區(qū)域，對(duì)設(shè)計(jì)區(qū)域進(jìn)行網(wǎng)格的劃分，給定目標(biāo)體積分?jǐn)?shù)及刪除率；

2) 對(duì)結(jié)構(gòu)進(jìn)行有限元方法分析；

3) 計(jì)算單元等效應(yīng)力，得到單元靈敏度數(shù)并進(jìn)行濾波處理[20]；

4) 根據(jù)單元靈敏度數(shù)大小，對(duì)單元進(jìn)行排序，按照刪除率對(duì)單元進(jìn)行刪、添；

5)重復(fù)以上過程，直至滿足收斂條件，得到最優(yōu)拓?fù)浣Y(jié)構(gòu)。

圖1給出了該算法的具體流程圖。

圖1 改進(jìn)后BESO算法流程圖Fig.1 The improved BESO method flowchart

改進(jìn)的BESO算法是通過聯(lián)合Abaqus有限元軟件與MATLAB編程語(yǔ)言來實(shí)現(xiàn)的，在算法實(shí)現(xiàn)過程中的關(guān)鍵步驟包括：

1)有限元軟件的調(diào)用。在MATLAB程序中，采用dos命令來調(diào)用Abaqus有限元軟件，命令格式為：status = dos ( 'abaqus job=file_name.inp' )，其中，status表示返回狀態(tài)；file_name表示有限元模型文件名。

2)有限元分析結(jié)果的讀取。通過讀取Abaqus的fil文件，將單元節(jié)點(diǎn)應(yīng)力數(shù)據(jù)保存為txt文件，然后再導(dǎo)入到MATLAB的workspace中。

3) 刪、添單元后新結(jié)構(gòu)的構(gòu)造?？赏ㄟ^修改inp文件來實(shí)現(xiàn)單元的刪、添，從而構(gòu)造新的結(jié)構(gòu)用于下一步有限元分析。

2 集材絞盤機(jī)摩擦卷筒設(shè)計(jì)

2.1 算法驗(yàn)證

為驗(yàn)證改進(jìn)的BESO算法，先來研究一個(gè)簡(jiǎn)單算例：旋轉(zhuǎn)飛輪的形狀優(yōu)化問題。在離心載荷作用下，盤型飛輪可看成二維軸對(duì)稱問題。假設(shè)飛輪由鋁合金制成，彈性模量為70 GPa、泊松比為0.3、密度為2 700 kg/m3；飛輪半徑為300 mm、軸向厚度為150 mm；飛輪以2 000 rad/s高速旋轉(zhuǎn)。設(shè)定BESO的進(jìn)化率ER=2%、懲罰系數(shù)p=3、濾波半徑rmin=3、目標(biāo)體積V*=50%。

圖2 (b)為優(yōu)化后飛輪的形狀及應(yīng)力分布。與初始狀態(tài)圖2(a)相比，最大應(yīng)力降低了43%，低應(yīng)力區(qū)域(深色)也明顯減少。圖3給出了優(yōu)化過程中飛輪結(jié)構(gòu)應(yīng)力均方差的變化曲線，從圖3可以看出，經(jīng)過優(yōu)化，飛輪內(nèi)部應(yīng)力變得更加均勻。此外，從圖2 (b)還可以看到，優(yōu)化后的飛輪形狀與等應(yīng)力飛輪的形狀十分相似。綜上可見，改進(jìn)的BESO方法是可行、有效的。下面采用該方法進(jìn)一步開展應(yīng)用研究。

圖2 旋轉(zhuǎn)飛輪優(yōu)化結(jié)果對(duì)比Fig.2 Omparison of flywheel optimization results

圖3 飛輪應(yīng)力均方與體積分?jǐn)?shù)變化差曲線Fig.3 Stress standard deviation convergence curve of flywheel

2.2 摩擦卷筒拓?fù)鋬?yōu)化設(shè)計(jì)

集材絞盤機(jī)是依靠摩擦卷筒與鋼索產(chǎn)生的摩擦力牽引鋼索運(yùn)行的，對(duì)摩擦卷筒進(jìn)一步優(yōu)化設(shè)計(jì)，有利于實(shí)現(xiàn)絞盤機(jī)輕量化和提高作業(yè)效率。摩擦卷筒主要承受的載荷是鋼索拉力對(duì)其作用的扭矩，設(shè)摩擦卷筒最大牽引力為19 208 N，根據(jù)文獻(xiàn)[7]，可以得到工作半徑為150 mm的摩擦卷筒所承受的最大扭矩約為2 881.2 N·m。摩擦卷筒的材質(zhì)一般為低合金高強(qiáng)鋼，彈性模量為210 GPa，泊松比為0.3，密度為7 900 kg/m3。摩擦卷筒是旋轉(zhuǎn)周期對(duì)稱結(jié)構(gòu)，圖4為1/4周期摩擦卷筒的初始結(jié)構(gòu)(旋轉(zhuǎn)角度為60°)，外半徑為150 mm，內(nèi)、外端徑向厚度都為20 mm，軸向最大厚度為70 mm，軸向最小厚度為50 mm；其中深色部分為非設(shè)計(jì)域，淺色部分為實(shí)際的設(shè)計(jì)域。其他優(yōu)化參數(shù)不變。

圖4 摩擦卷筒的初始結(jié)構(gòu)Fig.4 Initial model of friction drum

圖5給出了旋轉(zhuǎn)周期數(shù)為4、設(shè)計(jì)域材料體積分?jǐn)?shù)為40%時(shí)的摩擦卷筒最優(yōu)拓?fù)浣Y(jié)構(gòu)及其應(yīng)力分布。從圖中可看出，優(yōu)化后摩擦卷筒結(jié)構(gòu)邊界清晰，不存在中間過渡材料；且最大應(yīng)力為24.11MPa，遠(yuǎn)小于材料的許用應(yīng)力113 MPa，可見仍有進(jìn)一步優(yōu)化的空間。

圖5 旋轉(zhuǎn)周期數(shù)為4、材料體積分?jǐn)?shù)為40%時(shí)的最優(yōu)摩擦卷筒拓?fù)浣Y(jié)構(gòu)Fig.5 Optimal friction drum topology when the number of rotation periods is 4 and the material volume fraction is 40%

為進(jìn)一步研究材料體積分?jǐn)?shù)及旋轉(zhuǎn)周期數(shù)對(duì)摩擦卷筒結(jié)構(gòu)的影響，圖6、7、8分別給出了在體積分?jǐn)?shù)為20%、40%以及60%情況下，不同旋轉(zhuǎn)周期數(shù)對(duì)應(yīng)的摩擦卷筒的最優(yōu)拓?fù)浣Y(jié)構(gòu)。

圖6 材料體積分?jǐn)?shù)為20%時(shí)不同旋轉(zhuǎn)周期數(shù)的摩擦卷筒最優(yōu)拓?fù)浣Y(jié)構(gòu)Fig.6 Topological optimization results of the friction drum with a volume fraction of 20%

圖7 材料體積分?jǐn)?shù)為40%時(shí)不同旋轉(zhuǎn)周期數(shù)的摩擦卷筒最優(yōu)拓?fù)浣Y(jié)構(gòu)Fig.7 Topological optimization results of the friction drum with a volume fraction of 40%

圖8 材料體積分?jǐn)?shù)為60%時(shí)不同旋轉(zhuǎn)周期數(shù)的摩擦卷筒最優(yōu)拓?fù)浣Y(jié)構(gòu)Fig.8 Topological optimization results of the friction drum with a volume fraction of 60%

對(duì)比圖6-8可看出，不同約束條件，摩擦卷筒的最優(yōu)拓?fù)洳槐M相同。隨著旋轉(zhuǎn)周期數(shù)的增加，摩擦卷筒內(nèi)的孔洞也跟著增多，但結(jié)構(gòu)形式變化不大；隨著材料體積分?jǐn)?shù)減小，摩擦卷筒的孔洞逐漸增大，結(jié)構(gòu)形式變得明顯不同，結(jié)構(gòu)中 “桿件”增多，結(jié)構(gòu)趨于簡(jiǎn)單化。在實(shí)際設(shè)計(jì)生產(chǎn)中，應(yīng)綜合考慮使用的材料性能及加工成本，選擇經(jīng)濟(jì)實(shí)用的拓?fù)浣Y(jié)構(gòu)。

3 結(jié)論

1)改進(jìn)的BESO算法無(wú)需復(fù)雜的公式推導(dǎo)，具有簡(jiǎn)單高效、易于實(shí)現(xiàn)的特點(diǎn)，便于實(shí)際工程應(yīng)用。

2)采用改進(jìn)的BESO方法可得到結(jié)構(gòu)合理、邊界清晰的集材絞盤機(jī)摩擦卷筒最優(yōu)拓?fù)浣Y(jié)構(gòu)，為摩擦卷筒輕量化設(shè)計(jì)與開發(fā)提供了理論和技術(shù)支撐。

3) 不同約束條件，摩擦卷筒的最優(yōu)拓?fù)洳槐M相同；隨著旋轉(zhuǎn)周期數(shù)目增加，摩擦卷筒結(jié)構(gòu)形式變化不大，但隨著材料體積分?jǐn)?shù)減小，摩擦卷筒結(jié)構(gòu)形式變化明顯，摩擦卷筒中“桿件”增多，結(jié)構(gòu)趨于簡(jiǎn)單化。