微细铣削加工颤振系统动力学分析及铣削过程仿真
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
微细铣削加工在制造业中有着较强的竞争力,并且随着科学技术水平的提高和社会的不断进步,民用和军工等领域对复杂微小尺寸光机电元器件需求日益增大,微细铣削加工在这些领域扮演着十分重要的角色。但在微细铣削过程中常常有切削颤振现象发生,其极大地影响机床加工性能、刀具的使用寿命和元件的加工质量,甚至会导致噪声等环境污染。本文针对铣削加工中主要的颤振形式—再生型切削颤振进行了分析,结合铣削加工参数(N,r,aε,kt,kr)和系统动力学参数(kx,ky,(?)x,(?)y,ωnx,ωny)建立了两自由度铣削系统动力学模型,利用Matlab软件得到切削过程的稳定性曲线,根据稳定性曲线预测加工过程的极限切削条件,从而选择合理的切削加工参数。论文着重分析了系统动力学参数变化对稳定性极限切深的影响规律,并得知各参数的变化不仅使稳定性曲线在纵坐标方向上有偏移,在横坐标上也有变化,此外又分析了各因素对最小极限切深的影响,从中得到刀具齿数的变化对最小极限切深的影响最为敏感的结论。同时论文在所建立的两自由度动力学模型基础上,利用稳定性曲线中的切削系统参数,计算出了X、Y方向上的平均铣削力。
     另外根据上述稳定性曲线,选择一组低于临界颤振条件的切削参数,以此切削参数为基础,对铣削过程进行有限元仿真和切削实验,然后根据二者的切削过程是否发生颤振及二者所得铣削力数据的比对结果,验证动力学模型的合理性和仿真模型的可靠性。对于仿真,用Deform-3D软件分析刀具在接触工件、切入工件刀具直径1/4处、切入工件刀具直径3/4处、完全切入工件以及刀具进入稳态切削等状态的应力和应变,据此分析了工件内部的弹性流动和塑性流动以及当工件发生大幅度塑性变形时,材料从工件表面发生去除的情况。根据Deform-3D软件分析结果,论文还研究了铣削过程工件所受到的作用力,得知在上述的五个状态中铣削力逐渐增加,其中在第三个状态切削力显著增大,达到五个切削状态的最大值,而后又降低,说明在切削过程存在弹性积累,使得刀具的切削负载增大,作用力急剧上升。
     对于实验考证,论文利用立扬MVL1300型立式加工中心对铣削加工过程中X、Y、Z三个方向的切削力进行了测试,将实验得到的切削力数据和软件仿真分析结果进行了对比分析,比对结果表明:仿真预测的铣削力数据与实验结果存在一定误差,但总体趋势上相符,基本验证了仿真模型的可靠性,同时根据实验和仿真中切削过程是否发生颤振验证了动力学模型的合理性;在文章最后将实验结果与理论计算的平均铣削力进行了比对,所得数据虽然存在误差,但都在一个数量级上,也验证了理论模型的合理性。
     论文工作得到科技部国际合作与交流专项项目(2010DFA72000)和教育部高等学校博士学科点专项基金(200801831024)的资助。
Micro-milling has Strong competitiveness in the manufacturing. With the improvement of science and technology and the progress of society, the demand for complex Optoelectronic components with micro-size is growing no matter in civil use or in national defense industry. Thus micro-milling is playing a very important role in the milling process. However, the chatter which accompanies the milling process often affects the processing properties of machine tools, life of tool and machining quality of workpiece, even produce the noise pollution. Based on the above mentioned, this paper analyzed the main chatter form-regenerative chatter in milling process. Meanwhile the paper established the cutting process dynamic model for the two degrees of freedom based on the cutting parameters of the milling process (N, r, ae, kt, kr) and the dynamics parameters of the chatter system(kx, ky, (?)x, (?)y,ωnx,ωny),then the stable curve of the cutting process was obtained using MATLAB software. Thus the appropriate cutting parameters were selected according to the limit cutting conditions predicted by stable curve. This paper mainly analyzed the laws that how the changes of the system dynamics parameters effected the stability limit depth of cut. And then got to know the factors not only made the stable curve changes in the vertical coordinate direction; there were also changes in the horizontal axis. And the paper also analyzed the factors'effect on the minimum limit depth of cut, and the result showed that the change of the number of teeth for tools is most sensitive to the minimum limit depth of cut. At the same time, using the two degrees of freedom dynamic model and the cutting system parameters, it calculated the average milling force in X, Y directions.
     In addition, chose a group of cutting parameters which is below the critical cutting conditions according to the stability curve. Based on the cutting parameters, Did the finite element simulation and cutting experiments in cutting process and contrast analysis results.then tested the rationality of the dynamics model and the reliability of the simulation model. Using Deform-3D,got the diagram of the stress and strain in five states, when the cutter was getting ready to cut into the workpiece, cut into the workpiece at 1/4 diameter of the cutter, cut into the workpiece at 3/4 diameter of the cutter, exactly completely cut into the workpiece and the cutter got into stable cutting. Then analyzed the elastic flow、the plastic flow in the workpiece and the macroscopic situation that the material removed from the surface of the workpiece, according to the diagram and the situation that when there is a large number of plastic deformation in the workpiece. And it also analyzed the force on the workpiece during the milling process using the Deform-3D software. In the five above states, the cutting force showed a gradual increase. But in the third state, the cutting force has a sharp increase which was the largest one in the five states. And then the force decreased, indicating that there is a flexible accumulation process during the cutting process. This process made the load of the cutter became larger and the force increased dramatically.
     As to the experiments, this paper tested the cutting force on the X, Y, Z three directions in the machining centre using LiYang MVL1300. The data from the Deform software was compared with experimental cutting force data. Although there was deviation among the comparison, they basically correspond with each other in the trend. This proved that the simulation model was reliable. The curve is also reliable according to the experiment and the simulation (flutter didn't occur when it below the critical cutting conditions). When we compared the experimental data with the theoretic average milling force, we got to know that there was deviation, but they were in the same magnitude. This also proved the rationality of the theory.
     This paper had gotten help from the special project of ministry of science and technology for international cooperation and exchanges(2010DFA72000) and special funded projects of the doctoral program of higher education for ministry of education (200801831024)
引文
[1]国家自然科学基金委员会.学科发展战略研究报告[c].机械与制造科学,2006.
    [2]于骏一,吴博达.机械加工振动的诊断、识别与控制[M].北京:清华大学出版社,1994.
    [3]星铁太郎著机械加工颤振的分析与对策[M].上海:上海科学技术出版社,1984年.
    [4]杨橚 唐恒龄 廖伯瑜 机床动力学(Ⅰ)[M].北京:机械工业出版社,1983年.
    [5]杨橚 唐恒龄 廖伯瑜 机床动力学(Ⅱ)[M].北京:机械工业出版社,1983年.
    [6]Colding, Bertil N Intelligent selection of machining parameter for metal cutting operations [J] Robotics and Computer-Integrated Manufacturing,1992,9(4-5):407-412,
    [7]于骏一,邹青著。机械制造技术基础[M].北京:机械工业出版社,2004.1.
    [8]Arnold R N.The mechanism of tool vibration in the cutting of steals [J]. ARCHIVE:Proceed ing-s of the Institution of Mechanical Engineers,1946,154:261-284.
    [9]Sission TR, Kegg RL.Explanation of low speed chatter effects [J]. ASME J.Ind,1969,97(4):9 51-958.
    [10]Minis I, Yanushevsky R. New theoretical approach for the prediction of machine tool chatter in milling[J] J.of ASME,1993,115(1):1-8.
    [11]Heisel U,Milberg J. Vibrations and Surface Generation in Slab Milling[J]. Annals of the CIRP, 1994,43(1):337-340.
    [12]Elbestawi M A. Modeling Machining Dynamics Including Damping in the Tool-Workpiece in-terface[J]. ASME Prod Eng Div Publ PED,1991,54:253-263.
    [13]Baradie E I. Statistical analysis of the dynamic cutting coefficient and machine tool stability [J]. Journal of engineering for industry,1995,115(2):205-215.
    [14]Altintas Y,Budak E. Analytical prediction of stability lobes in milling [J]. CIRP Annals-Manuf acturing Technology,1995,44(1):357-362.
    [15]Altintas Y, Shamoto E, Lee P..Analytical prediction of stability lobes in ball end milling [J]. Journal of manufacturing science and engineering-transactions of the ASME,1999,121(4):586-592.
    [16]Jensen S A. Stability analysis in face milling operation,(part Ⅰ):theory of stability lobe predict-ion[J]. Journal of Manufacturing Science and Engineering,1999,121(4):600-605.
    [17]Jensen S A,Shin Y C. Stability analysis in face milling operation,part Ⅱ:experimental validation and influence factors[J]. Journal of Manufacturing Science and Engineering,1999,121(4):606-614.
    [18]Faassen R P H, Nijmeijer H,Oosterling J A J. Prediction of regenerative chatter by modelling and analysis of milling high-speed [J].International Journal of Machine Tools&Manufacture,2 003,43:1437-1446.
    [19]Solis E, Peres C R, Jimenez J E,Alique J R. A new analytical-experimental method for the identification of stability lobes in high-speed milling [J]. International Journal of Machine Too Is and Manufacture,2004,44(15):1591-1597.
    [20]Shamoto E,Akazawa K. Analytical prediction of chatter stability in ball end milling with tool inclination[J]. CIRP Annals-Manufacturing Technology,2009,58:351-354.
    [21]Altintas Y, Stepan G Chatter stability of milling in frequency and discrete time domain [J].C-IRP Journal of Manufacturing Science and Technology,2008, 1(1):35-44.
    [22]Ismail F, Soliman E. A new method for the identification of stability lobes in maching [J]. International Journal of Machine Tools and Manufacture,1997,37(6):763-774.
    [23]Shi H M, obias S A. Theory of finite amplitude machine tool instability [J]. International Jou-rnal of Machine Tool Design and Research,1984,24(1):45-69.
    [24]Mohammad Malekian,imon S Park,Martin B G Jun.Modeling of dynamic micro-milling cutting forces[J]. International Journal of Machine Tools and Manufacture,2009,49(7-8):586-598.
    [25]Eynian M, Altintas Y. Analytical Chatter Stability of Milling With Rotating Cutter Dynamics at Process Damping Speeds [J].Journal of Manufacturing Science and Engineering,2010,132 /021012:1-14.
    [26]Zhao M X, Balachandran B. Dynamics and stability of milling process [J].International Journa 1 of Solids and Structures,2001,38:2233-2248.
    [27]Miao J C,Chen G L, Lai X M. Review of dynamic issues in micro-end-milling [J].Int J Adv-Manuf Technol,2007,31:897-904.
    [28]Msaddek E 1 Bechir, Bouaziz Zoubeir, Baili Maher,Gilles Dessein Modeling and simulation of high-speed milling centers dynamics [J].Int J Adv Manuf Technol,2010,53(9-12):877-888.
    [29]Altintas Y,Weck M. Chatter Stability of Metal Cutting and Grinding [J].CIRP Annals-Manufac t-uring Technology,2004,53(2):619-642.
    [30]张廷贤.动态铣削过程的稳定性分析[N].天津大学学报,1985,2:84-92.
    [31]张廷贤.切削加工稳定图[N].天津大学学报,1986,3:1-7.
    [32]张廷贤.铣削力建模与仿真[N].天津大学学报,1991,3:73-80.
    [33]阎兵,徐安平.球头铣刀铣削过程的物理建模与仿真[N].天津职业技术师范学院学报,2001.12,11(4):1-7.
    [34]徐安平,张大卫.基于再生颤振的端铣动态铣削过程建模与仿真[N].河北工业大学学报,1998,1:63-70.
    [35]袁伟军.铣削加工过程稳定性计算机仿真及实验研究[D].北京:北京工业大学硕士学位论文,2003.5.
    [36]李志刚.铣削加工过程的稳定性仿真及有限元研究[D].四川:西南交通大学硕士学位论文,2007.4.
    [37]李志刚,刘桓龙,王国志.铣削加工过程的动态稳定性极限分析[N].金属加工机械,2007,1(34):45-47.
    [38]孙孟琴.球头铣刀铣削加工过程的动力学建模及仿真[D].浙江:浙江工业大学硕士学位论文,2-006.4
    [39]李忠伟.铣削系统动力学预测和稳定性分析[D].上海:上海交通大学硕士学位论文,2008.6
    [40]张孝峰.球头铣削动力学建模及仿真[D].山东:山东大学硕士学位论文,2009.4.
    [41]赵宏伟,王晓军,于骏一.再生型机床切削颤振系统稳定性极限预测[J].吉林:精密制造与自动化,2003,12:62-64+67-4.
    [42]赵宏伟.再生型机床切削颤振系统稳定性极限预测[D].吉林:吉林大学硕士学位论文,2003.2.
    [43]王晓军.车削加工系统稳定性极限预测的研究[D].吉林:吉林大学博士学位论文,2005.12.
    [44]成志清,金瑞琪.切削颤振研究的历史、现状及其发展动向[N].江苏工学院学报,1992,13(2):20-26.
    [45]付连宇,于骏一,王文才.变速切削系统稳定性的实验研究[J].机械制造,1998,11:21-22.
    [46]李沪曾,于信汇,张国红,呙清强.铣削振动的计算机仿真[N].振动工程学报,2001,14(3):48-53.
    [47]李彬,邓建新,张松.高速切削加工稳定性研究的现状和发展趋势[J].制造技术与机床,2008,4:33-36.
    [48]王跃辉,王民.金属切削过程颤振控制技术的研究进展[N].机械工程学报,2010.4,46(7):166-174.
    [49]张莉,李升军Deform在金属塑性成形中的应用[M].北京:机械工业出版社,2009.
    [50]陈明祥.弹塑性力学[M].北京:科学出版社,2007.
    [51]刘谨.机床动力学(一)[J].机械制造,1987,12:40-43.
    [52]刘谨.机床动力学(二)[J].机械制造,1988,1:43-46.
    [53]刘谨.机床动力学(三)[J].机械制造,1988,2:44-47.
    [54]刘谨.机床动力学(四)[J].机械制造,1988,3:34-37.
    [55]周敏,张之敬.微小型正交车铣单齿圆周刃理论切削力研究[J].中国机械工程,2009,(13):1527-1532.
    [56]刘胤,李志忠,张智海,郑力.球头铣刀加工仿真中几何边界参数的提取[J].清华大学学报(自然科学版),2005,45(2):254-257.
    [57]王苏东.基于DEFORM-3D的钛合金切削过程有限元仿真[J]装备制造技术,2009,(12):30-40.
    [58]董兆伟,张以都,刘胜永.铣削过程中残余应力仿真分析[J]北京航空航天大学学报,2005,31(7):762-765.
    [59]焦辉,汤文成.薄壁件铣削过程仿真与分析[J]中国制造业信息化,2007,19:121-125.
    [60]包杰,李亮,何宁,曹自洋,黄亮.国外微细铣削研究综述[J].机械科学与技术,2009,28(8):1018-1022.