雙環(huán)減速器環(huán)板內齒接觸分析及齒輪系統(tǒng)模態(tài)分析
3.1 雙環(huán)減速器環(huán)板內齒接觸有限元分析
因為在少齒差內嚙合中祠够,不處于嚙合位置的齒對在進入嚙合之前以及在脫離嚙合之后袍嬉,其內、外齒廓間的間隙非常小旦委,因此雙環(huán)減速器環(huán)板內齒輪嚙合傳動在載荷作用下,因彈性變形實際接觸齒對數(shù)遠大于理論重合度喧务,在傳遞載荷的過程中仰税,輪齒的變形量要大于一部分齒對的間隙,這些齒對就要接觸并同時分擔載荷割坠,這就大大提高了整個齒輪傳動裝置的承載能力齐帚。同時實際齒輪的接觸對數(shù)及接觸應力,隨著載荷增大而相應增加彼哼。因而少齒差內嚙合行星齒輪傳動的實際接觸齒對數(shù)的確定对妄,對齒輪承載能力的估算以及齒輪模數(shù)的正確確定具有重要意義。對于齒輪的接觸問題的敢朱,國內外學者對此已作了不少研究剪菱,取得一些研究成果。然而這些研究有的沒有考慮誤差的影響拴签,有的僅考慮內外嚙合齒輪副齒廓理論間隙渗某,沒有考慮輪齒制造誤差和輪齒的彈性變形的影響。本文對于環(huán)板內齒輪接觸問題進行了計算喇暴,計算時考慮內外嚙合齒輪副齒廓理論間隙牍系、制造誤差及輪齒彈性變形的影響,在此基礎上倍熄,建立了少齒差內嚙合行星齒輪傳動實際接觸齒對數(shù)及各齒間載荷分配的理論分析計算的實體模型户寺,利用I-DEAS軟件進行了環(huán)板內齒輪接觸有限元計算,并用該軟件分析計算了在不同載荷工況下的實際接觸齒對數(shù)混检、最大接觸應力以及兩環(huán)板上各接觸齒的載荷分配情況抢督。
3.1.1 誤差對直齒內嚙合傳動重合度的影響分析
通常,內嚙合漸開線齒輪的理論計算是按無側隙的情況設計的写并,但要保證齒輪傳動靈活喧雹,不發(fā)生卡滯,必須保留足夠的齒側間隙于嚼,在實際制造與安裝中齒輪系統(tǒng)都必然存在尺寸誤差越名,這些誤差直接影響著齒輪傳動的精度、齒輪的強度和振動噪聲。
眾所周知渺蒿,內嚙合直齒圓柱齒輪傳動的理論重合度為:
設內齒直齒輪傳動的齒輪模數(shù)為m痢士,分度圓壓力角為a,外齒輪齒數(shù)為Z1茂装,內齒輪齒數(shù)為Z2怠蹂,外、內齒輪變位系數(shù)分別為X1少态、X2城侧,外、內齒輪變位系數(shù)分別為S1′彼妻、S2′嫌佑,外、內齒輪節(jié)圓上的齒槽寬分別為e1′澳骤、e2′歧强,兩嚙合齒輪節(jié)圓上的圓周側隙為jt,兩嚙合齒輪非工作側齒廓沿公法線方向的側隙(法向側隙)為jn为肮。顯然有

(3.2)
聯(lián)系式3.6與式3.8,從兩工中可以看出肤京,計算中心距a′與嚙合角

′是一一對應關系果喊,這樣式3.6可以表達為三種形式:
(1)由確定的齒側間隙jn與齒輪變位系數(shù)X2-X1來求中心距a′;
(2)由確定的中心距a′與齒輪變位系數(shù)X2-X1來求齒側間隙jn腿逞;
(3)由確定的齒側間隙jn與中心距a′來求齒輪變位系數(shù)X2-X1实朗。
當齒輪傳動的幾何設計參數(shù)確定后,影響實際傳動的就只有制造和安裝誤差汤钻,對于中心距a′而言垮撇,影響它的因素是制造誤差、安裝誤差荡唾,而不是齒側間隙jn與齒輪變位系數(shù)X2-X1状堰。齒輪設計是按無側隙進行的,同時不考慮制造矩允、安裝誤差糯驯,但加工時,齒輪齒厚公差都是負偏差赵椰,即齒輪厚公差都是負偏差耙福,即齒輪必須減薄,這樣實際傳動存在齒側隙姜钳。如果不考慮中心距的誤差坦冠,由式3.6的第3種形式可知,齒側隙的變化必然影響到齒輪變位系數(shù)的變化哥桥。顯然辙浑,這樣的變化也就會使齒輪的齒頂圓產生變化激涤,相應的齒頂圓奢力角產生變化,進而影響齒輪嚙合重合度例衍。
由式3.5變化為
本文實驗樣機的環(huán)板齒輪設計參數(shù)見表3.1所示昔期,根據本文實驗樣機的環(huán)板齒輪設計參數(shù),假定嚙合角不變化佛玄,則對應齒側隙的變化硼一,應用MATLAB繪圖分析,得到了齒輪變位系數(shù)隨齒側隙變化時的變化規(guī)律圖梦抢,見圖3.1所示般贼。
表3.1 內嚙合齒輪參數(shù)表
|
外齒輪 |
環(huán)板內齒 |
模數(shù) |
2 |
齒數(shù) |
56 |
58 |
嚙合角 |
39.039° |
變位系數(shù) |
1.37 |
1.685 |
齒頂高系數(shù) |
0.7 |
中心距 |
2.42 |
齒頂圓直徑 |
119.519 |
120.774 |
計算重合度 |
1.124 |
由3.2圖可以看出,當齒側隙jn由0變化至.25毫米時(查資料取較大最小側隙參考值)奥吩,齒輪變位系數(shù)差X2-X1的變化量由近0.32變化到0.50哼蛆,變化了約0.18。
對于齒頂圓直徑按德國工業(yè)標準DIN的簡單方法有:

(3.10)(3.11)
由式3.10划搓、3.11可以看出沮念,當總變位系數(shù)

產生變化時,必然引起起嚙合齒輪齒頂圓直徑的變化验脐。按X
1澜碎、X
2變化最小原則,有三種可能性织活,即X
1不變嫡境,X
2變化0~0.18;或X
1變化0~0.18稼那,X
2不變褒饱;或X
1變小,X
2變大搞吱,變化范圍0~0.09援儡。假定齒輪嚙合角不變,按本文實驗樣機參數(shù)纸牌,當齒輪變位系數(shù)發(fā)生變化時瘾显,其變化量對重合度的影響狀況見圖3.2所示。從圖3.2所示可以看出涎才,隨外齒輪變位變小時鞋既,齒輪嚙合重合度減小,隨內齒輪變位變小時耍铜,齒輪嚙合重合度增大邑闺,從圖中的變化趨勢為看,外齒輪變位系數(shù)的變化對重合度的影響程度與內齒輪變位系數(shù)的變化對重合度的影響程度大約是相同的,只是趨向是反向的陡舅。
對于中心距而言抵乓,查資料表35.2-60可知中心距的極限偏差為±0.011mm,假定外齒輪齒頂圓壓力角不變,內齒輪齒頂圓壓力角隨變位系數(shù)變化靶衍,對于中心距變化對重合度的影響變化趨勢及變位系數(shù)變化對重合度的影響的變化趨勢比較情形見圖3.3灾炭。假定內齒輪齒頂圓壓力角不變,外齒輪齒頂圓壓力角隨變位系數(shù)變化颅眶,對于中心距變化對重合度的影響變化趨勢及變位系數(shù)變化對重合度的影響的變化趨勢比較情形見圖3.4所示蜈出。

從圖3.3、3.4明顯看出中心距變化對重合度的影響程度明顯小于變位系數(shù)變化對重合度的影響程度涛酗。
齒輪制造誤差中對齒廓間隙影響較大的有最小齒側間隙净处、基節(jié)偏差、齒形誤差居鸳、齒距偏差婉固、齒厚偏差等等,但這些因素相互關聯(lián)常彰、相互制約稼乃、相互影響,特別是基節(jié)偏差盼涵、齒形誤差砰吻、齒距偏差、齒厚偏差等最終影響齒側隙的變化彤擒,因此最終影響齒輪傳動的主要因素仍是中心距和齒側隙的變化。分析圖3.2垢雨、3.3锯叮、3.4,在允許的制造誤差范圍內及最小齒側隙的情況下辐践,重合度的理論計算值變化不是很大坦报。
3.1.2 接觸分析中制造安裝誤差的引入
從上述分析可知狂鞋,齒輪的制造片择、安裝誤差對齒輪嚙合的重合度有一定的影響,但就理論計算值而言骚揍,從上面的分析結果可以看出其影響程度是不大的字管,然而在少齒差內嚙合齒輪副中,嚙合線附近兩相鄰工作齒對齒側間隙非常小信不。假設內齒輪固定嘲叔,當輪齒承擔載荷時,兩接觸輪齒將分別產生彈性變形,外齒輪將繞其軸線旋轉一附加角硫戈,當這一附加轉角足夠大時锰什,除原有接觸齒對繼續(xù)發(fā)生接觸外,相鄰工作齒對的齒側間隙也會消失丁逝,這樣就產生了多齒承載接觸汁胆。
在本節(jié)中,引入制造霜幼、安裝誤差撰类,用有限元方法計算內齒傳動的接觸對數(shù),同時驗證上節(jié)中討論的影響重合度的關鍵參數(shù)潘谴。用I-DEAS軟件計算環(huán)板內齒接觸問題英嬉,其關鍵是齒輪輪齒的實體建模。對于相嚙合的齒輪份肠,因其兩工作齒廓側面相互位置誤差的變化導致了實際齒廓間隙的變化,在嚙合線附近只有少數(shù)幾對齒在我們討論范圍內贬嚷,這里我們只考慮對齒廓間隙影響較大的最小齒側間隙僧少、基節(jié)偏差、齒形誤差供撵、齒距偏差、齒厚偏差等旱谐。在用用I-DEAS軟件進行齒輪輪齒的實體建模時蔗括,對輪齒而言,沿圓周是均勻分布的撤防,因此,齒輪的制造誤差在實體造型時寄月,表現(xiàn)在齒輪輪齒的齒厚及基圓尺寸的誤差上辜膝。而經多元回歸分析,得出各變量對齒廓間隙的內在影響關系表明:齒輪的基節(jié)偏差是最主要的影響因素漾肮。
本文雙環(huán)減速器齒輪加工精度為7級厂抖,根據資料克懊,查得齒輪的基節(jié)極限偏差△fpb為±0.014忱辅,對于齒厚極限偏差,根據資料保檐,查得齒輪齒厚的極限偏差△Es的上偏差為-128μm崔梗,下偏差為-192μm,將這些影響因素的極限偏差垒在,通過幾何關系的轉換蒜魄,換算到基節(jié)上愚贩,通過參數(shù)方程繪出齒輪單個輪齒的輪廓線盯辅,再由單個輪齒的齒廓沿圓周接齒輪齒數(shù)均勻復制,這樣得到整個齒輪的齒廓線米搭,由齒廓線延伸成齒輪實體模型趴鹰。
由于零部件的制造誤差,同時因裝配原因避纤,在進行整體裝配時纠沉,產生零部件間的安裝誤差,這些誤差直接影響齒輪間的嚙合狀態(tài)岔雾。對于雙環(huán)減速器中的環(huán)板內齒輪嚙合傳動,由于環(huán)板內齒中心與輸出齒輪中心之間的偏心距尺寸較卸辍(本文雙環(huán)減速器的偏心距為2.42mm) 黑乞,因此,偏心軸的偏心距制造誤差及中心距制造安裝誤差是影響環(huán)板內齒嚙合傳動的關鍵膛姊。本文雙環(huán)減速器兩偏心軸偏心距的制造誤差為±0.011 拟祖。
綜合以上誤差分析,將環(huán)板與輸出齒輪軸按極限偏差尺寸所定的位置進行輪齒嚙合裝配擦酌,然后進行有限元接觸分析俱诸。
3.1.3 齒輪有限元模型的建立及邊界條件的確定
環(huán)板內齒輪與輸出軸齒輪之間的傳動是直齒輪傳動,為保證進行接觸計算時赊舶,有較多的有較接觸區(qū)域乙埃,齒輪輪齒部分的有限元網格采用映射網格劃分法(Mapped Meshing )對齒輪進行有限元網格劃分。同時锯岖,為保證兩對嚙合齒輪所劃分的網格接觸對相應的面能對應齊,因此齒輪嚙合的接觸面相應應當齊整甫何。用輸出軸齒輪的齒頂面所在的輪廓曲面對環(huán)板內齒進行劃分(partition 命令)出吹,同理用環(huán)板內齒的齒頂面所在的輪廓曲面對輸出軸齒輪進行劃分,這樣得到對應齊整的齒輪嚙合的接觸面辙喂,如圖3.5所示捶牢。根據圖3.5所示,采用映射網格劃分法( Mapped Meshing)進行網格劃分的輪齒部分即是圖中輸出齒輪齒頂圓與環(huán)板內齒齒頂圓之間所劃分的部分巍耗,而其余部分采用自由網格劃分法(Free Meshing)進行劃分秋麸,劃分網格后的部分網格圖見圖3.6所示渐排。考慮實際接觸齒對數(shù)不會太多筒臂,所以每塊環(huán)板上只給出6個輪齒進行接觸計算烤酌。環(huán)板嚙合的實際過程是環(huán)板齒輪帶動輸出齒輪運動,計算時為方便脚自,假設環(huán)板固定不動惫康,設定輸出軸繞軸線轉動,將減速器的輸出功率換算為輸出扭矩联缝,加在輸出軸的輸出端赢瘦。

3.1.4 接觸計算結果分析
按考慮偏心距誤差影響及考慮齒厚誤差影響的兩個極限狀態(tài),組合為十種情況進行計算苇葫。十種計算工況見表3.2所示汇泰。考慮齒厚誤差影響的兩個極限狀態(tài)分別是齒輪的輪齒最薄堰聪、最厚柄立,偏心距誤差的兩個極限狀態(tài)分別是偏心軸的偏心距最大、最小兩種情況曾蚊。種接觸情況的計算結果分別由表3.3~3.7列出赃阀。
3.2計算工況列表
序
號 |
工況說明 |
序
號 |
工況說明 |
齒厚狀態(tài) |
偏心距狀態(tài) |
兩環(huán)板狀態(tài) |
齒厚狀態(tài) |
偏心距狀態(tài) |
兩環(huán)板狀態(tài) |
1 |
最厚 |
無誤差 |
對稱 |
6 |
最薄 |
最大正偏差 |
對稱 |
2 |
最薄 |
無誤差 |
對稱 |
7 |
最厚 |
最大負偏差 |
不對稱 |
3 |
最厚 |
最大負偏差 |
對稱 |
8 |
最薄 |
最大負偏差 |
不對稱 |
4 |
最薄 |
最大負偏差 |
對稱 |
9 |
最厚 |
最大正偏差 |
不對稱 |
5 |
最厚 |
最大負偏差 |
對稱 |
10 |
最薄 |
最大正偏差 |
不對稱 |
3.4工況一接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
一 |
200 |
1 |
1 |
75.2 |
160 |
1 |
1 |
66.8 |
120 |
1 |
1 |
56.1 |
80 |
1 |
1 |
41.1 |
40 |
1 |
1 |
28.0 |
3.5工況二接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
二 |
200 |
1 |
1 |
75.2 |
160 |
1 |
1 |
66.8 |
120 |
1 |
1 |
56.1 |
80 |
1 |
1 |
41.1 |
40 |
1 |
1 |
28.0 |
3.6工況三接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
三 |
200 |
1 |
1 |
68.2 |
160 |
1 |
1 |
62.6 |
120 |
1 |
1 |
55.1 |
80 |
1 |
1 |
40.2 |
40 |
1 |
1 |
25.3 |
3.7工況四接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
四 |
200 |
2 |
2 |
74.6 |
160 |
2 |
2 |
62.1 |
120 |
2 |
2 |
49.0 |
80 |
2 |
2 |
41.8 |
40 |
1 |
1 |
26.6 |
3.8工況五接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
五 |
200 |
2 |
2 |
2.2 |
160 |
2 |
2 |
58.6 |
120 |
2 |
2 |
45.6 |
80 |
2 |
2 |
42.1 |
40 |
1 |
1 |
26.5 |
3.9工況六接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
六 |
200 |
1 |
1 |
76 |
160 |
1 |
1 |
63.2 |
120 |
1 |
1 |
49.7 |
80 |
1 |
1 |
40.5 |
40 |
1 |
1 |
28.5 |
3.10工況七接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
七 |
200 |
1 |
3 |
76 |
160 |
0 |
2 |
63.2 |
120 |
0 |
2 |
49.7 |
80 |
0 |
1 |
40.5 |
40 |
0 |
1 |
28.5 |
3.11工況八接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
八 |
200 |
1 |
2 |
130.0 |
160 |
0 |
2 |
114.0 |
120 |
0 |
2 |
93.4 |
80 |
0 |
2 |
70.2 |
40 |
0 |
2 |
39.0 |
3.12工況九接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
九 |
200 |
2 |
1 |
104.0 |
160 |
2 |
0 |
92.6 |
120 |
2 |
0 |
86.6 |
80 |
2 |
0 |
65.3 |
40 |
1 |
0 |
42.8 |
3.13工況十接觸情況計算表
誤差情況 |
載荷(N·m) |
接觸齒對數(shù) |
最大接觸應力(MPa) |
環(huán)板1
|
環(huán)板2 |
工
況
十 |
200 |
2 |
1 |
136.0 |
160 |
2 |
0 |
117. |
120 |
2 |
0 |
102.0 |
80 |
2 |
0 |
77.4 |
40 |
2 |
0 |
47.6 |
由表3.2~3.13可以看出,無論有無制造擎颖、安裝誤差的影響榛斯,隨著負荷的增加,齒輪實際嚙合接觸對數(shù)增加搂捧,但從計算結果分析證明驮俗,負荷對接觸對的影響小于誤差對接觸對的影響。對比表中的接觸對數(shù)和接觸應力允跑,數(shù)據表明當齒厚發(fā)生變化時王凑,實際接觸齒對數(shù)變化比較大,同時最大接觸應力也有明顯變化聋丝。從變化的趨勢上看出與前面的理論推導變化趨勢相同索烹,這里只是更加明顯。
工況六至工況十這四和情況弱睦,是指兩環(huán)板之間不對稱百姓,即偏心軸的兩個偏心拐的尺寸誤差不一致,這里考慮的是特殊狀況况木,從表3.10-3.13看出垒拢,當兩環(huán)板不對稱時,將產生嚴重的載荷不均現(xiàn)象靡勾,載荷偏向中心距小的環(huán)板內齒上北拔,甚至只有一個環(huán)板承受載荷换秧。由此說明,對于雙環(huán)減速器偏心軸有加工以搏,應特別注意兩曲捌的制造誤差辫航,應盡量使安裝好的兩環(huán)板具有較好的對稱性。圖3.7~圖3.12為六種工況下環(huán)板內齒的接觸齒對的應力云圖鲜附。





3.2 雙環(huán)減速器有限元模型建立及模態(tài)分析
近年來钙蕉,隨著齒輪振動噪聲研究的不斷深入,不僅需要考慮齒輪系統(tǒng)及其嚙合過程屋包,還必須考慮動態(tài)嚙合力在整個齒輪系統(tǒng)中的傳遞官孝,以及系統(tǒng)中各零部件的固有特性和動態(tài)響應的性質,因此厌圈,必須以整個齒輪系統(tǒng)為研究對象门俏,建立同時包括傳動系統(tǒng)和結構系統(tǒng)的完整系統(tǒng)的分析模型。本節(jié)將在雙環(huán)減速器運動特性分析的基礎上摩疑,以整個雙環(huán)減速器齒輪系統(tǒng)為對象危融,建立有限元動力分析模型,全面分析該系統(tǒng)的固有特性和動態(tài)響應雷袋。
齒輪系統(tǒng)的固有特性一般指系統(tǒng)的固有頻率和固有振型吉殃,是齒輪系統(tǒng)的動態(tài)特性之一,它對系統(tǒng)的動態(tài)響應楷怒,動載荷的產生與傳遞蛋勺,以及系統(tǒng)振動的形式等都具有重要意義。此外鸠删,固有特性還是用振型疊加法求解系統(tǒng)響應的基礎抱完。由于系統(tǒng)的固有特性表明了在哪些頻率下結構會產生共振以及在各階頻率下結構的相對變形,因此對于改善結構動態(tài)特性具有重要意義刃泡。通常巧娱,研究齒輪系統(tǒng)固有特性有理論計算法和實驗測試法兩種。實驗測試將在以后的章節(jié)中研究烘贴,本章將利用數(shù)值計算法計算雙環(huán)減速器齒輪系統(tǒng)的固有特性禁添。結構系統(tǒng)固有特性的數(shù)值計算通常采用有限元模態(tài)分析方法。在有限元分析中庸伏,結構的動力學特征用微分方程組表征度堤,則結構的固有頻率和固有振型相應地由動力學方程的特征值和特征向量所確定。所以可以說完簿,求系統(tǒng)的固有頻率和固有振型即為求系統(tǒng)的特征值和特征向量。特征值問題的解法很多桅蕊,主要有多項式割線迭代法聚灸、矢量逆迭代法娇皇、廣義雅可比法、子空間迭代法和行列式收縮法朱鹤。動力問題有限元法是求解復雜結構的大型動力學方程組的特征值問題的有效方法呛仁。本章將用動力有限元法對雙環(huán)減速器齒輪系統(tǒng)進行特征值分析,計算其固有頻率和固有振型延杯。
本文建立了齒輪系統(tǒng)的有限元動力學模型捅没,用I-Deas 集成化軟件的固有模態(tài)求解模塊求解了齒輪系統(tǒng)的固有頻率和振型,并用響應動力(Response Dynamios ) 分析模塊研究了齒輪傳動系統(tǒng)的動態(tài)響應特性玫斋。要進行有限元模態(tài)分析舀黄,必須先建立齒輪系統(tǒng)的有限元力學模型。
3.2.1 齒輪系統(tǒng)有限元模型的建立
通常蚯涮,實際工程結構都是質量治专、剛度和阻尼連續(xù)分布的連續(xù)系統(tǒng),其動力學特性常常用偏微分方程組描述遭顶,但由于實際結構的幾何形狀及邊界條件復雜张峰,難以用解析法求解偏微分方程組。于是棒旗,通常對實際工程結構離散化喘批,將連續(xù)分布參數(shù)系統(tǒng)簡化為具有有限多個自由度的集中參數(shù)系統(tǒng)。用矩陣表達式對每個單元進行單元分析铣揉,根據虛位移原理建立單元的動態(tài)平衡方程式饶深;然后,作整體分析老速,得到整個結構系統(tǒng)的動力學方程式粥喜。應用自由振動邊界條件,得到由有限元法計算整個結構系統(tǒng)固有振動特性的基本方程橘券。這樣额湘,通過離散化模型,可將復雜的偏微分方程組簡化為可以求解的常微分方程組荞诡,從而使問題的求解大大簡化忘哼。如前所述,雙環(huán)減速器結構運動特性復雜灸尾,既有旋轉齒輪運動锨穷,又有連桿平面運動,根據這一具體情況汗势,本文用集中參數(shù)法建立了該傳動系統(tǒng)的彎——扭耦合振動的動力學分析模型岗憨。對箱體建立了有限元分析模型,該有限元模型的建立是用SDRC公司的I-Deas大型集成化軟件完成锁熟,并在此基礎上進行動力學響應分析計算帝膊。當進行系統(tǒng)的模態(tài)計算時确奄,為了模型建立的方便,則將齒輪輪齒部分簡化為分度圓柱僻跳,將齒輪和輪毅及輪毅與軸的緊配合聯(lián)接看成剛性聯(lián)接闪铸,忽略鍵槽、倒角和退刀槽的影響茶月,利用自由網格劃分法(Free Meshing)對齒輪進行有限元網格劃分需五。圖3.17為輸出齒輪軸的結構示意圖,齒輪的有關參數(shù)見表2.1轧坎。物理參數(shù)為:彈性模量E = 2.06×l011N/m2宏邮,泊松比為0.3, 材料密度p=7.8×103kg/m3。

采用四節(jié)點四面體單元眶根,利用軟件自由網格劃分法(Free Meshing )對輸出齒輪軸進行有限元網格劃分蜀铲。減速箱的其它零件采用相同的方法進行有限元建模,圖3.18為輸出齒輪軸的有限元網格模型共劃分為3021個單元属百,835個節(jié)點记劝,圖3.19為偏心軸的有限元網格模型,共劃分為4098個單元族扰,1188個節(jié)點厌丑,圖3.20 為環(huán)板的有限元網格模型,共劃分為1505個單元渔呵,594個節(jié)點怒竿。
雙環(huán)減速器箱體是用鑄鐵鑄造而成,不考慮或對分析影響甚微的結構特征扩氢,如小孔耕驰、園角、倒角椎裕、螺栓孔等言丧,在進行減速箱實體建模時,忽略上述因素的影響皿进,同時將箱蓋與箱體合為一體辕寺,而齒輪箱與前后端蓋、軸承組合而成臊瞬,結構比較復雜误扯,建模時不考慮頂蓋與箱體、前后端蓋银景、軸承與箱體的接合部的影響骇瓦。軸承采用相同尺寸的圓環(huán)套,而其材料性質用具有與軸承等效剛度相同的材料參數(shù),這些處理不會對齒輪箱體的重量及剛度產生大的影響庶咨,完全能保證足夠的計算精度概丢。用自由網格劃分方法(Free Meshing),采用八節(jié)點四面體單元劃分網格民轴,圖3.21為所建立的齒輪箱有限元模型網格,節(jié)點數(shù)為8761球订,單元數(shù)為28286后裸。前面分別建立了齒輪系統(tǒng)主要零件和箱體的有限元網格,將簡化后的軸承用相同的方法進行有限元網格劃分冒滩,再用I-DEAS軟件的APPED命令將齒輪系統(tǒng)的各零部件集成為整個系統(tǒng)的有限元網格模型微驶,如圖3.22 所示。

齒輪系統(tǒng)整個模型的單元數(shù)為112072开睡,節(jié)點數(shù)為66120因苹。
劃分有限元網格后,還必須對整個網格模型進行檢查篇恒,從而保證計算結果的真實性和準確性扶檐。首先檢查自由單元邊,當單元的某一邊不在其它單元之內時胁艰,稱為自由單元邊款筑。在復雜模型的建立過程中,通過拉伸旋轉等操作產生的各個部件腾么,有時會沒有連接在一起浩出,這將導致有限元模型開裂,影響計算結果靡隔,嚴重時將使計算失敗扶蜻。其次檢查重復單元、重復節(jié)點谆纸,根據情況決定是否將它們合并在一起搔绿。最后檢查單元的形狀參數(shù),過度扭曲的單元將影響計算片迁,必須并將其修改為可以接受的形狀抖唧。
3.2.2 齒輪系統(tǒng)有限元模態(tài)分析
結構模態(tài)是振動系統(tǒng)特性的一種表征。一般來說煞肠,模態(tài)參數(shù)包括固有頻率瓤鬓、固有振型、模態(tài)質量乏束、模態(tài)剛度和模態(tài)阻尼比等四乱,但其中最重要的是前兩項。模態(tài)參數(shù)將表明在哪幾種頻率下結構會產生共振以及在各階頻率下結構的相對變形,對于改善結構的動態(tài)特性宣吱,這是最重要的基本參數(shù)窃这。
通常,結構的前幾階固有頻率較低征候,對結構的動力學特性影響較大杭攻,所以分析時只提取齒輪系統(tǒng)結構的前20階固有頻率及其對應的振型便能滿足工程需要。為了能反映箱體的實際情況疤坝,在進行模態(tài)分析時將箱體底座z方向約束兆解,地腳螺栓處用固定約束,研究系統(tǒng)的約束模態(tài)跑揉。利用I-DEAS軟件的固有模態(tài)求解模塊求解模型的固有頻率和固有振型的方法有:Lanczos法锅睛、Guyan Reduction法和同步向量迭代法(Simultaneous vecto: Iteraron,SVI)。Guyan法是決動力學問題的一種古典方法历谍,但計算結果的精度在很大程度上取決于模型中的主自由度的選取现拒。SVI法的精度不取決于用戶的判斷,但通常計算時間比Guyan法要長望侈。一般求結構的固有頻率和固有振型時用Lanczos法印蔬,Lanczos法是近年來新興的一種模態(tài)求解方法,其特征值惜肃、特征向量求解精度高喝壹。它的計算速度比另兩種方法快,不需要選取“主自由度”盲狈,且要求的輸入數(shù)據較少旱已。Lanczos法使用Stunn序列檢查,在用戶感興趣的頻率范圍內螺矮,在每個漂移點處如果不能找到所有的特征根殿漆,Lanczos法會給出提示信息,以避免漏根和引入不可能的根畅型。用Lanczos法求解特征值和特征向量時肾蕉,先根據載荷空間分布模式按一定規(guī)律生成一組相互正交的Lanezos向量,在將系統(tǒng)運動方程轉換到這組Lanczos向量空間以后乘澈,對運動方程進行縮減笤敞,然后通。過求解一次縮減了的運動方程的特征值問題尝鬓,再經過坐標系的變換瞻绝,進而就可得到原系統(tǒng)的部分特征解,從而避免了一般算法如反迭代法或子空間迭代法中的多次迭代步驟秸讹,因而檀咙,顯著提高了求解效率。Lanczos方法適應面廣泛,可適用于所有的大中型結構問題弧可,并都能給出很好的結果蔑匣。
與其它特征值問題求解方法一樣,Lanczos算法的主要目的是求解結構系統(tǒng)的動力學方程棕诵。

(3.12)
式中M裁良、C、K——系統(tǒng)的總質量校套、阻尼和剛度矩陣

——系統(tǒng)節(jié)點加速度趴久、速度和位移向量

——節(jié)點載荷向量
不考慮阻尼影響的多自由度系統(tǒng)自由振動方程為

(3.13)
設式(3.12)解的形式為

(3.14)
式中

——n階向量,位移向量

的幅值

——向量

振動的頻率
t0——由初始條件確定的時間常數(shù)
將式(3.14)代入式(3.13)得到特征方程為

(3.15)
或
KX=MXΩ2 (3.16)
式中搔确,

,

灭忠,X膳算、

分別為固有振型矩陣和固有頻率矩陣。
引入變換

(3.17)
式中攘体,q(t)=[q1好浆、q2 … qn]T,為廣義得位移向量办轮,qi(i=1浓朋,2,…枫欢,n)為廣義位移值磨搭。
將上述變換代入式(3.12)中,兩端前乘XT返引,則式(3.12)成為n個解耦的二階常微分方程

(3.18)
式中

——模態(tài)阻尼因子

——無阻尼固有頻率

——模態(tài)力缭亦,

,是載荷向量

在振型{X}
i上的投影隔每,若

按一定的空間分布模式而隨時間變化培穆,即

(3.19)
上工中的符號s表示空間坐標。
則求解

(3.20)
正則化得贼酵,

糙笛,其中
再求解

(3.21)
將式(3.20)的解正交化得

(3.22)
式中,
為了保證和提供計算的穩(wěn)定性烛愧,需要對由式(3.22)計算得到的

進行重正交處理油宜,即

(3.23)
式中,

屑彻,上式中的

即為式(3.22)中的

验庙,其迭計算終止于

滿足規(guī)定的誤差,或s達到某規(guī)定的值。然后將

恢復寫出為

粪薛,進行正則化計算得

悴了,其中

。
由式(3.21)违寿、(3.22)和

得到
(i=2湃交,3,…藤巢,且q0=0)

(3.24)
式中搞莺,A=K-1M。將上式寫成矩陣形式為
AQ=QT (3.25)
式中掂咒,Q=[q1 q2 … qr]抛现,矩陣T如下式所示
(3.26)
引入原特征向量和Lanczos向量間的變換X
r=QZ(Z為r×r階矩陣),將原廣義特征值方程KX=MX

兩邊用Q
TMK
-1前乘谊阐、用

后乘得到用Lanczos向量內三對角矩陣T表示的標準特征值問題贾反,即

(3.27)
求解標準特征值問題式(3.27)的特征解:

(3.28)
由此計算得到原問題的部分特征解

(3.29)
從而使對A的特征值求解轉化為對T 的求解。由此可見它蛔,Lanczos算法的最大優(yōu)點是其所得的三對角陣T的特征解直接近似于原矩陣A中r個特征根七萧。因此,對于僅需求解部分模態(tài)解的大型結構問題首袍,Lanczos法可以在子空間內非潮镉拢快速有效地求解。本課題對系統(tǒng)各主要零部件及整體系統(tǒng)固有特性的求解就采用了這種方法衰呢。
如圖3.22 所示的齒輪系統(tǒng)有限元模型夫凭,根據實際的安裝情況,為了能比較準確地反映實際狀態(tài)猪出,箱體的邊界條件為:箱體底座垂直方向約束顶赎,四個螺栓孔處全部固定約束。則研究齒輪系統(tǒng)的約束模態(tài)坤学,表3.14~3.15分別為用I-Deas軟件的Lanczos模態(tài)求解法計算得到的各重要零部件及齒輪系統(tǒng)前上10階固有頻率和響應的振型情況疯坤。限于篇幅,圖3.23給出了齒輪系統(tǒng)的前十階固有振型深浮。
表3.14 齒輪系統(tǒng)前10階固有頻率和振型
模型 |
固有頻率(Hz) |
振型特征 |
1 |
974.3714 |
箱體橫向(y向)彎曲 |
2 |
1219.1347 |
箱體橫向(x向)彎曲 |
3 |
1575.9162 |
輸入軸輸出軸彎曲+箱體扭擺 |
4 |
1690.4077 |
箱體底面向內凹陷 |
5 |
1733.7038 |
輸入軸輸出軸向下彎曲 |
6 |
1738.6459 |
輸出軸輸出軸彎向上彎曲 |
7 |
1871.7482 |
輸入軸輸出軸水平彎曲 |
8 |
2345.6277 |
箱體兩邊膨脹中間向下凹陷 |
9 |
2657.8228 |
環(huán)板扭擺 |
10 |
2706.0426 |
輸入軸輸出軸彎曲向上彎曲+箱體扭彎 |
表3.15 零件前10階固有頻率
模態(tài) |
固有頻率(Hz) |
環(huán)板 |
輸出軸 |
偏心軸 |
輸入軸 |
箱體 |
1 |
1097.70 |
3004.30 |
2037.79 |
4538.75 |
890.79 |
2 |
1490.23 |
3021.40 |
2039.17 |
4542.69 |
1317.35 |
3 |
2355.64 |
6231.67 |
5118.85 |
10006.96 |
1428.24 |
4 |
2638.10 |
7985.83 |
5126.75 |
10326.25 |
1510.36 |
5 |
2977.05 |
8009.92 |
5590.44 |
10329.09 |
1681.76 |
6 |
3034.42 |
9210.48 |
8503.90 |
15266.32 |
1725.52 |
7 |
3457.06 |
13634.45 |
8939.29 |
18645.97 |
1761.74 |
8 |
3754.52 |
13704.50 |
8961.50 |
18668.88 |
2010.17 |
9 |
4296.74 |
14793.46 |
11425.93 |
19068.83 |
2021.84 |
10 |
4781.41 |
15484.25 |
13581.56 |
27836.97 |
2059.41 |
3.3 本章小結
在分析誤差對直齒內嚙合傳動重合度的影響分析的基礎上压怠,考慮部分制造誤差及安裝的影響,對環(huán)板內齒的嚙合傳動飞苇,建立了考慮誤差的環(huán)板內齒接觸有限元模型菌瘫,并按不同的誤差狀態(tài),組合成十種不同的工況布卡,對不同的工況分別建立了有限元接觸分析計算模型雨让,通過計算結果分析討論了誤差及載荷變化的對齒輪接觸齒對的影響雇盖,指出了因載荷不均造成的減速器振動因素。
1)在雙環(huán)減速器齒輪系統(tǒng)的制造栖忠、安裝誤差中崔挖,齒輪側隙的誤差對齒輪重合度的影響大于中心距誤差對重合度的影響。而對于齒側隙的誤差大小重要是由于齒輪變位系數(shù)的變化造成的庵寞。即齒輪變位系數(shù)的誤差對齒輪重合度的影響大于中心距誤差對重合度的影響狸相。
2)對于兩環(huán)板與輸出齒輪的傳動,如果兩環(huán)板對于輸出齒輪中心不對稱朗溶,即當兩偏心軸上的兩偏心曲拐的誤差不一致時荡档,將產生兩環(huán)板嚴重的載荷不均現(xiàn)象。這是雙環(huán)減速器產生振動故障的重要因素彼窥。