Validation of computational models in biomechanics

Influence of landmark and surgical variability on virtual assessment of total knee arthroplasty

Given advances in recent years in imaging modalities and computational hardware/software, virtual analyses are increasingly valuable and practical for evaluating total knee arthroplasty (TKA). However, the influence of variabilities at each step in computational analyses on predictions of TKA performance for a population has not yet been thoroughly investigated, nor the relationship between these variabilities and expected variations in surgical practice. Understanding these influences is nevertheless essential for ensuring the clinical relevance of theoretical predictions. Here, a morphological analysis of proximal tibial resections within TKA is proposed and investigated. The goals of this analysis are to quantify the influence of variability in landmark detection on resection parameters and to evaluate this sensitivity relative to expected clinical variability in TKA resections. Results here are directly applicable to population-level computational analyses of morphological and functional TKA performance.

Modeling stented coronary arteries: where we are, where to go

In the last two decades, numerical models have become well-recognized and widely adopted tools to investigate stenting procedures. Due to limited computational resources and modeling capabilities, early numerical studies only involved simplified cases and idealized stented arteries. Nowadays, increased computational power allows for numerical models to meet clinical needs and include more complex cases such as the implantation of multiple stents in bifurcations or curved vessels. Interesting progresses have been made in the numerical modeling of stenting procedures both from a structural and a fluid dynamics points of view. Moreover, in the drug eluting stents era, new insights on drug elution capabilities are becoming essential in the stent development. Lastly, image-based methods able to reconstruct realistic geometries from medical images have been proposed in the recent literature aiming to better describe the peculiar anatomical features of coronary vessels and increase the accuracy of the numerical models. In this light, this review provides a comprehensive analysis of the current state-of-the-art in this research area, discussing the main methodological advances and remarkable results drawn from a number of significant studies.

FEBio: finite elements for biomechanics

In the field of computational biomechanics, investigators have primarily used commercial software that is neither geared toward biological applications nor sufficiently flexible to follow the latest developments in the field. This lack of a tailored software environment has hampered research progress, as well as dissemination of models and results. To address these issues, we developed the FEBio software suite (http://mrl.sci.utah.edu/software/febio), a nonlinear implicit finite element (FE) framework, designed specifically for analysis in computational solid biomechanics. This paper provides an overview of the theoretical basis of FEBio and its main features. FEBio offers modeling scenarios, constitutive models, and boundary conditions, which are relevant to numerous applications in biomechanics. The open-source FEBio software is written in C++, with particular attention to scalar and parallel performance on modern computer architectures. Software verification is a large part of the development and maintenance of FEBio, and to demonstrate the general approach, the description and results of several problems from the FEBio Verification Suite are presented and compared to analytical solutions or results from other established and verified FE codes. An additional simulation is described that illustrates the application of FEBio to a research problem in biomechanics. Together with the pre- and postprocessing software PREVIEW and POSTVIEW, FEBio provides a tailored solution for research and development in computational biomechanics.

Multiscale mechanics of articular cartilage: potentials and challenges of coupling musculoskeletal, joint, and microscale computational models

Articular cartilage experiences significant mechanical loads during daily activities. Healthy cartilage provides the capacity for load bearing and regulates the mechanobiological processes for tissue development, maintenance, and repair. Experimental studies at multiple scales have provided a fundamental understanding of macroscopic mechanical function, evaluation of the micromechanical environment of chondrocytes, and the foundations for mechanobiological response. In addition, computational models of cartilage have offered a concise description of experimental data at many spatial levels under healthy and diseased conditions, and have served to generate hypotheses for the mechanical and biological function. Further, modeling and simulation provides a platform for predictive risk assessment, management of dysfunction, as well as a means to relate multiple spatial scales. Simulation-based investigation of cartilage comes with many challenges including both the computational burden and often insufficient availability of data for model development and validation. This review outlines recent modeling and simulation approaches to understand cartilage function from a mechanical systems perspective, and illustrates pathways to associate mechanics with biological function. Computational representations at single scales are provided from the body down to the microstructure, along with attempts to explore multiscale mechanisms of load sharing that dictate the mechanical environment of the cartilage and chondrocytes.

Validation of radiocarpal joint contact models based on images from a clinical MRI scanner

This study was undertaken to assess magnetic resonance imaging (MRI)-based radiocarpal surface contact models of functional loading in a clinical MRI scanner for future in vivo studies, by comparison with experimental measures from three cadaver forearm specimens. Experimental data were acquired using a Tekscan sensor during simulated light grasp. Magnetic resonance (MR) images were used to obtain model geometry and kinematics (image registration). Peak contact pressures (PPs) and average contact pressures (APs), contact forces and contact areas were determined in the radiolunate and radioscaphoid joints. Contact area was also measured directly from MR images acquired with load and compared with model data. Based on the validation criteria (within 25% of experimental data), out of the six articulations (three specimens with two articulations each), two met the criterion for AP (0%, 14%); one for peak pressure (20%); one for contact force (5%); four for contact area with respect to experiment (8%, 13%, 19% and 23%), and three contact areas met the criterion with respect to direct measurements (14%, 21% and 21%). Absolute differences between model and experimental PPs were reasonably low (within 2.5 MPa). Overall, the results indicate that MRI-based models generated from 3T clinical MR scanner appear sufficient to obtain clinically relevant data.

Application and evaluation of biomechanical models and scores for the planning of total hip arthroplasty

Intimate knowledge of the biomechanics of a given individual hip joint provides a potential advantage during the planning of total hip arthroplasty, and would thus have a positive influence over the outcome of such an intervention. In current clinical practise, the surgical planning is based solely on the status of the individual hip and its radiographic appearance. However, additional information could be gathered from the radiography to be used as input data for biomechanical models aimed at calculating the resultant force FR within the hip joint.

An investigation of the biomechanical models by Pauwels, Debrunner and Iglič was performed, where the magnitude of FR calculated by the models showed a favourable comparison to the in-vivo data from instrumented prostheses by Bergmann. The Blumentritt model returned abnormally high results. The computational results showed large variations for FR orientation, which tends to depend more on the model used than on patient-specific parameters.

Furthermore, a discrepancy was found between the data gathered from instrumented prostheses and the Standing Human Model within the ‘AnyBody Modeling System™’ software by AnyBody Tech. Additionally, the variations in inter-rater and intra-rater errors made while localizing radiographic landmarks were analysed with respect to their influence on Babisch-Layher-Blumentritt (BLB)-scoring using the Blumentritt hip model.

Analysis of several biomechanical models for the simulation of lamb liver behaviour using similarity coefficients from medical image

In this study, six biomechanical models for simulating lamb liver behaviour are presented. They are validated using similarity coefficients from Medical Image on reconstructed volumes from computerised tomography images. In particular, the Jaccard and Hausdorff coefficients are used. Loads of 20 and 40 g are applied to the livers and their deformation is simulated by means of the finite element method. The models used are a linear elastic model, a neo-Hookean model, a Mooney-Rivlin model, an Ogden model, a linear viscoelastic model and a viscohyperelastic model. The model that provided a behaviour that is closest to reality was the viscohyperelastic model, where the hyperelastic part was modelled with an Ogden model.

MRI-based modeling for radiocarpal joint mechanics: validation criteria and results for four specimen-specific models

The objective of this study was to validate the MRI-based joint contact modeling methodology in the radiocarpal joints by comparison of model results with invasive specimen-specific radiocarpal contact measurements from four cadaver experiments. We used a single validation criterion for multiple outcome measures to characterize the utility and overall validity of the modeling approach. For each experiment, a Pressurex film and a Tekscan sensor were sequentially placed into the radiocarpal joints during simulated grasp. Computer models were constructed based on MRI visualization of the cadaver specimens without load. Images were also acquired during the loaded configuration used with the direct experimental measurements. Geometric surface models of the radius, scaphoid and lunate (including cartilage) were constructed from the images acquired without the load. The carpal bone motions from the unloaded state to the loaded state were determined using a series of 3D image registrations. Cartilage thickness was assumed uniform at 1.0 mm with an effective compressive modulus of 4 MPa. Validation was based on experimental versus model contact area, contact force, average contact pressure and peak contact pressure for the radioscaphoid and radiolunate articulations. Contact area was also measured directly from images acquired under load and compared to the experimental and model data. Qualitatively, there was good correspondence between the MRI-based model data and experimental data, with consistent relative size, shape and location of radioscaphoid and radiolunate contact regions. Quantitative data from the model generally compared well with the experimental data for all specimens. Contact area from the MRI-based model was very similar to the contact area measured directly from the images. For all outcome measures except average and peak pressures, at least two specimen models met the validation criteria with respect to experimental measurements for both articulations. Only the model for one specimen met the validation criteria for average and peak pressure of both articulations; however the experimental measures for peak pressure also exhibited high variability. MRI-based modeling can reliably be used for evaluating the contact area and contact force with similar confidence as in currently available experimental techniques. Average contact pressure, and peak contact pressure were more variable from all measurement techniques, and these measures from MRI-based modeling should be used with some caution.

On validation of multibody musculoskeletal models

We review the opportunities to validate multibody musculoskeletal models in view of the current transition of musculoskeletal modelling from a research topic to a practical simulation tool in product design, healthcare and other important applications. This transition creates a new need for justification that the models are adequate representations of the systems they simulate. The need for a consistent terminology and established standards is identified and knowledge from fields with a more progressed state-of-the-art in verification and validation is introduced. A number of practical steps for improvement of the validation of multibody musculoskeletal models are pointed out and directions for future research in the field are proposed. It is hoped that a more structured approach to model validation can help to improve the credibility of musculoskeletal models.

Considerations for reporting finite element analysis studies in biomechanics

Simulation-based medicine and the development of complex computer models of biological structures is becoming ubiquitous for advancing biomedical engineering and clinical research. Finite element analysis (FEA) has been widely used in the last few decades to understand and predict biomechanical phenomena. Modeling and simulation approaches in biomechanics are highly interdisciplinary, involving novice and skilled developers in all areas of biomedical engineering and biology. While recent advances in model development and simulation platforms offer a wide range of tools to investigators, the decision making process during modeling and simulation has become more opaque. Hence, reliability of such models used for medical decision making and for driving multiscale analysis comes into question. Establishing guidelines for model development and dissemination is a daunting task, particularly with the complex and convoluted models used in FEA. Nonetheless, if better reporting can be established, researchers will have a better understanding of a model's value and the potential for reusability through sharing will be bolstered. Thus, the goal of this document is to identify resources and considerate reporting parameters for FEA studies in biomechanics. These entail various levels of reporting parameters for model identification, model structure, simulation structure, verification, validation, and availability. While we recognize that it may not be possible to provide and detail all of the reporting considerations presented, it is possible to establish a level of confidence with selective use of these parameters. More detailed reporting, however, can establish an explicit outline of the decision-making process in simulation-based analysis for enhanced reproducibility, reusability, and sharing.

Mechanical characterization of fourth generation composite humerus

Mechanical data on upper extremity surrogate bones, supporting use as biomechanical tools, is limited. The objective of this study was to characterize the structural behaviour of the fourth-generation composite humerus under simulated physiologic bending, specifically, stiffness, rigidity, and mid-diaphysial surface strains. Three humeri were tested in four-point bending, in anatomically defined anteroposterior (AP) and mediolateral (ML) planes. Stiffness and rigidity were derived using load–displacement data. Principal strains were determined at the anterior, posterior, medial, and lateral surfaces in the humeral mid-diaphysial transverse plane of one specimen using stacked rosettes. Linear structural behaviour was observed within the test range. Average stiffness and rigidity were greater in the ML (918 ± 18 N/mm; 98.4 ± 1.9 Nm2) than the AP plane (833 ± 16 N/mm; 89.3 ± 1.6 Nm2), with little inter-specimen variability. The ML/AP rigidity ratio was 1.1. Surface principal strains were similar at the anterior (5.41 µε/N) and posterior (5.43 µε/N) gauges for AP bending, and comparatively less for ML bending, i.e. 5.1 and 4.5 µε/N, at the medial and lateral gauges, respectively. This study provides novel strain and stiffness data for the fourth-generation composite humerus and also adds to published construct rigidity data. The presented results support the use of this composite bone as a tool for modelling and experimentation.

3D video-based deformation measurement of the pelvis bone under dynamic cyclic loading

Background

Dynamic three-dimensional (3D) deformation of the pelvic bones is a crucial factor in the successful design and longevity of complex orthopaedic oncological implants. The current solutions are often not very promising for the patient; thus it would be interesting to measure the dynamic 3D-deformation of the whole pelvic bone in order to get a more realistic dataset for a better implant design. Therefore we hypothesis if it would be possible to combine a material testing machine with a 3D video motion capturing system, used in clinical gait analysis, to measure the sub millimetre deformation of a whole pelvis specimen.

Method

A pelvis specimen was placed in a standing position on a material testing machine. Passive reflective markers, traceable by the 3D video motion capturing system, were fixed to the bony surface of the pelvis specimen. While applying a dynamic sinusoidal load the 3D-movement of the markers was recorded by the cameras and afterwards the 3D-deformation of the pelvis specimen was computed. The accuracy of the 3D-movement of the markers was verified with 3D-displacement curve with a step function using a manual driven 3D micro-motion-stage.

Results

The resulting accuracy of the measurement system depended on the number of cameras tracking a marker. The noise level for a marker seen by two cameras was during the stationary phase of the calibration procedure ± 0.036 mm, and ± 0.022 mm if tracked by 6 cameras. The detectable 3D-movement performed by the 3D-micro-motion-stage was smaller than the noise level of the 3D-video motion capturing system. Therefore the limiting factor of the setup was the noise level, which resulted in a measurement accuracy for the dynamic test setup of ± 0.036 mm.

Conclusion

This 3D test setup opens new possibilities in dynamic testing of wide range materials, like anatomical specimens, biomaterials, and its combinations. The resulting 3D-deformation dataset can be used for a better estimation of material characteristics of the underlying structures. This is an important factor in a reliable biomechanical modelling and simulation as well as in a successful design of complex implants.