基于工作模态分析的铣削过程阻尼标定方法与流程

本发明涉及一种铣削过程阻尼标定方法，特别涉及一种基于工作模态分析的铣削过程阻尼标定方法。

背景技术：

切削过程中振动的能量被两种形式的阻尼吸收，一种是由工件、刀具-刀柄和机床系统产生的结构阻尼，另一种是由刀具后刀面与工件之间挤压与摩擦作用而产生的过程阻尼。在稳定性叶瓣图绘制过程中，如果只考虑结构阻尼而忽略过程阻尼的影响，其稳定性预测结果在较低转速范围内与实验观测值偏差很大。因此，对过程阻尼的识别具有重要理论及工程应用价值。

文献“l.tunc,e.budak,identificationandmodelingofprocessdampinginmilling,transactionsofasmejournalofmanufacturingscienceandengineering135(2013)021001.”公开了一种适用于铣削过程阻尼的标定方法。该方法首先利用颤振实验，通过在转速固定的情况下逐渐增大切深来获取稳定性极限切深值；然后根据稳定性求解的逆过程，从获得的极限切深中反推出过程阻尼系数。但这种方法只适用于稳定切削状态与不稳定切削状态界限清晰可辨的情况。

上述文献的典型特点是：在铣削过程阻尼的识别过程中，如果在稳定性极限切深附近不能清楚判定该切削状态是稳定切削还是非稳定切削，那么阻尼识别精度就会有所降低。

技术实现要素：

为了克服现有铣削过程阻尼标定方法实用性差的不足，本发明提供一种基于工作模态分析的铣削过程阻尼标定方法。该方法首先通过标准力锤冲击实验获得刀具-机床系统模态参数；然后进行无颤振-稳定铣削实验，利用位移传感器测量出铣刀刀杆x、y方向上、下两个位置处的振动信号；求出每个频率点下所测振动信号的功率谱密度矩阵，将该矩阵进行奇异值分解，然后剔除周期性切削力所对应的特征值，筛选出主导模态附近、mac值接近1、随机力激励所对应的较大特征值进行傅里叶逆变换，对逆变换结果的峰值和谷值应用指数衰减法来提取平均阻尼系数。得到平均阻尼系数之后，为表征瞬态下的过程阻尼，先利用切削力和刀尖频响得到刀尖位置振动幅值，再将其连同平均过程阻尼系数一起带入x、y两个方向上建立的能量平衡关系式中，同时求得径向和切向的犁切力系数；此时，利用所标定出的犁切力系数来表示瞬态下动态犁切力，该动态犁切力即为过程阻尼产生的本质原因。本发明方法能适用于绝对稳定切削与绝对不稳定切削之间界限不明确的情况，通过少量稳定铣削实验即可标定出过程阻尼系数，实用性好。

本发明解决其技术问题所采用的技术方案：一种基于工作模态分析的铣削过程阻尼标定方法，其特点是包括以下步骤：

步骤一、将加速度计贴在刀尖、刀杆下端和刀杆上端处，力锤在刀尖处进行敲击，从获得的频率响应曲线中提取x和y方向的主导模态参数mx，ζsr,x，ωn,x，和my，ζsr,y，ωn,y，其中mx和my是刀尖处的模态质量，ζsr,x和ζsr，y是刀尖处的结构阻尼比，ωn,x和ωn,y是主导模态频率，和是刀杆下端和刀杆上端处的模态振型向量。

步骤二、进行无颤振-稳定铣削实验，利用位移传感器测量铣刀刀杆x、y方向上刀杆下端和刀杆上端处的振动信号x(t)和y(t)，其中x(t)是x方向上刀杆下端和刀杆上端处振动信号时域值组成的向量，y(t)是y方向上刀杆下端和刀杆上端处振动信号时域值组成的向量。

步骤三、将时域中表示的振动信号x(t)进行快速傅里叶变换得到x(jω)，x(jω)是刀尖到奈奎斯特频率之间各个频率点x方向上刀杆下端和刀杆上端处的振动信号频域值组成的向量。

步骤四、通过下式计算振动信号的功率谱密度矩阵：

sxx(jω)＝x(jω)·x^h(jω)

式中，sxx(jω)是x方向上刀杆下端和刀杆上端处振动信号的功率谱密度矩阵，h表示对矩阵进行复共轭转置操作。

步骤五、对矩阵sxx(jω)进行奇异值分解得到特征值σx1，σx2和特征向量ux1，uy1，其中σx1是奇异值分解后得到的较大的特征值，σx2是奇异值分解后得到的较小的特征值，ux1是奇异值分解后σx1对应的特征向量，ux2是奇异值分解后σx2对应的特征向量。

步骤六、通过下式计算向量和ux1之间的mac值：

式中，mac表示特征值ux1与模态振型之间的相似程度。

步骤七、在主导模态附近选择mac值接近1、随机激励力所对应的σx1进行傅里叶逆变换，剔除周期性切削力所对应的σx1。

步骤八、对步骤七中傅里叶逆变换结果峰值与谷值的绝对值取对数运算后，进行线性拟合，得到拟合直线的斜率s，s表示拟合所得到直线的斜率。

步骤九、通过下式计算x方向上产生的总阻尼比：

式中，ζx是x方向上的总阻尼比。

步骤十、通过下式计算x方向上产生的平均过程阻尼系数：

cpr,x＝2mxωn,x(ζx-ζsr,x)

式中，cpr,x是x方向上的平均过程阻尼系数。

步骤十一、重复步骤三至步骤十，计算y方向上产生的平均过程阻尼系数cpr,y，cpr,y是y方向上的平均过程阻尼系数。

步骤十二、利用切削力和刀尖位置频响得到铣削过程中刀尖位置振动幅值x0，y0和振动频率f0，其中x0是x方向振动幅值，y0是y方向振动幅值，f0是振动频率。

步骤十三、通过下式计算x、y方向上刀尖位置振动速度：

式中，是x方向振动速度，是y方向振动速度。

步骤十四、通过下式计算动态犁切体积：

式中，vdy,ij是铣刀第i个刀齿上第j个单元参与切削时的动态犁切体积，v是切削线速度，lw是刀尖磨损带长度，θij是铣刀第i个刀齿上第j个单元参与切削时的切触角，zij是铣刀第i个刀齿上第j个单元的轴向高度。

步骤十五、通过联立下式求解切向和径向犁切力系数：

式中，kp,t是切向犁切力系数，kp,r是径向犁切力系数，t是主轴旋转周期，n是铣刀齿数，nel是铣刀轴向单元的数目。

本发明的有益效果是：该方法首先通过标准力锤冲击实验获得刀具-机床系统模态参数；然后进行无颤振-稳定铣削实验，利用位移传感器测量出铣刀刀杆x、y方向上、下两个位置处的振动信号；求出每个频率点下所测振动信号的功率谱密度矩阵，将该矩阵进行奇异值分解，然后剔除周期性切削力所对应的特征值，筛选出主导模态附近、mac值接近1、随机力激励所对应的较大特征值进行傅里叶逆变换，对逆变换结果的峰值和谷值应用指数衰减法来提取平均阻尼系数。得到平均阻尼系数之后，为表征瞬态下的过程阻尼，先利用切削力和刀尖频响得到刀尖位置振动幅值，再将其连同平均过程阻尼系数一起带入x、y两个方向上建立的能量平衡关系式中，同时求得径向和切向的犁切力系数；此时，利用所标定出的犁切力系数来表示瞬态下动态犁切力，该动态犁切力即为过程阻尼产生的本质原因。本发明方法能适用于绝对稳定切削与绝对不稳定切削之间界限不明确的情况，通过少量稳定铣削实验即可标定出过程阻尼系数，实用性好。

下面结合附图和具体实施方式对本发明作详细说明。

附图说明

图1是本发明基于工作模态分析的铣削过程阻尼标定方法所涉及铣刀的示意图。

图2是本发明方法实施例x方向特征值σx1。

图3是本发明方法实施例x方向特征值σx1傅里叶逆变换结果。

图4是本发明方法实施例y方向特征值σy1。

图5是本发明方法实施例y方向特征值σy1傅里叶逆变换结果。

图6是本发明实施例预测稳定性叶瓣图与实测颤振数据的对比图。

具体实施方式

参照图1-6。本发明基于工作模态分析的铣削过程阻尼标定方法具体步骤如下：

步骤一、实验选定半径为6mm、螺旋角为35°、法向前角为10°，后角为5°，刀尖圆弧半径20μm，齿数为4的硬质合金立铣刀在三坐标立式加工中心对铝合金7050-t7451进行顺铣切削。将加速度计贴在刀尖0、刀杆下端1和刀杆上端2处，力锤在刀尖0处进行敲击，从获得的频率响应曲线中提取x和y方向的主导模态参数mx＝0.3193kg，ζsr,x＝1.79％，ωn,x＝1060hz，和my＝0.3971kg，ζsr,y＝2.35％，ωn,y＝1041hz，其中mx和my是刀尖0处的模态质量，ζsr,x和ζsr,y是刀尖0处的结构阻尼比，ωn,x和ωn,y是主导模态频率，和是刀杆下端1和刀杆上端2处的模态振型向量。

步骤二、进行无颤振-稳定铣削实验，利用位移传感器测量铣刀刀杆x、y方向上刀杆下端1和刀杆上端2处的振动信号x(t)和y(t)，其中x(t)是x方向上刀杆下端1和刀杆上端2处振动信号时域值组成的向量，y(t)是y方向上刀杆下端1和刀杆上端2处振动信号时域值组成的向量。

步骤三、将时域中表示的振动信号x(t)进行快速傅里叶变换得到x(jω)，x(jω)是0到奈奎斯特频率之间各个频率点x方向上刀杆下端1和刀杆上端2处的振动信号频域值组成的向量。

步骤四、通过下式计算振动信号的功率谱密度矩阵：

sxx(jω)＝x(jω)·x^h(jω)

式中，sxx(jω)是x方向上刀杆下端1和刀杆上端2处振动信号的功率谱密度矩阵，h表示对矩阵进行复共轭转置操作。

步骤五、对矩阵sxx(jω)进行奇异值分解得到特征值σx1，σx2和特征向量ux1，uy1，其中σx1是奇异值分解后得到的较大的特征值，σx2是奇异值分解后得到的较小的特征值，ux1是奇异值分解后σx1对应的特征向量，ux2是奇异值分解后σx2对应的特征向量。

步骤六、通过下式计算向量和ux1之间的mac值：

式中，mac表示特征值ux1与模态振型之间的相似程度。

步骤七、在主导模态附近选择mac值接近1、随机激励力所对应的σx1进行傅里叶逆变换，剔除周期性切削力所对应的σx1。

步骤八、对步骤七中傅里叶逆变换结果峰值与谷值的绝对值取对数运算后，进行线性拟合，得到拟合直线的斜率s＝-0.2279，s表示拟合所得到直线的斜率。

步骤九、通过下式计算x方向上产生的总阻尼比：

式中，ζx是x方向上的总阻尼比。

步骤十、通过下式计算x方向上产生的平均过程阻尼系数：

cpr,x＝2mxωn,x(ζx-ζsr,x)

＝2×0.3193×2×π×1060×(7.24％-1.79％)

＝231.8n/(m/s)

式中，cpr,x是x方向上的平均过程阻尼系数。

步骤十一、重复步骤三至步骤十，计算y方向上产生的平均过程阻尼系数cpr,y＝407.8n/(m/s)，cpr,y是y方向上的平均过程阻尼系数。

步骤十二、利用切削力和刀尖位置频响得到铣削过程中刀尖位置振动幅值x0＝49μm，y0＝64μm和振动频率f0＝46.67hz，其中x0是x方向振动幅值，y0是y方向振动幅值，f0是振动频率。

步骤十三、通过下式计算x、y方向上刀尖位置振动速度：

式中，是x方向振动速度，是y方向振动速度。

步骤十四、通过下式计算动态犁切体积：

式中，vdy,ij是铣刀第i个刀齿上第j个单元参与切削时的动态犁切体积，v是切削线速度，lw是刀尖磨损带长度，θij是铣刀第i个刀齿上第j个单元参与切削时的切触角，zij是铣刀第i个刀齿上第j个单元的轴向高度。

步骤十五、通过联立下式求解切向和径向犁切力系数：

经联立求解得kp,r＝40059n/mm³，kp,t＝14386n/mm³，式中kp,t是切向犁切力系数，kp,r是径向犁切力系数，t是主轴旋转周期，n是铣刀齿数，nel是铣刀轴向单元的数目。

通过图6可以看出，采用实施例中所标定犁切力系数绘制的稳定性叶瓣图在较低的转速范围内能很好的与实测颤振结果相匹配，证明了本发明标定方法的准确性。