Patient-specific finite element modeling of bones

A comparative analysis of dimensionality reduction surrogate modeling techniques for full human body finite element impact simulations

Fast-running surrogate computational models (simpler computational models) have been successfully used to replace time-intensive finite element models. However, it is unclear how well they perform in accurately and efficiently replicating complex, full human body finite element models. Here we survey several surrogate modeling techniques and assess their accuracy in predicting full strain fields of tissues of interest during a highly dynamic behind armor blunt trauma impact to the liver. We found that coupling dimensionality reduction on the high-dimensional output space (principal component analysis or autoencoders) with machine learning techniques (Gaussian Process Regression or multi-output neural networks) provides a framework capable of accurately and efficiently replacing complex full human body models. It was found that these surrogate models can successfully predict the strain fields (<10% average strain error) of select tissues during a nonlinear impact event but careful consideration should be given to element parsing and modeling technique.

Finite element analysis of patient-specific additive-manufactured implants

Introduction: Bone tumors, characterized by diverse locations and shapes, often necessitate surgical excision followed by custom implant placement to facilitate targeted bone reconstruction. Leveraging additive manufacturing, patient-specific implants can be precisely tailored with complex geometries and desired stiffness, enhancing their suitability for bone ingrowth. Methods: In this work, a finite element model is employed to assess patient-specific lattice implants in femur bones. Our model is validated using experimental data obtained from an animal study (n = 9). Results: The results demonstrate the accuracy of the proposed finite element model in predicting the implant mechanical behavior. The model was used to investigate the influence of reducing the elastic modulus of a solid Ti6Al4V implant by tenfold, revealing that such a reduction had no significant impact on bone behavior under maximum compression and torsion loading. This finding suggests a potential avenue for reducing the endoprosthesis modulus without compromising bone integrity. Discussion: Our research suggests that employing fully lattice implants not only facilitates bone ingrowth but also has the potential to reduce overall implant stiffness. This reduction is crucial in preventing significant bone remodeling associated with stress shielding, a challenge often associated with the high stiffness of fully solid implants. The study highlights the mechanical benefits of utilizing lattice structures in implant design for enhanced patient outcomes.

The Potential of Deep Learning to Advance Clinical Applications of Computational Biomechanics

When combined with patient information provided by advanced imaging techniques, computational biomechanics can provide detailed patient-specific information about stresses and strains acting on tissues that can be useful in diagnosing and assessing treatments for diseases and injuries. This approach is most advanced in cardiovascular applications but can be applied to other tissues. The challenges for advancing computational biomechanics for real-time patient diagnostics and treatment include errors and missing information in the patient data, the large computational requirements for the numerical solutions to multiscale biomechanical equations, and the uncertainty over boundary conditions and constitutive relations. This review summarizes current efforts to use deep learning to address these challenges and integrate large data sets and computational methods to enable real-time clinical information. Examples are drawn from cardiovascular fluid mechanics, soft-tissue mechanics, and bone biomechanics. The application of deep-learning convolutional neural networks can reduce the time taken to complete image segmentation, and meshing and solution of finite element models, as well as improving the accuracy of inlet and outlet conditions. Such advances are likely to facilitate the adoption of these models to aid in the assessment of the severity of cardiovascular disease and the development of new surgical treatments.

Dual-zone material assignment method for correcting partial volume effects in image-based bone models

In image-based finite element analysis of bone, partial volume effects (PVEs) arise from image blur at tissue boundaries and as a byproduct of geometric reconstruction and meshing during model creation. In this study, we developed and validated a material assignment approach to mitigate partial volume effects. Our validation data consisted of physical torsion testing of intact tibiae from N = 20 Swiss alpine sheep. We created finite element models from micro-CT scans of these tibiae using three popular element types (10-node tetrahedral, 8-node hexahedral, and 20-node hexahedral). Without partial volume management, the models over-predicted the torsional rigidity compared to physical biomechanical tests. To address this problem, we implemented a dual-zone material model to treat elements that overlap low-density surface voxels as soft tissue rather than bone. After in situ inverse optimization, the dual-zone material model produced strong correlations and high absolute agreement between the virtual and physical tests. This suggests that with appropriate partial volume management, virtual mechanical testing can be a reliable surrogate for physical biomechanical testing. For maximum flexibility in partial volume management regardless of element type, we recommend the use of the following dual-zone material model for ovine tibiae: soft-tissue cutoff density of 665 mgHA/cm3 with a soft tissue modulus of 50 MPa (below cutoff) and a density-modulus conversion slope of 10,225 MPa-cm3/mgHA for bone (above cutoff).

A combined active shape and mean appearance model for the reconstruction of segmental bone loss

This study investigates the novel combination of an active shape and mean appearance model to estimate missing bone geometry and density distribution from sparse inputs simulating segmental bone loss of the femoral diaphysis. An active shape Gaussian Process Morphable model was trained on healthy right femurs of South African males to model shape. The density distribution was approximated based on the mean appearance of computed tomography images from the training set. Estimations of diaphyseal resections were obtained by probabilistic fitting of the active shape model to sparse inputs consisting of proximal and distal femoral data on computed tomography images. The resulting shape estimates of the diaphyseal resections were then used to map the mean appearance model to the patients' missing bone geometry, constructing density estimations. In this way, resected bone surfaces were estimated with an average error of 2.24 (0.5) mm. Density distributions were approximated within 87 (0.7) % of the intensity of the original target images before the simulated segmental bone loss. These results fall within the acceptable tolerances required for surgical planning and reconstruction of long bone defects.

A Density-Dependent Target Stimulus for Inverse Bone (Re)modeling with Homogenized Finite Element Models

Inverse bone (re)modeling (IBR) can infer physiological loading conditions from the bone microstructure. IBR scales unit loads, imposed on finite element (FE) models of a bone, such that the trabecular microstructure is homogeneously loaded and the difference to a target stimulus is minimized. Micro-FE (µFE) analyses are typically used to model the microstructure, but computationally more efficient, homogenized FE (hFE) models, where the microstructure is replaced by an equivalent continuum, could be used instead. However, also the target stimulus has to be translated from the tissue to the continuum level. In this study, a new continuum-level target stimulus relating relative bone density and strain energy density is proposed. It was applied using different types of hFE models to predict the physiological loading of 21 distal radii sections, which was subsequently compared to µFE-based IBR. The hFE models were able to correctly identify the dominant load direction and showed a high correlation of the predicted forces, but mean magnitude errors ranged from − 14.7 to 26.6% even for the best models. While µFE-based IBR can still be regarded as a gold standard, hFE-based IBR enables faster predictions, the usage of more sophisticated boundary conditions, and the usage of clinical images.

Computer-assisted preoperative planning of bone fracture fixation surgery: A state-of-the-art review

Background: Bone fracture fixation surgery is one of the most commonly performed surgical procedures in the orthopedic field. However, fracture healing complications occur frequently, and the choice of the most optimal surgical approach often remains challenging. In the last years, computational tools have been developed with the aim to assist preoperative planning procedures of bone fracture fixation surgery.

Objectives: The aims of this review are 1) to provide a comprehensive overview of the state-of-the-art in computer-assisted preoperative planning of bone fracture fixation surgery, 2) to assess the clinical feasibility of the existing virtual planning approaches, and 3) to assess their clinical efficacy in terms of clinical outcomes as compared to conventional planning methods.

Methods: A literature search was performed in the MEDLINE-PubMed, Ovid-EMBASE, Ovid-EMCARE, Web of Science, and Cochrane libraries to identify articles reporting on the clinical use of computer-assisted preoperative planning of bone fracture fixation.

Results: 79 articles were included to provide an overview of the state-of-the art in virtual planning. While patient-specific geometrical model construction, virtual bone fracture reduction, and virtual fixation planning are routinely applied in virtual planning, biomechanical analysis is rarely included in the planning framework. 21 of the included studies were used to assess the feasibility and efficacy of computer-assisted planning methods. The reported total mean planning duration ranged from 22 to 258 min in different studies. Computer-assisted planning resulted in reduced operation time (Standardized Mean Difference (SMD): -2.19; 95% Confidence Interval (CI): -2.87, -1.50), less blood loss (SMD: -1.99; 95% CI: -2.75, -1.24), decreased frequency of fluoroscopy (SMD: -2.18; 95% CI: -2.74, -1.61), shortened fracture healing times (SMD: -0.51; 95% CI: -0.97, -0.05) and less postoperative complications (Risk Ratio (RR): 0.64, 95% CI: 0.46, 0.90). No significant differences were found in hospitalization duration. Some studies reported improvements in reduction quality and functional outcomes but these results were not pooled for meta-analysis, since the reported outcome measures were too heterogeneous.

Conclusion: Current computer-assisted planning approaches are feasible to be used in clinical practice and have been shown to improve clinical outcomes. Including biomechanical analysis into the framework has the potential to further improve clinical outcome.

Model Construction for the Estimation of Healthy Bone Shape and Density Distribution

Statistical models are widely used within biomedical fields for automated segmentation and reconstruction of healthy geometry. In the absence of contralateral geometry, statistical models are a viable alternative for reconstructing healthy bone anatomy. Therefore, statistical models of shape and appearance were constructed from sample data based on the right femur of South African males, and their use in an automated segmentation and density estimation application was investigated. The models reproduced the shape and density distribution of the population with an average error of 1.3 mm and a 90% density fit. These results fall within the acceptable tolerance limits of reconstructive surgery and appear promising for practical use in implant design.Clinical Relevance- Constructing and validating statistical models and registration algorithms provides the groundwork for further investigation into automating the digital reconstruction of pathological bone for use in implant design.

A computational framework for crack propagation in spatially heterogeneous materials

This paper presents a mathematical formulation and numerical modelling framework for brittle crack propagation in heterogeneous elastic solids. Such materials are present in both natural and engineered scenarios. The formulation is developed in the framework of configurational mechanics and solved numerically using the finite-element method. We show the methodology previously established for homogeneous materials without the need for any further assumptions. The proposed model is based on the assumption of maximal dissipation of energy and uses the Griffith criterion; we show that this is sufficient to predict crack propagation in brittle heterogeneous materials, with spatially varying Young's modulus and fracture energy. Furthermore, we show that the crack path trajectory orientates itself such that it is always subject to Mode-I. The configurational forces and fracture energy release rate are both expressed exclusively in terms of nodal quantities, avoiding the need for post-processing and enabling a fully implicit formulation for modelling the evolving crack front and creation of new crack surfaces. The proposed formulation is verified and validated by comparing numerical results with both analytical solutions and experimental results. Both the predicted crack path and load-displacement response show very good agreement with experiments where the crack path was independent of material heterogeneity for those cases. Finally, the model is successfully used to consider the real and challenging scenario of fracture of an equine bone, with spatially varying material properties obtained from CT scanning. This article is part of a discussion meeting issue 'A cracking approach to inventing new tough materials: fracture stranger than friction'.

Finite element analysis of different locking plate fixation methods for the treatment of ulnar head fracture

Background

Ulnar head fractures are increasingly higher with the growing proportion of the elderly people. Failure to achieve a stable anatomic reduction of ulna head fracture may lead to a distal radioulnar joint (DRUJ) dysfunction and nonunion of the distal radius. Due to the lack of the postoperative reporting outcomes and the biomechanical studies, it has not been well established about the optimal management of the comminuted distal ulna head fracture. Hence, the purpose of this study is to use finite element analysis to explain the advantages and disadvantages of ulnar-side locking plate fixation compared with dorsal-side locking plate fixation and its screw arrangement in the treatment of ulnar head fractures.

Methods

FE models of the ulnar head fracture and the models of ulnar-side locking plate and dorsal-side plate with two or three distal screws was constructed. In order to simulate forces acting on the ulnar and the osteosynthesis material during daily-life activity in subjects who underwent reconstructive surgery, we applied three loading conditions to each model, viz. 20 N axial compression, 50 N axial compression, 1 N∙m torsion moment, 1 N∙m lateral bending moments, and 1 N∙m extension bending moments. Under these conditions, values of the von Mises stress (VMS) distribution of the implant, peak VMS, the relative displacement of the head and shaft fragments between the fracture ends and the displacement and its direction of the models were investigated.

Results

The stress values of ulnar-side plates were lower than those of dorsal-side plates. And the ulnar-plate fixation system also has smaller maximum displacement and relative displacement. When adding a screw in the middle hole of the ulnar head, the values of model displacement and the peak stress in fixation system are lower, but it may evidently concentrate the stress on the middle screw.

Conclusions

In conclusion, our study indicated that ulnar-side locking plates resulted in a lower stress distribution in the plate and better stability than dorsal-side locking plates for ulnar head fracture fixation. Adding an additional screw to the ulnar head could increase the stability of the fixation system and provide an anti-torsion function. This study requires clinical confirmation of its practicality in the treatment of ulnar head fractures. This study requires clinical confirmation as to its practicality in the treatment of ulnar head fracture.

A GPU-based caching strategy for multi-material linear elastic FEM on regular grids

In this study, we present a novel strategy to the method of finite elements (FEM) of linear elastic problems of very high resolution on graphic processing units (GPU). The approach exploits regularities in the system matrix that occur in regular hexahedral grids to achieve cache-friendly matrix-free FEM. The node-by-node method lies in the class of block-iterative Gauss-Seidel multigrid solvers. Our method significantly improves convergence times in cases where an ordered distribution of distinct materials is present in the dataset. The method was evaluated on three real world datasets: An aluminum-silicon (AlSi) alloy and a dual phase steel material sample, both captured by scanning electron tomography, and a clinical computed tomography (CT) scan of a tibia. The caching scheme leads to a speed-up factor of ×2-×4 compared to the same code without the caching scheme. Additionally, it facilitates the computation of high-resolution problems that cannot be computed otherwise due to memory consumption.

The biomechanical evaluation of magnetic forces to drive osteogenesis in newborn's with cleft lip and palate

This study examined the potential for dental magnets to act as a driving force for osteogenesis in the palate of newborns with a unilateral cleft lip and palate. In the first part of the study dental magnets were arranged in a set up mimicking a distraction device and the curves of the magnetic attraction force versus gap distance curves generated, with and without the presence of palatal rugae tissue in between both sides of the distraction device. The attraction forces ranged from 1 to 12 N depending on the gap distance and the presence of soft tissue in the gap. In the second part of the study these forces were used as input for a 3D finite element model of the palate of a newborn affected by unilateral cleft lip and palate. In the analysis of load transfer, it was found that the strains generated by a magnetically induced distraction exceed 1,500 µstrain suggesting that bone locally is submitted to mild overload leading to bone apposition.

Image-based biomechanical models of the musculoskeletal system

Finite element modeling is a precious tool for the investigation of the biomechanics of the musculoskeletal system. A key element for the development of anatomically accurate, state-of-the art finite element models is medical imaging. Indeed, the workflow for the generation of a finite element model includes steps which require the availability of medical images of the subject of interest: segmentation, which is the assignment of each voxel of the images to a specific material such as bone and cartilage, allowing for a three-dimensional reconstruction of the anatomy; meshing, which is the creation of the computational mesh necessary for the approximation of the equations describing the physics of the problem; assignment of the material properties to the various parts of the model, which can be estimated for example from quantitative computed tomography for the bone tissue and with other techniques (elastography, T1rho, and T2 mapping from magnetic resonance imaging) for soft tissues. This paper presents a brief overview of the techniques used for image segmentation, meshing, and assessing the mechanical properties of biological tissues, with focus on finite element models of the musculoskeletal system. Both consolidated methods and recent advances such as those based on artificial intelligence are described.

Image-based finite-element modeling of the human femur

Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.

Empirical relationships between bone density and ultimate strength: A literature review

INTRODUCTION

Ultimate strength-density relationships for bone have been reported with widely varying results. Reliable bone strength predictions are crucial for many applications that aim to assess bone failure. Bone density and bone morphology have been proposed to explain most of the variance in measured bone strength. If this holds true, it could lead to the derivation of a single ultimate strength-density-morphology relationship for all anatomical sites.

METHODS

All relevant literature was reviewed. Ultimate strength-density relationships derived from mechanical testing of human bone tissue were included. The reported relationships were translated to ultimate strength-apparent density relationships and normalized with respect to strain rate. Results were grouped based on bone tissue type (cancellous or cortical), anatomical site, and loading mode (tension vs. compression). When possible, the relationships were compared to existing ultimate strength-density-morphology relationships.

RESULTS

Relationships that considered bone density and morphology covered the full spectrum of eight-fold inter-study difference in reported compressive ultimate strength-density relationships for trabecular bone. This was true for studies that tested specimens in different loading direction and tissue from different anatomical sites. Sparse data was found for ultimate strength-density relationships in tension and for cortical bone properties transverse to the main loading axis of the bone.

CONCLUSIONS

Ultimate strength-density-morphology relationships could explain measured strength across anatomical sites and loading directions. We recommend testing of bone specimens in other directions than along the main trabecular alignment and to include bone morphology in studies that investigate bone material properties. The lack of tensile strength data did not allow for drawing conclusions on ultimate strength-density-morphology relationships. Further studies are needed. Ideally, these studies would investigate both tensile and compressive strength-density relationships, including morphology, to close this gap and lead to more accurate evaluation of bone failure.

Anisotropic computational modelling of bony structures from CT data: An almost automatic procedure

BACKGROUND AND OBJECTIVE

The use of modelling techniques that combine CT data and bone tissue micromechanics is spreading in computational biomechanics. Finite Element models show great potential in surgical planning of intervention and in prediction of stress and strain fields through a non-invasive method. The main challenge pertains to the reliable characterization of bone mechanical behaviour. An almost automatic procedure is here defined, which provides computational models of bony structures considering the actual anisotropy of bone tissue response. The innovative aspect resides on the automatic detection of the directions of anisotropy as the eigenvectors of a three-dimensional distribution matrix of HU values.

METHODS

The procedure combines CT data and micromechanics modelling techniques. Regarding a specific location, the procedure reports both the orthotropic elastic constants, by the analysis of the local HU value, and the anisotropic material directions, by the analysis of the HU values distribution around the specific location.

RESULTS

The procedure returns the distribution of bone tissue orthotropic elasticity tensor. The procedure proves to correctly respect the differentiation between cortical and trabecular bone. Principal directions show to be consistent with experimental data from ultrasound measurements. Regarding the material mapping from voxel to FE model, the developed strategies show to be reliable, leading to marginal errors (lower than 10%) for most of CT voxels (more than 90%). The computational analyses of typical structural loading conditions lead to strain values that are comparable with results from strain gauges experimentations. The development and the exploitation of FE models of different bony structures allow assessing the reliability of the procedure for cortical bone.

CONCLUSIONS

The results highlight the potentialities of the procedure in providing accurate patient-specific biomechanical models of bony structures starting from CT data. The accuracy and the automatism of the procedure are important factors for the development of real time clinical tools. The main limitations of this work remain the not fully automatism and the reliability assessment, which is based mainly on cortical bone regions only.

Effect of orthopedic implants on canine long bone compression stiffness: a combined experimental and computational approach

Osteosynthesis for canine long bones is a complex process requiring knowledge of biology, surgical techniques and (bio)mechanical principles. Subject-specific finite element analysis constitutes a promising tool to evaluate the effect of surgical intervention on the global properties of a bone-implant construct, but suffers from a lack of validation. In this study, the biomechanical behavior of 10 canine humeri was compared before and after creation of a 10 mm bone defect stabilized with an eight-hole locking compression plate (Synthes^®) and two locking screws on each fragment. The response under compression of both intact and plated samples was measured experimentally and reproduced with a finite element model. The experimental stiffness ratio between plated and intact bone was equal to 0.39 ± 0.06. A subject-specific finite element analysis including density-dependent elasto-plastic material properties for canine bone and automatic generation of orthopedic implants was then conducted to recover these experimental results. The stiffness of intact and plated samples could be predicted, with no significant differences with experimental data. The simulated stiffness ratio between plated and intact canine bone was equal to 0.43 ± 0.03. This study constitutes a first step toward the building of a virtual database of pre-computed cases, aiming at helping the veterinary surgeons to make decisions regarding the most suited orthopedic solution for a given dog and a given fracture.

In vitro experimental and numerical study on biomechanics and stability of a novel adjustable hemipelvic prosthesis

Hemipelvic prostheses are used to reconstruct the damaged pelvis due to bone tumors and traumas. However, biomechanical properties of the reconstructed pelvis remain unclear, causing difficulties to implant development and prediction of surgical outcome. In this study, a novel adjustable hemipelvic prosthesis for the Type 1-3 pelvis resection was used to reconstruct the intact pelvic ring. Two types of Pedicle Screw Rod Systems were proposed to improve the stability of fixation between the prosthesis and the bone. Finite Element models of the reconstructed pelvis were built to analyze the performance of the prosthesis and PSRS. Moreover, an in vitro experimental study was performed to measure the deformation of the human reconstructed pelvis. Numerical results agree well with the experimental data. It was found that displacements and stresses bilaterally transferred more evenly in the reconstructed pelvis enhanced by bilateral Pedicle Screw Rod System. The load-transfer function of the pelvis under double-leg standing stance could be recovered. The bilateral pedicle system has better biomechanical performance than the unilateral pedicle system.

Application of quality by design for 3D printed bone prostheses and scaffolds

3D printing is an emergent manufacturing technology recently being applied in the medical field for the development of custom bone prostheses and scaffolds. However, successful industry transformation to this new design and manufacturing approach requires technology integration, concurrent multi-disciplinary collaboration, and a robust quality management framework. This latter change enabler is the focus of this study. While a number of comprehensive quality frameworks have been developed in recent decades to ensure that the manufacturing of medical devices produces reliable products, they are centred on the traditional context of standardised manufacturing techniques. The advent of 3D printing technologies and the prospects for mass customisation provides significant market opportunities, but also presents a serious challenge to regulatory bodies tasked with managing and assuring product quality and safety. Before 3D printing bone prostheses and scaffolds can gain traction, industry stakeholders, such as regulators, clients, medical practitioners, insurers, lawyers, and manufacturers, would all require a high degree of confidence that customised manufacturing can achieve the same quality outcomes as standardised manufacturing. A Quality by Design (QbD) approach to custom 3D printed prostheses can help to ensure that products are designed and manufactured correctly from the beginning without errors. This paper reports on the adaptation of the QbD approach for the development process of 3D printed custom bone prosthesis and scaffolds. This was achieved through the identification of the Critical Quality Attributes of such products, and an extensive review of different design and fabrication methods for 3D printed bone prostheses. Research outcomes include the development of a comprehensive design and fabrication process flow diagram, and categorised risks associated with the design and fabrication processes of such products. An extensive systematic literature review and post-hoc evaluation survey with experts was completed to evaluate the likely effectiveness of the herein suggested QbD framework.

The use of subject-specific Finite Element analysis of L1-L4 vertebra to screening osteoporosis in postmenopausal women

The aim of this work is to study Computed Tomography (CT)-based Finite Element (FE) modelling of lumbar vertebra to classify osteoporotic patients using dual-energy x-ray absorptiometry (DXA) as reference. Cohort comprised 15 postmenopausal female patients with CT-DXA pairs. Sixty subject-specific CT-based FE models were constructed of L1 to L4 vertebrae. Load-displacement in each node of the models was analyzed. Mean kurtosis values of displacement distributions between osteoporosis and healthy groups at L1 level showed a statistically significant difference (p<;0.005). The CT-based FE model of L1 vertebra with patient specific material properties can be used to identify osteoporosis patients.

Biomechanical advantages of robot-assisted pedicle screw fixation in posterior lumbar interbody fusion compared with freehand technique in a prospective randomized controlled trial-perspective for patient-specific finite element analysis

BACKGROUND CONTEXT

There have been conflicting results on the surgical outcome of lumbar fusion surgery using two different techniques: robot-assisted pedicle screw fixation and conventional freehand technique. In addition, there have been no studies about the biomechanical issues between both techniques.

PURPOSE

This study aimed to investigate the biomechanical properties in terms of stress at adjacent segments using robot-assisted pedicle screw insertion technique (robot-assisted, minimally invasive posterior lumbar interbody fusion, Rom-PLIF) and freehand technique (conventional, freehand, open approach, posterior lumbar interbody fusion, Cop-PLIF) for instrumented lumbar fusion surgery.

STUDY DESIGN

This is an additional post-hoc analysis for patient-specific finite element (FE) model.

PATIENT SAMPLE

The sample is composed of patients with degenerative lumbar disease.

OUTCOME MEASURES

Intradiscal pressure and facet contact force are the outcome measures.

METHODS

Patients were randomly assigned to undergo an instrumented PLIF procedure using a Rom-PLIF (37 patients) or a Cop-PLIF (41), respectively. Five patients in each group were selected using a simple random sampling method after operation, and 10 preoperative and postoperative lumbar spines were modeled from preoperative high-resolution computed tomography of 10 patients using the same method for a validated lumbar spine model. Under four pure moments of 7.5 Nm, the changes in intradiscal pressure and facet joint contact force at the proximal adjacent segment following fusion surgery were analyzed and compared with preoperative states.

RESULTS

The representativeness of random samples was verified. Both groups showed significant increases in postoperative intradiscal pressure at the proximal adjacent segment under four moments, compared with the preoperative state. The Cop-PLIF models demonstrated significantly higher percent increments of intradiscal pressure at proximal adjacent segments under extension, lateral bending, and torsion moments than the Rom-PLIF models (p=.032, p=.008, and p=.016, respectively). Furthermore, the percent increment of facet contact force was significantly higher in the Cop-PLIF models under extension and torsion moments than in the Rom-PLIF models (p=.016 under both extension and torsion moments).

CONCLUSIONS

The present study showed the clinical application of subject-specific FE analysis in the spine. Even though there was biomechanical superiority of the robot-assisted insertions in terms of alleviation of stress increments at adjacent segments after fusion, cautious interpretation is needed because of the small sample size.

Development of a biomechanical model of the wrist joint for patient-specific model guided surgical therapy planning: Part 1

An enhanced musculoskeletal biomechanical model of the wrist joint is presented in this article. The developed computational model features the two forearm bones radius and ulna, the eight wrist bones, the five metacarpal bones, and a soft tissue apparatus. Validation of the model was based on information taken from the literature as well as own experimental passive in vitro motion analysis of eight cadaver specimens. The computational model is based on the multi-body simulation software AnyBody. A comprehensive ligamentous apparatus was implemented allowing the investigation of ligament function. The model can easily patient specific personalized on the basis of image information. The model enables simulation of individual wrist motion and predicts trends correctly in the case of changing kinematics. Therefore, patient-specific multi-body simulation models are potentially valuable tools for surgeons in pre- and intraoperative planning of implant placement and orientation.

The Role of Fluid Dynamics in Distributing Ankle Stresses in Anatomic and Injured States

BACKGROUND

In 1976, Ramsey and Hamilton published a landmark cadaveric study demonstrating a dramatic 42% decrease in tibiotalar contact area with only 1 mm of lateral talar shift. An increase in maximum principal stress of at least 72% is predicted based on these findings though the delayed development of arthritis in minimally misaligned ankles does not appear to be commensurate with the results found in dry cadaveric models. We hypothesized that synovial fluid could be a previously unrecognized factor that contributes significantly to stress distribution in the tibiotalar joint in anatomic and injured states.

METHODS

As it is not possible to directly measure contact stresses with and without fluid in a cadaveric model, finite element analysis (FEA) was employed for this study. FEA is a modeling technique used to calculate stresses in complex geometric structures by dividing them into small, simple components called elements. Four test configurations were investigated using a finite element model (FEM): baseline ankle alignment, 1 mm laterally translated talus and fibula, and the previous 2 bone orientations with fluid added. The FEM selected for this study was the Global Human Body Models Consortium-owned GHBMC model, M50 version 4.2, a model of an average-sized male (distributed by Elemance, LLC, Winston-Salem, NC). The ankle was loaded at the proximal tibia with a distributed load equal to the GHBMC body weight, and the maximum principal stress was computed.

RESULTS

All numerical simulations were stable and completed with no errors. In the baseline anatomic configuration, the addition of fluid between the tibia, fibula, and talus reduced the maximum principal stress computed in the distal tibia at maximum load from 31.3 N/mm² to 11.5 N/mm². Following 1 mm lateral translation of the talus and fibula, there was a modest 30% increase in the maximum stress in fluid cases. Qualitatively, translation created less high stress locations on the tibial plafond when fluid was incorporated into the model.

CONCLUSIONS

The findings in this study demonstrate a meaningful role for synovial fluid in distributing stresses within the ankle that has not been considered in historical dry cadaveric studies. The increase in maximum stress predicted by simulation of an ankle with fluid was less than half that projected by cadaveric data, indicating a protective effect of fluid in the injured state. The trends demonstrated by these simulations suggest that bony alignment and fluid in the ankle joint change loading patterns on the tibia and should be accounted for in future experiments.

CLINICAL RELEVANCE

Synovial fluid may play a protective role in ankle injuries, thus delaying the onset of arthritis. Reactive joint effusions may also function to additionally redistribute stresses with higher volumes of viscous fluid.

Quantitative Computed Tomography (QCT) derived Bone Mineral Density (BMD) in finite element studies: a review of the literature

BACKGROUND

Finite element modeling of human bone provides a powerful tool to evaluate a wide variety of outcomes in a highly repeatable and parametric manner. These models are most often derived from computed tomography data, with mechanical properties related to bone mineral density (BMD) from the x-ray energy attenuation provided from this data. To increase accuracy, many researchers report the use of quantitative computed tomography (QCT), in which a calibration phantom is used during image acquisition to improve the estimation of BMD. Since model accuracy is dependent on the methods used in the calculation of BMD and density-mechanical property relationships, it is important to use relationships developed for the same anatomical location and using the same scanner settings, as these may impact model accuracy. The purpose of this literature review is to report the relationships used in the conversion of QCT equivalent density measures to ash, apparent, and/or tissue densities in recent finite element (FE) studies used in common density-modulus relationships. For studies reporting experimental validation, the validation metrics and results are presented.

RESULTS

Of the studies reviewed, 29% reported the use of a dipotassium phosphate (K2HPO4) phantom, 47% a hydroxyapatite (HA) phantom, 13% did not report phantom type, 7% reported use of both K2HPO4 and HA phantoms, and 4% alternate phantom types. Scanner type and/or settings were omitted or partially reported in 31% of studies. The majority of studies used densitometric and/or density-modulus relationships derived from different anatomical locations scanned in different scanners with different scanner settings. The methods used to derive various densitometric relationships are reported and recommendations are provided toward the standardization of reporting metrics.

CONCLUSIONS

This review assessed the current state of QCT-based FE modeling with use of clinical scanners. It was found that previously developed densitometric relationships vary by anatomical location, scanner type and settings. Reporting of all parameters used when referring to previously developed relationships, or in the development of new relationships, may increase the accuracy and repeatability of future FE models.

Biomechanical analysis of a novel hemipelvic endoprosthesis during ascending and descending stairs

In this study, the biomechanical characteristic of a newly developed adjustable hemipelvic prosthesis under dynamic loading conditions was investigated using explicit finite element method. Both intact and reconstructed pelvis models, including pelvis, femur and soft tissues, were established referring to human anatomic data using a solid geometry of a human pelvic bone. Hip contact forces during ascending stairs and descending stairs were imposed on pelvic models. Results showed that maximum von Mises stresses in reconstructed pelvis were 421.85 MPa for prostheses and 109.12 MPa for cortical bone, which were still within a low and elastic range below the yielding strength of Ti-6Al-4V and cortical bone, respectively. Besides, no significant difference of load transferring paths along pelvic rings was observed between the reconstructed pelvis and natural pelvis models. And good agreement was found between the overall distribution of maximum principal stresses in trabecular bones of reconstructed pelvis and natural pelvis, while at limited stances, principal stresses in trabecular bone of reconstructed pelvis were slightly lower than natural pelvis. The results indicated that the load transferring function of pelvis could be restored by this adjustable hemipelvic prosthesis. Moreover, the prosthesis was predicted to have a reliable short- and long-term performance. However, due to the occurrence of slightly lower principal stresses at a few stances, a porous structure applied on the interface between the prosthesis and bone would be studied in future work to obtain better long-term stability.

Automatic deformable surface registration for medical applications by radial basis function-based robust point-matching

Deformable surface mesh registration is a useful technique for various medical applications, such as intra-operative treatment guidance and intra- or inter-patient study. In this paper, we propose an automatic deformable mesh registration technique. The proposed method iteratively deforms a source mesh to a target mesh without manual feature extraction. Each iteration of the registration consists of two steps, automatic correspondence finding using robust point-matching (RPM) and local deformation using a radial basis function (RBF). The proposed RBF-based RPM algorithm solves the interlocking problems of correspondence and deformation using a deterministic annealing framework with fuzzy correspondence and RBF interpolation. Simulation tests showed promising results, with the average deviations decreasing by factors of 21.2 and 11.9, respectively. In the human model test, the average deviation decreased from 1.72±1.88mm to 0.57±0.66mm. We demonstrate the effectiveness of the proposed method by presenting some medical applications.

An Integrated Musculoskeletal-Finite-Element Model to Evaluate Effects of Load Carriage on the Tibia During Walking

Prior studies have assessed the effects of load carriage on the tibia. Here, we expand on these studies and investigate the effects of load carriage on joint reaction forces (JRFs) and the resulting spatiotemporal stress/strain distributions in the tibia. Using full-body motion and ground reaction forces from a female subject, we computed joint and muscle forces during walking for four load carriage conditions. We applied these forces as physiological loading conditions in a finite-element (FE) analysis to compute strain and stress. We derived material properties from computed tomography (CT) images of a sex-, age-, and body mass index-matched subject using a mesh morphing and mapping algorithm, and used them within the FE model. Compared to walking with no load, the knee JRFs were the most sensitive to load carriage, increasing by as much as 26.2% when carrying a 30% of body weight (BW) load (ankle: 16.4% and hip: 19.0%). Moreover, our model revealed disproportionate increases in internal JRFs with increases in load carriage, suggesting a coordinated adjustment in the musculature functions in the lower extremity. FE results reflected the complex effects of spatially varying material properties distribution and muscular engagement on tibial biomechanics during walking. We observed high stresses on the anterior crest and the medial surface of the tibia at pushoff, whereas high cumulative stress during one walking cycle was more prominent in the medioposterior aspect of the tibia. Our findings reinforce the need to include: (1) physiologically accurate loading conditions when modeling healthy subjects undergoing short-term exercise training and (2) the duration of stress exposure when evaluating stress-fracture injury risk. As a fundamental step toward understanding the instantaneous effect of external loading, our study presents a means to assess the relationship between load carriage and bone biomechanics.

Interactive graph-cut segmentation for fast creation of finite element models from clinical ct data for hip fracture prediction

In this study, we propose interactive graph cut image segmentation for fast creation of femur finite element (FE) models from clinical computed tomography scans for hip fracture prediction. Using a sample of N = 48 bone scans representing normal, osteopenic and osteoporotic subjects, the proximal femur was segmented using manual (gold standard) and graph cut segmentation. Segmentations were subsequently used to generate FE models to calculate overall stiffness and peak force in a sideways fall simulations. Results show that, comparable FE results can be obtained with the graph cut method, with a reduction from 20 to 2-5 min interaction time. Average differences between segmentation methods of 0.22 mm were not significantly correlated with differences in FE derived stiffness (R² = 0.08, p = 0.05) and weakly correlated to differences in FE derived peak force (R² = 0.16, p = 0.01). We further found that changes in automatically assigned boundary conditions as a consequence of small segmentation differences were significantly correlated with FE derived results. The proposed interactive graph cut segmentation software MITK-GEM is freely available online at https://simtk.org/home/mitk-gem .

Etiology of Femoroacetabular Impingement in Athletes: A Review of Recent Findings

The relationship between hip deformities and osteoarthritis has recently received a lot of attention. In particular, it has been shown that both osteoarthritis and its precursors, such as the hip deformities that lead to femoroacetabular impingement (FAI), are more prevalent in elite athletes compared with the general population. However, the etiology of the above-mentioned types of hip deformity is not currently well understood. Many recent studies have attempted to shed light on the etiology of this disease. In this article, the main clinical, radiological, mechanobiological, and biomechanical findings of relevance to understanding the etiology of hip deformities leading to FAI are reviewed. Based on these findings, a consistent biomechanical theory explaining the development of hip deformities in athletes is then presented. According to the presented theory, the repetitive, impact-like musculoskeletal loads that athletes experience, particularly when they undertake extreme ranges of hip motion, cause the development of hip deformities. According to this theory, these musculoskeletal loads trigger abnormal growth patterns during the years of skeletal development and cause the formation of hip deformities. A number of hypotheses based on the proposed theory are then formulated that could be tested in future studies to ascertain whether the proposed theory could sufficiently describe the development of hip deformities in athletes.

Patient-specific geometrical modeling of orthopedic structures with high efficiency and accuracy for finite element modeling and 3D printing

We improved the geometrical modeling procedure for fast and accurate reconstruction of orthopedic structures. This procedure consists of medical image segmentation, three-dimensional geometrical reconstruction, and assignment of material properties. The patient-specific orthopedic structures reconstructed by this improved procedure can be used in the virtual surgical planning, 3D printing of real orthopedic structures and finite element analysis. A conventional modeling consists of: image segmentation, geometrical reconstruction, mesh generation, and assignment of material properties. The present study modified the conventional method to enhance software operating procedures. Patient's CT images of different bones were acquired and subsequently reconstructed to give models. The reconstruction procedures were three-dimensional image segmentation, modification of the edge length and quantity of meshes, and the assignment of material properties according to the intensity of gravy value. We compared the performance of our procedures to the conventional procedures modeling in terms of software operating time, success rate and mesh quality. Our proposed framework has the following improvements in the geometrical modeling: (1) processing time: (femur: 87.16 ± 5.90 %; pelvis: 80.16 ± 7.67 %; thoracic vertebra: 17.81 ± 4.36 %; P < 0.05); (2) least volume reduction (femur: 0.26 ± 0.06 %; pelvis: 0.70 ± 0.47, thoracic vertebra: 3.70 ± 1.75 %; P < 0.01) and (3) mesh quality in terms of aspect ratio (femur: 8.00 ± 7.38 %; pelvis: 17.70 ± 9.82 %; thoracic vertebra: 13.93 ± 9.79 %; P < 0.05) and maximum angle (femur: 4.90 ± 5.28 %; pelvis: 17.20 ± 19.29 %; thoracic vertebra: 3.86 ± 3.82 %; P < 0.05). Our proposed patient-specific geometrical modeling requires less operating time and workload, but the orthopedic structures were generated at a higher rate of success as compared with the conventional method. It is expected to benefit the surgical planning of orthopedic structures with less operating time and high accuracy of modeling.

The influence of bone density and anisotropy in finite element models of distal radius fracture osteosynthesis: Evaluations and comparison to experiments

Continuum-level finite element (FE) models can be used to analyze and improve osteosynthesis procedures for distal radius fractures (DRF) from a biomechanical point of view. However, previous models oversimplified the bone material and lacked thorough experimental validation. The goal of this study was to assess the influence of local bone density and anisotropy in FE models of DRF osteosynthesis for predictions of axial stiffness, implant plate stresses, and screw loads. Experiments and FE analysis were conducted in 25 fresh frozen cadaveric radii with DRFs treated by volar locking plate osteosynthesis. Specimen specific geometries were captured using clinical quantitative CT (QCT) scans of the prepared samples. Local bone material properties were computed based on high resolution CT (HR-pQCT) scans of the intact radii. The axial stiffness and individual screw loads were evaluated in FE models, with (1) orthotropic inhomogeneous (OrthoInhom), (2) isotropic inhomogeneous (IsoInhom), and (3) isotropic homogeneous (IsoHom) bone material and compared to the experimental axial stiffness and screw-plate interface failures. FE simulated and experimental axial stiffness correlated significantly (p<0.0001) for all three model types. The coefficient of determination was similar for OrthoInhom (R(2)=0.807) and IsoInhom (R(2)=0.816) models but considerably lower for IsoHom models (R(2)=0.500). The peak screw loads were in qualitative agreement with experimental screw-plate interface failure. Individual loads and implant plate stresses of IsoHom models differed significantly (p<0.05) from OrthoInhom and IsoInhom models. In conclusion, including local bone density in FE models of DRF osteosynthesis is essential whereas local bone anisotropy hardly effects the models׳ predictive abilities.

Computed tomography-based finite element analysis to assess fracture risk and osteoporosis treatment

Finite element analysis (FEA) is a computer technique of structural stress analysis and developed in engineering mechanics. FEA has developed to investigate structural behavior of human bones over the past 40 years. When the faster computers have acquired, better FEA, using 3-dimensional computed tomography (CT) has been developed. This CT-based finite element analysis (CT/FEA) has provided clinicians with useful data. In this review, the mechanism of CT/FEA, validation studies of CT/FEA to evaluate accuracy and reliability in human bones, and clinical application studies to assess fracture risk and effects of osteoporosis medication are overviewed.

ASSESSMENT OF OSTEOPOROTIC FEMORAL FRACTURE RISK: FINITE ELEMENT METHOD AS A POTENTIAL REPLACEMENT FOR CURRENT CLINICAL TECHNIQUES

Femoral fracture risk prediction is a necessary step preceding effective pharmacological intervention or pre-operative planning. Current clinical methods for fracture risk prediction rely on 2D imaging methods and have limited predictive value. Researchers are therefore trying to find improved methods for fracture prediction. During last few decades, many studies have focused on integration of 3D imaging techniques and the finite element (FE) method to improve the accuracy of fracture assessment techniques. In this paper, we review the recent advances in FE and other techniques for predicting the risk of femoral fractures. Based on a number of selected studies, the different steps that are involved in generation of patient-specific FE models are reviewed with particular emphasis on the fracture criteria. The inaccuracies that might arise due to the imperfections of the involved steps are also discussed. It is concluded that compared to image- and geometry-based techniques, FE is a more promising approach for prediction of fracture loads. However, certain technological advancements in FE modeling protocols are required before FE modeling can be recruited in clinical settings.

Additively Manufactured Open-Cell Porous Biomaterials Made from Six Different Space-Filling Unit Cells: The Mechanical and Morphological Properties

It is known that the mechanical properties of bone-mimicking porous biomaterials are a function of the morphological properties of the porous structure, including the configuration and size of the repeating unit cell from which they are made. However, the literature on this topic is limited, primarily because of the challenge in fabricating porous biomaterials with arbitrarily complex morphological designs. In the present work, we studied the relationship between relative density (RD) of porous Ti6Al4V EFI alloy and five compressive properties of the material, namely elastic gradient or modulus (E_s20–70), first maximum stress, plateau stress, yield stress, and energy absorption. Porous structures with different RD and six different unit cell configurations (cubic (C), diamond (D), truncated cube (TC), truncated cuboctahedron (TCO), rhombic dodecahedron (RD), and rhombicuboctahedron (RCO)) were fabricated using selective laser melting. Each of the compressive properties increased with increase in RD, the relationship being of a power law type. Clear trends were seen in the influence of unit cell configuration and porosity on each of the compressive properties. For example, in terms of E_s20–70, the structures may be divided into two groups: those that are stiff (comprising those made using C, TC, TCO, and RCO unit cell) and those that are compliant (comprising those made using D and RD unit cell).

Effects of anodizing parameters and heat treatment on nanotopographical features, bioactivity, and cell culture response of additively manufactured porous titanium

Anodizing could be used for bio-functionalization of the surfaces of titanium alloys. In this study, we use anodizing for creating nanotubes on the surface of porous titanium alloy bone substitutes manufactured using selective laser melting. Different sets of anodizing parameters (voltage: 10 or 20V anodizing time: 30min to 3h) are used for anodizing porous titanium structures that were later heat treated at 500°C. The nanotopographical features are examined using electron microscopy while the bioactivity of anodized surfaces is measured using immersion tests in the simulated body fluid (SBF). Moreover, the effects of anodizing and heat treatment on the performance of one representative anodized porous titanium structures are evaluated using in vitro cell culture assays using human periosteum-derived cells (hPDCs). It has been shown that while anodizing with different anodizing parameters results in very different nanotopographical features, i.e. nanotubes in the range of 20 to 55nm, anodized surfaces have limited apatite-forming ability regardless of the applied anodizing parameters. The results of in vitro cell culture show that both anodizing, and thus generation of regular nanotopographical feature, and heat treatment improve the cell culture response of porous titanium. In particular, cell proliferation measured using metabolic activity and DNA content was improved for anodized and heat treated as well as for anodized but not heat-treated specimens. Heat treatment additionally improved the cell attachment of porous titanium surfaces and upregulated expression of osteogenic markers. Anodized but not heat-treated specimens showed some limited signs of upregulated expression of osteogenic markers. In conclusion, while varying the anodizing parameters creates different nanotube structure, it does not improve apatite-forming ability of porous titanium. However, both anodizing and heat treatment at 500°C improve the cell culture response of porous titanium.

A clinically oriented introduction and review on finite element models of the human cochlea

Due to the inaccessibility of the inner ear, direct in vivo information on cochlear mechanics is difficult to obtain. Mathematical modelling is a promising way to provide insight into the physiology and pathology of the cochlea. Finite element method (FEM) is one of the most popular discrete mathematical modelling techniques, mainly used in engineering that has been increasingly used to model the cochlea and its elements. The aim of this overview is to provide a brief introduction to the use of FEM in modelling and predicting the behavior of the cochlea in normal and pathological conditions. It will focus on methodological issues, modelling assumptions, simulation of clinical scenarios, and pathologies.

Patient-specific bone modeling and analysis: the role of integration and automation in clinical adoption

Patient-specific analysis of bones is considered an important tool for diagnosis and treatment of skeletal diseases and for clinical research aimed at understanding the etiology of skeletal diseases and the effects of different types of treatment on their progress. In this article, we discuss how integration of several important components enables accurate and cost-effective patient-specific bone analysis, focusing primarily on patient-specific finite element (FE) modeling of bones. First, the different components are briefly reviewed. Then, two important aspects of patient-specific FE modeling, namely integration of modeling components and automation of modeling approaches, are discussed. We conclude with a section on validation of patient-specific modeling results, possible applications of patient-specific modeling procedures, current limitations of the modeling approaches, and possible areas for future research.