Simulation of the mechanical interlocking capacity of a rough bone implant surface during healing
BioMedical Engineering OnLine volume 14, Article number: 45 (2015)
When an implant is inserted in the bone the healing process starts to osseointegrate the implant by creating new bone that interlocks with the implant. Biomechanical interlocking capacity is commonly evaluated in in vivo experiments. It would be beneficial to find a numerical method to evaluate the interlocking capacity of different surface structures with bone. In the present study, the theoretical interlocking capacity of three different surfaces after different healing times was evaluated by the means of explicit finite element analysis.
The surface topographies of the three surfaces were measured with interferometry and were used to construct a 3D bone-implant model. The implant was subjected to a displacement until failure of the bone-to-implant interface and the maximum force represents the interlocking capacity.
The simulated ratios (test/control) seem to agree with the in vivo ratios of Halldin et al. for longer healing times. However the absolute removal torque values are underestimated and do not reach the biomechanical performance found in the study by Halldin et al. which might be a result of unknown mechanical properties of the interface.
Finite element analysis is a promising method that might be used prior to an in vivo study to compare the load bearing capacity of the bone-to-implant interface of two surface topographies at longer healing times.
When an implant is inserted into bone the healing process starts by creating new bone that interlocks with the implant surface, so called osseointegration. It has been found that an implant with a rough surface topography enhances osseointegration and load bearing capacity [1, 2]. A significant amount of research has been undertaken in order to identify the biomechanical response to a structured surface, which is characterized by a set of statistical surface roughness parameters [3–7]. The general result of in vivo experiments is that increased implant surface roughness (Sa value) of cylindrical implants results in increased interfacial shear strength . Whether this empirical correlation is an effect of enhanced mechanical strength of the bone caused by the biological response to the surface and/or enhanced interlocking capacity is unclear. In addition, it is unclear whether surface roughness parameters truly reflect the load bearing capacity . Hansson and Norton  developed artificial surfaces characterized by various pit size, pit shape and pit density which were used to simulate the theoretical interfacial load bearing capacity. They found that the topographical structure affects the interfacial load bearing capacity . However, the theoretical surfaces considered by Hansson and Hansson  do not represent the surface topography of implants. Besides the implant surface topographical structure, the interfacial load bearing capacity depends on the mechanical properties of the bone during healing. In brief, bone formation in a gap proceeds, regardless of the presence of an implant, in two steps: (1) woven bone formation (2) remodeling of woven bone to lamellar/Harversian bone [11, 12]. Woven bone is characterized by a random arrangement of collagen fibers with poor mechanical properties . Woven bone is thereafter remodeled to Haversian bone over time by basic multicellular units (BMU) . The mineralization process of a new Haversian system can be divided into two stages; primary mineralization and secondary mineralization . Primary mineralization is characterized by a rapid, constant rate of mineralization that proceeds until 50–60% of the mineralization maximum has been reached [15, 16]. Following primary mineralization, a decreased rate of mineralization progressively continues and typically stabilizes at around 90–95% of the maximum level [15, 16]. The current knowledge of the mechanical properties of bone during healing is scarce. However, it is known that the degree of mineralization greatly affects the mechanical properties of bone [17–20]. It is evident that newly formed bone has different mechanical properties compared to mature bone due to different degrees of mineralization [21, 22]. Nano indentation has been used to investigate the mechanical properties of newly formed bone. Leong and Morgan  measured the mechanical properties of fractured rat bones after 24 days of healing and found that Young’s modulus of woven and cortical bone was 36.2 MPa and 7.2 GPa respectively. In a study by Ishimoto et al.  the mechanical properties of regenerated (2 weeks) and mature rabbit bone were measured. They obtained an elastic modulus of 17.0 and 27.6 GPa respectively which is higher than the maximum Young’s modulus for mature rabbit bone (~8 GPa) found by Isaksson et al. . Few studies have measured the mechanical properties of bone in the vicinity of implants. Chang et al.  obtained a value of Young’s modulus of pig alveolar bone close to an implant surface, after one month of healing, of 6.17 GPa that gradually increased to 7.9 GPa at a distance of 150 µm from the implant surface. Finally, it reached 10.1 GPa at a distance of 1,500 µm from the implant surface. This is slightly lower than what Vayron et al.  found in rabbit bone after longer healing times (7 and 13 weeks). An alternative way to determine the mechanical properties during healing was used in a study by Warzen et al. . They measured the stiffness of the implant bone interface in vivo of mice hind legs after 0–6 days of healing and by reverse engineering, obtained values of Young’s modulus ranging from 3 to 17 MPa. The mechanical properties of the bone and the surface topography affect the bone implant interfacial shear strength . The objective of the present study was to estimate the theoretical bone-to-implant interfacial shear strength for different surface topographies during healing by means of finite element analysis (FEA).
Geometrical representation of a surface topography
The coordinates (x, y and z) from a representative patch of each of the three different (A, B and C) surfaces used in the study by Halldin et al.  were extracted from interferometry measurements. Group A were blasted with coarse titanium particles thereafter treated in oxalic acid and hydrofluoric acid sequentially (CB-AT-I). Group B were blasted with fine particles titanium oxide particles and thereafter treated in oxalic acid and hydrofluoric acid sequentially (FB-AT-I). Group C implants were blasted with coarse titanium particles thereafter treated in hydrofluoric acid (CB-HF). The representative patches were selected to have surface average height values (Sa) close to the mean Sa values of the implant in Halldin et al. . The x, y and z coordinates were imported into MATLAB (MATLAB 2013b, Mathworks Inc., USA) to calculate the average height values (Sa [μm]), root-mean-square (Sq [μm]), skewness (Ssk [μm/μm]) and the root-mean-square of the surface slope (Sdq [μm/μm]) of the selected patches . Thereafter, the surfaces were exported to STL file format followed by conversion with high accuracy, to IGES format using Geomagic studio® (Geomagic Solutions, Cary, NC, USA). The IGES surfaces were imported into Pro Engineer (PTC, Needham, MA, USA) for creation of a 3D geometry (Figure 1). Finally the geometries were imported to ANSYS® 14.5 (ANSYS INC, Canonsburg, PA, USA) for simulation of the interfacial load bearing capacity.
Finite element model
The bone-to-implant interface was modeled according to Figure 2. The implant surface was in contact with interfacial bone which in turn was adjacent to the fictitious bone. The fictitious bone represents the structural stiffness of the surrounding bone and was modeled as an isotropic material with a thickness (l 0) of 10 µm and a Young’s modulus represented by Es. Substituting the surrounding bone by a fictitious bone results in a less computationally heavy model. The implant was locked in y direction and was moved in x direction by a constant acceleration (ax) during the simulation time. The interface between interfacial bone and fictitious bone was fixed in x direction and free to move in y direction. The outer sides of the X–Y plane of the 3D model were restricted to move in z direction. The implant was modeled as a rigid component and the interface bone-to-implant interface was modeled as a frictional contact with an assumed coefficient of friction of 0.4 .
The Young’s modulus (Es) of the fictitious bone was determined by simulating the deformation (d) at a predefined internal pressure (P) of 100 MPa in a 2D plane strain model representing a disc of mature bone. The disc had an inner diameter of 3.5 mm similar to the outer diameter of the implant in the study by Halldin et al. . The structural stiffness was computed for different outer diameters (Do) of the disc ranging from 2.5 to 12.5 mm. The fictitious Young’s modulus Es, for each outer diameter (Do) with corresponding deformation (d), was calculated using Eq. (1).
Since the current knowledge of the mechanical properties of bone in the vicinity of an implant during healing is limited, the assumed mechanical properties of bone during healing were derived from the correlation between mechanical properties and mineral content  and from the correlation between mineral content and healing time . Currey  investigated the mineral content, expressed as milligrams of calcium (Ca) per gram dried defatted bone, for various species with corresponding Young’s modulus (E), yield strain (εy) and post yield strain (εpy). Using Currey´s data the correlation between Ca and E, εy and εpy can be identified (Figure 3a). Furthermore, by assuming 100% mineralization after 350 days of healing , and a corresponding Young’s modulus of 7,950 MPa , the Young’s modulus (E), yield strain (εy) and ultimate strain (εu = εy + εpy) at different healing times can be approximated (Figure 3b). In the present study it was assumed that the implant interfacial bone exhibits a bilinear isotropic material behavior without hardening. Thus the stress strain curves and the mechanical properties at different healing times can be estimated (Figure 4; Table 1). The strength and Young’s modulus of titanium are an order of magnitude higher than those of the bone. Therefore, to reduce computing time the implant was assumed to be rigid. In the present study the shear strength of a bone-to-implant interface was simulated with explicit finite element analysis (Figure 2) using the ultimate strain as a failure criterion. According to the simulations the maximum shear force occurred at implant displacements of less than 0.75 µm. Therefore, the simulations were discontinued after an implant displacement of 0.75 µm was reached. The displacement (0.75 µm) was achieved by a constant acceleration of the implant during the simulation time. Using the representative patch area the maximum reaction force between implant interfacial bone and the fictitious bone was converted, to interfacial shear strength (τp). In explicit FE simulation the computing time depends on the choice of simulation time, minimum element size and material density . Reduced computing time is obtained by decreased simulation time (increased acceleration), increased material density and increased minimum element size . Prior to the 3D simulation a parameter study (convergence test) with several combinations of different values of the acceleration (P1), material density increase (P2) and minimum element size (P3) (Table 2) was performed to obtain reasonable computing time without compromising the accuracy. The parameter study was performed by simulating the reaction force on a 2D profile. The 2D profile was extracted from an x–y plane of the 3D model (Figure 2). The meshes were auto generated in ANSYS® with different values of the minimum element size at the bone-to-implant interface. Otherwise the simulations were performed with the same settings. The simulated shear strength for the three surfaces were converted to removal torque value (RTQ) of an implant with the same design as in the Halldin et al.  study according to Eq. (2).
where A i represents the outer threaded area of the implant which was 78 mm2, r represents the radius to the implant surface, which was 1.75 mm, and BIC represents the bone to implant contact obtained in the in vivo study , which was set to 40% for all three surfaces.
The surface roughness (Sa value) of the selected representative patch (Figure 1) and the Sa value with the standard deviation (SD) of the implant in the study by Halldin et al.  are presented in Table 3. The fictitious Young’s modulus Es of the bone surrounding a 3.5 mm (Di) implant with different values of the outer diameter (D0) is presented in Figure 5. The fictitious Young’s modulus was chosen to 33 MPa. A reasonable computing time was obtained with an acceleration of 600 m/s2 (resulting in a simulation time of 5 × 10−5 s), material density increase by a factor of 1,000 and a minimum element size of 1 µm. The results of the 2D analysis indicate that these settings do not seem to significantly affect the reaction force compared to simulation results using lower parameters values (P1, P2 and P3) (Figure 6). Hence these settings were used in the 3D simulations which resulted in a computing time of 24 h for each simulation. The 3D models were meshed with ANSYS® default settings and a minimum element size of 1 µm resulting in meshes according to Table 4 and Figure 7. The interfacial shear strength of the three surfaces is presented in Figure 8a. The theoretical removal torques, according to Eq. (1), and the corresponding in vivo removal torques of Halldin et al.  are presented in Figure 8b. The ratio (test/control) of the simulated shear strength and the mean of the corresponding ratios of the in vivo study by Halldin et al.  for the surfaces are presented in Figure 8c.
In the present study the theoretical interlocking capacity of three different implant surfaces with bone after different healing times was evaluated. The simulated ratios (test/control) seem to agree with those of the Halldin et al.  study for longer healing times (Figure 8). However, the simulated absolute removal torque is underestimated compared to the biological removal torque as found in Halldin et al.  (Figure 8b). There are many factors that may influence the simulated results: (1) model setup (2) representative mechanical properties of the interfacial bone during healing (3) selection of the representative patch (4) pressure at the bone-to-implant interface (5) bone-to-implant contact. In simulation it is essential to use a model that captures the essence of what is intended to be analyzed. However, a model will always be a simplification of reality. In this study a simplified model of the interface was developed to analyze the interlocking capacity of a rough surface with bone.
In the present 3D-model the circumferential structural displacement was neglected. The circumferential structural displacement would affect the angle of rotation but is not assumed to affect the interfacial shear strength. This simplification reduces the size of the model and thus the computing time. The 2D analyses were performed to investigate how the reaction force is affected by different simulation settings. A reduced element size was assumed to provide a higher accuracy but would have implied a longer computing time . The results of the 2D analysis showed that a minimum element size of 1 µm could be used without negatively affecting the accuracy of the results. Furthermore to reduce computing time the simulation time was decreased by the use of increased acceleration. An increased simulated acceleration might affect the interfacial shear force. To reduce the inertia effects the simulation was performed with a constant acceleration of 600 m/s2, which ramped up the velocity from 0 to 30 mm/s during 5 × 10−5 s. To reduce the computing time even further the material density was increased in turn affecting the inertia and consequently the results. In the 2D analysis it was found that the material density could be increased by a factor of 1,000 without substantially affecting the interfacial shear strength.
The magnitude of the interfacial shear strength is affected by the mechanical properties of interfacial bone as well as the interface properties (i.e. friction) To the authors’ knowledge the coefficient of friction on micro level is unknown and the assumed coefficient of friction was based on a coefficient of friction between bone and titanium . The coefficient of friction is assumed to affect the magnitude of the interfacial shear strength to the same extent between the simulations.
The knowledge about the mechanical properties of bone, including failure behavior, in the vicinity of an implant surface during healing is limited. Therefore, an estimation of the mechanical properties of bone was derived from the mineral content measurements of osteoid under mineralization in an osteon during healing by Fuchs  and the relationship between mineral content and material properties found by Currey . Furthermore, it was assumed that the mineral content was 100% after 350 days of healing which corresponds to a Young’s modulus of 7.95 GPa . To estimate the degree of mineralization during the first weeks of healing, extrapolation was made of the experimental data obtained by Currey . Thus in the present study, Young’s modulus during healing was set in the range of 2–10 GPa. Other studies have found a variation in Young’s modulus during healing which could be caused by differences in experimental setup and choice of animals used [13, 23–25, 30]. The assumed mechanical properties used in this study might not reflect the true material properties during healing which in turn affect the absolute value of the removal torque. However, when the relative interfacial load bearing capacity of two surfaces is considered the actual material property during healing is assumed to have limited impact. Furthermore, the same bone material properties were used for all surfaces at coinciding healing times and a potential biological change in the mechanical properties caused by a specific surface topography was not considered. This might be an explanation of the consistency in ratios for longer healing times and the deviations during the early healing phase. This indicates that, for longer healing times the surface topography might be the main parameter influencing the interlocking capacity of an osseointegrated implant surface.
In this simulation the ultimate strain was used as a failure criterion based on available ultimate strain data. Principal strain failure can be used to represent brittle or ductile failure in materials . Other material failure criteria, such as Mohr–Coulomb, Tsai–Wu and Hill, have been discussed for bone . However to our knowledge the anisotropy and a suitable yield and failure criterion for the bone material during healing in the vicinity of the implant surface are unknown. When the failure criterion was reached in an element the element was removed and did no longer contribute to the interfacial shear strength. This might result in an underestimation of the simulated interfacial shear strength value. The magnitude of the reaction force depends on the surface topography of the selected patch. The standard deviation of the Sa value (Table 3) of the implants in the Halldin et al.  study, from which the representative patches were selected, indicates a variation. Even though the Sa value of the representative patch was selected to be in the range of the Sa value of the implant, the selected patch might not fully represent the surface topography with respect to load bearing capacity. When an implant is inserted it induces static strains in the bone that gradually decrease during healing [32, 33]. These initial residual stresses increase the removal torque value during the early healing period. The current simulations did not include the change of residual stresses during healing time. Therefore, the absolute removal torque value might be underestimated during the early healing phase. During osseintegration the biological processes start to integrate the implant with bone and thus increase the bone-to-implant contact and the removal torque over time. In this simulation, 40% bone-to-implant contact was assumed, for coinciding healing times, which might lead to an over- or underestimated simulated absolute removal torque, but is assumed to have a limited impact when analyzing the ratios.
In this simulation the theoretical interfacial shear strength of the three surface topographies was simulated and compared to the results in the study by Halldin et al. . The simulated ratios (test/control) seem to agree with those of the Halldin et al.  for longer healing times (Figure 8). Despite differences in the absolute removal torque values between simulations and the in vivo data, finite element analysis is a promising method, that can be used prior to in vivo study, to compare the theoretical load bearing capacities of the bone-to-implant interfaces especially for longer healing times. It was also concluded that the mechanical properties and the material model of bone in the vicinity of an implant affect the theoretical load bearing capacity and requires more research.
bone to implant contact
Albrektsson T, Wennerberg A. Oral implant surfaces: part 1–review focusing on topographic and chemical properties of different surfaces and in vivo responses to them. Int J Prosthodont. 2004;17(5):536–43.
Wennerberg A, Albrektsson T. Effects of titanium surface topography on bone integration: a systematic review. Clin Oral Implants Res. 2009;20(Suppl 4):172–84.
Johansson CB, Gretzer C, Jimbo R, Mattisson I, Ahlberg E. Enhanced implant integration with hierarchically structured implants: a pilot study in rabbits. Clin Oral Implants Res. 2012;23(8):943–53.
Buser D, Nydegger T, Oxland T, Cochran DL, Schenk RK, Hirt HP, et al. Interface shear strength of titanium implants with a sandblasted and acid-etched surface: a biomechanical study in the maxilla of miniature pigs. J Biomed Mater Res. 1999;45(2):75–83.
Wennerberg A, Ektessabi A, Albrektsson T, Johansson C, Andersson B. A 1-year follow-up of implants of differing surface roughness placed in rabbit bone. Int J Oral Maxillofac Implants. 1997;12(4):486–94.
Wennerberg A, Albrektsson T, Lausmaa J. Torque and histomorphometric evaluation of c.p. titanium screws blasted with 25- and 75-microns-sized particles of Al2O3. J Biomed Mater Res. 1996;30(2):251–60.
Wennerberg A, Albrektsson T, Johansson C, Andersson B. Experimental study of turned and grit-blasted screw-shaped implants with special emphasis on effects of blasting material and surface topography. Biomaterials. 1996;17(1):15–22.
Loberg JMI, Hansson S, Ahlberg E. Characterisation of titanium dental implants I: critical assessment of surface roughness parameters. TOBIOMTJ. 2010;2:1–8.
Hansson S, Hansson KN. The effect of limited lateral resolution in the measurement of implant surface roughness: a computer simulation. J Biomed Mater Res A. 2005;75(2):472–7.
Hansson S, Norton M. The relation between surface roughness and interfacial shear strength for bone-anchored implants. A mathematical model. J Biomech. 1999;32(8):829–36.
Brunski JB. In vivo bone response to biomechanical loading at the bone/dental-implant interface. Adv Dent Res. 1999;13:99–119.
Colnot C, Romero DM, Huang S, Rahman J, Currey JA, Nanci A, Brunski JB, Helms JA. Molecular analysis of healing at a bone-implant interface. J Dent Res. 2007;86(9):862–7.
Leong PL, Morgan EF. Correlations between indentation modulus and mineral density in bone-fracture calluses. Integr Comp Biol. 2009;49(1):59–68.
Parfitt AM. Osteonal and hemi-osteonal remodeling: the spatial and temporal framework for signal traffic in adult human bone. J Cell Biochem. 1994;55(3):273–86.
Fuchs RK, Allen MR, Ruppel ME, Diab T, Phipps RJ, Miller LM, et al. In situ examination of the time-course for secondary mineralization of Haversian bone using synchrotron Fourier transform infrared microspectroscopy. Matrix Biol. 2008;27(1):34–41.
Davison KS, Siminoski K, Adachi JD, Hanley DA, Goltzman D, Hodsman AB, et al. Bone strength: the whole is greater than the sum of its parts. Semin Arthritis Rheum. 2006;36(1):22–31.
Currey JD. Tensile yield in compact bone is determined by strain, post-yield behaviour by mineral content. J Biomech. 2004;37(4):549–56.
Currey JD. The effect of porosity and mineral content on the Young’s modulus of elasticity of compact bone. J Biomech. 1988;21(2):131–9.
Currey JD. Strain rate and mineral content in fracture models of bone. J Orthop Res. 1988;6(1):32–8.
Spatz HC, O’Leary EJ, Vincent JF. Young’s moduli and shear moduli in cortical bone. Proc Biol Sci. 1996;263(1368):287–94.
Isaksson H, Harjula T, Koistinen A, Iivarinen J, Seppanen K, Arokoski JP, et al. Collagen and mineral deposition in rabbit cortical bone during maturation and growth: effects on tissue properties. J Orthop Res. 2010;28:1626–33.
Isaksson H, Malkiewicz M, Nowak R, Helminen HJ, Jurvelin JS. Rabbit cortical bone tissue increases its elastic stiffness but becomes less viscoelastic with age. Bone. 2010;47(6):1030–8.
Ishimoto T, Nakano T, Yamamoto M, Tabata Y. Biomechanical evaluation of regenerating long bone by nanoindentation. J Mater Sci Mater M. 2011;22(4):969–76.
Chang MC, Ko CC, Liu CC, Douglas WH, DeLong R, Seong WJ, et al. Elasticity of alveolar bone near dental implant-bone interfaces after one month’s healing. J Biomech. 2003;36(8):1209–14.
Vayron R, Matsukawa M, Tsubota R, Mathieu V, Barthel E, Haiat G. Evolution of bone biomechanical properties at the micrometer scale around titanium implant as a function of healing time. Phys Med Biol. 2014;59(6):1389–406.
Wazen RM, Currey JA, Guo H, Brunski JB, Helms JA, Nanci A. Micromotion-induced strain fields influence early stages of repair at bone-implant interfaces. Acta Biomater. 2013;9(5):6663–74.
Halldin A, Jimbo R, Johansson CB, Gretzer C, Jacobsson M. Improved osseointegration and interlocking capacity with dual acid-treated implants: a rabbit study. Clin Oral Implants Res. 2014. doi:10.1111/clr.12507.
Mischler S, Pax G. Tribological behavior of titanium sliding against bone. Eur Cell Mater Suppl. 2002;1(3):28–9.
ANSYS: Manual ANSYS software 14.5. Canonsburg.
Clark PA, Clark AM, Rodriguez A, Hussain MA, Mao JJ. Nanoscale characterization of bone–implant interface and biomechanical modulation of bone ingrowth. Mat Sci Eng C. 2007;27(3):382–93.
Garcia D. Elastic plastic damage laws for cortical bone. Thesis EPFL, no 3435, LAUSANNE 2006.
Halldin A, Jimbo R, Johansson CB, Wennerberg A, Jacobsson M, Albrektsson T, Hansson S. The effect of static bone strain on implant stability and bone remodeling. Bone. 2011;49(4):783–9.
Halldin A, Jimbo R, Johansson CB, Wennerberg A, Jacobsson M, Albrektsson T, et al. Implant stability and bone remodeling after 3 and 13 days of implantation with an initial static strain. Clin Implant Dent Relat Res. 2014;16(3):383–93.
AH developed appropriate material model performed finite element simulations and drafted the article. MA participated in the design and development of the model and drafting the article. SH performed the research of the mechanical behavior of bone and developed appropriate material model and drafting the article. MJ contributed to develop biological interfacial model and drafting the article and contributed for intellectual content. All authors read and approved the final manuscript.
This work was supported by the Swedish Research Council (621-2010-4760). We would like to thank Klas Johansson and Mikael Lauth at EDR&Medeso for the valuable support during simulation.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
About this article
Cite this article
Halldin, A., Ander, M., Jacobsson, M. et al. Simulation of the mechanical interlocking capacity of a rough bone implant surface during healing. BioMed Eng OnLine 14, 45 (2015). https://doi.org/10.1186/s12938-015-0038-0
- Bone-to-implant interface shear strength
- Bone healing