High-Speed Cutting of Synthetic Trabecular Bone—A Combined Experimental–Computational Investigation

: Orthopaedic surgical cutting instruments are required to generate sufﬁcient forces to penetrate bone tissue while minimising the risk of thermal and mechanical damage to the surrounding environment. This study presents a combined experimental–computational approach to determine relationships between key cutting parameters and overall cutting performance of a polyurethane-based synthetic trabecular bone analogue under orthogonal cutting conditions. An experimental model of orthogonal cutting was developed, whereby an adaptable cutting tool ﬁxture driven by a servo-hydraulic uniaxial test machine was used to carry out cutting tests on Sawbone ® trabecular bone analogues. A computational model of the orthogonal cutting process was developed using Abaqus/Explicit, whereby an Isotropic Hardening Crushable Foam elastic-plastic model was used to capture the complex post-yield behaviour of the synthetic trabecular bone. It was found that lower tool rake angles resulted in the formation of larger discontinuous chips and higher cutting forces, while higher rake angles tended to lead to more continuous chip formation and lower cutting forces. The computational modelling framework provided captured features of both chip formation and axial cutting forces over a wide range of cutting parameters when compared with experimental observations. This experimentally based computational modelling framework for orthogonal cutting of trabecular bone analogues has the potential to be applied to more complex three-dimensional cutting processes in the future.


Introduction
Bone is a naturally occurring composite material whose constituent phases are hierarchically organised to provide a highly optimised structure that exhibits high stiffness and excellent resistance to fracture while also being lightweight for efficient movement. The complex structural arrangement of bone gives rise to anisotropic, rate-dependent failure behaviour, which varies significantly depending on tissue composition and architecture [1][2][3]. This presents significant challenges in the development of surgical cutting instruments, which are required to generate sufficient forces to penetrate bone tissue while minimising the risk of thermo-mechanical damage to the surrounding tissue. The biomechanics of cutting processes is further complicated by the range of complex interactions taking place at the tool-bone interface, particularly when three-dimensional cutting processes such as sawing (James et al., 2014), drilling [4] and burring [5] are considered. Currently, instrument designers rely heavily on empirical-based strategies to understand tool-bone interactions, with significant amounts of prototyping and validation experiments required throughout the design process. Both experimental and computational models of twodimensional orthogonal cutting are widely used as a platform to investigate the influence of key cutting parameters, such as tool geometry, cutting depth and/or feed rate, on overall cutting performance [6][7][8][9][10]. However, the vast majority of these studies have focussed on cutting process in cortical bone, and to date, there is a distinct lack of studies examining the biomechanics of cutting in trabecular bone [11].

Dynamic Uniaxial Compression Experiments
Uniaxial compression tests were performed using the hydraulic Dartec Series HC10 uniaxial test machine (Dartec, Cranfield, UK) to evaluate the rate-dependent properties of the PU foam. Cubic specimens of Sawbone ® grade 20 (ρ = 0.32 g/cm 3 ) solid-rigid closedcell PU foam (Pacific Research Labs, Malmö, Sweden) measuring 10 mm × 10 mm × 10 mm were tested at three different uniaxial strain rates, . ε = 0.1 s −1 , 1 s −1 and 10 s −1 , to a nominal strain of ε = 0.6 at room temperature (n = 3). Force-displacement data were recorded by the machine load cell and crosshead, from which nominal stress-nominal strain curves were constructed based on undeformed specimen dimensions. The apparent compressive Young's modulus (E) was computed as the slope of the initial linear-elastic region. The yield stress (σ y ) was determined through the 0.2% offset method. The results of these tests formed the basis of constitutive model parameters used in the computational model.

High-Speed Orthogonal Cutting Experiments
Similar to Malak and Anderson [8], an experimental model of the high-speed orthogonal cutting process was developed to evaluate cutting forces generated during high-speed orthogonal cutting of the PU foam analogue. This consisted of a custom-built cutting tool and hydraulic Dartec Series HC10 uniaxial testing machine (see Figure 1). The cutting tool was attached to the upper crosshead of the Dartec machine, with the PU foam samples fixed beneath using custom-built clamps. Using a removable pin, the cutting tool rake angle (α) was varied to test α = 0 • , 20 • , 45 • and 60 • , with corresponding clearance angles of 55 • , 35 • , 10 • and 5 • , respectively. The PU foam was machined to rectangular specimens of 30 mm length, 40 mm height and 7.5 mm width. A total of 24 high-speed cutting tests were performed on the material, using a constant cutting velocity of v = 200 mm/s, with two separate cutting depths of h = 1 mm and h = 2 mm tested. This meant that three samples were tested for each configuration considered (n = 3). Force-displacement data were recorded, with results from axial cutting force (parallel to the direction of cutting) normalised and presented as the cutting force per unit (mm) thickness of material in the z-direction ( Figure 1). A Sony DSC-RX100 III (Sony Group, Tokyo, Japan) high-speed camera was used to visualise deformation and chip formation during the cutting experiments, with images captured at 250 frames per second.

Orthogonal Cutting Geometry
A finite element model of the orthogonal cutting process was developed using the Abaqus/Explicit 2017 finite element code, with sample dimensions and boundary conditions replicating the experimental setup, as shown in Figure 2. Here, three-dimensional elements were used to represent the workpiece and cutting tool, with both components comprising a single element in the thickness direction, with boundary conditions applied to planar faces to represented plane strain conditions during simulations. A reference point constrained to the cutting tool was used to apply the cutting velocities, with reaction forces and displacement extracted to generate cutting force data. A surface-to-surface contact (Explicit) formulation was used, with 'Hard' contact used to represent normal behaviour and a penalty friction algorithm used to describe tangential behaviour, in which a friction coefficient of µ = 0.25 assumed [29]. The PU workpiece was represented as a homogenous continuum, and the cutting region was discretised using~146,000 8-node linear brick elements (C3D8R) with distortion control, enhanced hourglass control and element deletion activated. Wedge elements (C3D6) were used outside the tool-material contact region. The cutting tool was modelled as a discrete rigid part and discretised using 130 4-node 3D bilinear rigid quadrilateral elements (R3D4). A mesh dependency study was conducted to ensure that the results obtained from the finite element analysis were independent of the mesh quality. Similar to the orthogonal cutting experiment, the cutting force was normalised and is presented as the cutting force per unit (mm) thickness of the material.

Orthogonal Cutting Geometry
A finite element model of the orthogonal cutting process was developed using the Abaqus/Explicit 2017 finite element code, with sample dimensions and boundary conditions replicating the experimental setup, as shown in Figure 2. Here, three-dimensional elements were used to represent the workpiece and cutting tool, with both components comprising a single element in the thickness direction, with boundary conditions applied to planar faces to represented plane strain conditions during simulations. A reference point constrained to the cutting tool was used to apply the cutting velocities, with reaction forces and displacement extracted to generate cutting force data. A surface-to-surface contact (Explicit) formulation was used, with 'Hard' contact used to represent normal behaviour and a penalty friction algorithm used to describe tangential behaviour, in which a friction coefficient of μ = 0.25 assumed [29]. The PU workpiece was represented as a homogenous continuum, and the cutting region was discretised using ~146,000 8-node linear brick elements (C3D8R) with distortion control, enhanced hourglass control and element deletion activated. Wedge elements (C3D6) were used outside the tool-material contact region. The cutting tool was modelled as a discrete rigid part and discretised using ~130 4-node 3D bilinear rigid quadrilateral elements (R3D4). A mesh dependency study was conducted to ensure that the results obtained from the finite element analysis were independent of the mesh quality. Similar to the orthogonal cutting experiment, the cutting force was normalised and is presented as the cutting force per unit (mm) thickness of the material.

Material Constitutive Model Implementation
To capture the behaviour of the workpiece during the cutting process, two plasticity formulations were implemented: (i) Johnson-Cook plasticity, which has been widely used but does not account for pressure-dependent yielding and (ii) Isotropic Crushable Foam (CFI) plasticity, which does include a pressure-dependent yield criterion. A damage evolution law was also implemented in conjunction with both plasticity models to account for post-yield softening/failure behaviour.
(i) Johnson-Cook plasticity: The isothermal strain-rate dependent Johnson-Cook plasticity model has the following form and is essentially a strain-rate dependent version of the Von Mises yield function, given by, where is the material flow stress, is the effective plastic strain, and = ,

Material Constitutive Model Implementation
To capture the behaviour of the workpiece during the cutting process, two plasticity formulations were implemented: (i) Johnson-Cook plasticity, which has been widely used but does not account for pressure-dependent yielding and (ii) Isotropic Crushable Foam (CFI) plasticity, which does include a pressure-dependent yield criterion. A damage evolution law was also implemented in conjunction with both plasticity models to account for post-yield softening/failure behaviour.
(i) Johnson-Cook plasticity: The isothermal strain-rate dependent Johnson-Cook plasticity model has the following form and is essentially a strain-rate dependent version of the Von Mises yield function, given by, where σ is the material flow stress, ε p e f f is the effective plastic strain, and ε 0 is the reference strain rate used to determine the model constants A, B & n. The first bracketed term accounts for strain-hardening effects, while the second bracketed term accounts for strain-rate dependency. Model constants for the PU20 were determined from the results from uniaxial compression tests, described in Section 2.1, and are provided in Table 1.  Deshpande and Fleck (2000) and has been widely used to represent the post-yield behaviour of cellular materials. The CFI model phenomenologically represents the buckling of cell walls under compression and has a pressure-dependent yield criterion, which is given by and and The flow potential for the CFI model is given by In the above equations, α is the shape of the yield ellipse in the Von Mises-Pressure stress plane, B is the size of the yield ellipse, p c is the yield strength in hydrostatic compression, q is the Von Mises yield stress, K is the compression yield stress ratio, σ 0 c is the initial yield strength in uniaxial compression, p 0 c is the initial yield strength in hydrostatic compression, β is the ellipse for the potential flow and v p is the plastic Poisson's ratio. The yield ellipse progression is governed by the equivalent plastic strain [30]. Values for the compression yield stress ratio and the plastic Poisson's ratio and are summarised in Table 2, which were obtained from the results of Kelly and McGarry [27], who also tested PU20. The isotropic hardening data were obtained on experimental uniaxial compression data and are prescribed in the model as a piecewise linear curve. The strain-rate dependence of the yield surface is also input in tabular form in the Abaqus software as part of the material property definition, which defines a yield stress ratio in terms of equivalent plastic strain rate, which is determined from the uniaxial compressive behaviour of the PU20 (described in Section 2.1). Damage evolution law: Both the Johnson-Cook and CFI plasticity models were coupled to a damage evolution law (often termed the Johnson-Cook damage model) to describe softening behaviour of the material post-yield. This is described by three failure parameters (d 1 , d 2 and d 3 ), represented by the following equation where ε f is the equivalent strain at failure and η is the stress triaxiality parameter, which is given as the ratio of hydrostatic stress to Von Mises equivalent stress. This damage evolution law was calibrated through an iterative process, whereby the damage parameters (d 1 -d 3 ) and the effective plastic displacement at the failure parameter (u pl f ) were adjusted until the predicted chip formation provided good agreement with those observed in the experiments ( Table 3). The damage parameters (d 1 -d 3 ) and the effective plastic displacement at failure parameter (u pl f ) used for these simulations are provided in Table 3. The damage evolution law defines the element removal process. Damage is initiated once the threshold for equivalent strain at failure is exceeded. Beyond this point, a reduction in stiffness and strength occurs until the defined equivalent plastic strain at failure is reached; at this point, the element is deleted from the mesh. Table 3. Calibrated damage evolution law parameters for PU20.

High-Speed Orthogonal Cutting Simulations
Cutting simulations were performed using a dynamic, explicit step with automatic incrementation in the Abaqus Explicit finite element code. For all simulations, the cutting tool was prescribed a velocity of v = 200 mm/s. A variation of cutting parameters was carried out using this modelling framework, with rake angles of α = 0 • , 20 • , 45 • or 60 • and cutting depths of h = 1 mm and h = 2 mm simulated.

Categorisation of Chip Formation
This study categorises chip formation observed in cutting experiments according to the categories established by Malak and Anderson [8] in their experimental study on the orthogonal cutting of PU foams. Type 1, 2 and 3 chip formation categories are described in Table 4 and illustrated schematically in Figure 3, along with typical resulting cutting forcetool displacement curves that correspond to each type of chip formation. A new type of chip formation was observed in this current study, which has not been previously reported, and this has been defined as "Type 4" chip formation and is also described in Table 4. The images that we present in the results are those that were considered representative of the steady state cutting process observed. In some cases, key events such as a large crack occurred during the cutting process, and in these cases, such events are presented in the images in the results. Table 4. Categorisation of chip formation based on study by Malak and Anderson (2005).

Chip Type Chip Type
Surface fragmentation (Type 1) An accumulation of broken cell walls formed as the tool scrapes along the surface at low cut depths relative to cell size (not observed in this study, so not shown in Figure 3).

Continuous chip formation (Type 2)
A continuous plastic deformation of the material without fracture ahead of the cutting edge and results in the smooth flow of the chip up the rake face.

Discontinuous chip formation (Type 3)
Fragmented chips forming and may be further subdivided into (i) confined crack propagation and lifting of the chip (Type 3a), (ii) continuous crack propagation (Type 3b) and (iii) shearing (Type 3c).

Obliteration (Type 4)
Material ahead of the tool is obliterated, resulting in fine dust. Type 4 chips form a compression failure in the material, ahead of a relatively low rake angle, approaching 0 • . This forms an acceptable chip type in the form of fine dust and is a result of the brittle failure of the cell walls within the cellular material.
R PEER REVIEW 8 Figure 3. Schematic drawings of the types of chip formation are presented alongside the typical cutting force-tool displacement curves that correspond to each type of chip formation.

Uniaxial Dynamic Compression Tests
Experimentally measured nominal stress-nominal strain curves for the dynamic uniaxial compression testing of PU20 are shown in Figure 4, with associated Young's modulus and yield stress (mean ± SD (range)) shown in Table 5

Uniaxial Dynamic Compression Tests
Experimentally measured nominal stress-nominal strain curves for the dynamic uniaxial compression testing of PU20 are shown in Figure 4, with associated Young's modulus and yield stress (mean ± SD (range)) shown in Table 5. The measured values for Young's modulus are E = 155 ± 2.6 MPa and E = 197.9 ± 28.6 MPa at strain rates of . ε = 0.1 s −1 and . ε = 10 s −1 , respectively. The yield stress for PU20 was found to be σ y = 5.7 ± 0.1 MPa are σ y = 6.9 ± 0.3 MPa at strain rates of . ε = 0.1 s −1 and . ε = 10 s −1 , respectively. lus and yield stress (mean ± SD (range)) shown in Table 5. The measured values for Young's modulus are = 155 2.6 MPa and = 197.9 28.6 MPa at strain rates of = 0.1 s −1 and = 10 s −1 , respectively. The yield stress for PU20 was found to be = 5.7 0.1 MPa are = 6.9 0.3 MPa at strain rates of = 0.1 s −1 and = 10 s −1 , respectively.   Figure 4. Nominal stress-nominal strain curves for PU20 obtained from dynamic uniaxial compression testing at high strain rates.

Chip Formation
Following the analysis of the results from both the Johnson-Cook and CFI plasticity models, it was found that the CFI model was more accurate in predicting the experimental cutting forces (see Section 3.2.2) and better captured the chip formation. Therefore, Figure 5 compares experimental and computational (CFI model only) chip formation at two time points during the cutting process for all cutting parameters considered (i.e., rake angle and cutting depth). The computational model accurately predicts the formation and separation of chips when compared with the experimental observations. In general, it was found that high rake angles and low cutting depths tended towards more continuous (Type 2) chip formation while decreasing the rake angle and/or increasing the cutting depth led to the formation of more discontinuous chip types (Type 4, 3a, 3b, 3c). Figure 5a shows that for a α = 0 • and h = 2 mm, steady-state chip type was observed, with obliteration of the material (Type 4) and some discontinuous fragments occurring. Interestingly, a large crack developed in the sample towards the end of the test as the tool exited the workpiece, a feature that was also predicted in the computational model following its prediction of steady-state chip formation. Figure 5b shows that for a α = 0 • and h = 1 mm, material obliteration (Type 4) and discontinuous fragmentation occurred experimentally, with the computational model predicting both of those features. Figure 5c illustrates that for a α = 20 • and h = 2 mm, one discontinuous Type 3b chip formed, followed by a large crack that split the sample ahead of the tooltip. The model captured the features of discontinuous Type 3b chipping and the edge effect quite well. However, discontinuous Type 3b chipping also occurred computationally but not experimentally. Appl Figure 5d presents the chip formation for α = 20 • and h = 1 mm. The initial steadystate chip type was discontinuous, suggesting that these cutting parameters may be within a transition zone between the two chip types (Type 3a/3b). Similar to the previous case, as the tool exited the workpiece, a large crack occurred ahead of the tooltip. For this case, the model incorrectly predicted discontinuous Type 3c chipping but predicted similar cracking of the material due to the edge effect. Of the eight simulations, this was the only case in which the model failed to accurately predict the steady-state chip type. Figure 5e presents the chip formation for α = 45 • and h = 2 mm, in which discontinuous chip formation with confined crack propagation and lifting of the chips (Type 3a) occurred experimentally and was also captured by the computational model. Figure 5f presents the results for α = 45 • and h = 1 mm. Similar to Figure 5e, discontinuous chip formation with confined crack propagation occurred (Type 3a). However, in this case, the size of the chips formed was smaller due to the lower cutting depth. Again, these features were captured by the computational model. Figure 5g,h present the chip formation for α = 60 • at h = 2 mm and h = 1 mm, respectively. In both cases, discontinuous chip formation with confined crack propagation and lifting of the chips was observed (e.g., Type 3a). Again, for lower cutting depths, smaller chips were produced, with all these features being captured by the computational model. A summary of all chip types observed experimentally and those predicted by the computational model is provided in Table 6.

Cutting Force
The cutting force-tool displacement data measured from experimental tests (mean ± SD, n = 3) and predicted by computational models is presented in Figure 6. The largest experimental cutting forces were measured at α = 0 • and h = 2 mm, while the smallest were measured at α = 60 • and h = 1 mm. Figure 6a shows cutting force-displacement data for α = 0 • and h = 2 mm. The experimental curve shows an initial peak in cutting force followed by a large drop-off to zero, which is a result of the propagation of a large crack within the sample that was visible in chip formation analysis in Figure 5a. The computational CFI model predicts a similar cutting force to the experimental, with slight oscillations in the curve, with local maxima and minima caused by chip fragmentation. The crack that developed computationally (Figure 5a) did not cause complete separation of the sample, thus explaining why a large drop-off in the cutting force did not occur. While similar features are predicted by the J-C model, it over-predicts the cutting force compared with the measured experimental data for this case. Figure 6b shows the cutting force-displacement data for α = 0 • and h = 1 mm. The experimental data resembles a typical material obliteration force curve as characterised in Section 2.3. Figure 6c presents cutting force-displacement data for α = 20 • and h = 2 mm. The initial spike in cutting force is a result of the formation and separation of a large chip (Type 3b). The second spike in cutting force corresponds to a large crack propagating within the sample due to an edge effect. The CFI model predicted a combination of discontinuous chips (e.g., Type 3b & 3c; see Figure 5c), thus explaining the combination of cutting force-displacement curves (characterised in Section 2.3) displayed in Figure 6c. The J-C cutting force-displacement prediction has a similar profile to the CFI model but at a greater magnitude. Figure 6d shows cutting force-displacement results for α = 20 • and h = 1 mm. The experimental curve corresponds to discontinuous chip formation and is similar to that presented in Figure 6c. As previously mentioned (comments on Figure 5d), both computational models failed to accurately predict the experimental chip formation at these cutting parameters. This reflects the incorrect prediction of cutting force-displacement presented in Figure 6d. Appl

Wavelength of Discontinuous Chip Formation and Specific Cutting Energies
Discontinuous chip formation was observed during many of the cutting parameters tested, with cutting force-displacement curves showing fluctuations of distinct wavelengths. The size of these wavelengths is directly related to the size of the chips formed. These fluctuations were quantified based on the experimentally measured cutting force data, and Figure 7 shows the relationship between the wavelength of discontinuous chip formation and rake angle and cutting depth. Here, lower rake angles and larger cutting depths have higher measured wavelengths (larger chips were also observed in this case). No measurable wavelength was observed for a rake angle of 0°.  Figure 6e shows the cutting force-displacement data for α = 45 • and h = 2 mm. This type of curve corresponds to discontinuous chip formation with confined crack propagation, as characterised in Section 2.3. Discontinuous chip formation was found to fluctuate with a distinct wavelength, which was visually measured from the cutting force-displacement data. The wavelength is the distance between cutting force local maxima and corresponds to the length of the discontinuous chips formed. The cutting force local maxima and minima and curve wavelength are caused by confined crack propagation and lifting of the chips and is accurately predicted by the CFI model in this case (see Figure 5e). Figure 6f shows the cutting force-displacement data for α = 45 • and h = 1 mm. Similar to Figure 6e, this curve corresponds to discontinuous chip formation. The cutting force local maxima and curve wavelength are lower than in Figure 6e. The CFI model provides a good prediction of these results, with the J-C model overestimating the cutting forces in both Figure 6e,f. Figure 6g presents cutting force-displacement data for α = 60 • and h = 2 mm. The experimental data corresponds to discontinuous chip formation and is similar to Figure 6e,f. However, this case presents non-zero cutting force local minima. The CFI model slightly underestimates the cutting forces local minima but still predicts cutting force local maxima and curve wavelength. Figure 6h exhibits the cutting force-displacement data for α = 60 • and h = 1 mm. The experimental data could be misinterpreted as a corresponding curve to continuous chip formation. However, it, in fact, corresponds to discontinuous chip formation with a relatively small wavelength. In this case, the CFI model predicts lower cutting forces than those measured experimentally and predicted by the J-C model.

Wavelength of Discontinuous Chip Formation and Specific Cutting Energies
Discontinuous chip formation was observed during many of the cutting parameters tested, with cutting force-displacement curves showing fluctuations of distinct wavelengths. The size of these wavelengths is directly related to the size of the chips formed. These fluctuations were quantified based on the experimentally measured cutting force data, and Figure 7 shows the relationship between the wavelength of discontinuous chip formation and rake angle and cutting depth. Here, lower rake angles and larger cutting depths have higher measured wavelengths (larger chips were also observed in this case). No measurable wavelength was observed for a rake angle of 0 • . The specific cutting energy u, which is the energy consumed per unit volume of material removed (Shaw, 2005), was calculated using = ℎ (9) where is the cutting force, is the width of the sample and ℎ represents the cutting depth. Figure 7b presents the relationship between specific cutting energy, rake angle and cutting depth. Increasing rake angle and increasing cutting depth correspond to a decrease in specific cutting energy. The CFI model predictions for specific cutting energies are superimposed in Figure 7b, allowing comparison with the experimental results. The model provided a reasonable prediction of the experimental outcomes.  The specific cutting energy u, which is the energy consumed per unit volume of material removed (Shaw, 2005), was calculated using

Tool-Bone Interactions
where F is the cutting force, w is the width of the sample and h represents the cutting depth. Figure 7b presents the relationship between specific cutting energy, rake angle and cutting depth. Increasing rake angle and increasing cutting depth correspond to a decrease in specific cutting energy. The CFI model predictions for specific cutting energies are superimposed in Figure 7b, allowing comparison with the experimental results. The model provided a reasonable prediction of the experimental outcomes. Figure 8 shows contour plots of stress triaxiality (η), which is defined as the ratio of hydrostatic stress (p) to Von Mises stress (q), where η = p q . Here, the contour limits are chosen as η = −0.33 and η = 0.33, which represent values for unconfined compression and tension, respectively. For these simulations, large regions of the workpiece tended to have values of stress triaxiality that were outside these limits (−0.33 < η < 0.33), indicating that the hydrostatic pressure component ahead of the tooltip is substantial.

Tool-Bone Interactions
where is the cutting force, is the width of the sample and ℎ represents the cutting depth. Figure 7b presents the relationship between specific cutting energy, rake angle and cutting depth. Increasing rake angle and increasing cutting depth correspond to a decrease in specific cutting energy. The CFI model predictions for specific cutting energies are superimposed in Figure 7b, allowing comparison with the experimental results. The model provided a reasonable prediction of the experimental outcomes. Figure 8 shows contour plots of stress triaxiality ( ), which is defined as the ratio of hydrostatic stress (p) to Von Mises stress (q), where = . Here, the contour limits are chosen as = −0.33 and = 0.33, which represent values for unconfined compression and tension, respectively. For these simulations, large regions of the workpiece tended to have values of stress triaxiality that were outside these limits (−0.33 < < 0.33), indicating that the hydrostatic pressure component ahead of the tooltip is substantial.

Discussion
A combined experimental-computational approach was used to determine important relationships between key cutting parameters and overall cutting performance of synthetic trabecular bone analogue under high-speed orthogonal cutting conditions. The study presents a computational modelling framework for high-speed cutting of trabecular bone analogues, whose model parameters were determined experimentally. The modelling framework generally captured features of chip formation and axial cutting forces over a wide range of cutting parameters when compared with experimental measurements and observations and has the potential to be used to examine more complex three-dimensional cutting processes in the future. Through the development of a novel experimental cutting rig, this study represents the first investigation of the high-speed (v = 200 mm s −1 ) orthogonal cutting behaviour of PU foam. It was found that low rake angles (α = 0 • ) resulted in obliteration (Type 4) of the workpiece ahead of the tooltip and the highest cutting forces. Intermediate rake angles (20 • < α < 45 • ) typically produced discontinuous chip formation, either through confined crack propagation and lifting of the chip (Type 3a), which was usually associated with smaller cutting depths or continuous crack propagation ahead of the cutting tool (Type 3b) at larger cutting depths. Higher rake angles (α = 60 • ) resulted in the lowest cutting forces and tended towards more continuous chip formation (Type 2), providing a smoother finish on the cut surface. It is interesting that features of chip formation observed at each rake angle and cutting depth here are largely similar to the orthogonal cutting experiments carried out by Malak and Anderson [8], who tested a range of PU foams, however at much lower cutting speeds (v = 0.187 mm s −1 ), suggesting that the cutting rate does not have a drastic effect on visual chip formation. Our study also identified a new mode of chip formation at a rake angle of α = 0 • , defined as Type 4 chip formation, which describes compression failure ahead of the tooltip, resulting in the formation of fine dust and some material fragmentation. This observed behaviour could be relevant for the design of tools for burring operations, which tend to have very low (and sometimes negative) rake angles [16]. Our study also presents the efficiency of the cutting processes using the specific cutting energy (u), a parameter defined as the energy consumed per unit volume of material removed [31]. With the emergence of minimally invasive surgical techniques [32], the requirement for more compact tools and surgical equipment with high cutting efficiency is becoming more important. The results from this study show that higher rake angles (α = 20-60 • ) generate lower specific cutting energies during the cutting process, a feature that could be exploited to optimise the design and efficiency of smaller surgical cutting tools.
This study also developed a computational modelling framework for high-speed orthogonal cutting that generally captured the features of chip formation and cutting forces over a range of cutting parameters. This framework implements the strain rate-and pressure-dependent crushable foam yield criteria within the cutting model and directly compares its performance in capturing the experimental cutting process against the more widely used Johnson-Cook model (originally developed for metals), with material constants for both plasticity formulations determined experimentally. It is interesting that the CFI model provided a more accurate prediction of cutting forces, while the Johnson-Cook plasticity model consistently overestimated the cutting forces (see Figure 6) in almost all cases. Previous studies have demonstrated that the yield behaviour of PU foams is pressuredependent [27] or, more specifically, depends on stress triaxiality (η), which is defined as the ratio of hydrostatic stress (p) to Von Mises stress (q). However, the most widely used plasticity formulation in cutting simulations is the Johnson-Cook model, which is a strain rate-and temperature-dependent version of Von Mises plasticity. This implies that yielding in the Johnson-Cook model is independent of hydrostatic stress/triaxiality. This model should work well where deformation is relatively unconfined or where stress triaxiality values are between η = −0.33 and η = 0.33, respectively. For more confined stress states, or if the stress triaxiality value moves outside the limits (−0.33 < η < 0.33), it would be expected that a Johnson-Cook plasticity model may fail to capture yielding behaviour. Figure 8 shows contour plots of stress triaxiality (η = p q ) during the cutting process, with high (absolute) values of stress triaxiality in the workpiece ahead of the cutting tool, which leads to a general over-prediction of cutting forces by the J-C model within our computational framework. The implications of this finding are that the accuracy of cutting simulations involving porous, cellular structures relies heavily on choosing correct constitutive material models that account for the complex yield behaviour of such materials.
While the modelling framework has provided very good predictions of overall cutting behaviour, there are certain limitations within this study that should be acknowledged. Firstly, the out-of-plane thickness of the model was assumed to be represented by a single element. While this represents a simplification of the 3D geometry, the orthogonal cutting process presented here is fully planar, implying that out-of-plane forces and effects are negligible. Moreover, while the initial post-yield behaviour of the computational model was determined based on uniaxial compression testing, subsequent damage (or softening) behaviour was described by a damage evolution law (Equation (8)), which is formulated as a function of the equivalent plastic strain (PEEQ) at failure. Due to the difficulty in experimentally determining the individual damage parameters (d 1 -d 3 ) for this evolution law, our study used an iterative process to calibrate a set of damage constants until an appropriate representation of chip formation compared with the experiments was observed. This process was not trivial, as the computational model was required to capture many different categories of chip formation. Considering the complex dynamic failure processes taking place ahead of the tooltip during orthogonal cutting, the calibrated parameters correctly predicted chip formation in all but one of the cases examined. It should be noted that in the case where the model failed to capture chip formation (α = 20 • , h = 1 mm), several different types of chip formation were observed experimentally (Type 3a and 3b), suggesting that there may have been a transition in chip formation occurring for these set of parameters. An alternative reason could be the use of a continuum-based modelling approach for the PU foam material. While this study represents the first implementation of a CFI plasticity model for this type of foam structure, this constitutive formulation may have difficulty in capturing certain features occurring at particular length scales within this cellular solid. It is worth noting that the PU20 foam used is described by the manufacturer as a closed-cell rigid cellular foam. While this material has been shown to have similar bulk post-yield mechanical behaviour as trabecular bone, the characteristic cell size is much smaller (~0.17 mm) than trabecular bone, whose pore dimensions would be closer to the cutting depths used in this study (1-2 mm). At this scale, it may be the case that a discrete modelling approach in which the cell walls themselves are represented may be required to capture these intricate features of the cutting process.

Conclusions
In summary, the study presents a systematic experimental-computational investigation of how key cutting parameters govern orthogonal cutting performance. A novel finite element model was developed that combined a CFI plasticity model and a damage evolution law to predict chip formation and cutting forces in high-speed cutting simulations of trabecular bone analogues. This modelling framework predicts chip formation and cutting forces in a high-speed cutting model of trabecular analogues. The modelling framework developed has the potential to be extended to three-dimensional situations and could be used for the design of cutting tools in more operations such as burring and drilling.