Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

Comparison of theoretical fixation stability of three devices employed in medial opening wedge high tibial osteotomy: a finite element analysis.

BACKGROUND Medial open wedge high tibial osteotomy is a well-established procedure for the treatment of unicompartmental osteoarthritis and symptomatic varus malalignment. We hypothesized that different fixation devices generate different fixation stability profiles for the various wedge sizes in a finite element (FE) analysis. METHODS Four types of fixation were compared: 1) first and 2) second generation Puddu plates, and 3) TomoFix plate with and 4) without bone graft. Cortical and cancellous bone was modelled and five different opening wedge sizes were studied for each model. Outcome measures included: 1) stresses in bone, 2) relative displacement of the proximal and distal tibial fragments, 3) stresses in the plates, 4) stresses on the upper and lower screw surfaces in the screw channels. RESULTS The highest load for all fixation types occurred in the plate axis. For the vast majority of the wedge sizes and fixation types the shear stress (von Mises stress) was dominating in the bone independent of fixation type. The relative displacements of the tibial fragments were low (in μm range). With an increasing wedge size this displacement tended to increase for both Puddu plates and the TomoFix plate with bone graft. For the TomoFix plate without bone graft a rather opposite trend was observed.For all fixation types the occurring stresses at the screw-bone contact areas pulled at the screws and exceeded the allowable threshold of 1.2 MPa for at least one screw surface. Of the six screw surfaces that were studied, the TomoFix plate with bone graft showed a stress excess of one out of twelve and without bone graft, five out of twelve. With the Puddu plates, an excess stress occurred in the majority of screw surfaces. CONCLUSIONS The different fixation devices generate different fixation stability profiles for different opening wedge sizes. Based on the computational simulations, none of the studied osteosynthesis fixation types warranted an intransigent full weight bearing per se. The highest fixation stability was observed for the TomoFix plates and the lowest for the first generation Puddu plate. These findings were revealed in theoretical models and need to be validated in controlled clinical settings.

Finite element analysis in asymptomatic, symptomatic, and ruptured abdominal aortic aneurysms: in search of new rupture risk predictors.

OBJECTIVES To compare biomechanical rupture risk parameters of asymptomatic, symptomatic and ruptured abdominal aortic aneurysms (AAA) using finite element analysis (FEA). STUDY DESIGN Retrospective biomechanical single center analysis of asymptomatic, symptomatic, and ruptured AAAs. Comparison of biomechanical parameters from FEA. MATERIALS AND METHODS From 2011 to 2013 computed tomography angiography (CTA) data from 30 asymptomatic, 15 symptomatic, and 15 ruptured AAAs were collected consecutively. FEA was performed according to the successive steps of AAA vessel reconstruction, segmentation and finite element computation. Biomechanical parameters Peak Wall Rupture Risk Index (PWRI), Peak Wall Stress (PWS), and Rupture Risk Equivalent Diameter (RRED) were compared among the three subgroups. RESULTS PWRI differentiated between asymptomatic and symptomatic AAAs (p < .0004) better than PWS (p < .1453). PWRI-dependent RRED was higher in the symptomatic subgroup compared with the asymptomatic subgroup (p < .0004). Maximum AAA external diameters were comparable between the two groups (p < .1355). Ruptured AAAs showed the highest values for external diameter, total intraluminal thrombus volume, PWS, RRED, and PWRI compared with asymptomatic and symptomatic AAAs. In contrast with symptomatic and ruptured AAAs, none of the asymptomatic patients had a PWRI value >1.0. This threshold value might identify patients at imminent risk of rupture. CONCLUSIONS From different FEA derived parameters, PWRI distinguishes most precisely between asymptomatic and symptomatic AAAs. If elevated, this value may represent a negative prognostic factor for asymptomatic AAAs.

Cardiopulmonary Bypass has No Significant Impact on Survival in Patients Undergoing Nephrectomy and Level III-IV Inferior Vena Cava Thrombectomy. A Multi-Institutional Analysis.

PURPOSE The impact of cardiopulmonary bypass in level III-IV tumor thrombectomy on surgical and oncologic outcomes is unknown. We determine the impact of cardiopulmonary bypass on overall and cancer specific survival, as well as surgical complication rates and immediate outcomes in patients undergoing nephrectomy and level III-IV tumor thrombectomy with or without cardiopulmonary bypass. MATERIALS AND METHODS We retrospectively analyzed 362 patients with renal cell cancer and with level III or IV tumor thrombus from 1992 to 2012 at 22 U.S. and European centers. Cox proportional hazards models were used to compare overall and cancer specific survival between patients with and without cardiopulmonary bypass. Perioperative mortality and complication rates were assessed using logistic regression analyses. RESULTS Median overall survival was 24.6 months in noncardiopulmonary bypass cases and 26.6 months in cardiopulmonary bypass cases. Overall survival and cancer specific survival did not differ significantly in both groups on univariate analysis or when adjusting for known risk factors. On multivariate analysis no significant differences were seen in hospital length of stay, Clavien 1-4 complication rate, intraoperative or 30-day mortality and cancer specific survival. Limitations include the retrospective nature of the study. CONCLUSIONS In our multi-institutional analysis the use of cardiopulmonary bypass did not significantly impact cancer specific survival or overall survival in patients undergoing nephrectomy and level III or IV tumor thrombectomy. Neither approach was independently associated with increased mortality on multivariate analysis. Greater surgical complications were not independently associated with the use of cardiopulmonary bypass.

Finite element analysis predicts experimental failure patterns in vertebral bodies loaded via intervertebral discs up to large deformation

Vertebral compression fracture is a common medical problem in osteoporotic individuals. The quantitative computed tomography (QCT)-based finite element (FE) method may be used to predict vertebral strength in vivo, but needs to be validated with experimental tests. The aim of this study was to validate a nonlinear anatomy specific QCT-based FE model by using a novel testing setup. Thirty-seven human thoracolumbar vertebral bone slices were prepared by removing cortical endplates and posterior elements. The slices were scanned with QCT and the volumetric bone mineral density (vBMD) was computed with the standard clinical approach. A novel experimental setup was designed to induce a realistic failure in the vertebral slices in vitro. Rotation of the loading plate was allowed by means of a ball joint. To minimize device compliance, the specimen deformation was measured directly on the loading plate with three sensors. A nonlinear FE model was generated from the calibrated QCT images and computed vertebral stiffness and strength were compared to those measured during the experiments. In agreement with clinical observations, most of the vertebrae underwent an anterior wedge-shape fracture. As expected, the FE method predicted both stiffness and strength better than vBMD (R2 improved from 0.27 to 0.49 and from 0.34 to 0.79, respectively). Despite the lack of fitting parameters, the linear regression of the FE prediction for strength was close to the 1:1 relation (slope and intercept close to one (0.86 kN) and to zero (0.72 kN), respectively). In conclusion, a nonlinear FE model was successfully validated through a novel experimental technique for generating wedge-shape fractures in human thoracolumbar vertebrae.

Clinical Use of Quantitative Computed Tomography–Based Finite Element Analysis of the Hip and Spine in the Management of Osteoporosis in Adults: the 2015 ISCD Official Positions—Part II

The International Society for Clinical Densitometry (ISCD) has developed new official positions for the clinical use of quantitative computed tomography (QCT)-based finite element analysis of the spine and hip. The ISCD task force for QCT reviewed the evidence for clinical applications and presented a report with recommendations at the 2015 ISCD Position Development Conference. Here we discuss the agreed upon ISCD official positions with supporting medical evidence, rationale, controversy, and suggestions for further study. Parts I and III address the clinical use of QCT of the hip, and the clinical feasibility of existing techniques for opportunistic screening of osteoporosis using CT scans obtained for other diagnosis such as colonography was addressed.

Biomechanical validation of computer assisted planning of periacetabular osteotomy: A preliminary study based on finite element analysis

Periacetabular Osteotomy (PAO) is a joint preserving surgical intervention intended to increase femoral head coverage and thereby to improve stability in young patients with hip dysplasia. Previously, we developed a CT-based, computer-assisted program for PAO diagnosis and planning, which allows for quantifying the 3D acetabular morphology with parameters such as acetabular version, inclination, lateral center edge (LCE) angle and femoral head coverage ratio (CO). In order to verify the hypothesis that our morphology-based planning strategy can improve biomechanical characteristics of dysplastic hips, we developed a 3D finite element model based on patient-specific geometry to predict cartilage contact stress change before and after morphology-based planning. Our experimental results demonstrated that the morphology-based planning strategy could reduce cartilage contact pressures and at the same time increase contact areas. In conclusion, our computer-assisted system is an efficient tool for PAO planning.

Clay mineralogy and multi-element chemistry of surface sediments on the Aiberian-Arctic Shelf

Clay mineral and bulk chemical (Si, Al, K, Mg, Sr, La, Ce, Nd) analyses of terrigenous surface sediments on the Siberian-Arctic shelf indicate that there are five regions with distinct, or endmember, sedimentary compositions. The formation of these geochemical endmembers is controlled by sediment provenance and grain size sorting. (1) The shale endmember (Al, K and REE rich sediment) is eroded from fine-grained marine sedimentary rocks of the Verkhoyansk Mountains and Kolyma-Omolon superterrain, and discharged to the shelf by the Lena, Yana, Indigirka and Kolyma Rivers. (2) The basalt endmember (Mg rich) originates from NE Siberia's Okhotsk-Chukotsk volcanic belt and Bering Strait inflow, and is prevalent in Chukchi Sea Sediments. Concentrations of the volcanically derived clay mineral smectite are elevated in Chukchi fine-fraction sediments, corroborating the conclusion that Chukchi sediments are volcanic in origin. (3) The mature sandstone endmember (Si rich) is found proximal to Wrangel Island and sections of the Chukchi Sea's Siberian coast and is derived from the sedimentary Chukotka terrain that comprises these landmasses. (4) The immature sandstone endmember (Sr rich) is abundant in the New Siberian Island region and reflects inputs from sedimentary rocks that comprise the islands. (5) The immature sandstone endmember is also prevalent in the western Laptev Sea, where it is eroded from sedimentary deposits blanketing the Siberian platform that are compositionally similar to those on the New Siberian Islands. Western Laptev can be distinguished from New Siberian Island region sediments by their comparatively elevated smectite concentrations and the presence of the basalt endmember, which indicate Siberian platform flood basalts are also a source of western Laptev sediments. In certain locations grain size sorting noticeably affects shelf sediment chemistry. (1) Erosion of fines by currents and sediment ice rafting contributes to the formation of the coarse-grained sandstone endmembers. (2) Bathymetrically controlled grain size sorting, in which fines preferentially accumulate offshore in deeper, less energetic water, helps distribute the fine-grained shale and basalt endmembers. An important implication of these results is that the observed sedimentary geochemical endmembers provide new markers of sediment provenance, which can be used to track sediment transport, ice-rafted debris dispersal or the movement of particle-reactive contaminants.

Geochemical processes and rare earth element analysis in cold-seep pore waters of Hydrate Ridge, northeast Pacific Ocean

The concentrations of rare earth elements (REEs), sulphate, hydrogen sulphide, total alkalinity, calcium, magnesium and phosphate were measured in shallow (<12 cm below seafloor) pore waters from cold-seep sediments on the northern and southern summits of Hydrate Ridge, offshore Oregon. Downward-decreasing sulphate and coevally increasing sulphide concentrations reveal sulphate reductionas dominant early diagenetic process from ~2 cm depth downwards. A strong increase of total dissolved REE concentrations is evident immediately below the sediment-water interface, which can be related to early diagenetic release of REEs into pore water resulting from the remineralization of particulate organic matter. The highest pore water REE concentrations were measured close to the sediment-water interface at ~2 cm depth. Distinct shale normalized REE patterns point to particulate organic matter and iron oxides as main REE sources in the upper ~2-cm depth interval. In general, the pore waters have shalenormalized patterns reflecting heavy REE (HREE) enrichment, which suggests preferential complexation of HREEs with carbonate ions. Below ~2 cm depth, a downward decrease in REE correlates with a decrease in pore water calcium concentrations. At this depth, the anaerobic oxidation of methane (AOM) coupled to sulphate reduction increases carbonate alkalinity through the production of bicarbonate, which results in the precipitation of carbonate minerals. It seems therefore likely that the REEs and calcium are consumed during vast AOM-induced precipitation of carbonate in shallow Hydrate Ridge sediments. The analysis of pore waters from Hydrate Ridge shed new light on early diagenetic processes at cold seeps, corroborating the great potential of REEs to identify geochemical processes and to constrain environmental conditions.