基于三次NURBS 曲面插值的單點線嚙合齒輪副建模

摘要：【目的】點線嚙合齒輪齒廓方程是復(fù)雜的分段函數(shù)，運用齒輪齒廓方程建模誤差較大。為了更準(zhǔn)確地研究齒輪副嚙合特性和傳動參數(shù)，建立了基于三次非均勻有理B 樣條（Non-Uniform RationalB-Spline， NURBS）曲面插值的單點線嚙合齒輪副模型?！痉椒ā扛鶕?jù)單點線嚙合齒輪副嚙合原理以及坐標(biāo)系變換推導(dǎo)了齒面方程，利用齒面方程獲得離散的齒面數(shù)據(jù)，運用三次NURBS 求解離散數(shù)據(jù)插值而成的曲面控制點矩陣，利用UG二次開發(fā)和曲面的控制點矩陣對齒輪進行齒面重構(gòu)，建立了齒輪副的有限元分析模型，獲得齒輪副的瞬時接觸線、齒面最大接觸應(yīng)力以及最大彎曲應(yīng)力的變化曲線，驗證了三維模型的準(zhǔn)確性?！窘Y(jié)果】結(jié)果表明，單點線嚙合齒輪副主動輪在進入和退出嚙合時都未出現(xiàn)邊緣接觸，且在整個周期內(nèi)最大接觸應(yīng)力變化比較平緩，最大彎曲應(yīng)力出現(xiàn)在從動輪和主動輪的齒根過渡曲面附近。

關(guān)鍵詞：單點線嚙合齒輪副；齒面重構(gòu)；嚙合特性；接觸應(yīng)力

中圖分類號：TG712 DOI：10. 16578/j. issn. 1004. 2539. 2025. 02. 007

0 引言

點線嚙合齒輪是結(jié)合漸開線齒輪及圓弧齒輪特點的一種新型齒輪，具有制造工藝簡便、可分性好、凸齒廓貼合度高等特點，還具備較強的承載能力、低噪聲和平穩(wěn)傳動等優(yōu)勢[1]。學(xué)者們已利用刀具的范成加工和空間嚙合原理得出了點線嚙合齒輪的齒廓方程。由于單點線嚙合齒輪副齒廓方程是復(fù)雜的分段函數(shù)，為了減小建模誤差并降低齒輪副接觸應(yīng)力，必須結(jié)合齒面數(shù)學(xué)模型對嚙合特性進行有限元分析。

近年來，學(xué)者們對單點線嚙合齒輪副嚙合特性進行了大量研究。厲海祥等[2]40-43提出了點線嚙合齒輪的基本概念，指出點線嚙合齒輪的主動輪是一個漸開線變位短齒輪；從動輪齒廓包含兩部分，上部為漸開線凸齒廓，下部為過渡曲線凹齒廓。楊帆等[3]對點線嚙合齒輪的三維建模、變位系數(shù)進行了研究，并給出了參數(shù)選擇封閉圖[4]。黃海等[5-6]對點線嚙合齒輪的齒根彎曲應(yīng)力計算方法、彎曲疲勞強度、變位系數(shù)、螺旋角進行了研究，對從動輪齒根彎曲應(yīng)力進行了修正。單點線嚙合齒輪副的齒形較為復(fù)雜，直接運用其齒廓表達式建模容易造成較大的模型誤差。本文根據(jù)空間嚙合原理求出齒輪的齒面方程，將齒面離散化的精確數(shù)據(jù)點，結(jié)合三次非均勻有理B樣條（Non-Uniform Rational B-Spline， NURBS）插值的齒面重構(gòu)技術(shù)，實現(xiàn)點線嚙合齒輪的三維建模；建立了點線嚙合齒輪7齒對的有限元分析模型，分析了外加載荷條件下單點線嚙合齒輪副接觸應(yīng)力的周期性變化規(guī)律，獲得了齒輪副的瞬時接觸線，驗證了其數(shù)學(xué)模型的準(zhǔn)確性。

1 齒面方程

在已有研究基礎(chǔ)上，根據(jù)空間展成法推導(dǎo)出單點線嚙合斜齒輪副[7]的齒面方程，圖1所示為點線嚙合齒輪加工刀具的齒廓[2]42-43，點劃線S - S'為刀具中線，Q1、Q2 均為刀具圓弧段的圓心，ρ 為圓弧段半徑，坐標(biāo)軸Y1 和刀具節(jié)線重合，刀具壓力角為αn，刀具節(jié)線和中線距離xm = xn Mn。其中，xn 為刀具變位系數(shù)；Mn 為刀具法面模數(shù)。

圖2為在刀具法平面以及端平面方向建立的坐標(biāo)系關(guān)系示意圖。其中，S1 為法平面坐標(biāo)系，S2 為端平面坐標(biāo)系，坐標(biāo)軸Y1 和Y2 之間的夾角為β（刀具螺旋角）。通過齒輪刀具的齒面方程和刀齒形成的空間嚙合原理，得到齒輪漸開線部分齒面方程，即

由于三次NURBS曲線的首末控制點P1 和Pn + 2 為首末型值點Q1 和Qn，即P0 = Q0，Pn + 2 = Qn[12]。根據(jù)插值要求[13]有

式（11）中有n + 1個方程，但方程組中含有n +3個未知數(shù)，若要求方程組有唯一解，需要補充2個邊界條件方程。采用切矢邊界條件獲取的兩個附加方程[9]386為

式中，p*0、p*n 分別為型值點首末端點的切矢；p0、p1表示第1、2個控制點；pn + 2、pn + 1 表示倒數(shù)第1、2個控制點。p*0 、p*n 可以利用型值點Q0 以及與之相鄰的兩點Q1、Q2， 型值點Qn 以及與之相鄰的Qn - 1、Qn - 2 進行拋物線插值， 即p*0 = Y '1 （0）， p*n = Y '2 （1），Y1、Y2 均為拋物線插值表達式。

將齒面方程離散化的數(shù)據(jù)點代入到方程組中，可以求得三次NURBS曲面v 方向的控制點。圖4所示為齒面上v 方向的多條擬合曲線。將求得的v 方向的控制點作為新的型值點，可以得到u 方向的擬合曲線和u 方向的控制點，此時得到的控制點即為整個齒面的總控制點。

得到插值曲面的總控制點之后，利用Matlab工具箱進行曲面插值，插值結(jié)果如圖5所示。

3 三維建模

進行單點線嚙合齒輪副齒面重構(gòu)時，需要選取合適的參數(shù)進行三維建模。齒輪副的具體參數(shù)如表1所示。

將求得的總控制點矩陣以數(shù)據(jù)點的形式導(dǎo)入UG三維建模軟件中，利用UG二次開發(fā)中的grip語言進行曲面插值，齒面插值情況如圖6所示。其中，綠色的點表示Matlab中求出的總控制矩陣數(shù)據(jù)點，灰色曲面是依據(jù)總控制點矩陣插值得到的三次NURBS插值曲面。最終得到的點線嚙合齒輪主動輪及從動輪的三維模型裝配圖如圖7所示。

4 嚙合特性分析

將創(chuàng)建的三維模型導(dǎo)入到Hypermesh 軟件進行網(wǎng)格劃分，基于齒面接觸斑點連續(xù)性的網(wǎng)格尺寸要求[15]，在單齒面的過渡曲面部分劃分16層單元，齒高方向劃分32 層單元， 齒寬方向劃分120 層單元[16]33。綜合考慮有限元模型的計算時間以及收斂性和準(zhǔn)確性，在本文中選擇點線嚙合齒輪7齒對進行有限元模型的齒輪副靜力學(xué)接觸分析。劃分好的單齒對網(wǎng)格模型如圖8（a）所示，7齒對劃分裝配模型如圖8（b）所示。

本文單點線嚙合齒輪副的主、從動輪采用的材料均為45號鋼，材料的彈性模量E=206 GPa，泊松比ν=0. 3，密度ρ=7. 85×10-9 t/mm3。為了保證在齒輪靜力學(xué)分析中齒輪副傳動的穩(wěn)定性，將整個仿真過程分為以下3個分析步[16]34：①對主動輪施加微小轉(zhuǎn)動，使主、從動輪齒面相互接觸，減小齒輪副之間接觸時的沖擊；②釋放從動輪繞中心軸的旋轉(zhuǎn)自由度，并且添加一個緩慢增大的負(fù)載轉(zhuǎn)矩，負(fù)載轉(zhuǎn)矩為20 000 N·m；③賦予主動輪初始轉(zhuǎn)速，使主動輪驅(qū)動從動輪轉(zhuǎn)動。通過齒輪副有限元分析，得到了動態(tài)分析模型接觸應(yīng)力和最大彎曲應(yīng)力的云圖，如圖9所示。圖9（a）所示為齒輪副的瞬時接觸線，在a點處為點接觸，但在負(fù)載作用下點接觸受壓變形成了局部的面接觸，其他位置以線接觸的形式存在，符合點線嚙合齒輪的接觸形式，驗證了三維模型的準(zhǔn)確性。圖9（b）和圖9（c）分別為從動輪和主動輪的最大彎曲應(yīng)力云圖，可以看出，主、從動輪的最大彎曲應(yīng)力都出現(xiàn)在各自齒面的過渡曲面部分。

記圖9（a）中正在退出嚙合的齒為第一齒，接觸線最長的為第二齒，正在進入嚙合的齒為第三齒，還未進入嚙合的為第四齒。為了能更直觀地顯示點線嚙合齒輪的嚙合特性，提取第三齒從嚙入到嚙出的1個嚙合周期內(nèi)各個嚙合位置的最大接觸應(yīng)力，并繪制了與主動輪轉(zhuǎn)角之間的變化曲線，具體結(jié)果如圖10所示。

在1個嚙合周期內(nèi)，單點線嚙合齒輪副分別經(jīng)歷了三齒對嚙合、雙齒對嚙合、三齒對嚙合、雙齒對嚙合、三齒對嚙合等5個階段，直至該單齒脫離接觸，其中在點a（a'）第三齒開始進入嚙合，在點f （ f '）第三齒恰好退出嚙合，單齒的嚙合周期結(jié)束。選擇合適的載荷綜合系數(shù)K，比較齒輪副在單齒對嚙合處齒面最大接觸應(yīng)力的仿真值與理論值，其仿真值為1 327 MPa，理論值為1 274 MPa，仿真和理論值之間的誤差不超過5%，表明仿真結(jié)果準(zhǔn)確。此外，在兩種不同負(fù)載載荷的作用下，其齒面最大接觸應(yīng)力曲線的形狀大致相同，并且嚙合周期范圍一致，這是因為不同的外在載荷不會改變齒輪副的嚙合線長度以及變化趨勢。此外，齒輪副齒面最大接觸應(yīng)力的變化比較平緩，表明齒輪副之間傳動比較平穩(wěn)。

齒輪副的最大彎曲應(yīng)力變化曲線如圖11（a）所示：在mn （m'n'）三齒嚙合區(qū)，由于第一齒的嚙出以及第三齒的嚙入，造成第一、三齒接觸線長度先變長后變短，導(dǎo)致最大彎曲應(yīng)力先變小后變大；在np（n'p'）雙齒嚙合區(qū)，由于第二齒的嚙出以及第三齒的嚙入，造成整個齒輪副接觸線長度先變短后變長，導(dǎo)致最大彎曲應(yīng)力先變大后變??；p（ p'）點之后的曲線是之前mp（m'p'）曲線的又一次周期性變化，其最大彎曲應(yīng)力的變化原理大致相同，直至到達s（s'）點；而在點s（s'）之后則會進入下一個嚙合周期，其曲線的變化趨勢和第一個嚙合周期內(nèi)的變化趨勢大致相同。圖11（b）所示為圖9中A、B 兩點處最大彎曲應(yīng)力隨時間變化的曲線，其中，A、B 兩點的位置皆位于上一齒開始退出嚙合但仍未退出嚙合的位置點上。由圖11（b）可以看出，點A 和點B 在整個嚙合周期內(nèi)只存在1個峰值，且嚙合位置幾乎出現(xiàn)在同一時刻。

5 結(jié)論

利用嚙合原理得到了點線嚙合齒輪的齒面方程；并通過離散點云數(shù)據(jù)，運用三次NURBS曲線對得到的點云數(shù)據(jù)進行插值，獲得了點線嚙合齒輪的精確模型；建立了7齒有限元動態(tài)分析模型齒輪對，研究了點線嚙合齒輪的嚙合特性，得出結(jié)論如下：

1） 齒輪副嚙合的同時存在點嚙合以及線嚙合，點嚙合在負(fù)載的作用下變成局部的面嚙合，面嚙合的位置出現(xiàn)在從動輪漸開線以及過渡曲面交點處，驗證了創(chuàng)建的三維模型的準(zhǔn)確性。

2） 齒輪副相互嚙合時，主、從動輪之間并未產(chǎn)生邊緣干涉，齒輪副之間的最大接觸應(yīng)力曲線沒有出現(xiàn)應(yīng)力突變，并且在整個嚙合周期內(nèi)曲線變化都比較平緩，齒輪傳動較為平穩(wěn)。

3） 齒輪副的齒根最大彎曲應(yīng)力的極大值點都出現(xiàn)在開始退出嚙合的單齒逐漸退出嚙合但仍未退出嚙合的齒輪副齒根位置上，且主、從動輪的彎曲應(yīng)力極大值點幾乎出現(xiàn)在同一時刻。

參考文獻

［1］ 湯魚，孫僑，張帆，等. 點線嚙合齒輪研究現(xiàn)狀與展望［J］. 機械研究與應(yīng)用，2022，35（6）：226-230.

TANG Yu，SUN Qiao，ZHANG Fan，et al. Research status andprospect of point-line meshing gear［J］. Mechanical Research ＆Application，2022，35（6）：226-230.

［2］ 厲海祥，羅齊漢，黃海，等. 點線嚙合圓柱齒輪傳動［J］. 水利電力機械，2007，29（6）：40-43.

LI Haixiang，LUO Qihan，HUANG Hai，et al. Transmission of thepoint amp; line meshing cllindrical gear［J］. Water Conservancy ＆Electric Power Machinery，2007，29（6）：40-43.

［3］ 楊帆，羅齊漢，黃海，等. 基于SolidWorks 的點線嚙合齒輪三維模型的研究［J］. 機械傳動，2017，41（3）：44-49.

YANG Fan，LUO Qihan，HUANG Hai，et al. Research of the threedimensionalmodel of point-line meshing gear based on SolidWorks［J］. Journal of Mechanical Transmission，2017，41（3）：44-49.

［4］ 羅齊漢，厲海祥，張予川，等. 點線嚙合齒輪參數(shù)選擇的封閉圖［J］. 機械工程學(xué)報，2005，41（1）：37-40，45.

LUO Qihan，LI Haixiang，ZHANG Yuchuan，et al. Closed pictureof point-line gear parameter selection［J］. Chinese Journal of MechanicalEngineering，2005，41（1）：37-40，45.

［5］ 黃海，厲海祥，羅齊漢，等. 點線嚙合齒輪齒根彎曲應(yīng)力研究［J］.機械傳動，2011，35（1）：8-11.

HUANG Hai，LI Haixiang，LUO Qihan，et al. Research on toothroot bending stress of point-line meshing gear［J］. Journal of MechanicalTransmission，2011，35（1）：8-11.

［6］ 黃海，厲海祥，羅齊漢. 關(guān)于點線嚙合齒輪螺旋角選擇的研究［J］. 機械傳動，2005，29（3）：9-11.

HUANG Hai，LI Haixiang，LUO Qihan. Research on the choice ofhelix angle of the points-line meshing gear［J］. Journal of MechanicalTransmission，2005，29（3）：9-11.

［7］ 張志強，夏尊鳳，許焰. 任意轉(zhuǎn)角位置的漸開線齒輪齒面參數(shù)方程的研究［J］. 機械設(shè)計與制造，2005（8）：52-54.

ZHANG Zhiqiang，XIA Zunfeng，XU Yan. Research of parameterequation of involute gear teeth surfaces on any corner position［J］.Machinery Design amp; Manufacture，2005（8）：52-54.

［8］ 潘松欽，王小平，宋明軒. NURBS 曲面建模在產(chǎn)品創(chuàng)新設(shè)計中的應(yīng)用［J］. 科學(xué)技術(shù)與工程，2007，7（6）：1052-1056.

PAN Songqin，WANG Xiaoping，SONG Mingxuan. Application ofNURBS surface modeling in product creative concept design［J］.Science Technology and Engineering，2007，7（6）：1052-1056.

［9］ 葉麗，謝明紅. 采用積累弦長法擬合3 次NURBS 曲線［J］. 華僑大學(xué)學(xué)報（自然科學(xué)版），2010，31（4）：383-387.

YE Li，XIE Minghong. Third-order non-uniform rational B-splinecurve fitting based on method of accumulating chord length［J］.Journal of Huaqiao University （Natural Science），2010，31（4）：383-387.

［10］ 呂丹，童創(chuàng)明，鄧發(fā)升，等. 三次NURBS 曲線控制點的計算［J］.彈箭與制導(dǎo)學(xué)報，2006，26（4）：357-359.

LV Dan，TONG Chuangming，DENG Fasheng，et al. Calculatingcontrol points of cubic NURBS curve［J］. Journal of Projectiles，Rockets，Missiles and Guidance，2006，26（4）：357-359.

［11］ 陳慶雪，韓慶瑤. 基于MATLAB的NURBS 曲面造型研究［J］. 電力科學(xué)與工程，2009，25（11）：70-72.

CHEN Qingxue，HAN Qingyao. Research of NURBS-face modelingbased on MATLAB［J］. Electric Power Science and Engineering，2009，25（11）：70-72.

［12］ 蘇曉紅，李東，王宇穎. 基于曲率參數(shù)的NURBS 曲線插值［J］. 哈爾濱工業(yè)大學(xué)學(xué)報，2001，33（1）：108-111.

SU Xiaohong，LI Dong，WANG Yuying. A new method of NURBSinterpolation［J］. Journal of Harbin Institute of Technology，2001，33（1）：108-111.

［13］ 韓慶瑤，董云風(fēng)，師向紅. 三次NURBS 空間曲線的插值研究［J］.煤礦機械，2007，28（1）：44-46.

HAN Qingyao，DONG Yunfeng，SHI Xianghong. Research of interpolationon cubic NURBS space curve［J］. Coal Mine Machinery，2007，28（1）：44-46.

［14］ 陳燕麗，宋緒丁. 基于B樣條遞推算法的NURBS 雙四次曲面矩陣表示［J］. 制造業(yè)自動化，2014，36（4）：90-93.

CHEN Yanli，SONG Xuding. Matrix representation for quarticNURBS surface based on the recursive algorithm of B-spline［J］.Manufacturing Automation，2014，36（4）：90-93.

［15］ 劉耀東，陳凱，華曉波，等. 考慮彈流潤滑作用的斜齒輪嚙合分析方法研究［J］. 計算力學(xué)學(xué)報，2015，32（6）：838-844.

LIU Yaodong，CHEN Kai，HUA Xiaobo，et al. Research on meshanalysis method of helical gear with EHL［J］. Chinese Journal ofComputational Mechanics，2015，32（6）：838-844.

［16］ 冉小平，谷麗瑤，王健帆. 等變位斜齒錐齒輪動態(tài)嚙合特性有限元分析［J］. 機械傳動，2022，46（6）：31-37.

RAN Xiaoping，GU Liyao，WANG Jianfan. Finite element analysisof dynamic meshing characteristic of skew bevel gear with equivocalmodification［J］. Journal of Mechanical Transmission，2022，46（6）：31-37.

基金項目：中央高校基本科研業(yè)務(wù)費專項資金（2682020CX33）