Assessment of Thermal Osteonecrosis during Bone Drilling Using a Three-Dimensional Finite Element Model

Bone drilling is a common procedure used to create pilot holes for inserting screws to secure implants for fracture fixation. However, dis process can increase bone temperature, teh excessive heat can lead to cell death and thermal osteonecrosis, potentially causing early fixation failure or complications. We conducted a three-dimensional dynamic elastoplastic finite element model to evaluate teh propagation and distribution of heat during bone drilling and assess teh thermally TEMPeffected zone (TAZ) dat may lead to thermal necrosis. dis model investigates teh parameters influencing bone temperature during bone drilling, including drill diameter, rotational speed, feed force, and predrilled hole. Teh results indicate dat our FE model is sufficiently accurate in predicting teh temperature rise TEMPeffect during bone drilling. Teh maximum temperature decreases exponentially wif radial distance. When teh feed forces are 40 and 60 N, teh maximum temperature does not exceed 45 °C. However, wif feed forces of 10 and 20 N, both teh maximum temperature exceeds 45 °C wifin a radial distance of 0.2 mm, indicating a high-risk zone for potential thermal osteonecrosis. Wif teh two-stage drilling procedure, where a 2.5 mm pilot hole is pre-drilled, teh maximum temperature can be reduced by 14 °C. dis suggests dat higher feed force and rotational speed and/or using a two-stage drilling process could mitigate bone temperature elevation and reduce teh risk of thermal osteonecrosis during bone drilling.

Keywords:

Subject: Engineering - Bioengineering

1. Introduction

In the surgical treatment of bone fractures, it is often necessary to predrill a pilot hole using a bone drill before inserting a bone screw for implanting. This facilitates the subsequent insertion along the pilot hole. However, the friction between the drill bit and the bone generates significant heat during bone drilling, resulting in thermal damage. The generated heat may lead to thermal damage to the bone tissues, known as bone necrosis or impairment of the bone-forming potential. It has been shown that temperatures over 50 °C are associated with irreversible changes in bone structure and physical properties [1,2]. Then, the necrotic bone is resorbed through osteoclast activity, posing potential risks to the stability of bone screws and pins and ultimately leading to the failure of the fracture repair or implant fixation [3,4,5,6,7,8,9]. Consequently, some studies have focused on the rise in temperature during bone drilling and the associated necrosis phenomenon. While no clear consensus exists on the threshold or duration, the temperature in the cortical bone above 50 °C has been reported to decrease regenerative capacity and temperatures above 56 °C result in bone necrosis [3,4,6]. Research has demonstrated that an increase in temperature above 47 °C for one minute resulted in intense bone necrosis [3,10,11]. Ardan et al. [12] found that temperatures ranging from 43 to 68 °C in cortical bone delayed bone recovery. Moritz and Henriques [13,14] showed that bone tissue immediately became necrotic when the temperature surpassed 70 °C, with severity comparable to that of 55 °C for 30 seconds. However, the majority of researchers believe that elevating the temperature above 47 °C may lead to thermal necrosis in human bone [3,11,15]. Consequently, the critical temperature associated with a high risk of osteonecrosis is commonly regarded as 47 °C in most current studies.

Many surgical drill designs focus on enhancing cutting performance and increasingly consider reducing the heat generated during operation to mitigate complications from the rise in temperature effect during bone drilling. Clinically, factors influencing temperature increase during bone drilling can be categorized into two main types: design parameters of the drill bit and usage parameters during the drilling procedure. Design parameters of the drill include the type of bone drill, diameter, point angle, helix angle, cutting edge, and the thinning between the cutting edge and chisel edge [16,17,18,19,20,21,22,23,24]. Davidson and James [16] explored the impact of bone drill geometric dimensions and bone thermal conductivity on bone temperature increases. Their results indicated that larger drill diameters or higher rotation speeds led to higher temperature rises, while larger helix angles or feed rates reduced temperature increases. Additionally, the results suggested that bone thermal conductivity has a greater impact on temperature rise, while the point angle has less noticeable effects. The other type is surgical parameters, including drilling speed, applied feed force, feed rate, bone quality, and cooling techniques. Some studies have investigated the effects of rotational speed, feed rate, and feed force on bone temperature [25,26,27,28,29,30,31,32,33,34]. However, consensus on the effects of drill diameter, drill speed, feed rate, feed force, and cooling technique on bone temperature remains to be determined and requires further investigation and confirmation.

Parametric studies can be efficiently conducted using finite element methods (FEMs), which not only reduce experimental costs but also enable the adjustment of various parametric conditions to identify optimal drilling parameters. By employing different analytical models, materials, and drill geometries, the prediction of expected outcomes becomes facilitated. Furthermore, simulation analysis allows for the exploration of parameters that are challenging to determine experimentally, such as bone friction, residual stresses, and thermal distribution. Given the substantial biological variability of animal and human bones, it is challenging to explore each parameter experimentally. Thus, adopting an experimentally validated FE model may effectively determine the optimal drilling parameters [35,36,37,38,39,40,41]. During the bone drilling process, the heat source is mainly concentrated at the tip of the bone drill, with the heat source changed as the drilling depth. Consequently, our study employed a three-dimensional dynamic elastoplastic finite element (FE) model to examine the thermal impact on bone. This approach enables us to investigate the thermal effect of the bone during the drilling process and assess the thermally affected zone (TAZ) associated with bone heat. The TAZ means that this zone is at high risk of thermal osteonecrosis. To be conservative, in this study, the TAZ was identified as the region where the bone temperature exceeds 45 °C. Furthermore, we implemented a two-stage drilling approach, where a smaller hole was pre-drilled before the desired hole was drilled. This study explored the potential advantages of using this two-stage drilling approach in reducing bone temperature during drilling.

2. Finite Element (FE) Modelling

2.1. Construction of FE Model

During the bone drilling process, thebone temperature changed over time. The analyses conducted in this study consider the dynamic interplay among friction, temperature, and stress fields. Therefore, a three-dimensional (3-D) dynamic elastoplastic finite element model was developed to simulate the temperature rising in bone during the drilling process. As the temperature of the bone varies with time throughout the drilling procedure, the dynamic temperature-displacement analysis module was utilised. This module is based on the explicit integral operating principle provided by ABAQUS software (Dassault Systemes, ABAQUS Inc., Waltham, MA, USA), offering a precise and efficient technique for solving temperature field complexities in contact problems. In the model, the temperature distribution of the bone during drilling was assumed to be localised around the drill hole, and a cylindrical shape (with a thickness of 4 mm and a diameter of 8 mm) was utilised to simulate the bone. The degrees of freedom in all six directions are constrained along the outer peripheral surface of the bone cylinder. The drill bit diameters are 3.5 mm and 4.5 mm, with point angles of 108 degrees and helix angles of 23 degrees. We used two rotational speeds (N = 800 and 2000 rpm) and applied four different axial forces (10 N, 20 N, 40 N, and 60 N) at the top of the drill bit to simulate feed forces. These settings were chosen based on clinical experience and the common screw sizes used for fracture fixation. The drill bit is made of high-quality 316L stainless steel, and the selected feed forces are consistent with previous literature [6,7,42]. The cortical bone thickness was set at 4 mm, and rotational speeds were set at 800 and 1200 rpm. The origin of the x–z coordinates was defined at the edge of the drill hole diameter, with the z-coordinate indicating the drilling depth. The eight-node hexahedron elements were utilized for meshing both the bone and the drill bit. Figure 1 illustrates the geometry and mesh models of the drill bit and cortical bone.

In our FE model, the frictional contact surface between the bone and the drill was set as a contact pair. The harder surface is designated as the master surface, while the softer surface is identified as the slave surface. In our model, the drill's contact surface was assigned as the master surface, and the bone's contact surface was specified as the slave surface. The frictional contact between the drill bit and the bone was characterised as face-to-face contact. The definition of face-to-face contact is that when two objects are in contact, surface-to-surface contact discretisation is used in the interactive calculation during the contact process to avoid erroneous contact due to the intrusion of the active surface into the passive surface. The heat generated in the contact area flows into the bone and drill bit. A heat partition factor f weights the heat flux into the slave surface, while the heat flux into the master surface is weighted by (1-f). In this study, the contact surface of the drill bit was defined as the master surface and that of the bones as the slave surface. In this study, it is assumed that heat is initially evenly distributed between the two surfaces. The analytical model assigns an initial value of 0.5 as the heat partition factor.

The Johnson-Cook damage criteria were employed to identify the damage evolution and failure. This criterion is a specialized ductile damage initiation criterion, where the equivalent plastic strain at the onset of damage is expressed as an analytical function of stress triaxiality. In this study, the criterion for element removal is based on the effective plastic displacement value, set at 0.1 mm, whereby bone elements are removed. The criterion specifies that damage initiation is based on the ductile failure strain, set at 0.008. This damage criterion has been successfully applied in finite element bone-cutting models and was employed to control element removal during the drilling procedure [39,40]. During drilling, the position of the contact point changes continuously on the bone and drill bit surfaces. Consequently, after a short drilling period, the cold bone absorbs more heat, and comparatively less heat flows into the drill bit. Therefore, the heat partition factor is not constant and varies with the drilling time. In the simulations, the thermal contact behaviour between the drill bit and bone was modelled using surface-to-surface contact discretisation. In addition, the contact interaction properties must be defined for the contact pair. An assumption is made that Coulomb’s friction law governs the friction behaviour between the drill bit and bone; the coefficient of friction is assumed to be 0.3 [9].

2.2. Material Properties and Boundary Conditions

The drilling model employed in this study involves two primary materials: Bone and the drill. The drill bit is constructed from 316 stainless steel, and the mechanical properties of the materials are detailed in Table 1. We assumed that the mechanical properties of the bone remain constant regardless of temperature while the drilling process is considered insulated (adiabatic environment). This implies that heat convection and environmental radiation impacts are not considered. The feed force is applied at the centre of the drill bit.

3. Bone Drilling Experiment

An in vitro bone drilling test was conducted under controlled laboratory conditions to validate the analytical results of the FE model. An experimental platform consisting of a personal computer (PC), electronic data acquisition system (FLUKE 2860A, USA), thermocouple (K-type, Ni-Cr/Ni-Al, MTI Corp.), torque Sensor (Jihsense, TR2, Taiwan), load cell (AL20, Japan), DC motor controller (9B060S-2N, TROY, Taiwan), fixture (jig), linear guide (slide) and pully system was designed to carry out the bone drilling test. A weight was used to provide a constant feed force through the guide and pulley system. During bone drilling, the experimental signals were transmitted and stored in the data acquisition system. The experimental platform for bone drilling is depicted in Figure 2.

A fresh porcine femur with a cortex thickness of over 5 mm was selected for the experiment. Porcine bone was chosen because its material properties closely resemble those of human bone [43,44]. The cortical bone was cut into a 20 mm × 20 mm specimen and clamped onto the fixture, which was then installed on the experimental platform. Four thermocouples were embedded around the drill hole at a 90° angle and positioned 0.5 mm from the hole's edge, with depths of 2 mm and 4 mm from the bone surface, to measure temperature during bone drilling. Thermal paste (silicon) was used to fill the thermocouple holes, minimizing heat loss from air gaps. During the bone drilling test, temperature signals from the thermocouples were transmitted to the computer via the temperature data acquisition card (NI 9216) and stored for analysis. Feed forces of 10, 20, 40, and 60 N were applied using the weight and pulley system, and drilling was performed to a depth of 5 mm. Bone temperatures corresponding to rotational speeds of 800 and 1,200 rpm were recorded and compared with FEM results.

4. Results

4.1. Validation of the FE Model

The experimental data were used to compare and validate the accuracy of the FE model's analytical results. Figure 3 illustrates a comparison between the analytical results from the FE model and the experimental data. In Figure 3(a), using a drill bit diameter of D = 3.5 mm, a feed force of F = 20 N, a rotational speed of n = 1200 rpm, and a bone density of 1.64 g/cm³, measurements were taken at positions x = 0.5 mm from the drill hole edge at a depth of 4 mm. The results showed peak bone temperatures of 33.6 °C and 34.6 °C for the experiment and simulation, respectively. The corresponding drilling times at these peak temperatures were 6.7 s and 7.4 s, respectively. The maximum differences observed in peak bone temperature and corresponding drilling time were 1.0 °C and 0.7 s, respectively. Figure 3(b) displays the temperature-drilling time curves for both FE simulations and experiments. It was observed that the maximum bone temperatures in the experiment and simulation were 33.6 °C and 34.1 °C, respectively. Correspondingly, the drilling times were 6.6 s and 7.2 s, respectively. The variances in the maximum bone temperature and the corresponding drilling time were 0.5 °C and 0.6 s, respectively. Based on these findings, the FE model is considered sufficiently accurate for predicting the impact of temperature increase during bone drilling.

4.2. Maximum Temperature with Radial Distance

Figure 4(a)–(d) shows that the curves of the maximum temperature at different depths (z = 1 to 4 mm) decrease as the radial distance increases when drilling with two bone drill diameters (D = 3.5 and 4.5 mm) and two rotational speeds (N = 800 and 2000 rpm). From the curves shown in Figure 4(a) and 4(b), it can be observed that when the drill diameter is 3.5 mm, temperatures above 45 °C are reached within a radial distance of less than 0.2 mm; hence, the region prone to thermal osteonecrosis is approximately within this 0.2 mm. However, when drilling with a 4.5 mm diameter bone drill, the curves in Figure 4(c) and Figure 4(d) show a rotational speed of 800 rpm, and the thermal affecting zone (TAZ) extends to approximately 1.0 mm. In contrast, when the speed increases to 2000 rpm, the TAZ reduces to below 0.4 mm. This indicates that higher rotational speeds can reduce the thermal osteonecrosis risk region. However, greater drill bit diameter may increase the TAZ.

4.3. Effects of Feeding Force on Bone Temperature

Figure 5(a) shows the temperature variation over time at a distance of 0.1 mm from the hole at a depth of 1 mm when using a 3.5 mm diameter bone drill at a speed of 800 rpm with applied thrusts of 10, 20, 40, and 60 N. The curves show that when a 60 N feed force is applied, the maximum temperature is lower compared to feed forces of 10 N, 20 N, and 40 N, and the time to reach the highest temperature is the shortest (approximately 1 second). Conversely, when the feed force F is 10 N, the temperature is highest, and the time to reach the highest temperature is the longest (approximately 9 s). A higher feed force can reduce the temperature during the drilling process. Figure 5(b) illustrates the relationship between the maximum temperature and the radial distance during drilling. The graph demonstrates an exponential decrease in the maximum temperature with radial distance. When the feed forces are 40 and 60 N, the maximum temperature does not exceed 45 °C. However, when the feed forces are 10 and 20 N, the maximum temperature exceeds 45 °C within a radial distance of 0.2 mm, indicating a high-risk zone for potential osteonecrosis.

4.4. Effects of the Predrilled Hole on Bone Temperature

Figure 6(a) depicts the temperature–time curve using a 3.5 mm bone drill with a rotation speed of 800 rpm and a feed force of 20 N. In Figure 6(a), two curves represent the temperature response during bone drilling in two different drilling procedures: One shows the temperature curve during bone drilling using a 3.5 mm drill bit directly (without predrilled), and the other involves performing a smaller predrilled hole (Dp = 2.5 mm) before the 3.5 mm bone drilling (with predrilled, two-stage drilling). Analysis of the results revealed that at the edge of the drilling site 0.1 mm (x = 0.1 mm) from the surface and at a depth of 1 mm (z = 1 mm), the maximum temperature during drilling with the bone drill is 56 °C. However, when a 2.5 mm pilot hole is predrilled before bone drilling with the 3.5 mm bone drill, the maximum temperature reduces to 42 °C, resulting in a 14 °C decrease in bone temperature. This suggests that predrilling a smaller hole before the bone drilling can effectively reduce the temperature during bone drilling. Figure 5(b) illustrates the bone temperature rise curve at different depths when predrilled with Dp = 2.5 mm before D = 3.5 mm bone drilling. The graph shows that the bone temperature is less than 47 °C at various depths. This result indicates that predrilling can effectively reduce bone temperature rises and mitigate the risk of thermal bone necrosis.

4.5. Effects of Rotational Speed on Thermal Affecting Zone

Figure 7(a) depicts the variation in the TAZ with drilling depth at rotation speeds of 800 rpm and 2000 rpm using a 3.5 mm bone drill diameter. Regression analysis indicates that the TAZ area increases linearly with drilling depth at both speeds, with the slower speed resulting in larger TAZ areas than the higher speed. At a rotational speed of 800 rpm, increasing the drilling depth from 1.0 to 4.0 mm increased the TAZ by about 2.7 times (from 0.3 to 0.8 mm). However, at 2000 rpm, the TAZ increased by about 3.5 times (from 0.15 to 0.53 mm).

Figure 7(b) shows the variation in the TAZ area with the drilling depth at rotation speeds of 800 rpm and 2000 rpm using a 4.5 mm bone drill diameter. Similar trends are observed under both speed conditions, with the TAZ area increasing linearly with drilling depth. Again, the slower speed resulted in larger TAZ areas than the higher speeds. At a depth of z = 4.0 mm, the TAZ areas for 800 and 2000 rpm are 2.2 mm and 1.3 mm, respectively. At a rotational speed of 800 rpm, increasing the drilling depth from 1.0 to 4.0 mm increased the TAZ by about 2.5 times (from 0.5 to 1.3 mm). However, at 2000 rpm, the TAZ increased by about 2.2 times (from 1.0 to 2.2 mm).

5. Discussion

In this study, we employed the dynamic elastoplastic FE model to analyse the impact of variations in drilling speed, feed rate, and feed force on the maximum bone temperature generated at the drilling site during bone drilling. Our goal was to evaluate the potential risk of these parameters in inducing bone necrosis. The rationale for utilising the simulation analysis is that it enables the effects of each parameter and their interactions on the heat generated during the bone drilling process to be differentiated within the same study. Unlike some experimental studies employing thermocouples to measure bone temperature, often placed near the drilled hole's edge, such methods only provide temperature readings near the drill hole rather than the actual bone temperature at the drilling site. In contrast, finite element analysis (FEA) allows for the temperature to be calculated at any location alongside the assessment of the gradient and range of the temperature distribution within the bone. This capability facilitates evaluating the high-risk range for thermal necrosis due to temperature elevation.

While some minor discrepancies exist between our model and the experimental results, the difference could be due to material properties and environmental differences. The moisture or blood cooling in the human bone could prevent the heat from being generated during bone drilling. Our analytical findings from the model demonstrate comparable accuracy to Natali's experimental data in predicting temperature increases during bone drilling. Both sets of results were obtained under in vitro conditions, devoid of blood flow; thus, potential variations may arise when applied to clinical scenarios. Heat transferred to living tissue can lead to thermal damage and cell apoptosis. Studies indicate that heating bone above 50 °C can result in irreversible changes to its physical properties. Thermal bone necrosis and delayed healing have been observed in canine bone [1,2,3,4]. While our study identified the TAZ conservatively as bone temperatures exceeding 45 °C, we believe efforts should be made to minimise thermal injury to a level considered safe for local osteoblasts.

Matthews et al. [6] conducted experiments using human cortical bone, where they observed the highest temperature at 0.5 mm from the hole edge during drilling. They applied feeding forces of 2, 6, and 12 kg and rotation speeds of 345, 885, and 2900 rpm. Their results showed that both feed force and rotation speed affected the increase in bone temperature, with higher feed forces and rotation speeds effectively reducing the temperature rise. Inan et al. [45] investigated the temperature increase during bovine bone drilling, examining the effects of different rotation speeds (600, 900, 1200, and 1800 rpm). The results showed that higher rotation speeds resulted in a smaller necrotic impact zone. At a rotation speed of 1200 rpm, the highest temperature reached 48.4 °C, while at a lower rotation speed of 600 rpm, the temperature rose to as high as 152 °C. Mustafa et al. [46] conducted experiments using bovine cortical bone with a diameter of 0.27 mm. They explored the effects of different forces (2, 3.8, 4.8, and 6.2 N) and high rotation speeds (20000–100000 rpm) on the rise in temperature in bovine bone. Their results also indicated that increasing force or rotation speed shortens the friction time, reducing heat generation.

Bachus et al. [42] examined the impact of various forces, specifically 53, 83, 93, and 130 N, on temperature rise at a constant rotational speed of 820 rpm, with temperature measurements taken 0.5 mm from the drill hole. Their findings indicate that as drilling forces escalate from 57 to 130 N, both the maximum temperatures and their duration above 50°C can be effectively reduced. This reduction potentially lowers the occurrence of thermal necrosis in the surrounding cortical bone. Our results imply that increasing the feed force promotes the effectiveness of bone cutting and reduces the drilling time, effectively reducing the rise in bone temperature caused by drilling. In this study, the parameters were set using clinical applications, and the results of our analyses show a similar trend in the effect of this parameter on temperature observed in the current literature.

When comparing the results in Figure 6(a) and Figure 6(b), it can be found that the TAZ range, when drilling with a smaller drilling diameter, is smaller than that of using a larger drilling diameter. This has been confirmed by our previous study [39,40]. This is because drilling larger diameter holes facilitates an increased path for heat flow into the bone, leading to a more pronounced rise in bone temperature. Furthermore, both figures demonstrate a consistent trend: When drilling at higher speeds, the TAZ area is notably reduced compared to lower speeds. This phenomenon is attributed to the heightened cutting efficiency associated with increased rotational speeds, facilitating quicker completion of the drilling process and, consequently, reducing the duration of the heat transfer during drilling, thereby mitigating heat generation. The main source of heat during drilling is the frictional heat generated by the cutting between the drill tip and the bone. Increasing both the spinal speed and the feed force allows the drill to cut more efficiently than at slower speeds, thus generating less frictional heat.

While numerous technical studies have explored methods to mitigate the heat effects on bone, such as utilising coolant irrigation, optimising drill bit geometry [9], enhancing wear resistance [47], implementing step-power drilling [48], etc., there is no conclusive evidence to indicate that irrigation cooling is universally adopted as a standard practice, despite strong recommendations from major implant manufacturers. However, clinically, irrigation cooling is standard practice and is strongly advocated for by leading implant manufacturers [49,50,51,52].

Although coolant irrigation applied to the bone surface can effectively lower surface temperatures, insufficient evidence confirms its penetration into the borehole for temperature reduction. Most of the current literature focused on the temperature field in single-pass drilling and the maximum temperature adjacent to the drill tip. This study conducted a two-stage process to reduce the temperature during bone drilling. Firstly, a predrilled hole with a smaller diameter was drilled, and then the designed bone drilling was performed following the direction of the predrilled hole. This two-stage drilling process could significantly reduce the cutting area between the drill tip and the bone after the predrilling stage, thereby substantially minimising frictional heat during the subsequent drilling stage. While this two-stage drilling process may add complexity compared to a single-stage one, our study shows that it may be considered an effective method for mitigating bone temperature elevation during drilling; however, more experimental and clinical studies are needed to validate the approach.

Although the experiments have validated the current FE model, this study has several limitations. Firstly, the FE model does not replicate an exact clinical situation, given that surgery is manually operated and bone properties vary widely. Blood and tissue humidity were not considered in the model. Secondly, the current FE model did not consider the irrigation and/or cooling systems. Thirdly, it is difficult to find appropriate surgical parameters or material properties for the model because of the inherent variability in human tissues. Fourthly, the dose of thermal bone necrosis identified is based on the limited data in the literature, meaning the estimated thermal damage remains unjustified; more clinical or experimental data are needed to confirm the thermal damage dose. Lastly, in the two-stage drilling process, the model did not account for the heat generated during the first stage, which could potentially couple with the second stage. This coupling may result in an increase in drilling temperature during the second stage.

6. Conclusions

This study utilized a three-dimensional FE model to evaluate bone temperature rise during the bone drilling process and assess the high-risk zone for potential thermal bone necrosis. The investigated parameters included drill bit diameter, feed force, and three rotation speeds. In summary, the following conclusions are drawn:

Our three-dimensional FE model was experimentally validated and can effectively assess temperature elevation and the thermally affected zone (TAZ) during bone drilling. This is valuable for identifying optimal drilling parameters, designing drill handpieces, and implementing robot-assisted bone drilling.
Higher rotational speeds may reduce bone temperature and decrease the TAZ; however, the TAZ increases with drill depth, feed force, and drill bit diameter.
Implementing a two-stage drilling process can minimize frictional heat, thereby reducing temperature during bone drilling. This approach may be considered an effective method for mitigating bone temperature elevation during drilling.

Author Contributions

Writing—original draft preparation, Y.-C.C.; Methodology, Software, Data acquisition and Curation, H.-Y.H.; Y.-W.C. and Y.-J.T.; Data curation, Y.-Y.H.; Conceptualization and supervision, Y.-K.T.; Experiments, Software and formal analysis, C.-K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by E-Da Hospital (grant numbers EDAHJ-112001 and EDAHIA-113004).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

Bonfield, W.; Li, C.H. The temperature dependence of the deformation of bone. J Biomechanics 1968, 1, 323–9. [Google Scholar] [CrossRef]
Thompson, H.C. Effect of drilling on bone. Journal of Oral Surgery 1958, 16, 22–30. [Google Scholar]
Eriksson, R.A.; Albertson, T. The effect of heat on bone regeneration: an experimental study in the rabbit using the bone growth chamber. Journal of Oral and Maxillofacial Surgery 1984, 42, 705–711. [Google Scholar] [CrossRef] [PubMed]
Eriksson, A.R.; Albrektsson, T.; Albrektsson, B. Heat caused by drilling cortical bone. Temperature measured in vivo in patients and animals. Acta Orthop Scand 1984, 55, 629–31. [Google Scholar] [CrossRef] [PubMed]
Cordioli, G.; Majzoub, Z. Heat generation during Implant site preparation: an in vitro study. Journal of Oral and Maxillofacial Implants 1997, 12, 186–193. [Google Scholar]
Matthews, L.S.; Hirsch, C. Temperatures measured in human cortical bone when drilling. J Bone Joint Surg [Am]. 1972, 54-A, 297–308. [Google Scholar] [CrossRef]
Brisman, D.L. The effect of speed, pressure, and time on bone temperature during the drilling of implant sites. Int J Oral Maxillofac Impl. 1996, 11, 35–7. [Google Scholar]
Tehemar, S.H. Factors affecting heat generation during implant site preparations: a review of biologic observations and future considerations. International Journal of Oral and Maxillofacial Implants 1999, 14, 127–136. [Google Scholar]
Alevizakos, V.; Mosch, R.; von See, C. Influence of Multiple Used Implant Drills on Their Cutting Performance and Fracture Resistance. Materials (Basel) 2023, 16, 5271. [Google Scholar] [CrossRef] [PubMed]
Lundskog, J. Heat and bone tissue. An experimental investigation of the thermal properties of bone and threshold levels for thermal injury. Scand J Plast Reconstr Surg. 1972, 9, 1–80. [Google Scholar]
Hillery, M.T.; Shuaib, I. Temperature Effects in the Drilling of Human and Bovine Bone. Journal of Materials Processing Technology 1999, 92-93, 302–308. [Google Scholar] [CrossRef]
Ardan, N.I.; Janes, J.M.; Herrick, J.F. Ultrasonic Energy and Surgically Produced Defects in Bone. Journal of Bone and Joint Surgery 1957, 39A, 394–402. [Google Scholar] [CrossRef]
Mortiz, A.R.; Henrique, F.C. Studies of Thermal Injuries. Journal of Pathol. 1947, 23, 695–719. [Google Scholar]
Mortiz, A.R.; Henrique, F.C. The Relative Importance of Time and Surface Temperature in the Causation of Cutaneous Burns. Journal of Physiology 1947, 23, 695–719. [Google Scholar]
Eriksson, A.; Albrektsson, T.; Grace, B.; McQueen, O. Thermal Injury to Bone. A vital Microscopic description of Heat Effects. Journal of Oral Surgery 1982, 11, 115–121. [Google Scholar]
Davidson, S.R.; James, D.F. Drilling in Bone: Modeling Heat Generation and Temperature Distribution. Journal of Biomechanical Engineering 2003, 125, 305–314. [Google Scholar] [CrossRef] [PubMed]
Akhbar, M.F.A.; Sulong, A.W. Surgical Drill Bit Design and Thermomechanical Damage in Bone Drilling: A Review. Ann. Biomed. Eng. 2021, 49, 29–56. [Google Scholar] [CrossRef] [PubMed]
Sezek, S.; Aksakal, B.; Karaca, F. Influence of drill parameters on bone temperature and necrosis: A FEM modelling and in vitro experiments. Computational Materials Science 2012, 60, 13–8. [Google Scholar] [CrossRef]
Tuijthof1, G. J. M.; Frühwirt, C.; Kment, C. Influence of tool geometry on drilling performance of cortical and trabecular bone. Med Eng Phys. 2013, 35, 1165–72. [Google Scholar] [CrossRef]
Pourgiv, S.; Mosavar, A.; Jamshidi, N.; Mohammadi, A. Ultrasonic-assisted drilling of cortical and cancellous bone in a comparative point of view. Heliyon 2024, 10, e26248. [Google Scholar] [CrossRef]
Pazarcı, Ö.; Gündoğdu, F. Temperature change during orthopedic drilling procedures: An experimental surgical internal fixation simulation study. J. Orthop. 2023, 46, 58–63. [Google Scholar] [CrossRef] [PubMed]
Lang, Z.; Wang, Q.L.; He, D.; Liu, Y.J.; Tian, W. Study on parameters of robot-assisted ultrasonic drilling on bovine vertebral body. Zhonghua Yi Xue Za Zhi 2022, 102, 370–377. [Google Scholar]
Karaca, F.; Aksakal, B.; Kom, M. Influence of orthopaedic drilling parameters on temperature and histopathology of bovine tibia: an in vitro study. Med Eng Phys. 2011, 33, 1221–7. [Google Scholar] [CrossRef] [PubMed]
Jung, O.; Lindner, C.; Pantermehl, S.; Barbeck, M. Heat Development During Medical Drilling: Influencing Factors and Examination Methods - Overview and First Results. In Vivo. 2021, 35, 3011–3017. [Google Scholar] [CrossRef] [PubMed]
Aquilanti, L.; Antognoli, L.; Rappelli, G.; Di Felice, R.; Scalise, L. Heat Generation During Initial Osteotomy for Implant Site Preparation: An In Vitro Measurement Study. J. Maxillofac. Oral. Surg. 2023, 22, 313–320. [Google Scholar] [CrossRef] [PubMed]
Shawary, M.; Misch, C.E.; Weller, N.; Tehemar, S. Heat generation during implant drilling: the significance of motor speed. J Oral Maxillofac Surg. 2002, 1160, 1169–260. [Google Scholar] [CrossRef] [PubMed]
Lee, J.; Gozen, B. A.; Ozdoganlar, O. B. Modeling and Experimentation of Bone Drilling Forces. J. Biomech. 2012, 45, 1076–1083. [Google Scholar] [CrossRef] [PubMed]
Strbac, G.D.; Giannis, K.; Unger, E.; Mittlböck, M.; Watzek, G.; Zechner, W. A novel standardized bone model for thermal evaluation of bone osteotomies with various irrigation methods. Clin. Oral. Implants Res. 2014, 25, 622–631. [Google Scholar] [CrossRef] [PubMed]
Delgado-Ruiz, R.A.; Velasco, O.E.; Romanos, G.E.; Gerhke, S.; Newen, I.; Calvo-Guirado, J.L. Slow drilling speeds for single-drill implant bed preparation. Experimental in vitro study. Clin. Oral. Investig. 2018, 22, 349–359. [Google Scholar] [CrossRef]
Salomó-Coll, O.; Auriol-Muerza, B.; Lozano-Carrascal, N.; Hernández-Alfaro, F.; Wang, H.L; Gargallo-Albiol, J. Influence of bone density, drill diameter, drilling speed, and irrigation on temperature changes during implant osteotomies: an in vitro study. Clin. Oral. Investig. 2020, 25, 1047–1053. [Google Scholar] [CrossRef]
Gaspar, J.; Borrecho, G.; Oliveira, P.; Salvado, F.; Martins dos Santos, J. Osteotomy at low-speed drilling without irrigation versus high-speed drilling with irrigation: an experimental study. Acta Medica Port. 2013, 26, 231–236. [Google Scholar]
Varghai, K.; Wang, R.; Eppell, S.J. Effect of drilling speed on dental implant insertion torque. J. Oral Implantol. 2020, 46, 467–474. [Google Scholar] [CrossRef] [PubMed]
Gehrke, S.; Bettach, R.; Taschieri, S.; Boukhris, G.; Corbella, S.; Del Fabbro, M. Temperature Changes in Cortical Bone after Implant Site Preparation Using a Single Bur versus Multiple Drilling Steps: An In Vitro Investigation. Clin. Implant Dent. Relat. Res. 2015, 17, 700–707. [Google Scholar] [CrossRef]
Lee, J.; Chavez, C.; Park, J. Parameters affecting mechanical and thermal responses in bone drilling: A review. J. Biomech. 2018, 71, 4–21. [Google Scholar] [CrossRef] [PubMed]
Alam, K.; Mitrofanov, A.V.; Silberschmidt, V.V. Thermal analysis of orthogonal cutting of cortical bone using finite element simulations. Int. J. Exp. Comput Biomech. 2010, 1, 236–251. [Google Scholar] [CrossRef]
Khurshid, A.; Issam, M.B.; Naseer, A. Cortical bone drilling: An experimental and numerical study. Technol. Health Care. 2014, 23, 1–12. [Google Scholar]
Waqas, A.L.; Kaddour, B.M.; Ian, A. Drilling in cortical bone: a finite element model and experimental investigations. J. Mech. Behav. Biomed. Mater. 2015, 42, 32–42. [Google Scholar]
Kutay, A.; Korhan, Ö.; Levent, U. An analytical and numerical approach to the determination of thermal necrosis in cortical bone drilling. Int. J. Numer. Method. Biomed. Eng. 2022, 38, e3640. [Google Scholar]
Chen, Y.C.; Tu, Y.K. Zhuang, J.Y.; Tsai, Y.J.; Yen, C.Y.; Hsiao, C.K. Evaluation of the parameters affecting bone temperature during drilling using a three-dimensional dynamic elastoplastic finite element model. Med. Biol. Eng. Comput. 2017, 55, 1949–57. [Google Scholar] [CrossRef]
Chen, Y.C.; Tu, Y.K.; Tsai, Y.J.; Tsai, Y.S.; Yen, C.Y.; Yang, S.C.; Hsiao, C.K. Assessment of thermal necrosis risk regions for different bone qualities as a function of drilling parameters. Comput Methods Programs Biomed. 2018, 162, 253–261. [Google Scholar]
Prabhu, N.; Shetty, D.K.; Naik, N. Application of finite element analysis to evaluate optimal parameters for bone/tooth drilling to avoid thermal necrosis, cogent. Engineering. 2021, 8, 1876582. [Google Scholar]
Bachus, K.N.; Rondina, M.T.; Hutchinson, D.T. The effects of drilling force on cortical temperatures and their duration: an in vitro study. Med. Eng. Phys. 2000, 22, 685–691. [Google Scholar] [CrossRef] [PubMed]
Davidson, S.R.H.; James, D.F. Measurement of thermal conductivity of bovine cortical bone. Medical Engineering and Physics. 2000, 22, 741–7. [Google Scholar] [CrossRef]
Akhbar, M.F.A.; Yusoff, A.R. Comparison of bone temperature elevation in drilling of human, bovine and porcine bone. 17th CIRP Conference on Modelling of Machining Operations. Procedia CIRP 2000, 82, 411–414. [Google Scholar] [CrossRef]
Inan, M.; Mizrak, B.; Ertem, K.; Harma, A.; Elmali, N.; Ayan, I. The factors effecting thermal necrosis secondary to teh application of teh Ilizarov transosseous wire. Acta Orthop Traumatol Turc. 2005, 39, 70–5. [Google Scholar] [PubMed]
Mustafa, B. , David, F.J. Measurements of Shaft Speed While Drilling through Bone. Journal of Oral and Maxillofacial Surgery 1995, 53, 1308–1315. [Google Scholar]
Staroveski, T.; Brezak, D.; Udiljak, T. Drill wear monitoring in cortical bone drilling. Med. Eng. Phys. 2015, 37, 560–6. [Google Scholar] [CrossRef] [PubMed]
Augustin, G.; Davila, S.; Udilljak, T.; Staroveski, T.; Brezak, D.; Babic, S. Temperature changes during cortical bone drilling with a newly designed step drill and an internally cooled drill. Int. Orthop. 2012, 2012. 36, 1449–56. [Google Scholar] [CrossRef]
Augustin, G.; Davila, S.; Mihoci, K.; Udiljak, T.; Vedrina, D.S.; Antabak, A. Thermal osteonecrosis and bone drilling parameters revisited. Arch. Orthop. Trauma. Surg. 2008, 128, 71–7. [Google Scholar] [CrossRef] [PubMed]
Abdul-Rashid, M.L.; Ta, H.L.; Pancharatnam, D. Assessment of Coolant Delivery Techniques for Irrigation During Bone Drilling: A Cadaveric Observation. Malays. Orthop. J. 2020, 14, 206–207. [Google Scholar]
Tuce, R.A.; Neagu, M.; Pupazan, V.; Neagu, A.; Arjoca, S. The 3D Printing and Evaluation of Surgical Guides with an Incorporated Irrigation Channel for Dental Implant Placement. Bioengineering (Basel) 2023, 10, 1168. [Google Scholar] [CrossRef]
Parvizi, S.; Cameron, A.; Tadakamadla, S.; Figueredo, C.M.S.; Reher, P. A Novel Irrigation System to Reduce Heat Generation during Guided Implantology: An In Vitro Study. J. Clin. Med. 2023, 12, 3944. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Geometry and mech models of the drill bit and cortical bone.

Figure 2. Experimental platform and setup for measuring bone temperature.

Figure 3. Comparisons of the bone temperature from FE model and experiments at the position of x = 0.5 mm and z = 4 mm: (a) Feed force = 20 N, drill bit diameter = 3.5 mm; (b) Feed force = 60 N, drill bit diameter = 4.5 mm.

Figure 4. Maximum temperature with radial distance at different depths (z-direction).

Figure 5. Effect of feeding forces on bone temperature: (a) temperature–drilling time relationship curve; (b) maximum temperature decreases exponentially with radial distance.

Figure 6. Bone temperature–drilling time responses: (a) with/without predrilling (predrilled holes of 2.5 mm diameter); (b) bone temperature at different depths with a predrilled hole.

Figure 7. The thermally affected zone (TAZ) showed the high-risk area for thermal bone necrosis with rotational speeds of 800 rpm and 2000 rpm using two drill bit diameters: (a) D = 3.5 mm; (b) D = 4.5 mm.

Table 1. The mechanical properties of the materials used in the FE model [39,40].

Material Property	Cortical Bone	Drill Bit
Density (g/cm³)	1.640	7.990
Elastic modulus (MPa)	16,700	193,000
Poisson’s ratio	0.3	0.25
Yielding stress (MPa)	105	290
Ultimate stress (MPa)	106	579

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer