Development of a statistical shape-function model of the implanted knee for real-time prediction of joint mechanics

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.

Instantaneous Generation of Subject-Specific Finite Element Models of the Hip Capsule

Subject-specific hip capsule models could offer insights into impingement and dislocation risk when coupled with computer-aided surgery, but model calibration is time-consuming using traditional techniques. This study developed a framework for instantaneously generating subject-specific finite element (FE) capsule representations from regression models trained with a probabilistic approach. A validated FE model of the implanted hip capsule was evaluated probabilistically to generate a training dataset relating capsule geometry and material properties to hip laxity. Multivariate regression models were trained using 90% of trials to predict capsule properties based on hip laxity and attachment site information. The regression models were validated using the remaining 10% of the training set by comparing differences in hip laxity between the original trials and the regression-derived capsules. Root mean square errors (RMSEs) in laxity predictions ranged from 1.8° to 2.3°, depending on the type of laxity used in the training set. The RMSE, when predicting the laxity measured from five cadaveric specimens with total hip arthroplasty, was 4.5°. Model generation time was reduced from days to milliseconds. The results demonstrated the potential of regression-based training to instantaneously generate subject-specific FE models and have implications for integrating subject-specific capsule models into surgical planning software.

Robust automatic hexahedral cartilage meshing framework enables population-based computational studies of the knee

Osteoarthritis of the knee is increasingly prevalent as our population ages, representing an increasing financial burden, and severely impacting quality of life. The invasiveness of in vivo procedures and the high cost of cadaveric studies has left computational tools uniquely suited to study knee biomechanics. Developments in deep learning have great potential for efficiently generating large-scale datasets to enable researchers to perform population-sized investigations, but the time and effort associated with producing robust hexahedral meshes has been a limiting factor in expanding finite element studies to encompass a population. Here we developed a fully automated pipeline capable of taking magnetic resonance knee images and producing a working finite element simulation. We trained an encoder-decoder convolutional neural network to perform semantic image segmentation on the Imorphics dataset provided through the Osteoarthritis Initiative. The Imorphics dataset contained 176 image sequences with varying levels of cartilage degradation. Starting from an open-source swept-extrusion meshing algorithm, we further developed this algorithm until it could produce high quality meshes for every sequence and we applied a template-mapping procedure to automatically place soft-tissue attachment points. The meshing algorithm produced simulation-ready meshes for all 176 sequences, regardless of the use of provided (manually reconstructed) or predicted (automatically generated) segmentation labels. The average time to mesh all bones and cartilage tissues was less than 2 min per knee on an AMD Ryzen 5600X processor, using a parallel pool of three workers for bone meshing, followed by a pool of four workers meshing the four cartilage tissues. Of the 176 sequences with provided segmentation labels, 86% of the resulting meshes completed a simulated flexion-extension activity. We used a reserved testing dataset of 28 sequences unseen during network training to produce simulations derived from predicted labels. We compared tibiofemoral contact mechanics between manual and automated reconstructions for the 24 pairs of successful finite element simulations from this set, resulting in mean root-mean-squared differences under 20% of their respective min-max norms. In combination with further advancements in deep learning, this framework represents a feasible pipeline to produce population sized finite element studies of the natural knee from subject-specific models.

Development of an additively-manufactured functionally-graded expandable implant via the application of the adaptive response surface method: feasibility study on intramedullary humerus nail

The paper introduces the concept of the rational design of a deployable humeral intramedullary nail plug based on a honeycomb structure used for the surgical treatment of humeral shaft fractures. The concept serves for to restore the axial alignment of bone fragments and to maintain stability via bone-nail friction and locking screws. The design nail plug was gained by optimisation process the Latin Hypercube Sampling Design algorithm and Multi-Objective Genetic Algorithm. It was shown that we can use statistical shape function combined by 3 D printing for designing of a new rationally designed implants.

The impact of shape uncertainty on aortic-valve pressure-drop computations

Patient-specific image-based computational fluid dynamics (CFD) is widely adopted in the cardiovascular research community to study hemodynamics, and will become increasingly important for personalized medicine. However, segmentation of the flow domain is not exact and geometric uncertainty can be expected which propagates through the computational model, leading to uncertainty in model output. Seventy-four aortic-valves were segmented from computed tomography images at peak systole. Statistical shape modeling was used to obtain an approximate parameterization of the original segmentations. This parameterization was used to train a meta-model that related the first five shape mode coefficients and flowrate to the CFD-computed transvalvular pressure-drop. Consequently, shape uncertainty in the order of 0.5 and 1.0 mm was emulated by introducing uncertainty in the shape mode coefficients. A global variance-based sensitivity analysis was performed to quantify output uncertainty and to determine relative importance of the shape modes. The first shape mode captured the opening/closing behavior of the valve and uncertainty in this mode coefficient accounted for more than 90% of the output variance. However, sensitivity to shape uncertainty is patient-specific, and the relative importance of the fourth shape mode coefficient tended to increase with increases in valvular area. These results show that geometric uncertainty in the order of image voxel size may lead to substantial uncertainty in CFD-computed transvalvular pressure-drops. Moreover, this illustrates that it is essential to assess the impact of geometric uncertainty on model output, and that this should be thoroughly quantified for applications that wish to use image-based CFD models.

Combining statistical shape modeling, CFD, and meta-modeling to approximate the patient-specific pressure-drop across the aortic valve in real-time

BACKGROUND

Advances in medical imaging, segmentation techniques, and high performance computing have stimulated the use of complex, patient-specific, three-dimensional Computational Fluid Dynamics (CFD) simulations. Patient-specific, CFD-compatible geometries of the aortic valve are readily obtained. CFD can then be used to obtain the patient-specific pressure-flow relationship of the aortic valve. However, such CFD simulations are computationally expensive, and real-time alternatives are desired.

AIM

The aim of this work is to evaluate the performance of a meta-model with respect to high-fidelity, three-dimensional CFD simulations of the aortic valve.

METHODS

Principal component analysis was used to build a statistical shape model (SSM) from a population of 74 iso-topological meshes of the aortic valve. Synthetic meshes were created with the SSM, and steady-state CFD simulations at flow-rates between 50 and 650 mL/s were performed to build a meta-model. The meta-model related the statistical shape variance, and flow-rate to the pressure-drop.

RESULTS

Even though the first three shape modes account for only 46% of shape variance, the features relevant for the pressure-drop seem to be captured. The three-mode shape-model approximates the pressure-drop with an average error of 8.8% to 10.6% for aortic valves with a geometric orifice area below 150 mm² . The proposed methodology was least accurate for aortic valve areas above 150 mm² . Further reduction to a meta-model introduces an additional 3% error.

CONCLUSIONS

Statistical shape modeling can be used to capture shape variation of the aortic valve. Meta-models trained by SSM-based CFD simulations can provide an estimate of the pressure-flow relationship in real-time.