###
工程科学与技术:2021,53(2):178-186
←前一篇   |   后一篇→
本文二维码信息
码上扫一扫!
高速滚珠丝杠副非赫兹接触应力计算
(1.北京工业大学 材料与制造学部 先进制造技术重点实验室,北京 100124;2.电火花加工技术北京市重点实验室,北京 100191;3.中国铁道科学研究院 标准计量研究所,北京 100081)
Non-Hertzian Contact Stress Calculation of High-speed Ball Screw Mechanism
(1.Beijing Key Lab. of Advanced Manufacturing Technol., Faculty of Materials and Manufacturing, Beijing Univ. of Technol., Beijing 100124, China;2.Beijing Key Lab. of Electro-Machining Technol., Beijing 100191, China;3.Inst. of Standard Metrology, Chinese Academy of Railway Sciences, Beijing 100081, China)
摘要
图/表
参考文献
相似文献
本文已被:浏览 1421次   下载 544
投稿时间:2020-05-27    修订日期:2020-07-01
中文摘要: 在滚珠丝杠副运动过程中,滚珠与滚道之间的摩擦和磨损导致滚珠丝杠副精度衰退,降低其精度保持性。针对高速滚珠丝杠副螺旋线滚道接触区域的非对称性,采用传统的赫兹接触理论对丝杠和螺母滚道接触几何形状和关系的计算存在较大误差,提出了一种基于最小余能原理对高速滚珠丝杠副滚珠与滚道之间的法向接触应力进行分析计算的方法。首先,基于Frenet-Serret坐标转换公式建立滚珠丝杠副坐标系,对滚珠与滚道接触点处的接触情况进行几何描述,并通过2维插值算法重新进行网格划分以便于后续的接触力学分析;对滚珠与滚道之间的接触区域和控制方程进行离散化分析,提高数值计算的精确性;为了解决接触问题计算求解过程复杂、收敛速度缓慢,通过运用变分原理将接触问题转成极值问题,采用共扼梯度法进行循环迭代来实现快速收敛求解,提高接触数值解法的计算效率,求解出滚珠与滚道之间的接触应力分布;由于接触表面变形的计算相当于求影响系数矩阵与法向压应力之间的卷积,2维卷积运算和对应的2维快速傅里叶变换被用于求解接触问题的数值解。在数值求解中,主要的计算量集中在弹性变形的计算上,使用2维快速傅里叶变换(FFT)技术方便快速地计算接触区域弹性变形。将最小余能原理与赫兹接触理论分别用于求解光滑的球与平面之间的弹性接触应力分布,将两种计算方法得出的结果进行对比,计算数值重合度很好,验证了基于最小余能原理求解接触应力分布的正确性;基于最小余能原理分析螺旋升角对滚珠与滚道接触点处的接触面积和弹性变形的影响,并将该方法所得结果与赫兹接触理论计算结果进行对比分析,发现赫兹接触解与最小余能非赫兹接触精确解相比,随着滚珠丝杠副滚道螺旋升角的增大,丝杠与滚珠接触点A处的误差始终大于螺母与滚珠接触点B处的误差,丝杠滚道与滚珠接触点A处非赫兹解接触应力峰值逐渐变小,螺母滚道与滚珠接触点B处非赫兹解接触应力峰值逐渐增大,且接触点A处应力峰值始终大于接触点B处峰值,计算所得滚珠与滚道接触区域的面积和长短轴有所不同。最小余能原理非赫兹接触计算方法可以准确计算滚珠和滚道接触区域应力分布,并且可以全面和准确的计算滚道磨损带的宽度和深度。研究结果表明,采用最小余能非赫兹接触解可以有效提高滚珠丝杠副接触应力分布计算精度,保证滚珠丝杠副的精度保持性。
Abstract:During the movement of the ball screw mechanism (BSM), the friction and wear between the balls and raceway cause the accuracy to deteriorate and reduce its accuracy retention. For the asymmetry of the contact area of the helical raceway of the high-speed BSM, there is a big error in the calculation of the contact geometry and relationship between the screw and the nut raceway by using the traditional Hertzian contact theory. Based on the minimum excess principle, a method was proposed to analyze and calculate the normal contact stress between the ball and the raceway of the high-speed BSM. Firstly, a BSM coordinate system was established based on the Frenet-Serret coordinate conversion formula to geometrically describe the contact between the ball and the raceway at the contact point, and then re-gridded by a two-dimensional interpolation algorithm to facilitate subsequent contact mechanics analysis; the contact area and control equation between the ball and raceway were discretized to improve the accuracy of numerical calculations. To solve the complex computational solution of the contact problem and the slow convergence, the contact problem was transformed by applying the variational fraction principle. The polar value problem was solved by using the conjugate gradient method for cyclic iteration to achieve fast convergence and improve the computational efficiency of the contact problem numerical solution, solving the contact stress distribution between the ball and the raceway. Since the calculation of the contact surface deformation was equivalent to the convolution between the influence coefficient matrix and the normal compressive stress, the two-dimensional convolution calculation, and the corresponding two-dimensional fast Fourier transform were used to solve the numerical solution of the contact problem. In numerical solutions, the main computational effort was focused on the computation of elastic deformations using the two-dimensional fast Fourier transform (FFT) technique Easy and fast calculation of elastic deformation in the contact area. Applying the minimum excess principle and Hertzian contact theory to solve the elastic contact stress distribution between a smooth ball and a plane, respectively, The results of the two calculation methods are compared, and the calculated values overlap well, which verifies the correctness of solving the contact stress distribution based on the principle of minimum excess; Based on the principle of minimum excess to analyze the influence of helix angle on the contact area and elastic deformation at the contact point of ball and raceway, and the results obtained by this method and Hertzian contact theoretical calculations for comparative analysis, as the BSM raceway helix angle increases, the error at the screw and ball contact point A is always greater than the error at the contact point of the nut and ball contact point B. The non-Hertzian solution contact stress peak at the contact point of the screw raceway and ball contact point A gradually becomes smaller, the non-Hertzian solution contact stress peak at the contact point of the nut raceway and ball contact point B gradually increases, and the stress peak at the contact point A is always greater than the peak at the contact point B. The non-Hertzian solution contact stress peak at the contact point of the nut raceway and ball contact point B gradually increases. The minimum excess principle non-Hertzian contact calculation method accurately calculates the stress distribution in the contact area of the ball and raceway and allows comprehensive and accurately calculate the width and depth of the wear zone of the raceway. Therefore, using the minimum excess principle non-Hertzian contact solution can improve the calculation accuracy of the contact stress distribution of the ball screw mechanism, and ensure the precision prediction retention of its accuracy.
文章编号:202000443     中图分类号:TH132.1    文献标志码:
基金项目:国家自然科学基金项目(51975020;51575014;51875008);北京市自然科学基金项目(3202005)
作者简介:第一作者:王民(1972-),男,教授,博士.研究方向:机床动力学;制造系统智能监控.E-mail:wangm@bjut.edu.cn
引用文本:
王民,孙铁伟,董朝阳,孔德顺,高相胜.高速滚珠丝杠副非赫兹接触应力计算[J].工程科学与技术,2021,53(2):178-186.
WANG Min,SUN Tiewei,DONG Zhaoyang,KONG Deshun,GAO Xiangsheng.Non-Hertzian Contact Stress Calculation of High-speed Ball Screw Mechanism[J].Advanced Engineering Sciences,2021,53(2):178-186.