Prediction of micromotion initiation of an implanted femur under physiological loads and constraints using the finite element method

Gear-shaped carbonate apatite granules with a hexagonal macropore for rapid bone regeneration

Synthetic bone grafts are in high demand owing to increased age-related bone disorders in the global aging population. Here, we report fabrication of gear-shaped granules (G-GRNs) for rapid bone healing. G-GRNs possessed six protrusions and a hexagonal macropore in the granular center. These were composed of carbonate apatite, i.e., bone mineral, microspheres with ∼1-μm micropores in the spaces between the microspheres. G-GRNs formed new bone and blood vessels (both on the granular surface and within the macropores) 4 weeks after implantation in the rabbit femur defects. The formed bone structure was similar to that of cancellous bone. The bone percentage in the defect recovered to that in a normal rabbit femur at week-4 post-implantation, and the bone percentage remained constant for the following 8 weeks. Throughout the entire period, the bone percentage in the G-GRN-implanted group was ∼10% higher than that of the group implanted with conventional carbonate apatite granules. Furthermore, a portion of the G-GRNs resorbed at week-4, and resorption continued for the following 8 weeks. Thus, G-GRNs are involved in bone remodeling and are gradually replaced with new bone while maintaining a suitable bone level. These findings provide a basis for the design and fabrication of synthetic bone grafts for achieving rapid bone regeneration.

A Review of Biomaterials and Associated Performance Metrics Analysis in Pre-Clinical Finite Element Model and in Implementation Stages for Total Hip Implant System

Total hip replacement (THR) is a common orthopedic surgery technique that helps thousands of individuals to live normal lives each year. A hip replacement replaces the shattered cartilage and bone with an implant. Most hip implants fail after 10–15 years. The material selection for the total hip implant systems is a major research field since it affects the mechanical and clinical performance of it. Stress shielding due to excessive contact stress, implant dislocation due to a large deformation, aseptic implant loosening due to the particle propagation of wear debris, decreased bone remodeling density due to the stress shielding, and adverse tissue responses due to material wear debris all contribute to the failure of hip implants. Recent research shows that pre-clinical computational finite element analysis (FEA) can be used to estimate four mechanical performance parameters of hip implants which are connected with distinct biomaterials: von Mises stress and deformation, micromotion, wear estimates, and implant fatigue. In vitro, in vivo, and clinical stages are utilized to determine the hip implant biocompatibility and the unfavorable local tissue reactions to different biomaterials during the implementation phase. This research summarizes and analyses the performance of the different biomaterials that are employed in total hip implant systems in the pre-clinical stage using FEA, as well as their performances in in vitro, in vivo, and in clinical studies, which will help researchers in gaining a better understanding of the prospects and challenges in this field.

Application of finite element analysis to tissue differentiation and bone remodelling approaches and their use in design optimization of orthopaedic implants: A review

Post-operative bone growth and long-term bone adaptation around the orthopaedic implants are simulated using the mechanoregulation based tissue-differentiation and adaptive bone remodelling algorithms, respectively. The primary objective of these algorithms was to assess biomechanical feasibility and reliability of orthopaedic implants. This article aims to offer a comprehensive review of the developments in mathematical models of tissue-differentiation and bone adaptation and their applications in studies involving design optimization of orthopaedic implants over three decades. Despite the different mechanoregulatory models developed, existing literature confirm that none of the models can be highly regarded or completely disregarded over each other. Not much development in mathematical formulations has been observed from the current state of knowledge due to the lack of in vivo studies involving clinically relevant animal models, which further retarded the development of such models to use in translational research at a fast pace. Future investigations involving artificial intelligence (AI), soft-computing techniques and combined tissue-differentiation and bone-adaptation studies involving animal subjects for model verification are needed to formulate more sophisticated mathematical models to enhance the accuracy of pre-clinical testing of orthopaedic implants.

Multiscale design of artificial bones with biomimetic elastic microstructures

Cancellous bone is a highly porous, heterogeneous, and anisotropic material which can be found at the epiphyses of long bones and in the vertebral bodies. The hierarchical architecture makes cancellous bone a prime example of a lightweight natural material that combines strength with toughness. Better understanding the mechanics of cancellous bone is of interest for the diagnosis of bone diseases, the evaluation of the risk of fracture, and for the design of artificial bones and bone scaffolds for tissue engineering. A multiscale optimization method to maximize the stiffness of artificial bones using biomimetic cellular microstructures described by a finite set of geometrical micro-parameters is presented here. The most outstanding characteristics of its implementation are the use of: an interior point optimization algorithm, a precalculated response surface methodology for the evaluation of the elastic tensor of the microstructure as an analytical function of the micro-parameters, and the adjoint method for the computation of the sensitivity of the macroscopic mechanical response to the variation of the micro-parameters. The performance and effectiveness of the tool are evaluated by solving a problem that consists in finding the optimal distribution of the microstructures for a proximal end of a femur subjected to physiological loads. Two strategies for the specification of the solid volume fraction constraints are assessed. The results are compared with data of a computed tomography study of an actual human bone. The model successfully predicts the main features of the spatial arrangement of the trabecular and cortical microstructures of the natural bone.

Impact of screw thread shape on stress transfer in bone: a finite element study

Gradual screw loosening is a well-known failure mechanism in internal fixation. Loosening is primarily due to progressive bone loss caused by stress shielding, a phenomenon in which a medical device absorbs a disproportionate amount of load within the screw-bone construct. The proximity of elastic moduli of magnesium and bone presents the potential for alleviating screw loosening by allowing optimum stress to be transferred between screw and bone, and in turn, supporting bone remodeling around the screw. In this study, the effect of thread profile on stress transfer in a magnesium fixation was simulated using a 2-D finite element model. Modified stress parameters from a previous study were used to estimate stress transfer across three thread profiles. Results showed highest stress transfer in trapezoidal-shaped magnesium screw thread. In accordance, this study corroborates the potential for magnesium as an ultimate screw material to eliminate progressive screw loosening.

Finite element analysis of the effect of medullary contact on fracture healing and remodeling in the intramedullary interlocking nail-fixed tibia fracture

Intramedullary interlocking nail is an effective treatment for tibial diaphyseal fracture. The contact between medullary rod and diaphyseal cortex is able to enhance fracture stability. However, how and to what degree the contact affects fracture healing and subsequent bone remodeling is still unclear. To investigate this, fracture healing and remodeling algorithms were combined, improved, and used to simulate the healing and remodeling processes in a transverse tibial diaphyseal fracture fixed with an intramedullary interlocking nail device. Two different diaphyseal fracture statuses, three different initial loading levels, and two nail materials were considered. The results showed that the medullary contact could significantly enhance the fixation stability; the strain reduction was up to 80% in the initial granulation callus. However, low initial loading level was found to be a more potential risk factor for the insufficient loading-induced nonunion other than medullary contact and stiffer nail material. Furthermore, the stabilizing effect of medullary contact diminished when stiff bone tissue formed in the callus; thus, the remodeling in the long-term was not affected by medullary contact. Copyright © 2016 John Wiley & Sons, Ltd.

Biomechanical determinants of the stability of dental implants: influence of the bone-implant interface properties

Dental implants are now widely used for the replacement of missing teeth in fully or partially edentulous patients and for cranial reconstructions. However, risks of failure, which may have dramatic consequences, are still experienced and remain difficult to anticipate. The stability of biomaterials inserted in bone tissue depends on multiscale phenomena of biomechanical (bone-implant interlocking) and of biological (mechanotransduction) natures. The objective of this review is to provide an overview of the biomechanical behavior of the bone-dental implant interface as a function of its environment by considering in silico, ex vivo and in vivo studies including animal models as well as clinical studies. The biomechanical determinants of osseointegration phenomena are related to bone remodeling in the vicinity of the implants (adaptation of the bone structure to accommodate the presence of a biomaterial). Aspects related to the description of the interface and to its space-time multiscale nature will first be reviewed. Then, the various approaches used in the literature to measure implant stability and the bone-implant interface properties in vitro and in vivo will be described. Quantitative ultrasound methods are promising because they are cheap, non invasive and because of their lower spatial resolution around the implant compared to other biomechanical approaches.

Optimal bone density distributions: numerical analysis of the osteocyte spatial influence in bone remodeling

In this paper a control and optimization procedure for bone remodeling simulations was adopted to study the effect of the osteocyte influence range on the predicted density distribution. In order to reach this goal, the osteocyte network regulating bone remodeling process in a 2-D bone sample was numerically simulated. The assumed proportional-integral-derivative (PID) bone remodeling rule was related to the error signal between the strain energy density and a selected target. Furthermore the control parameters and the target were optimally determined minimizing a suitable cost index: the goal was to minimize the final mass and the energy thus maximizing the stiffness. The continuum model results show that the developed and adapted trabecular structure was consistent with the applied loads and only depended on the external forces, the value of the cost index, the maximum attainable elastic modulus value (hence, the maximum density value) and the value of the energy target. The remodeling phenomenon determined the number and thickness of the trabeculae which are formed from a uniform distribution of mass density in the considered domain; this number and these thicknesses are controlled by the values assigned to the parameters of the model. In particular, the osteocyte decay distance (D) of the influence range affected the trabecular patterns formation, showing an important effect in the adaptive capacity of the optimization numerical model.

Finite element analysis modelling of proximal femoral fractures, including post-fixation periprosthetic fractures

Biomechanical testing has been a cornerstone for the development of surgical implants used in fracture stabilisation. In a multi-disciplinary collaboration complex at the University of Wales, Swansea, novel computerised clinically relevant models were developed using advanced computational engineering. In-house software (developed initially for commercial aerospace engineering), allowed accurate finite element analysis (FEA) models of the whole femur to be created, including the internal architecture of the bone, by means of linear interpolation of greyscale images from multiaxial CT scans. This allowed for modelling the changing trabecular structure and bone mineral density as seen in progressive osteoporosis. Falls from standing were modelled in a variety of directions (with and without muscle action) using analysis programmes which resulted in fractures consistent with those seen in clinical practice. By meshing implants into these models and repeating the mechanism of injury in simulation, periprosthetic fractures were also recreated. Further development with simulated physiological activities (e.g. walking and rising from sitting) along with attrition in the bone (in the boundary zones where stress concentration occurs) will allow further known modes of failure in implants to be reproduced. Robust simulation of macro and micro-scale events will allow the testing of novel new designs in simulations far more complex than conventional biomechanical testing will allow.

Real-time controller for foot-drop correction by using surface electromyography sensor

Foot drop is a disease caused mainly by muscle paralysis, which incapacitates the nerves generating the impulses that control feet in a heel strike. The incapacity may stem from lesions that affect the brain, the spinal cord, or peripheral nerves. The foot becomes dorsiflexed, affecting normal walking. A design and analysis of a controller for such legs is the subject of this article. Surface electromyography electrodes are connected to the skin surface of the human muscle and work on the mechanics of human muscle contraction. The design uses real surface electromyography signals for estimation of the joint angles. Various-speed flexions and extensions of the leg were analyzed. The two phases of the design began with surface electromyography of real human leg electromyography signal, which was subsequently filtered, amplified, and normalized to the maximum amplitude. Parameters extracted from the surface electromyography signal were then used to train an artificial neural network for prediction of the joint angle. The artificial neural network design included various-speed identification of the electromyography signal and estimation of the angles of the knee and ankle joints by a recognition process that depended on the parameters of the real surface electromyography signal measured through real movements. The second phase used artificial neural network estimation of the control signal, for calculation of the electromyography signal to be stimulated for the leg muscle to move the ankle joint. Satisfactory simulation (MATLAB/Simulink version 2012a) and implementation results verified the design feasibility.

The effect of cup orientation and coverage on contact mechanics and range of motion of metal-on-metal hip resurfacing arthroplasty

Implant malpositioning has been identified as a factor associated with clinical failures of metal-on-metal hip resurfacings (MoMHRs). This study investigated the effect of cup orientation and cup coverage on the contact mechanics (incidence of edge-loading) and range of motion (ROM) of MoMHR. Three generic MoMHRs with differing amounts of cup coverage were considered at various orientations. Contact area and contact pressure at the bearing surface were predicted for each design using finite element (FE) method. The ROM was determined based on the geometry overlap. Edge contact was found at lower angles of inclination (65°) for lower coverage cup designs; however, they also provided the greatest ROM. Conversely, cups with greater coverage did not exhibit edge contact until the cup was more steeply positioned (75°), however ROM was reduced. This study enables both sets of variable to be considered in the design of metal-on-metal bearings in hip.

Experimental evaluation of new concepts in hip arthroplasty

In this thesis we evaluated two different hip arthroplasty concepts trough in vitro studies and numerical analyses. The cortical strains in the femoral neck area were increased by 10 to 15 % after insertion of a resurfacing femoral component compared to values of the intact femur, shown in an in vitro study on human cadaver femurs. There is an increased risk of femoral neck fracture after hip resurfacing arthroplasty. An increase of 10 to 15 % in femoral neck strains is limited, and cannot alone explain these fractures. Together with patient specific and surgical factors, however, increased strain can contribute to increased risk of fracture. An in vitro study showed that increasing the neck length in combination with retroversion or reduced neck shaft angle on a standard cementless femoral stem does not compromise the stability of the stem. The strain pattern in the proximal femur increased significantly at several measuring sites when the version and length of neck were altered. However, the changes were probably too small to have clinical relevance. In a validation study we have shown that a subject specific finite element analysis is able to perform reasonable predictions of strains and stress shielding after insertion of a femoral stem in human cadaver femurs. The usage of finite element models can be a valuable supplement to in vitro tests of femoral strain pattern around hip arthroplasty. Finally, a patient case shows that bone resorption around an implant caused by stress shielding can in extreme cases lead to periprosthetic fracture.