Debonding of coin-shaped osseointegrated implants: Coupling of experimental and numerical approaches

While cementless implants are now widely used clinically, implant debonding still occur and is difficult to anticipate. Assessing the biomechanical strength of the bone-implant interface can help improving the understanding of osseointegration phenomena and thus preventing surgical failures. A dedicated and standardized implant model was considered. The samples were tested using a mode III cleavage device to assess the mechanical strength of the bone-implant interface by combining experimental and numerical approaches. Four rough (Sa = 24.5 μm) osseointegrated coin-shaped implants were left in sheep cortical bone during 15 weeks of healing time. Each sample was experimentally rotated at 0.03°/sec until complete rupture of the interface. The maximum values of the torque were comprised between 0.48 and 0.72 N m, while a significant increase of the normal force from 7-12 N to 31-43 N was observed during the bone-implant interface debonding, suggesting the generation of bone debris at the bone-implant interface. The experimental results were compared to an isogeometric finite element model describing the adhesion and debonding phenomena through a modified Coulomb's law, based on a varying friction coefficient to represent the transition from an unbroken to a broken bone-implant interface. A good agreement was found between numerical and experimental torques, with numerical friction coefficients decreasing from 8.93 to 1.23 during the bone-implant interface rupture, which constitutes a validation of this model to simulate the debonding of an osseointegrated bone-implant interface subjected to torsion.

Using an Instrumented Hammer to Predict the Rupture of Bone Samples Subject to an Osteotomy

Osteotomies are common procedures in maxillofacial and orthopedic surgery. The surgeons still rely on their proprioception to control the progression of the osteotome. Our group has developed an instrumented hammer that was shown to provide information on the biomechanical properties of the tissue located around the osteotome tip. The objective of this study is to determine if this approach may be used to predict the rupture of a bone sample thanks to an instrumented hammer equipped with a force sensor. For each impact, an indicator τ is extracted from the signal corresponding to the variation of the force as a function of time. A linear by part regression analysis is applied to the curve corresponding to the variation of τ as a function of the distance d between the tip of the osteotome and the end of the sample. The experiments were conducted with plywood and bovine trabecular bone samples. The results show that τ starts increasing when the value of d is lower than 2.6 mm on average, which therefore corresponds to a typical threshold detection distance between the osteotome tip and the sample end. These findings open new paths for the development of this instrumented surgical hammer.

Anatomical subject validation of an instrumented hammer using machine learning for the classification of osteotomy fracture in rhinoplasty

Osteotomies during rhinoplasty are usually based on the surgeon's proprioception to determine the number and the strength of the impacts. The aim of this study is to determine whether a hammer instrumented with a force sensor can be used to classify fractures and to determine the location of the osteotome tip. Two lateral osteotomies were realized in nine anatomical subjects using an instrumented hammer recording the evolution of the impact force. Two indicators τ and λ were derived from the signal, and video analysis was used to determine whether the osteotome tip was located in nasal or frontal bone as well as the condition of the bone tissue around the osteotome tip. A machine-learning algorithm was used to predict the condition of bone tissue after each impact. The algorithm was able to predict the condition of the bone after the impacts with an accuracy of 83%, 91%, and 93% when considering a tolerance of 0, 1, and 2 impacts, respectively. Moreover, in nasal bone, the values of τ and λ were significantly lower (p < 10^-10) and higher (p < 10^-4) than in frontal bone, respectively. This study paves the way for the development of the instrumented hammer as a decision support system.

Validation of an Instrumented Hammer for Rhinoplasty Osteotomies: A Cadaveric Study

Background: Osteotomies during rhinoplasty are usually based on surgeon's proprioception to determine the number, energy, and trajectory of impacts. Objective: The first objective was to detect the occurrence of fractures. The second objective was to determine when the thicker frontal bone was encountered by the osteotome. Materials and Methods: An instrumented hammer was used to measure the impact force during lateral osteotomies on nine human anatomic specimens. A prediction algorithm was developed using machine learning techniques, to detect the occurrence of fractures, and the proximity of the osteotome to the frontal bone. Results: The algorithm was able to predict the occurrence of fractures and the proximity to the frontal bone with a prediction rate of 83%, 91%, and 93% when allowing for an error of 0, 1, and 2 impacts, respectively. The location of the osteotome in the frontal bone was predicted with an error of 7.7%. Conclusion: An osteotomy hammer measuring the impact force when performing lateral osteotomies can predict the occurrence of fractures and the proximity to the frontal bone, providing the surgeon with instant feedback.

Assessment of dental implant stability using resonance frequency analysis and quantitative ultrasound methods

Purpose Quantitative ultrasound (QUS) and resonance frequency analyses (RFA) are promising methods to assess the stability of dental implants. The aim of this in vivo preclinical study is to compare the results obtained with these two techniques with the bone-implant contact (BIC) ratio, which is the gold standard to assess dental implant stability.Methods Twenty-two identical dental implants were inserted in the tibia and femur of 12 rabbits, which were sacrificed after different healing durations (0, 4, 8 and 13 weeks). For each implant, the ultrasonic indicator (UI) and the implant stability quotient (ISQ) were retrieved just before the animal sacrifice using the QUS and RFA techniques, respectively. Histomorphometric analyses were carried out to estimate the bone-implant contact ratio.Results UI values were found to be better correlated to BIC values (R²=0.47) compared to ISQ values (R²=0.39 for measurements in one direction and R²=0.18 for the other direction), which were shown to be dependent on the direction of measurements. Errors realized on the UI were around 3.3 times lower to the ones realized on the ISQ.Conclusions QUS provide a better estimation of dental implant stability compared to RFA. This study paves the way for the future clinical development of a medical device aiming at assessing dental implant stability in a patient-specific manner. Clinical studies should confirm these results in the future.

Ultrasonic assessment of osseointegration phenomena at the bone-implant interface using convolutional neural network

Although endosseous implants are widely used in the clinic, failures still occur and their clinical performance depends on the quality of osseointegration phenomena at the bone-implant interface (BII), which are given by bone ingrowth around the BII. The difficulties in ensuring clinical reliability come from the complex nature of this interphase related to the implant surface roughness and the presence of a soft tissue layer (non-mineralized bone tissue) at the BII. The aim of the present study is to develop a method to assess the soft tissue thickness at the BII based on the analysis of its ultrasonic response using a simulation based-convolution neural network (CNN). A large-annotated dataset was constructed using a two-dimensional finite element model in the frequency domain considering a sinusoidal description of the BII. The proposed network was trained by the synthesized ultrasound responses and was validated by a separate dataset from the training process. The linear correlation between actual and estimated soft tissue thickness shows excellent R² values equal to 99.52% and 99.65% and a narrow limit of agreement corresponding to [ -2.56, 4.32 μm] and [ -15.75, 30.35 μm] of microscopic and macroscopic roughness, respectively, supporting the reliability of the proposed assessment of osseointegration phenomena.

Neutron microtomography to investigate the bone-implant interface-comparison with histological analysis

Bone properties and especially its microstructure around implants are crucial to evaluate the osseointegration of prostheses in orthopaedic, maxillofacial and dental surgeries. Given the intrinsic heterogeneous nature of the bone microstructure, an ideal probing tool to understand and quantify bone formation must be spatially resolved. X-ray imaging has often been employed, but is limited in the presence of metallic implants, where severe artifacts generally arise from the high attenuation of metals to x-rays. Neutron tomography has recently been proposed as a promising technique to study bone-implant interfaces, thanks to its lower interaction with metals. The aim of this study is to assess the potential of neutron tomography for the characterisation of bone tissue in the vicinity of a metallic implant. A standardised implant with a bone chamber was implanted in rabbit bone. Four specimens were imaged with neutron tomography and subsequently compared to non-decalcified histology to stain soft and mineralised bone tissues, used here as a ground-truth reference. An intensity-based image registration procedure was performed to place the 12 histological slices within the corresponding 3D neutron volume. Significant correlations (p < 0.01) were obtained between the two modalities for the bone-implant contact (BIC) ratio (R = 0.77) and the bone content inside the chamber (R = 0.89). The results indicate that mineralised bone tissue can be reliably detected by neutron tomography. However, theBICratio and bone content were found to be overestimated with neutron imaging, which may be explained by its sensitivity to non-mineralised soft tissues, as revealed by histological staining. This study highlights the suitability of neutron tomography for the analysis of the bone-implant interface. Future work will focus on further distinguishing soft tissues from bone tissue, which could be aided by the adoption of contrast agents.

Using an impact hammer to perform biomechanical measurements during osteotomies: Study of an animal model

Osteotomies are common surgical procedures used for instance in rhinoplasty and usually performed using an osteotome impacted by a mallet. Visual control being difficult, osteotomies are often based on the surgeon proprioception to determine the number and energy of each impact. The aim of this study is to determine whether a hammer instrumented with a piezoelectric force sensor can be used to (i) follow the displacement of the osteotome and (ii) determine when the tip of the osteotome arrives in frontal bone, which corresponds to the end of the osteotomy pathway. Seven New Zealand White rabbit heads were collected, and two osteotomies were performed on their left and right nasal bones using the instrumented hammer to record the variation of the force as a function of time during each impact. The second peak time τ was derived from each signal while the displacement of the osteotome tip D was determined using video motion tracking. The results showed a significant correlation between τ and D (ρ² = 0.74), allowing to estimate the displacement of the osteotome through the measurement of τ. The values of τ measured in the frontal bone were significantly lower than in the nasal bone (p<10^-10), which allows to determine the transition between the nasal and frontal bones when τ becomes lower than 0.78 its initial averaged value. Although results should be validated clinically, this technology could be used by surgeons in the future as a decision support system to help assessing the osteotome environment.

Quantitative ultrasound assessment of the influence of roughness and healing time on osseointegration phenomena

The evolution of bone tissue quantity and quality in contact with the surface of orthopedic and dental implants is a strong determinant of the surgical outcome but remains difficult to be assessed quantitatively. The aim of this study was to investigate the performance of a quantitative ultrasound (QUS) method to measure bone-implant interface (BII) properties. A dedicated animal model considering coin-shaped titanium implants with two levels of surface roughness (smooth, S_a = 0.49 µm and rough, S_a = 3.5 µm) allowed to work with a reproducible geometry and a planar interface. The implants were inserted in rabbit femurs and tibiae for 7 or 13 weeks. The ultrasonic response of the BII was measured ex vivo, leading to the determination of the 2-D spatial variations of bone in contact with the implant surface. Histological analysis was carried out to determine the bone-implant contact (BIC) ratio. The amplitude of the echo was significantly higher after 7 weeks of healing time compared to 13 weeks, for both smooth (p < 0.01) and rough (p < 0.05) implants. A negative correlation (R = − 0.63) was obtained between the ultrasonic response and the BIC. This QUS technique is more sensitive to changes of BII morphology compared to histological analyses.

Analytical modeling of the interaction of an ultrasonic wave with a rough bone-implant interface

Quantitative ultrasound can be used to characterize the evolution of the bone-implant interface (BII), which is a complex system due to the implant surface roughness and to partial contact between bone and the implant. The determination of the constitutive law of the BII would be of interest in the context of implant acoustical modeling in order to take into account the imperfect characteristics of the BII. The aim of the present study is to propose an analytical effective model describing the interaction between an ultrasonic wave and a rough BII. To do so, a spring model was considered to determine the equivalent stiffness K of the BII. The stiffness contributions related (i) to the partial contact between the bone and the implant and (ii) to the presence of soft tissues at the BII during the process of osseointegration were assessed independently. K was found to be comprised between 10¹³ and 10¹⁷ N/m³ depending on the roughness and osseointegration of the BII. Analytical values of the reflection and transmission coefficients at the BII were derived from values of K. A good agreement with numerical results obtained through finite element simulation was obtained. This model may be used for future finite element bone-implant models to replace the BII conditions.

Multimodal Evaluation of the Spatiotemporal Variations of Periprosthetic Bone Properties

Titanium implants are widely used in dental and orthopedic surgeries. However, implant failures still occur because of a lack of implant stability. The biomechanical properties of bone tissue located around the implant need to be assessed to better understand the osseointegration phenomena and anticipate implant failure. The aim of this study was to explore the spatiotemporal variation of the microscopic elastic properties of newly formed bone tissue close to an implant. Eight coin-shaped Ti6Al4V implants were inserted into rabbit tibiae for 7 and 13 weeks using an in vivo model allowing the distinction between mature and newly formed bone in a standardized configuration. Nanoindentation and micro-Brillouin scattering measurements were carried out in similar locations to measure the indentation modulus and the wave velocity, from which relative variations of bone mass density were extracted. The indentation modulus, the wave velocity and mass density were found to be higher (1) in newly formed bone tissue located close to the implant surface, compared to mature cortical bone tissue, and (2) after longer healing time, consistently with an increased mineralization. Within the bone chamber, the spatial distribution of elastic properties was more heterogeneous for shorter healing durations. After 7 weeks of healing, bone tissue in the bone chamber close to the implant surface was 12.3% denser than bone tissue further away. Bone tissue close to the chamber edge was 16.8% denser than in its center. These results suggest a bone spreading pathway along tissue maturation, which is confirmed by histology and consistent with contact osteogenesis phenomena.

Spatio-temporal evolution of hydroxyapatite crystal thickness at the bone-implant interface

A better understanding of bone nanostructure around the bone-implant interface is essential to improve longevity of clinical implants and decrease failure risks. This study investigates the spatio-temporal evolution of mineral crystal thickness and plate orientation in newly formed bone around the surface of a metallic implant. Standardized coin-shaped titanium implants designed with a bone chamber were inserted into rabbit tibiae for 7 and 13 weeks. Scanning measurements with micro-focused small-angle X-ray scattering (SAXS) were carried out on newly formed bone close to the implant and in control mature cortical bone. Mineral crystals were thinner close to the implant (1.8 ± 0.45 nm at 7 weeks and 2.4 ± 0.57 nm at 13 weeks) than in the control mature bone tissue (2.5 ± 0.21 nm at 7 weeks and 2.8 ± 0.35 nm at 13 weeks), with increasing thickness over healing time (+30 % in 6 weeks). These results are explained by younger bone close to the implant, which matures during osseointegration. Thinner mineral crystals parallel to the implant surface within the first 100 µm indicate that the implant affects the ultrastructure of neighbouring bone , potentially due to heterogeneous interfacial stresses, and suggest a longer maturation process of bone tissue and difficulty in binding to the metal. The bone growth kinetics within the bone chamber was derived from the spatio-temporal evolution of bone tissue's nanostructure, coupled with microtomographic imaging. The findings indicate that understanding mineral crystal thickness or plate orientation can improve our knowledge of osseointegration.

Ultrasonic Propagation in a Dental Implant

Ultrasound techniques can be used to characterize and stimulate dental implant osseointegration. However, the interaction between an ultrasonic wave and the implant-bone interface (IBI) remains unclear. This study-combining experimental and numerical approaches-investigates the propagation of an ultrasonic wave in a dental implant by assessing the amplitude of the displacements along the implant axis. An ultrasonic transducer was excited in a transient regime at 10 MHz. Laser interferometric techniques were employed to measure the amplitude of the displacements, which varied 3.2-8.9 nm along the implant axis. The results demonstrated the propagation of a guided wave mode along the implant axis. The velocity of the first arriving signal was equal to 2110 m.s^-1, with frequency components lower than 1 MHz, in agreement with numerical results. Investigating guided wave propagation in dental implants should contribute to improved methods for the characterization and stimulation of the IBI.

Reflection of an ultrasonic wave on the bone-implant interface: Comparison of two-dimensional and three-dimensional numerical models

Quantitative ultrasound is used to characterize osseointegration at the bone-implant interface (BII). However, the interaction between an ultrasonic wave and the implant remains poorly understood. Hériveaux, Nguyen, and Haiat [(2018). J. Acoust. Soc. Am. 144, 488-499] recently employed a two-dimensional (2D) model of a rough BII to investigate the sensitivity of the ultrasonic response to osseointegration. The present letter aimed at assessing the validity of the 2D assumption. The values of the reflection coefficient of the BII obtained with two and three-dimensional models were found not to be significantly different for implant roughness lower than 20 μm. 2D modeling is sufficient to describe the interaction between ultrasound and the BII.

A cadaveric validation of a method based on impact analysis to monitor the femoral stem insertion

The success of cementless hip arthroplasty depends on the primary stability of the femoral stem (FS). It remains difficult to assess the optimal impaction energy to guarantee the FS stability while avoiding bone fracture. The aim of this study is to compare the results of a method based on the use of an instrumented hammer to determine the insertion endpoint of cementless FS in a cadaveric model with two other methods using i) the surgeon proprioception and ii) video motion tracking. Different FS were impacted in nine human cadaveric femurs. For each configuration, the number of impacts realized when the surgeon felt that the FS was correctly inserted was noted N_surg. For each impact, the insertion depth E was measured and an indicator D was determined based on the time-variation of the force. The impact number N_vid (respectively N_d), corresponding to the end of the migration phase, was estimated analyzing the evolution of E (respectively D). The respective difference between N_surg, N_vid and N_d was similar and lower than 3 for more than 85% of the configurations. The results allow a validation of the use of an impact hammer to assess the moment when the surgeon should stop the impaction, paving the way towards the development of a decision support system to assist the surgeon.

Elastography of the bone-implant interface

The stress distribution around endosseous implants is an important determinant of the surgical success. However, no method developed so far to determine the implant stability is sensitive to the loading conditions of the bone-implant interface (BII). The objective of this study is to investigate whether a quantitative ultrasound (QUS) technique may be used to retrieve information on compressive stresses applied to the BII. An acousto-mechanical device was conceived to compress 18 trabecular bovine bone samples onto coin-shaped implants and to measure the ultrasonic response of the BII during compression. The biomechanical behavior of the trabecular bone samples was modeled as Neo-Hookean. The reflection coefficient of the BII was shown to decrease as a function of the stress during the elastic compression of the trabecular bone samples and during the collapse of the trabecular network, with an average slope of −4.82 GPa⁻¹. The results may be explained by an increase of the bone-implant contact ratio and by changes of bone structure occurring during compression. The sensitivity of the QUS response of the BII to compressive stresses opens new paths in the elaboration of patient specific decision support systems allowing surgeons to assess implant stability that should be developed in the future.

Biomechanical behaviours of the bone-implant interface: a review

In recent decades, cementless implants have been widely used in clinical practice to replace missing organs, to replace damaged or missing bone tissue or to restore joint functionality. However, there remain risks of failure which may have dramatic consequences. The success of an implant depends on its stability, which is determined by the biomechanical properties of the bone-implant interface (BII). The aim of this review article is to provide more insight on the current state of the art concerning the evolution of the biomechanical properties of the BII as a function of the implant's environment. The main characteristics of the BII and the determinants of implant stability are first introduced. Then, the different mechanical methods that have been employed to derive the macroscopic properties of the BII will be described. The experimental multi-modality approaches used to determine the microscopic biomechanical properties of periprosthetic newly formed bone tissue are also reviewed. Eventually, the influence of the implant's properties, in terms of both surface properties and biomaterials, is investigated. A better understanding of the phenomena occurring at the BII will lead to (i) medical devices that help surgeons to determine an implant's stability and (ii) an improvement in the quality of implants.

Reflection of an ultrasonic wave on the bone-implant interface: Effect of the roughness parameters

Quantitative ultrasound can be used to characterize the evolution of the bone-implant interface (BII), which is a complex system due to the implant surface roughness and to partial contact between bone and the implant. The aim of this study is to derive the main determinants of the ultrasonic response of the BII during osseointegration phenomena. The influence of (i) the surface roughness parameters and (ii) the thickness W of a soft tissue layer on the reflection coefficient r of the BII was investigated using a two-dimensional finite element model. When W increases from 0 to 150 μm, r increases from values in the range [0.45; 0.55] to values in the range [0.75; 0.88] according to the roughness parameters. An optimization method was developed to determine the sinusoidal roughness profile leading to the most similar ultrasonic response for all values of W compared to the original profile. The results show that the difference between the ultrasonic responses of the optimal sinusoidal profile and of the original profile was lower to typical experimental errors. This approach provides a better understanding of the ultrasonic response of the BII, which may be used in future numerical simulation realized at the scale of an implant.

Resonant frequency analysis of dental implants

Dental implant stability influences the decision on the determination of the duration between implant insertion and loading. This work investigates the resonant frequency analysis by means of a numerical model. The investigation is done numerically through the determination of the eigenfrequencies and performing steady state response analyses using a commercial finite element package. A peri-implant interface, of simultaneously varying stiffness, density and layer thickness is introduced in the numerical 3D model in order to probe the sensitivity of the eigenfrequencies and steady state response to an evolving weakened layer, in an attempt to identify the bone reconstruction around the implant. For the first two modes, the resonant frequency is somewhat insensitive to the healing process, unless the weakened layer is rather large and compliant, like in the very early stages of the implantation. A "Normalized Healing Factor" is devised in the spirit of the Implant Stability Quotient, which can identify the healing process especially at the early stages after implantation. The sensitivity of the resonant frequency analysis to changes of mechanical properties of periprosthetic bone tissue seems relatively weak. Another indicator considering the amplitude as well as the resonance frequency might be more adapted to bone healing estimations. However, these results need to be verified experimentally as well as clinically.

Effect of intracortical bone properties on the phase velocity and cut-off frequency of low-frequency guided wave modes (20-85 kHz)

The assessment of intracortical bone properties is of interest since early-stage osteoporosis is associated with resorption in the endosteal region. However, understanding the interaction between ultrasonic guided waves and the cortical bone structure remains challenging. The purpose of this work is to investigate the effect of intracortical bone properties on the ultrasonic response obtained at low-frequency (<100 kHz) using an axial transmission configuration. The semi-analytical finite element method was used to simulate the propagation of guided waves in a waveguide with realistic geometry and material properties. An array of 20 receivers was used to calculate the phase velocity and cut-off frequency of the excited modes using the two-dimensional Fourier transform. The results show that the position of the emitter around the circumference of the bone is an important parameter to control since it can lead to variations of up to 10 dB in the amplitude of the transmitted modes. The cut-off frequency of the high order modes was, however, only slightly affected by the circumferential position of the emitter, and was sensitive mainly to the axial shear modulus. The phase velocity and cut-off frequency in the 20-85 kHz range are promising parameters for the assessment of intracortical properties.

Monitoring cementless femoral stem insertion by impact analyses: An in vitro study

The primary stability of the femoral stem (FS) implant determines the surgical success of cementless hip arthroplasty. During the insertion, a compromise must be found for the number and energy of impacts that should be sufficiently large to obtain an adapted primary stability of the FS and not too high to decrease fracture risk. The aim of this study is to determine whether a hammer instrumented with a force sensor can be used to monitor the insertion of FS. Cementless FS of different sizes were impacted in four artificial femurs with an instrumented hammer, leading to 72 configurations. The impact number when the surgeon empirically felt that the FS was fully inserted was noted N_surg. The insertion depth E was assessed using video motion tracking and the impact number N_vid corresponding to the end of the insertion was estimated. For each impact, two indicators noted I and D were determined based on the analysis of the variation of the force as a function of time. The pull-out force F was significantly correlated with the indicator I (R² = 0.67). The variation of D was analyzed using a threshold to determine an impact number N_d, which is shown to be closely related to N_surg and N_vid, with an average difference of around 0.2. This approach allows to determine i) the moment when the surgeon should stop the impaction procedure in order to obtain an optimal insertion of the FS and ii) the FS implant primary stability. This study paves the way towards the development of a decision support system to assist the surgeon in hip arthroplasty.

Reflection of an ultrasonic wave on the bone-implant interface: A numerical study of the effect of the multiscale roughness

Quantitative ultrasound is used to characterize and stimulate osseointegration processes at the bone-implant interface (BII). However, the interaction between an ultrasonic wave and the implant remains poorly understood. This study aims at investigating the sensitivity of the ultrasonic response to the microscopic and macroscopic properties of the BII and to osseointegration processes. The reflection coefficient R of the BII was modeled for different frequencies using a two-dimensional finite element model. The implant surface roughness was modeled by a sinusoidal function with varying amplitude h and spatial frequency L. A soft tissue layer of thickness W was considered between bone tissue and the implant in order to model non-mineralized fibrous tissue. For microscopic roughness, R is shown to increase from around 0.55 until 0.9 when kW increases from 0 to 1 and to be constant for kW > 1, where k is the wavenumber in the implant. These results allow us to show that R depends on the properties of bone tissue located at a distance comprised between 1 and 25 μm from the implant surface. For macroscopic roughness, R is highly dependent on h and this dependence may be explained by phase cancellation and multiple scattering effects for high roughness parameters.

Ex Vivo Evaluation of Cementless Acetabular Cup Stability Using Impact Analyses with a Hammer Instrumented with Strain Sensors

The acetabular cup (AC) implant stability is determinant for the success of cementless hip arthroplasty. A method based on the analysis of the impact force applied during the press-fit insertion of the AC implant using a hammer instrumented with a force sensor was developed to assess the AC implant stability. The aim of the present study was to investigate the performance of a method using a hammer equipped with strain sensors to retrieve the AC implant stability. Different AC implants were inserted in five bovine samples with different stability conditions leading to 57 configurations. The AC implant was impacted 16 times by the two hammers consecutively. For each impact; an indicator I_S (respectively I_F) determined by analyzing the time variation of the signal corresponding to the averaged strain (respectively force) obtained with the stress (respectively strain) hammer was calculated. The pull-out force F was measured for each configuration. F was significantly correlated with I_S (R² = 0.79) and I_F (R² = 0.80). The present method has the advantage of not modifying the shape of the hammer that can be sterilized easily. This study opens new paths towards the development of a decision support system to assess the AC implant stability.

Assessment of the biomechanical stability of a dental implant with quantitative ultrasound: A three-dimensional finite element study

Dental implant stability is an important determinant of the surgical success. Quantitative ultrasound (QUS) techniques can be used to assess such properties using the implant acting as a waveguide. However, the interaction between an ultrasonic wave and the implant remains poorly understood. The aim of this study is to investigate the sensitivity of the ultrasonic response to the quality and quantity of bone tissue in contact with the implant surface. The 10 MHz ultrasonic response of an implant used in clinical practice was simulated using an axisymmetric three-dimensional finite element model, which was validated experimentally. The amplitude of the echographic response of the implant increases when the depth of a liquid layer located at the implant interface increases. The results show the sensitivity of the QUS technique to the amount of bone in contact with the implant. The quality of bone tissue around the implant is varied by modifying the bone biomechanical properties by 20%. The amplitude of the implant echographic response decreases when bone quality increases, which corresponds to bone healing. In all cases, the amplitude of the implant response decreased when the dental implant stability increased, which is consistent with the experimental results.

Ultrasound Speed of Sound Measurements in Trabecular Bone Using the Echographic Response of a Metallic Pin

Bone quality is an important parameter in spine surgery, but its clinical assessment remains difficult. The aim of the work described here was to demonstrate in vitro the feasibility of employing quantitative ultrasound to retrieve bone mechanical properties using an echographic technique taking advantage of the presence of a metallic pin inserted in bone tissue. A metallic pin was inserted in bone tissue perpendicular to the transducer axis. The echographic response of the bone sample was determined, and the echo of the pin inserted in bone tissue and water were compared to determine speed of sound, which was compared with bone volume fraction. A 2-D finite-element model was developed to assess the effect of positioning errors. There was a significant correlation between speed of sound and bone volume fraction (R(2) = 0.6). The numerical results indicate the relative robustness of the measurement method, which could be useful to estimate bone quality intra-operatively.

Assessment of in vitro dental implant primary stability using an ultrasonic method

Dental implants are used for oral rehabilitation. However, there remain risks of failure that depend on the implant stability. The objective of this study is to investigate whether quantitative ultrasound technique can be used to assess the amount of bone in contact with dental implants. Ten implants are first inserted in the bone samples. The 10 MHz ultrasonic response of each implant is measured using a dedicated device and an indicator I is derived based on the amplitude of the signal. Then, the implant is unscrewed by 2 π radians and the measurement is realized again. A statistical analysis of variance was carried out and revealed a significant effect of the amount of bone in contact with the implant on the values of I (p value < 10⁻⁵). The results indicates the feasibility of quantitative ultrasound techniques to assess implant primary stability in vitro.

Effects of biomechanical properties of the bone-implant interface on dental implant stability: from in silico approaches to the patient's mouth

Dental implants have become a routinely used technique in dentistry for replacing teeth. However, risks of failure are still experienced and remain difficult to anticipate. Multiscale phenomena occurring around the implant interface determine the implant outcome. The aim of this review is to provide an understanding of the biomechanical behavior of the interface between a dental implant and the region of bone adjacent to it (the bone-implant interface) as a function of the interface's environment. First, we describe the determinants of implant stability in relation to the different multiscale simulation approaches used to model the evolution of the bone-implant interface. Then, we review the various aspects of osseointegration in relation to implant stability. Next, we describe the different approaches used in the literature to measure implant stability in vitro and in vivo. Last, we review various factors affecting the evolution of the bone-implant interface properties.

Radial anatomic variation of ultrasonic velocity in human cortical bone

Quantitative ultrasound techniques can be used to retrieve cortical bone quality. The aim of this study was to investigate the anatomic variations in speed of sound (SOS) in the radial direction of cortical bone tissue. SOS measurements were realized in 17 human cortical bone samples with a 3.5-MHz transverse transmission device. The radial dependence of SOS was investigated in a direction perpendicular to the periosteum. For each sample, bone porosity was measured using an X-ray micro-computed tomography device. The mean SOS was 3586 ± 255 m/s. For 16 of 17 specimens, similar radial variations in SOS were observed. In the periosteal region, SOS first decreased in the direction of the endosteum and reached a minimum value approximately in the middle of the cortical bone. SOS then increased, moving to the endosteal region. A significant negative correlation was obtained between SOS and porosity (R = -0.54, p = 0.02).

Anatomical distribution of the degree of mineralization of bone tissue in human femoral neck: impact on biomechanical properties

Osteoporotic hip fractures represent a major public health problem associated with high human and economic costs. The anatomical variation of the tissue mineral density (TMD) and of the elastic constants in femoral neck cortical bone specimens is an important determinant of bone fragility. The purpose of this study was to show that a Synchrotron radiation microcomputed tomography system coupled with a multiscale biomechanical model allows the determination of the 3-D anatomical dependence of TMD and of the elastic constants (i.e. the mechanical properties of an anisotropic material) in human femoral neck. Bone specimens from the inferior femoral neck were obtained from 18 patients undergoing standard hemiarthroplasty. The specimens were imaged using 3-D synchrotron micro-computed tomography with a voxel size of 10.13 μm, leading to the determination of the anatomical distributions of porosity and TMD. The elastic properties of bone tissue were computed using a multiscale model. The model uses the experimental data obtained at the scale of several micrometers to estimate the components of the elastic tensor of bone at the scale of the organ. Statistical analysis (ANOVA) revealed a significant effect of the radial position on porosity and TMD and a significant effect of axial position on TMD only. Porosity was found to increase in the radial direction moving from the periosteum inwards (p<10(-5)). At any given distance from the periosteum, porosity does not vary noticeably along the bone axis. TMD was found to be significantly higher (p<10(-5)) in the periosteal region than in other bone locations and decreases from the periosteal to the endosteal region with an average slope of 10.05 g.cm(-3).m(-1), the decrease being faster in the porous part of the samples (average slope equal of 30.04 g.cm(-3).m(-1)) than in dense cortical bone. TMD was found to decrease from the distal to the proximal part of the femur neck (average slope of 6.5 g.cm(-3).m(-1)). Considering TMD variations in the radial direction induces weak changes of bone properties compared to constant TMD. TMD variations in the axial direction are responsible for a significant variation of elastic constants. These results demonstrate that the anatomical variations of TMD affect the bone elastic properties, which could be explained by the complex stress field in bone affecting bone remodeling. TMD spatial variations should be taken into account to properly describe the spatial heterogeneity of elastic coefficients of bone tissue at the organ scale.

Influence of healing time on the ultrasonic response of the bone-implant interface

The aim of the present study is to investigate the effect of bone healing on the ultrasonic response of coin-shaped titanium implants inserted in rabbit tibiae. The ultrasound response of the interface was measured in vitro at 15 MHz after 7 and 13 weeks of healing time. The average value of the ratio r between the amplitudes of the echo of the bone-implant interface and of the water-implant interface was determined. The bone-implant contact (BIC) was measured by histomorphometry and the degree of mineralisation of bone was estimated qualitatively by histologic staining. The significant decrease of the ultrasonic quantitative indicator r (p = 2.10⁻⁴) vs. healing time (from r = 0.53 to r = 0.49) is explained by (1) the increase of the BIC (from 27% to 69%) and (2) the increase of mineralization of newly formed bone tissue, both phenomena inducing a decrease of the gap of acoustical impedance.

Equivalent contributing depth investigated by a lateral wave with axial transmission in viscoelastic cortical bone

Cortical bone is a viscoelastic heterogeneous medium which may be assessed with axial transmission. This work aims at evaluating the average depth investigated by the lateral wave for radial variations of material properties in relatively thick cortical bone. The equivalent contributing depth (ECD) is derived from the finite element simulation results for spatial variations of a viscoelastic coefficient (η(11)) and of porosity. A value of ECD equal to around 1.6 mm is obtained for a spatial variation of η(11). The method fails to predict accurate values of the ECD for a spatial variation of porosity, because all parameters vary simultaneously.

Two-dimensional simulation of linear wave propagation in a suspension of polymeric microcapsules used as ultrasound contrast agents

A generation of tissue-specific stable ultrasound contrast agent (UCA) composed of a polymeric capsule with a perfluorocarbone liquid core has become available. Despite promising uses in clinical practice, the acoustical behavior of such UCA suspensions remains unclear. A simulation code (2-D finite-difference time domain, FDTD) already validated for homogeneous particles [Galaz Haiat, Berti, Taulier, Amman and Urbach, (2010). J. Acoust. Soc. Am. 127, 148-154] is used to model the ultrasound propagation in such UCA suspensions at 50 MHz to investigate the sensitivity of the ultrasonic parameters to physical parameters of UCA. The FDTD simulation code is validated by comparison with results obtained using a shell scatterer model. The attenuation coefficient (respectively, the sound velocity) increases (respectively, decreases) from 4.1 to 58.4 dB/cm (respectively, 1495 to 1428 m/s) when the concentration varies between 1.37 and 79.4 mg/ml, while the backscattered intensity increases non-linearly, showing that a concentration of around 30 mg/ml is sufficient to obtain optimal backscattering intensity. The acoustical parameters vary significantly as a function of the membrane thickness, longitudinal and transverse velocity, indicating that mode conversions in the membrane play an important role in the ultrasonic propagation. The results may be used to help manufacturers to conceive optimal liquid-filled UCA suspensions.

Ultrasonic evaluation of dental implant biomechanical stability: an in vitro study

Dental implants are widely used for oral rehabilitation. However, there remain risks of failure that are difficult to anticipate. The objective of this ex vivo study is to investigate the potentiality of quantitative ultrasound (QUS) to assess the amount of bone in contact with titanium prototype cylindrical implants. Four groups of 10 rabbit femurs each are considered, corresponding to different amounts of bone in contact with the implant. The 10 MHz ultrasonic response of the implant is processed to derive a quantitative indicator I, based on the temporal variation of the signal amplitude. Analysis of variance (ANOVA) (p < 10(-5)) tests revealed a statistical distribution of I significantly correlated with the amount of bone in contact with the cylinders. An analytical model considering the propagation of lateral waves allows the understanding of the physical origin of the echoes. QUS technique may be used to investigate the amount of bone in contact with a cylinder implant.

Determination of the heterogeneous anisotropic elastic properties of human femoral bone: from nanoscopic to organ scale

Cortical bone is a multiscale composite material. Its elastic properties are anisotropic and heterogeneous across its cross-section, due to endosteal bone resorption which might affect bone strength. The aim of this paper was to describe a homogenization method leading to the estimation of the variation of the elastic coefficients across the bone cross-section and along the bone longitudinal axis. The method uses the spatial variations of bone porosity and of the degree of mineralization of the bone matrix (DMB) obtained from the analysis of 3-D synchrotron micro-computed tomography images. For all three scales considered (the foam (100 nm), the ultrastructure (5 microm) and the mesoscale (500 microm)), the elastic coefficients were determined using the Eshelby's inclusion problem. DMB values were used at the scale of the foam. Collagen was introduced at the scale of the ultrastructure and bone porosity was introduced at the mesoscale. The pores were considered as parallel cylinders oriented along the bone axis. Each elastic coefficient was computed for different regions of interest, allowing an estimation of its variations across the bone cross-section and along the bone longitudinal axis. The method was applied to a human femoral neck bone specimen, which is a site of osteoporotic fracture. The computed elastic coefficients for cortical bone were in good agreement with experimental results, but some discrepancies were obtained in the endosteal part (trabecular bone). These results highlight the importance of accounting for the heterogeneity of cortical bone properties across bone cross-section and along bone longitudinal axis.

Influence of viscoelastic and viscous absorption on ultrasonic wave propagation in cortical bone: Application to axial transmission

Cortical bone and the surrounding soft tissues are attenuating and heterogeneous media, which might affect the signals measured with axial transmission devices. This work aims at evaluating the effect of the heterogeneous acoustic absorption in bone and in soft tissues on the bone ultrasonic response. Therefore, a two-dimensional finite element time-domain method is derived to model transient wave propagation in a three-layer medium composed of an inhomogeneous transverse isotropic viscoelastic solid layer, sandwiched between two viscous fluid layers. The model couples viscous acoustic propagation in both fluid media with the anisotropic viscoelastic response of the solid. A constant spatial gradient of material properties is considered for two values of bone thicknesses (0.6 and 4 mm). In the studied configuration, absorption in the surrounding fluid tissues does not affect the results, whereas bone viscoelastic properties have a significant effect on the first arriving signal (FAS) velocity. For a thin bone, the FAS velocity is governed by the spatially averaged bone properties. For a thick bone, the FAS velocity may be predicted using a one-dimensional model.

Experimental validation of a time domain simulation of high frequency ultrasonic propagation in a suspension of rigid particles

Ultrasonic propagation in suspensions of particles is a difficult problem due to the random spatial distribution of the particles. Two-dimensional finite-difference time domain simulations of ultrasonic propagation in suspensions of polystyrene 5.3 mum diameter microdisks are performed at about 50 MHz. The numerical results are compared with the Faran model, considering an isolated microdisk, leading to a maximum difference of 15% between the scattering cross-section values obtained analytically and numerically. Experiments are performed with suspensions in through transmission and backscattering modes. The attenuation coefficient at 50 MHz (alpha), the ultrasonic velocity (V), and the relative backscattered intensity (I(B)) are measured for concentrations from 2 to 25 mg/ml, obtained by modifying the number of particles. Each experimental ultrasonic parameter is compared to numerical results obtained by averaging the results derived from 15 spatial distributions of microdisks. alpha increases with the concentration from 1 to 17 dB/cm. I(B) increases with concentration from 2 to 16 dB. The variation of V versus concentration is compared with the numerical results, as well as with an effective medium model. A good agreement is found between experimental and numerical results (the larger discrepancy is found for alpha with a difference lower than 2.1 dB/cm).

Relationship between ultrasonic parameters and apparent trabecular bone elastic modulus: a numerical approach

The physical principles underlying quantitative ultrasound (QUS) measurements in trabecular bone are not fully understood. The translation of QUS results into bone strength remains elusive. However, ultrasound being mechanical waves, it is likely to assess apparent bone elasticity. The aim of this study is to derive the sensitivity of QUS parameters to variations of apparent bone elasticity, a surrogate for strength. The geometry of 34 human trabecular bone samples cut in the great trochanter was reconstructed using 3-D synchrotron micro-computed tomography. Finite-difference time-domain simulations coupled to 3-D micro-structural models were performed in the three perpendicular directions for each sample and each direction. A voxel-based micro-finite element linear analysis was employed to compute the apparent Young's modulus (E) of each sample for each direction. For the antero-posterior direction, the predictive power of speed of sound and normalized broadband ultrasonic attenuation to assess E was equal to 0.9 and 0.87, respectively, which is better than what is obtained using bone density alone or coupled with micro-architectural parameters and of the same order of what can be achieved with the fabric tensor approach. When the direction of testing is parallel to the main trabecular orientation, the predictive power of QUS parameters decreases and the fabric tensor approach always gives the best results. This decrease can be explained by the presence of two longitudinal wave modes. Our results, which were obtained using two distinct simulation tools applied on the same set of samples, highlight the potential of QUS techniques to assess bone strength.

Influence of a gradient of material properties on ultrasonic wave propagation in cortical bone: application to axial transmission

The aim of this work is to evaluate the effect of a spatial gradient of material properties (mass density and stiffness coefficients) of cortical bone on its ultrasonic response obtained with an axial transmission device. Therefore, a two-dimensional finite element time-domain method is derived to model transient wave propagation in a three-layer medium composed of an inhomogeneous transverse isotropic solid layer sandwiched between two acoustic fluid layers and excited by an acoustic linear source located in one fluid layer, delivering broadband ultrasonic pulses. The model couples the acoustic propagation in both fluid media with the elastodynamic response of the solid layer. A constant spatial gradient of material properties is considered for two values of bone thicknesses corresponding to relatively thick and thin bone widths. For a thin bone (0.6 mm) compared to wavelength (around 4 mm at 1 MHz), the results are in good agreement with a S(0) Lamb wave assuming a homogeneous material with spatially averaged material properties. For a thick bone (4 mm), the results are in agreement with the propagation of a lateral wave and allow the derivation of an equivalent contributing depth in the case of a transverse isotropic inhomogeneous solid layer.

Velocity dispersion in trabecular bone: influence of multiple scattering and of absorption

Speed of sound measurements are widely used clinically to assess bone strength. Trabecular bone is an attenuating composite material in which negative values of velocity dispersion have been measured, this behavior remaining poorly explained physically. The aim of this work is to describe the ultrasonic propagation in trabecular bone modeled by infinite cylinders immersed in a saturating matrix, and to derive the physical determinants of velocity dispersion. A homogenization model accounting for the coupling of multiple scattering and absorption phenomena allows the computation of phase velocity and of dispersion while varying bone properties. The present model is adapted from the generalized self-consistent method developed in the work of Yang and Mal [(1994). "Multiple-scattering of elastic waves in a fiber-reinforced composite," J. Mech. Phys. Solids 42, 1945-1968]. It predicts negative values of velocity dispersion, in agreement with experimental results obtained in phantoms mimicking trabecular bone. In trabecular bone, mostly negative and also positive values of velocity dispersion are predicted by the model, which span within the range of values measured experimentally. Scattering effects are responsible for the negative values of dispersion, whereas the frequency dependence of the absorption coefficient in bone marrow and/or in the trabeculae results in an increase in dispersion, which may then become positive.

Ultrasonic velocity dispersion in bovine cortical bone: an experimental study

Cortical bone quality is determinant in bone fragility and its ultrasonic evaluation has become possible in clinical practice. However, the interaction between a broadband ultrasonic pulse and this complex multiscale medium remains poorly understood. The frequency dependence of phase velocity, which may impact clinical measurements, has been sparsely investigated. Our objective is to evaluate the determinants of the frequency dependence of phase velocity in bovine femoral cortical bone samples using an in vitro ultrasonic transmission device. The apparent phase velocity varies quasilinearly on a 1 MHz restricted bandwidth around 4 MHz. After compensating for diffraction effects, significant differences in velocity dispersion are obtained according to the anatomical location. The microstructure of each sample is determined using an optical microscope, which allows assessing the dependence of dispersion on the type of bone microstructure. Mostly positive but also negative values of dispersion are measured. Negative dispersion is mainly obtained in samples constituted of mixed microstructure, which may be explained by phase cancellation effects due to the presence of different microstructures within the same sample. Dispersion is shown to be related to broadband ultrasonic attenuation values, especially in the radial direction. Results are compared with the local Kramers-Kronig relationships.

Fast wave ultrasonic propagation in trabecular bone: numerical study of the influence of porosity and structural anisotropy

Our goal is to assess the potential of computational methods as an alternative to analytical models to predict the two longitudinal wave modes observed in cancellous bone and predicted by the Biot theory. A three-dimensional (3D) finite-difference time-domain method is coupled with 34 human femoral trabecular microstructures measured using microcomputed tomography. The main trabecular alignment (MTA) and the degree of anisotropy (DA) were assessed for all samples. DA values were comprised between 1.02 and 1.9. The influence of bone volume fraction (BV/TV) between 5% and 25% on the properties of the fast and slow waves was studied using a dedicated image processing algorithm to modify the initial 3D microstructures. A heuristic method was devised to determine when both wave modes are time separated. The simulations (performed in three perpendicular directions) predicted that both waves generally overlap in time for a direction of propagation perpendicular to the MTA. When these directions are parallel, both waves are separated in time for samples with high DA and BV/TV values. A relationship was found between the least bone volume fraction required for the observation of nonoverlapping waves and the degree of anisotropy: The higher the DA, the lower the least BV/TV.