杜雅寧 王士軍 齊 娜 王 程 鄭蘇芮

（山東理工大學(xué)機(jī)械工程學(xué)院，山東 淄博 255000）

諧波傳動(dòng)因其體積小、傳動(dòng)精度高、承載能力強(qiáng)等特點(diǎn)，被廣泛應(yīng)用于航空航天、工業(yè)機(jī)器人、精密機(jī)床等多個(gè)領(lǐng)域[1-2]。諧波傳動(dòng)主要依靠柔輪在波發(fā)生器的作用下產(chǎn)生周期性的波動(dòng)變形與剛輪進(jìn)行嚙合運(yùn)動(dòng)來(lái)實(shí)現(xiàn)運(yùn)動(dòng)的傳遞[3]，柔輪在與剛輪的嚙合運(yùn)動(dòng)過(guò)程中極易發(fā)生嚙合干涉導(dǎo)致柔輪發(fā)生應(yīng)力集中、輪齒磨損等現(xiàn)象，因此柔輪齒與剛輪齒間的無(wú)干涉嚙合是研究的重點(diǎn)。張立勇等[4]根據(jù)諧波傳動(dòng)的包絡(luò)理論，分析了諧波傳動(dòng)徑向變形量對(duì)嚙合特性的影響，并通過(guò)有限元法進(jìn)行驗(yàn)證。張寧等[5]發(fā)現(xiàn)柔輪在裝配后實(shí)際齒廓線空間位置與理論位置不符，提出了一種輪齒齒廓線的修正方法。李鵬超等[6]通過(guò)線性法調(diào)整柔輪齒的軸向齒高，使柔輪齒與剛輪齒不發(fā)生嚙合干涉。但在實(shí)際工作工程中，由于柔輪齒圈的存在，使得柔輪沿軸向的徑向變形量呈非線性關(guān)系，導(dǎo)致線性法設(shè)計(jì)修形存在一定偏差，不滿足實(shí)際的傳動(dòng)要求。

本文以無(wú)公切線雙圓弧齒廓為研究對(duì)象，采用改進(jìn)運(yùn)動(dòng)學(xué)法求解出共軛剛輪齒廓。為消除柔輪在不同徑向變形量下與剛輪嚙合產(chǎn)生的干涉，基于有限元徑向變形量結(jié)合輪齒修形理論，對(duì)柔輪齒廓進(jìn)行三維修形，并對(duì)比分析了線性修形法與有限元修形法修形后的柔輪等效應(yīng)力分布情況。

1 雙圓弧諧波傳動(dòng)共軛齒廓設(shè)計(jì)

1.1 柔輪輪齒雙圓弧齒廓設(shè)計(jì)

諧波傳動(dòng)中，雙圓弧齒廓具有更大的包絡(luò)存在區(qū)間。本文研究采用無(wú)公切線雙圓弧齒廓，相較于公切線雙圓弧齒廓，其共軛存在區(qū)間較大且設(shè)計(jì)求解過(guò)程更加簡(jiǎn)單，便于加工生產(chǎn)，其基本齒形如圖1 所示，圖1 中各參數(shù)的意義見(jiàn)表1。

表1 雙圓弧齒廓參數(shù)及意義

圖1 無(wú)公切線雙圓弧齒廓坐標(biāo)系

將柔輪輪齒的對(duì)稱(chēng)軸定義為Y1軸，將Y1軸與柔輪中性層曲線的交點(diǎn)定義為坐標(biāo)原點(diǎn)O1，過(guò)原點(diǎn)O1作中性面的切線，切線定義為X1軸，建立柔輪輪齒局部動(dòng)坐標(biāo)系{O1,X1,Y1}。

雙圓弧齒廓由凸圓弧AB、凹圓弧BC以及過(guò)渡圓弧CD組成。以從齒頂A處起的弧長(zhǎng)s為參數(shù)，r為右側(cè)齒廓的失徑，n為其對(duì)應(yīng)的法向量。

凸圓弧齒廓AB段為

凹圓弧齒廓BC段為

將分段函數(shù)的弧長(zhǎng)s離散為n個(gè)點(diǎn)，采用改進(jìn)運(yùn)動(dòng)學(xué)法[7-8]求解不同弧長(zhǎng)s對(duì)應(yīng)的共軛轉(zhuǎn)角φl(shuí)i（i=1,2,···,n），通過(guò)諧波傳動(dòng)精確轉(zhuǎn)角關(guān)系將式（1）與式（3）進(jìn)行坐標(biāo)變換，再將φl(shuí)i代入到變換后的公式中即可以求得離散的剛輪齒廓數(shù)值坐標(biāo)點(diǎn)集合，最后通過(guò)最小二乘法擬合得到特定的剛輪實(shí)際齒廓。

1.2 雙圓弧柔輪齒廓修形原理

假設(shè)柔輪內(nèi)壁與波發(fā)生器的接觸狀態(tài)為線接觸，在波發(fā)生器裝入柔輪后，柔輪輪齒在柔輪筒底和波發(fā)生器的影響下會(huì)產(chǎn)生彈性變形和法向傾角，柔輪齒在波發(fā)生器長(zhǎng)軸附近會(huì)向外張開(kāi)一定角度，在短軸處會(huì)向內(nèi)收縮一定角度，柔輪齒沿軸向方向發(fā)生不同程度的徑向變形。變形示意圖如圖2 所示。

圖2 柔輪變狀態(tài)圖

根據(jù)理論線性變形原理[9]可以求得在柔輪軸向上距筒底為Zg的任意截面的理論徑向變形量為

式中：ω0為柔輪主截面最大徑向變形量；Z0為主截面到柔輪筒底的距離。

在實(shí)際工作中，由于柔輪內(nèi)壁與波發(fā)生器接觸狀態(tài)復(fù)雜以及柔輪齒圈的影響，導(dǎo)致線性計(jì)算的徑向變形量不滿足使用要求，故需要通過(guò)有限元法計(jì)算柔輪長(zhǎng)短軸實(shí)際徑向變形量。

通過(guò)改進(jìn)運(yùn)動(dòng)學(xué)嚙合理論與Matlab 參數(shù)化編程，在不改變?nèi)我饨孛纨X廓形狀的條件下，合理調(diào)整柔輪齒根壁厚，使柔輪齒廓在任意轉(zhuǎn)角下都位于剛輪齒槽內(nèi)，從而實(shí)現(xiàn)無(wú)干涉嚙合。線性法修形是根據(jù)柔輪壁厚的調(diào)整量與齒寬的線性關(guān)系，確定其余各截面的壁厚調(diào)整量的。有限元法修形是根據(jù)各截面的徑向變形量與輪齒間的嚙合狀態(tài)來(lái)確定壁厚調(diào)整量的。

線性法修形后與有限元法修形后的柔輪齒形如圖3 所示，圖中dq與dh分別為前端修形量與后端修形量，虛線為修行前的柔輪齒，實(shí)線為修形后的柔輪齒。

圖3 柔輪輪齒軸向修形圖

2 諧波傳動(dòng)的有限元分析

2.1 共軛空間齒廓計(jì)算實(shí)例

以模數(shù)m=0.4、柔輪齒數(shù)z1=160、剛輪齒數(shù)z2=162、徑向變形量ω0=0.42、剛輪固定、余弦波發(fā)生器作為輸入、柔輪作為輸出、傳動(dòng)比i=80 的無(wú)公切線雙圓弧諧波減速器為研究對(duì)象，采用改進(jìn)運(yùn)動(dòng)學(xué)法計(jì)算出柔輪齒廓和共軛剛輪齒廓參數(shù)，結(jié)果見(jiàn)表2 和表3。

表2 柔輪齒廓參數(shù)

表3 剛輪齒廓參數(shù)

2.2 仿真實(shí)體建模

在SolidWorks 中建立諧波減速器三維模型，此時(shí)柔輪輪齒在軸向?yàn)槠烬X，未修形。諧波減速器的波發(fā)生器由凸輪和柔性軸承組成，在有限元分析中可將凸輪與軸承進(jìn)行一體化建模[10]。建模過(guò)程中從短軸處將凸輪對(duì)稱(chēng)分為兩部分，并施加位移約束，約束方向與凸輪長(zhǎng)軸方向一致，約束大小為柔輪徑向變形量0.42 mm，這樣做可以避免裝配過(guò)程中柔輪與波發(fā)生器產(chǎn)生干涉。將模型導(dǎo)入到Ansys Workbench 進(jìn)行分析，具體分析過(guò)程如下。

（1）材料選取。諧波減速器中柔輪材料選用30CrMnSiA（密度為7 750 kg/m3，彈性模量為204 GPa，泊松比為0.29）；剛輪和余弦凸輪材料選用45#優(yōu)質(zhì)碳素結(jié)構(gòu)鋼（密度為7 850 kg/m3，彈性模量為210 GPa，泊松比為0.269）。

（2）接觸設(shè)計(jì)與網(wǎng)格劃分。波發(fā)生器與柔輪，柔輪齒與剛輪齒間的接觸方式為面與面的摩擦接觸，波發(fā)生器外表面、剛輪輪齒表面為目標(biāo)面（Target），柔輪齒圈內(nèi)表面與輪齒表面為接觸面（Contact），摩擦系數(shù)設(shè)定為0.2，法向剛度因數(shù)設(shè)定為0.1，采用八節(jié)點(diǎn)六面體單元對(duì)模型進(jìn)行網(wǎng)格劃分[11]。

（3）設(shè)置約束及載荷。有限元分析步驟設(shè)定為3。第1 步對(duì)分割后的波發(fā)生器施加位移約束，使波發(fā)生器恢復(fù)處理前的形狀，模擬波發(fā)生器的裝配；第2 步對(duì)剛輪設(shè)置位移約束，使剛輪從后端裝入柔輪，剛輪齒廓與柔輪齒廓接觸；第3 步釋放柔輪繞軸向的旋轉(zhuǎn)約束并對(duì)柔輪輸出端施加額定扭矩T，T=60 N·m。

2.3 柔輪變形分析

通過(guò)Ansys Workbench 進(jìn)行三維實(shí)體有限元分析，得到柔輪徑向位移分布云圖，如圖4 所示。

圖4 柔輪徑向位移分布云圖

提取圖4 中截面1（前端截面）、截面2（主截面）、截面3（后端截面）中性層最大徑向變形量見(jiàn)表4。表中理論徑向變形量由式（5）計(jì)算而得。

表4 柔輪各截面中性層最大徑向變形量 mm

由表4 可知，柔輪中性層長(zhǎng)短軸實(shí)際徑向變形量與理論值相差較大，前端截面和主截面柔輪最大徑向變形量小于理論值，而在后端截面柔輪最大徑向變形量大于理論值，故線性修形不滿足實(shí)際要求，需采用有限元法對(duì)柔輪齒廓進(jìn)行修形。

將有限元仿真結(jié)果帶入Matlab 仿真程序中，可得到單個(gè)柔輪輪齒相對(duì)固定剛輪齒槽的運(yùn)動(dòng)軌跡，如圖5 所示，柔輪齒與剛輪齒在前端截面整體產(chǎn)生干涉，在主截面和后端截面柔輪齒與剛輪齒在嚙合區(qū)產(chǎn)生干涉。

圖5 未修形時(shí)柔輪齒運(yùn)動(dòng)軌跡

3 柔輪有限元法齒廓修形

柔輪有限元法修形的具體方法是利用有限元分析得到的柔輪輪齒軸向各截面的徑向變形量結(jié)合Matlab 參數(shù)化編程，觀察單個(gè)柔輪齒與剛輪齒的運(yùn)動(dòng)軌跡，根據(jù)輪齒間的干涉情況，調(diào)整各截面柔輪齒根壁厚，使剛輪齒與柔輪齒無(wú)干涉嚙合。

圖6 所示為線性修形法與有限元修形法的修形量對(duì)比圖，對(duì)有限元法的修形量進(jìn)行多項(xiàng)式擬合得到修形量方程：

圖6 柔輪齒修形量對(duì)比圖

式中：h為壁厚修形量，b為各截面距后端截面3的距離。修形后各截面運(yùn)動(dòng)軌跡如圖7 所示，有限元法修形后各截面柔輪齒與剛輪齒均可實(shí)現(xiàn)無(wú)干涉嚙合。

圖7 修形后柔輪各截面嚙合軌跡

通過(guò)Ansys Workbench 分別對(duì)未修形、線性法修形、有限元法修形后的仿真模型進(jìn)行有限元分析，柔輪在不同修形法下的應(yīng)力云圖如圖8 所示。由圖8可知，未修形時(shí)柔輪最大等效應(yīng)力為592.6 MPa，柔輪齒前端和后端產(chǎn)生了應(yīng)力集中的現(xiàn)象。線性法修形后柔輪最大等效應(yīng)力為404.09 MPa，有限元法修形后柔輪最大等效應(yīng)力為344.89 MPa，有限元法修形后柔輪最大等效應(yīng)力主要集中在柔輪齒中間截面附近，且應(yīng)力沿圓周方向分布更加均勻。

圖8 不同修形法修形后柔輪等效應(yīng)力云圖

4 結(jié)語(yǔ)

本文基于裝配變形與改進(jìn)運(yùn)動(dòng)學(xué)理論，對(duì)諧波傳動(dòng)無(wú)公切線雙圓弧齒廓進(jìn)行了設(shè)計(jì)和修形，并對(duì)實(shí)體模型進(jìn)行有限元仿真，分析了柔輪的變形及應(yīng)力分布情況，得到以下結(jié)論：

（1）柔輪在實(shí)際工作過(guò)程中柔輪徑向變形量與理論變形量有較大的差異，在齒圈前端柔輪最大徑向變形量小于理論值，而在齒圈后端柔輪最大徑向變形量大于理論值。

（2）通過(guò)有限元法修形后柔輪齒廓與剛輪齒廓在各截面均可實(shí)現(xiàn)無(wú)干涉嚙合。

（3）未修形時(shí)柔輪齒圈處最大等效應(yīng)力集中在柔輪齒圈前端位置、齒圈與筒體過(guò)度位置。有限元法修形后最大等效應(yīng)力集中在齒圈中間截面附近且齒圈圓周方向的應(yīng)力分布更加均勻。

（4）有限元法修形后柔輪的最大等效應(yīng)力比線性法修形后減小14.65%，比未修形減小41.79%，說(shuō)明有限元法修形可以有效解決線性法修形不充分以及未修形時(shí)輪齒間的嚙合干涉引起應(yīng)力集中的問(wèn)題。