CN107238480A - 基于工作模态分析的铣削过程阻尼标定方法 - Google Patents

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

Info

Publication number
CN107238480A
CN107238480A CN201710491648.7A CN201710491648A CN107238480A CN 107238480 A CN107238480 A CN 107238480A CN 201710491648 A CN201710491648 A CN 201710491648A CN 107238480 A CN107238480 A CN 107238480A
Authority
CN
China
Prior art keywords
mrow
msub
directions
munderover
mover
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710491648.7A
Other languages
English (en)
Other versions
CN107238480B (zh
Inventor
万敏
冯佳
张卫红
马颖超
杨昀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710491648.7A priority Critical patent/CN107238480B/zh
Publication of CN107238480A publication Critical patent/CN107238480A/zh
Application granted granted Critical
Publication of CN107238480B publication Critical patent/CN107238480B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00—Vibration-testing of structures; Shock-testing of structures

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automatic Control Of Machine Tools (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于工作模态分析的铣削过程阻尼标定方法，用于解决现有铣削过程阻尼标定方法实用性差的技术问题。技术方案是首先获取刀具‑机床系统模态参数；然后进行无颤振‑稳定铣削实验，测量出铣刀刀杆x、y方向上、下两个位置处的振动信号；求出每个频率点下所测振动信号的功率谱密度矩阵，将该矩阵进行奇异值分解，然后剔除周期性切削力所对应的特征值，对所对应的较大特征值进行傅里叶逆变换，对逆变换结果的峰值和谷值提取平均阻尼系数。得到平均阻尼系数之后，先利用切削力和刀尖频响得到刀尖位置振动幅值，再将其连同平均过程阻尼系数一起带入x、y两个方向上建立的能量平衡关系式中，标定出铣削过程阻尼系数，实用性好。

Description

基于工作模态分析的铣削过程阻尼标定方法
技术领域
本发明涉及一种铣削过程阻尼标定方法，特别涉及一种基于工作模态分析的铣削过程阻尼标定方法。
背景技术
切削过程中振动的能量被两种形式的阻尼吸收，一种是由工件、刀具-刀柄和机床系统产生的结构阻尼，另一种是由刀具后刀面与工件之间挤压与摩擦作用而产生的过程阻尼。在稳定性叶瓣图绘制过程中，如果只考虑结构阻尼而忽略过程阻尼的影响，其稳定性预测结果在较低转速范围内与实验观测值偏差很大。因此，对过程阻尼的识别具有重要理论及工程应用价值。
文献“L.Tunc,E.Budak,Identification and modeling of process damping inmilling,Transactions of ASME Journal of Manufacturing Science and Engineering135(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ω)·XH(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ω)·XH(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/mm3，Kp,t＝14386N/mm3，式中Kp,t是切向犁切力系数，Kp,r是径向犁切力系数，T是主轴旋转周期，N是铣刀齿数，Nel是铣刀轴向单元的数目。
通过图6可以看出，采用实施例中所标定犁切力系数绘制的稳定性叶瓣图在较低的转速范围内能很好的与实测颤振结果相匹配，证明了本发明标定方法的准确性。

Claims (1)

1.一种基于工作模态分析的铣削过程阻尼标定方法，其特征在于包括以下步骤：
步骤一、将加速度计贴在刀尖、刀杆下端和刀杆上端处，力锤在刀尖处进行敲击，从获得的频率响应曲线中提取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ω)·XH(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方向上产生的总阻尼比：
<mrow> <msub> <mi>&amp;zeta;</mi> <mi>x</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>s</mi> </mrow> <msqrt> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> </mrow>
式中，ζ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方向上刀尖位置振动速度：
<mrow> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mn>0</mn> </msub> <msub> <mi>x</mi> <mn>0</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mn>0</mn> </msub> <mi>t</mi> <mo>)</mo> </mrow> </mrow>
<mrow> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mn>0</mn> </msub> <msub> <mi>y</mi> <mn>0</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mn>0</mn> </msub> <mi>t</mi> <mo>)</mo> </mrow> </mrow>
式中，是x方向振动速度，是y方向振动速度；
步骤十四、通过下式计算动态犁切体积：
<mrow> <msub> <mi>V</mi> <mrow> <mi>d</mi> <mi>y</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msubsup> <mi>l</mi> <mi>w</mi> <mn>2</mn> </msubsup> <mrow> <mn>2</mn> <mi>v</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>sin&amp;theta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>cos&amp;theta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <msub> <mi>z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow>
式中，Vdy,ij是铣刀第i个刀齿上第j个单元参与切削时的动态犁切体积，v是切削线速度，lw是刀尖磨损带长度，θij是铣刀第i个刀齿上第j个单元参与切削时的切触角，zij是铣刀第i个刀齿上第j个单元的轴向高度；
步骤十五、通过联立下式求解切向和径向犁切力系数：
<mrow> <munderover> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>T</mi> </munderover> <mo>&amp;lsqb;</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mrow> <mi>e</mi> <mi>l</mi> </mrow> </msub> </munderover> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mi>cos</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mo>,</mo> <mi>r</mi> </mrow> </msub> <mi>sin</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>V</mi> <mrow> <mi>d</mi> <mi>y</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>=</mo> <munderover> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>T</mi> </munderover> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mi>p</mi> <mi>r</mi> <mo>,</mo> <mi>x</mi> </mrow> </msub> <msup> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> </mrow>
<mrow> <munderover> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>T</mi> </munderover> <mo>&amp;lsqb;</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mrow> <mi>e</mi> <mi>l</mi> </mrow> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mi>sin</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mo>,</mo> <mi>r</mi> </mrow> </msub> <mi>cos</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>V</mi> <mrow> <mi>d</mi> <mi>y</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>=</mo> <munderover> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>T</mi> </munderover> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mi>p</mi> <mi>r</mi> <mo>,</mo> <mi>y</mi> </mrow> </msub> <msup> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> </mrow>
式中，Kp,t是切向犁切力系数，Kp,r是径向犁切力系数，T是主轴旋转周期，N是铣刀齿数，Nel是铣刀轴向单元的数目。
CN201710491648.7A 2017-06-20 2017-06-20 基于工作模态分析的铣削过程阻尼标定方法 Active CN107238480B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710491648.7A CN107238480B (zh) 2017-06-20 2017-06-20 基于工作模态分析的铣削过程阻尼标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710491648.7A CN107238480B (zh) 2017-06-20 2017-06-20 基于工作模态分析的铣削过程阻尼标定方法

Publications (2)

Publication Number Publication Date
CN107238480A true CN107238480A (zh) 2017-10-10
CN107238480B CN107238480B (zh) 2019-02-12

Family

ID=59987918

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710491648.7A Active CN107238480B (zh) 2017-06-20 2017-06-20 基于工作模态分析的铣削过程阻尼标定方法

Country Status (1)

Country Link
CN (1) CN107238480B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732995A (zh) * 2018-04-25 2018-11-02 西北工业大学 铣削过程工件动力学参数的快速获取方法
CN108776430A (zh) * 2018-05-08 2018-11-09 苏州科技大学 一种滚珠丝杠进给驱动系统位置环增益优化取值方法
CN114492013A (zh) * 2022-01-07 2022-05-13 西北工业大学 一种考虑金属死区与材料回弹的微铣削过程阻尼建模方法
CN116698974A (zh) * 2023-06-16 2023-09-05 中南林业科技大学 基于自由振动曲线全采样的木基材料阻尼测试方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070088456A1 (en) * 2005-04-07 2007-04-19 University Of Florida Research Foundation, Inc. System and method for tool point prediction using multi-component receptance coupling substructure analysis
CN101122541A (zh) * 2007-08-03 2008-02-13 东方电气集团东方汽轮机有限公司 汽轮机叶片振动试验方法及装置
CN101539493A (zh) * 2009-04-28 2009-09-23 中国地震局工程力学研究所 用于结构损伤诊断的对称信号方法
WO2012125837A1 (en) * 2011-03-15 2012-09-20 Purdue Research Foundation Weak bond detection
CN103217284A (zh) * 2013-04-18 2013-07-24 西北工业大学 一种旋转机械模态阻尼系数的测量方法
CN103419089A (zh) * 2012-05-22 2013-12-04 湖南湖机国际机床制造有限公司 一种圆锯床整机动态工作性能优化设计方法
CN103939082A (zh) * 2014-04-29 2014-07-23 西北工业大学 一种基于主动阻尼法抑制钻杆粘滑振动的系统及方法
CN104850713A (zh) * 2015-05-28 2015-08-19 西北工业大学 机械结构随机振动动态应力高精度计算方法
CN105973621A (zh) * 2016-05-02 2016-09-28 国家电网公司 一种基于异常振动分析的gis机械故障诊断方法和系统
US9534654B2 (en) * 2013-03-27 2017-01-03 Zf Friedrichshafen Ag Vibration damper with a clamp-shaped attachment part

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070088456A1 (en) * 2005-04-07 2007-04-19 University Of Florida Research Foundation, Inc. System and method for tool point prediction using multi-component receptance coupling substructure analysis
CN101122541A (zh) * 2007-08-03 2008-02-13 东方电气集团东方汽轮机有限公司 汽轮机叶片振动试验方法及装置
CN101539493A (zh) * 2009-04-28 2009-09-23 中国地震局工程力学研究所 用于结构损伤诊断的对称信号方法
WO2012125837A1 (en) * 2011-03-15 2012-09-20 Purdue Research Foundation Weak bond detection
CN103419089A (zh) * 2012-05-22 2013-12-04 湖南湖机国际机床制造有限公司 一种圆锯床整机动态工作性能优化设计方法
US9534654B2 (en) * 2013-03-27 2017-01-03 Zf Friedrichshafen Ag Vibration damper with a clamp-shaped attachment part
CN103217284A (zh) * 2013-04-18 2013-07-24 西北工业大学 一种旋转机械模态阻尼系数的测量方法
CN103939082A (zh) * 2014-04-29 2014-07-23 西北工业大学 一种基于主动阻尼法抑制钻杆粘滑振动的系统及方法
CN104850713A (zh) * 2015-05-28 2015-08-19 西北工业大学 机械结构随机振动动态应力高精度计算方法
CN105973621A (zh) * 2016-05-02 2016-09-28 国家电网公司 一种基于异常振动分析的gis机械故障诊断方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
EROL TURKES: "A new process damping model for chatter vibration", 《MEASUREMENT》 *
万敏等: "铣削加工工艺力学机理研究", 《航空制造技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732995A (zh) * 2018-04-25 2018-11-02 西北工业大学 铣削过程工件动力学参数的快速获取方法
CN108732995B (zh) * 2018-04-25 2019-05-31 西北工业大学 铣削过程工件动力学参数的快速获取方法
CN108776430A (zh) * 2018-05-08 2018-11-09 苏州科技大学 一种滚珠丝杠进给驱动系统位置环增益优化取值方法
CN108776430B (zh) * 2018-05-08 2021-02-02 苏州科技大学 一种滚珠丝杠进给驱动系统位置环增益优化取值方法
CN114492013A (zh) * 2022-01-07 2022-05-13 西北工业大学 一种考虑金属死区与材料回弹的微铣削过程阻尼建模方法
CN116698974A (zh) * 2023-06-16 2023-09-05 中南林业科技大学 基于自由振动曲线全采样的木基材料阻尼测试方法

Also Published As

Publication number Publication date
CN107238480B (zh) 2019-02-12

Similar Documents

Publication Publication Date Title
CN107238480A (zh) 基于工作模态分析的铣削过程阻尼标定方法
Zatarain et al. Analysis of the influence of mill helix angle on chatter stability
CN106965032B (zh) 薄壁件铣削颤振抑制方法
Honeycutt et al. A new metric for automated stability identification in time domain milling simulation
CN104484506B (zh) 一种基于可靠性叶瓣图的车削加工颤振预测方法
CN103345200B (zh) 一种基于广义区间的切削加工颤振辨识方法
CN107457609A (zh) 基于刚度变化的铣削颤振抑制方法及铣削颤振优化系统
Gao et al. Experimental study of micro-milling mechanism and surface quality of a nickel-based single crystal superalloy
CN106934172A (zh) 一种碳纤维复合材料的多刃铣削去除率计算方法
Liang et al. Influence of chisel edge axial rake angle on the drilling performance of helical point micro-drill
Zhang et al. A novel finite element investigation of cutting force in orthogonal cutting considering plough mechanism with rounded edge tool
Namlu et al. Combined use of ultrasonic-assisted drilling and minimum quantity lubrication for drilling of NiTi shape memory alloy
CN102922010B (zh) 一种自润滑麻花钻及其加工方法
CN202726640U (zh) 一种用于检测机床切削能力的测试试件
CN107609311A (zh) 基于排屑力模型的深孔钻削深度优化方法
Gao et al. Experimental research on micro-milling force of a single-crystal nickel-based superalloy
Du et al. A new dynamic boring force calculation method using the analytical model of time-varying toolpath and chip fracture
Zedan et al. Effects of lubrication modes on part quality during drilling 6061-T6 aluminium alloy
Ibañez et al. Workpiece material influence on stability lobe diagram
Ding et al. Two dimensional vibration-assisted micro-milling: kinematics simulation, chip thickness computation and analysis
Wu et al. Study on the mechanism of AL2024-T351/Ti-6Al-4V laminated materials by ultrasonic vibration drilling
CN108747589B (zh) 切削过程中加工阻尼系数的快速准确识别方法
CN110516340A (zh) 基于铣削系统刚度变化的过程阻尼统一模型建模方法
CN110941914A (zh) 一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法
CN203636030U (zh) 一种可转位锥度钻孔刀具

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant