石瑩,沈子強(qiáng)

(1.大連交通大學(xué) 機(jī)械工程學(xué)院,遼寧 大連 116028;2.中國(guó)鐵路沈陽(yáng)局集團(tuán)有限公司沈陽(yáng)動(dòng)車段,遼寧 沈陽(yáng) 110000)

輕量化技術(shù)是集產(chǎn)品設(shè)計(jì)、生產(chǎn)制造及原材料等于一體的綜合性應(yīng)用,其最終目的就是減小重量、增強(qiáng)性能以及降低成本,并使這三方面達(dá)到一個(gè)平衡,輕量化研究具有很高的經(jīng)濟(jì)和社會(huì)效益。齒輪作為機(jī)械設(shè)備中廣泛應(yīng)用的一種傳動(dòng)零件,其結(jié)構(gòu)輕量化研究對(duì)實(shí)現(xiàn)工程機(jī)械輕量化起到了重要的推動(dòng)作用。在保證強(qiáng)度要求的前提下,對(duì)齒輪進(jìn)行輕量化設(shè)計(jì)有利于降低制造成本、提高系統(tǒng)的動(dòng)力性、減小由于摩擦產(chǎn)生的噪聲、降低機(jī)械設(shè)備的動(dòng)力消耗、創(chuàng)造更高的經(jīng)濟(jì)效益,符合當(dāng)今社會(huì)發(fā)展所要求的綠色生產(chǎn)目標(biāo)[1]。

很多研究學(xué)者對(duì)齒輪輕量化進(jìn)行了卓有成效的研究。蔣進(jìn)科等[2]結(jié)合人字齒輪特點(diǎn),通過(guò)將齒面與修形面疊加,結(jié)合TCA、LTCA技術(shù),對(duì)齒面修形展開(kāi)研究,研究結(jié)果表明,合理的修形優(yōu)化可以達(dá)到減重、減振的效果。周新建等[3]應(yīng)用Romax軟件以螺旋線修形和齒廓修形為基礎(chǔ),對(duì)汽車1擋齒輪修形參數(shù)進(jìn)行優(yōu)化,并對(duì)修行參數(shù)設(shè)計(jì)分配,基于遺傳算法優(yōu)化求解,得到最優(yōu)參數(shù)組合。葉小芬等[4]以高鐵齒輪為研究對(duì)象,以減小體積和提高可靠性為目標(biāo),將強(qiáng)度作為約束條件,建立多目標(biāo)數(shù)學(xué)模型,應(yīng)用Matlab工具箱優(yōu)化求解,得到了體積最小且滿足約束要求的新傳動(dòng)機(jī)構(gòu)?，F(xiàn)階段,大部分齒輪輕量化設(shè)計(jì)研究仍停留在單一的尺寸優(yōu)化設(shè)計(jì)階段[5],通常只是改變齒輪的模數(shù)、齒數(shù)、齒寬等基本結(jié)構(gòu)參數(shù),將尺寸和拓?fù)鋬?yōu)化相結(jié)合對(duì)齒輪進(jìn)行設(shè)計(jì)的研究卻很少。

本文以某減速器齒輪為研究對(duì)象,以齒輪重量最小為目標(biāo)函數(shù),基于遺傳算法與Ansys聯(lián)合優(yōu)化,獲得最優(yōu)尺寸參數(shù);應(yīng)用變密度法對(duì)尺寸優(yōu)化后的齒輪輪輻進(jìn)行拓?fù)鋬?yōu)化,獲得材料最優(yōu)分布。靜力學(xué)分析結(jié)果表明,輕量化設(shè)計(jì)后的減速器齒輪完全滿足結(jié)構(gòu)強(qiáng)度要求。減速器齒輪優(yōu)化流程見(jiàn)圖1。

圖1 減速器齒輪優(yōu)化流程圖

1 減速器齒輪優(yōu)化仿真分析

1.1 齒輪參數(shù)化設(shè)計(jì)

由于后續(xù)要應(yīng)用Matlab遺傳算法工具箱和Ansys聯(lián)合對(duì)齒輪進(jìn)行尺寸參數(shù)優(yōu)化,應(yīng)用參數(shù)化建模便于創(chuàng)建齒輪有限元模型, 并使Ansys和Matlab之間數(shù)據(jù)傳輸成為可能,因此應(yīng)用Ansys中自帶的參數(shù)化設(shè)計(jì)語(yǔ)言(APDL)對(duì)減速器齒輪進(jìn)行建模,為后續(xù)的結(jié)構(gòu)優(yōu)化打下基礎(chǔ)。

本文采用3D四面體單元和映射網(wǎng)格劃分方式劃分網(wǎng)格。對(duì)齒輪副接觸區(qū)域進(jìn)行網(wǎng)格細(xì)化,單元尺寸設(shè)置為0.2,接觸區(qū)域外,網(wǎng)格單元尺寸設(shè)置為3。劃分完成之后齒輪有限元模型一共有79 702個(gè)單元,齒輪副有限元模型見(jiàn)圖2。

圖2 齒輪副有限元模型

對(duì)從動(dòng)齒輪中心孔節(jié)點(diǎn)施加全約束,對(duì)主動(dòng)齒輪中心孔節(jié)點(diǎn)施加徑向和軸向自由度約束,并施加軸向轉(zhuǎn)矩,取值為214 875 N·mm,轉(zhuǎn)矩按式(1)計(jì)算:

(1)

式中:T為轉(zhuǎn)矩,N·mm;P為功率,取值為18 kW;n為主動(dòng)輪轉(zhuǎn)速,取值為800 r/min。

1.2 仿真結(jié)果分析

對(duì)減速器齒輪進(jìn)行強(qiáng)度分析,齒面接觸應(yīng)力云圖和齒根彎曲應(yīng)力云圖見(jiàn)圖3、圖4。

圖3 齒輪接觸應(yīng)力云圖

圖4 齒根彎曲應(yīng)力云圖

圖3中單齒嚙合區(qū)最大齒面接觸應(yīng)力為443.355 MPa。齒輪傳動(dòng)過(guò)程中最大齒根彎曲應(yīng)力出現(xiàn)在單齒嚙合區(qū),這與文獻(xiàn)[6]中描述一致,其應(yīng)力值為201.58 MPa。該減速器齒輪許用接觸應(yīng)力為550 MPa,許用彎曲應(yīng)力為320 MPa,可知得到的應(yīng)力仿真結(jié)果均滿足結(jié)構(gòu)強(qiáng)度要求。

2 減速器齒輪尺寸優(yōu)化

利用MATLAB遺傳算法工具箱對(duì)減速器齒輪進(jìn)行尺寸優(yōu)化,在滿足約束條件的基礎(chǔ)上獲得最小重量時(shí)的齒輪結(jié)構(gòu)參數(shù)。

2.1 設(shè)計(jì)變量選取

設(shè)計(jì)變量必須是相互獨(dú)立的變量,選取齒輪模數(shù)m、主動(dòng)齒輪齒數(shù)z1、從動(dòng)齒輪齒數(shù)z2、齒寬系數(shù)φd4個(gè)參數(shù)為設(shè)計(jì)變量。即:X=[m、z1、z2、φd]。

2.2 目標(biāo)函數(shù)選取

以齒輪重量最小為目標(biāo)函數(shù),即齒輪材料密度ρ乘以齒輪體積V最小。

minf(X)=min(ρ·V)

(2)

2.3 約束函數(shù)選取

(1)模數(shù)

減速器齒輪的模數(shù)一般取2～6,因此齒輪模數(shù)在2、2.5、3、4、5、6中取值。

(2)齒數(shù)

為了防止齒輪發(fā)生根切或者造成主動(dòng)齒輪的尺寸過(guò)大,通常主動(dòng)齒輪齒數(shù)控制在17～24,因此主動(dòng)齒輪約束取[17,24]范圍內(nèi)的整數(shù)。

為了使齒輪傳動(dòng)效率不至于太低,所以傳動(dòng)比不能太小,通常從動(dòng)齒輪齒數(shù)控制在24～40,因此從動(dòng)齒輪約束取[24,40]范圍內(nèi)的整數(shù)。

(3)齒寬系數(shù)

本文研究的齒輪屬于開(kāi)式齒輪,開(kāi)式齒輪齒寬系數(shù)取值一般在0.1～0.3,因此齒寬系數(shù)約束為取[0.1,0.3]范圍的。

(4)強(qiáng)度約束條件

Ansys得到的齒面接觸應(yīng)力仿真值應(yīng)小于接觸疲勞許用應(yīng)力值, 齒根彎曲應(yīng)力仿真值應(yīng)小于彎曲疲勞許用應(yīng)力值。

2.4 Matlab與Ansys間的數(shù)據(jù)傳遞

利用遺傳算法思想,對(duì)所建立的減速器齒輪進(jìn)行優(yōu)化求解,種群設(shè)置為50,交叉概率為0.9、變異概率為0.05。

Matlab與Ansys之間數(shù)據(jù)傳遞流程如下:

(1) 將寫(xiě)有設(shè)計(jì)變量的APDL命令流寫(xiě)入數(shù)據(jù)文件.txt中。

(2) Ansys讀入文本文件中的數(shù)據(jù),進(jìn)行參數(shù)化建模并做分析計(jì)算,得到應(yīng)力及模型重量。

(3) Ansys將計(jì)算結(jié)果提取并寫(xiě)入數(shù)據(jù)的結(jié)果文件。

(4) 應(yīng)用Matlab遺傳算法工具箱讀取的文件數(shù)據(jù),對(duì)步驟(3)中的結(jié)果進(jìn)行計(jì)算[7]。

(5) 判斷是否滿足終止條件,如果不滿足條件,這些參數(shù)將被寫(xiě)入初始數(shù)據(jù)文件,替換原來(lái)的值,重復(fù)步驟(2)并依次向下進(jìn)行。

(6) 滿足終止條件時(shí),那么搜索終止,輸出結(jié)果,即為最優(yōu)解。

2.5 尺寸優(yōu)化結(jié)果

遺傳優(yōu)化求解后可以得到齒輪副的結(jié)構(gòu)參數(shù),遺傳優(yōu)化前后結(jié)果對(duì)比見(jiàn)表1。

表1 遺傳優(yōu)化前后結(jié)果對(duì)比

從表1可以看出,優(yōu)化后齒輪模數(shù)和齒數(shù)均減小,齒寬系數(shù)增大,齒輪傳動(dòng)比減小1.8%,對(duì)傳動(dòng)速度及效率基本無(wú)影響,保證了其動(dòng)力性能。齒輪副總重量減小了43.5%,減重明顯。對(duì)尺寸優(yōu)化后的減速器齒輪進(jìn)行靜力學(xué)分析,得到應(yīng)力云圖見(jiàn)圖5、圖6。

圖5 優(yōu)化后齒輪齒面接觸應(yīng)力云圖

圖6 優(yōu)化后齒輪齒根彎曲應(yīng)力云圖

可以看出,遺傳優(yōu)化后齒輪最大齒面接觸應(yīng)力由優(yōu)化前的443 MPa減小為412 MPa,最大齒根彎曲應(yīng)力由201 MPa減小為170 MPa,均小于齒輪結(jié)構(gòu)許用應(yīng)力,表明優(yōu)化后齒輪滿足其強(qiáng)度要求。

3 減速器齒輪拓?fù)鋬?yōu)化分析

3.1 齒輪拓?fù)鋬?yōu)化設(shè)計(jì)

本文應(yīng)用Ansys Workbench中的Topology Optimization模塊,對(duì)尺寸優(yōu)化后的主動(dòng)齒輪輪輻進(jìn)行拓?fù)鋬?yōu)化,對(duì)優(yōu)化區(qū)域進(jìn)行網(wǎng)格細(xì)化,得到的齒輪模型見(jiàn)圖7。

圖7 主動(dòng)齒輪拓?fù)鋬?yōu)化網(wǎng)格模型

采用變密度拓?fù)鋬?yōu)化方法,最大允許迭代500次,設(shè)計(jì)變量是輪輻區(qū)域網(wǎng)格單元的相對(duì)密度,相對(duì)密度值為0～1,以結(jié)構(gòu)柔度最小為目標(biāo)函數(shù),去除材料質(zhì)量為約束條件,對(duì)齒輪輪輻部分進(jìn)行拓?fù)鋬?yōu)化。

3.2 拓?fù)鋬?yōu)化分析結(jié)果

將約束條件設(shè)置為去除材料質(zhì)量60%、70%、80%這3種情況,分別進(jìn)行拓?fù)鋬?yōu)化計(jì)算,目標(biāo)函數(shù)迭代及優(yōu)化后模型見(jiàn)圖8。

(a) 去除材料質(zhì)量60%時(shí)

由圖8可以看出:

(1) 去除材料質(zhì)量為60%、70%、80%這3種情況下,目標(biāo)函數(shù)經(jīng)過(guò)多次迭代后均實(shí)現(xiàn)收斂。

(2) 隨著迭代次數(shù)的增加,結(jié)構(gòu)柔性呈下降趨勢(shì),最終趨于穩(wěn)定,結(jié)構(gòu)的柔度和剛度為此消彼長(zhǎng)的關(guān)系,柔度的下降說(shuō)明結(jié)構(gòu)剛度增加。

(3) 輪輻作為拓?fù)鋬?yōu)化區(qū)域,厚度明顯減薄,原因是該部分材料承載較小載荷,表現(xiàn)的冗余量較大;隨著去除材料的不斷增大,齒輪重量減小,但齒輪強(qiáng)度也隨之降低,所以后續(xù)應(yīng)進(jìn)行強(qiáng)度校核以滿足安全使用要求。

3.3 模型改進(jìn)設(shè)計(jì)

拓?fù)鋬?yōu)化得到的齒輪模型不能直接進(jìn)行加工,需要對(duì)其進(jìn)行改進(jìn)設(shè)計(jì),以滿足加工制造及生產(chǎn)裝配等要求。

以圖8(c)去除材料質(zhì)量80%時(shí)得到的拓?fù)鋬?yōu)化后齒輪為例,進(jìn)行改進(jìn)設(shè)計(jì)。傳統(tǒng)設(shè)計(jì)的減重孔多為圓柱孔,為保守起見(jiàn),取圖8(c)中最小減重孔靠近輪轂處尺寸為減重孔直徑,圖8(c)的減重孔數(shù)量為10,改進(jìn)設(shè)計(jì)的減重孔仍取同等數(shù)量并沿圓周均布。對(duì)圖8(c)中去除材料質(zhì)量80%時(shí)得到的拓?fù)浣Y(jié)果進(jìn)行改進(jìn)后的齒輪模型見(jiàn)圖9。拓?fù)鋬?yōu)化前主動(dòng)輪質(zhì)量為0.557 35 kg,改進(jìn)后主動(dòng)輪質(zhì)量變?yōu)?.426 34 kg,整體質(zhì)量大約減小23.5%。

圖9 拓?fù)鋬?yōu)化改進(jìn)設(shè)計(jì)齒輪模型

3.4 優(yōu)化后仿真分析

在相同工況下,對(duì)拓?fù)鋬?yōu)化前后齒輪進(jìn)行靜力學(xué)分析,得到的等效應(yīng)力對(duì)比見(jiàn)圖10。

圖10 拓?fù)鋬?yōu)化前、后等效應(yīng)力對(duì)比

由圖10可知,拓?fù)鋬?yōu)化前主動(dòng)輪等效應(yīng)力為219.81 MPa,最大應(yīng)力在齒根處,優(yōu)化后最大等效應(yīng)力值為240.57 MPa,最大應(yīng)力也在齒根處,拓?fù)鋬?yōu)化后相比優(yōu)化前應(yīng)力增大20.76 MPa,應(yīng)力較之前增大9.4%,最大應(yīng)力出現(xiàn)的位置相同,都集中在齒根部位。從應(yīng)力云圖還可以看出,優(yōu)化后輪齒處應(yīng)力集中情況有所改善。優(yōu)化后應(yīng)力值小于齒根彎曲許用應(yīng)力值,滿足齒輪設(shè)計(jì)強(qiáng)度要求,證明拓?fù)涓倪M(jìn)后齒輪結(jié)構(gòu)滿足設(shè)計(jì)要求。

齒輪在嚙合過(guò)程中,輪齒會(huì)產(chǎn)生變形,變形量過(guò)大將會(huì)對(duì)齒輪的運(yùn)行性能造成影響。圖11為拓?fù)鋬?yōu)化前后主動(dòng)齒輪變形云圖。

圖11 拓?fù)鋬?yōu)化前后主動(dòng)齒輪變形云圖

由圖11可知,主動(dòng)齒輪優(yōu)化前后最大變形量分別為0.145 54 、0.200 37 mm.由優(yōu)化前后齒輪變形分析結(jié)果可知,最大變形均出現(xiàn)在齒輪輪齒處,齒輪輻板變形相對(duì)較小,雖然優(yōu)化后的齒輪重量有所減小,最大變形量較優(yōu)化前有所增大,但齒輪輻板在受載時(shí)的最大變形量仍相對(duì)較低,仍能滿足剛度要求。

4 結(jié)論

本文提出了某減速器齒輪的尺寸優(yōu)化和輪輻拓?fù)鋬?yōu)化方法。利用遺傳算法和Ansys聯(lián)合優(yōu)化方法對(duì)齒輪結(jié)構(gòu)進(jìn)行了尺寸優(yōu)化,利用變密度對(duì)齒輪輪輻進(jìn)行了拓?fù)鋬?yōu)化,獲得了材料最佳分布形式。在滿足強(qiáng)度設(shè)計(jì)要求下,實(shí)現(xiàn)減速器齒輪結(jié)構(gòu)輕量化的目的,為齒輪結(jié)構(gòu)設(shè)計(jì)提供了參考。