双丝杠驱动进给单元轴向刚度模型及影响因素分析

双丝杠驱动进给单元轴向刚度模型及影响因素分析许向荣;鹿群鹏;宋现春;姜洪奎【摘要】双丝杠驱动进给单元在高性能数控机床中的应用较为广泛,研究其动态特性对于提高数控机床的加工性能具有重要的意义.文章基于关键滚动功能部件结合面刚度各个方向耦合的影响,围绕双丝杠驱动进给单元轴向刚度进行了研究分析,通过研究双丝杠驱动进给单元的轴向刚度模型,结合动力学和Hertz接触理论的相关知识,推导其各组成部分的刚度计算公式,探索其主要影响因素及关系曲线.结果表明:双丝杠驱动进给单元的轴向刚度主要受到滚珠丝杠副、支承轴承的结构参数,螺母位置以及两根丝杠跨距的影响;支承轴承结构参数的变化,会引起进给单元轴向刚度的明显变化;螺母位置是进给单元轴向刚度的主要影响因素,在两端固定的方式下其轴向刚度中间最小,在机床加工过程中,应尽量避免在丝杠的中间位置进行操作.【期刊名称】《山东建筑大学学报》【年(卷),期】2016(031)002【总页数】8页(P125-132)【关键词】结合面;双丝杠驱动;进给单元;轴向刚度【作者】许向荣;鹿群鹏;宋现春;姜洪奎【作者单位】山东建筑大学机电工程学院,山东济南250101;山东建筑大学机电工程学院,山东济南250101;山东建筑大学机电工程学院,山东济南250101;山东建筑大学机电工程学院,山东济南250101【正文语种】中文【中图分类】TH132.1与传统的单丝杠驱动进给系统相比，双丝杠驱动进给系统有很大的优势，如高精度、高刚度、低振动等，因此，双丝杠驱动进给系统日益成为高性能数控机床、加工中心研究领域的一个热点[1-3]。

目前，国内外学者纷纷围绕双丝杠驱动进给系统展开了理论和实验方面的研究。

Hiramoto等介绍了重心驱动原理，设计了工作台和主轴均采用重心驱动的机床模型和实际的双驱动进给系统的实验模型，通过实验发现重心驱动进给系统可以明显地减小机床的振动[1]，但是缺乏理论支持。

Gomand等采用理论建模的方法，建立了龙门数控机床双驱动进给系统的数学模型，推导了丝杠轴向振动方程，并分析了丝杠变形对其静动态特性的影响[2]。

而国内这方面的研究起点较晚，但随着高性能数控机床的快速发展，各科研院所也加快了对双驱动进给系统的研究。

何王勇等仅考虑双滚珠丝杠同步驱动轴的同步轴间负载耦合关系及丝杠、螺母的轴向扭转耦合关系的影响，建立了其动力学模型，推导了动力学方程，并对其固有频率进行了理论计算[3]。

唐余林和芮执元等以铣车复合加工中心双驱进给系统为研究对象，将双丝杠耦合转化到横梁的各个方向上的运动，建立了双驱进给系统动力学模型及其动力学方程，并进行了模态分析与谐响应分析，找到了系统的薄弱环节，并对其进行了优化[4-5]。

郭崇高等采用有限元分析和虚拟样机技术，建立了双驱动进给系统动力学模型，采用拉格朗日方程，推导了系统运动微分方程，获得了其主要模态参数，最后对其进行了模态分析和谐响应分析，得出两对导轨滑块与两根丝杠跨距等对系统动态性能的影响，并提出了相应的改进措施[6]。

丁喜合等建立了二自由度系统振动模型，采用单因素分析方法研究了丝杠跨距和导轨跨距对双丝杠驱动进给系统动态特性的影响[7-9]。

但是对于导轨跨距的影响研究中没有考虑导轨法向结合面的刚度。

白茹通过拉格朗日方程在广义坐标下的变换建立了双驱滚珠丝杠进给系统的实体模型，并设计了加工中心双驱进给系统动态性能实验平台，对双驱滚珠丝杠进给系统的动态特性进行了分析，提出了双驱进给系统的改善措施,为进一步研究奠定了理论基础[10]。

夏田等建立了双丝杠驱动工作台动力学模型，找到了位置响应的主要因素，对重心驱动工作台结构进行了改进，并建立了有限元分析模型，进行了模态分析和谐响应分析。

其研究发现进给轴刚度是影响位置精度的关键因素，恰当选择导轨的间距对提高工作台的动态响应性能有重要的作用[11]。

杨勇等采用理论建模的方法，对MCH63卧式加工中心“重心驱动”机构建立了动力学模型，通过实验获得了导轨结合面的刚度和阻尼值，并计算分析了导轨结合面刚度、阻尼值及“重心驱动”结构对机床动力学性能的影响[12]。

上述对于双丝杠驱动进给系统的研究,多数采用建立动力学模型，推导动力学方程，有限元模态分析和谐响应分析的思路对进给系统动态特性进行研究，没有专门针对系统刚度进行的研究。

因此，文章考虑结合面刚度各个方向耦合的影响，建立了双丝杠驱动进给单元的动力学模型，并在此基础上建立了其轴向刚度模型，分析研究了影响进给单元刚度的主要因素，获得了其相应的关系曲线，为后面分析双丝杠驱动进给单元的动态特性奠定了基础。

双丝杠驱动进给单元结构组成示意图，如图1(a)所示，主要包括：两台相同型号的伺服电机，两套相同型号的滚珠丝杠副，两套相同型号的支承轴承，两套相同型号的直线滚动导轨副以及拖板组成。

由两台伺服电机同步驱动两套滚珠丝杠副带动拖板沿着直线滚动导轨副做往复直线运动，实现进给驱动。

由于60%的机床振动产生于其结合部[13-14]，因此在研究双丝杠驱动进给单元的动态特性时，不能忽略结合面的影响，尤其是其主要滚动功能部件，如滚珠丝杠副、滚动轴承、直线滚动导轨副结合部的影响。

这三种部件作为双丝杠驱动进给单元的关键部件，其变形及刚度对整个系统的变形及刚度特性都有很重要的影响。

因此，在该部分里，文章依据动力学相关理论，建立了双丝杠驱动进给单元动力学模型，如图1(b)所示。

并在此基础上建立了进给单元的轴向刚度模型。

然后，基于Hertz 接触理论[15]，推导了滚珠丝杠副、支承轴承结合面的接触变形和接触刚度的计算公式，进而获得进给单元的轴向刚度计算公式；最后，通过实例计算分析研究了其各个设计参数对进给单元轴向刚度的影响。

依据图1(b),可以建立双丝杠驱动进给单元的轴向刚度模型，如图1(c)所示。

从图1(c)可以看出，由弹簧串、并联的刚度系数计算方法可得出双丝杠驱动进给单元的轴向刚度由式(1)、(2)表示为式中：Kfi(i=1、2)为第i根滚珠丝杠传动系统的轴向刚度，N/mm，可以根据式(2)进行计算，Ksi(i=1、2)是滚珠丝杠副结合面的轴向刚度，N/mm；kφsi(i=1、2)是滚珠丝杠扭转—轴向耦合刚度，N/mm；Kbbil(i=1、2)、Kbbir(i=1、2)分别是左、右端支承轴承结合面的轴向刚度，N/mm；Kφbbil(i=1、2)、Kbbir(i=1、2)分别是左端轴承扭转—轴向耦合刚度，N/mm；Ksgi(i=1、2)是由拖板耦合引起的滚珠丝杠轴向刚度，N/mm。

2.1 滚珠丝杠副结合面的轴向刚度计算相关研究表明，两个曲面物体互相接触、挤压，其在接触部位的应力分布与接触面的形状、尺寸及表面粗糙度等许多因素有关[14]。

1881年，Hertz指出两光滑弹性体在相互作用的力下发生的接触变形及应力分布应满足Hertz接触假设[15]：(1) 接触物体只产生弹性变形，服从胡克定律。

(2) 接触表面充分光滑，不考虑接触面的介质以及动摩擦力的影响。

(3) 接触面的尺寸与接触物体表面的曲率半径相比是小量。

满足上述假设的接触为赫兹接触。

在上述假设下，忽略工作载荷分布不均匀的影响，单螺母丝杠副在法向载荷Q作用下，螺母滚道面与丝杠滚道面间由于法向弹性接触变形所产生的弹性变形量[16-18]δn由式(3)表示为式中：E′为等效弹性模量，)；E1、E2分别为丝杠和螺母的弹性模量；μ1、μ2分别为丝杠和螺母的泊松比；为丝杠侧和螺母侧的无量纲参数，仅与接触椭圆的椭圆参数有关；δs p、δnp分别为滚珠与丝杠滚道面、与螺母滚道面的接触弹性变形量；∑ρs、∑ρn分别是丝杠侧和螺母侧滚道主曲率的和。

根据几何关系，可以得到丝杠相对于螺母的轴向弹性位移δa由式(4)表示为将式(3)带入式(4)，可得单螺母滚珠丝杠副的载荷变形关系由式(5)表示为式中：Fa为螺母受到的轴向推力，N，Fa=zQsinβcosλ；β、λ分别为滚珠丝杠副的接触角和螺旋升角，°；z为滚珠丝杠副的滚珠数目。

将式(5)两边分别对δa求导，可得单螺母丝杠副的轴向接触刚度式(6)为通常采用双螺母预紧的方法减小或消除滚珠丝杠副的轴向间隙，以提高其轴向接触刚度。

因此，在上述单螺母滚珠丝杠副研究上，可以得到双螺母滚珠丝杠副结合面的轴向接触刚度Ksi(i=1,2)[13]由式(7)表示为式中：PP为丝杠上承受的预紧法向力，N；F为轴向载荷，N。

由式(7)可以看出，在轴向载荷作用下，双螺母滚珠丝杠副结合面的轴向接触刚度与轴向载荷、预紧力、螺旋升角、接触角等几何参数有关。

2.2 拖板引起的滚珠丝杠轴向刚度计算工作过程中，除了轴向载荷会引起单根丝杠轴向变形外，拖板会使得两根丝杠产生耦合变形，进而影响到进给单元的轴向动态变形量。

在轴向载荷作用下，丝杠会发生拉压变形，其大小与两端的支承方式及螺母的工作位置有关。

设支承方式为两端固定方式，由于两根滚珠丝杠通过拖板连接在一起，假设两根丝杠不能实现同步运动，则势必会出现如图2所示的耦合作用。

图2中，lsi(i=1,2)为丝杠i支承与螺母间的距离，mm；lsgi(i=1,2)为丝杠i的长度，mm；ly为两根滚珠丝杠跨距，mm。

由图2可以看出几何关系式(8)为因此，对于单根滚珠丝杠来说，其轴向变形量式(9)、(10)为式中：δs1、δs2分别为丝杠1、2的轴向直线位移量，mm；F为丝杠的轴向载荷，N；d1为丝杠的螺纹内径，mm；θ为两根丝杠不同步时，拖板的扭转角度；lsi(i=1,2)为丝杠i支承与螺母间的距离，mm；lsgi(i=1,2)为丝杠i的长度，mm；ly为两根滚珠丝杠跨距，mm；E为材料的纵向弹性模量，钢的E=207 GPa。

因此，由拖板引起的滚珠丝杠轴向刚度由式(11)、(12)表示为从式(11)，(12)中可以看出，影响进给单元中拖板引起的滚珠丝杠轴向刚度的主要因素有：两根丝杠的跨距ly，螺母的位置ls1以及两根丝杠不同步时拖板的扭转角度θ。

2.3 滚珠丝杠的扭转—轴向耦合刚度计算在驱动力矩和轴向载荷作用下，滚珠丝杠时会产生扭转变形，从而引起螺母轴向位移，因而会影响丝杠轴向动态变形量。

根据材料力学相关知识，可得两端固定时单根滚珠丝杠在扭矩作用下产生的轴向变形量由式(13)表示为式中：θi(i=1,2)为丝杠i产生的扭转变形角，rad/mm；T为丝杠的驱动扭矩；G为材料剪切弹性模量，钢的G=8.5×104N/mm2； d1为丝杠螺纹内径，mm；F 为轴向工作载荷；η为传动效率；lswifti(i=1,2)为丝杠i受扭部分长度，一般为两端支承跨距lsgi(i=1,2)；Pl为导程。

因此，丝杠的扭转—轴向耦合刚度式由(14)表示为2.4 支承轴承结合面的轴向刚度计算鉴于滚动轴承与滚珠丝杠副结构和受力的相似性，之前求滚珠丝杠副轴向接触刚度计算时的假设条件，仍适用于支承系统轴向刚度的计算，且单滚珠丝杠副轴向接触刚度计算分析方法也同样适用于支承系统轴向刚度的计算。