The role of parameter identification in finite element contact analyses with reference to orthopaedic biomechanics applications

Developing and Validating a Model of Humeral Stem Primary Stability, Intended for In Silico Clinical Trials

In silico clinical trials (ISCT) can contribute to demonstrating a device's performance via credible computational models applied on virtual cohorts. Our purpose was to establish the credibility of a model for assessing the risk of humeral stem loosening in total shoulder arthroplasty, based on a twofold validation scheme involving both benchtop and clinical validation activities, for ISCT applications. A finite element model computing bone-implant micromotion (benchtop model) was quantitatively compared to a bone foam micromotion test (benchtop comparator) to ensure that the physics of the system was captured correctly. The model was expanded to a population-based approach (clinical model) and qualitatively evaluated based on its ability to replicate findings from a published clinical study (clinical comparator), namely that grit-blasted stems are at a significantly higher risk of loosening than porous-coated stems, to ensure that clinical performance of the stem can be predicted appropriately. Model form sensitivities pertaining to surgical variation and implant design were evaluated. The model replicated benchtop micromotion measurements (52.1 ± 4.3 µm), without a significant impact of the press-fit ("Press-fit": 54.0 ± 8.5 µm, "No press-fit": 56.0 ± 12.0 µm). Applied to a virtual population, the grit-blasted stems (227 ± 78µm) experienced significantly larger micromotions than porous-coated stems (162 ± 69µm), in accordance with the findings of the clinical comparator. This work provides a concrete example for evaluating the credibility of an ISCT study. By validating the modeling approach against both benchtop and clinical data, model credibility is established for an ISCT application aiming to enrich clinical data in a regulatory submission.

Finite element analysis in the optimization of posterior-stabilized total knee arthroplasty

Posterior-stabilized total knee arthroplasty (PS-TKA) is associated with high rates of satisfaction and functional recovery. This is notably attributed to implant optimization in terms of design, choice of materials, positioning and understanding of biomechanics. Finite elements analysis (FEA) is an assessment technique that contributed to this optimization by ensuring mechanical results based on numerical simulation. By close teamwork between surgeons, researchers and engineers, FEA enabled testing of certain clinical impressions. However, the methodological features of the technique led to wide variations in the presentation and interpretation of results, requiring a certain understanding of numerical and biomechanical fields by the orthopedic community. The present study provides an up-to-date review, aiming to address the following questions: what are the principles of FEA? What is the role of FEA in studying PS design in TKA? What are the key elements in the literature for understanding the role of FEA in PS-TKA? What is the contribution of FEA for understanding of tibiofemoral and patellofemoral biomechanical behavior? What are the limitations and perspectives of digital simulation and FEA in routine practice, with a particular emphasis on the "digital twin" concept? LEVEL OF EVIDENCE: V, expert opinion.

A Parameter Sensitivity Analysis on Multiple Finite Element Knee Joint Models

The reproducibility of computational knee joint modeling is questionable, with models varying depending on the modeling team. The influence of model variations on simulation outcomes should be investigated, since knowing the sensitivity of the model outcomes to model parameters could help determine which parameters to calibrate and which parameters could potentially be standardized, improving model reproducibility. Previous sensitivity analyses on finite element knee joint models have typically used one model, with a few parameters and ligaments represented as line segments. In this study, a parameter sensitivity analysis was performed using multiple finite element knee joint models with continuum ligament representations. Four previously developed and calibrated models of the tibiofemoral joint were used. Parameters of the ligament and meniscus material models, the cartilage contact formulation, the simulation control and the rigid cylindrical joints were studied. Varus-valgus simulations were performed, changing one parameter at a time. The sensitivity on model convergence, valgus kinematics, articulating cartilage contact pressure and contact pressure location were investigated. A scoring system was defined to categorize the parameters as having a “large,” “medium” or “small” influence on model output. Model outcomes were sensitive to the ligament prestretch factor, Young’s modulus and attachment condition parameters. Changes in the meniscus horn stiffness had a “small” influence. Of the cartilage contact parameters, the penalty factor and Augmented Lagrangian setting had a “large” influence on the cartilage contact pressure. In the rigid cylindrical joint, the largest influence on the outcome parameters was found by the moment penalty parameter, which caused convergence issues. The force penalty and gap tolerance had a “small” influence at most. For the majority of parameters, the sensitivity was model-dependent. For example, only two models showed convergence issues when changing the Quasi-Newton update method. Due to the sensitivity of the model parameters being model-specific, the sensitivity of the parameters found in one model cannot be assumed to be the same in other models. The sensitivity of the model outcomes to ligament material properties confirms that calibration of these parameters is critical and using literature values may not be appropriate.

Mechanical performance comparison of two surgical constructs for wrist four-corner arthrodesis via dorsal and radial approaches

BACKGROUND

Four-corner arthrodesis, which involves fusing four carpal bones while removing the scaphoid bone, is a standard surgery for the treatment of advanced stages of wrist arthritis. Nowadays, it can be performed using a dorsal approach by fixing a plate to the bones and a new radial approach is in development. To date, there is no consensus on the biomechanically optimal and most reliable surgical construct for four-corner arthrodesis.

METHODS

To evaluate them biomechanically and thus assist the surgeon in choosing the best implant orientation, radial or dorsal, the two different four-corner arthrodesis surgical constructs were virtually simulated on a 3D finite element model representing all major structures of the wrist. Two different realistic load sets were applied to the model, representing common tasks for the elderly.

FINDINGS

Results consistency was assessed by comparing with the literature the force magnitude computed on the carpal bones. The Von Mises stress distribution in the radial and dorsal plates were calculated. Stress concentration was located at the plate-screw interface for both surgical constructs, with a maximum stress value of 413 MPa for the dorsal plate compared to 326 MPa for the radial plate, meaning that the stress levels are more unfavourable in the dorsal approach.

INTERPRETATION

Although some bending stress was found in one load case, the radial plate was mechanically more robust in the other load case. Despite some limitations, this study provides, for the first time, quantified evidence that the newly developed radial surgical construct is mechanically as efficient as the dorsal surgical construct.

Automated Pipeline to Generate Anatomically Accurate Patient-Specific Biomechanical Models of Healthy and Pathological FSUs

State-of-the-art preoperative biomechanical analysis for the planning of spinal surgery not only requires the generation of three-dimensional patient-specific models but also the accurate biomechanical representation of vertebral joints. The benefits offered by computational models suitable for such purposes are still outweighed by the time and effort required for their generation, thus compromising their applicability in a clinical environment. In this work, we aim to ease the integration of computerized methods into patient-specific planning of spinal surgery. We present the first pipeline combining deep learning and finite element methods that allows a completely automated model generation of functional spine units (FSUs) of the lumbar spine for patient-specific FE simulations (FEBio). The pipeline consists of three steps: (a) multiclass segmentation of cropped 3D CT images containing lumbar vertebrae using the DenseVNet network, (b) automatic landmark-based mesh fitting of statistical shape models onto 3D semantic segmented meshes of the vertebral models, and (c) automatic generation of patient-specific FE models of lumbar segments for the simulation of flexion-extension, lateral bending, and axial rotation movements. The automatic segmentation of FSUs was evaluated against the gold standard (manual segmentation) using 10-fold cross-validation. The obtained Dice coefficient was 93.7% on average, with a mean surface distance of 0.88 mm and a mean Hausdorff distance of 11.16 mm (N = 150). Automatic generation of finite element models to simulate the range of motion (ROM) was successfully performed for five healthy and five pathological FSUs. The results of the simulations were evaluated against the literature and showed comparable ROMs in both healthy and pathological cases, including the alteration of ROM typically observed in severely degenerated FSUs. The major intent of this work is to automate the creation of anatomically accurate patient-specific models by a single pipeline allowing functional modeling of spinal motion in healthy and pathological FSUs. Our approach reduces manual efforts to a minimum and the execution of the entire pipeline including simulations takes approximately 2 h. The automation, time-efficiency and robustness level of the pipeline represents a first step toward its clinical integration.

Validated Finite Element Models of Premolars: A Scoping Review

Finite element (FE) models are widely used to investigate the biomechanics of reconstructed premolars. However, parameter identification is a complex step because experimental validation cannot always be conducted. The aim of this study was to collect the experimentally validated FE models of premolars, extract their parameters, and discuss trends. A systematic review was performed following Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Records were identified in three electronic databases (MEDLINE [PubMed], Scopus, The Cochrane Library) by two independent reviewers. Twenty-seven parameters dealing with failure criteria, model construction, material laws, boundary conditions, and model validation were extracted from the included articles. From 1306 records, 214 were selected for eligibility and entirely read. Among them, 19 studies were included. A heterogeneity was observed for several parameters associated with failure criteria and model construction. Elasticity, linearity, and isotropy were more often chosen for dental and periodontal tissues with a Young’s modulus mostly set at 18–18.6 GPa for dentine. Loading was mainly simulated by an axial force, and FE models were mostly validated by in vitro tests evaluating tooth strains, but different conditions about experiment type, sample size, and tooth status (intact or restored) were reported. In conclusion, material laws identified herein could be applied to future premolar FE models. However, further investigations such as sensitivity analysis are required for several parameters to clarify their indication.

Predicting calvarial morphology in sagittal craniosynostosis

Early fusion of the sagittal suture is a clinical condition called, sagittal craniosynostosis. Calvarial reconstruction is the most common treatment option for this condition with a range of techniques being developed by different groups. Computer simulations have a huge potential to predict the calvarial growth and optimise the management of this condition. However, these models need to be validated. The aim of this study was to develop a validated patient-specific finite element model of a sagittal craniosynostosis. Here, the finite element method was used to predict the calvarial morphology of a patient based on its preoperative morphology and the planned surgical techniques. A series of sensitivity tests and hypothetical models were carried out and developed to understand the effect of various input parameters on the result. Sensitivity tests highlighted that the models are sensitive to the choice of input parameter. The hypothetical models highlighted the potential of the approach in testing different reconstruction techniques. The patient-specific model highlighted that a comparable pattern of calvarial morphology to the follow up CT data could be obtained. This study forms the foundation for further studies to use the approach described here to optimise the management of sagittal craniosynostosis.

A Systematic Review of Continuum Modeling of Skeletal Muscles: Current Trends, Limitations, and Recommendations

Finite elasticity theory has been commonly used to model skeletal muscle. A very large range of heterogeneous constitutive laws has been proposed. In this review, the most widely used continuum models of skeletal muscles were synthetized and discussed. Trends and limitations of these laws were highlighted to propose new recommendations for future researches. A systematic review process was performed using two reliable search engines as PubMed and ScienceDirect. 40 representative studies (13 passive muscle materials and 27 active muscle materials) were included into this review. Note that exclusion criteria include tendon models, analytical models, 1D geometrical models, supplement papers, and indexed conference papers. Trends of current skeletal muscle modeling relate to 3D accurate muscle representation, parameter identification in passive muscle modeling, and the integration of coupled biophysical phenomena. Parameter identification for active materials, assumed fiber distribution, data assumption, and model validation are current drawbacks. New recommendations deal with the incorporation of multimodal data derived from medical imaging, the integration of more biophysical phenomena, and model reproducibility. Accounting for data uncertainty in skeletal muscle modeling will be also a challenging issue. This review provides, for the first time, a holistic view of current continuum models of skeletal muscles to identify potential gaps of current models according to the physiology of skeletal muscle. This opens new avenues for improving skeletal muscle modeling in the framework of in silico medicine.

Development and in vitro validation of a simplified numerical model for the design of a biomimetic femoral stem

BACKGROUND

Dense and stiff metallic femoral stems implanted into femurs for total hip arthroplasties produce a stress shielding effect since they modify the original load sharing path in the bony structure. Consequently, in the long term, the strain adaptive nature of bones leads to bone resorption, implant loosening, and the need for arthroplasty revision. The design of new cementless femoral stems integrating open porous structures can reduce the global stiffness of the stems, allowing them a better match with that of bones and provide their firm fixation via bone ingrowth, and, thus reduce the risk of implantation failure.

METHODS

This paper aims to develop and validate a simplified numerical model of stress shielding, which calculates the levels of bone resorption or formation by comparing strain distributions on the surface of the intact and the implanted femurs subjected to a simulated biological loading. Two femoral stems produced by laser powder-bed fusion using Ti-6Al-4V alloy are employed: the first is fully dense, while the second features a diamond cubic lattice structure in its core. The validation consists of a comparison of the numerically calculated force-displacement diagrams, and displacement and strain fields with their experimental equivalents obtained using the digital image correlation technique.

RESULTS AND CONCLUSIONS

The numerical models showed reasonable agreement between the force-displacement diagrams. Also, satisfactory results for the correlation analyses of the total displacement and equivalent strain fields were obtained. The stress shielding effect of the implant was assessed by comparing the equivalent strain fields of the implanted and intact femurs. The results obtained predicted less bone resorption in the femur implanted with the porous stem than with its dense counterpart.

Total ankle replacement design and positioning affect implant-bone micromotion and bone strains

•

A finite element model was developed to calculate micromotion of ankle implants.

•

Both optimally-positioned and malpositioned cases were considered.

•

Fixation nearer to the joint line relying on plural pegs improved implant stability.

•

Gaps between the implant and bone greatly increased micromotion and bone strains.

Implant loosening – commonly linked with elevated initial micromotion – is the primary indication for total ankle replacement (TAR) revision. Finite element modelling has not been used to assess micromotion of TAR implants; additionally, the biomechanical consequences of TAR malpositioning – previously linked with higher failure rates – remain unexplored. The aim of this study was to estimate implant-bone micromotion and peri-implant bone strains for optimally positioned and malpositioned TAR prostheses, and thereby identify fixation features and malpositioning scenarios increasing the risk of loosening. Finite element models simulating three of the most commonly used TAR devices (BOX^®, Mobility^® and Salto^®) implanted into the tibia/talus and subjected to physiological loads were developed. Mobility and Salto demonstrated the largest micromotion of all tibial and talar components, respectively. Any malpositioning of the implant creating a gap between it and the bone resulted in a considerable increase in micromotion and bone strains. It was concluded that better primary stability can be achieved through fixation nearer to the joint line and/or while relying on more than a single peg. Incomplete seating on the bone may result in considerably elevated implant-bone micromotion and bone strains, thereby increasing the risk for TAR failure.

Effects of biomechanical properties of the bone-implant interface on dental implant stability: from in silico approaches to the patient's mouth

Dental implants have become a routinely used technique in dentistry for replacing teeth. However, risks of failure are still experienced and remain difficult to anticipate. Multiscale phenomena occurring around the implant interface determine the implant outcome. The aim of this review is to provide an understanding of the biomechanical behavior of the interface between a dental implant and the region of bone adjacent to it (the bone-implant interface) as a function of the interface's environment. First, we describe the determinants of implant stability in relation to the different multiscale simulation approaches used to model the evolution of the bone-implant interface. Then, we review the various aspects of osseointegration in relation to implant stability. Next, we describe the different approaches used in the literature to measure implant stability in vitro and in vivo. Last, we review various factors affecting the evolution of the bone-implant interface properties.

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

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

The effect of fracture stability on the performance of locking plate fixation in periprosthetic femoral fractures

Periprosthetic femoral fracture (PFF) fixation failures are still occurring. The effect of fracture stability and loading on PFF fixation has not been investigated and this is crucial for optimum management of PFF. Models of stable and unstable PPFs were developed and used to quantify the effect of fracture stability and loading in a single locking plate fixation. Stress on the plate was higher in the unstable compared to the stable fixation. In the case of unstable fractures, it is possible for a single locking plate fixation to provide the required mechanical environment for callus formation without significant risk of plate fracture, provided partial weight bearing is followed. In cases where partial weight bearing is unlikely, additional biological fixation could be considered.

Evaluation of a new approach for modelling the screw–bone interface in a locking plate fixation: A corroboration study

Computational modelling of the screw–bone interface in fracture fixation constructs is challenging. While incorporating screw threads would be a more realistic representation of the physics, this approach can be computationally expensive. Several studies have instead suppressed the threads and modelled the screw shaft with fixed conditions assumed at the screw–bone interface. This study assessed the sensitivity of the computational results to modelling approaches at the screw–bone interface. A new approach for modelling this interface was proposed, and it was tested on two locking screw designs in a diaphyseal bridge plating configuration. Computational models of locked plating and far cortical locking constructs were generated and compared to in vitro models described in prior literature to corroborate the outcomes. The new approach led to closer agreement between the computational and the experimental stiffness data, while the fixed approach led to overestimation of the stiffness predictions. Using the new approach, the pattern of load distribution and the magnitude of the axial forces, experienced by each screw, were compared between the locked plating and far cortical locking constructs. The computational models suggested that under more severe loading conditions, far cortical locking screws might be under higher risk of screw pull-out than the locking screws. The proposed approach for modelling the screw–bone interface can be applied to any fixation involved application of screws.

Reduced stress shielding with limited micromotions using a carbon fibre composite biomimetic hip stem: a finite element model

Total hip arthroplasty (THA) enjoys excellent rates of success in older patients, but younger patients are still at risk of aseptic loosening and bone resorption from stress shielding. One solution to the stress shielding problem is to use a hip stem with mechanical properties matching those of cortical bone. The objective of the present study was to investigate numerically the biomechanical performance of such a biomimetic hip stem based on a hydroxyapatite (HA)-coated carbon fibre composite. A finite element model (FEM) of the biomimetic stem was constructed. Contact elements were studied to model the bone–implant interface in a non-osseointegrated and osseointegrated state in the best way. Three static load cases representing slow walking, stair climbing, and gait in a healthy individual were considered. Stress shielding and bone–implant interface micromotions were evaluated and compared with the results of a similar FEM based on titanium alloy (Ti–6Al–4V). The composite stems allowed for reduced stress shielding when compared with a traditional Ti–6Al–4V stem. Micromotions were slightly higher with the composite stem, but remained below 40 μm on most of the HA-coated surface. It is concluded that a biomimetic composite stem might offer a better compromise between stress shielding and micromotions than the Ti–6Al–4V stem with the same external geometry.

A numerical study of failure mechanisms in the cemented resurfaced femur: effects of interface characteristics and bone remodelling

Failure mechanisms of the resurfaced femoral head include femoral neck fracture in the short-term and stress shielding and implant loosening in the long-term. In this study, finite element simulations of the resurfaced femur considering a debonded implant-cement interface, variable stem-bone interface conditions, and bone remodelling were used to study load transfer within the resurfaced femur and to investigate its relationship with known failure mechanisms. Realistic three-dimensional finite element models of an intact and resurfaced femur were used. Various conditions at the interface between the stem of the prosthesis and the bone were considered. Loading conditions included normal walking and stair climbing. For all stem-bone contact conditions, the tensile stresses in the cement mantle varied between 1 MPa and 5.4 MPa, except near the distal rim of the resurfacing component where they reached 5.4-7MPa. In the case of full stem-bone contact, high von Mises stresses (114-121MPa) were generated in the implant at the stem-cup junction. These stresses were considerably reduced (maximum von Mises stress, 76 MPa) where a gap was present at the stem-bone interface. Resurfacing led to strain shielding of the bone of the femoral head (20-75 per cent strain reductions) and periprosthetic bone resorption (50-80 per cent bone density reductions) for all interface stem-bone contact conditions. In the lateral femoral head and the proximal femoral shaft around the trochantric region, bone density reductions varied between 10 per cent and 50 per cent. Bone apposition was observed in the inferior-medial part of the femoral head and proximal femoral neck region. For full stem-bone contact, more load was transferred through the stem to the surrounding bone, exacerbating strain shielding. Although femoral hip resurfacing conserves bone stock at the primary operation, strain shielding and periprosthetic bone resorption might lead to eventual loosening over time. Post-operatively, the resurfacing procedure generated elevated strains (0.50-0.75 per cent strain) in the proximal femoral neck-component junction irrespective of the variation in interface conditions, indicating an initial risk of femoral neck fracture. Subsequent to bone remodelling, this strain concentration was considerably reduced (0.35-0.50 per cent strain), lowering the risk of neck fracture. In order to reduce the potential risk of neck fracture, patients should avoid activities which might induce high loading of the hip during the early post-operative period to allow the bone around the proximal femoral neck to remodel and heal.

Subject specific finite element analysis of implant stability for a cementless femoral stem

BACKGROUND

The primary stability of a cementless implant is crucial to ensure long term stability through osseointegration. In the present study we have examined how subject specific finite element models can be used to evaluate the stability of a cementless femoral stem.

METHODS

Micromotion on the bone-implant interface of a cementless stem was measured experimentally in six human cadaver femurs. Subject specific finite element models were built from computed tomography of the same femurs, and used to simulate the same load scenario used experimentally.

FINDINGS

Both experimental measurements and numerical analyses showed a tendency of increased rotational stability for bigger implants. Good correlation was found between measurements and calculated values of axial rotation (R(2)=0.74, P<0.001). The finite element models produced interface micromotion of the same magnitude as measured experimentally, with micromotion generally below 40 microm. Bigger femoral stems were found to decrease the micromotion in the experimental measurements. This tendency could not be recognised in the interface micromotion from the finite element models.

INTERPRETATION

The finite element models showed limited success in predicting interfacial micromotion, but reproduced a similar pattern of rotational stability for the implants as seen experimentally. Since rotation in retroversion is often the main concern when studying implant stability, subject specific finite element models could be employed for pre-clinical evaluation of implants.

Subject specific finite element analysis of stress shielding around a cementless femoral stem

BACKGROUND

Stress shielding around a femoral stem is usually assessed experimentally using composite or human cadaver femurs. In the present study we have explored the feasibility of using subject specific finite element models to determine stress shielding in operated femurs.

METHODS

Cortical strain was measured experimentally on seven human cadaver femurs, intact and implanted with a straight cementless prosthesis. Two load configurations were considered: single leg stance and stair climbing. Subject specific finite element models derived from computed tomography of the same femurs were analysed intact and with an implant. Principal cortical strain was used to validate the finite element models. Stress shielding was defined as the change in equivalent (von Mises) strain between pre- and post-operative femurs.

FINDINGS

Cortical strain predicted by the finite element analyses showed to be close to unity with the experimental observations for both intact (R2=0.94, slope=0.99), operated femurs (R2=0.86, slope=0.86) and stress shielding (R2=0.70, slope=0.90). In the proximal calcar area, the region most prone to periprosthetic remodelling, the finite element models were found to successfully reproduce the stress shielding observed experimentally.

INTERPRETATION

The study shows that subject specific finite element models manage to describe the stress shielding pattern measured in vitro in the different femurs. Finite element models based on actual human femurs (cadaver and/or patient) could thus be a useful tool in the pre-clinical evaluation of new implants.

Finite element modelling of primary hip stem stability: the effect of interference fit

The most commonly reported complications related to cementless hip stems are loosening and thigh pain; both of these have been attributed to high levels of relative micromotion at the bone-implant interface due to insufficient primary fixation. Primary fixation is believed by many to rely on achieving a sufficient interference fit between the implant and the bone. However, attempting to achieve a high interference fit not infrequently leads to femoral canal fracture either intra-operatively or soon after. The appropriate range of diametrical interference fit that ensures primary stability without risking femoral fracture is not well understood. In this study, a finite element model was constructed to predict micromotion and, therefore, instability of femoral stems. The model was correlated with an in vitro micromotion experiment carried out on four cadaver femurs. It was confirmed that interference fit has a very significant effect on micromotion and ignoring this parameter in an analysis of primary stability is likely to underestimate the stability of the stem. Furthermore, it was predicted that the optimal level of interference fit is around 50 microm as this is sufficient to achieve good primary fixation while having a safety factor of 2 against femoral canal fracture. This result is of clinical relevance as it indicates a recommendation for the surgeon to err on the side of a low interference fit rather than risking femoral fracture.

Verification, validation and sensitivity studies in computational biomechanics

Computational techniques and software for the analysis of problems in mechanics have naturally moved from their origins in the traditional engineering disciplines to the study of cell, tissue and organ biomechanics. Increasingly complex models have been developed to describe and predict the mechanical behavior of such biological systems. While the availability of advanced computational tools has led to exciting research advances in the field, the utility of these models is often the subject of criticism due to inadequate model verification and validation (V&V). The objective of this review is to present the concepts of verification, validation and sensitivity studies with regard to the construction, analysis and interpretation of models in computational biomechanics. Specific examples from the field are discussed. It is hoped that this review will serve as a guide to the use of V&V principles in the field of computational biomechanics, thereby improving the peer acceptance of studies that use computational modeling techniques.

Stress analysis of osseointegrated transfemoral prosthesis: a finite element model

In this study a three dimensional finite element model was established to analyze the interfacial stress distribution in the osseointegrated transfemoral prosthesis. The effects of implant length and screw pitch on the stress and strain distributions were analyzed in terms of bonded and contact bone-implant interfaces. The influences of the contact parameters, contact stiffness and friction coefficient, of the model were evaluated. The results show the differences on the stress distributions at the bone-implant interface in the models with different implant lengths and thread pitches. The highest stress in the model with contact interface is higher than that in the bonded model.

Predicting the subject-specific primary stability of cementless implants during pre-operative planning: Preliminary validation of subject-specific finite-element models

Pre-operative planning help the surgeon in taking the proper clinical decision. The ultimate goal of this work is to develop numerical models that allow the surgeon to estimate the primary stability during the pre-operative planning session. The present study was aimed to validate finite-element (FE) models accounting for patient and prosthetic size and position as planned by the surgeon. For this purpose, the FE model of a cadaveric femur was generated starting from the CT scan and the anatomical position of a cementless stem derived by a skilled surgeon using a pre-operative CT-based planning simulation software. In-vitro experimental measurements were used as benchmark problem to validate the bone-implant relative micromotions predicted by the patient-specific FE model. A maximum torque in internal rotation of 11.4 Nm was applied to the proximal part of the hip stem. The error on the maximum predicted micromotion was 12% of the peak micromotion measured experimentally. The average error over the entire range of applied torques was only 7% of peak measurement. Hence, the present study confirms that it is possible to accurately predict the level of primary stability achieved for cementless stems using numerical models that account for patient specificity and surgical variability.

Primary stability of an anatomical cementless hip stem: A statistical analysis

The primary stability that the surgeon can achieve during surgery is a determinant of the clinical success of cementless implants. Thus, estimating what level of primary stability can be obtained with a new design is an important aspect of pre-clinical evaluation. The primary stability of a cementless hip stem is not only affected by the implant design, but also by other factors such as the mechanical quality of the host bone, the presence of gaps around the bone-implant interface, the body weight of the patient, and the size of the implant. Even the most extensive experimental study can only explore a small sub-set of all possible combinations found in vivo. To overcome this limitation, we propose a combination of experimental and numerical methods. The primary stability of a cementless anatomical stem is assessed in vitro. A finite element model is developed to accurately replicate the same experiment. The model is then parameterised over the various factors that affect the primary stability, and used in a Monte Carlo scheme to assess the primary stability over a simulated population. In this study, the method was used to investigate the mechanical stability of an anatomical cementless stem over more than 1000 simulated cases. Twenty cases were found macroscopically unstable, due to a combination of unfavourable conditions. The rest of the Monte Carlo sample showed on average a peak micromotion under stair climbing loading of 206 +/- 159 microm. The proposed method can be used to evaluate new designs in conditions more representative of the variability in clinical practice.

A Circumferentially Flanged Tibial Tray Minimizes Bone-Tray Shear Micromotion

Aseptic loosening of the tibial component is the major complication of total knee arthroplasty. There is an association between early excessive shear micromotion between the bone and the tray of the tibial component and late aseptic loosening. Using non-linear finite element analysis, whether a tibial tray with a circumferentially flanged rim and a mating cut in the proximal tibia could minimize bone-tray shear micromotion was considered. fifteen competing tray designs with various degrees of flange curvature were assessed with the aim of minimizing bone-tray shear micromotion. A trade-off was found between reducing micromotion and increasing peripheral cancellous bone stresses. It was found that, within the limitations of the study, there was a theoretical design that could virtually eliminate micromotion due to axial loads, with minimal bone removal and without the use of screws or pegs.

Numerical model to predict the longterm mechanical stability of cementless orthopaedic implants

The objective of this research was to develop a purely biomechanical model, intended to predict the long-term secondary stability of the implant starting from the biomechanical stability immediately after the operation. A continuous rule-based adaptation scheme was formulated as a dynamic system, and the work verified if such a model produced unique and clinically meaningful solutions. It also investigated whether this continuous model provided results comparable with those of a simpler, discrete-states model used in a previous study. The proposed model showed stable convergence behaviour with all investigated initial conditions, with oscillatory behaviour limited to the first steps of the simulation. The results obtained with the wide range of initial conditions support the hypothesis of the existence and uniqueness of the solution for all initial conditions. The differences between the continuous model and the simpler and more efficient finite-states model were found to be extremely modest (less than 4% over the predicted bonded area). Because of these minimal differences, the use of the much faster finite-states model is recommended to investigate asymptotic conditions, and the continuous model described should be used to investigate the evolution over time of the adaptive process.

Shape optimization of a cementless hip stem for a minimum of interface stress and displacement

The primary stem stability is an essential factor for success of cementless hip stems. A correct choice of the stem geometry can improve the stem stability and, consequently, increase the life time of a hip implant. In this work, it is proposed a computational model for shape optimization of cementless hip stems. The optimization problem is formulated by the minimization of relative displacement and stress on bone/stem interface using a multi-criteria objective function. Also multiple loads are considered to incorporate several daily life activities. Design variables are parameters that characterize the geometry of selected cross sections, which are subject to geometric constraints to ensure a clinically admissible shape. The stem/bone set is considered a structure in equilibrium with contact conditions on interface. The contact formulation allows us to analyze different lengths of porous coating. The optimization problem is solved numerically by a steepest descent method. The interface stress and relative displacement are obtained solving the contact problem by the finite element method. Numerical examples are presented for a two-dimensional model of a hip stem, however, the formulation is general and can be applied to the three-dimensional case. The model gives indications about the relation between shape, porous coating and prosthesis stability.

Influence of the stiffness of bone defect implants on the mechanical conditions at the interface--a finite element analysis with contact

The study focused on the influence of the implant material stiffness on stress distribution and micromotion at the interface of bone defect implants. We hypothesized that a low-stiffness implant with a modulus closer to that of the surrounding trabecular bone would yield a more homogeneous stress distribution and less micromotion at the interface with the bony bed. To prove this hypothesis we generated a three-dimensional, non-linear, anisotropic finite element (FE) model. The FE model corresponded to a previously developed animal model in sheep. A prismatic implant filled a standardized defect in the load-bearing area of the trabecular bone beneath the tibial plateau. The interface was described by face-to-face contact elements, which allow press fits, friction, sliding, and gapping. We assumed a physiological load condition and calculated contact pressures, shear stresses, and shear movements at the interface for two implants of different stiffness (titanium: E=110GPa; composite: E=2.2GPa). The FE model showed that the stress distribution was more homogeneous for the low-stiffness implant. The maximum pressure for the composite implant (2.1 MPa) was lower than for the titanium implant (5.6 MPa). Contrary to our hypothesis, we found more micromotion for the composite (up to 6 microm) than for the titanium implant (up to 4.5 microm). However, for both implants peak stresses and micromotion were in a range that predicts adequate conditions for the osseointegration. This was confirmed by the histological results from the animal studies.

The primary stability of a cementless stem varies between subjects as much as between activities

The rehabilitation program adopted immediately after a cementless total hip replacement is a very important factor, because of the known relationship between osseointegration and implant micromotion. The present study was aimed to evaluate which type of task is the most critical in terms of bone-implant relative micromotion. Both inter-task and inter-subject variability were taken into account to verify if the movement strategy could be determinant on this assessment. A previously validated finite element model was used to predict the peak total micromovements over the entire bone-implant contact surface in four different patients, performing nine different tasks, using published data on joint forces recorded by instrumented hip prostheses. The results predicted by the various simulations suggest that while stair climbing is surely a critical task for primary stability, for some subjects other tasks may be as critical as stair climbing. From a variance analysis for simple crossover design on the predicted peak micromotion, the inter-subject variability had much more influence on the primary stability of cementless implant than the inter-task variability. Even if the results of Patient IBL, who was reported to have difficulties to perform any activities in a normal way, were excluded from the statistical analysis, the inter-subject variability remained still higher than the inter-task variability. The results obtained from simulations suggest that the strategy the hip replacement patient adopts to perform a given motor task, may be, for the implant stability, equally or even more critical than the type of motor task performed.