考慮粗糙表面影響因素的齒輪齒根處應(yīng)力集中系數(shù)計(jì)算是齒輪彎曲疲勞壽命精準(zhǔn)預(yù)估的難題，以磨削噴丸后的直齒輪為研究對(duì)象，研究粗糙表面下的齒根應(yīng)力與應(yīng)力集中計(jì)算問題。采用白光干涉儀 Wyko NT9100 對(duì)磨削噴丸齒根粗糙表面進(jìn)行測(cè)量，得到粗糙表面形貌數(shù)據(jù)，基于空間坐標(biāo)變換原理，使用 Python 對(duì)有限元軟件進(jìn)行二次開發(fā)，通過調(diào)整齒根細(xì)化網(wǎng)格節(jié)點(diǎn)坐標(biāo)實(shí)現(xiàn)了齒根表面粗糙形貌的添加工作，建立了齒根過渡曲面三維粗糙表面有限元模型并進(jìn)行仿真分析。通過數(shù)值計(jì)算得到三維粗糙表面參數(shù)下的齒根應(yīng)力分布與應(yīng)力集中系數(shù)，對(duì)粗糙表面參數(shù)與應(yīng)力集中系數(shù)的關(guān)聯(lián)規(guī)律進(jìn)行非線性回歸分析，建立粗糙表面參數(shù)與應(yīng)力集中系數(shù)的關(guān)聯(lián)規(guī)律。結(jié)果表明，粗糙表面參數(shù) Sa 、Sv 、S10z 擬合應(yīng)力集中系數(shù)得到擬合公式的相關(guān)系數(shù)分別為 0.799，0.784，0.914，十點(diǎn)區(qū)域高度參數(shù) S10z 能較好地表征齒根表面的應(yīng)力集中。

　　齒輪作為傳動(dòng)系統(tǒng)的核心部件，被廣泛運(yùn)用于航空航天、汽車坦克和高精密儀器制造等領(lǐng)域，噴丸強(qiáng)化是高性能齒輪制造的關(guān)鍵工序，通過齒根噴丸強(qiáng)化引入殘余壓應(yīng)力可以有效提高齒輪使用壽命。噴丸與其他機(jī)械加工方法一樣，最終零件表面是凹凸不平的微觀形貌，這些微觀形貌的谷底區(qū)域在工件工作時(shí)會(huì)產(chǎn)生應(yīng)力集中，表面粗糙度越大，產(chǎn)生的應(yīng)力集中越嚴(yán)重，工件的疲勞壽命下降越快。因此，精確計(jì)算齒根噴丸粗糙表面的應(yīng)力與應(yīng)力集中系數(shù)對(duì)評(píng)估齒輪彎曲疲勞壽命極為關(guān)鍵。

　　對(duì)于粗糙表面與應(yīng)力集中系數(shù)的關(guān)聯(lián)性研究，NEUBER將粗糙表面輪廓近似為連續(xù)相鄰的缺口, 并由缺口深度和寬度是影響缺口應(yīng)力集中的主要因素，推廣到基于粗糙表面參數(shù) R_z 和谷底曲率半徑 ρ 的應(yīng)力集中系數(shù)計(jì)算公式。AROLA 等依據(jù) NEUBER 提出的經(jīng)驗(yàn)公式建立了 Arola-Ramulu 模型，將實(shí)際的表面形貌近似為連續(xù)相鄰的正弦狀缺口展開研究，并通過實(shí)驗(yàn)驗(yàn)證了該模型的準(zhǔn)確性。CHAN對(duì)電子束加工得到的板塊表面進(jìn)行測(cè)量，研究了具有凹痕特征的形貌的應(yīng)力集中，提出了加工制造過程中對(duì)表面形貌的要求;KANTZOS借助神經(jīng)網(wǎng)絡(luò)，將應(yīng)力集中和表面形貌進(jìn)行關(guān)聯(lián)分析，實(shí)現(xiàn)了裂紋萌生的預(yù)測(cè)。在數(shù)值計(jì)算方面，學(xué)者們通過對(duì)具有特定規(guī)則的形貌進(jìn)行簡(jiǎn)化研究，采用一階邊界攝動(dòng)法，基于 Hilbert 變換和能量守恒原理，預(yù)測(cè)了拉伸作用下二維表面形貌上的應(yīng)力分布并進(jìn)行實(shí)驗(yàn)驗(yàn)證。然而，目前可見文獻(xiàn)對(duì)于粗糙表面與應(yīng)力集中系數(shù)的關(guān)聯(lián)規(guī)律大多基于二維表面形貌展開研究，對(duì)于三維實(shí)測(cè)表面形貌引起的應(yīng)力集中研究較少，MEREUTA分別將 Neuber 模型和 Arola-Ramulu 模型的公式中二維粗糙表面參數(shù)轉(zhuǎn)換為對(duì)應(yīng)的三維粗糙表面參數(shù)進(jìn)行計(jì)算，并經(jīng)過試驗(yàn)驗(yàn)證了修正公式的可靠性。LI基于表面形貌重構(gòu)技術(shù)對(duì)重構(gòu)表面引起的應(yīng)力集中進(jìn)行分析，利用少數(shù)幾個(gè)重構(gòu)平面模型分析了粗糙表面參數(shù) S_a 疲勞壽命之間的關(guān)系。XU研究了腐蝕對(duì) Q235 鋼板疲勞壽命的影響，借助 ANSYS 軟件生成了三維粗糙形貌并分析了三種腐蝕坑與疲勞壽命的關(guān)系，最后與疲勞實(shí)驗(yàn)進(jìn)行對(duì)比驗(yàn)證。在 GB/T 3480—1997 中通過大量的疲勞實(shí)驗(yàn)提出了相對(duì)齒根表面狀況系數(shù)經(jīng)驗(yàn)公式,然而該經(jīng)驗(yàn)公式?jīng)]有給出粗糙度表面應(yīng)力集中系數(shù)的計(jì)算表達(dá)式，而是通過相對(duì)齒根應(yīng)力集中系數(shù)的比值擬合的公式間接獲得應(yīng)力集中對(duì)齒輪彎曲疲勞壽命的影響，且該經(jīng)驗(yàn)公式僅用 1 個(gè)二維粗糙度參數(shù) R_z 進(jìn)行表征，顯然無法精確評(píng)估三維實(shí)際表面產(chǎn)生的應(yīng)力集中。

　　隨著對(duì)三維表面的認(rèn)識(shí)，用于表征粗糙表面形貌特征的三維粗糙度參數(shù)多達(dá)幾十個(gè)，不同的粗糙度參數(shù)表征表面的不同性能，而對(duì)于表征應(yīng)力集中的三維粗糙度參數(shù)，目前尚未得到共識(shí)。且無論是二維還是三維粗糙表面問題，現(xiàn)階段所有研究對(duì)象均是以平板或者簡(jiǎn)化為平面模型進(jìn)行研究，解決齒根粗糙表面下實(shí)際應(yīng)力精確計(jì)算，得到表征粗糙表面應(yīng)力集中合適的形貌參數(shù)問題是高端齒輪、高性能齒輪設(shè)計(jì)制造的核心技術(shù)之一，具有較大的理論與工程應(yīng)用價(jià)值。

　　論文以磨削后噴丸的直齒輪為研究對(duì)象，提取齒輪齒根三維粗糙表面形貌信息，依據(jù)空間坐標(biāo)變換原理，研究齒根過渡曲面添加三維粗糙面的建模方法，從而突破現(xiàn)有齒輪強(qiáng)度分析只能創(chuàng)建光滑表面模型的限制，建立基于實(shí)測(cè)三維表面形貌的齒根區(qū)域有限元模型，對(duì)有限元計(jì)算得到的數(shù)據(jù)進(jìn)行非線性曲線擬合，得到齒根三維粗糙表面參數(shù)與應(yīng)力集中系數(shù) Kt 的關(guān)聯(lián)公式，為后續(xù)的齒輪彎曲疲勞壽命精準(zhǔn)預(yù)測(cè)和齒輪表面完整性參數(shù)設(shè)計(jì)提供相應(yīng)的參考。

　　一、齒根表面形貌測(cè)量

　　為研究齒根噴丸表面危險(xiǎn)區(qū)域的表面粗糙度與應(yīng)力集中系數(shù)的關(guān)系，首先確定齒根危險(xiǎn)區(qū)域的具體位置，基于磨削噴丸加工的齒輪(表 1)，借助 Matlab 編程生成直齒輪齒廓，將齒廓導(dǎo)入 CATIA，通過陣列、拉伸等操作生成直齒輪實(shí)體模型，將齒輪實(shí)體模型導(dǎo)入 ABAQUS 中進(jìn)行準(zhǔn)靜態(tài)仿真分析，查看仿真結(jié)果應(yīng)力云圖，確定了齒根危險(xiǎn)區(qū)域(應(yīng)力較大區(qū)域)。

　　其次，測(cè)量表面形貌。本文借助 ABAQUS 研究齒根危險(xiǎn)區(qū)域表面粗糙度對(duì)應(yīng)力集中的影響，通過調(diào)整網(wǎng)格節(jié)點(diǎn)坐標(biāo)值實(shí)現(xiàn)粗糙度的添加，而齒根危險(xiǎn)區(qū)域的網(wǎng)格數(shù)量取決于形貌測(cè)量的采樣點(diǎn)數(shù)，因此，選擇較為合適的采樣間隔減小網(wǎng)格數(shù)量、提高仿真效率是很有必要的。文中采用白光干涉儀 Wyko NT9100 對(duì)粗糙形貌進(jìn)行測(cè)量。對(duì)于表面形貌的測(cè)量，采樣間隔越小，表面形貌還原越完整?？紤]文中對(duì)有限元網(wǎng)格數(shù)量的要求，在保留形貌基本特征前提下，適當(dāng)選擇較大的采樣間隔，舍棄部分形貌信息。一般地，采樣頻率大于形貌頻率 4 倍時(shí)，采樣離散數(shù)據(jù)點(diǎn)可以基本復(fù)原形貌的大致輪廓特征。

　　文中以噴丸后的粗糙表面為研究對(duì)象，主要研究噴丸處理后產(chǎn)生的丸坑形貌對(duì)應(yīng)力集中的影響。一般地，表面形貌的應(yīng)力集中最有可能產(chǎn)生在粗糙形貌最深的幾個(gè)谷底區(qū)域之一，已知齒輪噴丸直徑約為 120 μm，噴丸后產(chǎn)生的幾個(gè)最深丸坑深度平均值絕大部分在5 μm以上(詳細(xì)數(shù)據(jù)見表2中的 S_10z )，則丸坑的表面直徑大于 40 μm。因此，在滿足采樣頻率大于 4 倍的形貌頻率(丸坑表面直徑)的基礎(chǔ)上，論文測(cè)量時(shí)選擇采樣間隔為 10 μm，采樣面積為1 mm ×1 mm。得到磨削后噴丸表面形貌如圖 1 所示。

　　二、考慮三維粗糙度的齒根應(yīng)力有限元計(jì)算

　　齒輪網(wǎng)格劃分及細(xì)化

　　齒根應(yīng)力有限元計(jì)算需要精確的齒根形貌與齒面形貌，漸開線齒面的計(jì)算按漸開線方程計(jì)算，齒根過渡曲面的構(gòu)建按課題組前期提出的建模方法構(gòu)建，詳細(xì)齒根過渡曲面構(gòu)建方法可參考文獻(xiàn)。

　　本文通過調(diào)整網(wǎng)格節(jié)點(diǎn)坐標(biāo)實(shí)現(xiàn)粗糙度的添加工作，在齒根區(qū)域的網(wǎng)格節(jié)點(diǎn)數(shù)目應(yīng)與采樣離散點(diǎn)的數(shù)目一致?？紤]直齒輪實(shí)體模型是由齒廓方程生成的，沿齒廓方向各點(diǎn)的受力狀態(tài)基本相同。因此在建立齒輪有限元模型時(shí)將表 1 直齒輪參數(shù)的齒寬設(shè)定為與采樣面積寬度一致的 1 mm，大大降低重復(fù)的網(wǎng)格數(shù)量。此外，為了驗(yàn)證齒寬為 1 mm 時(shí)齒輪有限元模型的可靠性，考慮模型邊緣效應(yīng)等問題，分別對(duì)齒寬為 1 mm、5 mm、10 mm 的有限元模型進(jìn)行仿真，保證其他邊界、材料等條件相同，載荷扭矩分別為T、5T 、10T 進(jìn)行仿真，并對(duì)仿真結(jié)果進(jìn)行對(duì)比分析，仿真結(jié)果表明三種齒寬條件下齒根彎曲應(yīng)力的最大值相差在 5%以內(nèi)，證實(shí) 1 mm 齒寬條件下齒輪模型的可靠性。其次將齒輪切分成兩部分，一部分包含齒根危險(xiǎn)區(qū)域，需要進(jìn)行網(wǎng)格細(xì)化，另一部分不需要調(diào)整網(wǎng)格節(jié)點(diǎn)，因此，網(wǎng)格數(shù)量不做要求。切分操作在 CATIA 中完成，具體切分過程如圖 2 所示。

　　將切分后的齒輪兩部分分別導(dǎo)入 Hypermesh 中進(jìn)行網(wǎng)格劃分。對(duì)齒根危險(xiǎn)區(qū)域網(wǎng)格進(jìn)行局部細(xì)化，齒根細(xì)化區(qū)域網(wǎng)格尺寸為10 μm ×10 μm，細(xì)化層數(shù)為六層。同時(shí)，為保證網(wǎng)格質(zhì)量，在齒輪嚙合區(qū)域網(wǎng)格與齒根危險(xiǎn)區(qū)域網(wǎng)格的尺寸相差不宜過大，未包含齒根危險(xiǎn)區(qū)域部分齒輪對(duì)仿真過程關(guān)注的彎曲應(yīng)力無影響，網(wǎng)格尺寸可盡量調(diào)大，最后劃分得到有限元網(wǎng)格數(shù)約為 25 萬，直齒輪的網(wǎng)格劃分模型如圖 3 所示。

　　齒根危險(xiǎn)區(qū)域表面粗糙形貌的生成

　　將 Hypermesh 劃分好的齒輪模型保存為 inp 文件導(dǎo)入 ABAQUS 中，借助 Python 對(duì) ABAQUS 進(jìn)行二次開發(fā)實(shí)現(xiàn)齒根粗糙形貌的添加，具體操作步驟如下所示。

　　(1) 刪除直齒輪有限元模型中非齒根細(xì)化區(qū)域網(wǎng)格，保留齒根網(wǎng)格細(xì)化區(qū)域。

　　(2) 對(duì)齒根網(wǎng)格細(xì)化區(qū)域表層創(chuàng)建節(jié)點(diǎn)集，創(chuàng)建節(jié)點(diǎn)集后刪除表層網(wǎng)格選擇第二層網(wǎng)格表面創(chuàng)建節(jié)點(diǎn)集，依次創(chuàng)建不同層數(shù)的節(jié)點(diǎn)集，此處論文創(chuàng)建六層節(jié)點(diǎn)集。

　　(3) 利用 Python 對(duì) ABAQUS 進(jìn)行二次開發(fā)，提取步驟(2)中創(chuàng)建的節(jié)點(diǎn)集的節(jié)點(diǎn)坐標(biāo)信息;

　　(4) 由步驟(3)獲得的節(jié)點(diǎn)坐標(biāo)信息是雜亂無序的，對(duì)節(jié)點(diǎn)坐標(biāo)依次按照節(jié)點(diǎn)沿齒廓方向以及節(jié)點(diǎn)高度值大小進(jìn)行排序，獲得規(guī)則的網(wǎng)格節(jié)點(diǎn)編號(hào)以及對(duì)應(yīng)節(jié)點(diǎn)坐標(biāo);

　　(5) 基于空間坐標(biāo)變換原理，將采樣離散點(diǎn)高度信息以垂直于齒根危險(xiǎn)區(qū)域曲面的方式分別對(duì)應(yīng)添加到 ABAQUS 提取的節(jié)點(diǎn)集中，即改變相應(yīng)節(jié)點(diǎn)的坐標(biāo)值。具體變換過程如圖 4 所示。

　　如圖 4 所示，直齒輪由齒廓曲線沿齒寬方向拉伸生成，因此用平面齒廓曲線代替直齒輪模型。其中，P₀ 為齒廓上齒根危險(xiǎn)區(qū)域的一點(diǎn)，該點(diǎn)坐標(biāo)由步驟(4)獲得，P₀ 點(diǎn)對(duì)應(yīng)曲線上的法向量由齒廓方程得到，Δ_Zi 為采樣離散點(diǎn)高度值，P₁ 為 P₀ 沿其對(duì)應(yīng)法向量 n 方向上添加形貌高度值 Δ_Zi 后的坐標(biāo)。其中，第一層網(wǎng)格節(jié)點(diǎn)垂直曲面高度增加值為 Δ_Zi 的 6/6 ，第二層網(wǎng)格節(jié)點(diǎn)垂直曲面高度增加值為 Δ_Zi 的 5/6 ，第三層增加值為 4/6 ，依次類推，直至離散點(diǎn)高度信息全部添加完成。

　　(6) 將修改完成的網(wǎng)格節(jié)點(diǎn)信息替換前面的齒輪模型 inp 文件對(duì)應(yīng)的網(wǎng)格節(jié)點(diǎn)，至此，齒根危險(xiǎn)區(qū)域表面粗糙形貌生成。

　　生成的齒輪形貌示意如圖 5 所示。

　　直齒輪的齒根應(yīng)力仿真計(jì)算

　　為了獲得粗糙表面參數(shù)與應(yīng)力集中系數(shù)的關(guān)聯(lián)規(guī)律，根據(jù) ISO 25178 標(biāo)準(zhǔn)，選擇其中幾個(gè)常用的空間粗糙度表征參數(shù)進(jìn)行計(jì)算分析：算術(shù)平均高度 S_a 、最大谷深 S_v 和十點(diǎn)區(qū)域高度 S_10z 。各參數(shù)的定義如下

　　式中，z 為采集形貌離散點(diǎn)的高度，i和 j 分別表示噴丸形貌的凸峰和凹谷。三維粗糙表面參數(shù) S_a 、S_v 、S_10z 與之對(duì)應(yīng)的二維粗糙表面參數(shù)為 R_a、R_v、R_z，二維粗糙表面參數(shù)是在三維粗糙表面上截取一條線段進(jìn)行計(jì)算，因此，二維粗糙表面參數(shù)會(huì)丟失部分形貌信息，為了更完整地研究粗糙表面產(chǎn)生的應(yīng)力集中，論文選擇三維粗糙表面參數(shù)進(jìn)行形貌表征。

　　完成齒根危險(xiǎn)區(qū)域添加表面粗糙度的工作后，論文選擇三組齒輪副的嚙合模型作為研究對(duì)象，其中小齒輪的中間齒輪為添加粗糙度的齒輪，兩邊齒輪則未添加。建立耦合和接觸對(duì)，將切分的包含粗糙度的中間齒輪用綁定約束進(jìn)行連接，添加載荷并設(shè)置約束條件，采用準(zhǔn)靜態(tài)分析方法進(jìn)行仿真分析，齒輪對(duì)的有限元模型如圖 6 所示。

　　為了研究齒根危險(xiǎn)區(qū)域表面粗糙度與應(yīng)力集中系數(shù)的關(guān)聯(lián)規(guī)律，本文采用的應(yīng)力集中系數(shù)定義如下

　　式中，σ_max 表示帶有缺口模型的缺口處應(yīng)力最大值，σ_n 表示光滑試樣無缺口的名義應(yīng)力。此處取σ_max 為添加粗糙度齒輪模型齒根彎曲應(yīng)力最大值，σ_n 為無粗糙度時(shí)齒根彎曲應(yīng)力最大值。由該定義可知，通過兩者的比值可消除 1 mm 齒寬與實(shí)際齒寬不同帶來的影響。

　　將添加粗糙度的齒輪仿真得到的齒根彎曲應(yīng)力最大值與未添加粗糙度的齒輪彎曲應(yīng)力最大值相比，即得到該粗糙表面參數(shù)下對(duì)應(yīng)的應(yīng)力集中系數(shù) K_t 。有無粗糙度的齒輪彎曲應(yīng)力分布對(duì)比如圖 7 所示。為了更清晰地看出齒根粗糙度表面的應(yīng)力分布，只選取齒輪切分包含齒根危險(xiǎn)區(qū)域部分的進(jìn)行展示。

　　如圖 7 所示，帶有表面粗糙形貌的齒根應(yīng)力分布與無粗糙度時(shí)的齒根應(yīng)力分布相差較大，粗糙表面的齒根應(yīng)力集中主要分布在表面形貌的波谷處，符合實(shí)際情況。取 8 個(gè)形貌樣本進(jìn)行仿真分析，得到各粗糙表面參數(shù)與應(yīng)力集中系數(shù)之間的仿真結(jié)果如表 2 所示。

　　三、粗糙表面參數(shù)與齒根應(yīng)力集中系數(shù)關(guān)聯(lián)規(guī)律

　　不同于現(xiàn)有對(duì)平面粗糙形貌應(yīng)力集中研究的計(jì)算方法，論文提出了一種考慮齒根過渡曲線宏觀幾何結(jié)構(gòu)上的微觀形貌產(chǎn)生的應(yīng)力集中求解方法。為了驗(yàn)證論文所研究的求解方法，將論文的仿真結(jié)果與 GB/T 3480—1997 中的相對(duì)齒根表面狀況系數(shù) Y_RrelT (通過試驗(yàn)獲得)進(jìn)行比較。

　　在 GB/T 3480—1997 中，齒根表面狀況系數(shù)的定義是“考慮齒輪根部的表面狀況，主要是齒根圓角(即齒根危險(xiǎn)區(qū)域)處的表面粗糙度對(duì)齒根彎曲強(qiáng)度的影響”，齒根表面狀況系數(shù)即為齒根粗糙表面產(chǎn)生的應(yīng)力集中系數(shù)。而相對(duì)齒根表面狀況系數(shù)定義是“所計(jì)算齒輪的齒根表面狀況系數(shù)與試驗(yàn)齒輪的齒根表面狀況系數(shù)的比值”，給出的計(jì)算經(jīng)驗(yàn)公式是

　　式中，R_z 為二維粗糙表面參數(shù)中的微觀不平度十點(diǎn)高度，R_z 對(duì)應(yīng)論文中的十點(diǎn)區(qū)域高度 S_10z，基于 S_10z 的相對(duì)齒根表面狀況系數(shù)公式為

　　論文仿真計(jì)算的表面粗糙度產(chǎn)生的應(yīng)力集中系數(shù)(齒根表面狀況系數(shù))，與相對(duì)齒根表面狀況系數(shù)的定義不同，因此，不能直接將粗糙表面參數(shù)代入求解，需要進(jìn)行轉(zhuǎn)換計(jì)算，轉(zhuǎn)換計(jì)算的原理方法如下：將國(guó)標(biāo)中定義的試驗(yàn)齒輪(粗糙表面參數(shù) S_10z =11.27 μm )對(duì)應(yīng)的應(yīng)力集中系數(shù)除以所需要計(jì)算齒輪的應(yīng)力集中系數(shù)即為該計(jì)算齒輪的相對(duì)齒根表面狀況系數(shù)。試驗(yàn)齒輪的粗糙表面參數(shù) S_10z 與樣本 8 基本相同。因此，應(yīng)力集中系數(shù)按樣本 8 數(shù)據(jù)進(jìn)行計(jì)算，經(jīng)轉(zhuǎn)換計(jì)算，表 2 中各粗糙度對(duì)應(yīng)的仿真計(jì)算獲得的相對(duì)齒根表面狀況系數(shù)為樣本 8 的應(yīng)力集中系數(shù)除以對(duì)應(yīng)樣本的應(yīng)力集中系數(shù)。將各個(gè)樣本對(duì)應(yīng)的十點(diǎn)區(qū)域高度 S_10z 代入式(6)中，即得到國(guó)標(biāo)經(jīng)驗(yàn)公式所計(jì)算的Y_RrelT 值。對(duì)比結(jié)果如圖 8 所示。

　　如圖 8 所示，仿真計(jì)算得到的相對(duì)齒根表面狀況系數(shù)與國(guó)標(biāo)經(jīng)驗(yàn)公式計(jì)算得到的系數(shù)整體趨勢(shì)相同，且兩者計(jì)算得到的系數(shù)值也十分接近。此外，在國(guó)標(biāo)中注釋了經(jīng)過強(qiáng)化處理(如噴丸)的齒輪，其相對(duì)齒根表面狀況系數(shù)要稍微大于計(jì)算經(jīng)驗(yàn)公式所確定的數(shù)值，從圖 8 整體趨勢(shì)不難看出仿真值稍微大于經(jīng)驗(yàn)公式計(jì)算的數(shù)值，與國(guó)標(biāo)中的注釋說法一致。圖 8 從側(cè)面論證了論文仿真求解方法的可信度。此外，國(guó)標(biāo)經(jīng)驗(yàn)公式間接表示了齒輪齒根處形貌引起的應(yīng)力集中，而本文則直接求解對(duì)應(yīng)形貌的應(yīng)力集中，表現(xiàn)更加直觀、具體。此外，與實(shí)驗(yàn)獲得應(yīng)力集中時(shí)耗時(shí)耗力相比，通過論文可以快速、低成本計(jì)算出應(yīng)力集中系數(shù)，為后續(xù)精確計(jì)算齒輪彎曲疲勞壽命提供技術(shù)支持。

　　為了得到噴丸表面三維粗糙表面參數(shù)與應(yīng)力集中系數(shù)的關(guān)聯(lián)規(guī)律，借助 Matlab 軟件的 cftool 工具包對(duì)表 2 中的數(shù)據(jù)進(jìn)行關(guān)系式擬合，得到擬合關(guān)系如圖 9 所示。

　　圖 9a、9b、9c 分別表示粗糙表面參數(shù) S_a 、S_v 、S_10z 與應(yīng)力集中系數(shù) Kt 的擬合曲線圖，對(duì)應(yīng)的擬合相關(guān)系數(shù) R 分別為 0.799，0.784，0.914。在對(duì) K_t 的擬合關(guān)系中，十點(diǎn)區(qū)域高度 S_10z 的擬合效果最好，擬合公式如式(7)所示

　　式(7)的擬合公式可以較好表征粗糙表面形貌與應(yīng)力集中系數(shù)的定量關(guān)系。需要說明的是，加工零部件圖紙上標(biāo)注的粗糙表面參數(shù) Ra 是粗糙表面的二維表征參數(shù)，表征三維粗糙表面具有一定的局限性，選擇三維粗糙表面參數(shù)更為合適。從論文的分析結(jié)果看，為更加科學(xué)的表征三維粗糙表面對(duì)應(yīng)力集中的影響，使用十點(diǎn)區(qū)域高度參數(shù) S_10z 更加科學(xué)，工程應(yīng)用中，可優(yōu)先使用 S_10z 對(duì)應(yīng)力集中系數(shù)進(jìn)行計(jì)算。

　　四、結(jié)論

　　(1) 本文提出了一種基于齒根三維粗糙表面形貌的有限元建模與分析方法，解決了三維粗糙齒根表面應(yīng)力集中系數(shù)計(jì)算問題，為齒輪表面完整性與抗疲勞設(shè)計(jì)研究提供了一種基礎(chǔ)方法。

　　(2) 論文采用有限元計(jì)算得到的表面粗糙度引起的應(yīng)力集中系數(shù)與 GB/T 3480—1997 中給出的相對(duì)齒根表面狀況系數(shù)經(jīng)驗(yàn)公式進(jìn)行間接對(duì)比，兩者計(jì)算結(jié)果基本相同。

　　(3) 對(duì)有限元計(jì)算得到的數(shù)據(jù)進(jìn)行擬合分析，結(jié)果表明粗糙表面參數(shù) S_a 、Sv 、S_10z 擬合應(yīng)力集中系數(shù)的擬合公式的相關(guān)系數(shù)分別為 0.799，0.784， 0.914，十點(diǎn)區(qū)域高度參數(shù) S_10z 能較好地表征齒根表面的應(yīng)力集中，工程應(yīng)用中可用十點(diǎn)區(qū)域高度 S_10z 對(duì)噴丸粗糙表面產(chǎn)生的應(yīng)力集中進(jìn)行計(jì)算與表征。

　　參考文獻(xiàn)略.