姜淑鳳 賈瑞超 王俊峰 何鑫林
(①齊齊哈爾大學(xué)機電工程學(xué)院,黑龍江 齊齊哈爾 161006;②齊重數(shù)控裝備股份有限公司,黑龍江 齊齊哈爾 161005;③黑龍江省智能制造裝備產(chǎn)業(yè)化協(xié)同創(chuàng)新中心,黑龍江 齊齊哈爾 161006)
目前,機床行業(yè)正處于轉(zhuǎn)型升級的重要時期,國內(nèi)機床勢必需打破高速化、高效化、高精密化等技術(shù)瓶頸[1]。機床需要升級創(chuàng)新,橫梁作為機床關(guān)鍵支撐件,其性能好壞直接影響著機床性能,為提高機床加工精度首先需對橫梁進(jìn)行優(yōu)化設(shè)計[2]。
近年來,無論國外還是國內(nèi)對機床橫梁優(yōu)化設(shè)計的學(xué)者相對較多。如Zhong H L[3]等人提出一種無支撐件設(shè)計方法,成功實現(xiàn)了機床橫梁的輕量化設(shè)計。Sun Q[4]等人基于拓?fù)鋬?yōu)化方法對焊接梁進(jìn)行優(yōu)化設(shè)計,設(shè)計結(jié)果顯示梁減重10%,焊接梁的綜合性能也得到提高。我國在機床橫梁設(shè)計方面也取得了一定的成績。于英華[5]等人采用新材料填充橫梁結(jié)構(gòu)使橫梁質(zhì)量減輕7.11%;高志來[6]等人基于正交實驗設(shè)計方法,優(yōu)化橫梁結(jié)構(gòu)成功實現(xiàn)橫梁輕量化設(shè)計;大多數(shù)設(shè)計者們主要集中在橫梁筋板結(jié)構(gòu)、截面尺寸上。采用方法主要為拓?fù)鋬?yōu)化、仿生設(shè)計等,對橫梁筋板厚度研究較少,往往忽略橫梁筋板厚度對橫梁特性的影響。
針對以上問題本文采用有限元分析軟件,分析橫梁靜、動態(tài)特性,再針對影響橫梁性能的筋板厚度及薄弱位置進(jìn)行局部結(jié)構(gòu)優(yōu)化,采用層次分析法篩選最優(yōu)設(shè)計方案,從而獲得較合理的橫梁優(yōu)化結(jié)果,為機床優(yōu)化設(shè)計提供新思路。
研究對象為DVT250系列,其結(jié)構(gòu)如圖1所示。橫梁外形尺寸長寬高分別為5 040 mm×100 mm×650 mm。內(nèi)部筋板結(jié)構(gòu)較為復(fù)雜,內(nèi)壁設(shè)有橫梁加強筋板和縱向筋板。中部由6個厚度為20 mm縱向筋板構(gòu)成,兩側(cè)分別由2個厚度為25 mm縱向筋板構(gòu)成。橫梁導(dǎo)軌壁厚為60 mm ,其余壁厚均為25 mm。
橫梁的數(shù)字化模型采用SolidWorks進(jìn)行三維實體建模,為避免計算時間過長,橫梁建模時已忽略非必要特征[7]。橫梁材料采用HT250,材料屬性為彈性模量E=150 GPa,密度ρ=7.25 g/cm3,泊松比μ=0.25。采用1節(jié)點Solid95單元結(jié)構(gòu),劃分方式采用四面體,劃分后模型單元數(shù)為71 440,節(jié)點數(shù)為134 213。橫梁有限元網(wǎng)格劃分模型如圖2所示。
DVT系列車床橫梁的外形由原始設(shè)計數(shù)據(jù)和裝配關(guān)系所決定,其外觀尺寸相對固定,無法采用橫梁外形尺寸為設(shè)計目標(biāo)點,本文采用內(nèi)部筋板厚度為設(shè)計目標(biāo)點[8]。如圖3所示,橫梁承載著滑板部件(包括滑塊、滑枕等零部件),沿橫梁進(jìn)行左右移動。
為研究橫梁受力變形情況,需研究橫梁工況最惡略時的變形即可,由圖可知,當(dāng)滑板部件移動到橫梁中部時工況最惡劣,如圖4所示。N為滑板部件總質(zhì)量;X為橫梁變形處距立柱的距離;m為橫梁質(zhì)量;Ym為中點撓度。
依據(jù)橫梁振動學(xué)原理將橫梁作為矩形截面的簡支梁進(jìn)行彎曲變形計算,梁振動系統(tǒng)的撓曲線近似微分方程為[9]:
(1)
式中:M代表彎矩;ρ代表曲率半徑;E代表彈性模量,I代表截面慣性矩。
對式(1)積分得轉(zhuǎn)角方程為:
(2)
式中:θ代表轉(zhuǎn)角;C代表常數(shù)項。
再將式(2)再積分一次得撓曲線方程為:
(3)
式中:P代表常數(shù)項。
上式計算過程較復(fù)雜,故利用有限元軟件求解,將橫梁與立柱接觸面作為約束,求解總變形云圖如5所示。
橫梁最大變形發(fā)生在橫梁中部,總變形量為5.23 μm。沿X軸變形量為0.96 μm;沿Y軸變形量為0.002 1 μm;沿Z軸變形量為1.73 μm。
模態(tài)分析主要用于確定橫梁振動特性,即固有頻率和振型。固有頻率越大橫梁剛性越大,其振動微分方程為[10]:
(4)
模態(tài)分析為橫梁自由振動狀態(tài),因此需忽略阻尼,式(4)變?yōu)椋?/p>
(5)
式中:[M]、[K]代表常數(shù),自由振動為簡諧運動即式(5)變?yōu)椋?/p>
(6)
式中:ωn代表n階固有頻率;φn代表第n階振型向量。
模態(tài)約束條件與靜態(tài)約束條件相同。求解模態(tài)存在無窮多階模態(tài),前幾階模態(tài)結(jié)果對動態(tài)分析有較大影響,本文對橫梁前3階模態(tài)進(jìn)行著重分析。通過仿真軟件求解得出前3階振型云圖如圖6所示。
圖6分別為1階、2階、3階模態(tài)振型云圖,固有頻率分別為189.17 Hz、202.77 Hz、208.46 Hz。
由上述分析可知,橫梁變形主要發(fā)生在橫梁中部,變形長度約為1 300 mm。原橫梁中部由6個縱向筋版支撐,每個筋板間隔距離均相等為276 mm。橫梁筋板的厚度、數(shù)量及筋板間隔距離都對橫梁性能有一定的影響。為優(yōu)化橫梁特性,需先對橫梁筋板進(jìn)行設(shè)計與分析,其結(jié)果如表1所示。
任何一種筋板厚度、數(shù)量、間隔距離都對于橫梁的不同特性,呈現(xiàn)的變化規(guī)律是不同的。為得到最優(yōu)設(shè)計方案,采用層次分析法對橫梁筋板厚度優(yōu)選方案進(jìn)行分析選取,最終獲得筋板厚度與橫梁不同特性間的權(quán)重值,再依據(jù)權(quán)重值得到關(guān)聯(lián)度值。
表1 不同筋板厚度仿真分析結(jié)果
層次分析法是把與決策有關(guān)的所有因素以目標(biāo)、準(zhǔn)則、方案等分層,分別確定各個因素的權(quán)重值?;疑P(guān)聯(lián)度是依據(jù)各因素與最優(yōu)結(jié)果接近程度進(jìn)行分析[11-13]。本文將橫梁質(zhì)量、形變量、應(yīng)力、前3階固有頻率作為分析因素。
計算橫梁厚度與各因素權(quán)重關(guān)系需將橫梁4個評斷指標(biāo)轉(zhuǎn)化為矩陣I:
再對矩陣I的每列進(jìn)行歸一化處理,得到矩陣U:
歸一化處理計算公式為:
(7)
將式(7)中得到的矩陣U每行進(jìn)行相加得到矩陣V:
(8)
將矩陣V進(jìn)行歸一化,即得到權(quán)重B=(b1,b2,b3,…,b24,b25)T
計算評價指標(biāo)時共有兩種計算公式,對成正比例關(guān)系的評價指標(biāo)如前3階固有頻率計算公式為:
(9)
對成反比例關(guān)系的評價指標(biāo)如質(zhì)量,最大變形量,最大應(yīng)力計算公式為:
(10)
選擇矩陣U中每一列的最大值作為關(guān)聯(lián)系數(shù)的參考數(shù)據(jù)列,其中H=[h1h2h3h4],hβ=max(h1βh2βh3βh4β)。
計算關(guān)聯(lián)系數(shù)矩陣W:
計算公式為:
(11)
計算關(guān)聯(lián)度γβ得:
γβ=(wβ1b1+wβ2b2+wβ3b3+wβ4b4)/4
(12)
在式(7)~(12)中ρ為0.5,α為(1,2,3,4);β為(1,2,3,4)。
橫梁評價指標(biāo)有兩種情況,一種如橫梁前3階固有頻率越大則越優(yōu);另一種情況如橫梁質(zhì)量、最大變形、應(yīng)力越小則越優(yōu)。根據(jù)式(7)~(8)得到權(quán)重B=[1.778 6 1.305 7 1.415 5 2.502 9 2.819 2 0.744 9 1.223 1 1.245 0 2.405 1 2.728 2 1.545 2 1.177 7 1.168 7 2.465 7 2.796 0 1.524 7 1.086 8 1.356 7 2.874 1 2.795 8 1.485 6 1.125 0 1.354 1 2.686 0 2.729 3 ]。根據(jù)式(9)~(11)得到關(guān)聯(lián)系數(shù)矩陣W:
根據(jù)式(12)計算得出每個方案的關(guān)聯(lián)度γ,γ=[0,1],數(shù)值越接近1關(guān)聯(lián)度越高,該橫梁厚度越影響性能指標(biāo),數(shù)值越接近0則相反[14]。為方便方案的選擇,將矩陣γ中數(shù)據(jù)以表格形式呈現(xiàn),如表2所示。
表2 筋板厚度與關(guān)聯(lián)度對比
由表2可知,DVT系列橫梁筋板數(shù)量為5個,厚度為10 mm,間隔距離340 mm時關(guān)聯(lián)度最大為0.929 2,因此選用此筋版作為新橫梁筋板,并與優(yōu)化前橫梁特性值進(jìn)行比較,大于優(yōu)化前橫梁數(shù)據(jù)值為正數(shù),反之為負(fù)數(shù),如表3所示。
表3 橫梁優(yōu)化前后對比
由表3可知,優(yōu)化后橫梁與原橫梁質(zhì)量減少117 kg,占優(yōu)化前橫梁質(zhì)量的2.16%;形變量僅減少0.01 μm,應(yīng)力減少0.250 5 MPa,占優(yōu)化前橫梁應(yīng)力的18%,前3階固有頻率均有所增加。綜上所述,DVT系列橫梁筋板厚度中部以5個縱向筋版筋版厚度10 mm每個筋版間隔距離340 mm厚時設(shè)計為最優(yōu)方案。
(1)采用有限元軟件對橫梁靜、動態(tài)特性進(jìn)行分析,發(fā)現(xiàn)橫梁中部受力變形較大,因此對橫梁中部進(jìn)行優(yōu)化設(shè)計,再結(jié)合DVT橫梁生產(chǎn)制造要求,最終以橫梁中部筋板厚度作為優(yōu)化設(shè)計目標(biāo)點。
(2)采用層次分析法確定筋板個數(shù)、厚度及相隔距離與橫梁不同特性的權(quán)重,再采用灰色關(guān)聯(lián)度法對比,分析各組數(shù)據(jù)的關(guān)聯(lián)度值,最終選用筋板個數(shù)為5厚度為10 mm筋板間距340 mm時為最優(yōu)設(shè)計方案。
(3)優(yōu)化后的橫梁在質(zhì)量、應(yīng)力、形變量方面均得到不同程度的減少,前3階固有頻率均得到一定提升,優(yōu)化設(shè)計效果較顯著。
本文以橫梁質(zhì)量、形變量、應(yīng)力及固有頻率為評估指標(biāo),以內(nèi)部筋板筋板數(shù)量、厚度和筋板間距為試驗因素,運用層次分析法和灰色關(guān)聯(lián)度法分析數(shù)據(jù),從而實現(xiàn)橫梁最優(yōu)設(shè)計方案,為今后機床優(yōu)化設(shè)計提供新的設(shè)計思路。