孫文輝，李天箭

（200082 上海市 上海理工大學(xué) 機(jī)械工程學(xué)院）

0 引言

隨著現(xiàn)代制造業(yè)的高速發(fā)展,機(jī)床各個部件都在往高精度與高性能方向發(fā)展[1]，不僅要求機(jī)床部件的精度要高，還要求其性能更好[2]。機(jī)床在工作時，立柱對整個機(jī)床的模態(tài)、剛度都起到很重要的作用，還能起到支撐整個機(jī)床結(jié)構(gòu)的作用，因此立柱的剛度以及固有頻率大小對機(jī)床來說很重要。

機(jī)床的性能固然重要，而對于機(jī)床部件的設(shè)計，大多數(shù)設(shè)計者還都只停留在經(jīng)驗設(shè)計與結(jié)構(gòu)類比的階段[3-4]，這可能不會使設(shè)計部件的性能提高更多，拓?fù)鋬?yōu)化技術(shù)可以目標(biāo)性地提高結(jié)構(gòu)的性能，能給零件的結(jié)構(gòu)設(shè)計提供有效的參考，在機(jī)床的零部件結(jié)構(gòu)的設(shè)計中已經(jīng)得到廣泛應(yīng)用[5]。拓?fù)鋬?yōu)化技術(shù)可以使構(gòu)型內(nèi)部的力流更加有效地傳遞，還可以減少材料屬性，縮減成本，所以從輕量化與經(jīng)濟(jì)效益方面著手[6]，拓?fù)鋬?yōu)化技術(shù)有一定的優(yōu)勢。

在本文中利用拓?fù)鋬?yōu)化技術(shù)，通過有限元軟件ANSYS Workbench 對立柱進(jìn)行拓?fù)鋬?yōu)化。在考慮到機(jī)床工作時，立柱與主軸箱的連接部位會有多種工況產(chǎn)生，這時要分3 種不同情況分析，從而在以提升一階固有頻率為目標(biāo)的前提下，得出理想的立柱構(gòu)型。

1 立柱有限元模型建模

機(jī)床由很多部件組成，其中包括床身、立柱、主軸箱、主軸、工作臺等幾大部分。本文主要研究機(jī)床立柱的性能。一般機(jī)床的內(nèi)部都是空心的，在本文中建立的有限元模型是實體的結(jié)構(gòu)，對建立的實體的模型進(jìn)行拓?fù)鋬?yōu)化，得出一定的構(gòu)型，對新的模型重新計算與重構(gòu)，從而設(shè)計出滿足一定要求的立柱模型。立柱的有限元模型實體結(jié)構(gòu)如圖1 所示。此模型為長360 mm，寬460 mm，高120 mm，且主軸箱位于上工況位置時的一種情況，連接部位簡化為上工況受力面。

圖1 立柱實體結(jié)構(gòu)Fig.1 Column structure

建立有限元模型后，將立柱實體結(jié)構(gòu)模型導(dǎo)入ANSYS Workbench 中，將材料設(shè)置為鑄鐵。鑄鐵的密度為ρ=7.2×103kg/m3，泊松比為μ=0.28，彈性模量為E=110 GPa。將立柱模型劃分單元網(wǎng)格，單元數(shù)25 254 個，節(jié)點數(shù)109 966 個。進(jìn)行模態(tài)分析，模態(tài)分析時將底面約束面固定住，經(jīng)過分析可以得出此模型的前6 階固有頻率。各階固有頻率如表1 所示。

表1 立柱前6 階模態(tài)的固有頻率Tab.1 Natural frequencies of the first six modes of column

其中1 階振型如圖2 所示。

圖2 1 階振型Fig.2 First mode

2 拓?fù)鋬?yōu)化

2.1 拓?fù)鋬?yōu)化數(shù)學(xué)模型建立

拓?fù)鋬?yōu)化的方法有很多種，本文采用的是變密度法對模型進(jìn)行拓?fù)鋬?yōu)化，以最大剛度為目標(biāo)，提升一階固有頻率。

材料宏觀物理量與材料的密度之間存在非線性關(guān)系

式中：Ei——單元i 彈性模量；μ0——材料的初始泊松比；μ——泊松比；p——密度懲罰因子，常取p=3；ρi——單元相對密度。

變密度法的數(shù)學(xué)模型如式（2）-式（3）：

式中：ρ——單元密度單位向量；n——結(jié)構(gòu)里面單元的數(shù)量；K——結(jié)構(gòu)總剛度矩陣；F——節(jié)點處載荷；U——位移產(chǎn)生的向量；C——結(jié)構(gòu)柔度值；ui——單元i 節(jié)點位移向量；k0——剛度矩陣常量；ρm——相對密度最小的變量；λ1——1 階固有頻率；λ*——原始1 階固有頻率；d1——最大變形位移；d*——原始最大的變形。

2.3 三種工況下拓?fù)鋬?yōu)化

機(jī)床工作時，立柱部分既要考慮底部與底座相接觸部分，也要考慮側(cè)面與主軸箱接觸部分。底面部分可以采用固定端約束，而側(cè)面則需考慮其受力情況。而在機(jī)床工作時往往會遇到多種加工工況，這樣就要考慮機(jī)床在多種工況下拓?fù)鋬?yōu)化的情形，這里采用立柱部分與主軸箱部分連接在一起的情況來進(jìn)行研究。由于機(jī)床運轉(zhuǎn)的多樣性，采用主軸箱位置在上工況位置、中工況位置與下工況位置3 種情況分別進(jìn)行拓?fù)鋄7]，探究立柱的拓?fù)浣Y(jié)構(gòu)變化與規(guī)律。3 種工況如圖3 所示。

圖3 上、中、下三種工況Fig.3 Upper,middle and lower working conditions

對3 個模型分別添加載荷與約束，如圖4 所示。C，D 四個面受力大小方向一致。A 為向下的重力作用，B 為主軸下端面模擬加工受的三項力，C 為多工況運行軌道受三項力，D 面為固定約束面。

圖4 模型受力Fig.4 Model stress

分別進(jìn)行35，35，31 次迭代之后得出其拓?fù)浣Y(jié)果，如圖5 所示。

圖5 三種工況拓?fù)浣Y(jié)果Fig.5 Topological results of three working conditions

2.4 立柱計算外力與單獨拓?fù)?/h3>

上述情況為多工況情況，是由立柱加主軸箱兩部分結(jié)合在一起進(jìn)行分析，現(xiàn)就立柱部分在拆開之后，需要考慮力通過主軸箱傳遞給立柱[7]，力的大小在立柱右側(cè)面的分布情況進(jìn)行研究，在其他面受力不變情況下，探討主軸箱與立柱接觸部分的力傳導(dǎo)情況。以上工況為例，即主軸箱與立柱接觸面在立柱右側(cè)的最上端，其受力情況如圖6 所示。

圖6 立柱受力情況Fig.6 Stress condition of column

圖6 中E 受力面為上工況時主軸箱與立柱的接觸面，機(jī)床工作時，主軸箱上的力對此面產(chǎn)生彎矩作用，此彎矩如圖6 所示，E 受力面添加一個逆時針的彎矩作用。而這個面的受力大小分布可能會不均勻，現(xiàn)對其D 受力線與C 受力線即上下兩部分進(jìn)行探討，研究兩個部分的力大小分布權(quán)重關(guān)系。

現(xiàn)對D 受力線與C 受力線的權(quán)重做比較。原先合力的大小為定值，現(xiàn)改變兩者大小比例值，如表2 所示。其中XD/YD 與XC/YC 分別為D 受力線與C 受力線在X 和Y 方向受力大小。

表2 計算外力Tab.2 Calculation of external force

表2 數(shù)據(jù)中，力的單位為kN。XD/YD與XC/YC力大小的總和為100 kN，現(xiàn)討論兩個部位的力占比情況與拓?fù)浣Y(jié)果做對比，其各個受力大小不同情況拓?fù)涑鼋Y(jié)果如圖7 所示。

圖7 立柱拓?fù)浣Y(jié)果Fig.7 Column topology results

對這5 種拓?fù)浣Y(jié)果分別進(jìn)行后處理與重新計算其變形情況，其5 種情況最大整體變形如表3所示。

表3 不同情況變形大小Tab.3 Deformation size under different conditions

由其結(jié)果可以看出，比例為0.2/0.8 時的變形為0.287 mm，其整體變形量比其他4 種結(jié)果變形較小，則此受力比例較其他4 種較優(yōu)。此拓?fù)淠Ｐ偷南旅婢€受力更大，上面線的受力較小。經(jīng)過后處理重新分析此模型，得出應(yīng)力與整體變形大小，其靜應(yīng)力與整體變形大小如圖8 所示。

圖8 靜應(yīng)力與變形Fig.8 Static stress and deformation

由結(jié)果可知，其最大等效應(yīng)力為22.64 MPa，整體變形為0.287 18 mm。

對此拓?fù)淠Ｐ瓦M(jìn)行模態(tài)分析，可以得出前6階固有頻率，如表4 所示。

表4 前6 階模態(tài)的固有頻率Tab.4 Natural frequencies of the first six modes

由表4 可以看出，拓?fù)淠Ｐ偷囊浑A固有頻率為177.1 Hz。拓?fù)淠Ｐ鸵浑A振型由圖9 所示。

圖9 拓?fù)淠Ｐ? 階振型Fig.9 First mode shape of topological model

3 模型重構(gòu)

根據(jù)立柱拓?fù)浣Y(jié)果初步重構(gòu)出其模型，如圖10 所示。

圖10 立柱重構(gòu)模型Fig.10 Column reconstruction model

對重構(gòu)出的模型重新進(jìn)行分析，得到其應(yīng)力、變形與應(yīng)變，如圖11 所示。

圖11 重構(gòu)模型的應(yīng)力與變形Fig.11 Stress and deformation of reconstructed model

由重新分析可以得出最大等效應(yīng)力大小為28.743 MPa，其位于立柱與底座接觸結(jié)合部右側(cè)；整體變形最大值為0.254 34 mm；應(yīng)變與應(yīng)力成正比，大小為0.000 265 64 mm。可以看出整體變形比拓?fù)浣Y(jié)果小，其剛度性能良好。

將重構(gòu)的模型底面固定后進(jìn)行模態(tài)分析，分析出重構(gòu)模型的前6 階固有頻率，如表5 所示。

表5 前6 階模態(tài)的固有頻率Tab.5 Natural frequencies of the first six modes

由模態(tài)分析得出，前6 階固有頻率，其中1階固有頻率為237.3 Hz。重構(gòu)模型一階振型如圖12 所示。

圖12 重構(gòu)模型1 階振型Fig.12 The first mode shape of reconstructed model

4 拓?fù)淠Ｐ团c重構(gòu)模型模態(tài)對比

原始實體模型為149.8 Hz，拓?fù)淠Ｐ蜑?77.1 Hz，而重構(gòu)模型的一階固有頻率為237.2 Hz，可以看出重構(gòu)模型與拓?fù)淠Ｐ偷? 階固有頻率都高于原來的拓?fù)渲暗? 階固有頻率?，F(xiàn)將拓?fù)淠Ｐ团c重構(gòu)模型的前6 階固有頻率進(jìn)行對比，如圖13 所示。

圖13 前6 階模態(tài)對比Fig.13 Comparison of the first six modes

重構(gòu)模型的前6 階固有頻率總體上要略高于拓?fù)淠Ｐ偷那? 階固有頻率。在機(jī)床運轉(zhuǎn)時，最主要看1 階固有頻率的大小，1 階固有頻率越大，那么它在工作時的性能就越好。拓?fù)淠Ｐ团c重構(gòu)模型1 階固有頻率都大于原始實體模型的1 階固有頻率，總體上性能得以提升。

5 結(jié)論

利用拓?fù)鋬?yōu)化方法探究了立柱結(jié)構(gòu)型與性能，對一個立柱的實體結(jié)構(gòu)進(jìn)行拓?fù)?，拓?fù)溥^程中考慮到多種加工工況，對上中下3 種工況進(jìn)行拓?fù)鋬?yōu)化，分別得到其構(gòu)型。對立柱單獨優(yōu)化考慮了力通過主軸箱部分傳遞給立柱的分布情況，分為5 種情況，最終以變形大小選擇較優(yōu)模型。通過此法得到的模型與多種加工工況得到的模型類比，重構(gòu)一個模型。將他們的性能對比，發(fā)現(xiàn)重構(gòu)模型與拓?fù)淠Ｐ偷?階固有頻率都有所提升，重構(gòu)模型的前6 階固有頻率更好，而且重構(gòu)模型的總體變形也比拓?fù)淠Ｐ偷牡?，則此模型較優(yōu)。