切削过程中薄壁件动力学参数的预测方法
【专利摘要】本发明公开了一种切削过程中薄壁件动力学参数的预测方法，用于解决现有薄壁件动力学参数的预测方法实用性差的技术问题。技术方案是通过有限元方法计算出未加工的初始工件的固有频率及模态振型；然后利用软化材料属性的方法软化工件质量及刚度矩阵，模拟材料去除过程，得到工件质量及刚度修改矩阵；最后通过结构动力修改方法计算切削过程中的工件固有频率及多点模态振型。由于利用结构动力修改方法，实现了计算加工过程中的工件的固有频率及多点模态振型，适用于带有平面和曲面薄壁件的工件动力学参数的计算，能够应用于多点接触的切削动力学模型，进行切削过程的稳定性预测，实用性强。
【专利说明】
切削过程中薄壁件动力学参数的预测方法
技术领域
[0001] 本发明涉及一种薄壁件动力学参数的预测方法，特别涉及一种切削过程中薄壁件 动力学参数的预测方法。
【背景技术】
[0002] 切削加工过程中的颤振会极大地影响零件加工效率和加工质量，加速刀具磨损和 刀具破损，因此铣削加工中应尽量避免颤振的发生。通过稳定性叶瓣图选取合适的切削加 工参数是避免切削颤振的重要途径之一。研究表明，在薄壁件切削过程中，工件的动力学特 性对切削稳定性叶瓣图预测准确性有很大的影响，故工件的动力学参数的准确获取对切削 稳定性叶瓣图的预测至关重要。工件的动力学参数一般用频响函数或者固有频率及模态振 型表示，通常可通过力锤冲击实验获得。由于在薄壁件切削工程中，去除材料占整个工件材 料的比例较大，工件的动力学参数将随着材料去除过程而显著变化。此时要准确预测切削 稳定性需要准确获得加工过程中的工件动力学参数。但若每一个时刻的工件动力学参数都 通过力锤冲击实验测量，将花费大量的时间。为此，研究者提出了计算工件的动力学参数的 方法。
[0003] 文南犬 1 "Q. Song，Z. Liu，Y.Wan，G. Ju，J. Shi，Application of Sherman-Morrison-Woodbury formulas in instantaneous dynamic of peripheral milling for thin-walled component，International Journal of Mechanical Sciences 96-97(2015)79-90."公开了一种计算工件的动力学参数的方法。考虑了工件动力学参数因材料去除的变化 及其在不同刀具位置处的变化，使用Sherman-Morrison-Woodbury公式计算加工过程中的 工件某点的频响函数，并将其代入最常用的单点刀具工件动力学模型中预测薄壁件周铣稳 定性。
[0004] 文献 2"E.Budak，L.T.Tunc，S.Alan，H.N.0zgven，Predictionofworkpiece dynamics and its effects on chatter stability in milling,CIRP Annals-Manufacturing Technology 61 (1)(2012)339-342·"公开了一种计算工件的动力学参数的 方法。基于结构动力修改方法计算了加工过程中的工件某点的频响函数，并用于五轴铣削 稳定性的计算。
[0005] 现有的计算工件的动力学参数的方法均是针对计算加工过程中的工件某单点频 响函数的方法，不能计算加工过程中的工件的固有频率及多点模态振型，不适用于多点接 触的切削动力学模型。
【发明内容】
[0006] 为了克服现有薄壁件动力学参数的预测方法实用性差的不足，本发明提供一种切 削过程中薄壁件动力学参数的预测方法。该方法通过有限元方法计算出未加工的初始工件 的固有频率及模态振型；然后利用软化材料属性的方法软化工件质量及刚度矩阵，模拟材 料去除过程，得到工件质量及刚度修改矩阵;最后通过结构动力修改方法计算切削过程中 的工件固有频率及多点模态振型。由于利用结构动力修改方法，实现了计算加工过程中的 工件的固有频率及多点模态振型，适用于带有平面和曲面薄壁件的工件动力学参数的计 算，可应用于多点接触的切削动力学模型，进行切削过程的稳定性预测，实用性强。
[0007] 本发明解决其技术问题所采用的技术方案:一种切削过程中薄壁件动力学参数的 预测方法，其特点是包括以下步骤：
[0008] 步骤一、建立未加工的初始工件的有限元模型，得到未加工的初始工件的质量矩 阵胸,0和刚度矩阵&,〇,并对有限元模型进行计算模态分析，得到初始工件的固有频率矩阵 ω w, 〇和模态振型矩阵Uw, 〇;
[0009] 步骤二、对刀具轨迹进行离散，得到nGS个刀具位置点；
[0010] 步骤三、根据刀具轨迹计算刀具从初始刀具位置点运动到第m刀具位置点时的刀 具扫掠轮廓，并判断这个过程中被去除材料包含的单元；
[0011] 步骤四、通过对被去除材料包含的单元乘以ΚΓ6对被去除材料进行软化，并得到质 量和刚度变化矩阵，即AMw,o和Δ&,ο;
[0012] 步骤五、根据步骤一和步骤四中的〇^,〇、1^()、八|^()和八&,()，将刀具运动到第111刀 具位置点时的加工过程中的工件的有限元模型写为：
[0013] (1 + U^i0AMw mtJ w>0 )f w 0(/) + (〇w0_ + U^0AKW mUw o )Γw o(t) = 0
[0014] 其中Γ w,o(t)为模态坐标向量，I和0分别为单位矩阵和零矩阵；
[0015] 步骤六、对步骤五中的有限元模型进行计算模态分析，得到加工过程中的工件的 固有频率矩阵cow,m和模态变换矩阵U m;
[0016] 步骤七、根据步骤一中的Uw,o,得到加工过程中的工件的模态振型矩阵Uw,m = Uw, 〇Um 〇
[0017] 本发明的有益效果是:该方法通过有限元方法计算出未加工的初始工件的固有频 率及模态振型；然后利用软化材料属性的方法软化工件质量及刚度矩阵，模拟材料去除过 程，得到工件质量及刚度修改矩阵;最后通过结构动力修改方法计算切削过程中的工件固 有频率及多点模态振型。由于利用结构动力修改方法，实现了计算加工过程中的工件的固 有频率及多点模态振型，适用于带有平面和曲面薄壁件的工件动力学参数的计算，能够应 用于多点接触的切削动力学模型，进行切削过程的稳定性预测，实用性强。
[0018] 下面结合附图和【具体实施方式】对本发明作详细说明。
【附图说明】
[0019] 图1是本发明方法实施例1中平板薄壁件的示意图。
[0020]图2是本发明方法实施例1中计算的工件的第二阶固有频率随加工过程的变化图。
[0021] 图3是本发明方法实施例1中计算的不同刀具位置处的工件-刀具接触区域的模态 振型。
[0022] 图4是本发明方法实施例2中带有曲面的薄壁件的示意图。
[0023] 图5是本发明方法实施例2中计算的工件的第一阶固有频率随加工过程的变化图。
[0024] 图6是本发明方法实施例2中计算的不同刀具位置处的工件-刀具接触区域的模态 振型。
[0025]图中，1、2、3、4分别表示实施例1中刀具运动到第1、10、19、29个刀具位置处的工 件-刀具接触区域的模态振型。5、6、7、8分别表示实施例2中刀具运动到第1、6、11、15个刀具 位置处的工件-刀具接触区域的模态振型。
【具体实施方式】
[0026]以下实施例参照图1-6。
[0027]实施例1:采用本发明进行平板薄壁件周铣过程中工件动力学参数的计算，其中平 板尺寸为115mm X 36mm X 3 · 5mm，径向切削深度为0 · 5mm，工件材料为铝合金6061-T6，工件通 过底部装夹在机床工作台上。
[0028] (1)建立未加工的初始工件的有限元模型，得到未加工的初始工件的质量矩阵Mw,o 和刚度矩阵&,0,并对有限元模型进行计算模态分析，得到初始工件的固有频率矩阵〇^,〇和 模态振型矩阵Uw, 〇;
[0029] (2)对刀具轨迹进行离散，得到29个刀具位置点；
[0030] (3)根据刀具轨迹计算刀具从初始刀具位置点运动到第m刀具位置点时的刀具扫 掠轮廓，并判断这个过程中被去除材料包含的单元；
[0031] (4)通过对被去除材料包含的单元乘以10-6对被去除材料进行软化，并得到质量和 刚度变化矩阵，即和Δ&,ο;
[0032] (5)根据步骤（1)和(4)中的ω w,〇、Uw, ο、Δ Mw,〇和Δ Kw,〇，将刀具运动到第m刀具位置 点时的加工过程中的工件的有限元模型写为：
[0033] (I + UWi0AMWmUWi() )f w 0(?) + (〇)w 〇 + U w 0AKw<mU w>0 )rwo (i) = 0
[0034] 其中rw,o(t)为模态坐标向量，I和0分别为单位矩阵和零矩阵；
[0035] (6)对步骤(5)中的有限元模型进行计算模态分析，得到加工过程中的工件的固有 频率矩阵co w,m和模态变换矩阵Um;
[0036] (7)根据步骤(1)中的Uw,Q，得到加工过程中的工件的模态振型矩阵Uw, m=Uw,0Um;
[0037] 通过上面的步骤，可以计算平板薄壁件的周铣过程中工件动力学参数，从图2、3可 以看出，本发明的本发明可以预测切削过程中工件的固有频率和多点模态振型，证明了方 法的有效性。
[0038] 实施例2:采用本发明进行带有曲面的薄壁件周铣过程中工件动力学参数的计算， 曲面方程为
[0039] x(u,v)=20+190u(1-u)2+40u2(1-u)+80u3 u,ve[0,l]
[0040] y(u,v) =5+80u2( l-u)+25u3
[0041] z(u, v) =43v
[0042] 该薄壁件通过对曲面的内外两个方向分别偏置0.3mm和2.2mm得到，径向切削深度 为0.3mm，工件材料为铝合金6061-T6，工件通过底部装夹在机床工作台上。
[0043] (1)建立未加工的初始工件的有限元模型，得到未加工的初始工件的质量矩阵Mw,o 和刚度矩阵&,〇,并对有限元模型进行计算模态分析，得到初始工件的固有频率矩阵〇^,〇和 模态振型矩阵Uw, 〇;
[0044] (2)对刀具轨迹进行离散，得到21个刀具位置点；
[0045] (3)根据刀具轨迹计算刀具从初始刀具位置点运动到第m刀具位置点时的刀具扫 掠轮廓，并判断这个过程中被去除材料包含的单元；
[0046] (4)通过对被去除材料包含的单元乘以10-6对被去除材料进行软化，并得到质量和 刚度变化矩阵，即和Δ&,ο;
[0047] (5)根据步骤（1)和(4)中的〇^,〇、1^()、八|%,()和八&,()，将刀具运动到第111刀具位置 点时的加工过程中的工件的有限元模型写为：
[0048] (I + U^0AMW,,"UW,0 )f WJ)(〇 + (?w,02 +U；,0AKw,",Uw,0)rw, 0(〇 = 0
[0049] 其中rw,o(t)为模态坐标向量，I和0分别为单位矩阵和零矩阵；
[0050] (6)对步骤(5)中的有限元模型进行计算模态分析，得到加工过程中的工件的固有 频率矩阵co w,m和模态变换矩阵Um;
[0051] (7)根据步骤(1)中的Uw,Q，得到加工过程中的工件的模态振型矩阵Uw, m=Uw,0Um;
[0052] 通过上面的步骤，可以计算带有曲面的薄壁件的周铣过程中工件动力学参数，从 图5、6可以看出，本发明的本发明可以预测带有曲面的薄壁件的切削过程中工件的固有频 率和多点模态振型，证明了方法的可适用于平面和曲面。
【主权项】
1. 一种切削过程中薄壁件动力学参数的预测方法，其特征在于包括W下步骤： 步骤一、建立未加工的初始工件的有限元模型，得到未加工的初始工件的质量矩阵Mw,0 和刚度矩阵&,〇,并对有限元模型进行计算模态分析，得到初始工件的固有频率矩阵ww,o和 模态振型矩阵Uw, 0; 步骤二、对刀具轨迹进行离散，得到nes个刀具位置点； 步骤Ξ、根据刀具轨迹计算刀具从初始刀具位置点运动到第m刀具位置点时的刀具扫 掠轮廓，并判断运个过程中被去除材料包含的单元； 步骤四、通过对被去除材料包含的单元乘Wl(T6对被去除材料进行软化，并得到质量和 刚度变化矩阵，即AMw,o和AKw,o; 步骤五、根据步骤一和步骤四中的ωw,o、Uw,o、ΔMw,o和ΔKw,o，将刀具运动到第m刀具位 置点时的加工过程中的工件的有限元模型写为：其中rw,o(t)为模态坐标向量，I和0分别为单位矩阵和零矩阵； 步骤六、对步骤五中的有限元模型进行计算模态分析，得到加工过程中的工件的固有 频率矩阵cow,m和模态变换矩阵Um; 步骤屯、根据步骤一中的Uw,0,得到加工过程中的工件的模态振型矩阵Uw,m = Uw,〇Um。
【文档编号】G06F17/50GK106096146SQ201610422362
【公开日】2016年11月9日
【申请日】2016年6月15日 公开号201610422362.9, CN 106096146 A, CN 106096146A, CN 201610422362, CN-A-106096146, CN106096146 A, CN106096146A, CN201610422362, CN201610422362.9
【发明人】张卫红, 杨昀, 刘辉, 万敏, 马颖超
【申请人】西北工业大学