一种直角切削颤振解析建模方法技术

本发明专利技术涉及一种直角切削颤振解析建模方法，包括以下步骤：步骤1，建立直角切削颤振的动力学模型；步骤2，计算动态切削力和切削力系数；步骤3，计算直角切削过程的稳定性SLD图，本发明专利技术将动态切削过程看作是在每一时刻的准静态切削过程，将工件材料特性、刀具几何、切削参数作为输入参数，其中动态切削力可通过等效变换的切削参数计算，进而理论推导出动态切削力系数的表达式，通过代数方程整理，得出动态切削力系数，避免了实验标定的繁琐性，并提高了准确度，另外，机床系统切削颤振的动态模型可以用时滞微分方程表达，通过进行切削稳定性分析，由时域半离散法获得直角切削颤振稳定性SLD图，提供一个比较真实的颤振稳定性预测。

An analytical modeling method for right angle cutting chatter

The invention relates to a cutting chatter analytical modeling method, which comprises the following steps: 1, to establish dynamic model of cutting chatter; step 2, the calculation of dynamic cutting force and cutting force coefficient; step 3, calculate the orthogonal cutting process stability SLD diagram, the dynamic cutting process as in the quasi static cutting process every moment, the workpiece material properties, tool geometry and cutting parameters as the input parameters, the dynamic cutting force by cutting parameter equivalent transform calculation, and then the theory of dynamic cutting force coefficients derived from the expression, by means of algebraic equations, the dynamic cutting force coefficients, avoiding the tedious calibration test. And improve the accuracy, in addition, the system dynamic model of cutting chatter of machine tools can be expressed by differential equation, the stability of cutting The stability of the Cartesian cutting chatter is obtained by the time domain semi discrete method, which provides a more realistic prediction of the flutter stability of SLD.

【技术实现步骤摘要】
一种直角切削颤振解析建模方法
本专利技术涉及金属材料切削加工
，尤其涉及一种直角切削颤振解析建模方法。
技术介绍
切削颤振是由于切削力的变化以及加工系统的柔性而产生的切屑再生现象，属于自激振动的一种形式。当不加以控制时，其会对加工工件表面质量，材料去除率，表面尺寸精度，以及刀具与机床的寿命产生破坏性影响。因此，建立有效的颤振模型是预测与避免切削颤振发生的有效措施。迄今为止，已经形成了很多加工颤振模型，其采用的方法主要有三种：解析法、数值法及实验法，其中切削力系数通常是通过实验方法确定的常数。很明显，对于给定范围的切削参数，获得可靠的实验数据是非常耗时的。此外，切削力系数在实际切削过程中是变化的，经验常数的假设将产生不准确的结果。除此之外，这些模型在缺乏切削参数(工件材料特性、刀具几何)对切削颤振影响的直观描述，且仍依赖于经验或实验的方法。针对现有技术方案提出的建模方法的技术问题，本专利技术提出了一种直角切削颤振的解析建模方法，该方法从金属切削的机理入手，将工件材料特性、刀具几何及切削参数作为输入参数，通过对切削参数的等效变换将动态切削过程看作是在每一时刻的准静态切削过程，动态切削力系数可通过等效变换来解析标定，避免了实验标定的繁琐性。
技术实现思路
针对上述现有技术的现状，本专利技术所要解决的技术问题在于提供一种直角切削颤振解析建模方法，该方法可以克服切削力系数由实验标定的准确度太低太繁琐的问题，避免了实验标定的繁琐性，提高了切削力系数的准确度，提供一个比较真实的切削颤振稳定性的预测。本专利技术解决上述技术问题所采用的技术方案为：一种直角切削颤振解析建模方法，包括以下步骤：步骤1，建立直角切削颤振的动态模型确定刀具几何参数：前角α，单位为deg，后角γ，单位为deg；选定切削参数，切削速度V，单位为m/min，进给量ft，单位为mm/r、切削宽度b，单位为mm；在金属切削过程中，切削厚度的变化会导致切削力的波动，切削力和切削厚度是按照特定的周期发生变化的，并形成一个闭环反馈系统，通过公式(1)计算出动态切削厚度h(t)；h(t)＝h0-[y(t)-y(t-T)](1)其中，h0为名义切削厚度，单位为mm，数值上等于刀具的进给量ft；y(t)和y(t-T)分别表示y方向的当前和前一周期的振动幅值，单位为mm，也称为内调制和外调制，[y(t)-y(t-T)]表示在t时刻由于刀具振动产生的切削厚度差，T为当前时刻与前一时刻的延迟，即主轴旋转周期；在工件被认为是刚性的而刀具是柔性的情况下，机床系统切削颤振的动态模型为在进给方向的单自由度系统，进给方向即y方向，该系统的振动方程由公式(2)表示；其中，m表示系统的等效质量，单位为kg，c表示系统的阻尼，单位为Ns/m，k表示系统的等效刚度，单位为N/m，Fy(t)为在进给方向由动态切削厚度h(t)引起的动态切削力；步骤2，计算动态切削力和切削力系数动态切削力由公式(3)表示；Fx(t)＝kx(t)bh(t),Fy(t)＝ky(t)bh(t)(3)其中，kx(t),ky(t)分别表示x，y方向的动态切削力系数，表示为名义切削速度V与刀具振动速度的函数，分别由公式(4)和公式(5)表示；其中，cxi，cyi，i＝0,1,2是名义切削速度V的函数，通过给定范围内的切削参数几何标定；将切屑的形成过程看作一个准静态的过程，其切削波纹面的斜率ξ即为准静态方向与有效瞬时切削方向之间的方位角，准静态方向为x方向，有效瞬时切削方向为x′方向，斜率ξ由公式(6)表示；其中，V为名义切削速度，单位为m/min，为y方向的刀具振动速度，单位为m/min，ξ的正负号与的一致；然后，t时刻的有效前角α′，单位为deg，有效后角γ′，单位为deg,，有效切削厚度h′(t)，单位为mm，通过准静态切削和动态切削之间的变化关系得到，由公式(7)和公式(8)表示；α′＝α-ξ，γ′＝γ+ξ(7)其中，名义前角α和名义后角γ分别对应于准静态切削过程的前角和后角，单位为deg；瞬时剪切角φ′、有效切削速度V′与有效切屑速度V′c分别由公式(9)和公式(10)计算；其中，为动态摩擦角，单位为deg，该动态摩擦角由公式(11)计算：其中，为名义摩擦系数，p为指数参数，V'c为有效切屑速度，单位为m/min；通过以上公式的换算，得到瞬时剪切角φ′的最终表达公式(12)；其中，A1，A2为材料特性参数，无量纲单位，瞬时剪切角φ′依赖于名义切削速度V和刀具振动速度通过Newton-Raphson迭代算法求解；有效平均流动剪切应力由切削材料的Johnson-Cook本构模型确定，有效剪切力Fs′，x，y方向的切削力分别由公式(13)和公式(14)计算；在每一时刻的准静态，计算出切削波纹面的斜率ξ、有效前角α′、有效后角γ′和有效切削厚度h′(t)、有效切削速度V′、有效切屑速度V′c、瞬时剪切角φ′，进而通过代数方程整理，得出x、y方向的动态切削力系数kx(t),ky(t)；步骤3，计算直角切削过程的稳定性SLD图首先，结合公式(2)和公式(5)，并在平衡位置处引入微小振动量u(t)，得到机床系统切削颤振的动态模型整理后的表达公式(15)；然后通过在公式(15)所示的平衡位置处引入微小振动量通过忽略公式(15)右端的非线性部分的影响，进行颤振的线性分析，得到公式(16)；令整理得到公式(15)的状态空间方程公式(17)；其中，通过数值求解公式(17)，将时间周期T离散为n个等分小区间，即T＝nτ，在每一个区间[iτ,(i+1)τ](i＝0,...,n-1)的中间矩阵Φ通过使用连续的离散映射Di，构造公式(18)如下；yn＝Φy0＝Dn-1Dn-2…D1D0y0(18)其中，Di为离散映射矩阵，yi为2(n+1)列向量；然后根据Floquet理论，即具有周期系数的线性常微分方程，通过约化成为一个常系数的常微分方程，判别颤振系统的稳定性，当矩阵Φ的所有特征值的模量小于单位1，则系统是稳定的，否则是不稳定的，进而计算出切削速度与切削宽度的关系图，即为切削颤振系统的稳定性SLD图。与现有技术相比，本专利技术的优点在于：本专利技术将动态切削过程看作是在每一时刻的准静态切削过程，将工件材料特性、刀具几何、切削参数作为输入参数，其中动态切削力可通过等效变换的切削参数计算，进而理论推导出动态切削力系数的表达式，通过代数方程整理，得出动态切削力系数，避免了实验标定的繁琐性，并提高了准确度，另外，机床系统切削颤振的动态模型可以用时滞微分方程表达，通过进行切削稳定性分析，由时域半离散法获得直角切削颤振稳定性SLD图，提供一个比较真实的颤振稳定性预测。附图说明图1是本专利技术实施例中直角切削颤振模型示意图；图2是本专利技术实施例中直角切削颤振原理图；图3是本专利技术实施例中等效参数变换图；图4是本专利技术实施例中直角颤振SLD对比图。具体实施方式下面结合附图并通过实施例对本专利技术作进一步的详细说明，以下实施例是对本专利技术的解释而本专利技术并不局限于以下实施例。如图1～3所示，一种直角切削颤振解析建模方法，包括以下步骤：步骤1，建立直角切削颤振的动态模型确定刀具几何参数，前角α，单位为deg，后角γ，单位为deg；选定切削参数，切削速度V，单位为m/min，进给量ft，本文档来自技高网...

【技术保护点】
一种直角切削颤振解析建模方法，其特征在于，包括以下步骤：步骤1，建立直角切削颤振的动态模型确定刀具几何参数：前角α，单位为deg，后角γ，单位为deg；选定切削参数，切削速度V，单位为m/min，进给量ft，单位为mm/r、切削宽度b，单位为mm；在金属切削过程中，切削厚度的变化会导致切削力的波动，切削力和切削厚度是按照特定的周期发生变化的，并形成一个闭环反馈系统，通过公式(1)计算出动态切削厚度h(t)；h(t)＝h0‑[y(t)‑y(t‑T)]    (1)其中，h0为名义切削厚度，单位为mm，数值上等于刀具的进给量ft；y(t)和y(t‑T)分别表示y方向的当前和前一周期的振动幅值，单位为mm，也称为内调制和外调制，[y(t)‑y(t‑T)]表示在t时刻由于刀具振动产生的切削厚度差，T为当前时刻与前一时刻的延迟，即主轴旋转周期；在工件被认为是刚性的而刀具是柔性的情况下，机床系统切削颤振的动态模型为在进给方向的单自由度系统，进给方向即y方向，该系统的振动方程由公式(2)表示；

【技术特征摘要】
1.一种直角切削颤振解析建模方法，其特征在于，包括以下步骤：步骤1，建立直角切削颤振的动态模型确定刀具几何参数：前角α，单位为deg，后角γ，单位为deg；选定切削参数，切削速度V，单位为m/min，进给量ft，单位为mm/r、切削宽度b，单位为mm；在金属切削过程中，切削厚度的变化会导致切削力的波动，切削力和切削厚度是按照特定的周期发生变化的，并形成一个闭环反馈系统，通过公式(1)计算出动态切削厚度h(t)；h(t)＝h0-[y(t)-y(t-T)](1)其中，h0为名义切削厚度，单位为mm，数值上等于刀具的进给量ft；y(t)和y(t-T)分别表示y方向的当前和前一周期的振动幅值，单位为mm，也称为内调制和外调制，[y(t)-y(t-T)]表示在t时刻由于刀具振动产生的切削厚度差，T为当前时刻与前一时刻的延迟，即主轴旋转周期；在工件被认为是刚性的而刀具是柔性的情况下，机床系统切削颤振的动态模型为在进给方向的单自由度系统，进给方向即y方向，该系统的振动方程由公式(2)表示；其中，m表示系统的等效质量，单位为kg，c表示系统的阻尼，单位为Ns/m，k表示系统的等效刚度，单位为N/m，Fy(t)为在进给方向由动态切削厚度h(t)引起的动态切削力；步骤2，计算动态切削力和切削力系数动态切削力由公式(3)表示；Fx(t)＝kx(t)bh(t),Fy(t)＝ky(t)bh(t)(3)其中，kx(t),ky(t)分别表示x，y方向的动态切削力系数，表示为名义切削速度V与刀具振动速度的函数，分别由公式(4)和公式(5)表示；其中，cxi，cyi，i＝0,1,2是名义切削速度V的函数，通过给定范围内的切削参数几何标定；将切屑的形成过程看作一个准静态的过程，其切削波纹面的斜率ξ即为准静态方向与有效瞬时切削方向之间的方位角，准静态方向为x方向，有效瞬时切削方向为x′方向，斜率ξ由公式(6)表示；其中，V为名义切削速度，单位为m/min，为y方向的刀具振动速度，单位为m/min，ξ的正负号与的一致；然后，t时刻的有效前角α′，单位为deg，有效后角γ′，单位为deg,，有效切削厚度h′(t)，单位为mm，通过准静态切削和动态切削之间的变化关系得到，由公式(7)和公式(8)表示；α′＝α-ξ，γ′＝γ+ξ(7)其中，名义前角α和名义后角γ分别对应于准静态切削过程的前角和后角，单位为deg；瞬时剪切角φ′、有效切削速度V′与有效切屑速度Vc′分别由公式(9)和公式(10)计算；其中，为动态摩擦角，单位为deg，该动态摩擦角由公式(11)计算：

【专利技术属性】
技术研发人员：付中涛，陈绪兵，戴耀南，毛金城，
申请(专利权)人：武汉工程大学，
类型：发明
国别省市：湖北,42