Open Access
Issue
JNWPU
Volume 42, Number 4, August 2024
Page(s) 634 - 642
DOI https://doi.org/10.1051/jnwpu/20244240634
Published online 08 October 2024

© 2024 Journal of Northwestern Polytechnical University. All rights reserved.

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

行星滚柱丝杠是一种可将旋转运动转化为直线运动的精密机械传动机构[1], 具有转速高、承载能力强、精度高等特点[2]。在航空、航天、精密机床、食品包装、能源化工、冶金及特种装备等领域有着广泛的应用前景[3]。面向机械装备高速度、高精度发展需求, 对行星滚柱丝杠的动态承载性能要求也愈来愈高, 然而行星滚柱丝杠在高速工作状态下, 不仅承受轴向载荷, 还会承受较大的惯性力, 这将导致滚柱与丝杠及螺母的接触力及载荷分布发生显著变化。因此, 求解高速工作状态下行星滚柱丝杠接触力及建立动载荷分布模型至关重要。

近年来, 针对行星滚柱丝杠螺纹牙载荷分布的理论研究已取得一定成果。杨家军等[4]基于赫兹弹性接触理论, 考虑了赫兹接触变形及轴向变形, 建立了行星滚柱丝杠轴向静刚度模型, 通过递推计算得到了行星滚柱丝杠螺纹牙载荷分布曲线。Jan等[5]考虑螺纹牙接触变形及丝杠、滚柱与螺母的螺纹轴段变形, 建立了螺纹牙载荷分布计算模型, 并建立有限元验证模型, 将数值解与有限元解对比验证, 但二者吻合度较差。Jones等[67]基于直接刚度法建立了行星滚柱丝杠轴向刚度理论模型, 为计算行星滚柱丝杠螺纹副载荷分布提供了一种有效的理论方法。Folly等[8]等基于滚柱间载荷均匀分配以及材料处于弹性阶段的假设, 考虑4种边界条件及螺母相对于丝杠不同位置, 建立了行星滚柱丝杠载荷分布模型。刘柱等[910]建立了考虑装配误差和滚柱弯曲的行星滚柱丝杠载荷分布模型, 通过求解有限元模型分析验证了计算模型的正确性。郑正鼎等[11]基于空间啮合理论, 分析了差动式行星滚柱丝杠载荷分布特点。王敏[12]建立了反向式行星滚柱丝杠静载荷分布模型, 并对静载荷分布进行了试验验证, 分析了螺纹牙变形量及螺纹牙啮合段静载刚度值。Du等[13]基于赫兹接触理论和螺纹累积变形计算所有螺纹的接触载荷, 建立了包含径向载荷和加工误差的行星滚柱丝杠载荷分布模型。刘荣荣等[14]根据赫兹接触及变形协调理论建立凸凹接触式行星滚柱丝杠载荷分布模型, 分析了牙侧角、螺距等参数对啮合点位置、轴向间隙和载荷分布的影响。

综上所述, 目前考虑惯性力对螺纹副载荷分布影响的研究较少。故本文基于螺旋曲面啮合原理及变形协调理论, 建立考虑惯性力的行星滚柱丝杠螺纹副动载荷分布模型。

1 动载荷分布模型

1.1 静态接触力

螺纹牙静态接触力计算方法基于以下假设:

1) 多个滚柱间载荷均匀分布;

2) 仅考虑弹性变形的影响;

3) 不考虑螺纹牙啮合间隙产生冲击载荷影响。

1.1.1 滚柱与丝杠侧静态接触力

基于行星滚柱丝杠螺旋曲面啮合原理[15], 在如图 1所示的局部坐标系oPR-xPRyPRzPR中对滚柱与丝杠接触侧进行受力分析。

图 1中, 坐标系XOY为整体坐标系, ϕSR为滚柱侧啮合偏角, ϕSR为丝杠侧啮合偏角, rSR为滚柱侧啮合半径, r′SR为丝杠侧啮合半径, oR为滚柱轴线在局部坐标系上的投影, oSR为丝杠与滚柱的接触点, FSRs为滚柱与丝杠的静态接触力。

根据丝杠和滚柱螺旋曲面在接触点处的法向量分析结果[15], 接触力FSRs可表示为:

式中, FSRs为接触力的幅值, 可表示为

式中:λSRβSR分别为丝杠与滚柱接触点处的螺旋升角和牙侧角;ϕSR为滚柱在丝杠侧的啮合偏角;FNZ为螺母上的负载;nroller为滚柱个数。

式中:LR为滚柱的导程;rSR为滚柱在丝杠侧的啮合半径;rR为滚柱的名义半径;rTR为滚柱的牙型轮廓半径;uTR为滚柱牙型轮廓圆心在螺纹牙截面uRvRwRuR坐标值, 螺纹牙截面和参数定义见文献[15]。

thumbnail 图1

滚柱-丝杠侧接触力

1.1.2 滚柱与螺母侧静态接触力

同理, 在滚柱与螺母侧进行受力分析如图 2所示。其中, 坐标系XOY为整体坐标系, ϕNR为滚柱侧啮合偏角, ϕNR为螺母侧啮合偏角, rNR为滚柱侧啮合半径, r′NR为螺母侧啮合半径, oR为滚柱轴线在局部坐标系上的投影, oNR为滚柱与螺母在局部坐标系oPR-xPRyPRzPR中的接触点, FNRs为滚柱与螺母的静态接触力。

根据丝杠和滚柱螺旋曲面在接触点处的法向量分析结果[15], 滚柱-螺母侧接触力FNRs可表示为

由于螺母在触点处的螺旋升角λNR和牙侧角βNR等于滚柱的螺旋升角λR和牙侧角βR, 则螺母与滚柱的啮合偏角ϕNR=0, 故FNRs可表示为

式中:FNRs为滚柱-螺母侧静态接触力幅值;λRβR分别为螺母的螺旋升角和牙侧角。可分别表示为

thumbnail 图2

滚柱-螺母侧接触力

1.2 考虑惯性力的动态接触力

当行星滚柱丝杠承受轴向载荷作用时, 滚柱与丝杠及螺母之间发生接触变形。滚柱与丝杠及螺母的法向接触力的接触变形为

式中:δXR为滚柱与丝杠或螺母的接触变形;δ*为Hertz接触变形系数;FXRs为滚柱与丝杠或螺母的静态接触力;μX为丝杠或螺母的泊松比;EX为丝杠或螺母的弹性模量;∑ρ为滚柱与丝杠或螺母接触曲面的曲率和函数。

图 3所示, 行星滚柱丝杠在承载后, 滚柱与丝杠的接触点从oSR移动至o′SR, 滚柱与螺母的接触点从oNR移动至o′NR, 达到平衡后, 滚柱与丝杠及螺母的接触角为αSRαNR, 接触变形的轴向分量为δSRaδNRa

根据变形协调几何关系有

式中:αSR0αNR0分别为滚柱与丝杠及螺母侧接触角初值;rS, rNrR为丝杠、螺母及滚柱的名义半径。

根据(12)~(13)式可得, 承载后接触角的计算公式为

行星滚柱丝杠在高速运转时, 滚柱绕丝杠公转将产生不可忽略的惯性力。考虑惯性力对滚柱与丝杠及螺母侧的受力影响, 如图 4所示。

将滚柱作为研究对象, 对其进行受力分析, 可得其轴向和径向力平衡方程

式中:Fc为惯性力;ωP为滚柱公转角速度;r为丝杠和滚柱的中径和;FSRFNR分别为滚柱与丝杠及螺母侧的动态接触力。

thumbnail 图3

变形协调几何关系

thumbnail 图4

行星滚柱丝杠螺纹受力示意图

1.3 动载荷分布

根据行星滚柱丝杠结构特点与承载原理, 其变形主要包含螺纹轴段变形、螺纹牙变形及螺纹接触点变形。螺纹轴段将承受压应力或拉应力, 进而发生压缩或拉伸变形; 螺纹牙承受法向接触载荷后, 将发生弯曲、剪切变形以及由于螺纹基体径向变形导致的变形; 螺纹牙接触点则在接触载荷的作用下发生法向接触变形。对应于3类变形, 分别定义其刚度为轴段刚度、螺纹牙刚度及螺纹牙接触刚度, 这3类刚度是建立行星滚柱丝杠螺纹牙载荷分布计算模型的基础。

根据螺纹牙变形协调关系[16], 以滚柱与螺母侧为例, 可得第i个螺纹闭环内螺母的变形总和为

式中:i表示螺纹牙序号;下标B代表轴段、T代表螺纹牙、C代表接触点Hertz变形, 不同下标组合表示Δl是由对应部分的变形引起的。

根据变形协调关系有

由于PN=PR, 则有

将轴段刚度、螺纹牙刚度及接触刚度代入(21)式得

同理, 滚柱与丝杠侧的变形协调方程为

滚柱与丝杠或螺母侧的各螺纹牙接触力之和等于滚柱与丝杠或螺母侧的整体动态接触力

2 模型验证

2.1 算例

表 1为一组行星滚柱丝杠结构参数。在螺旋传动机构中, 直线进给速度达60 m/min以上, 为高速[16]。结合表 1中丝杠直径, 可得丝杠的转速为83 rad/s。表 2为行星滚柱丝杠工况及滚柱结构参数。

表 3所示, 当丝杠转速ωS=0时, 滚柱-丝杠侧静态接触力大于滚柱-螺母侧, 与文献[17]计算结果对比, 本文计算结果更精确。当丝杠转速ωS=83 rad/s时, 滚柱-丝杠侧动态接触力小于滚柱-螺母侧, 此结果与静态接触力的结果相反, 更贴近行星滚柱丝杠实际工作中的载荷分布, 故本文计算方法相比于文献[17]更完善。

当丝杠转速设置为83 rad/s时, 滚柱的公转速度为31.12 rad/s, 螺纹序号如图 4所示。未考虑惯性力的行星滚柱丝杠螺纹牙静载荷分布结果如图 5所示, 根据本文所述方法求得滚柱与丝杠及螺母侧动载荷分布的结果如图 6所示。

图 5可知, 未考虑惯性力时, 滚柱-丝杠侧各螺纹牙接触力随着螺纹牙序号的增大而减小, 滚柱-螺母侧各螺纹牙的接触力随着螺纹牙的增大而增大。由图 6可知, 当丝杠转速为83 rad/s, 即滚柱公转速度为31.12 rad/s时, 滚柱两侧的动载荷分布曲线趋势与未考虑惯性力时的载荷分布趋势相似, 但滚柱-丝杠侧的动态接触力小于滚柱-螺母侧的动态接触载荷; 产生这一结果的原因是滚柱绕丝杠公转产生惯性力使得滚柱靠近螺母远离丝杠, 从而使得滚柱-螺母侧的变形增大, 滚柱-丝杠侧的变形减小, 最终导致滚柱-丝杠侧的动态接触力小于滚柱-螺母侧。

表1

行星滚柱丝杠主要结构参数

表2

行星滚柱丝杠工况及滚柱结构参数

表3

滚柱-丝杠及滚柱-螺母侧接触力计算结果

thumbnail 图5

行星滚柱丝杠螺纹副静载荷分布

thumbnail 图6

行星滚柱丝杠螺纹副动载荷分布

2.2 有限元仿真及模型验证

建立行星滚柱丝杠螺纹牙动载荷分布有限元模型, 分析其在丝杠转速为83 rad/s时滚柱-丝杠侧及滚柱-螺母侧的螺纹牙动载荷分布, 并与理论模型计算的结果进行对比。

2.2.1 模型前处理

根据表 1行星滚柱丝杠结构参数建立三维模型。为提高计算效率, 对三维模型进行简化处理[18]: 仅保留滚柱、丝杠、螺母、保持架及内齿圈等主要零件, 将原有的10个滚柱简化为均布的3个, 每个滚柱保留20个螺纹牙。单元类型为Tetrahedrons四面体网格, 在螺纹副接触区域进行网格密化, 单元数955 519, 节点数1 555 868。零件材料均选择GCr15, 密度ρ=7 810 kg/m3, 弹性模量E=212 GPa, 泊松比μ=0.29。

根据行星滚柱丝杠各零件的运动方式分别添加对应的运动副。定义滚柱和丝杠、螺母螺纹副接触,滚柱齿和内齿圈齿轮副接触,滚柱轴端和保持架接触。丝杠施加转速ωS=83 rad/s, 螺母施加与运动方向相反的轴向负载FNZ=1 500 N。惯性力是转速的函数, 因此, 在滚柱上定义径向惯性力。行星滚柱丝杠有限元模型如图 7所示。

thumbnail 图7

行星滚柱丝杠有限元接触模型

2.2.2 有限元结果分析

滚柱与丝杠及螺母之间的动态接触力曲线及动载荷分布曲线如图 8所示。

图 8可知, 当丝杠转速ωS=83 rad/s时, 滚柱与丝杠之间动态接触力FSR的均值为683 N, 滚柱与螺母之间动态接触力FNR的均值为734 N, 与理论值的相对误差不足5%。此结果验证了行星滚柱丝杠螺纹牙动载荷分布计算方法的正确性。

图 9可知, 滚柱与丝杠及螺母侧的动载荷分布规律与理论模型计算的结果相同: 滚柱与丝杠侧的动态接触力小于滚柱与螺母侧的动态接触力; 滚柱与丝杠侧各螺纹牙动载荷随螺纹牙序号增大而减小, 滚柱与螺母侧各螺纹牙动载荷随螺纹牙序号的增大而增大; 滚柱与丝杠侧的动载荷分布不均程度大于滚柱与螺母侧。有限元模型存在微量的径向位移, 导致结果出现轻微波动, 而理论模型忽略了径向的轻微变形。

thumbnail 图8

行星滚柱丝杠动态接触力数值解

thumbnail 图9

行星滚柱丝杠螺纹副动载荷分布结果对比图

3 参数影响规律

3.1 丝杠转速

将丝杠转速设置为26~133 rad/s, 则滚柱的公转角速度为10~50 rad/s, 负载为5 000 N, 则滚柱与丝杠侧和螺母侧的螺纹牙动载荷分布分别如图 10~11所示。

图 10~11可知, 随着丝杠转速的增大, 滚柱与丝杠侧的动态接触力增大, 滚柱与螺母侧的动态接触力减小; 滚柱与丝杠侧各螺纹牙上的载荷随螺纹牙序号增大而减小, 滚柱与螺母侧各螺纹牙上的载荷随螺纹牙序号的增大而增大。产生此现象的原因是丝杠转速较小时, 滚柱绕丝杠公转的惯性力较小, 而当丝杠转速大于某一临界值时, 惯性力则会使滚柱向螺母侧偏移, 导致螺母侧的接触力增大, 丝杠侧的接触力减小。由此可见, 丝杠转速对行星滚柱丝杠螺纹牙动载荷分布的影响较大。为提高行星滚柱丝杠承载特性, 确定丝杠转速时, 不仅要满足运动学要求, 还需考虑丝杠转速对动载荷分布的影响。

thumbnail 图10

丝杠转速对滚柱-丝杠侧动载荷分布的影响

thumbnail 图11

丝杠转速对滚柱-螺母侧动载荷分布的影响

3.2 轴向载荷

设置螺母上的负载为0~6 000 N, 丝杠转速为83 rad/s, 求解滚柱与丝杠及螺母侧的载荷分布, 结果如图 12~13所示。

图 12~13可知, 随着负载的增大, 滚柱与丝杠及螺母侧的动态接触力均随之增大; 滚柱与丝杠侧各螺纹牙上的载荷随螺纹牙序号增大而减小, 滚柱与螺母侧各螺纹牙上的载荷随螺纹牙序号的增大而增大; 滚柱与丝杠侧的载荷分布不均程度大于滚柱与螺母侧, 因此, 改善载荷不均现象应着重关注滚柱与丝杠侧。

thumbnail 图12

轴向载荷对滚柱-丝杠侧动载荷分布的影响

thumbnail 图13

轴向载荷对滚柱-螺母侧动载荷分布的影响

3.3 螺距

将螺距作为控制变量, 分析其对动载荷分布的影响。设置螺距P分别为2, 4和6, 丝杠转速为83 rad/s, 负载为5 000 N, 求解滚柱与丝杠及螺母侧的动载荷分布, 结果如图 14~15所示。

图 14~15可知, 螺距越大, 滚柱与丝杠及螺母侧的动载荷分布曲线均越陡, 表明滚柱与丝杠及螺母侧的动载荷分布不均程度随螺距的增大而增大。可见, 为提高行星滚柱丝杠的承载能力及运动稳定性, 应在满足传动效率的基础上尽可能选择较小的螺距。

thumbnail 图14

螺距对滚柱-丝杠侧动载荷分布的影响

thumbnail 图15

螺距对滚柱-螺母侧动载荷分布的影响

3.4 牙侧角

设置牙侧角分别为40°, 45°和50°, 丝杠转速为83 rad/s, 负载为5 000 N, 分析牙侧角对滚柱与丝杠及螺母侧动载荷分布影响, 结果如图 16~17所示。

图 16~17可知, 随着牙侧角的增大, 滚柱与丝杠及螺母侧的螺纹牙动载荷分布曲线均变缓, 表明滚柱与丝杠及螺母侧螺纹牙动载荷分布不均程度随着牙侧角的增大而减小。可见, 在选择牙侧角时不仅要以传动效率为设计原则, 还需考虑牙侧角对动载荷分布的影响。

thumbnail 图16

牙侧角对滚柱-丝杠侧动载荷分布的影响

thumbnail 图17

牙侧角对滚柱-螺母侧动载荷分布的影响

4 结论

1) 考虑滚柱高速运转产生的惯性力作用, 建立了行星滚柱丝杠螺纹副动载荷分布计算模型, 与有限元接触模型对比, 验证了计算模型的正确性;

2) 惯性力使滚柱产生靠近螺母、远离丝杠的运动趋势, 导致滚柱-丝杠侧接触力减小, 滚柱-螺母侧接触力增大;

3) 动载荷分布不均匀程度随着丝杠转速及轴向载荷的增大而增大, 且轴向载荷对于滚柱-丝杠侧动载荷分布影响更大;

4) 动载荷分布不均匀程度随着牙侧角及螺距的增大而减小。

References

  1. STRANDGREN C B. Roller Screw: US, 3182522[P]. 1965-05-11 [Google Scholar]
  2. LIU Geng, FU Xiaojun, MA Shangjun, et al. Review on meshing theory of the planetary roller screw mechanism[J]. Mechanical Drive, 2015, 39(7): 1–5 (in Chinese) [Google Scholar]
  3. CLAEYSSEN F, JANKER P, LELETTY R, et al. New actuators for aircraft, space and military applications[C]//12th International Conference on New Actuators, Bremen, Germany, 2010: 324–330 [Google Scholar]
  4. YANG Jiajun, WEI Zhenxing, ZHU Jisheng. Calculation of load distribution of planetary roller screws and static rigidity[J]. Journal of Huazhong University of Science and Technology, 2011, 39(4): 1–4 (in Chinese) [Google Scholar]
  5. JAN R, LISOWSKI F. The computational model of the load distribution between elements in planetary roller screw[C]//9th International Conference on Fracture & Strength of Solids, Jeju, Korea, 2013: 1–6 [Google Scholar]
  6. JONES M H, VELINSKY S A. Stiffness of the roller screw mechanism by the direct method[J]. Mechanics Based Design of Structures and Machines, 2014, 42: 17–34 [Article] [CrossRef] [Google Scholar]
  7. JONES M. H, VELINSKY S A. Contact kinematics in the roller screw Mechanism[J]. Mechanical Design, 2013, 135(1): 1–10 [Google Scholar]
  8. FOLLY A, ALAIN D, MICHEL C. Static load distribution and axial stiffness in a planetary roller screw mechanism[J]. Journal of Mechanical Design, 2016, 138(1): 012301 [Article] [CrossRef] [Google Scholar]
  9. LIU Zhu, FAN Yuanxun. Analysis and research on the load of planetary roller screw pairs considering assembly errors[J]. Mechanical Manufacturing and Automation, 2020, 49(6): 120–123 (in Chinese) [Google Scholar]
  10. LIU Zhu, FAN Yuanxun. Research on load distribution of planetary roller screw pairs considering roller bending[J]. Modular Machine Tools and Automated Processing Technology, 20207): 60–64 (in Chinese) [Google Scholar]
  11. ZHENG Zhengding, CHEN Bingkui, DU Xing, et al. Analysis of load-bearing characteristics of differential planetary roller screw[J]. Journal of Chongqing University, 2020, 43(12): 23–32 (in Chinese) [Google Scholar]
  12. WANG Min. Research on static load carrying capacity of reverse planetary roller screw[D]. Harbin: Harbin Institute of Technology, 2021 (in Chinese) [Google Scholar]
  13. DU X, CHEN B, ZHENG Z. Investigation on mechanical behavior of planetary roller screw mechanism with the effects of external loads and machining errors[J]. Tribology International, 2021, 154: 106689 [Article] [CrossRef] [Google Scholar]
  14. LIU Rongrong, CHEN Yonghong, DU Xing, et al. The meshing load bearing characteristics of convex concave contact planetary roller screw[J]. Journal of Harbin Institute of Technology, 2022, 54(1): 49–57 (in Chinese) [Google Scholar]
  15. FU Xiaojun, LIU Geng, MA Shangjun, et al. Research on the meshing mechanism of planetary roller screw pairs with helical surfaces[J]. Journal of Mechanical Engineering, 2016, 52(3): 26–33 (in Chinese) [Google Scholar]
  16. HUANG Yuquan, YU Zhongzhi. The development history and future trends of ball screw pairs in China[J]. Mechanical Drive, 200410): 60–62 (in Chinese) [Google Scholar]
  17. WENJIE Z, GENG L, RUITING T, et al. Load distribution of planetary roller screw mechanism and its improvement approach[J]. Journal of Mechanical Engineering Science, 2016, 230(18): 3304–3318 [Article] [CrossRef] [Google Scholar]
  18. MA Tenglong, MA Shangjun, WU Linping, et al. Dynamic contact load characteristics of synchronous meshing between planetary roller screw thread pair and gear pair[J]. Journal of Xi'an Jiaotong University, 202210): 1–12 (in Chinese) [Google Scholar]

All Tables

表1

行星滚柱丝杠主要结构参数

表2

行星滚柱丝杠工况及滚柱结构参数

表3

滚柱-丝杠及滚柱-螺母侧接触力计算结果

All Figures

thumbnail 图1

滚柱-丝杠侧接触力

In the text
thumbnail 图2

滚柱-螺母侧接触力

In the text
thumbnail 图3

变形协调几何关系

In the text
thumbnail 图4

行星滚柱丝杠螺纹受力示意图

In the text
thumbnail 图5

行星滚柱丝杠螺纹副静载荷分布

In the text
thumbnail 图6

行星滚柱丝杠螺纹副动载荷分布

In the text
thumbnail 图7

行星滚柱丝杠有限元接触模型

In the text
thumbnail 图8

行星滚柱丝杠动态接触力数值解

In the text
thumbnail 图9

行星滚柱丝杠螺纹副动载荷分布结果对比图

In the text
thumbnail 图10

丝杠转速对滚柱-丝杠侧动载荷分布的影响

In the text
thumbnail 图11

丝杠转速对滚柱-螺母侧动载荷分布的影响

In the text
thumbnail 图12

轴向载荷对滚柱-丝杠侧动载荷分布的影响

In the text
thumbnail 图13

轴向载荷对滚柱-螺母侧动载荷分布的影响

In the text
thumbnail 图14

螺距对滚柱-丝杠侧动载荷分布的影响

In the text
thumbnail 图15

螺距对滚柱-螺母侧动载荷分布的影响

In the text
thumbnail 图16

牙侧角对滚柱-丝杠侧动载荷分布的影响

In the text
thumbnail 图17

牙侧角对滚柱-螺母侧动载荷分布的影响

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.