·，·

(1.新疆大學(xué)數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 新疆 烏魯木齊 830046；2.新疆大學(xué)機(jī)械工程學(xué)院, 新疆 烏魯木齊 830047)

一對嚙合齒輪模型離散成SPH粒子方法研究

熱合買提江·依明1，買買提明·艾尼2

(1.新疆大學(xué)數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 新疆 烏魯木齊 830046；2.新疆大學(xué)機(jī)械工程學(xué)院, 新疆 烏魯木齊 830047)

利用連續(xù)體離散成均勻分布SPH粒子的思路，以漸開線直齒輪為例，通過其各部分曲線進(jìn)行研究，確定齒輪齒根過渡曲線以及其他部分的參數(shù)方程，提出了三維齒輪和一對嚙合齒輪模型離散成SPH粒子的方法和算法，編程實(shí)現(xiàn)了正確嚙合一對齒輪模型離散成SPH粒子的SPH前處理程序。通過用SPH方法對齒輪傳動中沖擊、接觸等問題進(jìn)行計(jì)算和模擬，提供了一種較好的前處理平臺。

SPH方法；齒輪過渡曲線；離散粒子；前處理

0 引言

在齒輪傳動的研究中，準(zhǔn)確繪制齒輪的三維模型是機(jī)械設(shè)計(jì)中的一項(xiàng)重要工作。其中，齒輪的三維模型造型中齒廓線的繪制是難點(diǎn)，主要原因是齒廓曲線較為復(fù)雜,尤其是齒根部分的過渡曲線與齒輪的加工方式有關(guān),需要具體分析[1]。

兩齒輪高速旋轉(zhuǎn)嚙合工作時,由于兩齒輪輪齒的相互進(jìn)入嚙合和脫離嚙合等因素，將會帶來較大的沖擊。這種沖擊會造成齒輪接觸強(qiáng)度或彎曲強(qiáng)度超出許用值，從而導(dǎo)致齒面點(diǎn)蝕或齒根斷裂，嚴(yán)重影響機(jī)械結(jié)構(gòu)正常運(yùn)行，甚至?xí)磔^大的事故。為此，對沖擊問題的正確評價是至關(guān)重要的。沖擊問題是一種比較復(fù)雜、求解難度大的非線性工程問題。用解析方法和簡單的數(shù)值方法計(jì)算和仿真時，一般都會對問題做大量的簡化，不能正確地反映實(shí)際情況[2]?；诰W(wǎng)格方法分析時，用固定網(wǎng)格很難實(shí)時追蹤齒輪傳動過程。特別是，當(dāng)面臨大變形和裂紋擴(kuò)散時，出現(xiàn)網(wǎng)格纏結(jié)和扭曲，在求解過程中需要實(shí)時進(jìn)行跟蹤和網(wǎng)格再劃分，這樣不僅計(jì)算費(fèi)用昂貴、費(fèi)時，而且會使計(jì)算精度受損[3-5]。近年來，光滑粒子流體動力學(xué)(SPH)方法被廣泛應(yīng)用在流體和固體力學(xué)領(lǐng)域。SPH方法的基本思想是將齒輪模型離散成有限多個有密度、體積、位置和速度等屬性的粒子[6-7]，計(jì)算過程中粒子之間不需要網(wǎng)格連接，而通過核函數(shù)將它們聯(lián)系起來，從而保證了物體的連續(xù)性，同時滿足了能量、動量和質(zhì)量的守恒性。由于無網(wǎng)格方法可以完全拋開網(wǎng)格并用拉格朗日方法來描述，因此，可直接對物體的實(shí)時運(yùn)動和變形過程進(jìn)行追蹤和坐標(biāo)更新，而且容易編程實(shí)現(xiàn)，從而避免網(wǎng)格奇異性帶來的計(jì)算誤差，保證動態(tài)變形和沖擊破碎的計(jì)算精度。用SPH方法對齒輪傳動過程進(jìn)行數(shù)值計(jì)算和模擬時，對齒輪傳動系統(tǒng)正確和合理離散成SPH粒子模型是非常重要的。

在此，根據(jù)連續(xù)體離散成均勻分布SPH粒子思路，以漸開線齒輪為例確定齒輪齒廓和輪齒過度曲線參數(shù)方程，從而確定齒輪上每一個離散粒子點(diǎn)的原始坐標(biāo)，根據(jù)齒輪嚙合原理，通過變換實(shí)現(xiàn)一對嚙合齒輪離散成SPH粒子，建立了三維SPH離散粒子模型，用RasMol軟件來顯示所作的模型[8]，并通過所采集齒根、齒廓和齒頂上點(diǎn)的粒子坐標(biāo)，將它們光滑連接，以檢驗(yàn)?zāi)Ｐ偷恼_性。

1 連續(xù)體離散成SPH粒子的思路

將連續(xù)體離散成均勻分布的SPH粒子的基本思路是：當(dāng)物體的形狀為長方形時，相互距離為粒子間距的平面束，先將物體截成若干個矩形截面，然后計(jì)算出每個方向能容納的粒子數(shù)，用粒子數(shù)、每個方向的長度和粒子間距來確定粒子的坐標(biāo)(如圖 1a 所示)。

當(dāng)物體的形狀為球形、圓柱形和錐形時，相互距離為粒子間距的平面束，先將物體截成若干個圓截面，然后將每個圓截面分解成以一個粒子為中心，半徑相差粒子間距的一系列同心圓面，其次計(jì)算每個圓的周長和能容納的粒子數(shù),最后用粒子數(shù)、每個圓的周長和粒子間距確定粒子的坐標(biāo)(如圖 1b所示)。

對于齒輪而言，相互距離為粒子間距的平面束，先將齒輪沿著其厚度截成若干個齒輪截面，用齒廓曲線方程計(jì)算齒廓線上的點(diǎn)與其中心線上點(diǎn)的相應(yīng)弧長，來確定每個粒子的坐標(biāo)(如圖 1c所示)，具體實(shí)現(xiàn)方法將在下面討論。

圖1 連續(xù)體離散SPH粒子

2 齒廓曲線

2.1 齒廓曲線的漸開線方程

由漸開線的形成原理可知，用rb表示基圓半徑，rk和θk分別表示漸開線上k點(diǎn)的極徑和極角[9]

(1)

αk為k點(diǎn)的壓力角。

2.2 齒根過渡曲線方程

用齒條型刀具加工齒輪時，齒條與齒輪的無間隙嚙合，刀具的直線部分切除被加工的齒輪齒廓的漸開線部分，過渡曲線部分由刀具的圓角部分切除，加工過程中，刀具的加工節(jié)線與齒輪加工節(jié)圓發(fā)生純滾動。

刀具齒廓的參數(shù)和齒根過渡曲線的參數(shù)方程[1]為：

(2)

rp為刀頂圓角半徑；a為刀頂圓角圓心Cp距中線的距離；b為刀頂圓角圓心Cp距刀具齒槽中心線的距離；f為齒高系數(shù)；c為徑向間隙系數(shù)。

(3)

2.3 齒輪漸開線與齒根過渡曲線的光滑連接

當(dāng)極軸與橫坐標(biāo)、極點(diǎn)與原點(diǎn)重合時，由式(1)確定漸開線上的點(diǎn)，旋轉(zhuǎn)θ2才能與由式(3)確定出的齒根過渡曲線，連接為一條完整的光滑齒廓曲線。旋轉(zhuǎn)角度可由式(4)確定，即

(4)

z為齒數(shù)。

3 輪齒截面的離散化

離散化輪齒截面時，如圖2所示，先對齒根過渡曲線與垂直方向之間的部分進(jìn)行離散化，再對漸開線與垂直方向之間的部分進(jìn)行離散化。處理此問題的關(guān)鍵，是要確定對應(yīng)齒廓曲線上點(diǎn)的半徑與垂直方向之間的夾角與弧長。

圖2 輪齒截面的離散化過程

3.1 齒根過度曲線與輪齒中心線之間的離散化

根據(jù)齒根過渡曲線參數(shù)方程(3)，可得齒根過渡曲線上對應(yīng)于α′的半徑為：

(5)

當(dāng)α′=π/2時，R(α′)的值最小，記為Rmin;當(dāng)α′=α?xí)r，R(α′)的值最大，記為Rmαx。

假設(shè)粒子之間的距離為d，則R(α′)從Rmin每一次增加d到Rmax時，每一個R(α′)對應(yīng)的α′由式(5)來確定。此時式(5)變成一個非線性方程。本研究用牛頓迭代法求解此非線性方程[10]，確定α′近似數(shù)值，計(jì)算齒根過渡曲線上的點(diǎn)，然后計(jì)算此點(diǎn)與輪齒中心線之間的夾角和弧長。

3.2 漸開線與輪齒中心線之間的離散化

漸開線上的點(diǎn)對應(yīng)的半徑，從齒根過渡曲線最后一個粒子對應(yīng)的半徑Rlast開始，每一次增加d，直到齒頂圓半徑Ra。當(dāng)半徑為Rd時，由式(1)可得對應(yīng)這個半徑的壓力角αd為：

αd=arccos(rb/Rd)

(6)

由極坐標(biāo)與直角坐標(biāo)的關(guān)系，對應(yīng)半徑Rd漸開線上點(diǎn)的直角坐標(biāo)為：

(7)

式(7)確定的點(diǎn)，逆時針旋轉(zhuǎn)θ2后，可計(jì)算此點(diǎn)與輪齒中心線之間的夾角和弧長，具體流程如圖3所示。

圖3 輪齒離散化流程

4 一對嚙合齒輪模型SPH離散粒子的建立

用以上所述的離散化思路，確定齒輪上每個SPH離散粒子點(diǎn)的原始坐標(biāo)。要實(shí)現(xiàn)一對嚙合齒輪的SPH離散粒子模型，有必要將兩輪的每個SPH離散粒子點(diǎn)的坐標(biāo)通過變換，平移和旋轉(zhuǎn)到嚙合位置。主動輪SPH離散粒子點(diǎn)的坐標(biāo)變換方程為：

(8)

從動輪離散坐標(biāo)的變換方程為：

(9)

Y為2個嚙合齒輪中心之間的距離；xi1,yi1,zi1為主動輪的SPH離散粒子點(diǎn)的原始坐標(biāo)；xi2,yi2,zi2為從動輪的SPH離散粒子點(diǎn)的原始坐標(biāo)；θ1=2π/(4z1)，θ2=2π/(4z2)，z1和z2分別為主動輪和從動輪的齒數(shù)。

在以上研究的基礎(chǔ)上，編程實(shí)現(xiàn)一對嚙合齒輪離散SPH粒子程序。在程序中輸入齒輪的模數(shù)、齒數(shù)等參數(shù)及粒子間距，將生成RasMol能識別的有每個SPH離散粒子點(diǎn)的坐標(biāo)信息文件，用RasMol打開此文件，就能看見如圖4a所示的模型(模數(shù)為10mm，齒數(shù)分別為15，30，粒子間距為0.01m)。

將主動輪和從動輪的SPH離散粒子點(diǎn)粒子的坐標(biāo)采集出來，將它們光滑連接，發(fā)現(xiàn)2個齒輪接觸在基圓上，如圖4b 所示。

圖4 嚙合齒輪模型

5 結(jié)束語

為了用SPH方法對齒輪傳動過程進(jìn)行計(jì)算和模擬，建立了一對嚙合齒輪的SPH離散粒子模型，確定了齒根過渡曲線、齒廓曲線等的參數(shù)方程，編制了實(shí)現(xiàn)一對嚙合齒輪離散粒子的SPH前處理程序。對模型進(jìn)行離散的結(jié)果顯示，此研究可以將一對嚙合齒輪正確離散成SPH粒子，并保證離散后正確嚙合。為今后用SPH方法分析齒輪傳動中沖擊、接觸等問題，提供了一種較好的前處理平臺，具體計(jì)算和分析方面的研究將在后續(xù)的文章中給出。

[1] 吳繼澤,王 統(tǒng).齒根過渡曲線與齒根應(yīng)力[M].北京:國防工業(yè)出版社, 1989.

[2] 陳 霞,夏巨諶,胡國安,等.直齒圓錐齒輪嚙合過程數(shù)值模擬[J].機(jī)械設(shè)計(jì), 2006, 23(4): 21-23.

[3] 賈光輝, 黃 海,胡震東.超高速撞擊數(shù)值模擬仿真結(jié)果分析[J].爆炸與沖擊, 2005,25(1): 47-54.

[4] Hwang S C, Lee J H,Lee D H,et al. Contact stress analysis for a pair of mating gears[J].Mathematical and Computer Modeling,2013, 57(1/2):40-49.

[5] 唐進(jìn)元,劉 欣,戴 進(jìn).基于ANSYS/LS-DYNA的齒輪傳動線外嚙合沖擊研究[J].振動與沖擊, 2007,26(9):40-44.

[6] Liu M B, Liu G R.Smoothed particle hydrodynamics (SPH):an overview and recent developments[J].Archives of Computational Methods in Engineering, 2010,17(1):25-76.

[7] Kikuchi M, Miyamoto M.Numerical simulation of impact crush/buckling of circular tube using SPH method[J].Key Engineering Materials, 2004, 261/263: 697-702.

[8] Bernstein H J. Manual of RasMol 2.7.5(molecular graphics visualization tool) [EB/OL]. [2014-06-15].http://www.openrasmol.org/doc/.

[9] 孫 恒，陳作模.機(jī)械原理[M].北京:高等教育出版社, 1999.

[10] Rahman K, Yimin R, Helil N.A fourth order iterative method for solving nonlinear equations[C]//IEEE International Conference on Computer Science and Information Technology, 2010, 8: 628-630.

Research on a Pair Meshing Gear Model Discrete into SPH Particle

RAHMATJANImin1,MAMTIMINGeni2

(1.College of Mathematics and System Science, Xinjiang University, Urumqi 830046, China；2.School of Mechanical Engineering, Xinjiang University, Urumqi 830047, China)

In this paper,` based on continuums discrete into SPH particles idea ,take involute spur gear as an example, the parameter equations of the transition curve of gear root and the other parts of gear are decided by studied every part curve of gear, provide a algorithm of discrete into SPH particles for 3D gear model and a pair of meshing gear model, programmed pre-processor program for correct meshing a pair of gear, to calculate and simulate impact ad contact in gear transmission to provide a good pre-processor platform.

SPH method;gear transition-curve; discrete particles; pre-processor

2014-06-30

國家自然科學(xué)基金資助項(xiàng)目(51075346)

TH112;TH132.4

A

1001-2257(2014)10-0008-04

熱合買提江·依明(1974-)，男，新疆伊犁人，講師，研究方向?yàn)闄C(jī)械優(yōu)化設(shè)計(jì)與數(shù)值模擬研究；買買提明·艾尼(1958-)，男，新疆烏魯木齊人，教授，博士，博士研究生導(dǎo)師，研究方向?yàn)闄C(jī)械設(shè)計(jì)及理論、計(jì)算力學(xué)。