考慮主平面夾角的行星滾柱絲杠副接觸特性

2021-12-23 11:19:22鄭繼貴張兆晶何理論

導彈與航天運載技術 2021年6期

陳中，鄭繼貴，侍威，張兆晶，何理論

（北京精密機電控制設備研究所，北京，100076）

0 引言

行星滾柱絲杠副（Planetary Roller Screw Mechanism，PRSM）作為將旋轉運動轉換為直線運動的重要傳動機構，具有高承載、高精度、高壽命的優異性能，已經在精密數控設備、光學測量儀器、航空航天、武器裝備等領域得到了廣泛的應用[1,2]。PRSM在傳動過程中通過螺紋曲面相互接觸來傳遞載荷，特別對于滾柱和絲杠接觸側，由于絲杠的螺旋升角和滾柱不一致，使得兩者之間的接觸更為復雜，而接觸區域的精確求解對于后續接觸特性分析具有重要意義。

目前，一些學者針對PRSM的接觸特性開展研究，Wang Shuiming等[3]通過建立接觸匹配方程對接觸點的位置進行精確求解，并研究了零件的輪廓誤差對接觸特性的影響規律。關于PRSM接觸區域的求解多采用等效球法來計算主曲率，劉艷強等[4]基于等效球法求解了PRSM的接觸參數，并研究了結構參數對接觸變形及應力分布的影響規律。然而，等效球法計算主曲率存在較大偏差，因而部分學者引入微分幾何原理來計算主曲率。Jones等[5]通過建立Frenet坐標系推導出PRSM各零件螺紋曲面方程，求解了接觸點的位置以及主曲率。Sebastian Sandu等[6,7]建立了螺紋法截面輪廓的螺旋曲面方程，并基于微分幾何原理對主曲率以及接觸橢圓的位置和形狀展開了研究。喬冠等[8]基于螺紋滾道法截面輪廓，建立了PRSM各零件螺紋曲面統一方程，求解了接觸點處的主曲率，同時引入主曲率和與主曲率差對接觸特性進行分析。綜上所述，利用微分幾何原理的方法可以更準確計算出主曲率，但現有研究在求解時都有所簡化，認為接觸體之間的主平面相互重合，這在簡化計算的同時也會帶來偏差。

本文以標準式PRSM為研究對象，結構如圖1所示。針對PRSM的傳動特點，基于微分幾何原理計算出主平面夾角，進而精確求解接觸橢圓的形狀，并與現有研究的簡化方法進行對比，同時研究了主平面夾角對接觸特性的影響規律。

圖1 PRSM結構Fig.1 Structure Drawing of PRSM

1 主曲率及主方向

滾柱絲杠側和滾柱螺母側的接觸可看做兩空間曲面之間的彈性接觸，而主曲率以及主方向的求解是PRSM接觸特性分析的基礎。以絲杠、滾柱和螺母的軸線為zi軸，xi軸經過螺紋起始點，分別建立各零件的坐標系，螺紋曲面上任意點空間坐標可表示為(rpi,θpi,zi)，下標i=s,r,n，分別表示絲杠、滾柱和螺母，如圖2所示。

圖2 零件坐標系Fig.2 Part Coordinate System

依據螺紋法截面輪廓可建立PRSM零件螺旋曲面統一方程[8]：

式中ξ可取-1和1，分別表示上曲面和下曲面；h(rpi)為各零件在法截面的輪廓方程；li和λi分別為零件的導程和螺旋升角。

螺旋曲面參數方程可表示為

零件輪廓曲面的單位法向量可表示為

基于微分幾何理論[9]，可依據零件曲面的參數方程求得曲面在任意一點的主曲率及主方向。螺旋曲面的第1類基本量E、F、G和第2類基本量L、M、N的計算公式如下：

式中fr和fθ為一階導數；frr，frθ和fθθ為二階導數。

曲面上任意一點的兩個主曲率可用κ1和κ2表示。

式中H為平均曲率；K為高斯曲率。

主方向用dr:dθ表示，滿足下列等式：

這是關于dr:dθ的二次方程，可解得主方向：

將主方向轉化為零件的笛卡爾坐標系中的矢量為

2 主平面夾角

主方向與接觸點處的法向量相互垂直，兩者共同組成主平面。由式（12）計算可得到各零件坐標系下接觸點處的主方向，通過建立全局坐標系，利用坐標變換關系可求得兩主平面夾角。由于主平面夾角不受主方向矢量的位置影響，因此只需對各零件在接觸點處的主方向進行旋轉變換即可。本文選取全局坐標系和絲杠坐標系重合，為確保螺紋能夠正確接觸，建立如圖3所示的坐標系。

圖3 正確接觸時坐標系示意（xOy平面）Fig.3 Schematic Diagram of Coordinate System in Correct Contact (xOy Plane)

以5頭螺紋的絲杠為例，其他頭數推導過程類似。假設絲杠側接觸點在理論中徑所在螺旋線上，選取絲杠和螺母的某一頭螺紋，當絲杠在角度為αs位置上與某一滾柱接觸時（αi為接觸點角度，即接觸點相對xi軸轉過的角度，用弧度制表示，i=s,rs,rn,n ，分別表示絲杠、滾柱在絲杠側、滾柱在螺母側以及螺母處的接觸點），絲杠側接觸點相對螺紋起始點高度：

式中p為螺距。

為保證滾柱與絲杠的準確接觸，滾柱與絲杠在該點處的高度應該相等，即：

αrs=5αs，此時滾柱在滾柱螺母側的接觸點角度為

同理，滾柱螺母側接觸點處應滿足：

計算得到螺母接觸點角度為

上述每個零件的曲面方程所建立的曲面坐標系是相互獨立的，且都是基于全局坐標系xyz（見圖4）。

圖4 曲面方程確定的零件坐標系（xOy平面）Fig.4 Part Coordinate System Determined by Surface Equation (xOy Plane)

因此，為達到圖3所示的接觸狀態，需要對此時的零件坐標系進行坐標變換。由圖2可知，全局坐標系變換到滾柱坐標系的旋轉角度為

則滾柱的旋轉變換矩陣為

全局坐標系變換到螺母坐標系的旋轉角度為

因此，螺母的旋轉變換矩陣為

對于滾柱絲杠側，由式（12），可以得到絲杠和滾柱的主方向，設為ts1、ts2、trs1和trs2，對應的單位法向量為ns和nrs，則全局坐標系下絲杠兩主平面可表示為[ts1;ns]和[ts2;ns]，滾柱兩主平面為[kr?trs1;kr?nrs]和[kr?ts2;kr?nrs]，利用matlab軟件，即可求得絲杠側兩接觸零件之間主平面的夾角ωs，同理可計算得到螺母側主平面夾角ωn。

3 Hertz接觸特性分析

依據Hertz接觸理論，兩點接觸的彈性體在受到外加載荷的情況下，接觸點處發生彈性變形變為橢圓區域。兩接觸體點接觸一般情況如圖5所示，其中R11表示物體1在主平面1內的曲率半徑，第1個下標表示接觸體代號，第2個下標表示主平面代號[10]。

圖5 兩接觸體點接觸一般情況Fig.5 General Situation of Point Contact between Two ContactBodies

接觸橢圓長短半軸計算公式分別為

式中Q為法向接觸載荷；E，μ分別為材料的彈性模量和泊松比；ma，mb分別為橢圓長短軸系數；為主曲率和。計算公式如下：

式中e為橢圓偏心率。

主曲率函數滿足下列等式：

式中k為橢圓短半軸和長半軸之比；K(e)，L(e)分別表示第1類和第2類完全橢圓積分。

依據式（22）~（26）即可求得橢圓偏心率e。

4 算例

本文采用PRSM的結構參數如表1所示。

表1 PRSM基本參數Tab.1 Basic Parameters of PMSM

依據表1中的基本參數，利用Matlab軟件繪制螺旋曲面如圖6所示。由圖6可看出在滾柱絲杠接觸側存在明顯干涉，導致在名義中徑處主平面[ts1;ns]和[kr?trs1;kr?nrs]的夾角與主平面[ts2;ns]和[kr?ts2;kr?nrs]的夾角不一致，因此取兩夾角的均值來近似表示滾柱絲杠側的主平面夾角。

ω=0的簡化計算與本文的計算結果對比如表2所示。利用matlab軟件，計算出接觸點角度αs?[0,2π]時滾柱絲杠側與滾柱螺母側主平面之間的夾角ωs和ωn，變化曲線如圖7所示，可以看出，滾柱在任意接觸位置絲杠側和螺母側主平面之間的夾角保持不變。

圖6 各零件曲面Fig.6 Contact Surface Diagram of Each Part

表2 兩種方法計算結果對比Tab.2 Comparison of Calculation Results of Two Methods

圖7 不同接觸位置主平面夾角Fig.7 The Principal Plane Angel at Different Contact Positions

螺距和牙型半角對絲杠側和螺母側主平面夾角的影響如圖8所示，星號標出數據為絲杠直徑ds=30 mm情況下Rollvis行星滾柱絲杠副RV系列產品所使用螺距的計算結果。可以看出，絲杠側和螺母側主平面夾角隨螺距的增大而增大，呈現出近似線性關系。并且在螺距較大時，絲杠側夾角增速放緩，使得兩者夾角的偏差增大。而隨螺紋牙型半角的增大，主平面夾角減小，且螺母側減小的速度更快。

橢圓偏心率e的取值范圍為[0,1]，e越大，橢圓越扁。圖9和圖10分別為絲杠側、螺母側螺距和牙型半角對偏心率的影響曲線。在考慮主平面夾角的情況下，接觸橢圓的偏心率并不是隨螺距和牙型半角線性變化的，而是波動改變的，但絲杠側橢圓偏心率相較于ω=0時普遍偏大，而螺母側則呈現出普遍偏小的趨勢，且兩種方法螺距隨偏心率的最大計算偏差能達到8.7%，牙型半角隨偏心率的最大計算偏差能達到7.2%。