Finite element analysis of a hemi-pelvis: the effect of inclusion of cartilage layer on acetabular stresses and strain

Assessment of bone ingrowth around beaded coated tibial implant for total ankle replacement using mechanoregulatory algorithm

The long-term performance of porous coated tibial implants for total ankle replacement (TAR) primarily depends on the extent of bone ingrowth at the bone-implant interface. Although attempts were made for primary fixation for immediate post-operative stability, no investigation was conducted on secondary fixation. The aim of this study is to assess bone ingrowth around the porous beaded coated tibial implant for TAR using a mechanoregulatory algorithm. A realistic macroscale finite element (FE) model of the implanted tibia was developed based on computer tomography (CT) data to assess implant-bone micromotions and coupled with microscale FE models of the implant-bone interface to predict bone ingrowth around tibial implant for TAR. The macroscale FE model was subjected to three near physiological loading conditions to evaluate the site-specific implant-bone micromotion, which were then incorporated into the corresponding microscale model to mimic the near physiological loading conditions. Results of the study demonstrated that the implant experienced tangential micromotion ranged from 0 to 71 μm with a mean of 3.871 μm. Tissue differentiation results revealed that bone ingrowth across the implant ranged from 44 to 96 %, with a mean of around 70 %. The average Young's modulus of the inter-bead tissue layer varied from 1444 to 4180 MPa around the different regions of the implant. The analysis postulates that when peak micromotion touches 30 μm around different regions of the implant, it leads to pronounced fibrous tissues on the implant surface. The highest amount of bone ingrowth was observed in the central regions, and poor bone ingrowth was seen in the anterior parts of the implant, which indicate improper osseointegration around this region. This macro-micro mechanical FE framework can be extended to improve the implant design to enhance the bone ingrowth and in future to develop porous lattice-structured implants to predict and enhance osseointegration around the implant.

A macro-micro FE and ANN framework to assess site-specific bone ingrowth around the porous beaded-coated implant: an example with BOX® tibial implant for total ankle replacement

The use of mechanoregulatory schemes based on finite element (FE) analysis for the evaluation of bone ingrowth around porous surfaces is a viable approach but requires significant computational time and effort. The aim of this study is to develop a combined macro-micro FE and artificial neural network (ANN) framework for rapid and accurate prediction of the site-specific bone ingrowth around the porous beaded-coated tibial implant for total ankle replacement (TAR). A macroscale FE model of the implanted tibia was developed based on CT data. Subsequently, a microscale FE model of the implant-bone interface was created for performing bone ingrowth simulations using mechanoregulatory algorithms. An ANN was trained for rapid and accurate prediction of bone ingrowth. The results predicted by ANN are well comparable to FE-predicted results. Predicted site-specific bone ingrowth using ANN around the implant ranges from 43.04 to 98.24%, with a mean bone ingrowth of around 74.24%. Results suggested that the central region exhibited the highest bone ingrowth, which is also well corroborated with the recent explanted study on BOX®. The proposed methodology has the potential to simulate bone ingrowth rapidly and effectively at any given site over any implant surface.

A numerical investigation for the development of functionally graded Ti/HA tibial implant for total ankle replacement: Influence of material gradation law and volume fraction index

Stress shielding is one of the major concerns for total ankle replacement implants nowadays, because it is responsible for implant-induced bone resorption. The bone resorption contributes to the aseptic loosening and failure of ankle implants in later stages. To reduce the stress shielding, improvements can be made in the implant material by decreasing the elastic mismatch between the implant and the tibia bone. This study proposes a new functionally graded material (FGM) based tibial implant for minimizing the problem of stress shielding. Three-dimensional finite element (FE) models of the intact tibia and the implanted tibiae were created to study the influence of material gradation law and volume fraction index on stress shielding and implant-bone micromotion. Different implant materials were considered that is, cobalt-chromium, titanium (Ti), and FGM with Ti at the bottom and hydroxyapatite (HA) at the top. The FE models of FGM implants were generated by using different volume fractions and the rule of mixtures. The rule of mixtures was used to calculate the FGM properties based on the local volume fraction. The volume fraction was defined by using exponential, power, and sigmoid laws. For the power and sigmoid law varying volume fraction indices (0.1, 0.2, 0.5, 1, 2, and 5) were considered. The geometry resembling STAR® ankle system tibial implant was considered for the present study. The results indicate that FGMs lower stress shielding but also marginally increase implant-bone micromotion; however, the values were within the acceptable limit for bone ingrowth. It is observed that the material gradation law and volume fraction index influence the performance of FGM tibial implants. The tibial implant composed of FGM using power law with a volume fraction index of 0.1 was the preferred option because it showed the least stress shielding.

Implementing Machine Learning approaches for accelerated prediction of bone strain in acetabulum of a hip joint

The Finite Element (FE) methods for biomechanical analysis involving implant design and subject parameters for musculoskeletal applications are extensively reported in literature. Such an approach is manually intensive and computationally expensive with longer simulations times. Although Artificial Intelligence (AI) based approaches are implemented to a limited extent in biomechanics, such approaches to predict bone strain in acetabulum of a hip joint, are hardly explored. In this context, the primary objective of this paper is to evaluate machine learning (ML) models in tandem with high-fidelity FEA data for the accelerated prediction of the biomechanical response in the acetabulum of the human hip joint, during the walking gait. The parameters used in the FEA study included the subject weight, number and distribution of fins on the periphery of the acetabular shell, bone condition and phases of the gait cycle. The biomechanical response has also been evaluated using three different acetabular liners, including pre-clinically validated HDPE-20% HA-20% Al2O3, highly-crosslinked ultrahigh molecular weight polyethylene (HC-UHMWPE) and ZrO2-toughened Al2O3 (ZTA). Such parametric variation in FEA analysis, involving 26 variables and a full factorial design resulted in 10,752 datasets for spatially varying bone strains. The bone condition, as opposed to subject weight, was found to play a statistically significant role in determining the strain response in the periprosthetic bone of the acetabulum. While utilising hyperparameter tuning, K-fold cross validation and statistical learning approaches, a number of ML models were trained on the FEA dataset, and the Random Forest model performed the best with a coefficient of determination (R2) value of 0.99/0.97 and Root Mean Square Error (RMSE) of 0.02/0.01 on the training/test dataset. Taken together, this study establishes the potential of ML approach as a fast surrogate of FEA for implant biomechanics analysis, in less than a minute.

Towards an optimal design of a functionally graded porous uncemented acetabular component using genetic algorithm

Generation of polyethylene wear debris and peri‑prosthetic bone resorption have been identified as potential causes of acetabular component loosening in Total Hip Arthroplasty. This study was aimed at optimization of a functionally graded porous acetabular component to minimize peri‑prosthetic bone resorption and polyethylene liner wear. Porosity levels (porosity values at acetabular rim, and dome) and functional gradation exponents (radial and polar) were considered as the design parameters. The relationship between porosity and elastic properties were obtained from numerical homogenization. The multi-objective optimization was carried out using a non-dominated sorting genetic algorithm integrated with finite element analysis of the hemipelvises subject to various loading conditions of common daily activities. The optimal functionally graded porous designs (OFGPs -1, -2, -3, -4, -5) exhibited less strain-shielding in cancellous bone compared to solid metal-backing. Maximum bone-implant interfacial micromotions (63-68 μm) for OFGPs were found to be close to that of solid metal-backing (66 μm), which might facilitate bone ingrowth. However, OFGPs exhibited an increase in volumetric wear (3-10 %) compared to solid metal-backing. The objective functions were found to be more sensitive to changes in polar gradation exponent than radial gradation exponent, based on the Sobol' method. Considering the common failure mechanisms, OFGP-1, having highly porous acetabular rim and less porous dome, appears to be a better alternative to the solid metal-backing.

Biomechanical Evaluation of 2 Endoscopic Spine Surgery Methods for Treating Lumbar Disc Herniation: A Finite Element Study

OBJECTIVE

This study aimed to evaluate the effects of 2 endoscopic spine surgeries on the biomechanical properties of normal and osteoporotic spines.

METHODS

Based on computed tomography images of a healthy adult volunteer, 6 finite element models were created. After validating the normal intact model, a concentrated force of 400 N and a moment of 7.5 Nm were exerted on the upper surface of L3 to simulate 6 physiological activities of the spine. Five types of indices were used to assess the biomechanical properties of the 6 models, range of motion (ROM), maximum displacement value, intervertebral disc stress, maximum stress value, and articular protrusion stress, and by combining them with finite element stress cloud.

RESULTS

In normal and osteoporotic spines, there was no meaningful change in ROM or disc stress in the 2 surgical models for the 6 motion states. Model N1 (osteoporotic percutaneous transforaminal endoscopic discectomy model) showed a decrease in maximum displacement value of 20.28% in right lateral bending. Model M2 (unilateral biportal endoscopic model) increased maximum displacement values of 16.88% and 17.82% during left and right lateral bending, respectively. The maximum stress value of L4-5 increased by 11.72% for model M2 during left rotation. In addition, using the same surgical approach, ROM, maximum displacement values, disc stress, and maximum stress values were more significant in the osteoporotic model than in the normal model.

CONCLUSION

In both normal and osteoporotic spines, both surgical approaches were less disruptive to the physiologic structure of the spine. Furthermore, using the same endoscopic spine surgery, normal spine biomechanical properties are superior to osteoporotic spines.

A combined FE-hybrid MCDM framework for improving the performance of the conical stem tibial design for TAR with the addition of pegs

BACKGROUND AND OBJECTIVES

The conical stemmed design of the tibial component for total ankle replacement (TAR) (example Mobility design) uses a single intramedullary stem for primary fixation. Tibial component loosening is a common mode of failure for TAR. Primary causes of loosening are lack of bone ingrowth due to excessive micromotion at the implant-bone interface and bone resorption due to stress shielding after implantation. The fixation feature of the conical stemmed design can be modified with the addition of small pegs to avoid loosening. The aim of the study is to select the improved design for conical stemmed TAR using a combined Finite Element (FE) hybrid Multi-Criteria Decision-Making (MCDM) framework.

METHODS

The geometry and material properties of the bone for FE modeling were extracted from the CT data. Thirty-two design alternatives with varying pegs in number (one, two, four, eight), location (anterior, posterior, medial, lateral, anterior-posterior, medial-lateral, equally spaced), and height (5 mm, 4 mm, 3 mm, 2 mm) were prepared. All models were analyzed for dorsiflexion, neutral, and plantarflexion loading. The proximal part of the tibia was fixed. The implant-bone interface coefficient of friction was taken as 0.5. The implant-bone micromotion, stress shielding, volume of bone resection, and surgical simplicity were the important criteria considered for evaluating the performance of TAR. The designs were compared using a hybrid MCDM method of WASPAS, TOPSIS, EDAS, and VIKOR. The weight calculations were based on fuzzy AHP and the final ranks on the Degree of Membership method.

RESULTS

The addition of pegs decreased the mean implant-bone micromotions and increased stress shielding. There was a marginal decrease in micromotion and a marginal increase in stress shielding when the peg heights were increased. The results of hybrid MCDM indicated that the most preferable alternative designs were two pegs of 4 mm height in the AP direction to the main stem, two pegs of 4 mm height in the ML direction, and one peg of 3 mm height in the A direction.

CONCLUSIONS

Outcomes of this study suggest that the addition of pegs can reduce the implant-bone micromotions. Modified three designs would be useful by considering implant-bone micromotions, stress shielding, volume of bone resection, and surgical simplicity.

Biomechanical performance of the cemented acetabular cup with combined effects of bone quality, implant material combinations and bodyweight

The objective of this study is to understand the combined effects of bone quality, implant materials and bodyweight on the biomechanical performance of cemented acetabular cup. Additionally, the performance of the cemented acetabular cup was evaluated for obesity cases or obese people. A total of 84 FE models (based on CT data) were developed based on combinations of three different cancellous bone material distributions to represent bone quality, four different implant material combinations and seven different bodyweights. The biomechanical performance of the acetabular cup was evaluated based on bone stress (both cortical and cancellous bone), cement mantle stress, micromotion and contact pressure between the acetabular cup and femoral head. Cortical bone stress, cancellous bone stress, cement stress, the contact pressure between implants and micromotion between implants are affected by different bone quality, implant material combinations and bodyweights. An increase in bodyweight would increase the cortical bone stress, cancellous bone stress, cement stress, contact pressure between implants and micromotion between implants. However, bodyweight affects the cortical and cancellous bone stress more (stiff rise of the bone stresses; nonlinear relation) as compared to other output parameters (mostly linear relation). Comparing cortical and cancellous bone stress, the stress versus bodyweight curve is much stiffer (stiff rise in the curve) for cortical bone than cancellous bone and that even further increases as bone quality decreases. Especially considering obesity cases or obese people (very high bodyweight), the performance of the cemented acetabular component is poor.

Graphical abstract

[Formula: see text]

Computer-aided automatic planning and biomechanical analysis of a novel arc screw for pelvic fracture internal fixation

BACKGROUND AND OBJECTIVE

The sacroiliac joint screw is a common fixation method for pelvic posterior ring fractures. The complex anatomical structure around the pelvis makes it impossible to find a suitable fixed path, which increases the difficulty of surgical operation. In this paper, we propose an automatic planning algorithm based on a computer-aided internal arc fixation channel for pelvic fractures for the first time.

METHODS

A channel generation algorithm based on seed derived points was designed, and the optimal channel was selected by scoring rules based on 3D erode algorithm for the generated channel. The biomechanical properties of the internal arc fixation screw and traditional internal straight fixation screw in three postures were compared using biomechanical finite element analysis.

RESULTS

The proposed algorithm verified the existence of a more adaptable internal arc fixation channel and can quantitatively plan a relatively optimal constant-curvature internal arc fixation channel in pelvises of ten adults. Significantly high stresses concentrated around the interaction region between the screws and bone may increase the risk of bone fractures and screw loosening in the long term. The experimental results show that the internal arc fixation screw has better strain and deformation performance than the internal straight fixation screw.

CONCLUSIONS

A novel arc internal fixation method for pelvic fractures was proposed to improve the safety and stability of screw fixation of pelvic fracture. The nonparametric test proved that the sacroiliac dislocation model repaired by internal arc fixation screw was significantly different from that repaired by internal straight fixation screw. The computer-aided automatic planning algorithm provides the possibility of robot-assisted pelvic fracture fixation.

Role of the Anterior Center-Edge Angle on Acetabular Stress Distribution in Borderline Development Dysplastic of Hip Determined by Finite Element Analysis

Background: The joint with hip dysplasia is more likely to develop osteoarthritis because of the higher contact pressure, especially in the socket. The lateral center-edge angle (LCEA) is the major indicator for hip dysplasia via radiography. However, the pathological conditions of LCEA angles in the range of 18°–25° are still controversial, which challenges precise diagnosis and treatment decision-making.

Objective: The purpose of this study is to investigate the influence of anterior center-edge angle (ACEA) on the mechanical stress distribution of the hip joint, via finite element analysis, to provide insights into the severity of the borderline development dysplasia.

Methods: From 2017 to 2019, there were 116 patients with borderline developmental dysplasia of the hip (BDDH) enrolled in this research. Based on the inclusion criteria, nine patients were involved and categorized into three LCEA groups with the maximal ACEA differences. Patient-specific hip joint models were reconstructed from computed tomography scans, and the cartilages, including the labrum, were established via a modified numerical method. The finite element analysis was conducted to compare the stress distributions due to the different ACEA.

Results: As ACEA decreased, the maximum stress of the acetabulum increased, and the high stress area developed toward the edge. Quantitative analysis showed that in the cases with lower ACEA, the area ratio of high stress increased, and the contact facies lunata area significantly affected the stress distribution.

Conclusion: For patients with BDDH, both the ACEA and the area of facies lunata played essential roles in determining the severity of hip dysplasia and the mechanical mechanism preceding osteoarthritis.

Prediction of fracture initiation and propagation in pelvic bones

The objective is developing an XFEM model that is capable of predicting different types of fracture in the pelvic bone under various loading conditions. Previously published mechanical and failure characteristics of cortical and cancellous tissues were implemented and assigned to an intact pelvic bone with specified cortical and cancellous tissues. Various loading conditions, including combined load directions, were applied to the acetabulum to model different types of fracture (e.g., anterior/posterior wall fracture and transverse fracture) in the pelvic bone. The predicated types of fracture and the maximum force at fracture were compared to those acquired from previously published experimental tests. Anterior/posterior wall fracture and transverse fracture were the most common types of fractures determined in the simulations. The XFEM simulations were able to predict similar fractures to those reported in the experimental tests. The maximum fracture force in the XFEM model was found to be 18.6 kN compared to 8.85 kN reported in the previous experimental tests. The results revealed that different types of fracture in the pelvic bones can be caused by the various loading conditions in unstable high-rate impact loads. Using proper mechanical and failure behaviors of cortical and cancellous tissues, XFEM modeling of pelvic bone is capable of predicting bone fracture. In future work, the XFEM models of cancellous and cortical tissues can be assigned to other bones in human body skeleton so that the failure mechanism in such bones can be investigated.

Heterogeneous material mapping methods for patient-specific finite element models of pelvic trabecular bone: A convergence study

Patient-specific finite element (FE) models of bone require the assignment of heterogeneous material properties extracted from the subject's computed tomography (CT) images. Though node-based (NB) and element-based (EB) material mapping methods (MMMs) have been proposed, the sensitivity and convergence of FE models to MMM for varying mesh sizes are not well understood. In this work, CT-derived and synthetic bone material data were used to evaluate the effect of MMM on results from FE analyses. Pelvic trabecular bone data was extracted from CT images of six subjects, while synthetic data were created to resemble trabecular bone properties. The numerical convergence of FE bone models using different MMMs were evaluated for strain energy, von-Mises stress, and strain. NB and EB MMMs both demonstrated good convergence regarding total strain energy, with the EB method having a slight edge over the NB. However, at the local level (e.g., maximum stress and strain), FE results were sensitive to the field type, MMM, and the FE mesh size. The EB method exhibited superior performance in finer meshes relative to the voxel size. The NB method converged better than did the EB method for coarser meshes. These findings may lead to higher-fidelity patient-specific FE bone models.

The Impact of an Open-Book Pelvic Ring Injury on Bone Strain: Validation of a Finite Element Model and Analysis Within the Gait Cycle

The threshold for surgical stabilization for an open-book pelvic fracture is not well defined. The purpose of this research was to validate the biomechanical behavior of a specimen-specific pelvic finite element (FE) model with an open-book fracture with the biomechanical behavior of a cadaveric pelvis in double leg stance configuration under physiologic loading, and to utilize the validated model to compare open book versus intact strain patterns during gait. A cadaveric pelvis was experimentally tested under compressive loading in double leg stance, intact, and with a simulated open-book fracture. An intact FE model of this specimen was reanalyzed with an equivalent simulated open-book fracture. Comparison of the FE generated and experimentally measured strains yielded an R2 value of 0.92 for the open-book fracture configuration. Strain patterns in the intact and fractured models were compared throughout the gait cycle. In double leg stance and heel-strike/heel-off models, tensile strains decreased, especially in the pubic ramus contralateral to the injury, and compressive strains increased in the sacroiliac region of the injured side. In the midstance/midswing gait configuration, higher tensile and compressive FE strains were observed on the midstance side of the fractured versus intact model and decreased along the superior and inferior pubic rami and ischium, with midswing side strains reduced almost to zero in the fractured model. Identified in silico patterns align with clinical understanding of open-book fracture pathology suggesting future potential of FE models to quantify instability and optimize fixation strategies.

Numerical analysis of the mechanical behaviour of intact and implanted lumbar functional spinal units: Effects of loading and boundary conditions

The objective of this study was to develop an improved finite element (FE) model of a lumbar functional spinal unit (FSU) and to subsequently analyse the deviations in load transfer owing to implantation. The effects of loading and boundary conditions on load transfer in intact and implanted FSUs and its relationship with the potential risk of vertebral fracture were investigated. The FE models of L1-L5 and L3-L4 FSUs, intact and implanted, were developed using patient-specific CT-scan dataset and segmentation of cortical and cancellous bone regions. The effect of submodelling technique, as compared to artificial boundary conditions, on the elastic behaviour of lumbar spine was examined. Applied forces and moments, corresponding to physiologic movements, were used as loading conditions. Results indicated that the loading and boundary conditions considerably affect stress-strain distributions within a FSU. This study, based on an improved FE model of a vertebra, highlights the importance of using the submodelling technique to adequately evaluate the mechanical behaviour of a FSU. In the intact FSU, strains of 200-400 µε were observed in the cancellous bone of vertebral body and pedicles. High equivalent stresses of 10-25 MPa and 1-5 MPa were generated around the pars interarticularis for cortical and cancellous regions, respectively. Implantation caused reductions of 85%-92% in the range of motion for all movements. Insertion of the intervertebral cage resulted in major deviations in load transfer across a FSU for all movements. The cancellous bone around cage experienced pronounced increase in stresses of 10-15 MPa, which indicated potential risk of failure initiation in the vertebra.

The elastoplastic numerical model and verification by macroindentation experiment of femoral head

For internal fixation of proximal femoral fractures, a screw is commonly placed into the femoral head; therefore, mechanical matching of the femoral head and screw is important. This article proposes an elastoplastic numerical model of the femoral head that takes nonlinear deformation and cancellous bone heterogeneity into account. Force-depth curves from finite element analysis based on the model were compared with those from macroindentation experiments. The maximum difference between the indentation depth shown by the finite element model and that found with macroindentation testing was 5.9%, which demonstrates that the model is valid.

Patient and surgical variability in the primary stability of cementless acetabular cups: A finite element study

Aseptic loosening is the most common indication for revision of cementless acetabular cups and often depends on the primary stability achieved following surgery. Cup designs must be capable of achieving primary stability for a wide variety of individuals and surgical conditions to be successful. Typically, preclinical finite element (FE) testing of cups involves assessing the performance in a single patient and under a limited set of idealized conditions. The aim of this study was to assess the effect of patient and surgical parameters on the primary stability of an acetabular cup design in a set of subject-specific FE models. Interference fit was varied in a representative set of 12 patient-specific models of the implanted hemipelvis. Linear mixed models showed a significant association with micromotion for interference fit (P < .0001), acetabular bone elastic modulus (P < .001), native acetabular diameter (P  = .03), and the interference fit-elastic modulus interaction (P = .01). There were no significant associations between the polar gap and any of the parameters considered. The significant interference fit-elastic modulus interaction suggests that increasing the interference fit in patients with low bone quality leads to a greater reduction in micromotion than in patients with higher bone quality. However, the significant association between percentage bone yielding and interference fit (P < .0001) suggests a higher periacetabular fracture risk at higher interference fits. This work supports the development of preclinical testing of cup designs for the broad range patients and surgical conditions a cup may face following surgery.

Effects of interfacial crack and implant material on mixed-mode stress intensity factor and prediction of interface failure of cemented acetabular cup

This study deals with the effects of interfacial crack and implant material on mixed-mode stress intensity factor and prediction of interface failure of the cemented acetabular cup. A three dimensional (3D) finite element (FE) model of implanted pelvic bone was developed based on the computed tomography (CT) scan data. Combinations of four materials were considered for implant material. To understand the influence of interfacial crack at bone-cement and cement-implant interfaces on failure, 2D cracked models were developed based on the FE model and solved using the element-free Galerkin method (EFGM) by considering a rectangular section in the superior, inferior, anterior, and posterior locations. Interface failure was predicted in terms of mixed-mode stress intensity factor (SIF). The stress values obtained from FE analysis were transferred at the cut boundary of the rectangular section and considered as a mixed-mode loading condition to determine the SIF in the superior, inferior, anterior, and posterior locations at bone-cement and cement-implant interfaces using EFGM. Location wise, anterior seems to have more chances of failure because SIF in the anterior location was found to be higher than other locations. The bone-cement interface has more SIF and indicated more chances of failure than the cement-implant interface. Less SIF was found for the ceramic-ceramic material combination than other material combinations.

Pelvic Construct Prediction of Trabecular and Cortical Bone Structural Architecture

The pelvic construct is an important part of the body as it facilitates the transfer of upper body weight to the lower limbs and protects a number of organs and vessels in the lower abdomen. In addition, the importance of the pelvis is highlighted by the high mortality rates associated with pelvic trauma. This study presents a mesoscale structural model of the pelvic construct and the joints and ligaments associated with it. Shell elements were used to model cortical bone, while truss elements were used to model trabecular bone and the ligaments and joints. The finite element (FE) model was subjected to an iterative optimization process based on a strain-driven bone adaptation algorithm. The bone model was adapted to a number of common daily living activities (walking, stair ascent, stair descent, sit-to-stand, and stand-to-sit) by applying onto it joint and muscle loads derived using a musculoskeletal modeling framework. The cortical thickness distribution and the trabecular architecture of the adapted model were compared qualitatively with computed tomography (CT) scans and models developed in previous studies, showing good agreement. The sensitivity of the model to changes in material properties of the ligaments and joint cartilage and changes in parameters related to the adaptation algorithm was assessed. Changes to the target strain had the largest effect on predicted total bone volumes. The model showed low sensitivity to changes in all other parameters. The minimum and maximum principal strains predicted by the structural model compared to a continuum CT-derived model in response to a common test loading scenario showed good agreement with correlation coefficients of 0.813 and 0.809, respectively. The developed structural model enables a number of applications such as fracture modeling, design, and additive manufacturing of frangible surrogates.

Finite Element Analysis to Probe the Influence of Acetabular Shell Design, Liner Material, and Subject Parameters on Biomechanical Response in Periprosthetic Bone

The implant stability and biomechanical response of periprosthetic bone in acetabulum around total hip joint replacement (THR) devices depend on a host of parameters, including design of articulating materials, gait cycle and subject parameters. In this study, the impact of shell design (conventional, finned, spiked, and combined design) and liner material on the biomechanical response of periprosthetic bone has been analyzed using finite element (FE) method. Two different liner materials: high density polyethylene-20% hydroxyapatite-20% alumina (HDPE-20%HA-20%Al2O3) and highly cross-linked ultrahigh molecular weight polyethylene (HC-UHMWPE) were used. The subject parameters included bone condition and bodyweight. Physiologically relevant load cases of a gait cycle were considered. The deviation of mechanical condition of the periprosthetic bone due to implantation was least for the finned shell design. No significant deviation was observed at the bone region adjacent to the spikes and the fins. This study recommends the use of the finned design, particularly for weaker bone conditions. For stronger bones, the combined design may also be recommended for higher stability. The use of HC-UHMWPE liner was found to be better for convensional shell design. However, similar biomechanical response was captured in our FE analysis for both the liner materials in case of other shell designs. Overall, the study establishes the biomechanical response of periprosthetic bone in the acetabular with preclinically tested liner materials together with new shell design for different subject conditions.

Bone remodelling around the tibia due to total ankle replacement: effects of implant material and implant-bone interfacial conditions

One of the major causes of implant loosening is due to excessive bone resorption surrounding the implant due to bone remodelling. The objective of the study is to investigate the effects of implant material and implant-bone interface conditions on bone remodelling around tibia bone due to total ankle replacement. Finite element models of intact and implanted ankles were developed using CT scan data sets. Bone remodelling algorithm was used in combination with FE analysis to predict the bone density changes around the ankle joint. Dorsiflexion, neutral, and plantar flexion positions were considered, along with muscle force and ligaments. Implant-bone interfacial conditions were assumed as debonded and bonded to represent non-osseointegration and fully osseointegration at the porous coated surface of the implant. To investigate the effect of implant material, three finite element models having different material combinations of the implant were developed. For model 1, tibial and talar components were made of Co-Cr-Mo, and meniscal bearing was made of UHMWPE. For model 2, tibial and talar components were made of ceramic and meniscal bearing was made of UHMWPE. For model 3, tibial and talar components were made of ceramic and meniscal bearing was made of CFR-PEEK. Changes in implant material showed no significant changes in bone density due to bone remodelling. Therefore, ceramic appears to be a viable alternative to metal and CFR-PEEK can be used in place of UHMWPE. This study also indicates that proper bonding between implant and bone is essential for long-term survival of the prosthetic components.

Effect of Femoral Head Size, Subject Weight, and Activity Level on Acetabular Cement Mantle Stress Following Total Hip Arthroplasty

In cases where cemented components are used in total hip arthroplasty, damage, or disruption of the cement mantle can lead to aseptic loosening and joint failure. Currently, the relationship between subject activity level, obesity, and prosthetic femoral head size and the risk of aseptic loosening of the acetabular component in cemented total hip arthroplasty is not well understood. This study aims to provide an insight into this. Finite element models, validated with experimental data, were developed to investigate stresses in the acetabular cement mantle and pelvic bone resulting from the use of three prosthetic femoral head sizes, during a variety of daily activities and one high impact activity (stumbling) for a range of subject body weights. We found that stresses in the superior quadrants of the cortical bone-cement interface increased with prosthetic head size, patient weight, and activity level. In stumbling, average von Mises stresses (22.4 MPa) exceeded the bone cement yield strength for an obese subject (143 kg) indicating that the cement mantle would fail. Our results support the view that obesity and activity level are potential risk factors for aseptic loosening of the acetabular component and provide insight into the increased risk of joint failure associated with larger prosthetic femoral heads. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:1771-1783, 2019.

Influence of stems and metaphyseal sleeve on primary stability of cementless revision tibial trays used to reconstruct AORI IIB defects

Metaphyseal augments, such as sleeves, have been introduced to augment the fixation of revision total knee replacement (rTKR) components, and can be used with or without a stem. The effect of sleeve size in combination with stems on the primary stability and load transfer of a rTKR implant in AORI type IIB defects where the defect involves both condyles are poorly understood. The aim of this study was to examine the primary stability of revision tibial tray augmented with a sleeve in an AORI type IIB defect which involves both condyles with loss of cortical and cancellous bone. Finite element models were generated from computed tomography (CT) scans of nine individuals. All the bones used in the study had an AORI type IIB defect. The cohort included eight females (mean weight: 64 kg, height: 1.6 m). Material properties were sampled from CT data and assigned to the FE model. Joint contact forces for level gait, stair descent, and squat were applied. Stemless sleeved implants under various loading conditions were shown to have adequate primary stability in all AORI type IIB defects investigated. Adding a stem only marginally improved the primary stability of the implant but reduced the strain in the metaphysis compared to stemless implants. Once good initial mechanical stability was established with a sleeve, there was no benefit, in terms of primary stability or bone strains, from increasing sleeve size. This study suggests that metaphyseal sleeves, without a stem, can provide the required primary stability required by a rTKR tibial implant, to reconstruct an AORI type IIB defect. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res.

In Silico Pelvis and Sacroiliac Joint Motion: Refining a Model of the Human Osteoligamentous Pelvis for Assessing Physiological Load Deformation Using an Inverted Validation Approach

Introduction. Computational modeling of the human pelvis using the finite elements (FE) method has become increasingly important to understand the mechanisms of load distribution under both healthy and pathologically altered conditions and to develop and assess novel treatment strategies. The number of accurate and validated FE models is however small, and given models fail resembling the physiologic joint motion in particular of the sacroiliac joint. This study is aimed at using an inverted validation approach, using in vitro load deformation data to refine an existing FE model under the same mode of load application and to parametrically assess the influence of altered morphology and mechanical data on the kinematics of the model. Materials and Methods. An osteoligamentous FE model of the pelvis including the fifth lumbar vertebra was used, with highly accurate representations of ligament orientations. Material properties were altered parametrically for bone, cartilage, and ligaments, followed by changes in bone geometry (solid versus 3 and 2 mm shell) and material models (linear elastic, viscoelastic, and hyperelastic isotropic), and the effects of varying ligament fiber orientations were assessed. Results. Elastic modulus changes were more decisive in both linear elastic and viscoelastic bone, cartilage, and ligaments models, especially if shell geometries were used for the pelvic bones. Viscoelastic material properties gave more realistic results. Surprisingly little change was observed as a consequence of altering SIJ ligament orientations. Validation with in vitro experiments using cadavers showed close correlations for movements especially for 3 mm shell viscoelastic model. Discussion. This study has used an inverted validation approach to refine an existing FE model, to give realistic and accurate load deformation data of the osteoligamentous pelvis and showed which variation in the outcomes of the models are attributed to altered material properties and models. The given approach furthermore shows the value of accurate validation and of using the validation data to fine tune FE models.

Effects of implant orientation and implant material on tibia bone strain, implant–bone micromotion, contact pressure, and wear depth due to total ankle replacement

The aim of this study is to investigate the effects of implant orientation and implant material on tibia bone strain, implant–bone micromotion, maximum contact pressure, and wear depth at the articulating surface due to total ankle replacement. Three-dimensional finite element models of intact and implanted ankle were developed from computed tomography scan data. Four implanted models were developed having varus and valgus orientations of 5° and 10°, respectively. In order to determine the effect of implant material combination on tibia bone strain, micromotion, contact pressure, and wear depth, three other finite element models were developed having a different material combination of the implant. Dorsiflexion, neutral, and plantarflexion positions were considered as applied loading condition, along with muscle force and ligaments. Implant orientation alters the strain distribution in tibia bone. Strain shielding was found to be less in the case of the optimally positioned implant. Apart from the strain, implant orientation also affects implant–bone micromotion, contact pressure, and wear depth. Implant materials have less influence on tibia bone strain and micromotion. However, wear depth was reduced when ceramic and carbon fibre–reinforced polyetheretherketone material combination was used. Proper orientation of the implant is important to reduce the strain shielding. The present result suggested that ceramic can be used as an alternative to metal and carbon fibre–reinforced polyetheretherketone as an alternative to ultra-high molecular weight polyethylene to reduce wear, which would be beneficial for long-term success and fixation of the implant.

Patient-Specific Static Structural Analysis of Femur Bone of different lengths

Background:

The femur bone is an essential part of human activity, providing stability and support in carrying out our day to day activities. The inter-human anatomical variation and load bearing ability of humans of different heights will provide the necessary understanding of their functional ability.

Objective:

In this study, femur bone of two humans of different lengths (tall femur and short femur) were subjected to static structural loading conditions to evaluate their load-bearing abilities using Finite Element Analysis.

Methods:

The 3D models of femur bones were developed using MIMICS from the CT scans which were then subjected to static structural analysis by varying the load from 1000N to 8000N. The von Mises stress and deformation were captured to compare the performance of each of the femur bones.

Results:

The tall femur resulted in reduced Von-Mises stress and total deformation when compared to the short femur. However, the maximum principle stresses showed an increase with an increase in the bone length. In both the femurs, the maximum stresses were observed in the medullary region.

Conclusion:

When the applied load exceeds 10 times the body weight of the person, the tall femur model exceeded 134 MPa stress value. The short femur model failed at 9 times the body weight, indicating that the tall femur had higher load-bearing abilities.

Cement interface and bone stress in total hip arthroplasty: Relationship to head size

The use of larger prosthetic femoral heads in total hip arthroplasty (THA) has increased considerably in recent years in response to the need to improve joint stability and reduce risk of dislocation. However, data suggests larger femoral heads are associated with higher joint failure rates. For cemented implants, ensuring the continued integrity of the cement mantle is key to long term fixation. This paper describes an investigation into the effect of variation in femoral head size on stresses in the acetabular cement mantle and pelvic bone. Three commonly used femoral head sizes: 28, 32, and 36 mm diameter were investigated. The study was undertaken using a finite element model validated using surface strains obtained from Digital Image Correlation (DIC) during experimentation on a composite hemipelvis implanted with a cemented all-polyethylene acetabular cup. Following validation, the models were used to investigate stresses in the pelvic bone and acetabular cement mantle resulting from two loading scenarios; an average weight subject (700 N) and an overweight subject (1,000 N) undertaking a single leg stand. We found that the highest peak stresses occurred in the anterosuperior and posterosuperior regions of the bone-cement interface, in the line of action of the load, where debonding usually initiates. Stress on the cortical bone-cement interface increased with femoral head diameter by up to 9% whilst stresses in the trabecular bone remained relatively invariant. Our findings may help to explain higher joint failure rates associated with larger femoral heads. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:2966-2977, 2018.

Effects of Cutting the Sacrospinous and Sacrotuberous Ligaments

The sacrospinous (SS) and sacrotuberous (ST) ligaments form a complex at the posterior pelvis, with an assumed role as functional stabilizers. Experimental and clinical research has yielded controversial results regarding their function, both proving and disproving their role as pelvic stabilizers. These findings have implications for strategies for treating pelvic injury and pain syndromes. The aim of the present simulation study was to assess the influence of altered ligament function on pelvis motion. A finite elements computer model was used. The two-leg stance was simulated, with the load of body weight applied via the fifth lumbar vertebra and both femora, allowing for nutation of the sacroiliac joint. The in-silico kinematics were validated with in-vitro experiments using the same scenario of load application following SS and ST transection in six human cadavers. Modeling of partial or complete ligament failure caused significant increases in pelvis motion. This effect was most pronounced if the SS and ST were affected with 164% and 182%, followed by the sacroiliac and iliolumbar ligaments with 123% and 147%, and the pubic ligaments with 113% and 119%, for partial and complete disruption, respectively. Simultaneous ligament transection multiplied the effects on pelvis motion by up to 490%. Unilateral ligament injury altered the motion at the pelvis contralaterally. The experiments presented here provide strong evidence for the stabilizing role of the SS and ST. A fortiori, the instability resulting from partial or complete SS and ST injury merits consideration in treatment strategies involving these ligaments as important stabilizers. Clin. Anat. 32:231-237, 2019. © 2018 Wiley Periodicals, Inc.

Biomechanical analysis of supra-acetabular insufficiency fracture using finite element analysis

BACKGROUND

Supra-acetabular insufficiency fractures (SAIFs) occur in the upper acetabulum and are rare compared with insufficiency sacral, femoral head, or ischial fractures. However, SAIFs are known to occur in low grade trauma, and the underlying mechanism is still remained unclear.

METHODS

We performed biomechanical analysis using finite element analysis to clarify the mechanisms underlying the development of SAIFs. Patient-specific models and bone mineral density (BMD) were derived from pelvic computed tomography data from two patients with SAIF (unaffected side) and two healthy young adults. The bone was assumed to be an isotropic, linearly elastic body. We assigned Young's modulus of each element to the pelvis based on the BMD, and reported the relationships for BMD-modulus. Clinically relevant loading conditions-walking and climbing stairs-were applied to the models. We compared the region of failure risk in each acetabulum using a maximum principal strain criterion.

RESULTS

The average supra-acetabular BMD was less than that of the hemi-pelvis and femoral head, but was higher than that of the femoral neck and greater trochanter. Greater minimum principal strain was concentrated in the supra-acetabular portion in both the SAIF and healthy models. In the SAIF models, the higher region of the failure risk matched the fracture site on the acetabulum.

CONCLUSIONS

Relative fragility causes compressive strain to concentrate in the upper acetabulum when walking and climbing stairs. When presented with a patient complaining of hip pain without apparent trauma or abnormal X-ray findings, physicians should consider the possibility of SAIF and perform magnetic resonance imaging for the diagnosis of SAIF.

The effect of cement mantle thickness on strain energy density distribution and prediction of bone density changes around cemented acetabular component

Cement mantle thickness is known to be one of the important parameters to reduce the failure of the cemented acetabular component. The thickness of the cement mantle is also often influenced by the positioning of the acetabular cup. The aim of this study is to determine the effect of uniform and non-uniform cement mantle thickness on strain energy density distribution and prediction of the possibility of bone remodelling around the acetabular region. Furthermore, tensile stress distribution in the cement mantle due to non-uniform cement mantle thickness was also investigated. Three-dimensional finite element models of intact and 17 implanted pelvic bone were developed based on computed tomography data sets. Results indicate that implantation with non-uniform cement thickness variation in the anterior-posterior direction has a significant influence on strain energy density distribution around the acetabulum as compared to thickness variation in the superior-inferior direction. Increase in density is predicted at the anterior part of the acetabulum, whereas density decrease is predicted at the posterior, inferior and superior part of the acetabulum. The non-uniform cement mantle thickness affected the tensile stress distribution in the cement mantle, in particularly superiorly placed acetabular cup. This study concludes that uniform cement thickness is desired for the longer success of the cemented acetabular component.

Influence of varying stem and metaphyseal sleeve size on the primary stability of cementless revision tibial trays used to reconstruct AORI IIA defects. A simulation study

Traditionally, diaphyseal stems have been utilized to augment the stability of revision total knee replacement (rTKR) implants. More recently metaphyseal augments, such as sleeves, have been introduced to further augment component fixation. The effect of augments such as stems and sleeves have on the primary stability of a rTKR implant is poorly understood, however it has important implications on the complexity, costs and survivorship of the procedure. Finite element analysis was used to investigate the primary stability and strain distribution of various size stems and sleeves used in conjunction with a cementless revision tibial tray. The model was built from computer tomography images of a single healthy tibia obtained from an 81-year-old patient to which an Anderson Orthopaedic Research Institute (AORI) IIA defect was virtually added. The influences of varying body mass index (BMI) and bone modulus were also investigated. Stemless sleeves were found to provided adequate primary implant stability (average implant micro-motion <50 μm) for the studied defect. Addition of a stem did not enhance the primary stability. Furthermore, this study found that varying BMI and bone modulus had a considerable effect on strain distribution but negligible effect on micro-motion in the sleeve area. In conclusion, the addition of diaphyseal stem to a metaphyseal sleeve had little benefit in enhancing the primary stability of tibial trays augmented when simulating reconstructions of AORI IIA tibial defects. Additional studies are required to determine the relative benefit of the diaphyseal stem when using metaphyseal sleeves defects with more extensive bone loss. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1876-1886, 2018.

Combined Bone Ingrowth and Remodeling Around Uncemented Acetabular Component: A Multiscale Mechanobiology-Based Finite Element Analysis

Bone ingrowth and remodeling are two different evolutionary processes which might occur simultaneously. Both these processes are influenced by local mechanical stimulus. However, a combined study on bone ingrowth and remodeling has rarely been performed. This study is aimed at understanding the relationship between bone ingrowth and adaptation and their combined influence on fixation of the acetabular component. Based on three-dimensional (3D) macroscale finite element (FE) model of implanted pelvis and microscale FE model of implant–bone interface, a multiscale framework has been developed. The numerical prediction of peri-acetabular bone adaptation was based on a strain-energy density-based formulation. Bone ingrowth in the microscale models was simulated using the mechanoregulatory algorithm. An increase in bone strains near the acetabular rim was observed in the implanted pelvis model, whereas the central part of the acetabulum was observed to be stress shielded. Consequently, progressive bone apposition near the acetabular rim and resorption near the central region were observed. Bone remodeling caused a gradual increase in the implant–bone relative displacements. Evolutionary bone ingrowth was observed around the entire acetabular component. Poor bone ingrowth of 3–5% was predicted around the centro-inferio and inferio-posterio-superio-peripheral regions owing to higher implant–bone relative displacements, whereas the anterio-inferior and centro-superior regions exhibited improved bone ingrowth of 35–55% due to moderate implant–bone relative displacement. For an uncemented acetabular CoCrMo component, bone ingrowth had hardly any effect on bone remodeling; however, bone remodeling had considerable influence on bone ingrowth.

EVALUATION AND PREDICTION OF HUMAN GAIT PARAMETERS USING UNIVARIATE, MULTIVARIATE AND STEPWISE STATISTICAL METHODS

This research was carried out to establish the relationship between human anthropometric data and corresponding gait variables. A group comprising 35 participants (18 male and 17 female) was selected for the current study. The study consisted of trials in which each participant was asked to walk the length of the instrumented walkway (Kistler’s force platform inset) at a self-selected speed. Using a four-camera motion analysis system, the kinematic and kinetic parameters of each trial were calculated. The peak values obtained from the data curves were used to generate the necessary regression fits. In order to establish the correlation between the anthropometric data of human and the gait parameters, the univariate, multivariate and stepwise fits were generated. Further, the statistical methods were employed to evaluate the [Formula: see text], [Formula: see text] and [Formula: see text]-values for each fit. The current multivariate study indicates an increasing trend in [Formula: see text] values and decreasing trend for [Formula: see text]-values when compared with the univariate fits and the results follow the expected line.

A numerical study on stress distribution across the ankle joint: Effects of material distribution of bone, muscle force and ligaments

OBJECTIVE

The goal of this study is to develop a realistic three dimensional FE model of intact ankle joint.

METHODS

Three dimensional FE model of the intact ankle joint was developed using computed tomography data sets. The effect of muscle force, ligaments and proper material property distribution of bone on stress distribution across the intact ankle joint was studied separately.

RESULTS

Present study indicates bone material property, ligaments and muscle force have influence on stress distribution across the ankle joint.

CONCLUSION

Proper bone material, ligaments and muscle must be considered in the computational model for pre-clinical analysis of ankle prosthesis.

Influence of Implant Surface Texture Design on Peri-Acetabular Bone Ingrowth: A Mechanobiology Based Finite Element Analysis

The fixation of uncemented acetabular components largely depends on the amount of bone ingrowth, which is influenced by the design of the implant surface texture. The objective of this numerical study is to evaluate the effect of these implant texture design factors on bone ingrowth around an acetabular component. The novelty of this study lies in comparative finite element (FE) analysis of 3D microscale models of the implant-bone interface, considering patient-specific mechanical environment, host bone material property and implant-bone relative displacement, in combination with sequential mechanoregulatory algorithm and design of experiment (DOE) based statistical framework. Results indicated that the bone ingrowth process was inhibited due to an increase in interbead spacing from 200 μm to 600 μm and bead diameter from 1000 μm to 1500 μm and a reduction in bead height from 900 μm to 600 μm. Bead height, a main effect, was found to have a predominant influence on bone ingrowth. Among the interaction effects, the combination of bead height and bead diameter was found to have a pronounced influence on bone ingrowth process. A combination of low interbead spacing (P = 200 μm), low bead diameter (D = 1000 μm), and high bead height (H = 900 μm) facilitated peri-acetabular bone ingrowth and an increase in average Young's modulus of newly formed tissue layer. Hence, such a surface texture design seemed to provide improved fixation of the acetabular component.

A Computational Efficient Method to Assess the Sensitivity of Finite-Element Models: An Illustration With the Hemipelvis

Assessing the sensitivity of a finite-element (FE) model to uncertainties in geometric parameters and material properties is a fundamental step in understanding the reliability of model predictions. However, the computational cost of individual simulations and the large number of required models limits comprehensive quantification of model sensitivity. To quickly assess the sensitivity of an FE model, we built linear and Kriging surrogate models of an FE model of the intact hemipelvis. The percentage of the total sum of squares (%TSS) was used to determine the most influential input parameters and their possible interactions on the median, 95th percentile and maximum equivalent strains. We assessed the surrogate models by comparing their predictions to those of a full factorial design of FE simulations. The Kriging surrogate model accurately predicted all output metrics based on a training set of 30 analyses (R2 = 0.99). There was good agreement between the Kriging surrogate model and the full factorial design in determining the most influential input parameters and interactions. For the median, 95th percentile and maximum equivalent strain, the bone geometry (60%, 52%, and 76%, respectively) was the most influential input parameter. The interactions between bone geometry and cancellous bone modulus (13%) and bone geometry and cortical bone thickness (7%) were also influential terms on the output metrics. This study demonstrates a method with a low time and computational cost to quantify the sensitivity of an FE model. It can be applied to FE models in computational orthopaedic biomechanics in order to understand the reliability of predictions.

Does the optimal position of the acetabular fragment should be within the radiological normal range for all developmental dysplasia of the hip? A patient-specific finite element analysis

Background

The success of Bernese periacetabular osteotomy depends significantly on how extent the acetabular fragment can be corrected to its optimal position. This study was undertaken to investigate whether correcting the acetabular fragment into the so-called radiological “normal” range is the best choice for all developmental dysplasia of the hip with different severities of dysplasia from the biomechanical view? If not, is there any correlation between the biomechanically optimal position of the acetabular fragment and the severity of dysplasia?

Methods

Four finite element models with different severities of dysplasia were developed. The virtual periacetabular osteotomy was performed with the acetabular fragment rotated anterolaterally to incremental center-edge angles; then, the contact area and pressure and von Mises stress in the cartilage were calculated at different correction angles.

Results

The optimal position of the acetabular fragment for patients 1, 2, and 3 was when the acetabular fragment rotated 17° laterally (with the lateral center-edge angle of 36° and anterior center-edge angle of 58°; both were slightly larger than the “normal” range), 25° laterally following further 5° anterior rotation (with the lateral center-edge angle of 31° and anterior center-edge angle of 51°; both were within the “normal” range), and 30° laterally following further 10° anterior rotation (with the lateral center-edge angle of 25° and anterior center-edge angle of 40°; both were less than the “normal” range), respectively.

Conclusions

The optimal corrective position of the acetabular fragment is severity dependent rather than within the radiological “normal” range for developmental dysplasia of the hip. We prudently proposed that the optimal correction center-edge angle of mild, moderate, and severe developmental dysplasia of the hip is slightly larger than the “normal” range, within the “normal” range, and less than the lower limit of the “normal” range, respectively.

Electronic supplementary material

The online version of this article (doi:10.1186/s13018-016-0445-3) contains supplementary material, which is available to authorized users.

The effects of musculoskeletal loading regimes on numerical evaluations of acetabular component

The importance of clinical studies notwithstanding, the failure assessment of implant–bone structure has alternatively been carried out using finite element analysis. However, the accuracy of the finite element predicted results is dependent on the applied loading and boundary conditions. Nevertheless, most finite element–based evaluations on acetabular component used a few selective load cases instead of the eight load cases representing the entire gait cycle. These in silico evaluations often suffer from limitations regarding the use of simplified musculoskeletal loading regimes. This study attempts to analyse the influence of three different loading regimes representing a gait cycle, on numerical evaluations of acetabular component. Patient-specific computer tomography scan-based models of intact and resurfaced pelvises were used. One such loading regime consisted of the second load case that corresponded to peak hip joint reaction force. Whereas the other loading regime consisted of the second and fifth load cases, which corresponded to peak hip joint reaction force and peak muscle forces, respectively. The third loading regime included all the eight load cases. Considerable deviations in peri-acetabular strains, standard error ranging between 115 and 400 µε, were observed for different loading regimes. The predicted bone strains were lower when selective loading regimes were used. Despite minor quantitative variations in bone density changes (less than 0.15 g cm−3), the final bone density pattern after bone remodelling was found to be similar for all the loading regimes. Underestimations in implant–bone micromotions (40–50 µm) were observed for selective loading regimes after bone remodelling. However, at immediate post-operative condition, such underestimations were found to be less (less than 5 µm). The predicted results highlight the importance of inclusion of eight load cases representing the gait cycle for in silico evaluations of resurfaced pelvis.

Assessment of failure of cemented polyethylene acetabular component due to bone remodeling: A finite element study

The aim of the study is to determine failure of the cemented polyethylene acetabular component, which might occur due to excessive bone resorption, cement-bone interface debonding and fatigue failure of the cement mantle. Three-dimensional finite element models of intact and implanted pelvic bone were developed and bone remodeling algorithm was implemented for present analysis. Soderberg fatigue failure diagram was used for fatigue assessment of the cement mantle. Hoffman failure criterion was considered for prediction of cement-bone interface debonding. Results indicate fatigue failure of the cement mantle and implant-bone interface debonding might not occur due to bone remodeling.

Bone ingrowth around porous-coated acetabular implant: a three-dimensional finite element study using mechanoregulatory algorithm

Fixation of uncemented implant is influenced by peri-prosthetic bone ingrowth, which is dependent on the mechanical environment of the implant-bone structure. The objective of the study is to gain an insight into the tissue differentiation around an acetabular component. A mapping framework has been developed to simulate appropriate mechanical environment in the three-dimensional microscale model, implement the mechanoregulatory tissue differentiation algorithm and subsequently assess spatial distribution of bone ingrowth around an acetabular component, quantitatively. The FE model of implanted pelvis subjected to eight static load cases during a normal walking cycle was first solved. Thereafter, a mapping algorithm has been employed to include the variations in implant-bone relative displacement and host bone material properties from the macroscale FE model of implanted pelvis to the microscale FE model of the beaded implant-bone interface. The evolutionary tissue differentiation was observed in each of the 13 microscale models corresponding to 13 acetabular regions. The total implant-bone relative displacements, averaged over each region of the acetabulum, were found to vary between 10 and 60 μm. Both the linear elastic and biphasic poroelastic models predicted similar mechanoregulatory peri-prosthetic tissue differentiation. Considerable variations in bone ingrowth (13-88%), interdigitation depth (0.2-0.82 mm) and average tissue Young's modulus (970-3430 MPa) were predicted around the acetabular cup. A progressive increase in the average Young's modulus, interdigitation depth and decrease in average radial strains of newly formed tissue layer were also observed. This scheme can be extended to investigate tissue differentiation for different surface texture designs on the implants.