Subject specific finite element mesh generation of the pelvis from biplanar x-ray images: application to 120 clinical cases

Several Finite Element (FE) models of the pelvis have been developed to comprehensively assess the onset of pathologies and for clinical and industrial applications. However, because of the difficulties associated with the creation of subject-specific FE mesh from CT scan and MR images, most of the existing models rely on the data of one given individual. Moreover, although several fast and robust methods have been developed for automatically generating tetrahedral meshes of arbitrary geometries, hexahedral meshes are still preferred today because of their distinct advantages but their generation remains an open challenge. Recently, approaches have been proposed for fast 3D reconstruction of bones based on X-ray imaging. In this study, we adapted such an approach for the fast and automatic generation of all-hexahedral subject-specific FE models of the pelvis based on the elastic registration of a generic mesh to the subject-specific target in conjunction with element regularity and quality correction. The technique was successfully tested on a database of 120 3D reconstructions of pelvises from biplanar X-ray images. For each patient, a full hexahedral subject-specific FE mesh was generated with an accurate surface representation.     

Mesh considerations for finite element blast modelling in biomechanics

Finite element (FE) modelling is a popular tool for studying human body response to blast exposure. However, blast modelling is a complex problem owing to more numerous fluid–structure interactions (FSIs) and the high–frequency loading that accompanies blast exposures. This study investigates FE mesh design for blast modelling using a sphere in a closed-ended shock tube meshed with varying element sizes using both tetrahedral and hexahedral elements. FSI was consistent for sphere-to-fluid element ratios between 0.25 and 4, and acceleration response was similar for both element types (R² = 0.997). Tetrahedral elements were found to become increasingly volatile following shock loading, causing higher pressures and stresses than predicted with the hexahedral elements. Deviatoric stress response was dependent on the sphere mesh size (p < 0.001), while the pressure response was dependent on the shock tube mesh size (p < 0.001). The results of this study highlight the necessity for mesh sensitivity analysis in blast models.     

Biomechanical measurements on scaphoid bone screws in an experimental model

A number of screws commonly used for internal fixation in scaphoid bone fractures and nonunions are compared regarding biomechanical properties and clinical applicability. The experiments were carried out on models made of ash-wood, representing a reconstruction and fixation as is performed in a cortico-cancellous inlay bone graft for scaphoid non-union. For fixation use was made of 2.7 and 3.5 AO/ASIF cortical screws respectively, 4.0 AO/ASIF cancellous screws, Herbert screws, and a newly designed screw called the three components screw (D.K.S.). The models with implanted screws were tested for bending strength, tensile strength and torsion stability. No large differences between the various screws were found regarding the measured parameters, so that a small intra-osteal implant such as the Herbert screw and the D.K.S., which can be inserted easily and which gives a certain amount of interfragmentary compression, will be sufficient for osteosynthesis of the scaphoid bone. In case an intra-osteal implant is not available a single 3.5 AO/ASIF cortical screw, inserted following lag-screw principles, is recommended.     

An approximate model for cancellous bone screw fixation

This paper presents a finite element (FE) model to identify parameters that affect the performance of an improved cancellous bone screw fixation technique, and hence potentially improve fracture treatment. In cancellous bone of low apparent density, it can be difficult to achieve adequate screw fixation and hence provide stable fracture fixation that enables bone healing. Data from predictive FE models indicate that cements can have a significant potential to improve screw holding power in cancellous bone. These FE models are used to demonstrate the key parameters that determine pull-out strength in a variety of screw, bone and cement set-ups, and to compare the effectiveness of different configurations. The paper concludes that significant advantages, up to an order of magnitude, in screw pull-out strength in cancellous bone might be gained by the appropriate use of a currently approved calcium phosphate cement.     

Use of a statistical model of the whole femur in a large scale,multi-model study of femoral neck fracture risk

Interpatient variability is often overlooked in orthopaedic computational studies due to the substantial challenges involved in sourcing and generating large numbers of bone models. A statistical model of the whole femur incorporating both geometric and material property variation was developed as a potential solution to this problem. The statistical model was constructed using principal component analysis, applied to 21 individual computer tomography scans. To test the ability of the statistical model to generate realistic, unique, finite element (FE) femur models it was used as a source of 1000 femurs to drive a study on femoral neck fracture risk. The study simulated the impact of an oblique fall to the side, a scenario known to account for a large proportion of hip fractures in the elderly and have a lower fracture load than alternative loading approaches. FE model generation, application of subject specific loading and boundary conditions, FE processing and post processing of the solutions were completed automatically. The generated models were within the bounds of the training data used to create the statistical model with a high mesh quality, able to be used directly by the FE solver without remeshing. The results indicated that 28 of the 1000 femurs were at highest risk of fracture. Closer analysis revealed the percentage of cortical bone in the proximal femur to be a crucial differentiator between the failed and non-failed groups. The likely fracture location was indicated to be intertrochantic. Comparison to previous computational, clinical and experimental work revealed support for these findings.     

Stress analysis of a centrally fractured rib fixated by an intramedullary screw

The stress on an intramedullary screw rib fixation device holding together a centrally fractured human rib under in vivo force loadings was studied using finite element analysis (FEA). Validation of the FEA modelling using pullout from porcine ribs proved FEA to be suitable for assessing the structural integrity of screw/bone systems such as rib fixated by a screw. In the human rib fixation investigation, it was found that intramedullary bioresorbable Bioretec screws can fixate centrally fractured human ribs under normal breathing conditions. However, under coughing conditions, simulation showed Bioretec fixating screws to bend substantially. High stresses in the screw are mainly the result of flexion induced by the force loading, and are restricted to thin regions on the outside of the screw shaft. Stiffer screws result in less locally intense stress concentrations in bone, indicating that bone failure in the bone/screw contact regions can be averted with improvements in screw stiffness.     

An efficient and accurate prediction of the stability of percutaneous fixation of acetabular fractures with finite element simulation

Posterior wall fracture is one of the most common fracture types of the acetabulum and a conventional approach is to perform open reduction and internal fixation with a plate and screws. Percutaneous screw fixations, on the other hand, have recently gained attention due to their benefits such as less exposure and minimization of blood loss. However their biomechanical stability, especially in terms interfragmentary movement, has not been investigated thoroughly. The aims of this study are twofold: (1) to measure the interfragmentary movements in the conventional open approach with plate fixations and the percutaneous screw fixations in the acetabular fractures and compare them; and (2) to develop and validate a fast and efficient way of predicting the interfragmentary movement in percutaneous fixation of posterior wall fractures of the acetabulum using a 3D finite element (FE) model of the pelvis. Our results indicate that in single fragment fractures of the posterior wall of the acetabulum, plate fixations give superior stability to screw fixations. However screw fixations also give reasonable stability as the average gap between fragment and the bone remained less than 1 mm when the maximum load was applied. Our finite element model predicted the stability of screw fixation with good accuracy. Moreover, when the screw positions were optimized, the stability predicted by our FE model was comparable to the stability obtained by plate fixations. Our study has shown that FE modeling can be useful in examining biomechanical stability of osteosynthesis and can potentially be used in surgical planning of osteosynthesis.     

LISS钢板治疗股骨远端骨折有限元建模及分析

目的建立LISS-DF治疗股骨远端骨折近端螺钉不同单双皮质固定的三维有限元模型,并进行初步生物力学分析。方法提取CT图片相关数据,利用自行编写程序生成命令流文件,建立完整股骨以及16个不同LISS-DF治疗股骨远端AO分型33-A3型骨折的实体模型（钢板和股骨不接触、螺钉分别固定于钢板和股骨）,进行网格划分。分析不同载荷作用下完整股骨和LISS钢板近端螺钉全双皮质固定治疗骨折的模型受力状况。结果建立了相关的有限元模型。不同载荷作用下,LISS钢板近端螺钉全双皮质固定模型和完整股骨的应力集中均位于股骨颈内侧和股骨干外侧中下1/3处。相同载荷作用下,LISS钢板近端螺钉全双皮质固定模型的股骨颈部最大等效应力值略减小,股骨干最大等效应力值明显减小。结论研究建立的三维有限元模型,为应用LISS治疗股骨骨折的生物力学分析提供了良好的实验平台和基础。从生物力学角度而言,LISS-DF近端螺钉全双皮质固定为治疗股骨远端复杂骨折的有效方法。     

Micro-CT and micro-FE analysis of pedicle screw fixation under different loading conditions

Anchorage of pedicle screw instrumentation in the elderly spine with poor bone quality remains challenging. In this study, micro finite element (µFE) models were used to assess the specific influence of screw design and the relative contribution of local bone density to fixation mechanics. These were created from micro computer tomography (µCT) scans of vertebras implanted with two types of pedicle screws, including a full region-or-interest of 10 mm radius around each screw, as well as submodels for the pedicle and inner trabecular bone of the vertebral body. The local bone volume fraction (BV/TV) calculated from the µCT scans around different regions of the screw (pedicle, inner trabecular region of the vertebral body) were then related to the predicted stiffness in simulated pull-out tests as well as to the experimental pull-out and torsional fixation properties mechanically measured on the corresponding specimens. Results show that predicted stiffness correlated excellently with experimental pull-out strength (R² > 0.92, p < .043), better than regional BV/TV alone (R2 = 0.79, p = .003). They also show that correlations between fixation properties and BV/TV were increased when accounting only for the pedicle zone (R² = 0.66–0.94, p ≤ .032), but with weaker correlations for torsional loads (R² < 0.10). Our analyses highlight the role of local density in the pedicle zone on the fixation stiffness and strength of pedicle screws when pull-out loads are involved, but that local apparent bone density alone may not be sufficient to explain resistance in torsion.     

Morphing the feature-based multi-blocks of normative/healthy vertebral geometries to scoliosis vertebral geometries: development of personalized finite element models

Personalized Finite Element (FE) models and hexahedral elements are preferred for biomechanical investigations. Feature-based multi-block methods are used to develop anatomically accurate personalized FE models with hexahedral mesh. It is tedious to manually construct multi-blocks for large number of geometries on an individual basis to develop personalized FE models. Mesh-morphing method mitigates the aforementioned tediousness in meshing personalized geometries every time, but leads to element warping and loss of geometrical data. Such issues increase in magnitude when normative spine FE model is morphed to scoliosis-affected spinal geometry. The only way to bypass the issue of hex-mesh distortion or loss of geometry as a result of morphing is to rely on manually constructing the multi-blocks for scoliosis-affected spine geometry of each individual, which is time intensive. A method to semi-automate the construction of multi-blocks on the geometry of scoliosis vertebrae from the existing multi-blocks of normative vertebrae is demonstrated in this paper. High-quality hexahedral elements were generated on the scoliosis vertebrae from the morphed multi-blocks of normative vertebrae. Time taken was 3 months to construct the multi-blocks for normative spine and less than a day for scoliosis. Efforts taken to construct multi-blocks on personalized scoliosis spinal geometries are significantly reduced by morphing existing multi-blocks.     

Growing multiblock structures: a semi-automated approach to block placement for multiblock hexahedral meshing

Finite element (FE) analysis is a cornerstone of orthopaedic biomechanics research. Three-dimensional medical imaging provides sufficient resolution for the subject-specific FE models to be generated from these data-sets. FE model development requires discretisation of a three-dimensional domain, which can be the most time-consuming component of a FE study. Hexahedral meshing tools based on the multiblock method currently rely on the manual placement of building blocks for mesh generation. We hypothesise that angular analysis of the geometric centreline for a three-dimensional surface could be used to automatically generate building block structures for the multiblock hexahedral mesh generation. Our algorithm uses a set of user-defined points and parameters to automatically generate a multiblock structure based on a surface's geometric centreline. This significantly reduces the time required for model development. We have applied this algorithm to 47 bones of varying geometries and successfully generated a FE mesh in all cases. This work represents significant advancement in automatically generating multiblock structures for a wide range of geometries.     

Growing multiblock structures: a semi-automated approach to block placement for multiblock hexahedral meshing

Finite element (FE) analysis is a cornerstone of orthopaedic biomechanics research. Three-dimensional medical imaging provides sufficient resolution for the subject-specific FE models to be generated from these data-sets. FE model development requires discretisation of a three-dimensional domain, which can be the most time-consuming component of a FE study. Hexahedral meshing tools based on the multiblock method currently rely on the manual placement of building blocks for mesh generation. We hypothesise that angular analysis of the geometric centreline for a three-dimensional surface could be used to automatically generate building block structures for the multiblock hexahedral mesh generation. Our algorithm uses a set of user-defined points and parameters to automatically generate a multiblock structure based on a surface's geometric centreline. This significantly reduces the time required for model development. We have applied this algorithm to 47 bones of varying geometries and successfully generated a FE mesh in all cases. This work represents significant advancement in automatically generating multiblock structures for a wide range of geometries.     

Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

Trabecular deformations during screw pull-out: a micro-CT study of lapine bone

The mechanical fixation of endosseous implants, such as screws, in trabecular bone is challenging because of the complex porous microstructure. Development of new screw designs to improve fracture fixation, especially in high-porosity osteoporotic bone, requires a profound understanding of how the structural system implant/trabeculae interacts when it is subjected to mechanical load. In this study, pull-out tests of screw implants were performed. Screws were first inserted into the trabecular bone of rabbit femurs and then pulled out from the bone inside a computational tomography scanner. The tests were interrupted at certain load steps to acquire 3D images. The images were then analysed with a digital volume correlation technique to estimate deformation and strain fields inside the bone during the tests. The results indicate that the highest shear strains are concentrated between the inner and outer thread diameter, whereas compressive strains are found at larger distances from the screw. Tensile strains were somewhat smaller. Strain concentrations and the location of trabecular failures provide experimental information that could be used in the development of new screw designs and/or to validate numerical simulations.     

Cortical screw pullout strength and effective shear stress in synthetic third generation composite femurs

BACKGROUND: The use of artificial bone analogs in biomechanical testing of orthopaedic fracture fixation devices has increased, particularly due to the recent development of commercially available femurs such as the third generation composite femur that closely reproduce the bulk mechanical behavior of human cadaveric and/or fresh whole bone. The purpose of this investigation was to measure bone screw pullout forces in composite femurs and determine whether results are comparable to cadaver data from previous literature. METHOD OF APPROACH: The pullout strengths of 3.5 and 4.5 mm standard bicortical screws inserted into synthetic third generation composite femurs were measured and compared to existing adult human cadaveric and animal data from the literature. RESULTS: For 3.5 mm screws, the measured extraction shear stress in synthetic femurs (23.70-33.99 MPa) was in the range of adult human femurs and tibias (24.4-38.8 MPa). For 4.5 mm screws, the measured values in synthetic femurs (26.04-34.76 MPa) were also similar to adult human specimens (15.9-38.9 MPa). Synthetic femur results for extraction stress showed no statistically significant site-to-site effect for 3.5 and 4.5 mm screws, with one exception. Overall, the 4.5 mm screws showed statistically higher stress required for extraction than 3.5 mm screws. CONCLUSIONS: The third generation composite femurs provide a satisfactory biomechanical analog to human long-bones at the screw-bone interface. However, it is not known whether these femurs perform similarly to human bone during physiological screw "toggling."     

Effect of pedicle screw diameter on screw fixation efficacy in human osteoporotic thoracic vertebrae

The selection of an ideal screw size plays a crucial role in the success of spinal instrumentation as larger diameter screws are thought to provide better fixation strength but increase the risk of pedicle failure during insertion. On the other hand, smaller diameter screws are with lesser risk of pedicle breakage but are thought to compromise the stability of the instrumentation. By investigating the relationship between screw diameter and the pullout strength of pedicle screws after fatigue loading, this study seeks to find quantitative biomechanical data for surgeons in determining the most ideal diameter size screws when performing surgical implementations on osteoporotic vertebrae.Twenty-seven osteoporotic (BMD ranged: 0.353–0.848 g/cm²) thoracic vertebrae (T3-T8) were harvested from 5 human cadavers. Two sizes of poly-axial screws (5.0 mm × 35 and 4.35 mm × 35) were implanted into each pedicles of the vertebrae by an experienced surgeon. Specimens were randomly distributed into control group, fatigue group of 5000 and 10,000 cycles with peak-to-peak loadings of 10–100 N at 1 Hz. Each specimen was then axial pullout tested at a constant rate of 5 mm/min. The ultimate pullout strength (N) & stiffness (N/mm) were obtained for analysis.The results showed that although the larger diameter screws achieved superior pullout strength immediately after the implantation, both sizes of screws exhibited comparable pullout strengths post fatigue loading. This indicates that the smaller diameter screws may be considered for surgical techniques performed on osteoporotic vertebrae for reduced risk of pedicle breakage without sacrificing fixation strength.     

Statistical shape and appearance models for fast and automated estimation of proximal femur fracture load using 2D finite element models

Estimating the risk of osteoporotic fractures is an important diagnostic step that needs to be taken before medicinal treatment. Densitometry-based criteria are normally used in clinical practice for this purpose. However, densitometry-based techniques could not explain all low-energy fractures. As patient-specific finite element (FE) models allow for consideration of other parameters (e.g. load conditions) that are known to be associated with fracture, they are considered promising candidates for more accurate fracture risk estimation. Nevertheless, they are often time consuming, expensive, and complex to build and may need the type of expertise that is not normally available in clinical settings. In this study, we report the development of an automated platform for estimating proximal femur fracture loads using patient-specific 2D FE models generated using dual-energy x-ray absorptiometry (DXA) scans. First, a statistical shape and appearance model (SSAM) is built using DXA scans of patients screened for osteoporosis following a low energy fracture. SSAM is then used together with Active Appearance Models (AAM) for automated segmentation of the proximal femur from new unseen DXA scans. The mean point-to-curve error of the automated procedure, i.e. 1.2–1.4 mm, is shown to be only slightly larger than the intra-observer variability of manual segmentation, i.e. 1.0 mm. Moreover, the developed platform automatically meshes the segmented shape, assigns density-based mechanical properties, assigns loads and boundary conditions, submits the 2D FE model for solution, and performs post-processing of the 2D FE simulation data to determine fracture loads. The fracture loads predicted using the manually generated and automatically generated 2D FE models are shown to be very close with a mean difference of around 8.8%. Repeated measures ANOVA showed no significant differences between the fracture loads calculated using FE models manually generated by three independent observers and those calculated using the automatically generated FE models (p>0.05).     

Enhanced bone screw fixation with biodegradable bone cement in osteoporotic bone model

Purpose: The purpose of this study was to study the potential of novel biodegradable PCL bone cement to improve bone screw fixation strength in osteoporotic bone. Methods: The biomechanical properties of bone cement (ε-polycaprolactone, PCL) and fixation strength were studied using biomechanical tests and bone screws fixed in an osteoporotic bone model. Removal torques and pullout strengths were assessed for cortical, self-tapping, and cancellous screws inserted in the osteoporotic bone model (polyurethane foam blocks with polycarbonate plate) with and without PCL bone cement. Open cell and cellular rigid foam blocks with a density of 0.12 g/cm3 were used in this model. Results: Removal torques were significantly (more than six-fold) improved with bone cement for cancellous screws. Furthermore, the bone cement improved pullout strengths three to 12 times over depending on the screw and model material.?Conclusions: Biodegradable bone cement turned out to be a very potential material to stabilize screw fixation in osteoporotic bone. The results warrant further research before safe clinical use, especially to clarify clinically relevant factors using real osteoporotic bone under human body conditions and dynamic fatigue testing for long-term performance.     

A novel methodology for generating 3D finite element models of the hip from 2D radiographs

Finite element (FE) modelling has been proposed as a tool for estimating fracture risk and patient-specific FE models are commonly based on computed tomography (CT). Here, we present a novel method to automatically create personalised 3D models from standard 2D hip radiographs. A set of geometrical parameters of the femur were determined from seven a–p hip radiographs and compared to the 3D femoral shape obtained from CT as training material; the error in reconstructing the 3D model from the 2D radiographs was assessed. Using the geometry parameters as the input, the 3D shape of another 21 femora was built and meshed, separating a cortical and trabecular compartment. The material properties were derived from the homogeneity index assessed by texture analysis of the radiographs, with focus on the principal tensile and compressive trabecular systems. The ability of these FE models to predict failure load as determined by experimental biomechanical testing was evaluated and compared to the predictive ability of DXA. The average reconstruction error of the 3D models was 1.77 mm (±1.17 mm), with the error being smallest in the femoral head and neck, and greatest in the trochanter. The correlation of the FE predicted failure load with the experimental failure load was r²=64% for the reconstruction FE model, which was significantly better (p<0.05) than that for DXA (r²=24%). This novel method for automatically constructing a patient-specific 3D finite element model from standard 2D radiographs shows encouraging results in estimating patient-specific failure loads.