A Novel Finite Element Method Approach in the Modelling of Edge Trimming of CFRP Laminates

: Nowadays, the development of robust ﬁnite element models is vital to research cost-effectively the optimal cutting parameters of a composite machining process. However, various factors, such as the high computational cost or the complicated nature of the interaction between the workpiece and the cutting tool signiﬁcantly hinder the modelling of these types of processes. For these reasons, the numerical study of common machining operations, especially in composite machining, is still minimal. This paper presents a novel approach comprising a mixed multidirectional composite damage mode with composite edge trimming operation. An ingenious ﬁnite element framework which infer the cutting edge tool wear assessing the incremental change of the machining forces is developed. This information is essential to replace tool inserts before the tool wear could cause severe damage in the machined parts. Two unidirectional carbon ﬁbre specimens with ﬁbre orientations of 45 ◦ and 90 ◦ manufactured by pre-preg layup and cured in an autoclave were tested. Excellent machining force predictions were obtained with errors below 10% from the experimental trials. A consistent 2D FE composite damage model previously performed in composite machining was implemented to mimic the material failure during the machining process. The simulation of the spring back effect was shown to notably increase the accuracy of the numerical predictions in comparison to similar investigations. Global cutting forces simulated were analysed together with the cutting tool tooth forces to extract interesting conclusions regarding the forces received by the spindle axis and the cutting tool tooth, respectively. In general terms, vertical and normal forces steadily increase with tool wear, while tangential to the cutting tool, tooth and horizontal machining forces do not undergo a notable variation. and editing, F.C.-M.; visualization, F.C.-M.; supervision, J.L.C.-S., K.K. and V.A.P.; project administra-tion, F.C.-M.; funding acquisition, J.L.C.-S. and K.K. All authors have read and agreed the version of the manuscript.


Introduction
Composite materials have been often used in high-performance applications because of their excellent structural capabilities such as high strength-to-weight ratio, excellent damage tolerance, or good fatigue and corrosion resistance. These advanced properties make composite components incredibly useful for highly demanding high-tech industries such as the automotive [1], the sportive [2,3], or marine applications [4], among others. The aeronautical sector, where weight reduction is a primary issue, clearly benefits the most from the use of this material. The progressive replacement of typical aluminium components in conventional designs by carbon fibre reinforced polymer (CFRP) components has achieved notable improvements in weight reduction for the aeronautical sector. For example, the implementation in the Boeing 787 of 50% of its own weight in CRFP components has achieved significant weight reductions, which lead to a considerable decrease in maintenance costs, increase fuel efficiency, and improved passenger comfort [5,6].
Although composites are manufactured in a near-net shape, subsequent machining operations are needed to guarantee the strict dimensional tolerances and surface quality demands in aerospace. Commonly, turning, drilling, and milling are the machining operations used to achieve this purpose. Turning is predominately used to adjust the dimensional tolerances of shafts. Thus, the use of cutting parameters to achieve a good surface quality is essential [7]. In general, the selection of high cutting speeds, non-wear cutting tools, and cutting the plies parallel to the fibre direction are the most significant factors to obtain a good surface quality [8,9]. Drilling is used to make bolted or riveted assemblies in composite parts. The main defects observed in this operation are the high surface delamination at the entry "push-out" and the exit "peel-up" plies caused by the high thrust forces experienced during the machining process [10][11][12]. Finally, milling is considered as a finishing operation to obtain high-quality surfaces in complex shapes. However, delamination and burrs with uncut fibres are often observed due to the complex interaction between the cutting tool and the laminate [13][14][15]. These resultant machining defects could lead to the nucleation of cracks in the surface with a drastic reduction in the fatigue life of the components [16]. Therefore, the study of the optimum machining settings to reduce or mitigate these defects has become essential to enhance the structural integrity of parts loaded in service.
Haddad et al. [17] studied the effect of cutting tool type and high speed machining parameters on the edge trimming effect of an M21 CFRP laminate. A four flute end mill was found to generate surfaces with the lowest surface roughness. Cutting temperatures which were measured using a thermal camera were found to be in the region of 80 • C to 115 • C, below the glass transition temperature. The machining length was found to have the highest impact on cutting temperatures compared to other parameters.
Prakash et al. [18] measured acoustic emissions signals from a CFRP edge trimming process to measure the effects of cutting tool and machining parameters on surface roughness, delamination and cutting forces. Artificial neural networks have then been compared against regression model with slightly improved results. The authors found that burr style cutting tools were able to reduce delamination better than the helical shaped tools and that the mean acoustic emissions signal was found to increase with wear of the cutting tools.
Recently, Nguyen-Dinh could be consider one of the most relevant researcher in this topic. Initially, this author used new methods to assess the quality of machining using different experimental methods, scanning electron microscopy (SEM), 3D optical topography and X-ray tomography [19]. They have used areal surface roughness metrics (S a ) and X-ray tomography to assess surface quality and crater volume (C v ), respectively. The applied methods showed improvements in characterising damage using the Altisurf 3D non-contact scanner to take volumetric surface roughness measurements. The cutting temperature and machining induced surface roughness were both found to increase significantly with cutting distance. Low feed and high speed also significantly enhanced tool wear and cutting temperatures.
Lately, the same author [20] investigated the effect of machining parameters and tooling on harmful particles and dust which are produced during machining, measured by laser spectropic dust monitor. It was shown that a high feed speed and lower cutting speed optimised the chip thickness to increase the size of harmful dust particles, which therefore did not become as easily airborne. Additionally, machining at a 3 mm depth of cut compared to 2 mm was shown to improve the size of chip fragments and reduce the quantity of harmful particles. Three types of cutting tool were applied to find their effects on chip and dust production. It was found that the serrated straight flute tool showed some improvement in the emission of dust generated.
Finally, Nguyen-Dinh et al. [21] used new surface quality metrics crater volume and depth of damage and related these to the mechanically tested compressive behaviour of a machined CFRP laminate. These two metrics (which were measured using X-ray tomography and Altisurf techniques) have shown to be good indicators of damage; it was shown that matrix induced machining damage had a more dominating effects on material compressive strength than fibre dominated machining damage. The author's experiments showed that a decreasing trend in the composite compressive strength was correlated with an increasing damage depth.
Several trials are required to fully understand the effect of different cutting parameters, such as cutting speed, feed rate, helix angle, or tool wear, have on machining responses such as sub-surface damage, cutting forces, or surface quality. The complex tool/workpiece interactions, which substantially complicate the design of effective trials, in conjunction with the high cost of composites, turn experimental research into a cumbersome and expensive solution. Therefore, finite element models of composite machining that simulate the underlying physics of the problem and are experimentally validated have been proven to be an excellent cost-effective virtual alternative.
This alternative also provides detailed information on laminate failure. These failures occur mainly in the interfaces (delamination) and in the laminate constituents (fibre and matrix). These failures are essential to predict the debonding between various pieces or plies [22] as well as fibre and matrix characteristic failures in common impact or machining problems [23].
To date, most finite element (FE) investigations in composite machining have investigated orthogonal cutting operations, which shares the most relevant physical aspects of the cutting mechanism characteristic of turning operations [24]. FE orthogonal cutting investigations are focused on the study of the influence that cutting parameters have on machining responses such as cutting forces, sub-surface damage, or chip formation. A brief revision of various relevant papers in this field is collected in the following. Zenia et al. [25] found that increments in the depth of cut produce a significant increase of machining forces and damages underlying the machined surface, while rake angle variations do not significantly affect the machining responses. Wang et al. [26] concluded that, because of the high fibre bending experience in the machining of 90 • , laminates can drastically increase the sub-surface damage in the machining of laminates with similar fibre orientations. Finally, Cepero-Mejias et al. [27] achieved, for the first time, the successful simulation of the characteristic chip fracture parallel to the fibre direction of laminates with positive fibre orientations.
Drilling modelling investigations mostly assess the impact that machining factors have on the surface delamination of the outer plies [26,28]. Phadnis et al. [29] found that the selection of high cutting speeds and low feed rates reduce the delamination levels, torque, and thrust force considerably. Later, Feito et al. [30] developed a simplified drilling FE model, which drastically reduced the computational time required for a complex drilling model in a simulation without a significant decrease in the accuracy of the predictions. The morphology of the drill bit plays a relevant role in reducing the machining-induced damage in drilling operations. For instance, Isbilir et al. [31] concluded that the use of high stage ratios in step drill bits is essential to decrease the push-out delamination. The development of more sophisticated FE models in this field have successfully studied the drilling process in more complicate composite types, such as woven or biocomposites [32,33].
In general, the FE modelling of metal milling operations has been developed in numerous investigations in recent years. However, these models still offer important limitations. For instance, the development of advanced machining friction models is required to improve the quality of feed force predictions [34]. The Johnson-Cook model is the most used to define the plastic behaviour of these metals with linear or exponential progressive damage models [35]. Many advances in this field have been achieved with these investigations. For example, Wojciechowski et al. [36] found that the cut-off transition between the burnishing regime and the chip formation regime is in a feed per tool range of 1-1.4 microns. Finally, Gao et al. [37] showed that the use of advanced modelling techniques such as the Coupled Eulerian-Lagrangian approach significantly improves the predictions quality in terms of chip shape and machining forces.
Although composite milling is one of the most common operations to achieve high-quality surface finishing in the industry, they have not been widely modelled so far. To study these operations, simplified FE orthogonal cutting models are implemented because of the high difficulty of simulating explicitly the characteristic material failure and chip separation of the milling process [38]. Therefore, the development of sophisticated FE models capable of solving the doubts about the physical mechanism of failure concerning composite milling is necessary for future investigations. Furthermore, the implementation of damage algorithms to mimic the characteristic progressive loss of mechanical properties and crack propagation observed in the damaged areas is vital to accomplish this matter [39].
The orientation of the cracks plays a fundamental role in the appearance of defects in the mechanical surface, such as surface roughness or sub-surface damage. A clear increase in these variables is observed in fibre orientations greater than 90 • due to the fact that the great bending effect causes the fibre breakage to occur in a very irregular way below the cutting surface as well as a significant increase in the fibre matrix debonding extension into the sub-surface damage [40,41].
However, the orientation of the fibres is not the only factor that influences the increase in sub-surface damage; the wear of the tool tip also contributes significantly to the increase in the length of the cracks propagating into the machined surface. The impact of blunt tool tips (large tool edge radius) interacting with the fibre tends to push the fibres, generating a significant bending effect rather than shearing them. This mechanism increases the size of the cracks produced by the fibre-matrix debonding below the machined surface, considerably reducing the structural integrity of the part [42]. Therefore, this tool wear must be quickly detected by evaluating variations in the mechanical responses. This information can be used to know when the worn tool inserts should be replaced in order to avoid the massive rejection of parts in a manufacturing production system. The most effective and easy-to-obtain indicator that shows the emergence of severe tool wear is the visualisation of an increment in the cutting forces. Therefore, the investigation of the variations of the cutting forces caused by the tool wear during machining is necessary to enhance the final structural quality of the machined components.
This paper address this matter, developing a consistent FE study of a characteristic milling operation called edge trimming to analyse the effect of tool wear on the machining forces. For this purpose, the machining of two unidirectional (UD) CFRP laminates with 45 • and 90 • fibre orientations were assessed. Tool wear was accounted for by increasing, in several simulations, the cutter tool edge radius, from 10 µm to 50 µm. Laminate behaviour was modelled using a damage algorithm based on a hybrid Puck-Maximum stress failure criteria to determine composite damage onset and subsequently reduce the mechanical properties linearly [43]. The numerical implementation of the spring back phenomenon (partial thickness recovery of the machined surface) was found to be essential to improve the accuracy of the numerical predictions. Numerical results are validated using the experimental data collected in the trials developed by Duboust et al. [44]. Interesting conclusions from the evolution of global forces to the spindle axis (F X = Horizontal force and F Y = Vertical force) and related forces to the cutting tool tooth (F T = Tangential force and F N = Normal force) are collected in this manuscript.

Experimental Set-Up
The milling process differs from orthogonal cutting in that there is a rotating tool with multiple cutting edges and intermittent contact of each cutting edge with the workpiece. Additionally, there is also a constantly changing thickness in the size of the removed chip as the tool is fed into the workpiece. A conventional edge trimming approach was applied in this experiment, which means that the size of the chip will increase as the tool increases engagement with the workpiece. In a composite edge trimming process, the milling kinematics mean that there will be a new fragmented chip generated for every new pass of a cutting tooth. The effect of fibre orientation also plays a role, which will change in relative orientation to the cutting edge with the rotation of the milling tool, as shown in Figure 1. The edge trimming machining experiments were completed on a uni-directional laminate with two different fibre orientations: −45 • and 90 • fibre orientations. The carbon fibre laminates were manufactured by pre-preg hand layup and were cured in an autoclave. A fibre type T700G and an epoxy matrix named Hexply M21 were employed. A five axis Computer Numerical Control (CNC) machine tool, Cincinnati FTV5-2500, Cincinnati machines, Ohio (EEUU), was used to machine the carbon fibre material, see Figure 2. A polycrystalline diamond (PCD) cutting tool was selected. Unidirectional CFRP laminates of 36 plies were vacuum bagged to be cured in an autoclave. The curing cycle selected is specified in the below list.

3.
Reduce vacuum to a safety value of −0.2 bar when the autoclave pressure reaches 1 bar gauge 4.
Vent autoclave pressure when the component reaches 60 • C (140 • F) or below.
This laminate was cut into samples (at the required fibre orientations), which had length, width and thickness-160 mm, 63 mm and 6 mm, respectively. These samples were fixed using a clamping system as shown in Figure 3 and attached to the CNC machine tool. The cutting forces were recorded during the experiment by using a Kistler dynamometer, which was attached to the base of the CNC machine as shown in Figure 3. The effect of cutting tool wear on the machining forces was recorded by measuring the F X (feed force) and F Y (thrust force), depicted in Figure 3. The mean cutting forces were calculated and compared with numerical models, which were then used to validate FE cutting simulations. The PCD cutting tool had a 0 helix angle and 10 mm diameter, as shown in Figure 4. Two cutting tools were used in two different conditions in new and worn condition, with a 3 µm and 10 µm edge radius, respectively. The tool tooth had a clearance angle of 22.4 • , and a rake angle of 29.8 • with a negative 1 mm offset from the centre line. Each test was repeated 2 times with each cutting tool condition. The samples were cut with a full thickness axial depth of cut (ap) of 10 mm and a radial depth of cut of 2 mm. A spindle speed of 6000 rpm and 1200 mm/min feed rate was applied using the three flute PCD cutting tool.

FE Model Characteristics
This work focused on analysing the tool wear influence on machining forces in composite edge trimming. It was achieved by developing 2D FE plane strain simulations in Abaqus/Explicit. In the FE model, just the cutting tool teeth of the three flute, zero helix, and milling tool were simulated to reduce the computational cost. The initial laminate geometry was modelled to mimic the tool path when one cutting edge starts the material removal, as shown in Figure 5. CFRP mechanical and strength properties are shown in Tables 1 and 2, respectively.  Material CFRP [44] 2375 1465 51 119 112 In the above tables, E 11 and E 22 represent the composite Young modulus in directions parallel and perpendicular to the fibre, respectively. The terms G 12 and ν 12 are the composite shear modulus and the Poisson's ratio, respectively. X T and X C denote the longitudinal traction and longitudinal compression laminate strengths, respectively, whereas variables Y T and Y C indicate the transversal traction and transversal compression laminate strengths, respectively. Finally, the term S 12 is the shear composite strength.
The boundary conditions were selected to reproduce the cutting condition observed in real trials reliably. The displacement at the bottom of the laminate was fixed to simulate the fixtures' clamping conditions. The horizontal displacement was restricted to zero on the lateral sides to avoid the simulated laminate bending while the cutter is in contact with the workpiece. A total of 82,376 quadrilateral CPSR4 elements were employed in this research. The mesh employed was designed to guarantee the right balance between the accuracy and the model's computational cost. The element size of 10 µm in the cutting area was selected because machining forces converge around this element size; the element sizes gradually increased until achieving a maximum size of 200 µm at the laminate borders.  The cutting tool was treated as a solid rigid body to reduce the FE model's computational cost. This approach is plausible because of the high rigidity and hardness of the Polycrystalline Diamond (PCD) milling tool and the low deformations levels suffered at the cutting tool tip compared to the machined composite laminate. A constant Coulomb friction coefficient of 0.1 was used to model the friction between the PCD cutter and CFRP laminates accurately [45].

FEM Damage Algorithm Basics
The composite damage model was implemented via VUMAT Fortran user-subroutine in Abaqus/Explicit. Four damage types were employed to reduce the mechanical properties linearly: (1) fibre traction, (2) fibre compression, (3) matrix traction, and (4) matrix compression. These damage types are included in the stiffness matrix to degrade the mechanical properties separately, as shown in Equation (1). The formulation of the stiffness matrix was first proposed by Lapczyk and Hurtado [46] and later applied in composite machining by Cepero-Mejias et al. [43]. 1] and I = ( f t, f c, mt, mc). A hybrid composite damage algorithm was used to predict the onset of the damage modes that were studied. To determine the damage initiation in fibre damage modes, a simple maximum stress failure criteria [47] detailed in the following equations was used.
-Fibre traction (σ 11 > 0) The accurate modelling of matrix damage modes is vital to predict the machining forces reliably [48]; therefore, the use of robust composite failure criteria is essential for excellent predictions. Puck's failure criterion was selected because of its high capability of determining the failure with a large number of off-axis loading scenarios [49]. For the purpose of brevity, Puck's criterion insights are not explained in detail here. This criterion is thoroughly explained in [50]. Three different matrix damage modes were distinguished: Mode A, Mode B, and Mode C. Mode A predicts matrix failure in tensile states, while Mode B and Mode C were utilised to predict matrix failure in compressive loads. These damage modes are formulated below. - In the above equations, the term F I denotes the damage activation functions, where I is whichever of the four investigated damage modes (ft, fc, mt, and mc). All terms with R represent the strength components related to the fracture plane. Term R (+)A ⊥ is generally associated to the transverse matrix traction strength (Y T ), R A ⊥ is considered the composite shear strength (S 12 ). R A ⊥⊥ represents the transverse/transverse shear strength. In general terms, an accurate measurement of this variable presents a cumbersome venture. Therefore, this variable is usually inferred from Equation (7). Term p (+) ⊥ denotes the Puck failure envelope slope of positive matrix stresses (σ 22 > 0) when σ 22 = 0. Finally, the value of the expression p R is extracted using the Equation (8).
Subsequently, a linear energy-based softening is applied to mimic the loss of mechanical properties that each damaged area experiences. This softening has been previously used in previous investigations by these authors [43]. To reduce the mesh size sensitivity, the crack band theory is used [51]. In this theory, the stress and displacement are recalculated, accounting for the characteristic element length, which are renamed as equivalent stress (σ eq ) and equivalent displacement (δ eq ), respectively. This action minimises the strain localisation problems to slow down the decrease of energy dissipated in refinement meshes. This damage evolution is calculated with the use of these equivalent displacements as follows, where the terms δ f I,eq and δ 0 I,eq represent the final and initial equivalent displacements, respectively. As can be seen in Equations (10) and (11), both terms depends on the critical fracture energy of the studied failure mode (G c I ). This fact allows the dissipation of all this fracture energy when the element damage d I reaches a value of 1. In this work, the critical fracture energies used by Zhang et al. [52] were used, see Table 3. Table 3. Critical fracture toughnesses of the studied damage modes extracted from [52].

N/mm
For the 90 • laminate, a maximum of 95% of reduction of mechanical properties in matrix damage modes (d m ) was permitted, whereas in the shear component of the stiffness (d s ), the maximum degradation allowed was 80%. In the case of the 45 • laminate, both features were degraded to a 50% of their initial stiffness. These limitations in the damage reduction were implemented to mitigate the distortional problems that occur to elements when their stiffness is close to zero [53]. We would like to note that the characteristic hardening of matrix stiffness in high strain rates scenarios is not considered in this research. Hence, this variable should be included in further investigations to emulate the behaviour of the material to increase the quality of the numerical predictions.
The chip separation is also not studied in detail in this manuscript. In order to remove the machined material, a strict strain based criteria was chosen. All elements are removed when the element strain achieves a value higher of 1.2 in one of its components (ε 11 , ε 22 and ε 12 ). This value exceeds 100 times the order of magnitude of the strain levels required for an epoxy matrix to fail [54]. Nevertheless, it is maintained because the element erosion produces, in a discrete model, a total contact loss with the cutting tool. This fact significantly reduces the forces that the cutting tool received because an important area in contact with the cutting tool is drastically removed. This factor significantly decreases the accuracy of the prediction of machining forces because the forces coming from the spring back phenomenon are not accounted for [43]. Therefore, this element deletion criterion was used in this research to simulate these spring back forces.
Finally, all the constitutive damage model phases described in this section is represented in Figure 7.

Results and Discussion
The cutting mechanism notably varies with tool wear, as shown in Figure 8. A new tool cut the laminate shearing small chips away. In order to model this problem, the fibre and matrix fracture mechanism should be included in the model to achieve a good quality in the numerical results. However, this configuration was not studied in this manuscript because this matter is out of the scope of this research.
A worn tool edge radius considerably reduces its effectiveness to cut fibres. As a result, the tool edge removes the material due to the high rubbing forces with the laminate, as shown in Figure 8. The increment of tool/workpiece contact area increases abruptly the ploughing forces and sub-surface damage. The proposed FE model takes into account this phenomenon and check how the cutting forces increase because of the tool wear. These results can be employed in a real test to check the tool wear, analysing the machining forces. Therefore, the inserts can be replaced before the worn edges induce excessive laminate damage. The simulations focus on the study of the action of the tool teeth on the material. Therefore, the simulation time selected was the time it takes for one of the tool teeth to be in contact with the laminate. A feed rate ( f r ) of 1200 mm/min, spindle speed (w) of 6000 rpm, and a radial depth of cut of 2 mm were modelled in this research. Five cutting edges, with a radius of 10 µm, 20 µm, 30 µm, 40 µm, and 50 µm, were investigated to assess the tool wear influence on the machining forces. Four distinct machining forces were analysed F X , F Y , F T , and F N , which are represented in Figure 9. F X and F Y denotes the horizontal and vertical forces, respectively. These variables allow us to analyse the machining forces experienced by the cutting tool in the spindle axis. In the case of F T and F N , they are the tangent and normal forces to the cutting tool edge trajectory, respectively. Henceforth, F T is named tangent force and F N corresponds to the normal force for the purpose of clarity. These variables reliably represent the forces that the tool tooth undergoes during the machining process.
In the above figure, θ represents the angle between the line that connects the tip and the rotational centre of the cutter and the horizontal axis. Variables w, f r , and R denote the cutting speed, feed rate, and radius of the cutting tool, respectively. V t is the tangent velocity of the tip of the cutting tooth in contact with the laminate. Finally, β and α correspond to the angles that indicate the orientation of V t and the coordinate system to analyse F t and F n with respect to the global coordinate system, respectively.
To transfer the machining forces obtained in the global coordinates F x and F y to the relative coordinates F t and F n , it is necessary to develop a vector base change. For this purpose, the angle that forms the tangent velocity at the tool tip (V t ) with the horizontal axis of the coordinate system (β) is calculated. The angle between the two vector bases (α) is established then to transfer the forces obtained in the global coordinate system to the new coordinate axes. The relevant equations to complete the change of vector base are set out below. In this work, the spring back phenomenon, which is the partial thickness recovery after the cutter travels through the machined surface, is included. The insertion of this phenomenon in simulations was achieved, avoiding the deletion of the elements in contact with the cutting tool edge. The interaction between the laminate and the cutting tool was observed to displace the tool tooth to the workpiece surface, as shown in. This interaction inserts new forces not accounted for in previous studies, which improves the reliability of the predictions obtained as discussed below.
Another relevant factor that influences the results is the machining relative angle. This factor is the angle between the tangent velocity of the cutting tool tip and the fibre orientation. Along with the advance of the cutting tool tooth, this angle is progressively reduced, notably influencing the machining forces obtained. For 45 • laminates, the machining relative angle drops from 45 • to −15 • , whereas the 90 • laminates go from 90 • to 30 • , as shown in Figure 10.

Numerical FE Model Validation
For the validation of this model, the average force data provided in the research of Duboust et al. [44] were used. All the machining configurations tested in the machining trials developed by Duboust et al. [44] with a cutting tool edge radius of 10 µm were simulated. In total, eight simulations with feed rates of 800 mm/min and 1200 mm/min, spindle speeds of 6000 rpm and 8000 rpm, and fibre orientations of 45 • and 90 • were performed. The inclusion of the spring back phenomenon, together with the limitation in the damage of shear and matrix applied, can simulate trends in machining forces very similar to those obtained experimentally. Finally, a definite improvement was found when compared to the predictions made by Duboust et al. [44], as shown in Figures 11-14. The small error percentages simulated show the reliability of this FE model in predicting accurate machining forces in any cutting configuration. These errors do not exceed 10% in most of the average forces simulated, as shown in Table 4. Only the configuration with f r = 1200 mm/min, w = 8000 rpm, and 90 • obtained an error higher than 10% in F Y forces. However, a glance to Figure 14 reveals that the standard deviation calculated for these trials is considerably bigger than the standard deviation obtained in the rest of the trials. Furthermore, the model predicts a similar trend observed in the experimental findings for F X and F Y forces. From these factors it is concluded that a small error in the measurements of these forces could be the cause of this discrepancy. Table 4. Numerical errors obtained with respect to the experimental average forces measured by Duboust et al. [44].

Parameters
Average

Fibre Orientation of 45 •
Initially, the vertical forces F Y experience a steady growth mainly motivated by the spring back phenomenon. This takes place because the range of machining relative angles given in the edge trimming of 45 • laminates, from 45 • to −15 • , generate similar thrust forces in machining [55]. A parabolic decrease of the forces was observed because the spring back forces have a larger contribution to the X axis, which decreases notably the quantity measured in the vertical component. Horizontal forces F X experience a parabolic increase due to the contribution of the spring back phenomenon to the horizontal axis. An analysis of this result reveals that cutting forces are not important in comparison to spring back forces, which are the leading cause of the F X and F Y force distribution. The global machining forces simulated for cutting tool edge radius of 10 µm, 20 µm, 40 µm, and 50 µm are presented in Figure 15.
Note that the cutting forces recorded in the simulations do not follow a smooth distribution with the advancement of the tool, but rather an irregular distribution full of saw teeth. These fluctuations are explained due to the changes in the normal contact stiffness between the cutting tool and the laminate that occurs when the tool comes into contact with the different elements of the meshed laminate. When the cutting tool tip contacts with a new element, there is a slight rebound in the contact forces because of the fast increment of the contact rigidity that progressively diminishes as the tool advances its truck over the aforementioned element, as shown in Figure 16. The magnitude of these fluctuations increase notably with the implementation of higher contact stiffness and could create element distortion problems with the use of high values. In this research, a linear contact rigidity of 10 GPa was employed to avoid this problem.
Regarding the forces F T and F N , two steady rises with distinct slopes were observed. Therefore, it is concluded that the variation of the effect of the changes in machining relative angles is not crucial for this configuration. As discussed above, it was observed that spring back forces, which are responsible for the increment of F N forces, are higher than the cutting forces that govern the F T forces evolution. The representation of all graphs measuring these forces is showcased in Figure 17.  From the analysis of the results, it is inferred that tool wear does not increase significantly (less than 20%) the machining forces in the edge trimming of 45 • laminates. The most significant increments are observed in forces F N and F Y with maximum rises of 19.90% and 13.76%, respectively. The force F X keeps stable with small increments below 5%. Finally, forces F T are calculated to be low in comparison to the other forces analysed, as they experience a tight decrease from 23.63 N to 17.78 N. The average machining forces simulated in this work are presented in Figure 18, whereas the force increments obtained with respect to the unworn cutting tool, with a tool edge radius of 10 µm, are shown in Table 5.

Fibre Orientation of 90 •
An observation of the vertical forces obtained reveals that these initially undergo a steep rise due to the effect of the spring back that pushes the cutting tool in the normal direction to the machined surface. This steep rise slope occurs because of the cutting tool being initially pushed with the characteristic high rigidity of laminates that are loaded parallel to the fibre. After this sharp rise, there is a decrease in the slope of this curve until it reaches a significant drop attributable to two factors. The first factor that causes this drop is that the contribution of forces produced by the spring back phenomenon in the vertical component decreases appreciably with the advance of the cutting tool tooth. The second factor is that the cutting forces generated during the machining process are added in the negative component of the vertical axis.
Horizontal force steady rises because of the horizontal contributions given during the simulation of the spring back and cutting forces. The horizontal cutting forces are primarily increased because of the increase in chip thickness as the tool rotates. Subsequently, the forces coming from the spring back increase their contribution in the horizontal axis while the cutting force contribution becomes less relevant with the advance of the cutter tooth. Figure 19 visualises the global machining forces collected for the cutting tool edge radius of 10 µm, 20 µm, 40 µm, and 50 µm, which is the subject of the above statements. Regarding the forces measured in the tooth of the cutting tool F N and F T , they reach a stable level of force over which they fluctuate. Different reasons explain this trend for both investigated forces. In the case of F N , this force is obtained as a consequence of the spring back phenomenon. Initially, a rapid force increase is observed because the penetration stiffness is high when the tool pushes the laminate in a direction parallel to the fibre. As the tool advances, the relative machining angle decreases in a range from 90 • to 30 • , softening the workpiece/cutting tool contact and substantially reducing the contact stiffness. The level of force previously obtained is maintained due to the deeper penetration of the cutting tool tooth (with an increase in chip thickness), however this is compensated by the contact stiffness reduction when the cutting tool tooth advances. An identical mechanism occurs with the tangential forces F T , where the cutting forces decrease with relative cutting angles lower than 90 • . However, in this analysis, this factor is compensated by the larger amount of material removed when the simulation goes ahead. As seen in Figure 20, the reached force levels of F N are significantly higher than those of F T because each pass of the tool tooth has a small depth-of-cut. The average forces simulated in this investigation highlight several assertions. The machining forces tend to experience a linear incremental increase produced by the effect of the cutting tool wear. A gradual and significant increase in the forces F N and F Y caused by tool wear was observed to increase in maximum forces of 38.79% and 47.45%, respectively. This rise is less pronounced for F X forces, achieving a maximum increment of 16.12%, while the F T forces do not show a significant variation due to the wear of the cutting tool. This effect is due to the increased contact area between the workpiece and the worn cutting tool. All the data obtained in the simulations performed are illustrated in Figure 21, while the force increases obtained in comparison with the forces obtained with the unworn tool with a tool edge radius of 10 µm, are shown in Table 6. Table 6. Average forces and increment percentage because of the tool wear simulated.  A deeper analysis of the results simulated in 45 • and 90 • laminates yield the following conclusions. The increment of the machining forces occurs because the tool edge wear decrease the fibre cutting effectiveness. A worn tool, instead of shearing the fibres away, tends to bend the fibres, increasing the machining forces and sub-surface damage. Generally, the increment of machining forces occurs because the tool/workpiece contact area is increased. This factor increases the spring back and frictional forces. The spring back forces increase significantly because a bigger area of the tool penetrates inside the laminates, adding a substantial force in the normal component (F N ). The tangential force (F T ) changes are smaller because, despite the friction forces increasing due to contact area increments, the tangential component of the spring back forces is contrary to the frictional force. Therefore, both effects are counteracted, avoiding a clear trend in the tangential forces. Horizontal and vertical force (F X and F Y ) increments are the result of the projections of the significant spring back force increment. In general, the increments of F Y caused by tool wear are more severe than the force increments observed in F X .

Conclusions
In this research, we developed a novel study for the modelling of edge trimming of CFRP laminates using a worn cutting tool. The tool wear effect on machining forces was successfully modelled in 45 • and 90 • laminates. A composite damage algorithm based on the hybrid Puck-Maximum stress criteria was used to calculate damage initiation with a linear energy-based degradation of the mechanical properties and to mimic composite behaviour. Validation of average machining forces predictions was achieved by obtaining relatively small errors (within 10%) compared to experimental data. These predictions notably enhance the numerical results achieved in previous investigations in this matter. These improvements are obtained because of the numerical implementation of the spring back phenomenon, which adds to the simulation machining forces. Several relevant conclusions extracted from this study are summarised in the below bullet points.

•
In all developed simulations, the normal forces (F N ) simulated are much higher than the tangential forces (F T ); meaning that spring back forces are predominant to the cutting forces in these analyses. • The rise observed in F N , F X , and F Y is result of the increment of spring back forces, while F T keeps stable because the increment of the spring back and frictional forces is compensated for. • Increments in F Y caused by tool wear were observed to be more severe than increments in F X . Therefore, the evolution of vertical forces should be investigated to determine the appropriate timing to replace worn inserts. • A steady increase in normal (F N ) and vertical (F Y ) average forces because the effect of the tool wear was observed in both studied laminates. These increments are more severe in 90 • laminates. • Normal and tangential machining forces depict a steady rise in 45 • laminates. This happens because the range of machining relative angles from 45 • to −15 • does not modify the spring back force significantly, nor cutting forces. Therefore, this rise is governed by the progressive increment of the depth-of-cut when the cutting tool tooth advances. • Normal machining forces of 90 • laminates have an initial steep rise to later fluctuate between a stable value. In the normal force, the cutting tool tooth is pushing down the laminate parallel to the fibre, finding a high contact rigidity that decreases significantly when the machining relative angle is far from this 90 • .
Note that interesting technical applications could be achieved using the FE model developed in this investigation. For instance, it could detect the point where the tool should be replaced, thus saving manufacturing time and optimising the tool's use. Apart from this point, another analysis could be developed using this framework. The FE model developed in this research offers an excellent capability to analyse how other relevant cutting parameters such as feed rate, cutting speed, fibre orientation, or rake and relief angles affect the machining forces. Therefore, this FE model or similar FE models can be used in oncoming composite milling investigations to quantify the effects of cutting edge rounding and CFRP machining failure mechanisms.