Search

Material property assignment is a critical step in developing subject-specific finite element models of bone. Inhomogeneous material properties are often applied using an equation relating density and elastic modulus, with the density information coming from CT scans of the bone. Very few previous studies have investigated which density-elastic modulus relationships from the literature are most suitable for application in long bone. No such studies have been completed for the ulna. The purpose of this study was to investigate six such density-modulus relationships and compare the results to experimental strains from eight cadaveric ulnae. Subject-specific finite element models were developed for each bone using micro-CT scans. Six density-modulus equations were trialed in each bone, resulting in a total of 48 models. Data from a previously completed experimental study in which each bone was instrumented with twelve strain gauges were used for comparison. Although the relationship that best matched experimental strains was somewhat specimen and location dependent, there were two relations which consistently matched the experimental strains most closely. One of these under-estimated and one over-estimated the experimental strain values, by averages of 15% and 31%, respectively. The results of this study suggest that the ideal relationship for the ulna may lie somewhere in between these two relations.

Diffuse idiopathic skeletal hyperostosis (DISH) is a non‐inflammatory spondyloarthropathy identified radiographically by calcification of the ligaments and/or entheses along the anterolateral aspect of the vertebral column. The etiology and pathogenesis of calcifications are unknown, and the diagnosis of DISH is currently based on radiographic criteria associated with advanced disease. To characterize the features of calcifications associated with DISH, we used micro‐computed tomographic imaging to evaluate a cohort of 19 human cadaveric vertebral columns. Fifty‐three percent of the cohort (n = 10; 3 females, 7 males, mean age of death = 81 years, range 67–94) met the radiographic criteria for DISH, with calcification of four or more contiguous vertebral segments. In almost all cases, the lower thoracic regions (T8‐12) were affected by calcifications, consisting primarily of large, horizontal outgrowths of bony material. In contrast, calcifications localized to the upper thoracic regions demonstrated variability in their presentation and were categorized as either "continuous vertical bands" or "discontinuous‐patchy" lesions. In addition to the variable morphology of the calcifications, our analysis demonstrated remarkable heterogeneity in the densities of calcifications, ranging from internal components below the density of cortical bone to regions of hyper‐dense material that exceeded cortical bone. These findings establish that the current radiographic criteria for DISH capture heterogeneous presentations of ectopic spine calcification that can be differentiated based on morphology and density. These findings may indicate a naturally heterogenous disease, potential stage(s) in the natural progression of DISH, or distinct pathologies of ectopic calcifications. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res [ABSTRACT FROM AUTHOR]

Objective: To assess whether the transient stresses of foot impact with the ground are similar to those found during midstance loading and if the location of high stress correlate with the sites most commonly associated with mechanically induced osteoarthritis (OA). We compared impact stresses in subchondral bone between two subject-specific, three-dimensional, finite-element models of the equine metacarpophalangeal (MCP) joint—one with advanced OA and one healthy, and with similar published data on the stresses that occur at midstance. Methods: Two right MCP joints (third metacarpal and proximal phalanx) were scanned using micro-computed tomography (μCT). Images were segmented, and meshed using modified 10-node quadratic tetrahedral elements. Bone material properties were assigned based on the bone density. An impact velocity of 3.55 m/s was applied to each model and contact pressures and stress distribution were calculated for each. In a separate iteration, the third metacarpal was loaded statically. A sampling grid of 160 equidistant points was superimposed over selected slices, and average peak stresses were calculated for 6 anatomical regions. Within-region maximal peak and average von Mises stresses were compared between healthy and OA bones in both midstance and impact loading. Results: Average impact stresses across all regions, in both locations (palmar and dorsal) were greater in the OA model. Highest impact stresses were located in the dorsal medial condyle in the healthy (12.8 MPa) and OA (14.1MPa) models, and were lowest in the palmar medial and lateral parasagittal grooves in the healthy (5.94 MPa) and OA (7.07 MPa) models. The healthy static model had higher peak (up to 49.7% greater) and average (up to 38.6% greater) stresses in both locations and across all regions compared to the OA static model. Conclusions: Under simulated footfall a trot, loading on the dorsal aspect of the third metacarpal at impact created stresses similar to those found during midstance. The high accelerations that occur under impact loading are likely responsible for creating the high stresses, as opposed to midstance loading where the high stresses are the result of high mass loading. Although the stress magnitudes were found to be similar among the two loading conditions, the location of the high stress loading occurred in sites that are not typically associated with osteoarthritic changes. [ABSTRACT FROM AUTHOR]

Articular cartilage degeneration is the most consistently observed feature of osteoarthritis (OA). Animal and human studies have shown that various forms of exercise influence the course of the disease in different ways. In addition, early changes in articular cartilage that influence the progression of OA, such as the expression of cytokines, require further investigation. We have used a surgically induced experimental model of knee OA to address these questions. Here, we discuss our recent studies investigating the effects of an exercise paradigm in surgically induced OA, which determined that the destabilized knee joint is susceptible to enhanced degeneration when subjected to low-intensity, low-impact exercise. Further, we investigated early global changes in gene expression in articular chondrocytes from degenerating cartilage. Identified candidate genes including genes involved in chemokine, endothelin, and transforming growth factor-α signaling are discussed in the context of articular cartilage degeneration in early OA. [ABSTRACT FROM AUTHOR]

Background Quantitative computed tomography based subject-specific finite element modeling has potential to clarify the role of subchondral bone alterations in knee osteoarthritis initiation, progression, and pain. However, it is unclear what density-modulus equation(s) should be applied with subchondral cortical and subchondral trabecular bone when constructing finite element models of the tibia. Using a novel approach applying neural networks, optimization, and back-calculation against in situ experimental testing results, the objective of this study was to identify subchondral-specific equations that optimized finite element predictions of local structural stiffness at the proximal tibial subchondral surface. Methods Thirteen proximal tibial compartments were imaged via quantitative computed tomography. Imaged bone mineral density was converted to elastic moduli using multiple density-modulus equations (93 total variations) then mapped to corresponding finite element models. For each variation, root mean squared error was calculated between finite element prediction and in situ measured stiffness at 47 indentation sites. Resulting errors were used to train an artificial neural network, which provided an unlimited number of model variations, with corresponding error, for predicting stiffness at the subchondral bone surface. Nelder-Mead optimization was used to identify optimum density-modulus equations for predicting stiffness. Findings Finite element modeling predicted 81% of experimental stiffness variance (with 10.5% error) using optimized equations for subchondral cortical and trabecular bone differentiated with a 0.5 g/cm 3 density. Interpretation In comparison with published density-modulus relationships, optimized equations offered improved predictions of local subchondral structural stiffness. Further research is needed with anisotropy inclusion, a smaller voxel size and de-blurring algorithms to improve predictions. [ABSTRACT FROM AUTHOR]

Background Quantitative computed tomography based subject-specific finite element modeling has potential to clarify the role of subchondral bone alterations in knee osteoarthritis initiation, progression, and pain initiation. Calculation of bone elastic moduli from image data is a basic step when constructing finite element models. However, different relationships between elastic moduli and imaged density (known as density–modulus relationships) have been reported in the literature. The objective of this study was to apply seven different trabecular-specific and two cortical-specific density–modulus relationships from the literature to finite element models of proximal tibia subchondral bone, and identify the relationship(s) that best predicted experimentally measured local subchondral structural stiffness with highest explained variance and least error. Methods Thirteen proximal tibial compartments were imaged via quantitative computed tomography. Imaged bone mineral density was converted to elastic moduli using published density–modulus relationships and mapped to corresponding finite element models. Proximal tibial structural stiffness values were compared to experimentally measured stiffness values from in-situ macro-indentation testing directly on the subchondral bone surface (47 indentation points). Findings Regression lines between experimentally measured and finite element calculated stiffness had R 2 values ranging from 0.56 to 0.77. Normalized root mean squared error varied from 16.6% to 337.6%. Interpretation Of the 21 evaluated density–modulus relationships in this study, Goulet combined with Snyder and Schneider or Rho appeared most appropriate for finite element modeling of local subchondral bone structural stiffness. Though, further studies are needed to optimize density–modulus relationships and improve finite element estimates of local subchondral bone structural stiffness. [ABSTRACT FROM AUTHOR]