Advanced material modelling in numerical simulation of primary acetabular press-fit cup stability

Primary stability of artificial acetabular cups, used for total hip arthroplasty, is required for the subsequent osteointegration and good long-term clinical results of the implant. Although closed-cell polymer foams represent an adequate bone substitute in experimental studies investigating primary stability, correct numerical modelling of this material depends on the parameter selection.

Material parameters necessary for crushable foam plasticity behaviour were originated from numerical simulations matched with experimental tests of the polymethacrylimide raw material. Experimental primary stability tests of acetabular press-fit cups consisting of static shell assembly with consecutively pull-out and lever-out testing were subsequently simulated using finite element analysis.

Identified and optimised parameters allowed the accurate numerical reproduction of the raw material tests. Correlation between experimental tests and the numerical simulation of primary implant stability depended on the value of interference fit. However, the validated material model provides the opportunity for subsequent parametric numerical studies.     

Influence of anisotropic bone properties on the biomechanical behavior of the acetabular cup implant: a multiscale finite element study

Although the biomechanical behavior of the acetabular cup (AC) implant is determinant for the surgical success, it remains difficult to be assessed due to the multiscale and anisotropic nature of bone tissue. The aim of the present study was to investigate the influence of the anisotropic properties of peri-implant trabecular bone tissue on the biomechanical behavior of the AC implant at the macroscopic scale. Thirteen bovine trabecular bone samples were imaged using micro-computed tomography (μCT) with a resolution of 18 μm. The anisotropic biomechanical properties of each sample were determined at the scale of the centimeter based on a dedicated method using asymptotic homogenization. The material properties obtained with this multiscale approach were used as input data in a 3D finite element model to simulate the macroscopic mechanical behavior of the AC implant under different loading conditions. The largest stress and strain magnitudes were found around the equatorial rim and in the polar area of the AC implant. All macroscopic stiffness quantities were significantly correlated (R² > 0.85, p < 6.5 e-6) with BV/TV (bone volume/total volume). Moreover, the maximum value of the von Mises stress field was significantly correlated with BV/TV (R² > 0.61, p < 1.6 e-3) and was always found at the bone-implant interface. However, the mean value of the microscopic stress (at the scale of the trabeculae) decrease as a function of BV/TV for vertical and torsional loading and do not depend on BV/TV for horizontal loading. These results highlight the importance of the anisotropic properties of bone tissue.     

The influence of uncemented femoral stem length and design on its primary stability: a finite element analysis

One of the crucial factors for short- and long-term clinical success of total hip arthroplasty cementless implants is primary stability. Indeed, motion at the bone–implant interface above 40 μm leads to partial bone ingrowth, while motion exceeding 150 μm completely inhibits bone ingrowth. The aim of this study was to investigate the effect of two cementless femoral stem designs with different lengths on the primary stability. A finite element model of a composite Sawbones^® fourth generation, implanted with five lengths of the straight prosthesis design and four lengths of the curved prosthesis design, was loaded with hip joint and abductor forces representing two physiological activities: fast walking and stair climbing. We found that reducing the straight stem length from 146 to 54 mm increased the average micromotion from 17 to 52 μm during fast walking, while the peak value increased from 42 to 104 μm. With the curved stem, reducing length from 105 to 54 mm increased the average micromotion from 10 to 29 μm, while the peak value increased from 37 to 101 μm. Similar findings are obtained for stair climbing for both stems. Although the present study showed that femoral stem length as well as stem design directly influences its primary stability, for the two femoral stems tested, length could be reduced substantially without compromising the primary stability. With the aim of minimising surgical invasiveness, newer femoral stem design and currently well performing stems might be used with a reduced length without compromising primary stability and hence, long-term survivorship.     

Comparative analysis of bone remodelling models with respect to computerised tomography-based finite element models of bone

Subject-specific finite element models are an extensively used tool for the numerical analysis of the biomechanical behaviour of human bones. However, bone modelling is not an easy task due to the complex behaviour of bone tissue, involving non-homogeneous and anisotropic mechanical properties. Moreover, bone is a living tissue and therefore its microstructure and mechanical properties evolve with time in a known process called bone remodelling. This phenomenon has been widely studied, many being the numerical models that have been formulated to predict density distribution and its evolution in several bones. The aim of the present study is to assess the capability of a bone remodelling model to predict the bone density distribution of different types of human bone (femur, tibia and mandible) comparing the obtained results with the bone density estimated by means of computerised tomography. Good accuracy was observed for the bone remodelling predictions including the thickness of the cortical layer.     

Prediction of globe rupture caused by primary blast: a finite element analysis

Although a human eye comprises less than 0.1% of the frontal body surface area, injuries to the eye are found to be disproportionally common in survivors of explosions. This study aimed to introduce a Lagrangian–Eulerian coupling model to predict globe rupture resulting from primary blast effect. A finite element model of a human eye was created using Lagrangian mesh. An explosive and its surrounding air domain were modelled using Eulerian mesh. Coupling the two models allowed simulating the blast wave generation, propagation and interaction with the eye. The results showed that the peak overpressures caused by blast wave on the corneal apex are 2080, 932.1 and 487.3 kPa for the victim distances of 0.75, 1.0 and 1.25 m, respectively. Higher stress occurred at the limbus, where the peaks for the three victim distances are 25.5, 14.1 and 6.4 MPa. The overpressure threshold of globe rupture was determined as 2000 kPa in a small-scale explosion. The findings would provide insights into the mechanism of primary blast-induced ocular injuries.     

Mechanical significance of femoral head trabecular bone structure in Loris and Galago evaluated using micromechanical finite element models

Work on the interspecific and intraspecific variation of trabecular bone in the proximal femur of primates demonstrates important architectural variation between animals with different locomotor behaviors. This variation is thought to be related to the processes of bone adaptation whereby bone structure is optimized to the mechanical environment. Micromechanical finite element models were created for the proximal femur of the leaping Galago senegalensis and the climbing and quadrupedal Loris tardigradus by converting bone voxels from high-resolution X-ray computed tomography scans of the femoral head to eight-noded brick elements. The resulting models had approximately 1.8 million elements each. Loading conditions representing takeoff phase of a leap and more generalized load orientations were applied to the models, and the models were solved using the iterative "row-by-row" matrix-vector multiplication algorithm. The principal strain and Von Mises stress results for the leaping model were similar for both species at each load orientation. Similar hip joint reaction forces in the range of 4.9 x to 12 x body weight were calculated for both species under each loading condition, but the hip reaction values estimated for Loris were higher than predicted based on locomotor behavior. These results suggest that functional adaptation to hip joint loading may not fully explain the differences in femoral head trabecular bone structure in Galago and Loris. The finite element method represents a unique and useful tool for analyzing the functional adaptation of trabecular bone in a diversity of animals and for reconstructing locomotor behavior in extinct taxa.     

Using 3D finite element models verified the importance of callus material and microstructure in biomechanics restoration during bone defect repair

Background: There is lack of further observations on the microstructure and material property of callus during bone defect healing and the relationships between callus properties and the mechanical strength. Methods: Femur bone defect model was created in rabbits and harvested CT data to reconstruct finite element models at 1 and 2 months. Three types of assumed finite element models were compared to study the callus properties, which assumed the material elastic property as heterogeneous (R-model), homogenous (H-model) or did not change from 1 to 2 months (U-model). Results: The apparent elastic moduli increased at 2 months (from 355.58 ± 132.67 to 1139.30 ± 967.43 MPa) in R-models. But there was no significant difference in apparent elastic moduli between R-models (355.58 ± 132.67 and 1139.30 ± 967.43 MPa) and H-models (344.79 ± 138.73 and 1001.52 ± 692.12 MPa) in 1 and 2 months. A significant difference of apparent elastic moduli was found between the R-model (1139.30 ± 967.43 MPa) and U-model group (207.15 ± 64.60 MPa) in 2 months. Conclusions: This study showed that the callus structure stability remodeled overtime to achieve a more effective structure, while the material quality of callus tissue is a very important factor for callus strength. At the meantime, this study showed an evidence that the material heterogeneity maybe not as important as it is in bone fracture model.     

Factors affecting subject-specific finite element models of implant-fitted rat bone specimens: critical analysis of a technical protocol

The authors propose a protocol to derive finite element (FE) models from micro computer tomography scans of implanted rat bone. A semi-automatic procedure allows segmenting the images using specimen-specific bone mineral density (BMD) thresholds. An open-source FE model generator processes the segmented images to a quality tetrahedral mesh. The material properties assigned to each element are integrated from the BMD field. Piecewise, threshold-dependent density–elasticity relationships are implemented to limit the effects of metal artefacts. A detailed sensitivity study highlights the coherence of the generated models and quantifies the influence of the modelling parameters on the results. Two applications of the protocol are proposed. The stiffness of bare and implanted rat tibiae specimens is predicted by simulating three-point bending and inter-implant displacement, respectively. Results are compared with experimental tests. The mean value and the variability between the specimens are well captured in both tests.     

Evaluation of liver tissue damage and grasp stability using finite element analysis

Minimizing tissue damage and maintaining grasp stability are essential considerations in surgical grasper design. Most past and current research analyzing graspers used for tissue manipulation in minimally invasive surgery is based on in vitro experiments. Most previous work assessed tissue injury and grasp security by visual inspection; only a few studies have quantified it. The goal of the present work is to develop a methodology with which to compute tissue damage magnitude and grasp quality that is appropriate for a wide range of grasper–tissue interaction. Using finite element analysis (FEA), four graspers with varying radii of curvature and four graspers with different tooth sizes were analyzed while squeezing and pulling liver tissue. All graspers were treated as surgical steel with linear elastic material properties. Nonlinear material properties of tissue used in the FEA as well as damage evaluation were derived from previously reported in vivo experiments. Computed peak stress, integrated stress, and tissue damage were compared. Applied displacement is vertical and then horizontal to the tissue surface to represent grasp and retraction. A close examination of the contact status of each node within the grasper–tissue interaction surface was carried out to investigate grasp stability. The results indicate less tissue damage with increasing radius of curvature. A smooth wave pattern reduced tissue damage at the cost of inducing higher percentage of slipping area. This methodology may be useful for researchers to develop and test various designs of graspers. Also it could improve surgical simulator performance by reflecting more realistic tissue material properties and predicting tissue damage for the student.     

Accuracy of specimen-specific nonlinear finite element analysis for evaluation of radial diaphysis strength in cadaver material

The feasibility of a user-specific finite element model for predicting the in situ strength of the radius after implantation of bone plates for open fracture reduction was established. The effect of metal artifact in CT imaging was characterized. The results were verified against biomechanical test data. Fourteen cadaveric radii were divided into two groups: (1) intact radii for evaluating the accuracy of radial diaphysis strength predictions with finite element analysis and (2) radii with a locking plate affixed for evaluating metal artifact. All bones were imaged with CT. In the plated group, radii were first imaged with the plates affixed (for simulating digital plate removal). They were then subsequently imaged with the locking plates and screws removed (actual plate removal). Fracture strength of the radius diaphysis under axial compression was predicted with a three-dimensional, specimen-specific, nonlinear finite element analysis for both the intact and plated bones (bones with and without the plate captured in the scan). Specimens were then loaded to failure using a universal testing machine to verify the actual fracture load. In the intact group, the physical and predicted fracture loads were strongly correlated. For radii with plates affixed, the physical and predicted (simulated plate removal and actual plate removal) fracture loads were strongly correlated. This study demonstrates that our specimen-specific finite element analysis can accurately predict the strength of the radial diaphysis. The metal artifact from CT imaging was shown to produce an overestimate of strength.     

Validation of an efficient method of assigning material properties in finite element analysis of pelvic bone

Bone in the pelvis is a composite material with a complex anatomical structure that is difficult to model computationally. Rather than assigning material properties to increasingly smaller elements to capture detail in three-dimensional finite element (FE) models, properties can be assigned to Gauss points within larger elements. As part of a validation process, we compared experimental and analytical results from a composite beam under four-point load to FE models with material properties assigned to refined elements and Gauss points within larger elements. Both FE models accurately predicted deformation and the analytical predictions of internal shear stress.     

A finite element analysis of the vibration behaviour of a cementless hip system

An early diagnosis of aseptic loosening of a total hip replacement (THR) by plain radiography, scintigraphy or arthography has been shown to be less reliable than using a vibration technique. However, it has been suggested that it may be possible to distinguish between a secure and a loose prosthesis using a vibration technique. In fact, vibration analysis methods have been successfully used to assess dental implant stability, to monitor fracture healing and to measure bone mechanical properties. Several studies have combined the vibration technique with the finite element (FE) method in order to better understand the events involved in the experimental technique. In the present study, the main goal is to simulate the change in the resonance frequency during the osseointegration process of a cementless THR (Zweymüller). The FE method was used and a numerical modal analysis was conducted to obtain the natural frequencies and mode shapes under vibration. The effects were studied of different bone and stem material properties, and different contact conditions at the bone–implant interface. The results were in agreement with previous experimental and computational observations, and differences among the different cases studied were detected. As the osseointegration process at the bone–implant interface evolved, the resonance frequency values of the femur–prosthesis system also increased. In summary, vibration analysis combined with the FE method was able to detect different boundary conditions at the bone–implant interface in cases of both osseointegration and loosening.     

Three-dimensional finite element analysis of the composite and compomer onlays in primary molars

Resin onlay restoration is an esthetic alternative technique used for restoring extensively damaged primary molars. Understanding the behavior of materials under repeated functional stress and how the stress is transmitted to the remaining tooth structure is important. The aim of this study was to compare stresses in primary molars restored with indirect composite and compomer onlay. 3D frame models of the right mandibular and maxillary primary molars and the alveolar bone were created using computerized tomography images of a six-year-old girl. The enamel and dentine layers above the cement layer were unified to generate onlay restoration, and composite and compomer were used as restorative materials. The vertical occlusal load (100?N) was applied to the teeth in the occlusal contact areas. The von Mises stress distributions and normal stress distributions of the y-axis (parallel to the long axis of tooth) were evaluated. The occlusal stress is transmitted to the cervical part of healthy teeth by spreading it through the enamel layer. The composite and compomer restorative materials exhibited similar stress distribution patterns. An indirect technique creates a structure similar to the original morphological form, and it allows restorations to distribute high occlusal stresses and to minimize possible breakages.     

Comparison between mechanical stress and bone mineral density in the femur after total hip arthroplasty by using subject-specific finite element analyses

The mechanism underling bone mineral density (BMD) loss that occurs in the femur after total hip arthroplasty (THA) remains unknown. We compared the equivalent stress and strain energy density (SED) to BMD in the femur after THA using subject-specific finite element analyses. Twenty-four patients who had undergone primary cementless THA were analysed. BMD was measured using dual-energy X-ray absorptiometry (DEXA) at 1 week and 3, 6 and 12 months after THA. Seven regions of interest (ROIs) were defined in accordance with Gruen's system (ROIs 1–7). Computed tomography images of the femurs were acquired pre- and postoperatively, and the images were converted into three-dimensional finite element (FE) models. Equivalent stress and SED were analysed and compared with DEXA data. BMD was maintained 1 year after THA in ROIs 3, 4, 5 and 6, whereas BMD decreased in ROIs 1, 2 and 7. FE analysis revealed that equivalent stress in ROIs 3, 4, 5 and 6 was much higher than that in ROIs 1, 2 and 7. A significant correlation was observed between the rate of changes in BMD and equivalent stress. Reduction of equivalent stress may contribute to decrease in BMD in the femur after THA.     

The effect of implant design of linked total elbow arthroplasty on stability and stress: a finite element analysis

Several linked total elbow arthroplasty designs exist, which function similar to a loose hinge joint. Constraint behaviour is an important design consideration, as it affects joint stability, or how much secondary [e.g. varus–valgus (VV)] motion is permitted. Implant durability is also a concern, as bearing failures have been reported. This finite element analysis investigates the constraint characteristics and ultra high molecular weight polyethylene bearing stresses of three linked elbow design concepts [cylindrical (CY), hourglass (HG) and concave cylinder (CC)]. The bearing of the CY design was subjected to elevated Von Mises stresses (2.1–5.4 times higher than the HG and CC designs) due to edge loading. The HG design maintained low stresses, but was unable to provide consistent VV stability. The CC design also maintained low stresses while providing consistent VV stability. These results suggest that CC designs may provide better stability characteristics and durability in vivo, compared to the other two designs.     

Statistical factorial analysis approach for parameter calibration on material nonlinearity of intervertebral disc finite element model

In the biomechanics field, material parameters calibration is significant for finite element (FE) model to ensure a legit estimation of biomechanical response. Determining an appropriate combination of calibration factors is challenging as each constitutive component responds differently. This study proposes a statistical factorial analysis approach using L₁₆(4⁵) orthogonal array to evaluate material nonlinearity and applicable calibration factor of the intervertebral disc FE model in pure moment. The calibrated model exhibits improved agreement to the experimental findings for all directions. Appropriate combination of calibration parameter reduces the estimation gap to the experimental findings, ensuring agreeable biomechanical responses.     

Quantitative comparison of ligament formulation and pre-strain in finite element analysis of the human lumbar spine

Data has been published that quantifies the nonlinear, anisotropic material behaviour and pre-strain behaviour of the anterior longitudinal, supraspinous (SSL), and interspinous ligaments of the human lumbar spine. Additionally, data has been published on localized material properties of the SSL. These results have been incrementally incorporated into a previously validated finite element model of the human lumbar spine. Results suggest that the effects of increased ligament model fidelity on bone strain energy were moderate and the effects on disc pressure were slight, and do not justify a change in modelling strategy for most clinical applications. There were significant effects on the ligament stresses of the ligaments that were directly modified, suggesting that these phenomena should be included in FE models where ligament stresses are the desired metric.     

Investigation of effective intrusion and extrusion force for maxillary canine using finite element analysis

Abstract

Orthodontic tooth movement is mainly regulated by the biomechanical responses of loaded periodontal ligament (PDL). We investigated the effective intervals of orthodontic force in pure maxillary canine intrusion and extrusion referring to PDL hydrostatic stress and logarithmic strain. Finite element analysis (FEA) models, including a maxillary canine, PDL and alveolar bone, were constructed based on computed tomography (CT) images of a patient. The material properties of alveolar bone were non-uniformly defined using HU values of CT images; PDL was assumed to be a hyperelastic–viscoelastic material. The compressive stress and tensile stress ranging from 0.47 to 12.8?kPa and 18.8 to 51.2?kPa, respectively, were identified as effective for tooth movement; a strain 0.24% was identified as the lower limit of effective strain. The stress/strain distributions within PDL were acquired in canine intrusion and extrusion using FEA; root apex was the main force-bearing area in intrusion–extrusion movements and was more prone to resorption. Owing to the distinction of PDL biomechanical responses to compression and tension, the effective interval of orthodontic force was substantially lower in canine intrusion (80–90?g) than in canine extrusion (230–260?g). A larger magnitude of force remained applicable in canine extrusion. This study revised and complemented orthodontic biomechanical behaviours of tooth movement with intrusive–extrusive force and could further help optimize orthodontic treatment.