Cortical bone plate properties assessment using inversion of axially transmitted low frequency ultrasonic guided waves.

Over the past few decades, early osteoporosis detection using ultrasonic bone quality evaluation has gained prominence. Specifically, various studies focused on axial transmission using ultrasonic guided waves and have highlighted this technique's sensitivity to intrinsic properties of long cortical bones. This work aims to demonstrate the potential of low-frequency ultrasonic guided waves to infer the properties of the bone inside which they are propagating. A proprietary ultrasonic transducer, tailored to transmit ultrasonic guided waves under 500 kHz, was used for the data collection. The gathered data underwent two-dimensional fast Fourier transform processing to extract experimental dispersion curves. The proposed inversion scheme compares experimental dispersion curves with simulated dispersion curves calculated through the semi-analytical iso-geometric analysis (SAIGA) method. The numerical model integrates a bone phantom plate coupled with a soft tissue layer on its top surface, mimicking the experimental bone phantom plates. Subsequently, the mechanical properties of the bone phantom plates were estimated by reducing the misfit between the experimental and simulated dispersion curves. This inversion leaned heavily on the dispersive trajectories and amplitudes of ultrasonic guided wave modes. Results indicate a marginal discrepancy under 5% between the mechanical properties ascertained using the SAIGA-based inversion and those measured using bulk wave pulse-echo measurements.

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.

Ultrasonic Evaluation of the Bone-Implant Interface.

While implant surgical interventions are now routinely performed, failures still occur and may have dramatic consequences. The clinical outcome depends on the evolution of the biomechanical properties of the bone-implant interface (BII). This chapter reviews studies investigating the use of quantitative ultrasound (QUS) techniques for the characterization of the BII.First, studies on controlled configurations evidenced the influence of healing processes and of the loading conditions on the ultrasonic response of the BII. The gap of acoustical properties at the BII increases (i) during healing and (ii) when stress at the BII increases, therefore inducing a decrease of the reflection coefficient at the BII.Second, an acoustical model of the BII is proposed to better understand the parameters influencing the interaction between ultrasound and the BII. The reflection coefficient is shown to decrease when (i) the BII is better osseointegrated, (ii) the implant roughness decreases, (iii) the frequency of QUS decreases and (iv) the bone mass density increases.Finally, a 10 MHz device aiming at assessing dental implant stability was validated in vitro, in silico and in vivo. A comparison between QUS and resonance frequency analysis (RFA) techniques showed a better sensitivity of QUS to changes of the parameters related to the implant stability.

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 R2 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.

Using an Impact Hammer to Estimate Elastic Modulus and Thickness of a Sample During an Osteotomy.

Performing an osteotomy with a surgical mallet and an osteotome is a delicate intervention mostly based on the surgeon proprioception. It remains difficult to assess the properties of bone tissue being osteotomized. Mispositioning of the osteotome or too strong impacts may lead to bone fractures which may have dramatic consequences. The objective of this study is to determine whether an instrumented hammer may be used to retrieve information on the material properties around the osteotome tip. A hammer equipped with a piezo-electric force sensor was used to impact 100 samples of different composite materials and thicknesses. A model-based inversion technique was developed based on the analysis of two indicators derived from the analysis of the variation of the force as a function of time in order to (i) classify the samples depending on their material types, (ii) determine the materials stiffness, and (iii) estimate the samples thicknesses. The model resulting from the classification using support vector machines (SVM) learning techniques can efficiently predict the material of a new sample, with an estimated 89% prediction performance. A good agreement between the forward analytical model and the experimental data was obtained, leading to an average error lower than 10% in the samples thickness estimation. Based on these results, navigation and decision-support tools could be developed and allows surgeons to adapt their surgical strategy in a patient-specific manner.

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.

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.

Micromechanical modeling of the contact stiffness of an osseointegrated bone-implant interface.

BACKGROUND: The surgical success of cementless implants is determined by the evolution of the biomechanical properties of the bone-implant interface (BII). One difficulty to model the biomechanical behavior of the BII comes from the implant surface roughness and from the partial contact between bone tissue and the implant. The determination of the constitutive law of the BII would be of interest in the context of implant finite element (FE) modeling to take into account the imperfect characteristics of the BII. The aim of the present study is to determine an effective contact stiffness [Formula: see text] of an osseointegrated BII accounting for its micromechanical features such as surface roughness, bone-implant contact ratio (BIC) and periprosthetic bone properties. To do so, a 2D FE model of the BII under normal contact conditions was developed and was used to determine the behavior of [Formula: see text]. RESULTS: The model is validated by comparison with three analytical schemes based on micromechanical homogenization including two Lekesiz's models (considering interacting and non-interacting micro-cracks) and a Kachanov's model. [Formula: see text] is found to be comprised between 1013 and 1015 N/m3 according to the properties of the BII. [Formula: see text] is shown to increase nonlinearly as a function of the BIC and to decrease as a function of the roughness amplitude for high BIC values (above around 20%). Moreover, [Formula: see text] decreases as a function of the roughness wavelength and increases linearly as a function of the Young's modulus of periprosthetic bone tissue. CONCLUSIONS: These results open new paths in implant biomechanical modeling since this model may be used in future macroscopic finite element models modeling the bone-implant system to replace perfectly rigid BII conditions.

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.

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.

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.

Comparison of Resonance Frequency Analysis and of Quantitative Ultrasound to Assess Dental Implant Osseointegration.

Dental implants are widely used in the clinic. However, there remain risks of failure, which depend on the implant stability. The aim of this paper is to compare two methods based on resonance frequency analysis (RFA) and on quantitative ultrasound (QUS) and that aim at assessing implant stability. Eighty-one identical dental implants were inserted in the iliac crests of 11 sheep. The QUS and RFA measurements were realized after different healing times (0, 5, 7, and 15 weeks). The results obtained with the QUS (respectively RFA) method were significantly different when comparing two consecutive healing time for 97% (respectively, 18%) of the implants. The error made on the estimation of the healing time when analyzing the results obtained with the QUS technique was around 10 times lower than that made when using the RFA technique. The results corresponding to the dependence of the ISQ versus healing time were significantly different when comparing two directions of RFA measurement. The results show that the QUS method allows a more accurate determination of the evolution of dental implant stability when compared to the RFA method. This study paves the way towards the development of a medical device, thus providing a decision support system to dental surgeons.

Simulation of acoustic guided wave propagation in cortical bone using a semi-analytical finite element method.

Axial transmission techniques have been extensively studied for cortical bone quality assessment. However, the modeling of ultrasonic guided waves propagation in such a complex medium remains challenging. The aim of this paper is to develop a semi-analytical finite element method to simulate the propagation of guided waves in an irregular, multi-layer, and heterogeneous bone cross-section modeled with anisotropic and viscoelastic material properties. The accuracy of the simulations was verified against conventional time-domain three-dimensional finite element. The method was applied in the context of axial transmission in bone to investigate the feasibility of first arrival signal (FAS) to monitor degradation of intracortical properties at low frequency. Different physiopathological conditions for the intracortical region, varying from healthy to osteoporotic, were monitored through FAS velocity using a 10-cycle tone burst excitation centered at 32.5 kHz. The results show that the variation in FAS velocity is mainly associated with four of the eight modes supported by the waveguide, varying with velocity values between 550 and 700 m/s along the different scenarios. Furthermore, the FAS velocity is shown to be associated with the group velocity of the mode with the highest relative amplitude contribution at each studied scenario. However, because of the evolution of the mode with the highest contribution, the FAS velocity is shown to be limited to discriminate intracortical bone properties at low frequency.

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 IS (respectively IF) 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 IS (R² = 0.79) and IF (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.

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.

In vitro evaluation of the acetabular cup primary stability by impact analysis.

The implant primary stability of the acetabular cup (AC) is an important parameter for the surgical success of press-fit procedures used for the insertion of cementless hip prostheses. In previous studies by our group (Mathieu, V., Michel, A., Lachaniette, C. H. F., Poignard, A., Hernigou, P., Allain, J., and Haiat, G., 2013, "Variation of the Impact Duration During the in vitro Insertion of Acetabular Cup Implants," Med. Eng. Phys., 35(11), pp. 1558-1563) and (Michel, A., Bosc, R., Mathieu, V., Hernigou, P., and Haiat, G., 2014, "Monitoring the Press-Fit Insertion of an Acetabular Cup by Impact Measurements: Influence of Bone Abrasion," Proc. Inst. Mech. Eng., Part H, 228(10), pp. 1027-1034), the impact momentum and duration were shown to carry information on the press-fit insertion of the AC within bone tissue. The aim of the present study is to relate the impact momentum recorded during the AC insertion to the AC biomechanical primary stability. The experimental protocol consisted in testing 13 bovine bone samples that underwent successively series of 15 reproducible mass falls impacts (5 kg, 5 cm) followed by tangential stability testing. Each bone sample was tested with different hole sizes in order to obtain different stability configurations. The impact momentum and the tangential primary stability reach a maximum value for an interference fit equal to around 1 mm. Moreover, a correlation between the impact momentum and the stability was obtained with all samples and all configuration (R2 = 0.65). The implant primary stability can be assessed through the measurement of the impact force signal analysis. This study opens new paths for the development of a medical device which could be used as a decision support system to assist the surgeon during the insertion of the AC implant.

Effect of collagen fibril orientation on the anisotropic properties of peri-implant bone.

In orthopedic and dental surgery, the implantation of biomaterials within the bone to restore the integrity of the treated organ has become a standard procedure. Their long-term stability relies on the osseointegration phenomena, where bone grows onto and around metallic implants, creating a bone-implant interface. Bone is a highly hierarchical material that evolves spatially and temporally during this healing phase. A deeper understanding of its biomechanical characteristics is needed, as they are determinants for surgical success. In this context, we propose a multiscale homogenization model to evaluate the effective elastic properties of bone as a function of the distance from the implant, based on the tissue's structure and composition at lower scales. The model considers three scales: hydroxyapatite foam (nanoscale), ultrastructure (microscale), and tissue (mesoscale). The elastic properties and the volume fraction of the elementary constituents of bone matrix (mineral, collagen, and water), the orientation of the collagen fibril relative to the implant surface, and the mesoscale porosity constitute the input data of the model. The effect of a spatiotemporal variation in the collagen fibrils' orientation on the bone anisotropic properties in the proximity of the implant was investigated. The findings revealed a strong variation of the components of the effective elasticity tensor of the bone as a function of the distance from the implant. The effective elasticity appears to be primarily sensitive to the porosity (mesoscale) rather than to the collagen fibrils' orientation (sub-micro scale). However, the orientation of the fibrils has a significant influence on the isotropy of the bone. When analyzing the symmetry properties of the effective elasticity tensor, the ratio between the isotropic and hexagonal components is determined by a combination of the porosity and the fibrils' orientation. A decrease in porosity leads to a decrease in bone isotropy and, in turn, an increase in the impact of the fibrils' orientation. These results demonstrate that the collagen fibril orientation should be taken into account to properly describe the effective elastic anisotropy of bone at the organ scale.

An Instrumented Hammer to Detect the Rupture of the Pterygoid Plates.

PURPOSE: Craniofacial osteotomies involving pterygomaxillary disjunction are common procedures in maxillofacial surgery. Surgeons still rely on their proprioception to determine when to stop impacting the osteotome, which is important to avoid complications such as dental damage and bleeding. Our group has developed a technique consisting in using an instrumented hammer that can provide information on the mechanical properties of the tissue located around the osteotome tip. The aim of this study is to determine whether a mallet instrumented with a force sensor can be used to predict the crossing of the osteotome through the pterygoid plates. METHODS: 31 osteotomies were carried out in 16 lamb skulls. For each impact, the force signal obtained was analysed using a dedicated signal processing technique. A prediction algorithm based on an SVM classifier and a cost matrix was applied to the database. RESULTS: We showed that the device could always detect the crossing of the osteotome, sometimes before its occurrence. The prediction accuracy of the device was 94.7%. The method seemed to be sensitive to the thickness of the plate and to crack apparition and propagation. CONCLUSION: These results pave the way for the development of a per-operative decision support system in maxillofacial surgery.

Characterization of the concentration of agar-based soft tissue mimicking phantoms by impact analysis.

In various medical fields, a change of soft tissue stiffness is associated with its physio-pathological evolution. While elastography is extensively employed to assess soft tissue stiffness in vivo, its application requires a complex and expensive technology. The aim of this study is to determine whether an easy-to-use method based on impact analysis can be employed to determine the concentration of agar-based soft tissue mimicking phantoms. Impact analysis was performed on soft tissue mimicking phantoms made of agar gel with a mass concentration ranging from 1% to 5%. An indicator Δt is derived from the temporal variation of the impact force signal between the hammer and a small beam in contact with the sample. The results show a non-linear decrease of Δt as a function of the agar concentration (and thus of the sample stiffness). The value of Δt provides an estimation of the agar concentration with an error of 0.11%. This sensitivity of the impact analysis based method to the agar concentration is of the same order of magnitude than results obtained with elastography techniques. This study opens new paths towards the development of impact analysis for a fast, easy and relatively inexpensive clinical evaluation of soft tissue elastic properties.