Numerical Simulation of Rotary Ultrasonic Machining of the Nomex Honeycomb Composite Structure

: Nomex honeycomb composite (NHC) cores have seen significant growth in recent years, particularly in the aeronautics, aerospace, naval and automotive industries. This development presents significant challenges in terms of improving machining quality, requiring the use of specialized cutting tools and favorable cutting techniques. In this context, experimental studies have been carried out to highlight the characteristics of the milling of NHCs by rotary ultrasonic machining (RUM). However, the rapid motion of the cutting tool and the inaccessibility of the tool/part interface prevent the visualization of the chip formation process. For this purpose, a three-dimensional numerical model for milling the NHC structure using RUM technology was developed by Abaqus Explicit software. On the basis of this model, the components of the cutting force, the quality of the machined surface and the chip accumulation in front of the cutting tool were analyzed. The numerical results agree with the experimental tests, demonstrating that the use of RUM technology effectively reduces the cutting force components. An in-depth analysis of the influence of feed component F y on the quality of the generated surface was carried out, revealing that the surface quality improved with low values of feed component F y . Furthermore, the impact of ultrasonic vibrations on the accumulation of chips in front of the cutting tool is particularly optimized, in particular for large amplitudes.


Introduction
A Nomex honeycomb composite (NHC) core is made from aramid fibers impregnated with phenolic resin.This combination gives Nomex honeycomb composites superior mechanical strength, increased toughness and excellent chemical resistance [1][2][3].Due to their remarkable mechanical characteristics, NHC cores are widely used in the fields of aerospace and aeronautics, mainly in the manufacture of aircraft wings and tails [4][5][6][7].However, the machining of NHC structures represents a major challenge for researchers and engineers due to the geometric complexity of thin-walled structures and the heterogeneous characteristics of Nomex paper.To meet these challenges, in-depth experimental studies are needed to understand and analyze the machining characteristics of NHC structures.These studies allow the analysis of the cutting process, while providing specific data to develop advanced numerical models aimed at analyzing cutting behavior and optimizing machining conditions [8,9].An NHC core is generally machined using conventional methods, involving the use of cutting tools such as mills and lathes to remove material chips.However, despite the optimization of specially designed machining tools, the conventional machining process has significant shortcomings in terms of machining accuracy and the quality of the machined surface.In the fact, the problems with the machined surface of the NHC core manifest themselves in the form of burrs, cell deformation and the tearing of aramid fibers [10].Consequently, conventional machining fails to deal effectively with these peculiarities, resulting in a poor machined surface.To overcome these problems, it is necessary to develop more advanced machining methods adapted specifically to the NHC structure, as well as to optimize the cutting conditions to minimize cell deformation and fiber pull-out.In response to these challenges, RUM technology was proposed to overcome the limitations of conventional machining process when processing NHCs.It is distinguished by the integration of ultrasonic vibrations along the axis of the tool, thus allowing better management of cell deformation and optimization of the cutting of the aramid fiber [11][12][13][14].RUM technology has gained increasing interest in the machining of a variety of conductive and dielectric materials.These materials include titanium alloys [15,16], aluminum alloys [17][18][19][20], nickel alloys [21,22], glass [23,24], ceramics [25,26] and composites [27][28][29][30].A study by Xia et al. [31] aimed to improve the structural design of an ultrasonic circular disc tool and to evaluate the effect of its geometrical characteristics on the energy density and stiffness of the tool.Their results indicated that the effectiveness of the NHC core treatment is highly dependent on the cutting angle, tool radius and vibration amplitude.Research carried out by Sun et al. [32] highlighted the advantages of applying RUM technology when machining aluminum honeycomb structures.In this regard, the authors observed that the use of this technology led to a reduction in the forces exerted by the disc cutter, thus reducing the plastic deformation of the structure's walls.A study discussed by Ahmed et al. [33] focused on the analysis of various machining methods applied to NHC cores.The main aim of this research was to compare and evaluate the performance of three distinct machining approaches, conventional machining, rotary ultrasonic machining (RUM) and laser machining.The implementation of an experimental setup entails considerable costs, which gives numerical modeling the status of a reliable tool capable of complementing experimental tests.The choice of the constitutive law represents a major challenge for numerical modeling, given that the simulation of the walls of the NHC structure is mainly based on the use of elastoplastic constitutive laws, which can be either isotropic or orthotropic.The isotropic elastoplastic approach is generally simpler to implement, but it tends to neglect the composite architecture of the material.On the other hand, the orthotropic approach offers increased precision by taking into account the specific mechanical properties linked to the orientation of the fibers.Although experimental testing remains the primary method used to machine NHC cores, research on numerical models is evolving to enrich and deepen the understanding of this process.Thus, numerical modeling offers the possibility to take into account conditions that are difficult to measure experimentally, to study extreme conditions and to deepen our understanding of them.It should be emphasized that the development of numerical milling models of NHC structures assisted by RUM technology remains a relatively recent area of research.Until now, few numerical studies have been carried out in the literature to analyze the characteristics of the NHC core milling process using a UCK cutting tool.To achieve this objective, a 3D numerical approach toward the milling process of NHC structures, using RUM technology, was developed by Abaqus Explicit software based on the finite element method.The numerical model was validated by experimental tests relating to the influence of the spindle speed, the feed rate and the amplitude of ultrasonic vibrations on the components of the cutting force and the quality of the machined surface.Subsequently, an in-depth analysis was undertaken to examine the impact of vibration amplitude on chip accumulation in front of the cutting tool.The comparison between the results obtained by simulation and those resulting from experimentation validate the reliability of the numerical model developed.

Design of the Cutting Tool and the Structure
Experimental studies were carried out on the milling of NHC structures using RUM technology.These studies were carried out using the ultrasonic machine tool (THU Ultrasonic 850) developed by Tsinghua University in China (see Figure 1) [34].The experimental setup included the following components: BP4610 high-speed bipolar ultrasonic power supply, a high-frequency compensation capacitor, an ultrasonic spindle system with a BT40compliant tool holder, a KISTLER-9256C2 dynamometer, a control panel for the machine tool, a fixture and a clamping unit to hold the test piece.These components were used together to collect the necessary experimental data when milling the NHC structure.
undertaken to examine the impact of vibration amplitude on chip accumulation in front of the cu ing tool.The comparison between the results obtained by simulation and those resulting from experimentation validate the reliability of the numerical model developed.

Design of the Cu ing Tool and the Structure
Experimental studies were carried out on the milling of NHC structures using RUM technology.These studies were carried out using the ultrasonic machine tool (THU Ultrasonic 850) developed by Tsinghua University in China (see Figure 1) [34].The experimental setup included the following components: BP4610 high-speed bipolar ultrasonic power supply, a high-frequency compensation capacitor, an ultrasonic spindle system with a BT40-compliant tool holder, a KISTLER-9256C2 dynamometer, a control panel for the machine tool, a fixture and a clamping unit to hold the test piece.These components were used together to collect the necessary experimental data when milling the NHC structure.The ultrasonic machine tool (THU Ultrasonic 850) has key performance conditions summarized in Table 1.However, se ing up the experimental procedure remains a lengthy and costly process.Thus, numerical modeling is a very interesting tool for simulating machining problems in various configurations, under different cu ing conditions and at lower cost.For this purpose, a 3D numerical model of the milling process of NHC The ultrasonic machine tool (THU Ultrasonic 850) has key performance conditions summarized in Table 1.However, setting up the experimental procedure remains a lengthy and costly process.Thus, numerical modeling is a very interesting tool for simulating machining problems in various configurations, under different cutting conditions and at lower cost.For this purpose, a 3D numerical model of the milling process of NHC structures, using RUM technology, was developed by Abaqus Explicit software.In this paper, the NHC structure studied was composed of single walls and of double walls represented in the W and L directions, respectively, as shown in Figure 2a.The geometric dimensions of the hexagonal cell are described in Figure 2b [34].structures, using RUM technology, was developed by Abaqus Explicit software.In this paper, the NHC structure studied was composed of single walls and of double walls represented in the W and L directions, respectively, as shown in Figure 2a.The geometric dimensions of the hexagonal cell are described in Figure 2b [34].
Table 1.Performance conditions of the THU Ultrasonic 850.According to the experimental procedure, the NHC core was machined using a specific cu ing tool called UCK, which is a circular milling cu er made of ultrasonic highspeed steel (HSS-W18Cr4V) (Figure 3b).For this purpose, a cu ing tool was designed respecting the geometric dimensions used in the experimental phase (see Figure 4) [34].According to the experimental procedure, the NHC core was machined using a specific cutting tool called UCK, which is a circular milling cutter made of ultrasonic high-speed steel (HSS-W18Cr4V) (Figure 3b).For this purpose, a cutting tool was designed respecting the geometric dimensions used in the experimental phase (see Figure 4) [34].

Presentation of the Numerical Model
The present work proposes a numerical study of the milling of an NHC structure using RUM technology, where the structure to be machined is composed of thin walls with low thickness (Figure 5a).In this regard, the walls of the NHC structure were meshed using four-node classic shell elements (S4R), with each node having six degrees of freedom (three translations and three rotations) with reduced integration.During the numerical simulation, the cu ing tool was considered rigid, which means that it did not undergo any deformation throughout the machining operation.In order to represent it correctly in the numerical model, it was discretized using rigid quadrangular elements with 4 nodes (R3D4) commonly used for three-dimensional analysis (Figure 5b).

Setting Up of the Numerical Model 2.2.1. Presentation of the Numerical Model
The present work proposes a numerical study of the milling of an NHC structure using RUM technology, where the structure to be machined is composed of thin walls with low thickness (Figure 5a).In this regard, the walls of the NHC structure were meshed using four-node classic shell elements (S4R), with each node having six degrees of freedom (three translations and three rotations) with reduced integration.During the numerical simulation, the cutting tool was considered rigid, which means that it did not undergo any deformation throughout the machining operation.In order to represent it correctly in the numerical model, it was discretized using rigid quadrangular elements with 4 nodes (R3D4) commonly used for three-dimensional analysis (Figure 5b).Before the implementation of the numerical code, a main step was carried out to analyze the nonlinear parameters likely to influence the convergence of the numerical results.Several factors needed to be taken into account during numerical simulation, including the size mesh, the type of contact, and the friction coefficient.It should be noted that using a small mesh size can lead to a significant increase in computation time (CPU), especially in three-dimensional setups.Thus, to find the right compromise between the size of the mesh and the precision of the results, a mesh of around 0.4 mm was adopted, corresponding to a total number of elements of 83,649.During the milling simulation, general contact was used to model the interaction between the cu ing tool and the walls of the NHC core.Due to the small thickness of the walls forming the NHC core, the contact be- Before the implementation of the numerical code, a main step was carried out to analyze the nonlinear parameters likely to influence the convergence of the numerical results.Several factors needed to be taken into account during numerical simulation, including the size mesh, the type of contact, and the friction coefficient.It should be noted that using a small mesh size can lead to a significant increase in computation time (CPU), especially in three-dimensional setups.Thus, to find the right compromise between the size of the mesh and the precision of the results, a mesh of around 0.4 mm was adopted, corresponding to a total number of elements of 83,649.During the milling simulation, general contact was used to model the interaction between the cutting tool and the walls of the NHC core.Due to the small thickness of the walls forming the NHC core, the contact between the cutting tool and the NHC structure was treated as a point contact, verified by Coulomb's law.In this work, the value of the contact coefficient chosen for the numerical simulations was 0.1.In order to ensure full contact between the cutting tool and the NHC structure from the start of the simulation, an initial engagement was planned, taking into account the particulars geometries of the structure and the cutting tool, thus reducing the CPU calculation time (Figure 5a).In the numerical simulation, the boundary conditions were defined in accordance with the established experimental procedure (Figure 6b).In this context, the bottom surface of the NHC structure remained stationary, indicating that it was unable to move or rotate in any direction.For this purpose, in the numerical simulation, we applied total embedding to the lower surface of the NHC core, which had the effect of eliminating all translations (U x = U y = U z = 0) along X, Y and Z and all rotations (U Rx = U Ry = U Rz = 0) around X, Y and Z (Figure 6a).Rotary ultrasonic machining (RUM) is a mechanical manufacturing process for material removal that relies on the coordination of three movements of the cutting tool: the feed motion along the OY axis, characterized by the feed rate V f ; the rotational motion around the OZ axis, characterized by a spindle speed n; and the vibratory motion along the OZ axis, applying a sinusoidal ultrasonic wave (Figure 6a).In order to follow the evolution of the machining precisely during the numerical simulation, a reference point, RP, was assigned on the axis of rotation of the cutting tool (Figure 5b).This reference point was essential for establishing cutting conditions and evaluating the cutting forces during the machining process.The UCK tool motion can be expressed in terms of global coordinates xyz by the following equations: Machines 2024, 12, x FOR PEER REVIEW 7 of 20

𝑉 = 𝑉 𝑠𝑖𝑛 𝑡
(2) represents ultrasonic vibration amplitude,  is cu ing speed,  is feed rate, n is spindle speed,  is time and  is the vibration frequency which is 21.260 KHz in the present paper.

Law of Behavior and Criterion of Damage
The adoption of an appropriate and representative constitutive law remains essential in order to effectively simulate the chip formation process during machining.In this paper, the Nomex paper forming the NHC structure exhibited bri le elastic mechanical behavior, meaning that it reacted differently along the stress axes.This feature resulted in high tensile strength but relatively lower compressive strength [35,36].In the context of A represents ultrasonic vibration amplitude, V c is cutting speed, V f is feed rate, n is spindle speed, t is time and f is the vibration frequency which is 21.260 KHz in the present paper.

Law of Behavior and Criterion of Damage
The adoption of an appropriate and representative constitutive law remains essential in order to effectively simulate the chip formation process during machining.In this paper, the Nomex paper forming the NHC structure exhibited brittle elastic mechanical behavior, meaning that it reacted differently along the stress axes.This feature resulted in high tensile strength but relatively lower compressive strength [35,36].In the context of the numerical modeling of NHC core machining, multiple layers were not included in the model, as their impact on this process was considered negligible.Indeed, the single-layer orthotropic approach is preferred for modeling due to its simplicity of implementation.The mathematical equations describing the stress-strain relationship of orthotropic material can be summarized as follows: where σ ij and ε ij characterize the stress tensor and strain tensor, respectively.C designates the stiffness matrix, with i, j = 1, 2, 3,. .., 6 [37].
The orthotropic mechanical properties of Nomex paper were determined from the study conducted by Kilchert and are presented concisely in Table 2 [38].To evaluate the damage to the Nomex material, a Tsai-Wu failure criterion is used according the following equation: [39]: where F i and F j represent the second-and fourth-order force tensors, respectively, defined as follows: where X and Y represent the limiting stresses of the material in the longitudinal and transverse directions, respectively, while S denotes the limiting shear stress, as shown in Table 2.
When the failure criterion is reached, this means that a progressive damage process occurs, leading to the deterioration of the previously mentioned stiffness matrix.
The Tsai-Wu criterion constitutes a tool for assessing damage within each stressed element of a structure when it meets the conditions set out in Equation ( 6).This method allows you to determine if an element has suffered damage by comparing the stresses applied to the resistances of the material in different directions.To resolve the direction challenge and calculate the decrease in element stiffness, the Tsai-Wu criterion was decomposed into several stress components, according to Equation (8).
More precisely, once the Tsai-Wu failure criterion is reached, the H k components will be calculated and compared to determine their maximum value.The highest H k component is selected as the predominant damage mode, leading to a decrease in the stiffness of the associated elastic modulus and a deterioration in the corresponding properties, as shown in Table 3 [40,41].In this analysis, the damage method used aims to decrease the stiffness of the individual elastic modules as soon as the breaking point is reached.For the Tsai-Wu criterion, the damage assessment is carried out based on several variables denoted as dk (with k ranging from 1 to 6), as illustrated by Equations ( 9) and (10).The formula describing the reduction in stiffness depends on the type of material and its mechanical behavior.For composite materials, the formula for stiffness reduction is calculated as follows [42]: with m: softening parameter; e: Euler number; ε f ij : deformation of rupture, which are determined as follows:

Cutting Force Calculation
During the experimental phase, the components of the cutting force (F y , the feed component along Y, and F x , the cutting component along X) are measured using the KISTLER-9256C2 dynamometer.This method makes it possible to calculate the average value in both directions, according to the following formulas: F x and F y represent the averages of the components of the cutting force in the X and Y directions, while t 1 and t 2 represent, respectively, the start and end of the cutting process.

Results and Discussion
The deformation of the NHC core during the milling process can have significant consequences on the cutting force and the machined surface quality.The main cause of cell deformation is the low rigidity of the Nomex paper that forms the NHC core.However, the NHC structure has high rigidity in the direction of the thickness T along the z axis [43,44].Therefore, the crushing component F z has no significant influence on the machining process.This section of research focuses on three key aspects.Firstly, the influence of cutting conditions on the F x and F y components was carefully analyzed.In parallel, an in-depth evaluation of the impact of feed component F y on the quality of the machined surface was carried out.Simultaneously, careful study focused on examining the influence of vibration amplitude on chip accumulation in front of the cutting tool.

Influence of the Feed Rate on the Cutting Force
In order to validate the numerical model of the milling of the NHC structure by RUM technology using the UCK cutting tool, numerical simulations were carried out according to the experimental protocol presented above.The aim of these simulations was to evaluate the influence of the feed rates such as 500 mm/min, 1000 mm/min, 1500 mm/min, 2000 mm/min, 2500 mm/min and 3000 mm/min on components F x and F y , in which the simulations were carried out both with and without the use of ultrasonic vibrations.The other cutting conditions were kept constant throughout the numerical modeling, such as the spindle speed at 5000 rpm, the vibration amplitude at 25 µm and the depth of cut at 2 mm.An in-depth comparison was carried out between the obtained results by the numerical model and those obtained by the experiment [34] (see Figure 7).mm/min, 2500 mm/min and 3000 mm/min on components Fx and Fy, in which the simulations were carried out both with and without the use of ultrasonic vibrations.The other cu ing conditions were kept constant throughout the numerical modeling, such as the spindle speed at 5000 rpm, the vibration amplitude at 25 µm and the depth of cut at 2 mm.An in-depth comparison was carried out between the obtained results by the numerical model and those obtained by the experiment [34] (see Figure 7).The obtained results show that components Fx and Fy increase as the feed rate of the cu ing tool increases, both for the simulation and for the experiment, whether or not the ultrasonic vibrations are used.The increase in components Fx and Fy is mainly due to the The obtained results show that components F x and F y increase as the feed rate of the cutting tool increases, both for the simulation and for the experiment, whether or not the ultrasonic vibrations are used.The increase in components F x and F y is mainly due to the increase in the amount of material removed per unit time.Furthermore, it is possible that the elastic deformation of the Nomex paper forming the NHC structure is responsible for this phenomenon by delaying the removal of the elements, leading to an accumulation of material in front of the tool and an overestimated resistance of the material, hence the increase in components F x and F y .When cutting the NHC structure using RUM technology, it is evident that the use of ultrasonic vibrations considerably reduces components F x and F y .In this context, the strong rotation and vibration of the cutting tool favorize the creation of cracks in the walls of the NHC structure, making it easier for the cutting tool to penetrate without resistance from the material forming the NHC structure.In contrast, when machining the NHC core without ultrasonic vibration, there is no vibration of the cutting tool, which is also reflected in the absence of crack propagation in the contact zone between the tool and the thin walls of the NHC structure, thus increasing components F x and F y .Additionally, when observing the process, it was noted that the feed force F y was consistently greater than cutting component F x .This results from the resistance of the material that builds up in front of the cutting tool during the milling process.On the other hand, the cutting force F x in the OX cutting direction is quite low due to the specific characteristic of Nomex paper, which is considered fragile with a low density.The values of the obtained components F x and F y by the numerical simulation present a significant correspondence with the experimental results, demonstrating the robustness of the developed numerical model.

Influence of the Spindle Speed on the Cutting Force
This paragraph presents a series of rotary ultrasonic machining (RUM) simulations of the NHC core, in which the study was carried out both with and without the use of ultrasonic vibrations.The main objective was to evaluate the influence of different spindle speeds, such as 500 rpm, 1000 rpm, 2000 rpm, 3000 rpm, 4000 rpm and 5000 rpm, on components F x and F y generated by the cutting tool.It essential to note that the feed rate used to perform the numerical simulations was 500 mm/min, the vibration amplitude was 25 µm and the depth of cut was 2 mm.The numerical simulations were carried out for a duration of 0.24 s, which is equivalent to a cutting width of 2 mm.The results calculated by the FE model are compared with the experimental results and are presented in Figure 8 [34].
ponents Fx and Fy generated by the cu ing tool.It essential to note that the feed rate used to perform the numerical simulations was 500 mm/min, the vibration amplitude was 25 µm and the depth of cut was 2 mm.The numerical simulations were carried out for a duration of 0.24 s, which is equivalent to a cu ing width of 2 mm.The results calculated by the FE model are compared with the experimental results and are presented in Figure 8 [34].According to the results presented in Figure 8, a decrease in components Fx and Fy is observed with increasing spindle speed, both for the simulation and for the experiment, According to the results presented in Figure 8, a decrease in components F x and F y is observed with increasing spindle speed, both for the simulation and for the experiment, whether or not ultrasonic vibrations are used.As the spindle speed increases, the kinetic energy of the cutting tool increases significantly.This results in an increase in the material removal rate per unit time, which systematically leads to a decrease in the F x and F y components.On the contrary, at low spindle speeds, the cutting tool exerts insufficient pressure on the walls of the hexagonal structure, which can lead to the bending or deformation of the walls.Consequently, the inappropriate cutting of Nomex paper results in uncut fibers, weakening the overall structure of the material and thus hindering the efficient evacuation of the formed chips.These chips can accumulate around the cutting zone, causing local overloading of the tool and the disruption of the cutting process, thus increasing the components of the cutting force in the X and Y directions.Furthermore, for all spindle speeds examined, feed component F y is higher than cutting component F x .This is due to the accumulation of chips in front of the cutting tool UCK, generating additional force in the direction of advance and thus causing an increase in the advance component F y .On the other hand, the fragility of Nomex paper can lead to less resistance of the cutting tool in the cutting direction X, resulting in relatively low values of component F x .Generally, the application of ultrasonic vibrations during the milling process results in a significant reduction in the F x and F y components compared to conventional cutting, thereby helping to prevent premature wear of the cutting tool.By comparing the presented results, a satisfactory correlation was observed between the components of the cutting force calculated by the numerical model and those resulting from the experiment, thus demonstrating the reliability of the developed numerical model.

Influence of the Vibration Amplitude on the Cutting Force
As part of this study, numerical simulations were carried out to examine the impact of the vibration amplitude of the cutting tool on the F x and F y components using RUM technology.This research was carried out by changing the vibration amplitudes of the cutting tool, namely to 5 µm, 10 µm, 15 µm, 25 µm and 27 µm.Meanwhile, the machining parameters, such as feed rate of 500 mm/min, spindle speed of 3000 rpm and cutting depth of 2 mm, were kept constant throughout the machining process.It should be noted that the numerical simulations were carried out for a period of 0.36 s, corresponding to a cutting width of 3 mm.The results obtained by the numerical model are compared with the experimental results and are presented in Figure 9.
cu ing tool, namely to 5 µm, 10 µm, 15 µm, 25 µm and 27 µm.Meanwhile, the machining parameters, such as feed rate of 500 mm/min, spindle speed of 3000 rpm and cu ing depth of 2 mm, were kept constant throughout the machining process.It should be noted that the numerical simulations were carried out for a period of 0.36 s, corresponding to a cutting width of 3 mm.The results obtained by the numerical model are compared with the experimental results and are presented in Figure 9.According to the results shown in Figure 9, there is a decrease in components F x and F y as a function of increasing vibration amplitude.The reduction in the F x component is explained by the creation of incisions resulting from the rotation and vibration of the cutting tool.This action promotes the propagation of cracks in the phenolic resin, which facilitates the penetration of the cutting tool into the walls of the NHC structure.For vibration amplitudes between 5 µm and 10 µm, a slight decrease in the F y component is noted.In fact, when the tool vibrates slightly, this reduces the impact on the part, thus limiting the propagation of cracks.In this case, the space between the cutting tool and the walls of the NHC structure is very small, resulting in closer contact between the two.As a result, friction during the cutting process increases, resulting in a higher cutting force.By increasing the values of the vibration amplitude, a clear decrease in power supply component F y is observed.In this case, the contact between the cutting tool and the walls of the NHC structure becomes less significant, which leads to a reduction in the force generated.In general, a correlation between the simulation and the experiment is clearly observed, thus demonstrating the reliability of the numerical model developed.

Analysis of Machined Surface Quality
The quality of the machined surface of NHC cores plays a crucial role in the construction of Nomex honeycomb sandwich composite structures.Therefore, optimizing the surface quality of this type of structure is essential to prevent delamination between the core and the skins that form the sandwich structure.According to experimental results, the most common machining defects in NHC core milling are burrs and tears to fibers.Therefore, optimizing the cutting tool edge is of paramount importance in the machining process.In order to validate the numerical model of the milling of the NHC structure, numerical simulations were carried out to highlight the influence of component F y on the quality of the machined surface.These simulations were carried out while strictly respecting the previously established cutting conditions for two ranges of the values of F y : F y ≤ 3N and F y > 3N.During the experimental phase, the machined surface quality defects were revealed using an optical microscope (model 55XA, SOIF, Shanghai, China) [34].However, the machining defects resulting from the numerical model were identified on the basis of visual examination with the naked eye.Figures 10 and 11 present the machined surfaces obtained by the numerical model and the experiment.
previously established cu ing conditions for two ranges of the values of Fy: Fy ≤ 3N and Fy > 3N.During the experimental phase, the machined surface quality defects were revealed using an optical microscope (model 55XA, SOIF, Shanghai, China) [34].However, the machining defects resulting from the numerical model were identified on the basis of visual examination with the naked eye.Figures 10 and 11 present the machined surfaces obtained by the numerical model and the experiment.The main observation is that machining defects are closely linked to the val component Fy, whether with UV or without.The results of the experiment hav that the main machining defects observed on the NHC core are burrs and tears fiber [34].On the other hand, the machining defects revealed by the numerical si are deformation and tears to the walls.It should be noted that the experimental not take into account the potential deformation of the walls during the machining Furthermore, burrs, which manifest themselves as excess material on the cell wa not be detected by the numerical model, as they were represented by S4R shell without thickness.In general, it was observed that the tears and deformations of were more pronounced for the high values of the advance force Fy > 3N in bo cases, whether with or without UV.In fact, the high values of Fy were obtained the low spindle speed of the cu ing tool, thus causing the elastic deformation o walls until they tore.For the low spindle speed, the cu ing tool does not apply pressure to the walls of the hexagonal structure, which can cause elastic deform lead to their folding.As a result, some fibers in Nomex paper may remain intact uncut fibers that compromise the overall strength of the structure and cause dete in the quality of the machined surface.On the other hand, the numerical and exp studies have clearly demonstrated the benefits of ultrasonic vibrations during th process.These vibrations reduce contact between the tool cu ing edge and the t of the NHC structure, which reduces friction and improves heat dissipation du ting.As a result, the cu ing tool can penetrate the walls of the NHC structure mo The main observation is that machining defects are closely linked to the value of feed component F y , whether with UV or without.The results of the experiment have shown that the main machining defects observed on the NHC core are burrs and tears to uncut fiber [34].On the other hand, the machining defects revealed by the numerical simulation are deformation and tears to the walls.It should be noted that the experimental tests did not take into account the potential deformation of the walls during the machining process.Furthermore, burrs, which manifest themselves as excess material on the cell walls, could not be detected by the numerical model, as they were represented by S4R shell elements without thickness.In general, it was observed that the tears and deformations of the walls were more pronounced for the high values of the advance force F y > 3N in both tested cases, whether with or without UV.In fact, the high values of F y were obtained thanks to the low spindle speed of the cutting tool, thus causing the elastic deformation of the thin walls until they tore.For the low spindle speed, the cutting tool does not apply sufficient pressure to the walls of the hexagonal structure, which can cause elastic deformation or lead to their folding.As a result, some fibers in Nomex paper may remain intact, creating uncut fibers that compromise the overall strength of the structure and cause deterioration in the quality of the machined surface.On the other hand, the numerical and experimental studies have clearly demonstrated the benefits of ultrasonic vibrations during the cutting process.These vibrations reduce contact between the tool cutting edge and the thin walls of the NHC structure, which reduces friction and improves heat dissipation during cutting.As a result, the cutting tool can penetrate the walls of the NHC structure more easily, thereby reducing machining defects and improving the quality of the machined surface.By analyzing the results of the numerical and comparing them with those obtained by the experiment, we found a strong correlation between the simulation and experiment, which reinforces the confidence of the numerical model developed.

Analysis of Chip Distribution in Front of the Cutting Tool
The accumulation of chips in front of the cutting tool represents a technical challenge for engineers and researchers in certain fields, particularly machining and manufacturing.Chip accumulation in front of the cutting tool can cause various problems.It disrupts the cutting movement, generating unwanted vibrations likely to deteriorate the quality of the surface and accelerate the wear of the cutting tool.As part of this study, numerical simulations were carried out in order to highlight the influence of the vibration amplitude on the accumulation of the chips in front of the cutting tool.For this purpose, three amplitudes were tested: 5 µm, 15 µm and 27 µm.The numerical simulations were performed using the same cutting conditions, including a spindle speed of 5000 rpm and a feed rate of 3000 mm/min.The duration of each simulation was 0.06 s, which corresponded to a cutting width of the 3 mm.The obtained results are shown in Figures 12-14    the accumulation of the chips in front of the cu ing tool.For this purpose, three amplitudes were tested: 5 µm, 15 µm and 27 µm.The numerical simulations were performed using the same cu ing conditions, including a spindle speed of 5000 rpm and a feed rate of 3000 mm/min.The duration of each simulation was 0.06 s, which corresponded to a cu ing width of the 3 mm.The obtained results are shown in Figures 12-14.The obtained results show that the amplitude of ultrasonic vibrations directly influences the accumulation of chips in front of the cu ing tool.In this context, it is obvious that high amplitudes significantly reduce the accumulation of chips in front of the cu ing tool.For low amplitudes, the contact between the edge of the tool and the walls of the structure becomes very significant.In this case, the cu ing tool does not easily penetrate the material, creating short-length cracks on the walls of the NHC structure.Thus, the poor cu ing of the walls of the NHC structure promotes their elastic deformation and leads to the accumulation of chips in front of the cu ing tool.By increasing the amplitude of ultrasonic vibrations beyond 5 µm, a notable reduction in chip accumulation in front of the cu ing tool is observed.In this case, the high amplitudes of vibrations and the high rotation speed reduce the contact between the edge of the cu ing tool and the walls of the NHC structure, which facilitates the penetration of the cu ing tool into the thin wall, consequently intensifying the propagation of cracks.Thus, the walls of the structure are cut The obtained results show that the amplitude of ultrasonic vibrations directly influences the accumulation of chips in front of the cutting tool.In this context, it is obvious that high amplitudes significantly reduce the accumulation of chips in front of the cutting tool.For low amplitudes, the contact between the edge of the tool and the walls of the structure becomes very significant.In this case, the cutting tool does not easily penetrate the material, creating short-length cracks on the walls of the NHC structure.Thus, the poor cutting of the walls of the NHC structure promotes their elastic deformation and leads to the accumulation of chips in front of the cutting tool.By increasing the amplitude of ultrasonic vibrations beyond 5 µm, a notable reduction in chip accumulation in front of the cutting tool is observed.In this case, the high amplitudes of vibrations and the high rotation speed reduce the contact between the edge of the cutting tool and the walls of the NHC structure, which facilitates the penetration of the cutting tool into the thin wall, consequently intensifying the propagation of cracks.Thus, the walls of the structure are cut optimally, which favors the formation of small chips and reduces their accumulation in front of the cutting tool.However, it should be noted that chip accumulation persists for all amplitudes examined.This can be attributed to the mechanical behavior of Nomex paper, which is considered elastic in this model.In this regard, the elastic deformation of Nomex paper delays the removal of elements, which causes an accumulation of material in front of the cutting tool forces and an overestimation of the material's strength.By analyzing the obtained results, we noticed that the conclusions drawn from the simulation provided a coherent explanation of everything we had discussed previously.As a result, the numerical model developed is represented as a numerical platform capable of handling all special cases.

Conclusions
In this work, a three-dimensional numerical model was developed to analyze the machining process of NHC structures using RUM technology both with and without the use of ultrasonic vibrations.This model has been rigorously validated by experimental tests comparing the F x and F y components, as well as the quality of the machined surface, under different cutting conditions.Furthermore, the numerical model offers the possibility of evaluating the influence of the amplitude of ultrasonic vibrations on the accumulation of chips in front of the cutting tool.The results of the study led to the following conclusions: • The influence of the feed rate on components F x and F y was carefully examined; it was observed that these components increase with the increase in the feed rate, both for the simulation and just for the experiment.Furthermore, the results indicate that the use of ultrasonic vibrations helps minimize the negative effects of components F x and F y for both directions.In the end, concrete agreement is observed between the obtained results by the numerical model and the results determined by the experiment.• The analysis of components F x and F y as a function of the spindle speed reveals similarities between the orthotropic elastic behavior and that of composite materials.The significant reduction in components F x and F y with high spindle speeds is observed for the two cases tested, whether with or without UV.Interestingly, the simulation based on the orthotropic behavior with the Tsai-Wu failure criterion follows a similar trend as the experimental results.

•
By analyzing the influence of feed component F y on the quality of the machined surface, it becomes clear that this relationship is firmly established.The combination between ultrasonic vibrations and the reduction in feed component F y demonstrates a significant improvement in the quality of the machined surface.A significant correlation is noted in the obtained results by the numerical model and those obtained by the experiment, which confirms the credibility of the proposed numerical model.• The developed numerical model makes it possible to study the influence of vibration amplitude on components F x and F y .The obtained results clearly indicate that higher vibration amplitudes result in a significant reduction in these components, thus helping to avoid premature wear of the cutting tool.In the end, concrete agreement is noted between the obtained results by the numerical model and those determined by the experiment.• The numerical model shows that feed component F y prevails over F x , due to the accumulation of material in front of the cutting tool and the low density of the Nomex paper constituting the NHC structure.• The numerical model makes it possible to analyze the influence of vibration amplitude on chip accumulation in front of the cutting tool.The higher amplitudes appear to help reduce unwanted chip buildup, helping to improve machining efficiency, optimize machined surface quality and prevent premature chip wear of the cutting tool.
• The numerical model allows the analysis of complex machining simulations, which are difficult to perform experimentally.It thus offers researchers and engineers the possibility of improving cutting performance while avoiding costly experimental tudies.

Figure 1 .
Figure 1.The experimental configuration used for cu ing the NHC core by RUM technology.

Figure 1 .
Figure 1.The experimental configuration used for cutting the NHC core by RUM technology.

FrequencyFigure 2 .
Figure 2. (a) Geometric dimensions of the NHC core; (b) Geometric dimensions of the hexagonal cell.

Figure 2 .
Figure 2. (a) Geometric dimensions of the NHC core; (b) Geometric dimensions of the hexagonal cell.

Machines 2024 , 20 Figure 3 .
Figure 3. Presentation UCK cu ing tool used to cut the NHC core.Figure 3. Presentation of the UCK cutting tool used to cut the NHC core: (a) UCK designed by simulation; (b) UCK used in the experiment.

Figure 3 .
Figure 3. Presentation UCK cu ing tool used to cut the NHC core.Figure 3. Presentation of the UCK cutting tool used to cut the NHC core: (a) UCK designed by simulation; (b) UCK used in the experiment.

Figure 3 .
Figure 3. Presentation UCK cu ing tool used to cut the NHC core.

Figure 4 .
Figure 4. Dimensions of the UCK cu ing tool.

Figure 4 .
Figure 4. Dimensions of the UCK cutting tool.

Figure 5 .
Figure 5. (a) Mesh used in the FE model; (b) cu ing conditions defined at reference point RP: n is the spindle speed and Vf is the feed rate.

Figure 5 .
Figure 5. (a) Mesh used in the FE model; (b) cutting conditions defined at reference point RP: n is the spindle speed and V f is the feed rate.

Figure 6 .
Figure 6.Boundary conditions used in the numerical simulation: (a) used in the simulation; (b) used in the experience.

Figure 6 .
Figure 6.Boundary conditions used in the numerical simulation: (a) used in the simulation; (b) used in the experience.

Figure 7 .
Figure 7. Evolution of Fx and Fy according to the feed rates.

Figure 7 .
Figure 7. Evolution of F x and F y according to the feed rates.

Figure 8 .
Figure 8. Evolution of Fx and Fy according to the spindle speeds.

Figure 8 .
Figure 8. Evolution of F x and F y according to the spindle speeds.

Figure 9 .
Figure 9. Evolution of Fx and Fy according to the vibration amplitudes.Figure 9. Evolution of F x and F y according to the vibration amplitudes.

Figure 9 .
Figure 9. Evolution of Fx and Fy according to the vibration amplitudes.Figure 9. Evolution of F x and F y according to the vibration amplitudes.

Figure 10 .
Figure 10.Surface quality resulting from RUM technology with UV: (a) Surface quality resulting from the numerical model; (b) surface quality resulting from the experiment.

Figure 10 .Figure 11 .
Figure 10.Surface quality resulting from RUM technology with UV: (a) Surface quality resulting from the numerical model; (b) surface quality resulting from the experiment.Machines 2024, 12, x FOR PEER REVIEW

Figure 11 .
Figure 11.Surface quality resulting from RUM technology without UV: (a) Surface quality resulting from the numerical model; (b) surface quality resulting from the experiment. .

Machines 2024 ,
12,  x FOR PEER REVIEW 15 of 20 the accumulation of the chips in front of the cu ing tool.For this purpose, three amplitudes were tested: 5 µm, 15 µm and 27 µm.The numerical simulations were performed using the same cu ing conditions, including a spindle speed of 5000 rpm and a feed rate of 3000 mm/min.The duration of each simulation was 0.06 s, which corresponded to a cu ing width of the 3 mm.The obtained results are shown in

Figure 12 .
Figure 12.Accumulation of the chips in the front of the UCK cu ing tool, for A= 5µm.

Figure 12 .
Figure 12.Accumulation of the chips in the front of the UCK cutting tool, for A = 5 µm.

Figure 12 .
Figure 12.Accumulation of the chips in the front of the UCK cu ing tool, for A= 5µm.

Figure 13 .
Figure 13.Accumulation of the chips in the front of the UCK cu ing tool, for A= 15µm.

Figure 13 .
Figure 13.Accumulation of the chips in the front of the UCK cutting tool, for A = 15 µm.

Figure 14 .
Figure 14.Accumulation of the chips in the front of the UCK cu ing tool, for A= 27µm.

Figure 14 .
Figure 14.Accumulation of the chips in the front of the UCK cutting tool, for A = 27 µm.

Table 1 .
Performance conditions of the THU Ultrasonic 850.

Table 2 .
The mechanical characteristics of Nomex paper.

Table 3 .
Degradation of mechanical properties and corresponding variables according to the Tsai-Wu criterion.