skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Finite element modeling with subject‐specific mechanical properties to assess knee osteoarthritis initiation and progression
Finite element models of the knee can be used to identify regions at risk of mechanical failure in studies of osteoarthritis. Models of the knee often implement joint geometry obtained from magnetic resonance imaging (MRI) or gait kinematics from motion capture to increase model specificity for a given subject. However, differences exist in cartilage material properties regionally as well as between subjects. This paper presents a method to create subject-specific finite element models of the knee that assigns cartilage material properties from T2 relaxometry. We compared our T2-refined model to identical models with homogeneous material properties. When tested on three subjects from the Osteoarthritis Initiative data set, we found the T2-refined models estimated higher principal stresses and shear strains in most cartilage regions and corresponded better to increases in KL grade in follow-ups compared to their corresponding homogeneous material models. Measures of cumulative stress within regions of a T2-refined model also correlated better with the region's cartilage morphology MRI Osteoarthritis Knee Score as compared with the homogeneous model. We conclude that spatially heterogeneous T2-refined material properties improve the subject-specificity of finite element models compared to homogeneous material properties in osteoarthritis progression studies. Statement of Clinical Significance: T2-refined material properties can improve subject-specific finite element model assessments of cartilage degeneration.  more » « less
Award ID(s):
2149946
PAR ID:
10329835
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Orthopaedic Research
ISSN:
0736-0266
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. BackgroundHealthy articular cartilage presents structural gradients defined by distinct zonal patterns through the thickness, which may be disrupted in the pathogenesis of several disorders. Analysis of textural patterns using quantitative MRI data may identify structural gradients of healthy or degenerating tissue that correlate with early osteoarthritis (OA). PurposeTo quantify spatial gradients and patterns in MRI data, and to probe new candidate biomarkers for early severity of OA. Study TypeRetrospective study. SubjectsFourteen volunteers receiving total knee replacement surgery (eight males/two females/four unknown, average age ± standard deviation: 68.1 ± 9.6 years) and 10 patients from the OA Initiative (OAI) with radiographic OA onset (two males/eight females, average age ± standard deviation: 57.7 ± 9.4 years; initial Kellgren‐Lawrence [KL] grade: 0; final KL grade: 3 over the 10‐year study). Field Strength/Sequence3.0‐T and 14.1‐T, biomechanics‐based displacement‐encoded imaging, fast spin echo, multi‐slice multi‐echoT2mapping. AssessmentWe studied structure and strain in cartilage explants from volunteers receiving total knee replacement, or structure in cartilage of OAI patients with progressive OA. We calculated spatial gradients of quantitative MRI measures (eg, T2) normal to the cartilage surface to enhance zonal variations. We compared gradient values against histologically OA severity, conventional relaxometry, and/or KL grades. Statistical TestsMultiparametric linear regression for evaluation of the relationship between residuals of the mixed effects models and histologically determined OA severity scoring, with a significance threshold atα = 0.05. ResultsGradients of individual relaxometry and biomechanics measures significantly correlated with OA severity, outperforming conventional relaxometry and strain metrics. In human explants, analysis of spatial gradients provided the strongest relationship to OA severity (R2 = 0.627). Spatial gradients of T2 from OAI data identified variations in radiographic (KL Grade 2) OA severity in single subjects, while conventional T2 alone did not. Data ConclusionSpatial gradients of quantitative MRI data may improve the predictive power of noninvasive imaging for early‐stage degeneration. Evidence Level1 Technical EfficacyStage 1 
    more » « less
  2. 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. 
    more » « less
  3. In the medical sector, three-dimensional (3D) images are commonly used like computed tomography (CT) and magnetic resonance imaging (MRI). The 3D MRI is a non-invasive method of studying the soft-tissue structures in a knee joint for osteoarthritis studies. It can greatly improve the accuracy of segmenting structures such as cartilage, bone marrow lesion, and meniscus by identifying the bone structure first. U-net is a convolutional neural network that was originally designed to segment the biological images with limited training data. The input of the original U-net is a single 2D image and the output is a binary 2D image. In this study, we modified the U-net model to identify the knee bone structures using 3D MRI, which is a sequence of 2D slices. A fully automatic model has been proposed to detect and segment knee bones. The proposed model was trained, tested, and validated using 99 knee MRI cases where each case consists of 160 2D slices for a single knee scan. To evaluate the model’s performance, the similarity, dice coefficient (DICE), and area error metrics were calculated. Separate models were trained using different knee bone components including tibia, femur, patella, as well as a combined model for segmenting all the knee bones. Using the whole MRI sequence (160 slices), the method was able to detect the beginning and ending bone slices first, and then segment the bone structures for all the slices in between. On the testing set, the detection model accomplished 98.79% accuracy and the segmentation model achieved DICE 96.94% and similarity 93.98%. The proposed method outperforms several state-of-the-art methods, i.e., it outperforms U-net by 3.68%, SegNet by 14.45%, and FCN-8 by 2.34%, in terms of DICE score using the same dataset. 
    more » « less
  4. Finite element analysis is used to study brain axonal injury and develop Brain White Matter (BWM) models while accounting for both the strain magnitude and the strain rate. These models are becoming more sophisticated and complicated due to the complex nature of the BMW composite structure with different material properties for each constituent phase. State-of-the-art studies focus on employing techniques that combine information about the local axonal directionality in different areas of the brain with diagnostic tools such as Diffusion-Weighted Magnetic Resonance Imaging (Diffusion-MRI). The diffusion-MRI data offers localization and orientation information of axonal tracks which are analyzed in finite element models to simulate virtual loading scenarios. Here, a BMW biphasic material model comprised of axons and neuroglia is considered. The model’s architectural anisotropy represented by a multitude of axonal orientations, that depend on specific brain regions, adds to its complexity. During this effort, we develop a finite element method to merge micro-scale Representative Volume Elements (RVEs) with orthotropic frequency domain viscoelasticity to an integrated macro-scale BWM finite element model, which incorporates local axonal orientation. Previous studies of this group focused on building RVEs that combined different volume fractions of axons and neuroglia and simulating their anisotropic viscoelastic properties. Via the proposed model, we can assign material properties and local architecture on each element based on the information from the orientation of the axonal traces. Consecutively, a BWM finite element model is derived with fully defined both material properties and material orientation. The frequency-domain dynamic response of the BMW model is analyzed to simulate larger scale diagnostic modalities such as MRI and MRE. 
    more » « less
  5. Abstract Background We aimed to determine if composite structural measures of knee osteoarthritis (KOA) progression on magnetic resonance (MR) imaging can predict the radiographic onset of accelerated knee osteoarthritis. Methods We used data from a nested case-control study among participants from the Osteoarthritis Initiative without radiographic KOA at baseline. Participants were separated into three groups based on radiographic disease progression over 4 years: 1) accelerated (Kellgren-Lawrence grades [KL] 0/1 to 3/4), 2) typical (increase in KL, excluding accelerated osteoarthritis), or 3) no KOA (no change in KL). We assessed tibiofemoral cartilage damage (four regions: medial/lateral tibia/femur), bone marrow lesion (BML) volume (four regions: medial/lateral tibia/femur), and whole knee effusion-synovitis volume on 3 T MR images with semi-automated programs. We calculated two MR-based composite scores. Cumulative damage was the sum of standardized cartilage damage. Disease activity was the sum of standardized volumes of effusion-synovitis and BMLs. We focused on annual images from 2 years before to 2 years after radiographic onset (or a matched time for those without knee osteoarthritis). To determine between group differences in the composite metrics at all time points, we used generalized linear mixed models with group (3 levels) and time (up to 5 levels). For our prognostic analysis, we used multinomial logistic regression models to determine if one-year worsening in each composite metric change associated with future accelerated knee osteoarthritis (odds ratios [OR] based on units of 1 standard deviation of change). Results Prior to disease onset, the accelerated KOA group had greater average disease activity compared to the typical and no KOA groups and this persisted up to 2 years after disease onset. During a pre-radiographic disease period, the odds of developing accelerated KOA were greater in people with worsening disease activity [versus typical KOA OR (95% confidence interval [CI]): 1.58 (1.08 to 2.33); versus no KOA: 2.39 (1.55 to 3.71)] or cumulative damage [versus typical KOA: 1.69 (1.14 to 2.51); versus no KOA: 2.11 (1.41 to 3.16)]. Conclusions MR-based disease activity and cumulative damage metrics may be prognostic markers to help identify people at risk for accelerated onset and progression of knee osteoarthritis. 
    more » « less