譚 峰，殷國(guó)富，方 輝，姬坤海，王 亮

（四川大學(xué)制造科學(xué)與工程學(xué)院，成都 610065）

基于ANSYS Workbench的微型數(shù)控車床主軸動(dòng)靜態(tài)性能分析*

譚 峰，殷國(guó)富，方 輝，姬坤海，王 亮

（四川大學(xué)制造科學(xué)與工程學(xué)院，成都 610065）

數(shù)控機(jī)床主軸的結(jié)構(gòu)特性對(duì)其加工精度具有重要影響。根據(jù)自行設(shè)計(jì)的微型數(shù)控車床的主軸結(jié)構(gòu)特征，應(yīng)用有限元分析軟件ANSYS Workbench對(duì)其主軸進(jìn)行了動(dòng)靜態(tài)性能分析。通過對(duì)微型數(shù)控車床主軸的靜力分析、疲勞分析的和模態(tài)分析，得到其動(dòng)靜態(tài)特性參數(shù)。在靜態(tài)特性分析中將軸承約束等效為剛性約束，在模態(tài)分析中將軸承約束等效為彈簧單元的彈性約束，求主軸約束模態(tài)，進(jìn)而求其臨界轉(zhuǎn)速。求得主軸動(dòng)靜態(tài)特性參數(shù)均在合理的范圍內(nèi)，驗(yàn)證了主軸設(shè)計(jì)的合理性，為后續(xù)進(jìn)行微型數(shù)控車床整機(jī)的動(dòng)靜態(tài)特性研究與優(yōu)化設(shè)計(jì)奠定了基礎(chǔ)。

微型數(shù)控車床；主軸；靜力分析；疲勞分析；模態(tài)分析

0 引言

數(shù)控機(jī)床的高精度化是其主要的發(fā)展趨勢(shì)。主軸部件是數(shù)控機(jī)床最為關(guān)鍵的部件之一，其動(dòng)靜態(tài)性能對(duì)機(jī)床的最終加工性能有非常重要的影響。這一影響在刀具切削工件造成的綜合位移影響中所占的比重可達(dá)到60%～80%［1］。因此，隨著機(jī)床速度和精度的提高，對(duì)其主軸部件的動(dòng)靜態(tài)性能提出了更高的設(shè)計(jì)和加工制造要求。各科研院所和研究機(jī)構(gòu)對(duì)主軸部件的動(dòng)靜態(tài)性能進(jìn)行了深入、廣泛的研究［2-4］。

本文以自行設(shè)計(jì)的微型數(shù)控車床為研究對(duì)象，其主要用于加工小、微型零件，具有加工精度高，加工效率高，節(jié)約能源，占地面積少等優(yōu)點(diǎn)。在微型數(shù)控車床的設(shè)計(jì)中，保證主軸部件具有較好的動(dòng)靜態(tài)特性是十分重要的，其靜態(tài)特性（包括強(qiáng)度、剛度和疲勞特性等）和動(dòng)態(tài)特性（模態(tài)特性等）的優(yōu)劣都將直接影響到整臺(tái)機(jī)床的使用性能。靜態(tài)特性直接決定了主軸的使用性能與壽命，動(dòng)態(tài)特性將直接影響主軸抵抗自激振動(dòng)與受迫振動(dòng)的能力，影響主軸的加工精度和工件的表面質(zhì)量，從而影響微型數(shù)控車床的加工性能。為使微型數(shù)控車床主軸具有剛度高、振動(dòng)小等良好性能，需研究其主軸的動(dòng)靜態(tài)特性，改善其薄弱環(huán)節(jié)，減小其對(duì)機(jī)床整機(jī)的動(dòng)靜態(tài)性能的影響。因此，在設(shè)計(jì)階段需對(duì)主軸的動(dòng)靜態(tài)特性進(jìn)行合理而準(zhǔn)確的分析，以提高設(shè)計(jì)效率，減少試驗(yàn)成本，進(jìn)而提高其使用性能。從而對(duì)提高微型數(shù)控車床的設(shè)計(jì)水平具有非常重要的理論和現(xiàn)實(shí)意義。應(yīng)用ANSYS Workbench有限元分析軟件對(duì)微型數(shù)控車床的主軸進(jìn)行動(dòng)靜態(tài)特性分析。Workbench提供了Windows風(fēng)格的友好界面，與Solid-Works等CAD軟件的無縫接口技術(shù)、新一代的參數(shù)化建模工具和領(lǐng)先的優(yōu)化技術(shù)使用戶能夠方便快捷地進(jìn)行CAE分析。

1 微型數(shù)控車床及其主軸結(jié)構(gòu)

以自行設(shè)計(jì)的微型數(shù)控車床的主軸為研究對(duì)象。微型數(shù)控車床結(jié)構(gòu)簡(jiǎn)圖如圖1所示，長(zhǎng)700mm，寬288mm，高233mm。

該機(jī)床的主軸部件結(jié)構(gòu)簡(jiǎn)圖如圖2所示。主軸部件設(shè)計(jì)的轉(zhuǎn)速范圍為400～6000r/min，可實(shí)現(xiàn)無級(jí)調(diào)速，總長(zhǎng)度為216mm。主軸通過左右兩組軸承安裝在主軸箱內(nèi)，主軸軸承采用了洛陽軸承研究所的高精密角接觸球軸承，它具有膨脹系數(shù)小、彈性模量大、極限轉(zhuǎn)速高和抗振動(dòng)性能好的特點(diǎn)。左側(cè)為單列安裝的面向卡盤的角接觸球軸承，型號(hào)為7003AC（α=25°），起到徑向支承的作用；右側(cè)為背靠背安裝的角接觸球軸承，型號(hào)為7004AC/DB（α=25°），承受軸向力和徑向力。這樣的組合保證主軸具有足夠的剛度和回轉(zhuǎn)精度。

圖1 微型數(shù)控車床結(jié)構(gòu)圖

圖2 微型數(shù)控車床主軸結(jié)構(gòu)圖

2 靜態(tài)特性分析

2.1 靜力分析

主軸的靜力分析主要包括強(qiáng)度和剛度的計(jì)算。主軸采用40Cr合金結(jié)構(gòu)鋼，調(diào)質(zhì)處理，屈服強(qiáng)度σs=785MPa。主軸材料屬性如表1所示。

表1 主軸材料屬性

施加約束和載荷時(shí)，應(yīng)盡量按照實(shí)際工況進(jìn)行，這樣才能保證計(jì)算結(jié)果更準(zhǔn)確。根據(jù)實(shí)際工況，加載位置有三處：一是，主軸左側(cè)與帶輪連接的鍵槽，施加的是固定約束；二是，軸承和主軸的接觸，將軸承支撐簡(jiǎn)化為剛性支撐，左側(cè)軸承用圓柱約束（Cylindrical Support）提供徑向支撐，右側(cè)軸承用圓柱約束（Cylindrical Support）提供軸向和徑向支撐［5］；三是，切削過程中產(chǎn)生的切削力和轉(zhuǎn)矩經(jīng)過轉(zhuǎn)化后加載在主軸和卡盤的連接處，但直接施加比較麻煩，因此將這幾處等效為剛性單元，將車刀切削工件產(chǎn)生的切削力直接施加在剛性單元上，主軸輸入端的轉(zhuǎn)矩由負(fù)載決定，因此在主軸和卡盤結(jié)合面施加遠(yuǎn)程力（Remote Force）模擬實(shí)際車刀切削力。這樣的約束和加載方式是符合實(shí)際工況的。

根據(jù)車床切削力計(jì)算公式以及一般外圓切削實(shí)際工況［6］，計(jì)算出主切削力Fc=210N，背向力Fp=126N，進(jìn)給力Ff=105N。在Workbench中通過便利的遠(yuǎn)程力（Remote Force）來模擬三個(gè)正交的切削力，設(shè)置其大小為（-105，210，-126），該遠(yuǎn)程力施力位置為實(shí)際車刀切削位置，相對(duì)主軸左端面坐標(biāo)為（350，0，20）。

有限元分析的仿真精度取決于有限元模型，建立模型時(shí)，忽略倒角、倒圓、螺紋等小特征，簡(jiǎn)化為多階梯空心圓柱體。利用SolidWorks軟件完成主軸建模后，按照Parasolid標(biāo)準(zhǔn)輸出“.x_t”文件，導(dǎo)入到Workbench環(huán)境下進(jìn)行網(wǎng)格劃分。在網(wǎng)格尺寸設(shè)置中將關(guān)聯(lián)中心（Relevance Center）設(shè)置為密網(wǎng)格（Fine），設(shè)置單元尺寸（Element Size）為4mm，并采用默認(rèn)的四面體網(wǎng)格進(jìn)行自由網(wǎng)格劃分。主軸模型網(wǎng)格劃分結(jié)果如圖3所示，共有8791個(gè)單元，15177個(gè)節(jié)點(diǎn)。

圖3 網(wǎng)格劃分后的主軸有限元模型

主軸有限元靜力分析結(jié)果如圖4所示，從圖4a可知主軸最大等效應(yīng)力為73MPa，40Cr合金結(jié)構(gòu)鋼的屈服強(qiáng)度為785MPa，安全系數(shù)大于10，即使考慮應(yīng)力集中的情況，根據(jù)第四強(qiáng)度理論，主軸強(qiáng)度依然滿足要求［7］。從圖4b可知主軸最大等效應(yīng)變?yōu)?.35μm。說明主軸結(jié)構(gòu)在強(qiáng)度與剛度上均達(dá)到了較好的效果。

圖4 主軸有限元靜力分析結(jié)果圖

2.2 疲勞分析

微型數(shù)控車床在工作過程中，主軸旋轉(zhuǎn)，所受的切削力是一個(gè)交變應(yīng)力，在交變應(yīng)力的作用下，雖然主軸所承受的應(yīng)力低于材料的屈服極限，但經(jīng)過較長(zhǎng)時(shí)間的工作后主軸會(huì)產(chǎn)生裂紋或突然發(fā)生完全斷裂，這種現(xiàn)象稱為金屬疲勞，故需對(duì)主軸進(jìn)行疲勞分析。

在對(duì)主軸靜力分析的基礎(chǔ)上，進(jìn)一步對(duì)其進(jìn)行疲勞分析。S/N（應(yīng)力/壽命）曲線是材料疲勞失效時(shí)應(yīng)力幅值S與對(duì)應(yīng)的疲勞壽命N的關(guān)系曲線［8］。設(shè)置主軸材料的S/N曲線，如圖5所示。

圖5 主軸材料的S/N曲線

在Workbench的靜力分析結(jié)果中添加一個(gè)Fatigue Tool（疲勞工具），設(shè)計(jì)壽命設(shè)為1e6，進(jìn)行疲勞分析，得到疲勞結(jié)果如圖6所示，從圖6a可知主軸最低疲勞壽命為1e6，從圖6b可知主軸最低疲勞安全系數(shù)為1.181，最高為15，說明主軸設(shè)計(jì)滿足疲勞壽命要求。

圖6 主軸疲勞分析結(jié)果圖

3 動(dòng)態(tài)特性分析

3.1 模態(tài)分析理論

根據(jù)有限元理論，主軸的動(dòng)力學(xué)方程如下：式中，［M］為主軸質(zhì)量矩陣，［K］為主軸剛度矩陣，［C］為主軸阻尼矩陣，分別為節(jié)點(diǎn)的位移、速度和加速度向量，｛F（t）｝為節(jié)點(diǎn)所受外力向量。固有頻率只與系統(tǒng)本身的特性（質(zhì)量、剛度和阻尼）有關(guān)，模態(tài)分析即是求解振動(dòng)系統(tǒng)的固有頻率和振型［9］。當(dāng)彈性體的動(dòng)力學(xué)基本方程中的外力向量｛F（t）｝=｛0｝時(shí)，略去阻尼，便可得系統(tǒng)的自由振動(dòng)方程：

解得，其特征方程為：

式中，f即為系統(tǒng)的固有頻率。

機(jī)床動(dòng)態(tài)特性是影響機(jī)床性能的重要因素，將直接影響機(jī)床最后的加工性能，是評(píng)定機(jī)床性能的重要指標(biāo)。對(duì)微型數(shù)控車床，其主軸的動(dòng)態(tài)特性對(duì)機(jī)床的加工性能影響很大。因此，對(duì)微型數(shù)控車床主軸的動(dòng)態(tài)特性分析，研究其對(duì)機(jī)床性能的影響就尤為重要。機(jī)床抵抗振動(dòng)能力的大小是評(píng)價(jià)主軸動(dòng)態(tài)性能的重要指標(biāo)。振動(dòng)的幅值與激振力的頻率關(guān)系很大，對(duì)于微型數(shù)控車床，激振力的頻率和振幅隨著轉(zhuǎn)速的提高而增大，很容易接近于系統(tǒng)的固有頻率，當(dāng)主軸的某階固有頻率與激振頻率相等或相近時(shí)，將使振幅劇增，產(chǎn)生共振［10］，因此主軸的各階固有頻率應(yīng)作為主軸不可忽視的一項(xiàng)評(píng)價(jià)指標(biāo)。

主軸的振動(dòng)可以表達(dá)為各階固有振型的線性組合，其中低階固有振型要比高階固有振型對(duì)軸的振動(dòng)影響大，越是低階影響就越大，因此低階振型對(duì)軸的動(dòng)態(tài)特性起決定作用，故在進(jìn)行主軸的模態(tài)分析時(shí)取前6階低階模態(tài)［11］。

3.2 模態(tài)分析

為了更精確地求解主軸在實(shí)際工況約束下的模態(tài)，采用固定剛度的彈簧模擬軸承對(duì)主軸的約束，對(duì)主軸進(jìn)行約束模態(tài)分析。所用求解方法為Workbench默認(rèn)的Block lanczos法，該方法計(jì)算結(jié)果較精確，收斂較快，且計(jì)算速度也快。

軸承的預(yù)緊剛度計(jì)算十分復(fù)雜。預(yù)緊剛度越大，軸的剛性越高，但軸承壽命和最大轉(zhuǎn)速減少。單個(gè)軸承預(yù)緊后的徑向剛度Kr可采用式（5）進(jìn)行計(jì)算［12］。

式中Kr—軸承徑向剛度，N/m；

Z—軸承滾動(dòng)體數(shù)目；

Db—軸承滾動(dòng)體直徑，m；

α—軸承接觸角，°；

Fa0—軸承預(yù)緊力，N。

主軸軸承采用定位預(yù)緊方式，根據(jù)所選軸承的參數(shù)，預(yù)緊力為輕預(yù)緊，左右軸承預(yù)緊力大小分別為50N、80N。由洛陽軸承研究所提供的左右支撐軸承相關(guān)參數(shù)如表2所示。

表2 軸承參數(shù)

根據(jù)軸承參數(shù)和預(yù)緊力以及剛度計(jì)算公式，求算得左側(cè)軸承徑向剛度為6.359e7 N/m，即6.359e4 N/mm；右側(cè)軸承徑向剛度值為7.621e7N/m，即7.621e4N/mm。

用4個(gè)沿圓周方向上分布的彈簧模擬軸承支撐［13］。左側(cè)為單列角接觸球軸承，模擬為單組彈簧，約束在軸頸中點(diǎn)處；右側(cè)為背靠背安裝的雙列角接觸球軸承，模擬為雙組彈簧，分別約束在軸承中點(diǎn)處。軸承外圈全約束，即為彈簧固定端；軸承內(nèi)圈提供徑向支撐，即為彈簧游動(dòng)端。彈簧布置圖如圖7所示。

圖7 彈簧布置圖

在Workbench中用Connections中的Body-Ground中的Spring模擬軸承支撐，Ground對(duì)應(yīng)軸承外圈，全約束，Body對(duì)應(yīng)軸承內(nèi)圈，提供徑向支撐。再在右側(cè)軸承處用圓柱約束（Cylindrical Support）提供軸向約束。輸入彈簧剛度值，對(duì)主軸進(jìn)行約束模態(tài)分析，得到前6階固有頻率如表3所示。

表3 約束模態(tài)主軸前6階頻率

主軸前6階約束模態(tài)振型云圖，如圖8所示。

圖8 主軸約束模態(tài)振型云圖

從圖8可以看出，主軸為軸對(duì)稱結(jié)構(gòu)，第1階振型為主軸徑向伸縮，是因?yàn)橹鬏S大端具有卡盤定位結(jié)構(gòu)，為主軸高階振型。主軸徑向伸縮對(duì)主軸與軸承的配合、間隙以及預(yù)緊具有重要影響。第2、3階頻率很接近，振型為沿正交的兩個(gè)徑向方向的一次彎曲振動(dòng)。第4、5階頻率也很接近，振型為沿正交的兩個(gè)徑向方向的二次彎曲振動(dòng)。第6、7階振型應(yīng)為沿正交的兩個(gè)徑向方向的三次彎曲振動(dòng)。但主軸各階頻率太高，實(shí)際主軸達(dá)不到如此高的頻率。

3.3 臨界轉(zhuǎn)速

主軸在運(yùn)轉(zhuǎn)中都會(huì)發(fā)生振動(dòng)，主軸的振幅隨轉(zhuǎn)速的增大而增大，到某一轉(zhuǎn)速時(shí)振幅達(dá)到最大值（共振），超過這一轉(zhuǎn)速后振幅隨轉(zhuǎn)速增大逐漸減少，且穩(wěn)定于某一范圍內(nèi)，這一主軸振幅最大的轉(zhuǎn)速稱為主軸的臨界轉(zhuǎn)速。這個(gè)轉(zhuǎn)速與主軸的固有頻率相關(guān)。

由約束模態(tài)分析結(jié)果可知，第1階固有頻率為3.3704e-3Hz，約等于0，為剛體運(yùn)動(dòng)，可以忽略；第2階與第3階頻率值很接近，并且振型表現(xiàn)為正交；第4階與第5階，第6階與第7階頻率值也很接近，并且為正交振型。由于主軸設(shè)計(jì)的工作轉(zhuǎn)速在6000r/min以下，屬于中低速范圍，轉(zhuǎn)速對(duì)主軸固有頻率影響不大，所以忽略轉(zhuǎn)速對(duì)主軸臨界轉(zhuǎn)速的影響［14］。根據(jù)主軸模態(tài)分析得到的固有頻率由式（6）計(jì)算主軸的臨界轉(zhuǎn)速［15］。

式中：n—臨界轉(zhuǎn)速，r/min；f—固有頻率，Hz。

主軸約束模態(tài)第2階臨界轉(zhuǎn)速n2=60×3261.3=195678r/min，遠(yuǎn)遠(yuǎn)高于該主軸的最高工作轉(zhuǎn)速6000 r/min，表明設(shè)計(jì)的主軸工作轉(zhuǎn)速在安全范圍內(nèi)，說明該主軸設(shè)計(jì)是合理的，能有效地避開共振區(qū)，可保證主軸的加工精度。

4 結(jié)論

主軸系統(tǒng)是機(jī)床最為關(guān)鍵的系統(tǒng)之一，對(duì)其進(jìn)行動(dòng)靜態(tài)特性研究對(duì)提高微型數(shù)控車床整機(jī)的性能至關(guān)重要。本文利用ANSYS Workbench有限元分析軟件建立了主軸的動(dòng)靜態(tài)性能分析模型，對(duì)微型數(shù)控車床的主軸進(jìn)行了靜力分析、疲勞分析以及模態(tài)分析，在模態(tài)分析中采用彈簧單元模擬軸承支承，得到更精確的模態(tài)分析結(jié)果。驗(yàn)證了主軸建模及設(shè)計(jì)的合理性，為進(jìn)一步優(yōu)化主軸系統(tǒng)結(jié)構(gòu)設(shè)計(jì)提供了理論依據(jù)，為深入研究微型數(shù)控車床整機(jī)的動(dòng)靜態(tài)特性奠定了基礎(chǔ)，同時(shí)也為實(shí)際試驗(yàn)提供了參考和依據(jù)。

［1］張耀滿，劉春時(shí)，謝志坤，等.高速機(jī)床主軸部件有限元分析［J］.東北大學(xué)學(xué)報(bào)：自然科學(xué)版，2008，29（10）：1474-1477.

［2］杜官將，李東波.基于ANSYS的機(jī)床主軸結(jié)構(gòu)優(yōu)化設(shè)計(jì)［J］.組合機(jī)床與自動(dòng)化加工技術(shù)，2012（12）：22-24.

［3］張學(xué)玲，滿佳，柴樹峰，等.基于模態(tài)分析的內(nèi)孔磨床主軸動(dòng)剛度改善措施［J］.組合機(jī)床與自動(dòng)化加工技術(shù)，2013（11）：95-98.

［4］孫偉，汪博，聞邦椿.高速主軸系統(tǒng)靜止及運(yùn)轉(zhuǎn)狀態(tài)下動(dòng)力學(xué)特性對(duì)比分析［J］.機(jī)械工程學(xué)報(bào)，2012，48（11）：146-152.

［5］張耀滿，王偉，劉永賢.數(shù)控車床主軸部件有限元分析及其驗(yàn)證［J］.東北大學(xué)學(xué)報(bào) （自然科學(xué)版），2009，30（5）：720-723.

［6］王杰，李方信，肖素梅.機(jī)械制造工程學(xué)［M］.北京郵電大學(xué)出版社，2004.［7］方鵬，李健，韋遼.基于ANSYS Workbench機(jī)床主軸有限元分析［J］.裝備制造技術(shù)，2013（4）：28-30.

［8］王海軍，殷國(guó)富，何波，等.基于SolidWorks Simulation的球磨機(jī)轉(zhuǎn)軸有限元分析［J］.機(jī)械，2014，41（1）：44-47.

［9］凌桂龍，丁金濱，溫正.ANSYS Workbench 13.0從入門到精通［M］.北京：清華大學(xué)出版社，2012.

［10］張耀滿，王旭東，蔡光起，等.高速機(jī)床有限元分析及其動(dòng)態(tài)性能試驗(yàn)［J］.組合機(jī)床與自動(dòng)化加工技術(shù)，2004（12）：15-17.

［11］楊明亞，楊濤，湯本金，等.應(yīng)用有限元分析系統(tǒng)計(jì)算車床主軸的動(dòng)態(tài)特性［J］.機(jī)械工程與自動(dòng)化，2007，1（3）：41-44.

［12］戴曙.機(jī)床滾動(dòng)軸承應(yīng)用手冊(cè)［Z］.北京：機(jī)械工業(yè)出版社，1993.

［13］巫少龍，張?jiān)?基于ANSYS Workbench的高速電主軸動(dòng)力學(xué)特性分析［J］.組合機(jī)床與自動(dòng)化加工技術(shù)，2010（9）：20-22.

［14］穆楊，羅平，劉恒，等.車銑復(fù)合加工中心電主軸的穩(wěn)定性分析和實(shí)驗(yàn)研究［J］.組合機(jī)床與自動(dòng)化加工技術(shù)，2011（3）：4-7.

［15］梁雙翼，陳晨，尹輝俊.基于有限元方法的電主軸模態(tài)分析［J］.制造業(yè)自動(dòng)化，2012（2）：57-60.

The Finite Element Analysis of the Desktop CNC Lathe’s Spind le Based on the ANSYS Workbench

TAN Feng，YIN Guo-fu，F(xiàn)ANG Hui，JIKun-hai，WANG Liang
（School of Manufacturing Science and Engineering，Sichuan University，Chengdu 610065，China）

The structure characteristics of the CNC machine tool’s spindle has important influence on its machining accuracy.According to the spindle structure characteristics of our own designed micro CNC lathe，the spindle’s static and dynamic performance was carried out through the FEA software ANSYS Workbench. Through the static analysis，fatigue analysis and modal analysis of the micro CNC lathe’s spindle，the static and dynamic characteristic parameters were obtained.The bearing constraints are equivalent to rigid constraints and elastic constraints respectively during the static characteristic analysis and modal analysis.After got the spindle’s constraint modal，the critical speed was calculated.The spindle’s static and dynamic characteristic parameters are in reasonable range，so it verifies the rationality of the spindle design and lays a foundation for subsequent static and dynamic characteristics study and optimal design of the overall micro CNC lathe.

micro CNC lathe；spindle；static analysis；fatigue analysis；modal analysis

TH16；TG659

A

1001-2265（2015）04-0029-04 DOI：10.13462/j.cnki.mmtamt.2015.04.007

2014-07-04

四川省科技支撐計(jì)劃項(xiàng)目資助（2014GZX0001）

譚峰（1990—），男，重慶萬州人，四川大學(xué)碩士研究生，研究方向?yàn)闄C(jī)床動(dòng)力學(xué)，（Email）306413664@qq.com；通訊作者：殷國(guó)富（1956—），男，成都人，四川大學(xué)教授，博士生導(dǎo)師，研究方向?yàn)橹圃熳詣?dòng)化，智能設(shè)計(jì)技術(shù)。CAD/CAM/CIMS，（E-mail）gfyin@scu. edu.cn。