Experimental study on machining deformation of large scale slender beam with weak stiffness of Ti6Al4V

Large-scale slender beam structures with weak stiffness are widely used in the aviation field. There will be a great deformation problem in machining because the overall stiffness of slender beam parts is lower. Firstly, the cutting mechanism and stability theory of the Ti6Al4V material are analyzed, and then the auxiliary support is carried out according to the machining characteristics of the slender beam structure. The feasibility of the deformation suppression measures for the slender beam is verified by experiments. The experimental analysis shows that on the basis of fulcrum auxiliary support, the filling of paraffin melt material is capable of increasing the damping of the whole system, improving the overall stiffness of the machining system, and inhibiting the chatter effect of machining. This method is effective to greatly improve the accuracy and efficiency during machining of slender beam parts. On the premise of the method of processing support with the combination of fulcrum and paraffin, if the tool wear is effectively controlled, the high precision machining of large-scale slender beams can be realized effectively, and the machining deformation of slender beams can be reduced. Although high speed milling has excellent machining effect on the machining accuracy of titanium alloy materials, severe tool wear is observed during high-speed milling of titanium alloy materials. Therefore, high-speed milling of titanium alloy slender beam is suitable to be carried out in the finishing process, which can effectively control tool wear and improve the machining accuracy of parts. Finally, the process verification of typical weak stiffness slender beam skeleton parts is carried out. Through the theoretical and technical support of the experimental scheme, the machining of large-scale slender beam structure parts with weak stiffness is realized.


Introduction
Titanium alloy is one of the main materials of the whole structure of aviation. With the development of aviation technology in recent years, the application of titanium alloy materials in high maneuvering aircraft has been widely popularized, especially Ti6Al4V. It is a heat-resistant alloy, its elastic modulus is 70.3 GPa, the compression modulus is 73.8 GPa, the density is 4.44 g/cm3, the Poisson's ratio is 0.33, and the melting point range is 1560~1680℃. Its chemical composition is shown in Table 1. Ti6Al4V material is mainly used to make mechanical structure of aviation products, such as frame, impeller, wing plate and beam structure. However, there are still some technical problems in the application of titanium alloy materials for processing of large-scale slender beam structure at home and abroad and, for example, serious deformation is the most common problem encountered during machining [1][2]. Although the design of slender beam structure has many advantages, when the workpiece is removed from the fixture, deformation such as bending, twisting, or combined bending and torsion often occurs during machining due to its large aspect ratio, high material removal rate and poor rigidity, which makes it difficult for the parts to meet the design requirements.
The Nervi S of the United States established a mathematical prediction model of machining deformation caused by initial residual stress of blank, and studied the relationship [3] between structural deformation and the distribution state of initial blank stress [3]. Keith studied the influence of residual stress on machining deformation. Results indicated that the machining stress was closely related to the radius of arc tip and the radius of blunt circle of cutting edge. For aeronautical thin-walled structure, the machining residual stress had a certain effect on the machining deformation of workpiece [4]. By using the neural network theory, Ratchev and others established the flexible prediction model of the milling deformation of thin-walled parts, and then established the prediction model of the machining error of the thin-walled workpiece by combining it with the finite element analysis method [5][6][7]. Tlusty proposed the effective use of unprocessed parts as support, and explored the method of milling path optimization using the overall rigidity of aeronautical structural parts [8].
Hiroyasu studied the deformation of thin-walled parts caused by cutting force, and designed the parallel double spindle machining method [9].  -walled parts for the whole impeller from the aspects of blank and tool selection, fixture design,   machining process flow of integral impeller, NC program compilation and simulation, specimen processing and measurement, ensuring feasible mass milling of thin-walled parts with high efficiency, high quality and low cost [10]. Dr. Wang Ming from Harbin University of Technology selected the cutting parameters of stable point by MATLAB simulation. The results showed that the stability analysis method made the reasonable selection of cutting parameters possible which could avoid the occurrence of chatter, thus improving the machining quality and efficiency of thin-walled parts [11].
Dong Guojun, Harbin University of Technology, has carried out molecular dynamics simulation and experimental research on tool wear mechanism, and studied tool wear mechanism according to the change of chemical composition of tool tip, which provides technical support for cutting tool wear [12].
Although many scholars at home and abroad have done a lot of academic research in thin walls with weak stiffness, impellers and frame parts, there is little research on the machining deformation of slender beam parts with large scale aspect ratio, and a breakthrough is needed in this direction.
Processing a slender beam structure becomes difficult because it will be in the cantilever state and deformation and chatter will occur during machining. Therefore, the paper takes the slender beam structure with weak stiffness as a research object and studies a method to suppress the deformation occurred during processing. In addition, experiments are performed to further study the deformation behavior of a cantilever structure with weak stiffness during processing by considering the properties of Ti6Al4V. The clamping method of the slender beam structure is studied to improve the overall machining stiffness and damping to make it stable. Based on the optimal clamping mode, the influence of different tool wear and milling parameters on the machining deformation of weak stiffness slender beam is explored, which provides technical support for the machining technology and clamping method of typical slender beam parts with weak stiffness.

Cutting Mechanism Analysis
The chip formation process of metal machining is actually the deformation process of material cutting layer. The chip separation process can be roughly divided into the first deformation zone (shear slip), the second deformation zone (fibrosis), and the third deformation zone (fibrosis and work hardening).
As shown in fig.1a, the first deformation zone interacts with the second deformation zone. When the friction force of the knife surface is large, the chip discharge is not smooth, the extrusion deformation intensifies, and the shear slip deformation in the first deformation zone increases. In the contact area between the machined surface and the rear tool surface, due to the extrusion friction between the tool edge circle and the rear tool surface, the machined surface is fibrosis and work hardening, thus forming the third deformation zone [1].

Fig. 1 Principle of Metal Cutting
Burr is a part of the edge area of the workpiece, which is a part of the chip removal in theory. It is a part of the machining process. The definition of burr parameters is an important standard for burr study and evaluating the energy consumption of deburring. The basic parameters of burr are usually studied by arbitrary cross section of burr. The main parameters of burr are height, root thickness and root radius.
As shown in fig.1a, burr height H refers to the maximum distance between the ideal end face of the workpiece terminal surface measured on the cross section and the profile of the burr cross section. Burr root thickness W refers to the distance from the plastic deformation convex starting point to the ideal surface of the workpiece measured on the end surface of the workpiece. The radius of the burr root circle is the radius of curvature of the root geometric curve of the burr cross section measured on the cross section. Obviously, the size of burr height H directly affects the dimensional accuracy and shape accuracy of the workpiece, while the thickness of burr root W has a relation to the efforts of deburring.
By using characteristic parameters H and W, the specific size and shape of burr cross section can be expressed [13].
The cutting temperature produced by the milling cutter in the machining process can simplify the three-dimensional cutting into a simple two-dimensional plane oblique angle cutting, that is, the milling cutter can be separated into countless units along the spiral cutting line. The cutting process of a single 5 unit is simplified diagonal cutting. Figure 1b is a simplified diagram of milling temperature model. The milling temperature is mainly from three regions: (1) The cutting temperature region of the front cutter surface caused by friction between the chip and the front cutter surface; (2) The cutting temperature region of the rear cutter surface produced by the friction between the rear cutter surface and the workpiece and (3) Temperature region of plastic deformation of shear surface produced by shear slip of metal plastic deformation [14].
As shown in Fig. 1c, the action of the arc radius of the tool resulted in a triangular region where tensile stress and extrusion stress combined on the workpiece when cutting. The grains in this region were elongated along the Y-axis due to stretching, and plastic shortened along the X-axis because of compression, finally forming a residual tensile stress under the action of these two behaviors.
Consequently, the third deformation zone was formed in the workpiece, and the friction force Fx and normal force Fy acted on the surface layer of the workpiece. The normal force caused plastic deformation on the surface layer in the direction of Y-axis while the friction force Fx induced the surface layer to produce plastic deformation along the X-axis direction, which led to residual compressive stress on the machined surface [15].

Dynamic Milling Force Model
The mathematical model of dynamic milling stability is the basis for the study of machining process stability and dynamic machining error. The difficulty of dynamic milling stability modeling lies in the modeling of cutting thickness. As shown in Fig. 2, the milling system can be simplified into two "spring-damping" systems with vertical degrees of freedom.

Fig. 2 Milling Vibration Systems
The milling dynamics equation can be expressed as a differential equation where mx、my -Quality of milling system in X, Y direction(kg); 6 cx，cy -Damping of milling system in X, Y directions(Ns/m); kx，ky -Stiffness of milling system in X, Y direction(N/m) ; Fxj，Fyj -Cutting force component (N) in X and Y direction acting on cutter teeth j. Figure 3 shows the force distribution of milling cutter blade, in which the Ft is tangential cutting force and the Fr is radial cutting force. In the milling process, the machining system is excited by the milling force in the X and Y directions, resulting in the corresponding dynamic displacement x and y.

Fig. 3 Force distribution of milling cutter
Coordinate transformation of dynamic displacement in cutting thickness direction of cutter teeth vj=-xsinφj-ycosφj. As the two-dimensional static cutting thickness has nothing to do with the mechanism of regenerative chatter, it can be ignored in the flutter stability analysis. So the two-dimensional dynamic cutting thickness can be expressed in the following formula [10].
where (x, y) represents the dynamic displacement of the current tool tooth period, and (x0, y0) represents the dynamic displacement of the former tool tooth period. The rotation angular velocity of the spindle is Ω(rad/s). Then the radial contact angle φj is transformed into Ωt with time. φj=(j-1)φp+φ is the instantaneous radial contact angle of cutter tooth j. Since the number of cutting edges j is 4, the cutting angles of the four cutting edges are：φ0 = φ，φ1=φ+π/2，φ2=φ+π，φ3=φ+3π/2。 The tangential cutting force Ftj and the radial cutting force Frj of acting on the cutter tooth j are [10] ( ) where ap is axial tangential depth(mm); Kt，Kr-tangential and radial cutting force coefficients, respectively(N/m).
The total cutting force acting on the tool is Finally, we substitute formula (2) and (3) into formula (4) and express them in the form of matrix Considering the variation of these parameters with time and angular velocity, equation (6) is expressed in matrix form in time domain as follows where A(t) is a periodic function of ω=NΩ and T=2π/ω. The expansion form of Fourier series is Ignore the influence of higher harmonics and retain the DC component By substituting equation (8) into equation (6), the expression of dynamic milling force coefficient can be simplified as

Flutter stability limit
As shown in Fig. 2, a tool-workpiece system can be simplified to X and Y two-degree-of-freedom vibration systems in two vertical directions [13]. Let the frequency response function of the tool workpiece system be (10) Where Gxx(iω) and Gyy(iω) represent the direct frequency response functions of the tool workpiece system in X and Y directions respectively. Gxy(iω) and Gyx(iω) represent the cross frequency response functions of the tool workpiece system in X and Y directions respectively. Let the relative vibration vectors of the tool and the workpiece at the current time (t) and the previous cutting time (t -T) be The frequency domain of the relative vibration function of the tool and workpiece at the current time and the cutting time of the previous cutter tooth at flutter frequency ωc is expressed as [16] ) Let the determinant be zero, find the special solution of the equation, and the characteristic equation of the tool workpiece closed-loop dynamic system is The characteristic equations obtained by further simplification are as follows The eigenvalues can be obtained by neglecting the influence of the cross functions Gxy and Gyx： The upper imaginary part is zero because the axial cutting depth ap is a real number, and the final expression of ultimate axial depth can be simplified It can be seen that if the flutter frequency ωc is given, the limit axial cutting depth can be obtained according to the above formula, and Where Ψ=arctanκ is the phase shift of eigenvalue. ε=π-2Ψ is the phase shift between internal and external modulation. Therefore, if the vibration ripple left by the cutting arc is an integer k, ωcT =ε+2kπ. The spindle speed can be obtained by calculating the cutting cycle T It can be concluded that in order to improve the cutting stability, the high precision parts can be processed by increasing the damping and the overall stiffness of the machining system as well as reducing the cutting force.

Influencing factors and restraining methods of processing deformation of slender beam with weak stiffness
For the deformation problems of weak stiffness parts, the influencing factors of deformation in the processing of suspension beam parts with weak rigidity need to be investigated first, as shown in Fig. 4.
It is found that the problems of low stiffness and deformation can be effectively solved by increasing the overall machining stiffness of weak stiffness parts and the damping of the system. At the same time, tool wear and machined surface stress can be improved.

Fig. 4 Factors influencing deformation of slender beam parts with weak stiffness
The clamping mode of weak stiffness parts is the most important factor that directly influence the machining accuracy of such parts. When clamping slender beam parts with weak stiffness, pressure plate clamping is not recommended because the special suspension state of the parts will cause the parts to be subjected to cutting force during the machining process, and hence the suspension deformation will occur because of the weak stiffness. Upon unloading, even if the condition of the parts after machining is satisfied, the machined parts will suffer deformation due to the release of the external clamping force. Therefore, it is necessary to ensure the overall processing stiffness of the parts during processing, so that it not only can enable the stable clamping of the parts, but also improve the stiffness of the overall workpiece processing system. To achieve this, the auxiliary support is necessary for the cantilever parts aiming to inhibit the chattering or deformation phenomenon occurs during machining. On the other hand, paraffin filling is carried out on the basis of auxiliary support to indirectly increase the damping of the workpiece processing system, which is capable of further improving the machining accuracy. The filling of paraffin melt can make it possible to suppress chatter in the process of machining, thereby greatly improving process efficiency. Figure 5 shows the effect of deformation with or without auxiliary support.

Experimental scheme of deformation of slender beam with weak stiffness
In this experiment, the clamping methods of the slender beam are designed with no auxiliary support, fulcrum auxiliary support, paraffin auxiliary support, and fulcrum + paraffin auxiliary support. The number of fulcrum is divided and increased according to the length of suspension beam. The processing deformation scheme of slender beam with weak stiffness is shown in Table 2. The experimental environment is shown in Figure 6.

Influence of processing support mode on deformation of slender beam with weak stiffness
The cutting parameters are n=4000r/min, F=800mm/min, and ap=0.02mm. At the end of the experiment, the deformation of the workpiece in Z and Y directions was measured online by 3D edge finder. A micrometer was used to measure machining dimension of parts. When measuring the deformation of slender beam, mark the maximum deformation and measure the surface residual stress with residual stress measuring instrument. Table 3 is the experimental data of processing deformation of the slender beam with weak stiffness, and Figure 7 is the influence curve of processing method on the deformation of the slender beam.

Fig. 7 Influence of processing method on deformation of slender beam
It can be seen in Figure 7 and table 3 that chatter occurred and sharp squeal was produced if the slender beam was machined without auxiliary support, indicating that the chatter was severe during cutting. The reason for the large deviation is that the deformation of slender beams without auxiliary support occurred. From the data table 4, it is observed that as compared to the chatter phenomenon of the slender beam processed without no auxiliary support, this behavior was gradually improved when processing the beam with two-fulcrum auxiliary support, the chatter pattern on the surface was gradually reduced, and the deformation and dimension deviation occurred during machining of the slender beam were gradually improved. Moreover, the residual stress at the maximum deformation becomes smaller after machining. It shows that the unstable chatter phenomenon can be improved and the residual stress is reduced and the machining deformation of slender beam parts can be reduced by upgrading the clamping mode and using the fulcrum support method. When the number of fulcrum was increased to 4, the degree of deformation and machining accuracy of slender beams were controlled within 0.03mm. Although the fulcrum auxiliary support is an effective option to improve the machining accuracy of the slender beam, the combination of fulcrum and paraffin clamping processing can further optimize the machining accuracy. An increase in the number of paraffin and fulcrum provides an opportunity to reduce the deformation and improve the processing accuracy of the slender beam and the residual stress at the maximum deformation of the slender beam after processing is also significantly reduced. When the processing mode was molten material combined with four fulcrums, the dimensional accuracy of beam was 0.016 mm, and the deformation was less than 0.01 mm. From the comparative experiments of 10 clamping methods, it is found that in order to improve the overall processing stiffness and system damping, using the clamping method of the combination of molten 14 material and fulcrum is a breakthrough to realize the machining of large-scale slender beams.

Influence of tool wear on machining deformation of slender beam with weak stiffness
The optimal processing support mode was selected from the previous experimental results: four fulcrums + paraffin auxiliary support. The machining experiments of the slender beam with weak stiffness were carried out by using different worn tools. The cutting parameters are n=4000r/min, F=800mm/min, and ap=0.02mm. At the end of the experiment, the deformation of the workpiece in Z and Y directions was measured online by 3D edge finder. A micrometer was employed to measure machining dimension of parts. Figure 8 shows the curve of machining deformation of slender beam.

Fig. 8 Influence curve of processing deformation of slender beam
It is seen from the data in Figure 8 that tool wear had a great influence on the machining accuracy of the slender beam. The increasing tool wear caused the machining deformation and dimensional deviation of the slender beam to increase. The deformation and machining accuracy of the slender beam were within 0.05mm. It can also be observed in the Table and the Figure that if the tool wear was effectively controlled, the machining accuracy and deformation of the slender beam could be effectively improved. The reason to achieve such a high-precision processing is that the experimental processing support method improves the overall processing rigidity of the slender beam, and the paraffin also indirectly improves the damping of the processing system and suppresses the chatter phenomenon. Therefore, the experimental results show that under the premise of four fulcrum and molten material processing support mode, the high-precision processing of large-scale slender beam can be effectively realized if the tool wear is effectively controlled, and hence the processing deformation of the slender beam can be reduced.

Influence of milling parameters on machining deformation of slender beam with weak stiffness
The best processing support mode was selected four fulcrums + paraffin auxiliary support. At the 15 end of the experiment, the deformation of the workpiece in Z and Y directions was measured online by 3D edge finder. A micrometer was used to measure machining dimension of parts. Figure 9 shows the influence of milling parameters on machining deformation of the slender beam. Fig. 9 Influence curve of milling parameters on machining deformation of slender beam It can be seen in Figure 9 that with the increase of the spindle speed, the machining deformation of the slender beam was gradually reduced, and the machining dimensional accuracy was also gradually improved. However, an increase in feed speed and cutting depth caused the machining dimensional accuracy to become worse. The graph illustrates that the cutting depth has a great impact on the machining deformation of the slender beam while the spindle speed and feed speed have a certain regular influence on the deformation, the change was relatively flat. Although high-speed milling has excellent machining effect on the surface quality of titanium alloy materials, severe tool wear can be 16 observed when high speed milling of titanium alloy. Therefore, high speed milling is suggested to be used for the finishing process of the slender beam, which not only can effectively control the tool wear, but also improve the machining accuracy of parts. Figure 10 shows the 3D model of the typical slender beam with weak stiffness to be machined. The model parts are verified by the previous experiments. After the processing, the dimension of the slender beam was measured to evaluate the machining accuracy. Being similar to the processing of ordinary structural parts, it is necessary to analyze the drawings of the parts and determine the process flow. Figure 11 is the process flow chart. According to the clamping experiment of titanium alloy structure parts with weak stiffness, process and fixture design was carried out for typical weak stiffness parts to meet the dimensional accuracy requirements, and repeated clamping operations can be avoided. As shown in Fig. 12, the fixture is provided with a pouring port and an air outlet for pouring liquid materials. The contact surface between the blank and the fixture is bonded with AB glue to avoid external force as much as possible, and reduce the clamping deformation caused by the release of the external force. The process is shown in Table 4. CAM cutting strategy is as follows: the outer profile milling path is from the outer edge to the inner ring milling, and the depth is first. The milling path of the inner cavity is from the inner to the outer edge and the depth is first. The common cutting depth of each tool is constant; the milling tool needs to cut smoothly during machining to make the cutting state stable, and all tool paths of the tool must follow the arc feed path to prevent the deformation of parts due to overcutting. Avoid tool path repetition. Avoid unnecessary damage to workpiece surface caused by cutting tool during machining. The tool path extends a step on the edge to remove the edge burr. The CAM tool feeding strategy is shown in Figure 13 below.

Fig. 13 CAM cutting strategy
As shown in Figure 14, the surface quality of the processed weak rigid parts is in good condition, the maximum deformation of the longest slender beam in the middle is 0.019mm, the size deviation of the 420mm×8.6mm×8.6mm slender beam in the middle is -0.017mm, the machined surface is bright, and the roughness is less than 0.5μm.

Conclusions
Aiming at the deformation problem of large-scale slender beam with weak stiffness, this paper has 20 carried out theoretical analysis on milling stability, and has performed experiments by adopting some methods to restrain the deformation. In the end, the processing of typical large-scale parts with weak stiffness was verified. The conclusions are as follows: (1) From the theoretical analysis, it can be concluded that improving cutting stability can be achieved by increasing the damping, improving the overall stiffness of the processing system, and reducing the cutting force so as to effectively achieve high-precision processing of parts.
(2) When processing the slender beam parts, the fulcrum and paraffin combined with auxiliary support is an optimal way to increase the overall system damping, improve the overall stiffness of the processing system, and suppress the processing chatter effect. This method is capable of greatly improving the accuracy and efficiency during processing of slender beam parts.
(3) Under the premise of fulcrum and paraffin combined with auxiliary support method, if the tool wear is effectively controlled, the high-precision machining of large-scale slender beam can be effectively realized, and the machining deformation of the slender beam can be reduced.
(4) Although high-speed milling has excellent machining effect on the surface quality of titanium alloy materials, severe tool wear is observed in the process of high-speed milling of the titanium alloy.
Therefore, high speed milling is recommended for the finishing process of the slender beam, which not only can effectively control the tool wear, but also improve the machining accuracy of parts. (5) Through the theoretical and technical support of the experimental scheme, the machining of typical large-scale slender beams with weak stiffness can be realized.

Author contribution
The authors discussed each reference paper together and contributed useful ideas for this manuscript. The authors sincerely thank Professor Hua-dong Yu of Changchun University of Science and Technology for his critical discussion and reading during manuscript preparation.
Funding This work was supported by "Supported by national basic scientific research program" (JCKY2017208B006) and "Stability support project"(WDZC2019JJ016).

Declarations
Competing interests The authors declare no competing interests.