Stiffness modeling technique for five-axis milling process system with thin-walled parts
ZHAO Tong1, BIAN Pengxi1, WANG Yongfei2, ZHANG Yibo2
1. Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China; 2. Beijing Power Machinery Research Institute, Beijing 100074, China
Abstract:[Objective] The manufacturing industry is witnessing increasing precision requirements for thin-walled components, such as turbine blades and casings, in the field of aerospace engines. This is because these components suffer high material removal rates and poor machining performance, resulting in extreme changes in cutting forces during machining and the formation of chatter marks on the machined surface, leading to decreased surface quality. The relative vibration between the cutting edge and workpiece surface is mainly influenced by the milling forces and stiffness of various parts of the process system. Therefore, modeling of stiffness and analysis of such process systems are crucial for increasing the machining accuracy. [Methods] Based on the research on machine tool modeling using cutting tools, casings, and fixtures available domestically and internationally, simulation calculations and process system analyses were conducted under typical working conditions. A large amount of data was generated through hammering and vibration excitation experiments, and an eXtreme Gradient Boosting (XGBoost) algorithm was employed to establish a stiffness model for a cradle-type AC rotary table structure machine tool; furthermore, this model includes pose information. The stiffness model of the machine tool was integrated with the tool, workpiece, and fixture models obtained via the finite element method using a stiffness series formula and the spatial coordinate transformation method for developing the stiffness model of the process system. The stiffness model’s accuracy was increased through iterative updates. First, the harmonic response of the tool, casing, and fixture was calculated through the finite element simulation. Second, a suspended measurement device for measuring the dynamic stiffness of the machine tool was proposed using an electromagnetic exciter. The device offers the advantages of high excitation power, high measurement signal-to-noise ratio, and measurement ability in different machine tool poses. Subsequently, the dynamic stiffness values of the machine tool in different poses were measured using the exciter. After obtaining the stiffness models of the tool, fixture, and workpiece systems through the finite element method, the overall stiffness modeling of the process system was completed through the stiffness series theory and spatial coordinate transformation. Finally, the calibration and iteration mechanism of the stiffness model of the process system were established. An XGBoost algorithm was established to estimate the dynamic stiffness of the process system using milling force and surface topography information measured through milling experiments. An adaptive undersampling technique for discrete data with equal chord height was proposed to address the severe imbalance between the amounts of dynamic stiffness data obtained through the abovementioned model and vibration excitation test data. [Results] 1) After comparing the measurement results of the hammering experiment with simulation result, the maximum error of the dynamic stiffness obtained through the finite element simulation was <10% and the average error was <6%, indicating that relatively accurate tool and casing models can be obtained through this simulation. 2) After measuring the dynamic stiffness of the machine tool at different poses by using the electromagnetic exciter, the error rate of the machine tool stiffness model established using XGBoost algorithm meets the 3σ criterion. 3) The average error of the XGBoost algorithm for measuring the process system stiffness through milling force and surface topography is less than 13%, with a maximum value of no more than 18%. [Conclusions] The undersampled vibration excitation test data and stiffness data obtained through the cutting experiments were used for calibrating and iterating the stiffness model of the process system, ensuring that the error rate of the model meets the 3σ criterion.
[1] KAM M, DEM AI· RTAŞM. Experimental analysis of the effect of mechanical properties and microstructure on tool vibration and surface quality in dry turning of hardened AISI 4340 steels[J]. Surface Review and Letters, 2021, 28(6):2150041. [2] SANTO D R, MENCIK J M, TANG B, et al. On the dynamic behavior of rods and beams with nonlinear boundary stiffness[C]//Latin American Congress on Computational Methods in Engineering. Paris, France:Université de Technologie de Compiègne, 2018. [3] 孙孟琴,王立臣.球头铣刀动力学模型及其铣削加工稳定性的研究[J].机床与液压, 2012, 40(9):52-54. SUN M Q, WANG L C. Research on dynamic model and milling stability of ball end milling cutter[J]. Machine Tool&Hydraulics, 2012, 40(9):52-54.(in Chinese) [4] 邓聪颖,刘蕴,殷国富,等.基于响应面方法的数控机床空间动态特性研究[J].工程科学与技术, 2017, 49(4):211-218. DENG C Y, LIU Y, YIN G F, et al. Research on machine tool spatial dynamic characteristics based on response surface method[J]. Advanced Engineering Sciences, 2017, 49(4):211-218.(in Chinese) [5] 葛双好,张蔚波,王佩凤,等.基于模态实验的机床进给系统的分析[J].机床与液压, 2015, 43(1):37-40. GE S H, ZHANG W B, WANG P F, et al. Analysis of machine tool feed system based on modal experiment[J]. Machine Tool&Hydraulics, 2015, 43(1):37-40.(in Chinese) [6] CHOWDHURY A R, PRASAD B, KUMAR V, et al. Design, modeling and open-loop control of a BCF mode bio-mimetic robotic fish[C]//2011 International Siberian Conference on Control and Communications. Krasnoyarsk, Russia:IEEE, 2011:87-92. [7] 姜越友.五轴机床工作空间动刚度测试方法研究[D].成都:电子科技大学, 2019. JIANG Y Y. Research on experimental method for evaluating dynamic stiffness of five-axis machine tool in workspace[D]. Chengdu:University of Electronic Science and Technology of China, 2019.(in Chinese) [8] APPRICH S, WULLE F, LECHLER A, et al. Approach for a general pose-dependent model of the dynamic behavior of large lightweight machine tools for vibration reduction[J]. Procedia CIRP, 2016, 41:812-817. [9] 韩松林.动态切削力模拟加载及误差在线测量装置的研究与开发[D].南京:南京航空航天大学, 2016. HAN S L. Research and development of dynamic cutting force loading and online erro measureing device[D]. Nanjing:Nanjing University of Aeronautics and Astronautics, 2016.(in Chinese) [10] MATSUBARA A, TSUJIMOTO S, KONO D. Evaluation of dynamic stiffness of machine tool spindle by non-contact excitation tests[J]. CIRP Annals, 2015, 64(1):365-368. [11] KONO D, MORIYA Y, MATSUBARA A. Influence of rotary axis on tool-workpiece loop compliance for five-axis machine tools[J]. Precision Engineering, 2017, 49:278-286. [12] GUSEV A A. Time domain finite element estimates of dynamic stiffness of viscoelastic composites with stiff spherical inclusions[J]. International Journal of Solids and Structures, 2016, 88-89:79-87. [13] MOHAMED N, CHONG P L, FARAJI F, et al. Stress analysis of various designs of centrifugal pump impellers using finite element method[J]. Journal of Engineering and Technological Advances, 2022, 7(1):36-46. [14] 刘宗昊.双柔性铣削系统动力学建模与叶片加工稳定性研究[D].北京:北京交通大学, 2020. LIU Z H. Research on dynamic modeling of double flexible milling system and stability of blade processing[D]. Beijing:Beijing Jiaotong University, 2020.(in Chinese)