季彬彬，楊玉萍，莫亞梅,陳寶國(guó)，董在俊

(1.南通大學(xué) 機(jī)械工程學(xué)院，江蘇 南通 226019；2.江蘇萬(wàn)達(dá)特種軸承有限公司，江蘇 南通 226500)

浪形保持架的疲勞斷裂已成為深溝球軸承在使用過(guò)程中失效的重要形式之一。理論上在旋轉(zhuǎn)一周的運(yùn)動(dòng)中，鋼球和保持架兜孔間的作用力改變一次方向，即在承載區(qū)內(nèi)鋼球推動(dòng)保持架運(yùn)動(dòng)，而在非承載區(qū)內(nèi)鋼球被保持架推動(dòng)，這使鋼球與保持架的梁之間不僅產(chǎn)生摩擦、磨損，還至少遭遇一次周向的變向沖擊[1],所以在保持架內(nèi)部產(chǎn)生應(yīng)力的載荷主要來(lái)自保持架與鋼球的碰撞、套圈的約束以及離心載荷(高速場(chǎng)合)[2]。

關(guān)于浪形保持架疲勞斷裂的研究已經(jīng)取得了一定成果。這些成果都是圍繞保持架在不同工況運(yùn)行過(guò)程中所承受的最大載荷展開(kāi)，通過(guò)控制保持架的最大載荷來(lái)控制最大應(yīng)力，從而提高保持架的疲勞壽命。但保持架的最大應(yīng)力除了與最大載荷有關(guān)外，還與保持架的結(jié)構(gòu)參數(shù)如保持架兜孔與梁連接處的相對(duì)彎曲半徑等有關(guān)。

下文擬采用有限元法分析深溝球軸承浪形保持架兜孔與梁連接處取不同相對(duì)彎曲半徑時(shí)圓角處的應(yīng)力分布狀況。同時(shí)，考慮到浪形保持架是沖壓成形件，其在成形過(guò)程中截面厚度已不是嚴(yán)格意義上的均勻分布，故在應(yīng)力分析時(shí)，將保持架過(guò)渡圓角處的截面厚度分布數(shù)據(jù)融合進(jìn)保持架應(yīng)力分析的有限元模型中，從而得到最接近保持架圓角處實(shí)際應(yīng)力的應(yīng)力分布數(shù)值，為保持架的結(jié)構(gòu)設(shè)計(jì)提供依據(jù)。

1 保持架載荷的確定

作用在保持架上的力主要是鋼球與保持架兜孔之間的作用力以及保持架引導(dǎo)面與引導(dǎo)擋邊之間的作用力。一個(gè)鋼球可與兜孔兩邊任意一邊接觸，這取決于是保持架引導(dǎo)還是鋼球引導(dǎo)[3]。文獻(xiàn)[4]分析了徑向載荷下高速軸承保持架兜孔與鋼球間作用力，其大小隨時(shí)間而呈周期變化，且峰值為10～30 N。下文以中低速深溝球軸承為例，在分析保持架圓角處應(yīng)力分布時(shí)，取兜孔與鋼球間的作用力峰值為10 N，著重比較保持架兜孔與梁處不同相對(duì)彎曲半徑對(duì)圓角處應(yīng)力分布的影響(兜孔與鋼球間作用力峰值的取值對(duì)分析結(jié)果影響不大)。

2 保持架沖壓成形模擬

保持架的沖壓成形模擬可用多種軟件實(shí)現(xiàn)，如基于靜力隱式的AutoForm和基于動(dòng)力顯式的LS-DYNA3D，eta/DYNAFORM等。

現(xiàn)采用eta/DYNAFORM進(jìn)行保持架的沖壓成形模擬。eta/DYNAFORM是用于板料成形數(shù)值模擬的專用軟件，可顯著減少模具開(kāi)發(fā)設(shè)計(jì)時(shí)間及試模周期，不但易用性良好，而且包括大量的智能化自動(dòng)工具，可方便地求解各類板料成形問(wèn)題。采用eta/DYNAFORM的多工步成形模擬功能完成保持架從板料到壓彎以及最后整形的全部工藝流程[5]。模擬分析分兩組進(jìn)行，板料初始截面厚度為1.5 mm，一組保持架兜孔與梁處的相對(duì)彎曲半徑取0.8，即過(guò)渡圓角半徑r=1.2 mm；另一組相對(duì)彎曲半徑取1.0，即過(guò)渡圓角半徑r=1.5 mm。保持架材料為08鋼，整形工序完成后保持架的截面厚度分布分別如圖1、圖2所示。

從圖1、圖2可知，當(dāng)保持架過(guò)渡圓角半徑r=1.2mm時(shí)，圓角處最小厚度值為1.137 mm，最大減薄率為24.2%；當(dāng)過(guò)渡圓角半徑r=1.5 mm時(shí)，圓角處最小厚度值為1.145 mm，最大減薄率為23.6%。

圖1 保持架截面厚度分布(r=1.2 mm)

圖2 保持架截面厚度分布(r=1.5 mm)

3 保持架應(yīng)力分析

3.1 應(yīng)力分析的有限元模型

保持架應(yīng)力分析采用通用的應(yīng)力應(yīng)變分析軟件ANSYS完成，在ANSYS中建立保持架有限元網(wǎng)格模型。保持架網(wǎng)格采用Thin shell 181單元，網(wǎng)格大小取1 mm，并在過(guò)渡圓角處作網(wǎng)格細(xì)化[6-7]。為簡(jiǎn)化分析，取保持架其中一個(gè)兜孔及其相鄰的梁作加載分析。分析過(guò)程中在兜孔兩邊梁上施加位移約束，在兜孔靠近梁一側(cè)且處于中間位置的一網(wǎng)格單元上施加壓力載荷，壓力總值為10 N。施加約束和載荷的兜孔及梁有限元網(wǎng)格模型如圖3所示。

圖3 保持架有限元模型施加載荷和約束示意圖

3.2 截面厚度均勻的保持架應(yīng)力分析

不考慮保持架沖壓成形過(guò)程中截面厚度的變化，將板料作等厚度處理，分別分析過(guò)渡圓角半徑r=1.2 mm和1.5 mm時(shí)保持架某一兜孔承受壓力載荷時(shí)過(guò)渡圓角處的Von Mises應(yīng)力。保持架梁與兜孔過(guò)渡圓角處的應(yīng)力分布如圖4、圖5所示。為便于統(tǒng)計(jì)分析，將過(guò)渡圓角處的應(yīng)力分布以網(wǎng)格為單位劃分成不同的應(yīng)力單元，網(wǎng)格不同顏色代表該應(yīng)力單元處的應(yīng)力大小。

圖4 保持架截面厚度均勻時(shí)的圓角處應(yīng)力分布(r=1.2 mm)

圖5 保持架截面厚度均勻時(shí)的圓角處應(yīng)力分布(r=1.5 mm)

從圖4、圖5可統(tǒng)計(jì)出過(guò)渡圓角區(qū)域的應(yīng)力分布狀況，如表1所示。

表1 截面厚度均勻的保持架圓角處應(yīng)力分布狀況

3.3 考慮截面厚度變化的保持架應(yīng)力分析

保持架在加工成形后其截面厚度分布實(shí)際是不均勻的，所以在分析保持架應(yīng)力時(shí)應(yīng)考慮到截面厚度不均勻?qū)?yīng)力分布的影響。

目前可通過(guò)以下途徑將其截面厚度變化融合到保持架應(yīng)力分析中。

(1)在ANSYS/LS-DYNA中完成保持架沖壓成形模擬，然后通過(guò)顯式-隱式序列求解將顯式分析的結(jié)果導(dǎo)入隱式分析中，從而完成保持架的應(yīng)力分析；但是在ANSYS/LS-DYNA中做沖壓成形不如eta/DYNAFORM方便[8]。

(2)在eta/DYNAFORM中完成保持架沖壓成形模擬，將分析得到的保持架截面厚度數(shù)據(jù)映射成nstran結(jié)果文件輸出，然后在ANSYS中讀取nstran文件并在此基礎(chǔ)上完成保持架的應(yīng)力分析，但是ANSYS中讀取nstran文件需對(duì)ANSYS作二次開(kāi)發(fā)[9]，比較繁瑣。

(3)在eta/DYNAFORM中完成保持架沖壓成形模擬，得到包含保持架截面厚度數(shù)據(jù)的dynain結(jié)果文件，再將該數(shù)據(jù)映射到ANSYS做應(yīng)力分析的保持架網(wǎng)格模型中，然后在ANSYS中加載分析，得到保持架應(yīng)力分布[10]。

現(xiàn)采用第3種方法得到考慮保持架截面厚度變化的保持架圓角處應(yīng)力分布，如圖6、圖7所示。

圖6 保持架截面厚度不均勻時(shí)的圓角處應(yīng)力分布(r=1.2 mm)

圖7 保持架截面厚度不均勻時(shí)的圓角處應(yīng)力分布(r=1.5 mm)

從圖6、圖7可統(tǒng)計(jì)出過(guò)渡圓角區(qū)域的應(yīng)力分布狀況，如表2所示。

表2 截面厚度變化的保持架圓角處應(yīng)力分布狀況

3.4 分析結(jié)果對(duì)比

將以上4種分析結(jié)果進(jìn)行對(duì)比，可得到如圖8所示的應(yīng)力分布狀況曲線圖。

圖8 應(yīng)力分布狀況曲線圖

由表1和圖8可知，在其截面厚度均勻的情況下，浪形保持架過(guò)渡圓角半徑r=1.2mm時(shí)，該圓角區(qū)域處占比最大的應(yīng)力區(qū)域?yàn)?2.2～29.6 MPa；而r=1.5 mm時(shí)，圓角區(qū)域處占比最大的應(yīng)力區(qū)域降至14.8～22.2 MPa。且模擬結(jié)果顯示，兩種不同圓角半徑時(shí)各圓角處的最大應(yīng)力值區(qū)域雖然都為29.6～37.0 MPa，但在r=1.2 mm的圓角區(qū)域，29.6～37.0 MPa的最大應(yīng)力值區(qū)域占26.8%；而在r=1.5 mm的圓角區(qū)域，29.6～37.0 MPa的最大應(yīng)力值區(qū)域僅占1.6%，故r=1.5 mm即浪形保持架相對(duì)圓角半徑為1.0時(shí)，其圓角處的應(yīng)力狀況有較大改善。

由表2和圖8可知，考慮由沖壓引起的圓角處截面厚度減薄的影響，對(duì)于r=1.2 mm的浪形保持架，其圓角處占比最大的應(yīng)力區(qū)域從均勻截面厚度時(shí)的22.2～29.6 MPa增至29.6～37.0 MPa；而最大應(yīng)力值區(qū)域也從均勻截面厚度時(shí)的29.6～37.0 MPa增至大于44.4 MPa，增幅約33%；而對(duì)于r=1.5 mm的浪形保持架，截面厚度均勻與考慮截面厚度減薄兩種情況下的應(yīng)力分布狀況大體相似。

綜上所述，在其截面厚度減薄的影響下，r=1.2 mm的浪形保持架其圓角區(qū)域占比最大的應(yīng)力區(qū)域?yàn)?9.6～37.0 MPa，最大應(yīng)力值區(qū)域?yàn)榇笥?4.4 MPa；而r=1.5 mm的浪形保持架其圓角區(qū)域占比最大的應(yīng)力區(qū)域?yàn)?4.8～22.2 MPa，最大應(yīng)力值區(qū)域?yàn)?9.6～37.0 MPa，故r=1.5 mm的保持架其圓角區(qū)域占比最大的應(yīng)力區(qū)域的應(yīng)力值相對(duì)于r=1.2 mm的保持架降低了約40%，整個(gè)圓角區(qū)域的最大應(yīng)力值也降低了約40%。

4 結(jié)論

(1)eta/DYNAFORM是適用于各種沖壓成形分析的專用軟件，其專門(mén)優(yōu)化了沖壓成形的模擬設(shè)置過(guò)程。ANSYS在做靜力隱式分析方面有其特有的優(yōu)勢(shì)，將這兩種軟件結(jié)合起來(lái)可使沖壓成形件得到更接近于真實(shí)狀況的應(yīng)力分析結(jié)果。

(2)在不同的相對(duì)彎曲半徑情況下，由于沖壓成形引起的保持架截面厚度減薄的影響，保持架過(guò)渡圓角處的最大應(yīng)力值與采用均勻截面厚度分析所得到的最大應(yīng)力值最多相差33%左右。

(3)浪形保持架梁與兜孔處相對(duì)彎曲半徑由0.8增至1.0有利于改善圓角處應(yīng)力分布狀況，提高保持架的疲勞壽命。