針對螺旋錐齒輪加工過程中因無法避免誤差而導(dǎo)致齒面加工精度難以保證的問題，提出螺旋錐齒輪加工誤差控制模型及加工參數(shù)修正方法。首先，基于實際加工中的刀具與工件相對位置與相對運動關(guān)系，依據(jù)坐標(biāo)齊次變換與嚙合原理，確定加工參數(shù)與加工曲面之間的函數(shù)關(guān)系，建立螺旋錐齒輪精確齒面模型;然后，計算實際加工齒面與理論齒面的法向距離，從而建立由加工參數(shù)驅(qū)動的齒面幾何誤差控制模型;接著，對加工參數(shù)進行敏感性分析，選取敏感性較高的加工參數(shù)作為誤差補償模型優(yōu)化變量，以提高優(yōu)化效率;最后，將齒面誤差最小化問題轉(zhuǎn)化為最小二乘法問題，基于改進的 L-M 算法進行求解，得到加工參數(shù)補償量，以此對加工參數(shù)進行修正達到減小齒面加工誤差的目的。采用一對由雙重螺旋法磨削加工得到的螺旋錐齒輪副作為應(yīng)用實例，對該方法進行實際加工驗證，結(jié)果表明：加工參數(shù)調(diào)整后，螺旋錐齒輪齒面加工誤差降低了65%以上，實際測量的齒面絕對誤差均不超過0.005mm，能夠滿足工程實際需求，證明該方法能夠有效提升齒面加工精度。該方法可為螺旋錐齒輪乃至其他復(fù)雜曲面零件加工提供一種加工誤差補償思路。

　　螺旋錐齒輪是廣泛應(yīng)用于航空動力傳動及運載工具等高端裝備中的核心基礎(chǔ)零部件，其高精密制造、加工誤差補償一直以來都是熱點難點問題。在實際加工中，機床的幾何誤差、跟隨誤差等因素導(dǎo)致螺旋錐齒輪加工誤差產(chǎn)生，對其接觸性能、使用壽命、傳動精度產(chǎn)生諸多不良影響。Litvin等系統(tǒng)研究了加工參數(shù)變化對于螺旋錐齒輪嚙合特性等性質(zhì)的影響，提出了基于線性回歸算法的齒面優(yōu)化模型，但其求解精度穩(wěn)定性低、求解效率有待提高。為了增強求解穩(wěn)定性，Shih提出了序列二次歸納法，解決了雅克比矩陣的病態(tài)問題，但加工參數(shù)耦合作用導(dǎo)致的數(shù)值 求解不穩(wěn)定問題仍被忽略。Gabiccini及蔣進科等基于齒面ease-off修正理論以準(zhǔn)雙曲面齒輪為研究對象，分別進行了傳動誤差變 化規(guī)律的研究及多目標(biāo)優(yōu)化設(shè)計的研究。李麗霞等從數(shù)學(xué)角度研究了加工參數(shù)與大輪齒面形狀的關(guān)系，但實際加工中通常不修正大輪，這一研究成果具有局限性。精密測量儀器與技術(shù)的廣泛應(yīng)用，使螺旋錐齒輪加工誤差補償有了進一步發(fā)展。唐進元等基于三坐標(biāo)測量儀，研究了銑削加工的螺旋錐齒輪齒面加工精度控制方法。丁撼等提出了考慮齒面微觀形貌的拓撲優(yōu)化方法，其提出了帶 Dog-leg步的置信域算法但其求解效率仍舊有待提高。彭山東等基于6δ原則探討了螺旋錐齒輪形性協(xié)同參數(shù)反調(diào)問題。以上研究雖然取得了一定進展，但對螺旋齒錐齒輪雙重螺旋法加工齒面的高精度反調(diào)修正方法研究工作文獻較少、尚待深入研究。

　　本文通過建立螺旋齒錐齒輪雙重螺旋法加工精確數(shù)學(xué)模型，提出一種參數(shù)驅(qū)動的雙重螺旋加工參數(shù)修正方法，并通過一對工程齒輪實際加工試驗，驗證了這一方法的可行性。本文方法可為螺旋錐齒輪加工參數(shù)反調(diào)提供幫助。

　　一、螺旋錐齒輪雙重螺旋法加工齒面建模

　　螺旋錐齒輪齒面建模是建立加工參數(shù)與齒面映射關(guān)系的過程。雙重螺旋法通過刀盤模擬“假想平頂齒輪”運動，由刀具切削刃在齒坯上形成齒形完成加工，如圖1所示。

　　基于坐標(biāo)變換，建立刀具與工件的坐標(biāo)轉(zhuǎn)換關(guān)系 Mbc。螺旋錐齒輪加工過程可以表示為

　　式中：F(ξ)為刀具在加工過程中形成的一系列曲面;r(u，θ)為刀具數(shù)字化表達，u為高度方向變量，θ 為回轉(zhuǎn)變量;φ為機床運動基本參數(shù)，表示刀具旋轉(zhuǎn)步長;ξ為螺旋錐齒輪雙重螺旋法的加工參數(shù)，ξ=[q，S_R，E_M ，X_D，X_B，γ_m，R_ai，m，j，χ_HL];q 為角向刀位;S_R 為徑向刀位;E_M 為垂直輪位;X_D 為水平輪位;X_B 為床位;γ_m 為安裝角;R_ai為刀轉(zhuǎn)角;m_c 為滾比;d_q 為刀具傾角;χ_HL為螺旋修正系數(shù)。根據(jù)嚙合原理，聯(lián)立方程組，螺旋錐齒輪雙重螺旋齒面參數(shù)化模型可以表示為

　　式中：n 為加工表面上任意一點的法矢;v為刀具與工件的相對速度;· 為標(biāo)量積符號;× 表示矢量積符號。式(2)為螺旋主齒輪齒面隱式表達。將齒面投影范圍設(shè)置為邊界條件，通過合理的齒面均勻化手段，進行齒面網(wǎng)格化處理，逐點求解。

　　二、參數(shù)反調(diào)修正優(yōu)化算法

　　齒面誤差控制模型

　　加工誤差是指刀具實際形成的包絡(luò)面與工件理論輪廓之間的偏差，圖2所示為曲面零件幾何誤差模型?；诰軠y量技術(shù)對實際曲面進行離散網(wǎng)格測量，得到實際加工曲面采樣點 PCMM ，根據(jù)采樣點精確擬合出實際加工曲面?；邶X面參數(shù)化模型，離散化求解得到理論齒面點r₁，r₂，…，r_m，理論齒面點的擬合結(jié)果即為理論曲面，實際加工曲面與目標(biāo)曲面的法向距離h_m 即為加工誤差。

　　建立以機床加工參數(shù)為設(shè)計變量的曲面加工誤差模型

　　實際生產(chǎn)中，通常約束式 (3)中|h_i|≤0.005 mm。式(3)中：i表示第i 個采樣點;r((μ，θ，φ)_i，ξ)、n((μ，θ，φ)_i，ξ)為采樣點對應(yīng)的理論曲面離散點徑矢及法矢

　　誤差補償策略

　　理論上，加工誤差是機床加工參數(shù)微小變化導(dǎo)致誤差的疊加，因此考慮加工誤差的真實齒面可以表示為

　　式中：ξ_e =ξ+Δξ= [q₀ +Δq，S_R0 +ΔS_R，E_M0 + ΔE_M0，X_D0 +ΔX_D，X_B0 +ΔX_B，γ_m0 +Δγ_m，R_ai0 + ΔR_ai，m₀+Δm，j₀+Δj，χ_HL0+Δχ_HL]。通過調(diào)整加工參數(shù)可以實現(xiàn)齒面加工誤差的補償，然而優(yōu)化變量過多會降低誤差補償效率。采用敏感性分析方法定量描述加工參數(shù)對加工誤差的影響，制定合理的優(yōu)化策略非常必要。這一定量描述通常被稱為敏感性系數(shù)，參數(shù)敏感性系數(shù)越大，對模型輸出影響越大。基于 Morris法評估加工參數(shù)敏感性，通過對加工參數(shù)進行多次隨機變化采樣，計算并求其平均值，用來評估其敏感性。本文采用齒面誤差控制模型輸出對輸入變量的一階偏導(dǎo)數(shù)與對應(yīng)法矢的點乘結(jié)果作為敏感系數(shù)

　　式中：i∈[1，n]，n 為曲面離散點數(shù)量。對式(6)進行計算，得到全部曲面點關(guān)于各項機床參數(shù)的敏感系數(shù)組成的雅克比矩陣，稱為曲面采樣點隨各項加工參數(shù)變化的敏感系數(shù)矩陣

　　式中：j∈[1，m]，m 為機床設(shè)置參數(shù)數(shù)量。幾何誤差模型可改寫為

　　為便于排序，求每一個機床參數(shù)的綜合敏感系數(shù)，可得

　　對綜合敏感系數(shù)進行排序，選取敏感性最大的 3～5個的加工參數(shù)作為優(yōu)化變量。由于雙重螺旋法加工過程中，凹凸面具有較強的耦合效應(yīng)，在制定優(yōu)化策略時優(yōu)先考慮工作面(通常為小輪凹面與大輪凸面)的參數(shù)敏感性。

　　參數(shù)反調(diào)優(yōu)化算法

　　式(3)具有強烈的非線性方程性與超越穩(wěn)定性，為了克服求解 min(h)時非線性特性導(dǎo)致雅克比矩陣病態(tài)或奇異引起的求解不穩(wěn)定問題，將式(3) 表達成非線性最小二乘問題

　　考慮到螺旋錐齒輪誤差補償問題的強烈非線性特性，為了提高結(jié)果穩(wěn)定性，以改進的 L-M 方法作為優(yōu)化算法?；舅悸肥窃诿恳淮蔚鷷r求解

　　式中：Δξ_k+1、Δξ_k 分別為第k+1次循環(huán)和第k次循環(huán)求得的反調(diào)量;F_k 為 第k次循環(huán)的目標(biāo)函數(shù)值(即預(yù)測誤差);J_k 為第k次循環(huán)的敏感系數(shù)矩陣;λ_k 是為了克服J_k^TJ_k 所帶來不良矩陣的非負參數(shù)。

　　L-M 方法有線搜索及信賴域搜索兩種搜索思想。本文基于信賴域搜索思想，在每次迭代時，在當(dāng)前迭代點的信賴域內(nèi)找到一個試探解，這一解將用于反饋調(diào)節(jié)信賴域。值得注意的是，L-M 方法通過調(diào)節(jié)迭代參數(shù)λ_k間接調(diào)節(jié)信賴域半徑，選取合理的迭代參數(shù) λ_k有利于其發(fā)揮作用。Yamashita 和 Fukshima提及過可以選取λ_k=‖F(xiàn)_k‖² 的方式，而Fan和 Yuan則選取λ_k=‖F(xiàn)_k‖。固定的信賴域調(diào)節(jié)系數(shù)會給求解帶來問題，當(dāng)加工精度調(diào)整到接近最優(yōu)解集時，λ_k 容易因為過小而失去作用。引入阻尼系數(shù)μ，增強 λ_k 的調(diào)節(jié)作用。當(dāng)加工精度接近目標(biāo)值，借助阻尼系數(shù)放大迭代參數(shù)，使求解更穩(wěn)定。從實際下降量與模型下降量的匹配程度及當(dāng)前精度與目標(biāo)值的距離兩方面來更新阻尼系數(shù)，以完成迭代參數(shù)λ_k 的評價和控制，定義如下

　　式中：r_m 為實際下降量與模型下降量之比;h_max = max(h_i)，表示曲面殘余誤差最大值;hgoal為優(yōu)化目標(biāo)，即預(yù)設(shè)曲面最大誤差。在進行優(yōu)化參數(shù)計算時還考慮以下3個問題：一是確保修正后的加工參數(shù)可用于實際加工，需要在避免傳統(tǒng) L-M 方法導(dǎo)致的過度調(diào)整，采取2個舉措(①識別越界參數(shù)，將其賦值為上一次迭代結(jié)果;② 在優(yōu)化策略中去除該參數(shù));二是在計算過程中，傾向于選擇誤差最大值減小的結(jié)果，所以對于信賴域半徑的調(diào)節(jié)要進行綜合考慮;三是為了平衡控制精度、效率與穩(wěn)定性，需要根據(jù)初始誤差設(shè)置優(yōu)化精度。

　　加工參數(shù)修正誤差補償流程

　　加工參數(shù)修正基本流程如圖3所示。第一步，輸入原始加工參數(shù)。第二步，進行算法初始化，根據(jù)誤差模型進行加工參數(shù)敏感性分析，得到每一個加工參數(shù)的綜合敏感性系數(shù)ε_j，進行排序;接著將敏感性強的s(通常取3～5)個參數(shù)以變量形式計入優(yōu)化方案，保留敏感性較弱的參數(shù);制定出優(yōu)化策略。第三步，利用參數(shù)反調(diào)算法進行精確求解，求解完成后需要對當(dāng)前優(yōu)化結(jié)果進行驗證，需要檢查當(dāng)前加工參數(shù)是否溢出邊界，能否滿足加工要求;在此基礎(chǔ)上，還需要檢查是否存在局部誤差補償能力不足的情況，若存在，則需進行局部敏感性分析，調(diào)整優(yōu)化策略;最后將優(yōu)化后的參數(shù)與保留參數(shù)以矩陣形式統(tǒng)一輸出，稱為補償后的加工參數(shù)。在反復(fù)調(diào)整加工參數(shù)以達到補償加工誤差的過程中，可以選擇更多的加工參數(shù)參與誤差補償提高算法的誤差補償能力。

　　三、結(jié)果與討論

　　算例

　　以工程實際中一對高精度螺旋錐齒輪副作為應(yīng)用實例，表1為該對齒輪的初始加工參數(shù)，包含刀盤參數(shù)和機床加工參數(shù)。通過修正機床加工參數(shù)完成補償加工誤差。

　　圖4為小輪齒面離散點，包含工作面和根尖圓角兩部分，藍色部分為工作齒面區(qū)域，紅色部分為齒根圓角區(qū)域。圖中顯示雙重螺旋法加工一次加工形成的完整齒槽的凹凸兩面。圖 5 為齒輪副三維模型：圖5a、圖5b分別為小輪、大輪三維模型;圖5c為齒輪副裝配模型圖。

　　圖6為初始加工誤差測量結(jié)果，其中 N_H 為沿齒高方向齒面點編號，N_W 為沿齒寬方向齒面點編號。圖6a、6b分別為小輪凸面、凹面誤差測量結(jié)果;圖6c、6d分別為大輪凸面、大輪凹面誤差測量結(jié)果。

　　優(yōu)化結(jié)果及討論

　　對初始加工參數(shù)進行敏感性分析，選擇優(yōu)化變量。由于螺旋錐齒輪齒面形狀較為復(fù)雜，為了讓最終的敏感性分析結(jié)果充分體現(xiàn)出該參數(shù)對全齒面的影響力，選擇齒面上5×9個齒面點求敏感系數(shù)，構(gòu)成敏感系數(shù)矩陣。

　　圖7為部分加工參數(shù)的敏感性分析結(jié)果，雙重螺旋法同時加工兩側(cè)齒面，需要考慮兩側(cè)齒面對于加工參數(shù)的敏感性。由圖5可知，對于小輪，加工參數(shù)對齒面不同點的影響能力有所區(qū)別，凹凸面兩側(cè)的敏感性也不同。徑向刀位增大會導(dǎo)致凹面全齒面誤差增大，凸面誤差減小，按照優(yōu)先修正凹面加工誤差的原則，在修正時選擇徑向刀位增大的方向;角向刀位增大，凹面加工誤差均增大，凸面加工誤差趨勢相反;安裝角增大，凹凸面加工誤差均增大，且其影響能力對凹面加工誤差更強;滾比對齒頂與齒根兩個方向的齒面點齒面誤差敏感性更強。

　　按照式(8)計算綜合敏感性系數(shù)ε_j，繪制散點圖如圖8所示，可以觀察到徑向刀位(S_R)、角向刀位(q)、安裝角(γ_m)、滾比(m)對應(yīng)的敏感系數(shù)偏離零點更遠，具有更大的敏感性，故作為本次優(yōu)化過程中的優(yōu)化變量。經(jīng)過調(diào)整優(yōu)化后的加工參數(shù)如表 2 所示。

　　本文分別從累計誤差的下降量及最大誤差下降量來評估誤差補償效果，預(yù)測殘余誤差的統(tǒng)計結(jié)果如表3所示。由表3可見：第一次反調(diào)后，累計誤差與最大誤差均減小，但最大誤差大于0.005mm 不符合實際需求;兩次反調(diào)后小輪凹、凸面累計誤差分別減少了85.29%和60.59%，最大誤差分別減少了 83.7%和68.75%;大輪凹、凸面累計誤差分別減少 了 82.66% 和 88.02%，最大誤差分別減小了 78.79%和75%，說明該方法有效。圖9為螺旋錐齒輪齒面理論加工誤差擬合圖，可以看到擬合后的殘差曲面均勻連續(xù)，齒面殘余誤差均能控制在 0.005mm 內(nèi)。

　　用補償后的加工參數(shù)進行加工驗證，如圖10所示。基于格里森測量機對加工誤差進行測量，實際齒面誤差測量結(jié)果如圖11所示。圖11為網(wǎng)格劃分示意圖，表4為加工參數(shù)修正后小輪、大輪的實際誤差測量結(jié)果。與優(yōu)化前測量結(jié)果對比，可知小輪凹面、凸面加工誤差最大值分別為3.5μm、4.8μm，降低了81.5%、67.8%;大輪凹面、凸面加工誤差最大值分別為4.3μm、4.7μm，降低了78.1%、73.6%。齒面誤差得到了有效減小，滿足實際工程需求。

　　四、結(jié)論

　　本文提出了針對螺旋錐齒輪雙重螺旋加工的一種加工誤差控制方法，通過加工參數(shù)反調(diào)來補償加工誤差;基于改進的 L-M 算法建立加工參數(shù)優(yōu)化算法，克服了傳統(tǒng)方法中雅克比矩陣病態(tài)或奇異的問題;根據(jù)加工誤差對不同加工參數(shù)的靈敏度，選擇合適的加工參數(shù)進行優(yōu)化，有效提高了優(yōu)化效率。從優(yōu)化結(jié)果分析可以看出：本文方法能夠高效地實現(xiàn)螺旋錐齒輪磨削加工精度的控制與優(yōu)化，能在較少的反調(diào)次數(shù)內(nèi)，有效提高實際加工齒面與理論設(shè)計齒面的貼合程度;通過調(diào)整加工參數(shù)對加工誤差進行補償后，最大殘余誤差不超過5μm，得到滿足實際磨齒加工需求的加工參數(shù)。本文工作為雙重螺旋法螺旋錐齒輪齒面高精度加工提供了一種誤差補償新思路，在理論方法與工程應(yīng)用上均有一定參考價值。

　　參考文獻略.