一种基于有限元思想的进给系统刚柔耦合动力学建模方法技术方案

本发明专利技术公开了一种基于有限元思想的进给系统刚柔耦合动力学建模方法，包括如下步骤：(1)确定进给系统滚珠丝杠的几何参数和材料属性；(2)将进给系统滚珠丝杠依据划分方法划分为若干单元；(3)通过各单元的变形势能和各单元横向移动动能和转动动能，计算各单元的刚度矩阵和质量矩阵；(4)获取滚珠丝杠系统质量矩阵和刚度矩阵，构建无阻尼的滚珠丝杠系统的动力学微分方程；(5)根据进给系统的位置和运动关系，利用集中质量法通过拉格朗日方程建立进给系统的动力学模型。该方法充分考虑滚珠丝杠系统的弹性振动和柔性基础，降低了进给系统的自由度数目，提高了动力学微分方程的求解效率，对于进给系统动态特性的研究具有理论指导意义。

A Rigid-Flexible Coupled Dynamics Modeling Method for Feed System Based on Finite Element Method

The invention discloses a rigid-flexible coupling dynamic modeling method for feed system based on finite element method, which includes the following steps: (1) determining the geometric parameters and material properties of the feed system ball screw; (2) dividing the feed system ball screw into several units according to the partition method; (3) calculating each single unit through the deformation potential energy of each unit and the lateral movement kinetic energy and rotation kinetic energy of each unit. The element stiffness matrix and mass matrix; (4) Obtain the mass matrix and stiffness matrix of the ball screw system, and construct the dynamic differential equation of the undamped ball screw system; (5) According to the position and motion relationship of the feed system, the dynamic model of the feed system is established by using the lumped mass method through the Lagrange equation. The method takes full account of the elastic vibration and flexible foundation of the ball screw system, reduces the number of degrees of freedom of the feed system, improves the efficiency of solving dynamic differential equations, and has theoretical guiding significance for the study of dynamic characteristics of the feed system.

【技术实现步骤摘要】
一种基于有限元思想的进给系统刚柔耦合动力学建模方法
本专利技术涉及一种进给系统刚柔耦合动力学建模方法，特别是涉及一种基于有限元思想的进给系统刚柔耦合动力学建模方法。
技术介绍
近年来,高速度、高加速度和高精度已经成为数控机床的主要发展方向。进给系统作为数控机床的核心功能部件，其性能的高低直接影响和制约着数控机床的发展。进给系统的动态特性是造成数控机床动态误差的主要因素，因此建立精准的进给系统动力学模型是进行后续动力学特性研究和动态误差影响规律研究的基础。目前，进给系统的动力学建模方法主要有集中参数法和有限元法。其中，集中参数法将进给系统的各零部件简化为集中质量块，将各零部件结合面间的刚度利用弹簧阻尼单元进行等效，该方法基于刚体参数化动力学建模，无法真实反应各关键结构件的变形。有限元法是在集中参数法基础上，将滚珠丝杠结构件等效为离散单元进行分析，但该方法建立的模型结构复杂，自由度数目偏多，因此计算困难。
技术实现思路
本专利技术的目的在于克服现有方法的问题，提供了一种基于有限元思想的进给系统刚柔耦合动力学建模方法，在进给系统动态特性研究阶段搭建进给系统刚柔耦合动力学模型，充分考虑滚珠丝杠系统的弹性振动和柔性基础，提高计算效率。本专利技术的一种基于有限元思想的进给系统刚柔耦合动力学建模方法，包括以下步骤：步骤一、确定进给系统滚珠丝杠的几何参数和材料属性，滚珠丝杠的几何参数包括各光轴部分轴外径直径和长度，带螺纹部分公称直径和长度；滚珠丝杠的材料属性包括弹性模量、泊松比和密度；步骤二、将进给系统滚珠丝杠划分为若干单元，划分方法如下：将滚珠丝杠轴外径变化处作为节点，将位于两个相邻节点之间且轴外径一致的轴段作为一个单元；当滚珠丝杠某一单元轴段长度大于该单元截面直径设定倍数时，将该轴段划分为若干单元；将轴外径相同的滚珠丝杠光轴部分和带有螺纹部分划分为两个单元；将滚珠丝杠与螺母结合部分划分为一个单元，滚珠丝杠与螺母结合处的左右两端为节点；将滚珠丝杠两端的轴承支撑点取为节点，该节点设置在支撑轴承宽度的中点处；步骤三、根据各相邻单元的节点处横向位移和弹性转角的相互作用，计算各单元的变形势能和各单元横向移动动能和转动动能，通过各单元的变形势能计算各单元刚度矩阵，通过各单元横向移动动能和转动动能计算各单元质量矩阵；步骤四、将各单元的质量矩阵和刚度矩阵分别进行叠加，得到滚珠丝杠系统质量矩阵和刚度矩阵，构建无阻尼的滚珠丝杠系统的动力学微分方程；步骤五、将伺服电机、联轴器、螺母、工作台视为具有集中质量的刚体，基于有限元思想将滚珠丝杠视为柔体，根据进给系统的位置和运动关系，在伺服电机的连接轴和联轴器接触部分建立第一结合面，在联轴器和滚珠丝杠接触部分建立第二结合面，在滚珠丝杠和联轴器支撑轴承接触部分建立第三结合面，在螺母和滚珠丝杠接触部分建立第四结合面，在螺母和工作台接触部分建立第五结合面，在工作台和导轨接触部分建立第六结合面，各结合面处通过假想的弹簧阻尼单元连接，实现结合部的等效建模，最终利用集中质量法通过拉格朗日方程建立进给系统的动力学模型。与现有方法相比，本专利技术的基于有限元思想的进给系统刚柔耦合动力学建模方法将滚珠丝杠系统根据滚珠丝杠轴径的变化和滚珠丝杠轴承安装位置的不同划分为若干单元，通过各单元的变形势能和各单元横向移动动能和转动动能，分别计算各单元的刚度矩阵和质量矩阵，将各个单元的矩阵进行叠加形成无阻尼的滚珠丝杠系统的动力学微分方程，最后基于集中参数法构建进给系统动力学模型。该方法充分考虑滚珠丝杠系统的弹性振动和柔性基础，极大限度的降低了进给系统的自由度数目，提高了动力学微分方程的求解效率，对于进给系统动态特性的研究具有一定的理论指导意义。附图说明图1为基于有限元思想的进给系统刚柔耦合动力学建模方法的整体流程图；图2为滚珠丝杠系统单元划分示意图；图3为第i个单元节点处横向位移和弹性转角的相互作用示意图；图4为进给系统动力学模型示意图。具体实施方式下面结合实施图和附图对本专利技术做出详细说明。本专利技术的一种基于有限元思想的进给系统刚柔耦合动力学建模方法，整体流程图如图1所示。该流程包括以下步骤：步骤一、确定进给系统滚珠丝杠的几何参数和材料属性。滚珠丝杠的几何参数包括各光轴部分轴外径直径和长度，带螺纹部分公称直径和长度；滚珠丝杠的材料属性包括弹性模量、泊松比和密度。步骤二、将进给系统滚珠丝杠划分为若干单元，划分方法如下：(1)将滚珠丝杠轴外径变化处作为节点，将位于两个相邻节点之间且轴外径一致的轴段作为一个单元；(2)当滚珠丝杠某一单元轴段长度大于该单元截面直径设定倍数(1.5倍)时，将该轴段划分为若干单元；(3)将轴外径相同的滚珠丝杠光轴部分和带有螺纹部分划分为两个单元；(4)将滚珠丝杠与螺母结合部分划分为一个单元，滚珠丝杠与螺母结合处的左右两端为节点；(5)将滚珠丝杠两端的轴承支撑点取为节点(轴承安装处不必划分为一个单元，只需要根据轴径变化划分，但轴承支承点为一个节点)，且节点设置在支撑轴承宽度的中点处。作为本专利技术的一种实施方式，如图2所示，在该滚珠丝杠子结构中，将其划分为14个单元，分别为轴承安装Ⅰ处的单元1、单元2，轴肩部分的单元3，光轴部分单元4、单元5，螺纹部分Ⅰ的单元6、单元7、单元8，滚柱丝杠与螺母结合部分的单元9，螺纹部分Ⅱ的单元10、单元11、单元12和轴承安装Ⅱ处的单元13、单元14。步骤三、根据各相邻单元的节点处横向位移和弹性转角的相互作用，计算各单元的变形势能和各单元横向移动动能和转动动能，通过各单元的变形势能计算各单元刚度矩阵，通过各单元横向移动动能和转动动能计算各单元质量矩阵。各单元刚度矩阵和各单元质量矩阵的具体计算方法如下：第一步，设滚珠丝杠系统单元个数为N，单元间节点数目为Nn，单元个数与节点个数之间存在以下关系：Nn＝N+1将各单元和节点自左向右进行编号，每个节点处建立两个广义坐标为横向位移和弹性转角。在第i个节点处建立的横向位移广义坐标编号为U2i-1，建立的弹性转角广义坐标编号为U2i，最终形成广义坐标列阵U＝[U1U2…U2i-1U2i]T。广义坐标数目Nu＝2Nn。第二步，分析第i个单元的节点处横向位移和弹性转角的变化，如图3所示。第i个单元包括有两个节点，将左节点处的横向位移和弹性转角设为u1(t)和u2(t)，将右节点处的横向位移和弹性转角设为u3(t)和u4(t)，利用u(x,t)表示在该单元任意径向截面处的横向位移，该单元在受到载荷作用时变形较为复杂，因此假定横向振动单元的位移取为三次多项式u(x,t)＝c0+c1x+c2x2+c3x3，利用节点两侧的横向位移和弹性转角u1(t)、u2(t)、u3(t)和u4(t)的边界条件，该边界条件u1(t)、u2(t)、u3(t)和u4(t)的计算方法为现有方法，具体参见文献(王世斌,亢一澜,王燕群,等.普通高等教育“十一五”国家级规划教材,材料力学[M].高等教育出版社,2008.)，获取三次多项式中的待定系数c0、c1、c2和c3，至此，求的任意截面y处的横向位移u(x,t)。该横向位移u(x,t)的计算方法为现有方法，具体参见文献(张策.机械动力学,(第二版)[M].高等教育出版社,2008.)第三步，计算第i个单元的变形势能，并通过第i个单本文档来自技高网...

【技术保护点】
1.一种基于有限元思想的进给系统刚柔耦合动力学建模方法，其特征在于包括以下步骤：步骤一、确定进给系统滚珠丝杠的几何参数和材料属性，滚珠丝杠的几何参数包括各光轴部分轴外径直径和长度，带螺纹部分公称直径和长度；滚珠丝杠的材料属性包括弹性模量、泊松比和密度；步骤二、将进给系统滚珠丝杠划分为若干单元，划分方法如下：将滚珠丝杠轴外径变化处作为节点，将位于两个相邻节点之间且轴外径一致的轴段作为一个单元；当滚珠丝杠某一单元轴段长度大于该单元截面直径设定倍数时，将该轴段划分为若干单元；将轴外径相同的滚珠丝杠光轴部分和带有螺纹部分划分为两个单元；将滚珠丝杠与螺母结合部分划分为一个单元，滚珠丝杠与螺母结合处的左右两端为节点；将滚珠丝杠两端的轴承支撑点取为节点，该节点设置在支撑轴承宽度的中点处；步骤三、根据各相邻单元的节点处横向位移和弹性转角的相互作用，计算各单元的变形势能和各单元横向移动动能和转动动能，通过各单元的变形势能计算各单元刚度矩阵，通过各单元横向移动动能和转动动能计算各单元质量矩阵；步骤四、将各单元的质量矩阵和刚度矩阵分别进行叠加，得到滚珠丝杠系统质量矩阵和刚度矩阵，构建无阻尼的滚珠丝杠系统的动力学微分方程；步骤五、将伺服电机、联轴器、螺母、工作台视为具有集中质量的刚体，基于有限元思想将滚珠丝杠视为柔体，根据进给系统的位置和运动关系，在伺服电机的连接轴和联轴器接触部分建立第一结合面，在联轴器和滚珠丝杠接触部分建立第二结合面，在滚珠丝杠和联轴器支撑轴承接触部分建立第三结合面，在螺母和滚珠丝杠接触部分建立第四结合面，在螺母和工作台接触部分建立第五结合面，在工作台和导轨接触部分建立第六结合面，各结合面处通过假想的弹簧阻尼单元连接，实现结合部的等效建模，最终利用集中质量法通过拉格朗日方程建立进给系统的动力学模型。...

【技术特征摘要】
1.一种基于有限元思想的进给系统刚柔耦合动力学建模方法，其特征在于包括以下步骤：步骤一、确定进给系统滚珠丝杠的几何参数和材料属性，滚珠丝杠的几何参数包括各光轴部分轴外径直径和长度，带螺纹部分公称直径和长度；滚珠丝杠的材料属性包括弹性模量、泊松比和密度；步骤二、将进给系统滚珠丝杠划分为若干单元，划分方法如下：将滚珠丝杠轴外径变化处作为节点，将位于两个相邻节点之间且轴外径一致的轴段作为一个单元；当滚珠丝杠某一单元轴段长度大于该单元截面直径设定倍数时，将该轴段划分为若干单元；将轴外径相同的滚珠丝杠光轴部分和带有螺纹部分划分为两个单元；将滚珠丝杠与螺母结合部分划分为一个单元，滚珠丝杠与螺母结合处的左右两端为节点；将滚珠丝杠两端的轴承支撑点取为节点，该节点设置在支撑轴承宽度的中点处；步骤三、根据各相邻单元的节点处横向位移和弹性转角的相互作用，计算各单元的变形势...

【专利技术属性】
技术研发人员：罗茹楠，牛文铁，王晨升，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津,12