張國慶，李玉光*，陳殿華

（大連大學 機械工程學院，遼寧 大連 116622）

WN齒輪副耦合振動及動態(tài)應力仿真的研究

張國慶，李玉光*，陳殿華

（大連大學 機械工程學院，遼寧 大連 116622）

基于WN弧齒輪傳動特征建立了多自由度切-徑-軸及系統(tǒng)誤差耦合振動動力學模分析型。利用PROE建立齒輪三維實體模型，通過PROE軟件和Ansys軟件數據接口將WN弧齒輪三維模型導入Ansys軟件中進行了3維沖擊動力接觸模擬分析，利用施加載荷進行了WN弧齒輪的應力分析，仿真出WN齒輪齒面接觸的等效應力云圖。在ANSYS/LS-DYNA中構建了齒輪鋼柔性體動態(tài)接觸模型，以用于齒輪轉動的動態(tài)分析，結果表明WN弧齒輪的切向與軸向方向的振動明顯大于徑向的振動。

WN齒輪；動態(tài)接觸；應力仿真；耦合振動

0 引言

WN齒輪也名為雙圓弧齒輪，以其良好的跑合性能、高承載能力、有效壽命長等優(yōu)點和獨特的多點嚙合方式和凸、凹齒接觸，其接觸強度比漸開線齒輪在相同情況下更大。因此在石油化工、汽輪機、航空直升機等領域得到廣泛應用。然而 WN齒輪也具有嚙合狀態(tài)復雜、載荷反復交變等特點[1]。為了改善齒輪傳動的靜動態(tài)特性和運轉性能，對 WN齒輪的嚙合性能、接觸狀態(tài)進行仿真分析是迫切需要的，這樣也能更好確定齒面壓力分布、綜合嚙合剛度、齒間載荷分配等。

1 WN弧齒輪傳動系統(tǒng)的動力學模型

1.1 WN弧齒輪齒面的基本方程

用GB/12759－1991型號的WN齒輪作為研究對象，WN齒輪的齒廓是由凸弧、凹弧、凸凹弧的過渡圓弧、齒根圓4段圓弧構成。每個圓弧都由各自的半徑、圓心位置確定。WN齒輪齒廓的公式方程為：

式(1)中：r為齒輪節(jié)圓半徑，ρ為相對曲率半徑，β為螺旋角，e為凹齒接觸點處槽寬，N為齒廓圓心相對所選坐標軸的偏移量，φ為齒輪轉角。

1.2 傳動系統(tǒng)動力學方程

[2、3]，考慮切-徑-軸方向及系統(tǒng)誤差的耦合作用，利用PROE建立齒輪三維實體模型建立動力學模型。如圖1所示，為x軸、y軸、z軸方向的平移振動位移；為主從動齒輪繞 z軸方向的扭轉振動位移；為主從動齒輪沿y軸、x軸方向的擺動振動位移。WN齒輪傳動系統(tǒng)的廣義位移列陣可表示為列向量：

圖1 WN齒輪副嚙合的動力學模型

由此列出 WN齒輪傳動系統(tǒng)各振動方向的運動微分方程，得到該傳動系統(tǒng)的振動方程為：

以上齒輪傳動系統(tǒng)中的阻尼系數和剛度均由試驗方法獲得。用矩陣表示上述運動方程為：

2 模型求解方法及數值求解

2.1 模型求解方法

用直接積分法（也稱時間變量離散法）求解（4）中的動力學方程。即把時間區(qū)間劃分成n個時間均勻步長 tΔ。首先用t0時刻的推算出時刻的位移值、速度值和加速度值，然后再逐步積分求解出每個瞬時時間的，本文采用Newmark-β法，它是線加速度作為整個區(qū)間的實際加速度，即

則（7）變?yōu)?/p>

(8)

2.2 模型數值求解

由2.1的解析可知，齒輪耦合振動系統(tǒng)的動態(tài)響應求解的程序用Newmark-β法和數值方法進行編寫[4]，利用這個程序，對齒輪接觸的動力學運動微分方程的進行計算，得到 WN齒輪耦合振動系統(tǒng)中各自由度的時間響應歷程，并且計算出其振動加速度、速度和各坐標軸方向的動態(tài)力的相應值。如圖2所示，其計算程序框圖為。

圖2 齒輪動態(tài)響應的計算程序框圖

3 WN齒輪傳動系統(tǒng)的動態(tài)仿真

3.1 WN齒輪建模

由于 WN齒輪齒面復雜，使用齒面齒廓方程式構建完整的齒輪齒面較麻煩。故通過PROE軟件 和Ansys數據接口將WN弧齒輪三維模型導入Ansys軟件中，并對模型劃分網格，施加載荷，在Ansys軟件中進行WN弧齒輪的應力分析[5]。分析模型如圖3：

圖3 WN弧齒輪導入的分析模型

3.2 齒對動態(tài)接觸模型

單元選擇：針對剛柔體模型，選用SOLD164單元和SHELL163單元。SOLD164單元是三維實體單元，為8節(jié)點六面體單元，各節(jié)點擁有9個自由度，將用于齒輪輪體。而此單元沒有旋轉自由度，不能施加轉速和轉矩。SHELL163單元為二維實體單元，彌補了這一缺陷，用于齒輪內圈表面的單元劃分。將SHELL163定義為剛形體約束其旋轉自由度，這樣有利于齒輪扭矩和轉速的施加。

材料屬性和實常數的設置：對于 SHELL163單元，設置為S/R Co-rotional Hughes-liu，即所謂的面內多積分改進型，以避免沙漏變形在分析工程中出現(xiàn)。而對于SHELL163單元，其均勻厚度設置為0.1。對于齒輪材料屬性，其密度設置為7.81-6 kg/mm3，彈性模量設為2.05e8Mpa，泊松比為0.25。在材料的設置中還有必要對旋轉和平移設置約束。兩齒輪的參數均設置為：translational constraint parameter-7(設置x,y，z方向的位移約束)，rotational constraint parameter-4 (設置x和y方向的旋轉約束)。

劃分網格和創(chuàng)建PART：分別對大小輪體采用彈性材料1、3，而其內圈內圈則采用剛性材料2、4，分別進行網格劃分(如圖4)，然后創(chuàng)建PART 1、2、3、4；其中2、4是剛體。在前面定義材料是已經將剛體定義為除了繞 Z軸方向旋轉，其他所有約束進行了約束。

定義接觸類型：定義柔體PART1和PART3為接觸對，選擇自動單面接觸，即物體自身或者物體之間在任意時間內的接觸情況將由程序自動判斷。

加載：給剛體PART2加載恒定的剛體角速度，剛體 PART4加載恒定的剛體轉矩[6]，這里假設輸入軸的轉速和輸出軸的轉矩都是恒定不變。

圖4 WN弧齒輪的網格劃分模型

3.3 多齒對齒輪的動態(tài)仿真

利用ANSYS的非線性動力分析模塊進行WN齒輪的動力學接觸仿真分析[5.6]，得到了齒輪輪齒在一個嚙合過程中等效等值應力的變化情況，仿真的動畫在0.001 s內共有22幀，這里僅列出4幅連續(xù)的應力應變云圖（如圖5所示)。同時利用Hertz應力公式計算的相應理論值，將ANSYS仿真出的輪齒等效等值應力與相應的理論值相比較，發(fā)現(xiàn)其偏差不大，而且趨勢相同。

在LS-POST中選取齒根單元4101和單元70397，一個節(jié)點11364（如圖6），分別仿真單元的等效應力隨時間的變化曲線，齒根等效應力隨嚙合時間的增加，先是增加然后又減小，與事實嚙合情況相符。

圖5 嚙合齒輪的動態(tài)等效等值應力圖

圖6 嚙合中節(jié)點4101和70397的應力變化

圖7所示，在主動齒輪加入轉速和轉矩，被動齒輪則加入相應的轉速和轉矩，得到的等值應力分布圖，如圖我們可以看出，在齒根凹弧處，其齒根應力最大，故而會產生斷齒等齒輪失效情況，這與實際情況相符。

如圖8所示，節(jié)點11364在X、Y、Z（X為齒輪切向方向，Y為齒輪徑向方向，Z為齒輪軸向方向）方向的加速度隨時間的變化曲線，再看節(jié)點11364分別在切向、徑向、軸向方向的加速度隨時間增加的變化情況，其中切向與軸向方向的加速度明顯大于徑向方向的加速；從加速度的情況我們可以得出在切向與軸向方向的振動明顯大于徑向的振動。

圖7 嚙合中耦合動載瞬時作用的應力分布

圖8 嚙合周期內各方向的加速度變化

4 結論

(1) 采用切-徑-軸 3方向及系統(tǒng)誤差耦合建立WN弧齒輪動力學模分析型。與工程實際更加接近。通過PROE軟件和Ansys軟件數據接口將WN弧齒輪三維模型導入Ansys軟件中, 建模方便。

(2) 基于ANSYS/LS-DYNA顯式動力程序，建立了齒輪動態(tài)剛柔性體計算模型，真實地模擬齒輪嚙合情況，得到了直觀清楚的仿真結果。

利用 history模塊，成功的仿真出單元等效應力曲線與節(jié)點的加速度情況。得到了圓弧齒輪嚙合過程中的切向、軸向、徑向方向振動情況，為齒輪參數設計和減少齒輪傳動的噪聲提供了理論依據。

參考文獻：

[1] 陳殿華, 李成超. 基于負載與嚙合特征的WN齒輪傳動接觸強度分析[J]. 工程設計學報, 2012, 19(4): 294-297.

[2] 王立華. 汽車螺旋錐齒輪傳動耦合非線性振動研究[D].重慶: 重慶大學, 2003.

[3] Dianhua Chen, Guoqing Zhang, Xi Chen. Research on nonlinear vibration and bend strength character of WN gears pairs in meshing [J]. Advanced Materials Research, 2013, 706-708: 1537-1540.

[4] 陳殿華, 李玉光. 基于彈流潤滑與動力學耦合的WN齒輪傳動的強度設計[J]. 工程設計學報, 2009, 16(5): 364-368.

[5] 屈文濤等. 雙圓弧齒輪的有限元綜合分析方法[J]. 機械強度, 2006, 28(3): 415-418.

[6] 李成超, 陳殿華, 李玉光, 等. WN齒輪傳動的三維建模與嚙合仿真析[J]. 大連大學學報, 2012, 33(6): 5-8.

The Study of Coupling Vibration and Dynamic Stress Simulation for WN Gear Pair

ZHANG Guo-qing, LI Yu-guang*, CHEN Dian-hua
（College of Mechanical Engineering, Dalian University, Dalian 16622, China）

Based on double arc gear characteristics and considering system error, a multi-degree-freedom (Tangential -Radial-Axial) coupled vibration dynamics analysis model was established. In this article, we established the three-dimensional solid model by using PROE, imported the arc gear three-dimensional model into Ansys software through the data interface of Ansys software and PROE software and conducted a three-dimensional simulation anasys of the impact of dynamic contact. By applying load, it conducted a stress analysis for double arc gear and gain the double arc gear's the effective stress clouds. Meanwhile, the rigid-flexible body of gear dynamic contact model based on ANSYS / LS-DYNA was established and the dynamic simulation anasys for WN gear was analyzed. The results demonstrated that the tangential and axial vibration of double arc gear is significantly greater than the radial vibration.

WN gear; dynamic contact; stress simulation; coupling vibration

TH132.41TP311.1

：A

：1008-2395(2013)06-0001-05

2013-08-30

國家自然科學基金項目（51075047）。

張國慶（1988-)，男，碩士研究生，研究方向：摩擦磨損及數值仿真技術；陳殿華（1953-），男，教授，研究方向：機械設計理論及CAD/CAE。

李玉光（1963-），男，教授，研究方向：機械設計。