Numerical Simulation of an Intramedullary Elastic Nail: Expansion Phase and Load-Bearing Behavior

Determination of intramedullary nail based on centerline adaptive registration

Objective: Internal fixation with intramedullary nails is a gold standard for the treatment of femoral shaft fractures. However, both the mismatch between intramedullary nails and the medullary cavity and inaccurate positioning of entry points will lead to deformation of intramedullary nails after implantation. The study aimed to determine a suitable intramedullary nail with an optimal entry point for a specific patient based on centerline adaptive registration.

Method: A homotopic thinning algorithm is employed to extract centerlines of the femoral medullary cavity and the intramedullary nail. The two centerlines are registered to obtain a transformation. The medullary cavity and the intramedullary nail are registered based on the transformation. Next, a plane projection method is employed to calculate the surface points of the intramedullary nail laid outside the medullary cavity. According to the distribution of compenetration points, an iterative adaptive registration strategy is designed to decide an optimal position of the intramedullary nail in medullary cavity. The isthmus centerline is extended to the femur surface, where the entry point of the intramedullary nail is located. The suitability of an intramedullary nail for a specific patient was calculated by measuring the geometric quantities reflecting the interference between the femur and nail, and the suitability values of all nails are compared and the most suitable one is determined.

Results: The growth experiment indicated that the bone to nail alignment is indeed affected by the extension of the isthmus centerline, including the extension direction and velocity. The geometrical experiment showed that this method could find the best registration position of intramedullary nails and select the optimal intramedullary nail for a specific patient. In the model experiments, the determined intramedullary nail could be successfully placed into the medullary cavity through the optimal entry point. A pre-screening tool to determine nails which can be successfully used has been given. In addition, the distal hole was accurately located within 14.28 s.

Conclusion: These results suggest that the proposed method can select a suitable intramedullary nail with an optimal entry point. The position of the intramedullary nail can be determined in the medullary cavity, while deformation is avoided. The proposed method can determine the largest diameter intramedullary nail with as little damage to the intramedullary tissue as possible. The proposed method provides preparation aid for internal fixation with intramedullary nails guided by navigation systems or extracorporeal aimers.

Biomechanical Analysis of a Novel Double-Point Fixation Method for Displaced Intra-Articular Calcaneal Fractures

The development of minimally invasive procedures and implant materials has improved the fixation strength of implants and is less traumatic in surgery. The purpose of this study was to propose a novel “double-point fixation” for calcaneal fractures and compare its biomechanical stability with the traditional “three-point fixation.” A three-dimensional finite element foot model with a Sanders type IIIAB calcaneal fracture was developed based on clinical images comprising bones, plantar fascia, ligaments, and encapsulated soft tissue. Double-point and three-point fixation resembled the surgical procedure with a volar distal radius plate and calcaneal locking plate, respectively. The stress distribution, fracture displacement, and change of the Böhler angle and Gissane’s angle were estimated by a walking simulation using the model, and the predictions between the double-point and three-point fixation were compared at heel-strike, midstance, and push-off instants. Double-point fixation demonstrated lower bone stress (103.3 vs. 199.4 MPa), but higher implant stress (1,084.0 vs. 577.9 MPa). The model displacement of double-point fixation was higher than that of three-point fixation (3.68 vs. 2.53 mm). The displacement of the posterior joint facet (0.127 vs. 0.150 mm) and the changes of the Böhler angle (0.9° vs. 1.4°) and Gissane’s angle (0.7° vs. 0.9°) in double-point fixation were comparably lower. Double-point fixation by volar distal radius plates demonstrated sufficient and favorable fixation stability and a lower risk of postoperative stress fracture, which may potentially serve as a new fixation modality for the treatment of displaced intra-articular calcaneal fractures.

Stochastic PCA-Based Bone Models from Inverse Transform Sampling: Proof of Concept for Mandibles and Proximal Femurs

Principal components analysis is a powerful technique which can be used to reduce data dimensionality. With reference to three-dimensional bone shape models, it can be used to generate an unlimited number of models, defined by thousands of nodes, from a limited (less than twenty) number of scalars. The full procedure has been here described in detail and tested. Two databases were used as input data: the first database comprised 40 mandibles, while the second one comprised 98 proximal femurs. The “average shape” and principal components that were required to cover at least 90% of the whole variance were identified for both bones, as well as the statistical distributions of the respective principal components weights. Fifteen principal components sufficed to describe the mandibular shape, while nine components sufficed to describe the proximal femur morphology. A routine has been set up to generate any number of mandible or proximal femur geometries, according to the actual statistical shape distributions. The set-up procedure can be generalized to any bone shape given a sufficiently large database of the respective 3D shapes.

Mechanical Behavior of Elastic Self-Locking Nails for Intramedullary Fracture Fixation: A Numerical Analysis of Innovative Nail Designs

Intramedullary nails constitute a viable alternative to extramedullary fixation devices; their use is growing in recent years, especially with reference to self-locking nails. Different designs are available, and it is not trivial to foresee the respective in vivo performances and to provide clinical indications in relation to the type of bone and fracture. In this work a numerical methodology was set up and validated in order to compare the mechanical behavior of two new nailing device concepts with one already used in clinic. In detail, three different nails were studied: (1) the Marchetti-Vicenzi's nail (MV₁), (2) a revised concept of this device (MV₂), and (3) a new Terzini-Putame's nail (TP) concept. Firstly, the mechanical behavior of the MV₁ device was assessed through experimental loading tests employing a 3D-printed component aimed at reproducing the bone geometry inside which the device is implanted. In the next step, the respective numerical model was created, based on a multibody approach including flexible parts, and this model was validated against the previously obtained experimental results. Finally, numerical models of the MV₂ and TP concepts were implemented and compared with the MV₁ nail, focusing the attention on the response of all devices to compression, tension, bending, and torsion. A stability index (SI) was defined to quantify the mechanical stability provided to the nail-bone assembly by the elastic self-locking mechanism for the various loading conditions. In addition, results in terms of nail-bone assembly stiffness, computed from force/moment vs. displacement/rotation curves, were presented and discussed. Findings revealed that numerical models were able to provide good estimates of load vs. displacement curves. The TP nail concept proved to be able to generate a significantly higher SI (27 N for MV₁ vs. 380 N for TP) and a greater stiffening action (up to a stiffness difference for bending load that ranges from 370 Nmm/° for MV₁ to 1,532 Nmm/° for TP) than the other two devices which showed similar performances. On the whole, a demonstration was given of information which can be obtained from numerical simulations of expandable fixation devices.

Prosthetic Hip ROM from Multibody Software Simulation

The pre-operative planning of a hip arthroplasty entails the choice of the prosthetic hip model and of the position of both joint components with reference to bone. Assessing the impact of geometrical factors on the final hip range of motion (ROM) is not trivial, since it requires performing 3D evaluations. Nonetheless, it deserves to be studied since hip impingement and dislocation are still relevant complications in hip arthroplasty. This work pertains a numerical model for the assessment of the hip ROM in relation to cotyle position. External/internal rotation is considered as a benchmark, and multiple combinations of acetabular anteversion/inclination are considered. According to results, over two hundred different geometric configurations can be examined in few minutes, and the cotyle position can be so optimized with relevant benefits in term of hip ROM.

Failure Analysis of Retrieved Osteosynthesis Implants

Failure of osteosynthesis implants is an intricate matter with challenging management that calls for efficient investigation and prevention. Using implant retrieval analysis combined with standard radiological examination, we evaluated the main causes for osteosynthesis implant breakdown and the relations among them for a series of cases. Twenty-one patients diagnosed with implant failure were assessed for this work. For metallurgical analysis, microscopy techniques such as scanning electron microscopy (SEM), stereomicroscopy, and optical microscopy were employed. The results showed that material structural deficiencies (nine patients) and faulty surgical techniques (eight patients) were the main causes for failure. An important number of patients presented with material structural deficiencies superimposed on an imperfect osteosynthesis technique (six patients). Consequently, the importance of failure retrieval analysis should not be overlooked, and in combination with other investigational techniques, must provide information for both implant manufacturing and design improvement, as well as osteosynthesis technique optimization.

Determining the Online Measurable Input Variables in Human Joint Moment Intelligent Prediction Based on the Hill Muscle Model

Introduction: Human joint moment is a critical parameter to rehabilitation assessment and human-robot interaction, which can be predicted using an artificial neural network (ANN) model. However, challenge remains as lack of an effective approach to determining the input variables for the ANN model in joint moment prediction, which determines the number of input sensors and the complexity of prediction. Methods: To address this research gap, this study develops a mathematical model based on the Hill muscle model to determining the online input variables of the ANN for the prediction of joint moments. In this method, the muscle activation, muscle-tendon moment velocity and length in the Hill muscle model and muscle-tendon moment arm are translated to the online measurable variables, i.e. muscle electromyography (EMG), joint angles and angular velocities of the muscle span. To test the predictive ability of these input variables, an ANN model is designed and trained to predict joint moments. The ANN model with the online measurable input variables is tested on the experimental data collected from ten healthy subjects running with the speeds of 2, 3, 4 and 5 m/s on a treadmill. The variance accounted for (VAF) between the predicted and inverse dynamics moment is used to evaluate the prediction accuracy. Results: The results suggested that the method can predict joint moments with a higher accuracy (mean VAF = 89.67±5.56 %) than those obtained by using other joint angles and angular velocities as inputs (mean VAF = 86.27±6.6%) evaluated by jack-knife cross-validation. Conclusions: The proposed method provides us with a powerful tool to predict joint moment based on online measurable variables, which establishes the theoretical basis for optimizing the input sensors and detection complexity of the prediction system. It may facilitate the research on exoskeleton robot control and real-time gait analysis in motor rehabilitation.

Multibody Models for the Analysis of a Fall From Height: Accident, Suicide, or Murder?

The final subject position is often the only evidence in the case of the fall of a human being from a given height. Foreseeing the body trajectory and the respective driving force may not be trivial due to the possibility of rotations and to an unknown initial position and momentum of the subject. This article illustrates how multibody models can be used for this aim, with specific reference to an actual case, where a worker fell into a stair well, prior to stair mounting, and he was found in an unexpected posture. The aim of the analysis was establishing if this worker was dead in that same place, if he had been pushed, and which was his initial position. A multibody model of the subject has been built ("numerical android"), given his stature and his known mass. Multiple simulations have been performed, following a design of experiments where various initial positions and velocity as well as pushing forces have been considered, while the objective function to be minimized was the deviation of the numerical android position from the actual worker position. At the end of the analysis, it was possible to point how a very limited set of conditions, all including the application of an external pushing force (or initial speed), could produce the given final posture with an error on the distance function equal to 0.39 m. The full analysis gives a demonstration of the potentiality of multibody models as a tool for the analysis of falls in forensic inquiries.