Influence of the Temperature-Dependent Characteristics of CFRP Mechanical Properties on the Critical Axial Force of Drilling Delamination

Previous studies have often assumed that the mechanical properties of Carbon Fibre Reinforced Plastics (CFRP) remain unchanged during drilling. In fact, due to the increase in drilling temperature, the mechanical properties of the composites change greatly, and this then affects the critical force. In addition, previous studies have often assumed that the failure mode of CFRP drilling was a type I crack failure. In fact, due to the complexity of the CFRP drilling process, the failure modes are often coupled with different failure modes, so type I cracks alone cannot reflect the actual cracking situation. Therefore, a three-dimensional drilling Finite Element Modeling (FEM) was established to analyze the failure modes of CFRP drilling delamination, and the I/III mode was determined; then, a new drilling critical force mechanics model, which considers the temperature dependence of CFRP mechanical properties and the failure modes of CFRP drilling delamination, was established based on the classical drilling critical force mechanics model; the results of the mechanics model were validated by drilling critical force experiments under different temperatures. The effects of the temperature dependence of CFRP mechanical properties on the drilling critical force were investigated and analyzed.


Introduction
In recent years, carbon fiber-reinforced polymer composites (CFRP) have been widely used in aviation, aerospace, subway tracks, and other fields. However, due to their anisotropy, multiphase, and heterogeneous nature, they are prone to drilling delamination damage in the drilling process, which seriously affects the reliability of CFRP parts [1]. To reduce the delamination damage in CFRP drilling, scholars have conducted a lot of research on the delamination damage in the drilling process around the drilling force. It is generally believed that there is a critical axial force (P C ) [2]. When the drilling axial force is lower than P C , delamination damage will not occur during drilling; on the contrary, delamination damage will occur when the force is higher than Pc [3]. Based on these studies, a series of analytical models for predicting the P C has been established [4]. For example, Hocheng et al. [5] proposed the concept of P C in 1990 and established the mechanical analytical model of P C based on linear elastic fracture mechanics and classical beam-plate bending theory. They believe that delamination defects do not occur when the axial force is less than a certain critical value in the process of composite drilling. However, their model assumes that the composite material is isotropic, which is inconsistent with the actual situation of the material, resulting in large errors when predicting the P C .
Since then, many scholars have improved the model of P C based on the Hocheng model. For example, Lachuad et al. [6] established an analytical model of P C with anisotropic mechanical properties of composites based on the classical laminated plate theory and Hocheng model, and they proved that the analytical model established by distributed load has higher accuracy than the concentrated load model for the first time. Jain et al. [7] first introduced the idea that the shape of the drilling delamination area is elliptic and characterized the delamination area by a coefficient, which improved the prediction accuracy of the P C analytical model. However, neither the model of Lachuad [6] nor Jain [7] considers the bending torsion coupling effect of laminates. Therefore, their models can only predict the P C of unidirectional laminated plates. Zhang et al. [8] considered the bending torsion coupling effect of laminates based on Jain [7] and established an analytical model of P C that can analyze any laminate stack sequence. To further improve the accuracy of P C , Ojo et al. [4] subdivided the chisel edge and the main cutting edge of a drill. They assume that the force caused by the chisel edge is a concentrated force and that the force caused by the main cutting edge is a uniformly distributed force. In their model, they put forward the hypothesis that the concentrated force and the uniformly distributed force both exist, and they confirmed that the ratio of concentrated force and distributed force has an impact on the P C . In addition, the model also puts forward the assumption that the drilling delamination crack propagation form is I/II mixed type and analyzes it. Saoudi et al. [9] first studied the analytical model of the P C under mechanical-thermal coupling. The results showed that the drilling P C would change under the influence of drilling temperature, but the mechanical properties of the CFRP materials were not considered in the model. Table 1 is a summary of the P C model in the above references, where P C is the critical axial force, G IC is a mode I crack, and C 3, K, D 11 , D 22 , and D are the stiffness coefficients of undrilled materials. Table 1. Summary of P C Prediction Models.

3−(1−υ 2 )
Lachuad [6] P C = 8π G IC D 1/3−(D /8D) Jain [7] P C = 3π D 22 Saoudi [9] P C = π(K * +G IC ) ξ(C 3 /3−K) Through the above analysis of the P C model, it can be found that the "fracture toughness" is an important factor that affects the P C . However, due to the complexity of the CFRP drilling process, it is difficult to directly observe the failure mode of drilling delamination from the experimental results. Most of the above models assume that the failure mode of drilling delamination crack is the mode I crack failure form. In fact, due to the complexity of force in the drilling process of composite materials, the failure mode is often a combination of different failure modes [1], and thus, using the mode I crack alone cannot reflect the actual cracking situation. In addition, for CFRP, because the temperature has a great impact on the mechanical properties of the resin matrix, some mechanical properties related to the resin, especially the fracture toughness G IC or G IIC , will change greatly with an increase in temperature, but the above models do not consider the influence of the temperature-dependent characteristics of CFRP mechanical properties on P C . This research found that G IC , C 3 , K, D 11 , D 22 , and D in the classical calculation equation of P C are affected by temperature changes and that the influence trends are different. The drilling process of composite materials produces a lot of heat, which increases the drilling temperature. Therefore, the influence of drilling temperature on the P C should also be considered when analyzing the P C of CFRP drilling delamination. Therefore, we established a three-dimensional finite element drilling model of CFRP, including an interface phase, to study the interlaminate damage mechanism of CFRP drilling and obtain the interlaminate damage mode of CFRP caused by drilling; then, based on the P C model of Zhang [8], we deduced a new P C model that considers the temperaturedependent characteristics of CFRP mechanical properties and the damage mode of CFRP interlaminate; finally, the P C model was verified by static compression experiments at different temperatures.

Material Constitutive Definition
CFRPs are defined as anisotropic linear elastic constitutive materials, and the failure criterion adopts the three-dimensional Hashin failure criterion [10,11]. As the damage criterion that comes with Abaqus does not have a three-dimensional Hashin failure criterion, the FORTRAN language must be used to write the three-dimensional Hashin failure criterion and damage evolution into the user material subprogram VUMAT for analysis and calculation. The interlamination interface is established by cohesive elements. Because the interlamination thickness of the composite material is very thin, the zero-thickness cohesive element is used to model it in the simulation analysis. The damage criterion adopts the traction-separation law with the I/II/III mixed modes [12][13][14]. The initial failure mode adopts the secondary stress failure mode, and the damage evolution process adopts the energy-based Benzeggagh-Kenane (BK) damage evolution criterion [15]. The damage evolution process after cracking adopts the energy-based BK damage evolution criterion: where η is the energy index; G IC , G IIC , and G IIIC are the fracture toughness (N/mm) of the three cracking modes, respectively, and G I , G II , and G III are the energy release rates during calculation (N/mm).

Finite Element Modeling and Parameter Setting
A carbide twist drill with a diameter of 8 mm is introduced into Abaqus, and the elastic modulus of the drill material is set to 640,000 Mpa and the Poisson's ratio to 0.22. As the elastic modulus of the tool material is much larger than that of CFRP, to save calculations during the analysis, we set the drill as a rigid body. To define the material properties and directional properties of the composite material, CFRP adopts T300/epoxy resin, the fiber volume fraction is about 62% according to Hengshen Co. Ltd., Zhenjiang, China, the layering method is [0/45/90/−45]s, and the thickness of a single layer is 0.125 mm, totaling 1 mm. The element type is 3D stress. The material properties are shown in Table 2, which were provided by Hengshen Co. Ltd., China. The interlamination interface adopts a zero-thickness cohesive element, and the material properties are shown in Table 3.
During the simulation calculation, the X-axis, Y-axis, and Z-axis movement freedom around the CFRP plate is constrained, and the downward feed of the drill bit and the rotational speed of the spindle rotating around the drill bit are defined. To reduce the calculation time and increase the degree of delamination damage, a larger feed is selected in the simulation. The feed rate and spindle speed are 0.27 mm/r and 9000 rpm, respectively. The downward velocity along the Y-axis is set to 40.5 mm/s. The rotation speed around the Y-axis is set to 942 rad/s. At the same time, the movement freedom of the rigid body reference point of the drill's X-axis and Z-axis is constrained. The cutting edge of the drill is in contact with the composite material. Due to the anisotropy of CFRP, its friction coefficient will vary with the cutting angle, ranging from 0.3 to 0.8 [16]. As the purpose of finite element analysis in this study is to observe the three-dimensional drilling process, the intermediate value of 0.5 is selected as the friction coefficient between the drill bit and the composite material. The established CFRP three-dimensional drilling finite element model is shown in Figure 1.  During the simulation calculation, the X-axis, Y-axis, and Z-axis movement freedom around the CFRP plate is constrained, and the downward feed of the drill bit and the rotational speed of the spindle rotating around the drill bit are defined. To reduce the calculation time and increase the degree of delamination damage, a larger feed is selected in the simulation. The feed rate and spindle speed are 0.27 mm/r and 9000 rpm, respectively. The downward velocity along the Y-axis is set to 40.5 mm/s. The rotation speed around the Y-axis is set to 942 rad/s. At the same time, the movement freedom of the rigid body reference point of the drill's X-axis and Z-axis is constrained. The cutting edge of the drill is in contact with the composite material. Due to the anisotropy of CFRP, its friction coefficient will vary with the cutting angle, ranging from 0.3 to 0.8 [16]. As the purpose of finite element analysis in this study is to observe the three-dimensional drilling process, the intermediate value of 0.5 is selected as the friction coefficient between the drill bit and the composite material. The established CFRP three-dimensional drilling finite element model is shown in Figure 1.    Figure 2 shows the simulation results of CFRP drilling. The white part is the CFRP element, the blue part is the interlamination interface element, and the red part represents the damage degree of the interface layer. To more clearly show the position and failure mode of the interlamination element in the laminated plate during the drilling process, some composite elements are hidden in Figure 2. Figure 2 shows the simulation results of CFRP drilling. The wh element, the blue part is the interlamination interface element, and sents the damage degree of the interface layer. To more clearly sho failure mode of the interlamination element in the laminated plate process, some composite elements are hidden in Figure 2. In order to more clearly show the damage and failure process terface elements during the drilling process, part of the CFRP eleme terlaminate element layer of the middle part are hidden in Figures 3 first layer of CFRP, the first layer of interlaminate interface element terface), and the last layer of interlaminate interface elements (short f reserved for location reference. Figure 3 is a comprehensive diagram is a representative figure of the drill bit passing through every layer layer of the interface at different stages. Figure 4a1-d1 is the front v a4), (b2-b4), (c2-c4), and (d3,d4), respectively. In order to more clearly show the damage and failure process of interlaminate interface elements during the drilling process, part of the CFRP element layer and the interlaminate element layer of the middle part are hidden in Figures 3 and 4, and only the first layer of CFRP, the first layer of interlaminate interface elements (short for first interface), and the last layer of interlaminate interface elements (short for last interface) are reserved for location reference. Figure 3 is a comprehensive diagram of Figure 4.       It can be seen from Figure 4(a1-a4) that, when the drill bit is drilled to stage I, the chisel edge makes contact with the first layer of the CFRP. Although the removal amount of CFRP material in the first layer is very small Figure 4(a2), the interlaminar element in the first layer also begins to be damaged due to the downward bending load. When the drill bit is drilled to stage II Figure 4(b1-b4), the main cutting edge starts to work, the contact delamination area between the first layer of the CFRP and the drill increases Figure 4(b2), and the damage delamination area of the first interface increases Figure 4(b3). At this time, the last interface is still not affected. When the drill bit is drilled to stage III Figure 4(c1-c4), the main cutting edge starts working, the drilling axial force increases, the contact delamination area between the first layer of the CFRP and the drill increases Figure 4(c2), and the damage delamination area of the first interface increases Figure 4(c3). At this time, the last interface begins to be damaged Figure 4(c4). When the drill bit is drilled to stage IV Figure 4(d1-d4), the main cutting edge is fully working, and the drilling axial force increases to its maximum. Accordingly, the damage delamination area of both the first layer of the CFRP and the first interface increase to their maxima Figure 4(d2,d3). As the drill bit is not in contact with the last interface, the damage evolution area of the last interface at this stage does not expand significantly. The above results show that the [0/45/90/−45]s layers can cover all the drilling processes of this bit, and all CFRP layers and interface layers need to go through stages I to IV.

Analysis of Simulation Results
Due to the laminated structure characteristics of CFRP, the crack form often occurs between layers during drilling, accompanied by a coupling effect. To judge the cracking form of the laminated plate in the drilling process, it is necessary to analyze the interlamination failure form in FEM results combined with different drilling stages. Figure 5 shows the section view of the finite element simulation results when the chisel edge just contacts the material (drilling stage I). It can be seen from Figure 5 that due to the width of the chisel edge, it will exert a pressing force on the vertical feeding direction of the first layer of the CFRP material of the laminate, which will cause the material of this layer to produce a shear force (κ) in the vertical extrusion direction and then cause a type II crack. At the same time, the downward feeding movement will produce a downward extrusion force (σ) on the uncut layer and then cause a type I crack. As the shear force and extrusion force exist simultaneously, the CFRP at the entrance is affected by the chisel edge, resulting in the coupling of I/II cracks ( Figure 5). lution area of the last interface at this stage does not expand sults show that the [0/45/90/−45]s layers can cover all the dril all CFRP layers and interface layers need to go through stage Due to the laminated structure characteristics of CFRP between layers during drilling, accompanied by a coupling form of the laminated plate in the drilling process, it is ne lamination failure form in FEM results combined with diffe shows the section view of the finite element simulation resu contacts the material (drilling stage I). It can be seen from Fi of the chisel edge, it will exert a pressing force on the vertica layer of the CFRP material of the laminate, which will cause produce a shear force (κ) in the vertical extrusion direction an At the same time, the downward feeding movement will pro force (σ) on the uncut layer and then cause a type I crack. A sion force exist simultaneously, the CFRP at the entrance is resulting in the coupling of I/II cracks ( Figure 5). When the drill bit continues to feed along the axial dire begins to work (drilling stage II). Any point on the main shear force (κ) in the direction of the cutting speed, which w the same time, the chisel edge feeds downwards, and the do will produce a downward extrusion force (σ) on the uncut When the drill bit continues to feed along the axial direction, the main cutting edge begins to work (drilling stage II). Any point on the main cutting edge will generate a shear force (κ) in the direction of the cutting speed, which will form a type III crack. At the same time, the chisel edge feeds downwards, and the downward feeding movement will produce a downward extrusion force (σ) on the uncut layer, causing a type I crack. Because the shear force in the direction of the cutting speed and the downward extrusion force exist simultaneously, this results in the coupling of I/III cracks ( Figure 6). When the main cutting edge is fully working (drilling stage IV), t moval form and the delamination form caused by the chisel edge and th edge are the same, type I/II cracks are caused by the chisel edge, and type caused by the main cutting edge. However, when drilling to the exit, as uncut layers decreases, the residual stiffness also decreases, and the I/II cr the chisel edge are removed by the movement of the main cutting edg shear force (κ) in the direction of the cutting speed and the downward ext exist simultaneously, I/III cracks couple at the drilling exit. Therefore, th nation form of the drilling exit is the I/III mixed crack ( Figure 7).

Modeling of CFRP Drilling Delamination PC
In this section, based on the drilling axial force model of Zhang [8 ture-dependent characteristics of CFRP mechanical properties and the lamination failure mode at the drilling exit, a prediction model of the dril lished. The model derivation process is as follows: When the main cutting edge is fully working (drilling stage IV), the material removal form and the delamination form caused by the chisel edge and the main cutting edge are the same, type I/II cracks are caused by the chisel edge, and type I/III cracks are caused by the main cutting edge. However, when drilling to the exit, as the number of uncut layers decreases, the residual stiffness also decreases, and the I/II cracks caused by the chisel edge are removed by the movement of the main cutting edge. Because the shear force (κ) in the direction of the cutting speed and the downward extrusion force (σ) exist simultaneously, I/III cracks couple at the drilling exit. Therefore, the main delamination form of the drilling exit is the I/III mixed crack (Figure 7). Because the shear force in the direction of the cutting speed and the downward extrusion force exist simultaneously, this results in the coupling of I/III cracks ( Figure 6). When the main cutting edge is fully working (drilling stage IV), the material removal form and the delamination form caused by the chisel edge and the main cutting edge are the same, type I/II cracks are caused by the chisel edge, and type I/III cracks are caused by the main cutting edge. However, when drilling to the exit, as the number of uncut layers decreases, the residual stiffness also decreases, and the I/II cracks caused by the chisel edge are removed by the movement of the main cutting edge. Because the shear force (κ) in the direction of the cutting speed and the downward extrusion force (σ) exist simultaneously, I/III cracks couple at the drilling exit. Therefore, the main delamination form of the drilling exit is the I/III mixed crack (Figure 7).

Modeling of CFRP Drilling Delamination PC
In this section, based on the drilling axial force model of Zhang [8], the temperature-dependent characteristics of CFRP mechanical properties and the mixed I/III delamination failure mode at the drilling exit, a prediction model of the drilling PC is established. The model derivation process is as follows: First, assuming that the shape of the CFRP drilling exit layer is elliptical [5], the force state of the delamination damage at the drilling exit is shown in Figure 8. In Figure 8a, h is the thickness of the undrilled layers, a is the longitudinal radius of the ideal damaged ellipse, and b is the transverse radius of the ideal damaged ellipse (Figure 8b).

Modeling of CFRP Drilling Delamination P C
In this section, based on the drilling axial force model of Zhang [8], the temperaturedependent characteristics of CFRP mechanical properties and the mixed I/III delamination failure mode at the drilling exit, a prediction model of the drilling P C is established. The model derivation process is as follows: First, assuming that the shape of the CFRP drilling exit layer is elliptical [5], the force state of the delamination damage at the drilling exit is shown in Figure 8. In Figure 8a, h is the thickness of the undrilled layers, a is the longitudinal radius of the ideal damaged ellipse, and b is the transverse radius of the ideal damaged ellipse (Figure 8b). Assuming that all the work accomplished via the drilling axial force at the drilling exit ( Figure 8) is converted into energy released by delamination and energy required by material strain, then, based on the linear elastic fracture mechanics and the law of energy conservation, the balance relation equation of energy required for delamination can be established: where Pc is the critical axial force (N); dω0 is the differential cross-section deflection of uncut laminates; GC is the critical energy release rate (N/mm); dA is the differential layered area (mm); and dU is the differential of the strain energy required for the elastic strain of the material (N). The process of solving the Pc is also the process of solving dω0, GC, dA, and dU in the equation. In previous studies, the influence of cutting temperature has not been considered when solving the above parameters. This study analyzes the equation based on the temperature-dependent characteristics of CFRPs' mechanical properties. The detailed process is as follows.

The Solution of Section Deflection dw0 of Undrilled Layer Material
According to the classical theory of laminated plates, the constitutive relationship of CFRPs considering the mechanical and thermal coupling properties is shown in Equation    Assuming that all the work accomplished via the drilling axial force at the drilling exit ( Figure 8) is converted into energy released by delamination and energy required by material strain, then, based on the linear elastic fracture mechanics and the law of energy conservation, the balance relation equation of energy required for delamination can be established: where Pc is the critical axial force (N); dω 0 is the differential cross-section deflection of uncut laminates; G C is the critical energy release rate (N/mm); dA is the differential layered area (mm); and dU is the differential of the strain energy required for the elastic strain of the material (N). The process of solving the Pc is also the process of solving dω 0 , G C , dA, and dU in the equation. In previous studies, the influence of cutting temperature has not been considered when solving the above parameters. This study analyzes the equation based on the temperature-dependent characteristics of CFRPs' mechanical properties. The detailed process is as follows.

The Solution of Section Deflection dw 0 of Undrilled Layer Material
According to the classical theory of laminated plates, the constitutive relationship of CFRPs considering the mechanical and thermal coupling properties is shown in Equation (3).   [8].
Assuming that the displacement is very small during drilling, the relationship between strain, curvature, and displacement is: where ε is the in-plane strain; u and v are the displacements in the x and y directions, respectively; τ is the curvature of the laminate midplane; and ω is the displacement perpendicular to the laminate. The drilling force model of undrilled CFRP laminate can be simplified as a plate model with a concentrated force in the middle. According to the basic bending equilibrium differential equation of plate shell theory: ∂N y ∂y Assuming that the load is uniformly distributed on the surface of the failure zone: q = P C πab = ξ πa 2 P C , where ξ = a b . The relationship between the internal force, bending moment, and displacement can be obtained by substituting the strain component and stiffness matrix of [A i ], [B i ], and [D i ] into Equation (3). Then, the internal force and bending moment are substituted into Equation (12) to obtain the relationship between the displacement and P C (Equations (13)- (17)). C j is the material performance coefficient of the undrilled material.
The deflection of the undrilled layer material is obtained by deriving the diameter of the elliptical long axis in Equation (15):

Solution of Elastic Strain Energy dU of Undrilled Material
According to classical plate theory, the strain energy of the plate model can be obtained via the following equation: Here, the relationship between stress and strain under thermal-mechanical coupling is: where ε x is the total strain ε m x is the strain caused by mechanical force; and α x is the coefficient of thermal expansion.
Equation (21) is obtained by expanding Equation (20) and substituting the expanded Equation (20) into Equation (19): Assuming that the temperature difference between the drilling temperature and the drilling ambient temperature is 0, ε t = 0, and the strain energy generated by pure mechanical strain is: Pure mechanical strain can be obtained by substituting Equations (13)-(17) into Equations (4) and (9): By substituting Equations (23)-(28) into Equation (22), the strain energy of undrilled laminates caused by mechanical force can be obtained when the drilling temperature difference is 0: where: Referring to the solution of thermal strain in Saoudi [9], when the temperature difference between the drilling temperature and the drilling ambient temperature is not zero, the total strain energy including thermal strain caused by the undrilled laminate is as follows: The basic assumptions of classical laminated plate theory are: where: where ∆T is the temperature difference between the drilling temperature and the ambient temperature.
Then, the calculation equation of undrilled materials' strain energy under mechanicalthermal coupling can be deduced: From the derivation of Equation (39), we can obtain:

Solution of Fracture Randomness
Through the analysis of the results of a three-dimensional drilling finite element analysis, it can be seen that, when the drill bit is drilled at the exit, the delamination failure mode is type I/III mixed mode, so the coupling of I and III cracks should be considered in the analytical model of P C . According to the BK damage criterion: Then, where r is the mixing coefficient of the fracture toughness of type I/III cracks and r = ( G I I I G I +G I I I ) η . This can be solved via the conjugate gradient method according to inverse problem theory.
When substituting the Gc, dA, dw, and dU obtained above into the critical layered energy balance relationship, the analytical equation of the P C can be obtained as follows: The temperature-dependent parameters in the equation are G IC , G IIIC , K, C j , and D*. The mechanical parameters of G IC that vary with temperature are obtained via the ASTM D5528-01 standard [17], the experimental results of G IC are shown in Figure 9, and a detailed experimental process is shown in [18]; the mechanical parameters of G IIIC that vary with temperature are obtained from the measurement results of G IIC [19], the mechanical parameters of G IIC that vary with temperature are obtained via the ASTM D7905M-14 standard [20], the experimental results of G IIC are shown in Figure 10, and the detailed experimental process is shown in Appendix A; the values of K are calculated via Equation (30); the calculation methods of C j are from [8].
(1 ) where r is the mixing coefficient of the fracture toughness of type I/III crac . This can be solved via the conjugate gradient method according verse problem theory. When substituting the Gc, dA, dw, and dU obtained above into the critical energy balance relationship, the analytical equation of the PC can be obtained as fo The temperature-dependent parameters in the equation are GIC, GIIIC, K, Cj, a The mechanical parameters of GIC that vary with temperature are obtained via the D5528-01 standard [17], the experimental results of GIC are shown in Figure 9, an tailed experimental process is shown in [18]; the mechanical parameters of GIIIC th with temperature are obtained from the measurement results of GIIC [19], the mec parameters of GIIC that vary with temperature are obtained via the ASTM D79 standard [20], the experimental results of GIIC are shown in Figure 10, and the d experimental process is shown in Appendix A; the values of K are calculated via tion (30); the calculation methods of Cj are from [8].
Both K and Cj are calculated by the combination of the [Ai], [Bi], and [Di] s matrices. As the mechanical properties of carbon fibers hardly change with a ch drilling temperature in this temperature range, the parameters affecting the value [Ai], [Bi], and [Di] stiffness matrices are calculated using the modulus of the resin with temperature, which is obtained via the tensile test of the resin at different a temperatures ( Figure 11). The detailed experimental process is shown in [18]; in a to the stiffness matrix of [Ai], [Bi], and [Di], the value of D* at different tempera also affected by the thermal expansion coefficient. The thermal expansion coeffic the unidirectional CFRP were measured using a dilatometer (NETZSCH DIL 402C).

The Solution of the Layered Area
Assuming that the shape of the exit layer is an ellipse, the area of the layer is the area of the ellipse minus the area of the ideal hole. The assumptions are:

The Solution of the Layered Area
Assuming that the shape of the exit layer is an ellipse, the area of the layer is the area of the ellipse minus the area of the ideal hole. The assumptions are: where ξ is an ellipse ratio. Then, the differential of the stratified area is:

Experimental Setup
To verify the accuracy of the drilling P C model derived above, a temperature-controlled P C experiment must be carried out. However, due to the influence of the process parameters on the drilling temperature and drilling force, it is difficult to obtain a P C that only produces delamination at different drilling temperatures through drilling tests. In addition, the experimental data from past studies on drilling P C cannot be referenced because the influence of drilling temperature is not considered, but the experimental method can be referenced. In previous studies, a static compression experiment was used to measure the P C .
To verify the theoretical model introduced in the previous section, this section draws on this test method to design a P C equivalent test at different ambient temperatures. However, because the P C experiment can only approximately simulate the axial force generated by the downward feed of the drill bit and cannot simulate the torque caused by the rotation of the drill bit, the failure mode of the exit stratification of the P C experiment is the only mode I open failure.
The CFRP composite used in this study was a carbon T300/epoxy unidirectional prepreg with a ply thickness of 0.125 mm. The paving sequence of materials used in the experiment is [0/45/90/−45] 4S . The reason why we use [0/45/90/−45] 4S to verify the analytical mode is that a thick layer can form better-quality blind holes. After laying the plate preform manually, we put it into the autoclave for heating and curing. The curing conditions are heating to 80 • C, holding for 30 min, then pressurizing to 0.5 MPa, heating to 120 • C, holding for 90 min, and finally cooling in a furnace. During the experiment, non-drilled conical blind holes with a thickness of 1-6 layers are prefabricated on the tested sample. Due to the taper of the drill, when the number of layers is small, the horizontal edge of the drill drills the blind holes. To ensure that drilling stratification does not occur during the prefabrication of blind holes, a back plate is placed under the laminated plate to be drilled during drilling to increase the exit stiffness. The spindle speed is 4000 rpm, and the feed rate is 0.03 mm/r. Figure 13 shows the P C testing bench at different ambient temperatures. During the experiment, the sample was fixed on the platform of the universal testing machine. The exit temperature of the blind hole was locally heated by a silica gel heater. The temperature of the heater was set to 23 • C, 60 • C, 90 • C, and 120 • C, respectively. The experimental pressure head was a carbide drill bit with a diameter of 8 mm. The compression speed was set to 2 mm/min during the test.  Figure 14 shows a comparison between the PC results predicted by Equation (45) and the experimental results under different drilling ambient temperatures. It can be seen from Figure 14 that the theoretically predicted PC value is in good agreement with the experimental value when the drilling temperature is 120 °C, and there are some errors when the drilling temperature is 23 °C, 60 °C, or 90 °C. There are two reasons for this error. Firstly, the influence of the geometry of the drill bit on the drilling axial force distribution is not considered in this PC model. Secondly, when measuring the fracture toughness, the influence of different ply angles is not considered.

Result Discussion
However, both the prediction results and the experimental results show that the drilling temperature has a great impact on the critical axial force, and the critical axial force increases with an increase in the drilling temperature. Additionally, the change trend of the curve is consistent when the drilling temperature is 23 °C, 60 °C, and 90 °C. From the general trend, it can be seen that the PC of CFRP drilling delamination is greatly influenced by the mechanical properties of the CFRP. When the drilling temperature is not greater than the glass transition temperature range of the material itself (about 120 °C), due to the increase in fracture toughness, the PC increases with an increase in drilling temperature.  Figure 14 shows a comparison between the P C results predicted by Equation (45) and the experimental results under different drilling ambient temperatures. It can be seen from Figure 14 that the theoretically predicted P C value is in good agreement with the experimental value when the drilling temperature is 120 • C, and there are some errors when the drilling temperature is 23 • C, 60 • C, or 90 • C. There are two reasons for this error. Firstly, the influence of the geometry of the drill bit on the drilling axial force distribution is not considered in this P C model. Secondly, when measuring the fracture toughness, the influence of different ply angles is not considered.   Figure 14 shows a comparison between the PC results predicted by Equation (45) and the experimental results under different drilling ambient temperatures. It can be seen from Figure 14 that the theoretically predicted PC value is in good agreement with the experimental value when the drilling temperature is 120 °C, and there are some errors when the drilling temperature is 23 °C, 60 °C, or 90 °C. There are two reasons for this error. Firstly, the influence of the geometry of the drill bit on the drilling axial force distribution is not considered in this PC model. Secondly, when measuring the fracture toughness, the influence of different ply angles is not considered.

Result Discussion
However, both the prediction results and the experimental results show that the drilling temperature has a great impact on the critical axial force, and the critical axial force increases with an increase in the drilling temperature. Additionally, the change trend of the curve is consistent when the drilling temperature is 23 °C, 60 °C, and 90 °C. From the general trend, it can be seen that the PC of CFRP drilling delamination is greatly influenced by the mechanical properties of the CFRP. When the drilling temperature is not greater than the glass transition temperature range of the material itself (about 120 °C), due to the increase in fracture toughness, the PC increases with an increase in drilling temperature.  However, both the prediction results and the experimental results show that the drilling temperature has a great impact on the critical axial force, and the critical axial force increases with an increase in the drilling temperature. Additionally, the change trend of the curve is consistent when the drilling temperature is 23 • C, 60 • C, and 90 • C. From the general trend, it can be seen that the P C of CFRP drilling delamination is greatly influenced by the mechanical properties of the CFRP. When the drilling temperature is not greater than the glass transition temperature range of the material itself (about 120 • C), due to the increase in fracture toughness, the P C increases with an increase in drilling temperature.

Conclusions
This paper focuses on the influence of the temperature-dependent characteristics of CFRPs' mechanical properties on the P C of drilling delamination damage. First, the failure mode of drilling delamination is qualitatively analyzed through a three-dimensional finite element model. Second, a new drilling P C model, which considers the temperaturedependence of CFRPs' mechanical properties and the failure modes of CFRP drilling delamination, is established based on the classical drilling critical force mechanics model. The influence of the change in CFRPs' mechanical properties on the P C under different drilling temperatures is analyzed, and the model is verified at different temperatures. The main findings are as follows: 1.
Through the finite element simulation results, it can be found that the failure mode of drilling exit delamination is mainly the I/III mixed crack failure mode; 2.
Due to the influence of the temperature-dependent characteristics of CFRP mechanical properties, when the drilling temperature changes, the P C is not a fixed value.
The drilling temperature has a great impact on P C when the temperature is lower than the glass transition temperature. The fracture toughness increases with an increase in drilling temperature, and the value of P C increases with an increase in drilling temperature.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.

Conflicts of Interest:
The authors declare no conflict of interest.

Mode II Interlaminar Fracture Toughness Testing
The influence of temperature on G IIC was tested using the three-point bending test, according to the standard of ASTM D7905M-14 [20]. The sample dimensions were 160 mm × 22 mm × 3.5 mm. The paving method is [0] 28 , and the initial delamination length was 45 mm. All tests were performed at a speed of 2 mm/min. Set the test temperature as 23 • C, 60 • C, 90 • C, and 120 • C, respectively. Five tests were performed at each temperature. Detailed experiments and sample settings are shown in Appendix A Figure A1. where p is the load (N); δ is the load point displacement (mm); b is the specimen width (mm); a is the delamination length (mm); and L is half the span (mm). Combined with the force/displacement curve and Formula (2)  The test samples measured at different ambient temperatures are shown in Appendix A Figure A2. The force/displacement curves measured at different ambient temperatures are shown in Appendix A Figure A3. The following equation was used to calculate G IIC : where p is the load (N); δ is the load point displacement (mm); b is the specimen width (mm); a is the delamination length (mm); and L is half the span (mm). Combined with the force/displacement curve and Formula (2), the Mode II interlaminar fracture toughness at different ambient temperatures can be calculated. The results are shown in Figure 10. where p is the load (N); δ is the load point displacement (mm); b is the specimen width (mm); a is the delamination length (mm); and L is half the span (mm). Combined with the force/displacement curve and Formula (2)  where p is the load (N); δ is the load point displacement (mm); b is the specimen width (mm); a is the delamination length (mm); and L is half the span (mm). Combined with the force/displacement curve and Formula (2)