Next Article in Journal
Conceptualizing Scaffold Guided Breast Tissue Regeneration in a Preclinical Large Animal Model
Previous Article in Journal
Orthobiologic Management Options for Degenerative Disc Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Yung-Chuan Chen
1,
Yi-Jung Tsai
2,†,
Hao-Yuan Hsiao
3,4,†,
Yen-Wei Chiu
2,3,†,
You-Yao Hong
1,
Yuan-Kun Tu
4,* and
Chih-Kun Hsiao
2,4,5,*
1
Department of Vehicle Engineering, National Pingtung University of Science and Technology, Pingtung 912301, Taiwan
2
Department of Medical Research, E-Da Hospital, I-Shou University, Kaohsiung 82445, Taiwan
3
Institute of Medical Science and Technology, National Sun Yat-sen University, Kaohsiung 804201, Taiwan
4
Department of Orthopedics, E-Da Hospital, I-Shou University, Kaohsiung 82445, Taiwan
5
Department of Mechanical Engineering/Graduate Institute of Mechatronics Engineering, Cheng Shiu University, Kaohsiung 833301, Taiwan
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Bioengineering 2024, 11(6), 592; https://doi.org/10.3390/bioengineering11060592
Submission received: 1 May 2024 / Revised: 19 May 2024 / Accepted: 5 June 2024 / Published: 10 June 2024
(This article belongs to the Section Biomedical Engineering and Biomaterials)

Abstract

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

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]. The necrotic bone is resorbed through osteoclast activity, posing potential risks to the stability of bone screws and pins and ultimately leads 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 the severity comparable to that of 55 °C for 30 s. 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 the complications from the rise in temperature effect during bone drilling. Clinically, factors influencing the temperature increase during bone drilling can be categorised into two main types: the design parameters of the drill bit and the usage parameters during the drilling procedure. The 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 the temperature rise, while the point angle has less noticeable effects. The other type is the 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 is 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 changing with 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 on the bone during 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 predrilled 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. Temperature Distribution during Bone Drilling

The temperature field distributions in bone are governed by the conservation of heat. The temperature function (Tbone) of the bone can be expressed by the three-dimensional heat transfer equation as follows [42]:
k b [ 1 r r ( r T b o n e r ) + 1 r 2 2 T b o n e θ 2 + 2 T b o n e z 2 ] = ρ b c b T b o n e t
where k b , ρ b , and c b are thermal conductivity, bone density, and specific heat of bone, respectively. r, θ, and z express the cylindrical coordinates.
Thermo-mechanical equations have been developed from machining theory to predict the heat generated during drilling. The amount of heat generated when the bit is in frictional contact with the bone can be considered as the heat generation rate [43]:
G = q f ˙ + q ˙ p
where q ˙ f is the heat generated by the friction between the bit and bone in N-m/s, and q ˙ p is the heat generated by the deformation of the plastic during cutting (N-m/s), which can be expressed as [44]:
q ˙ p = η σ e ε ˙ p l
q ˙ f = 2 π R n μ F n
where η is the inelastic heat fraction, σ e is the effective stress (MPa), ε ˙ p l is the plastic strain rate (1/s), R is the drill radius (mm), n is the rotational speed (rpm), μ is the friction coefficient, and F n is the normal force applied during drilling (mm). The heat flux generated by the friction between the bone and drill bit can be expressed as [43]:
q f = q b o n e + q d r i l l
where q b o n e is the heat flux transmitted to the bones in W/m2, and q d r i l l is the heat flux transmitted to the drill bit in W/m2. It is assumed that the entire frictional power will be converted into heat and transferred to the drill bit and bones. If the heat partition factor is η 1 , then the heat flux transferred to the drill bit and bone is expressed as [44]:
q d r i l l = η 1 q f
q b o n e = ( 1 η 1 ) q f

2.2. Construction of FE Model

During the bone drilling process, the bone temperature changes 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 rise in bone during the drilling process. As the temperature of the bone varies with time throughout the drilling procedure, a dynamic temperature–displacement analysis module was utilised. This module is based on the explicit integral operating principle provided by ABAQUS (2023) 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,45]. 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. 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 contact surface of the drillwas assigned as the master surface, and the contact surface of the bone 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, 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.
Ductile damage criteria were employed to model the damage evolution and predict failure. In this study, one criterion specifies that damage initiation occurs at a ductile failure strain of 0.008 [31]. The criterion for element removal is based on an effective plastic displacement value of 0.3 mm, at which point bone elements are removed. This value was determined from the chip morphology of bone drilling [32]. This damage criterion has been successfully applied in finite element bone drilling models and has been used to control the removal of elements during the drilling procedure [39,40]. During drilling, the position of the contact point changes continuously on the bone and the 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]. The thermal conductivity for cortical bone and drill bit are 0.452 and 16.2 (W/m·k), respectively [39,40].

2.3. Material Properties and Boundary Conditions

The drilling model employed in this study involves two primary materials: bone and the drill bit. The drill bit is constructed from 316 stainless steel (type 310.25 and 310.31 for diameters of 2.5 and 3.5 mm, respectively), and the mechanical properties of the materials are detailed in Table 1. Feed forces of 10, 20, 40, and 60 N were applied to the centre of the drill bit. The measured room temperatures (22 and 25.5 °C) were input as the initial bone temperature for analysis. In the FE model, we assumed that almost all of the energy used in material removal is converted into heat. Bone is a poor conductor of heat; the external surface of the cortical bone is less affected by the heat flow during hole drilling, and the temperature is the same as the surrounding environment. Thus, the external surfaces of the cortical bone are considered adiabatic boundary conditions, i.e., the effects of environmental heat convection and radiation are not taken into account [12,13,14,15,16].

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. A bone drilling platform was designed to install the measurement instrumentation. It was equipped with an infrared positioning system to determine the drilling position. The experimental platform, consisting of a personal computer (PC), electronic data acquisition system (FLUKE 2860A, WA, USA. Frequency: 0.01 s), thermocouple (K-type, Ni-Cr/Ni-Al, MTI Corp. Accuracy: ±0.5 °C), torque sensor (RT25, UAS. Accuracy: 0.01%), load cell (AL50, Japan, Accuracy: 0.2%), DC motor controller (9B060S-2N, TROY, Shihlin, Taipei, Taiwan. Accuracy: ±5 rpm), fixture (jig), linear guide (slide), and pulley system, was designed to carry out the bone drilling test. The rotational speed was controlled by the DC motor, and a weight (1, 2, 4, or 6 kg) was used to provide the constant feed force (10, 20, 40, or 60 N) 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 2a.
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 [46,47]. The porcine bone was dried for 1 h to avoid any possible influence on the measured temperature of the cooling effect of the water on the bone surface. Then, the cortical bone was cut into a 20 mm × 20 mm specimen, clamped onto the fixture, and installed on the experimental platform. For bone temperature measurement, four thermocouples were embedded into the drill hole at an angle of 90° and at a distance of 0.5 mm from the edge of the drilling hole with depths of 2 mm and 4 mm from the bone surface (Figure 2b).
To validate the FE model, experiments were performed to measure bone temperature. During the experiments, the measured environmental temperatures 22.0 °C and 25.5 °C and the relative humidity (78%RH) corresponded to the initial temperatures and humidity of the bone specimens. As drilling progressed, the bone temperature changed, and the thermocouple signals were continuously transmitted to a computer system through an electronic data acquisition system. The electronic data collector can simultaneously capture temperature data from 20 thermocouples with a maximum frequency of 0.01 s. Upon completion of drilling, the drill bit was retracted to its original position, its rotation was stopped, and it was allowed to cool naturally. In the validation of the FE model, the measured temperatures, 22.0 °C and 25.5 °C, were used as the initial bone temperatures for the drilling process. Thermal paste (silicon) was used to fill the thermocouple holes, minimizing heat loss from air gaps. During the bone drilling test, the 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 1200 rpm were recorded and compared with FEM results. The results of this experiment can be used to validate the finite element model, confirming the accuracy of the simulation 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 (mean, maximum, and minimum temperature) which were obtained from 30 boreholes. In Figure 3a, 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 3b 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.5 °C and 34.0 °C, respectively. Correspondingly, the drilling times were 6.9 s and 7.5 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 the temperature increase during bone drilling.

4.2. Maximum Temperature with Radial Distance

Figure 4a–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 4a,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 4c,d show a rotational speed of 800 rpm, and the thermally affected 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, a greater drill bit diameter may increase the TAZ.

4.3. Effects of Feeding Force on Bone Temperature

Figure 5a 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 s). 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 5b 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

In this study, a two-stage drilling procedure was used to investigate the effect of drilling on bone temperature. The procedure involved initially drilling to the desired depth with a smaller diameter drill bit, followed by drilling with a specified larger diameter drill bit. In the analytical model, predrilling was performed using a 2.5 mm diameter drill bit. This was followed by drilling with the specified bone drilling diameters of 3.5 mm or 4.5 mm to predict the effect of the second drilling stage on temperature increase. Figure 6a 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 6a, 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 predrilling), and the other involves making 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 5b 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 rise and mitigate the risk of thermal bone necrosis.

4.5. Effects of Rotational Speed on Thermally Affecting Zone (TAZ)

Figure 7a 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 7b 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 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 a 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. Some experimental studies using thermocouples to measure bone temperature, often placed near the drilled hole’s edge, 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 human bone could prevent heat from being generated during bone drilling. Our analytical findings from the model demonstrate comparable accuracy to the 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. [48] 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 was 48.4 °C, while at a lower rotation speed of 600 rpm, the temperature rose to as high as 152 °C. Mustafa et al. [49] 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 (20,000–100,000 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. [45] 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 as observed in the current literature.
When comparing the results in Figure 6a,b, it is found that the TAZ range, when drilling with a smaller drilling diameter, is smaller than when 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. Clinically, the affected zone of thermal damage is difficult to measure experimentally; the FE model provides a reliable approach to assess the potential risk region of thermal damage.
An ideal experiment for determining the bone temperature during drilling is difficult because real bone is a complex and anisotropic biological tissue. The development of 3D models of real machining operations in surgery, including anisotropic constitutive modelling, is one of the challenges because the accuracy of the simulated model depends greatly on the thermo-mechanical properties and the constitutive behaviour of the bone tissue. The anisotropic and heterogeneous structure of the cortical bone could be approximated as an isotropic homogeneous material due to the relatively low mechanical anisotropy of the cortical bone tissue [50].
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 [51], implementing step-power drilling [52], etc., there is no conclusive evidence to indicate that irrigation cooling been 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 [53,54,55,56]. 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 has 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.
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 level 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 interact with the second stage. This interaction could lead to an increase in drilling temperature during the second stage. However, further research or experiments are needed to confirm this effect.

6. Conclusions

This study utilised 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 minimise 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

Methodology, Writing—original draft preparation, Y.-C.C.; 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

  1. Bonfield, W.; Li, C.H. The temperature dependence of the deformation of bone. J. Biomech. 1968, 1, 323–329. [Google Scholar] [CrossRef] [PubMed]
  2. Thompson, H.C. Effect of drilling on bone. J. Oral Surg. 1958, 16, 22–30. [Google Scholar] [PubMed]
  3. Eriksson, R.A.; Albertson, T. The effect of heat on bone regeneration: An experimental study in the rabbit using the bone growth chamber. J. Oral Maxillofac. Surg. 1984, 42, 705–711. [Google Scholar] [CrossRef] [PubMed]
  4. 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–631. [Google Scholar] [CrossRef] [PubMed]
  5. Cordioli, G.; Majzoub, Z. Heat generation during Implant site preparation: An in vitro study. J. Oral Maxillofac. Implants 1997, 12, 186–193. [Google Scholar]
  6. Matthews, L.S.; Hirsch, C. Temperatures measured in human cortical bone when drilling. J. Bone Jt. Surg. 1972, 54-A, 297–308. [Google Scholar] [CrossRef]
  7. Brisman, D.L. The effect of speed, pressure, and time on bone temperature during the drilling of implant sites. Int. J. Oral Maxillofac. Implants 1996, 11, 35–37. [Google Scholar] [PubMed]
  8. Tehemar, S.H. Factors affecting heat generation during implant site preparations: A review of biologic observations and future considerations. Int. J. Oral Maxillofac. Implants 1999, 14, 127–136. [Google Scholar]
  9. Alevizakos, V.; Mosch, R.; von See, C. Influence of Multiple Used Implant Drills on Their Cutting Performance and Fracture Resistance. Materials 2023, 16, 5271. [Google Scholar] [CrossRef]
  10. 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]
  11. Hillery, M.T.; Shuaib, I. Temperature Effects in the Drilling of Human and Bovine Bone. J. Mater. Process. Technol. 1999, 92–93, 302–308. [Google Scholar] [CrossRef]
  12. Merchant, M.E. Mechanics of the Metal Cutting Process. J. Appl. Phys. 1945, 16, 267324. [Google Scholar]
  13. Lee, J.; Rabin, Y.; Ozdoganlar, O.B. A new thermal model for bone drilling wif applications to orthopaedic surgery. Med. Eng. Phys. 2011, 33, 1234–1244. [Google Scholar] [CrossRef] [PubMed]
  14. Maani, N.; Farhang, K.; Hodaei, M. A Model for the Prediction of Thermal Response of Bone in Surgical Drilling. J. Thermal Sci. Eng. Appl. 2014, 6, 041005. [Google Scholar] [CrossRef]
  15. Hu, Y.; Yan, Z.; Li, X.; Zhang, C.; Zheng, Q. Prediction model of bone drilling temperature based on heat source method in surgical rehabilitation. CIRP BoMainufacturing Conference 2019. Procedia CIRP 2020, 89, 263–269. [Google Scholar] [CrossRef]
  16. Davidson, S.R.; James, D.F. Drilling in Bone: Modeling Heat Generation and Temperature Distribution. J. Biomech. Eng. 2003, 125, 305–314. [Google Scholar] [CrossRef] [PubMed]
  17. 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]
  18. Sezek, S.; Aksakal, B.; Karaca, F. Influence of drill parameters on bone temperature and necrosis: A FEM modelling and in vitro experiments. Comput. Mater. Sci. 2012, 60, 13–18. [Google Scholar] [CrossRef]
  19. 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–1172. [Google Scholar] [CrossRef]
  20. 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]
  21. 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]
  22. 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]
  23. 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–1227. [Google Scholar] [CrossRef]
  24. 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]
  25. 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]
  26. 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–1260. [Google Scholar] [CrossRef] [PubMed]
  27. 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]
  28. 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]
  29. 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]
  30. 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]
  31. Liaoa, Z.; Axinte, D.A. On monitoring chip formation, penetration depth and cutting malfunctions in bone micro-drilling via acoustic emission. J. Mater. Process Technol. 2016, 229, 82–93. [Google Scholar] [CrossRef]
  32. Libonati, F.; Vergani, L. Understanding the structure property relationship in cortical bone to design a biomimetic composite. Compos. Struct. 2016, 139, 188–198. [Google Scholar] [CrossRef]
  33. 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] [PubMed]
  34. 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]
  35. 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]
  36. Khurshid, A.; Issam, M.B.; Naseer, A. Cortical bone drilling: An experimental and numerical study. Technol. Health Care 2014, 23, 1–12. [Google Scholar]
  37. 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]
  38. 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]
  39. 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–1957. [Google Scholar] [CrossRef]
  40. 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] [CrossRef]
  41. 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]
  42. Shen, Q.; Lee, T.C.; Lau, W.S. A Finite-Element Analysis of Temperature Distributions in Spade Drilling. J. Mater. Process. Technol. 1997, 66, 112–122. [Google Scholar] [CrossRef]
  43. Scott, F.M.; Albert, L.S. Thermo-mechanical finite element modelling of the friction drilling process. J. Manuf. Sci. Eng. 2007, 129, 531–538. [Google Scholar]
  44. Bonoa, M.; Ni, J. The location of the maximum temperature on the cutting edges of a drill. Int. J. Mach. Tools Manuf. 2006, 46, 901–907. [Google Scholar] [CrossRef]
  45. 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]
  46. Davidson, S.R.H.; James, D.F. Measurement of thermal conductivity of bovine cortical bone. Med. Eng. Phys. 2000, 22, 741–747. [Google Scholar] [CrossRef] [PubMed]
  47. 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]
  48. Inan, M.; Mizrak, B.; Ertem, K.; Harma, A.; Elmali, N.; Ayan, I. The factors affecting thermal necrosis secondary to the application of the Ilizarov transosseous wire. Acta Orthop. Traumatol. Turc. 2005, 39, 70–75. [Google Scholar]
  49. Mustafa, B.; David, F.J. Measurements of Shaft Speed While Drilling through Bone. J. Oral Maxillofac. Surg. 1995, 53, 1308–1315. [Google Scholar]
  50. Li, S.; Demirci, E.; Silberschmidt, V.V. Variability and anisotropy of mechanical behavior of cortical bone in tension and compression. J. Mech. Behav. Biomed. Mater. 2013, 21, 109–120. [Google Scholar] [CrossRef]
  51. Staroveski, T.; Brezak, D.; Udiljak, T. Drill wear monitoring in cortical bone drilling. Med. Eng. Phys. 2015, 37, 560–566. [Google Scholar] [CrossRef] [PubMed]
  52. 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, 36, 1449–1456. [Google Scholar] [CrossRef] [PubMed]
  53. 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–77. [Google Scholar] [CrossRef] [PubMed]
  54. 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]
  55. 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 2023, 10, 1168. [Google Scholar] [CrossRef]
  56. 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]
Figure 1. Geometry and mech models of the drill-bit and cortical bone.
Figure 1. Geometry and mech models of the drill-bit and cortical bone.
Bioengineering 11 00592 g001
Figure 2. (a) Experimental setup for bone temperature measurement and (b) schematic illustration of the embedded positions of thermocouples at a distance of 0.5 mm from the edge of the drilling hole with depths of 2 mm and 4 mm from the bone surface (x = 0.5 mm; z = 2 mm and 4 mm).
Figure 2. (a) Experimental setup for bone temperature measurement and (b) schematic illustration of the embedded positions of thermocouples at a distance of 0.5 mm from the edge of the drilling hole with depths of 2 mm and 4 mm from the bone surface (x = 0.5 mm; z = 2 mm and 4 mm).
Bioengineering 11 00592 g002
Figure 3. Comparisons of the bone temperature from the FE model and the experiments at the position x = 0.5 mm and z = 4 mm. The rotational speed is n = 1200 rpm and the bone density is d = 1.6 g/cc. The conditions are as follows: (a) feed force = 20 N, drill bit diameter D = 3.5 mm; (b) feed force = 60 N, drill bit diameter D = 4.5 mm. The upper and lower whiskers (⊤ and ⊥) display the experimental (green) maximum and minimum temperatures from 30 drilling boreholes, respectively.
Figure 3. Comparisons of the bone temperature from the FE model and the experiments at the position x = 0.5 mm and z = 4 mm. The rotational speed is n = 1200 rpm and the bone density is d = 1.6 g/cc. The conditions are as follows: (a) feed force = 20 N, drill bit diameter D = 3.5 mm; (b) feed force = 60 N, drill bit diameter D = 4.5 mm. The upper and lower whiskers (⊤ and ⊥) display the experimental (green) maximum and minimum temperatures from 30 drilling boreholes, respectively.
Bioengineering 11 00592 g003
Figure 4. Maximum temperature with radial distance at different depths (z-direction).
Figure 4. Maximum temperature with radial distance at different depths (z-direction).
Bioengineering 11 00592 g004
Figure 5. Effect of feeding forces on bone temperature: (a) temperature–drilling time relationship curve; (b) maximum temperature decreases exponentially with radial distance.
Figure 5. Effect of feeding forces on bone temperature: (a) temperature–drilling time relationship curve; (b) maximum temperature decreases exponentially with radial distance.
Bioengineering 11 00592 g005
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 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.
Bioengineering 11 00592 g006
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.
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.
Bioengineering 11 00592 g007
Table 1. The mechanical properties of the materials used in the FE model [39,40].
Table 1. The mechanical properties of the materials used in the FE model [39,40].
Material PropertyCortical BoneDrill Bit
Density (g/cm3)1.6407.990
Elastic modulus (MPa)16,700193,000
Poisson’s ratio0.30.25
Yielding stress (MPa)105290
Ultimate stress (MPa)106579
Ultimate strain0.0080.003
Specific heat (J/kg·°C)1640500
Thermal conductivity (W/m·k)0.45216.2
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.

Share and Cite

MDPI and ACS Style

Chen, Y.-C.; Tsai, Y.-J.; Hsiao, H.-Y.; Chiu, Y.-W.; Hong, Y.-Y.; Tu, Y.-K.; Hsiao, C.-K. Assessment of Thermal Osteonecrosis during Bone Drilling Using a Three-Dimensional Finite Element Model. Bioengineering 2024, 11, 592. https://doi.org/10.3390/bioengineering11060592

AMA Style

Chen Y-C, Tsai Y-J, Hsiao H-Y, Chiu Y-W, Hong Y-Y, Tu Y-K, Hsiao C-K. Assessment of Thermal Osteonecrosis during Bone Drilling Using a Three-Dimensional Finite Element Model. Bioengineering. 2024; 11(6):592. https://doi.org/10.3390/bioengineering11060592

Chicago/Turabian Style

Chen, Yung-Chuan, Yi-Jung Tsai, Hao-Yuan Hsiao, Yen-Wei Chiu, You-Yao Hong, Yuan-Kun Tu, and Chih-Kun Hsiao. 2024. "Assessment of Thermal Osteonecrosis during Bone Drilling Using a Three-Dimensional Finite Element Model" Bioengineering 11, no. 6: 592. https://doi.org/10.3390/bioengineering11060592

APA Style

Chen, Y. -C., Tsai, Y. -J., Hsiao, H. -Y., Chiu, Y. -W., Hong, Y. -Y., Tu, Y. -K., & Hsiao, C. -K. (2024). Assessment of Thermal Osteonecrosis during Bone Drilling Using a Three-Dimensional Finite Element Model. Bioengineering, 11(6), 592. https://doi.org/10.3390/bioengineering11060592

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop