Subject-specific Finite Element Models Research Articles

Biomechanical analysis of various internal fracture fixation devices used for treating femoral neck fractures: A comparative finite element analysis

IntroductionSeveral internal fixation devices are available for treating Pauwels type I, II and III femoral neck fractures. The present study compared various fixation implants for all Pauwels fracture types using a CT-based subject-specific finite element model of the femur and determined the most effective implant for each fracture type. Materials and MethodsThe analysis included four different configurations of cannulated screw models, Femoral Neck System, Dynamic Hip Screw and Dynamic Condylar Screw (with and without anti-rotational screw). Ti-alloy was considered as the implant material. Heterogeneous bone material property was assigned based on CT grey value. Frictional contact was assumed in the contact interfaces. Peak loading corresponding to normal walking and stair-climbing were considered. Equivalent strain in bone, equivalent stress in the implants, femoral head deformation and rotation, micromotion in the contact interfaces, and strain-shielding in bone were evaluated for each implanted model. ResultsStresses generated in the implants were within the yield limit of the implant material. In Pauwels I and II, the micromotion predicted at the contact regions in all the implanted models was within 100 µm, which is suitable for bone integration. However, in Pauwels III fracture, most of the implanted models other than DHS with AR-screw model exhibited micromotion of more than 150 µm in the contact regions, which is expected to inhibit bone growth. ConclusionsThe DHS with AR-screw implanted model was identified as the most effective in treating Pauwels I and III fractures. However, for Pauwels type II, DCS with an AR-screw implant was deemed superior to the other configurations.

Validation and Estimation of Obesity-Induced Intervertebral Disc Degeneration through Subject-Specific Finite Element Modelling of Functional Spinal Units.

(1) Background: Intervertebral disc degeneration has been linked to obesity; its potential mechanical effects on the intervertebral disc remain unknown. This study aimed to develop and validate a patient-specific model of L3-L4 vertebrae and then use the model to estimate the impact of increasing body weight on disc degeneration. (2) Methods: A three-dimensional model of the functional spinal unit of L3-L4 vertebrae and its components were developed and validated. Validation was achieved by comparing the range of motions (RoM) and intradiscal pressures with the previous literature. Subsequently, the validated model was loaded according to the body mass index and estimated stress, deformation, and RoM to assess disc degeneration. (3) Results: During validation, L3-L4 RoM and intradiscal pressures: flexion 5.17° and 1.04 MPa, extension 1.54° and 0.22 MPa, lateral bending 3.36° and 0.54 MPa, axial rotation 1.14° and 0.52 MPa, respectively. When investigating the impact of weight on disc degeneration, escalating from normal weight to obesity reveals an increased RoM, by 3.44% during flexion, 22.7% during extension, 29.71% during lateral bending, and 33.2% during axial rotation, respectively. Also, stress and disc deformation elevated with increasing weight across all RoM. (4) Conclusions: The predicted mechanical responses of the developed model closely matched the validation dataset. The validated model predicts disc degeneration under increased weight and could lay the foundation for future recommendations aimed at identifying predictors of lower back pain due to disc degeneration.

Load transfer in bone after partial, multi-compartmental, and total knee arthroplasty.

Introduction: Arthroplasty-associated bone loss remains a clinical problem: stiff metallic implants disrupt load transfer to bone and, hence, its remodeling stimulus. The aim of this research was to analyze how load transfer to bone is affected by different forms of knee arthroplasty: isolated partial knee arthroplasty (PKA), compartmental arthroplasty [combined partial knee arthroplasty (CPKA), two or more PKAs in the same knee], and total knee arthroplasty (TKA). Methods: An experimentally validated subject-specific finite element model was analyzed native and with medial unicondylar, lateral unicondylar, patellofemoral, bi-unicondylar, medial bicompartmental, lateral bicompartmental, tricompartmental, and total knee arthroplasty. Three load cases were simulated for each: gait, stair ascent, and sit-to-stand. Strain shielding and overstraining were calculated from the differences between the native and implanted states. Results: For gait, the TKA femoral component led to mean strain shielding (30%) more than three times higher than that of PKA (4%-7%) and CPKA (5%-8%). Overstraining was predicted in the proximal tibia (TKA 21%; PKA/CPKA 0%-6%). The variance in the distribution for TKA was an order of magnitude greater than for PKA/CPKA, indicating less physiological load transfer. Only the TKA-implanted femur was sensitive to the load case: for stair ascent and gait, almost the entire distal femur was strain-shielded, whereas during sit-to-stand, the posterior femoral condyles were overstrained. Discussion: TKA requires more bone resection than PKA and CPKA. These finite element analyses suggest that a longer-term benefit for bone is probable as partial and multi-compartmental knee procedures lead to more natural load transfer compared to TKA. High-flexion activity following TKA may be protective of posterior condyle bone resorption, which may help explain why bone loss affects some patients more than others. The male and female bone models used for this research are provided open access to facilitate future research elsewhere.

Custom orthotic design by integrating 3D scanning and subject-specific FE modelling workflow

The finite element (FE) foot model can help estimate pathomechanics and improve the customized foot orthoses design. However, the procedure of developing FE models can be time-consuming and costly. This study aimed to develop a subject-specific scaled foot modelling workflow for the foot orthoses design based on the scanned foot surface data. Six participants (twelve feet) were collected for the foot finite element modelling. The subject-specific surface-based finite element model (SFEM) was established by incorporating the scanned foot surface and scaled foot bone geometries. The geometric deviations between the scaled and the scanned foot surfaces were calculated. The SFEM model was adopted to predict barefoot and foot-orthosis interface pressures. The averaged distances between the scaled and scanned foot surfaces were 0.23 ± 0.09 mm. There was no significant difference for the hallux, medial forefoot, middle forefoot, midfoot, medial hindfoot, and lateral hindfoot, except for the lateral forefoot region (p = 0.045). The SFEM model evaluated slightly higher foot-orthoses interface pressure values than measured, with a maximum deviation of 7.1%. These results indicated that the SFEM technique could predict the barefoot and foot-orthoses interface pressure, which has the potential to expedite the process of orthotic design and optimization.Graphical abstract

Accuracy and precision of internal displacement and strain measurements in long human bones using HR-pQCT and digital volume correlation

ABSTRACT Digital volume correlation (DVC) is a technique for measuring 3D, internal displacements and strains in loaded structures. One application of DVC is the study of internal bone mechanics and the validation of subject-specific finite element (FE) models. While micro-computed tomography (µCT) is a common choice for DVC studies of small samples of bone, long human bones such as the tibia may exceed the spatial limitations of conventional µCT. High resolution peripheral quantitative CT (HR-pQCT) scanners, with their large bore and open-ended design, may be a viable option for long-bone DVC. However, HR-pQCT is afflicted by stitching artefacts, caused by the stacking of scan blocks to form the large scan volume, resulting in erroneous steps in DVC displacements and bands of increased strain error. This study proposed a modified HR-pQCT scanning protocol and stitching methodology to mitigate the effects of stitching artefacts on DVC measurements of displacement and strain. With the application of the proposed scanning/stitching methodology, displacements greater than 11.7µm (0.29 voxels) and strains greater than 1633µε could be repeatedly measured by DVC. These results show that HR-pQCT combined with DVC is suitable for measuring internal bone displacements in long human bones with sub-voxel precision and strains greater than 1633µε.

Integration of Swin UNETR and statistical shape modeling for a semi-automated segmentation of the knee and biomechanical modeling of articular cartilage

Simulation studies, such as finite element (FE) modeling, provide insight into knee joint mechanics without patient involvement. Generic FE models mimic the biomechanical behavior of the tissue, but overlook variations in geometry, loading, and material properties of a population. Conversely, subject-specific models include these factors, resulting in enhanced predictive precision, but are laborious and time intensive. The present study aimed to enhance subject-specific knee joint FE modeling by incorporating a semi-automated segmentation algorithm using a 3D Swin UNETR for an initial segmentation of the femur and tibia, followed by a statistical shape model (SSM) adjustment to improve surface roughness and continuity. For comparison, a manual FE model was developed through manual segmentation (i.e., the de-facto standard approach). Both FE models were subjected to gait loading and the predicted mechanical response was compared. The semi-automated segmentation achieved a Dice similarity coefficient (DSC) of over 98% for both the femur and tibia. Hausdorff distance (mm) between the semi-automated and manual segmentation was 1.4 mm. The mechanical results (max principal stress and strain, fluid pressure, fibril strain, and contact area) showed no significant differences between the manual and semi-automated FE models, indicating the effectiveness of the proposed semi-automated segmentation in creating accurate knee joint FE models. We have made our semi-automated models publicly accessible to support and facilitate biomechanical modeling and medical image segmentation efforts (https://data.mendeley.com/datasets/k5hdc9cz7w/1).

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.

Analysis of the pattern of microstructural changes in the brain after mTBI with diffusion tensor imaging and subject-specific FE models

Traumatic brain injury (TBI) is a major public health challenge. Up to 90 % of TBIs are on the mild spectrum of TBI (mTBI), where diagnosis is a major challenge. Majority of studies in this field have been conducted on human subjects, which inherently suffer from the lack of appropriate control group, selection bias, and individual differences in patients. To overcome these limitations, animal studies have been used as an alternative approach to provide deeper insights into the underlying mechanism related to the injury. Therefore our aim is to investigate various quantitative imaging biomarkers acquired from T1-W and diffusion tensor imaging (DTI) sequences to provide more information about the microstructural changes in the brain after mTBI. We then use this to generate subject-specific finite element models of the brain and examine how the changes in the brain material properties reflected in MR images affects strain distribution patterns on a subsequent head hit. Our study revealed a decrease in FA and an increase in diffusivity indices (MD, AD, RD) in the white matter tracts of the brain. This finding may represent the axonal damage, demyelination and gliosis after mild TBI, which have been shown in other animal and human studies. Moreover, our FE analysis showed that microstructural changes in the brain after mTBI might have weakened the structural integrity of the brain as the subsequent head hit led to wider and more severe brain deformations. SignificanceAnimal models have been used to investigate biomechanical and pathophysiological aspects of mild traumatic brain injuries in the past. Still, most of them used small animals such as rats and mice. These models provided valuable insight into the pathophysiology of mTBI, but their findings have limitations due to their inherent differences to human brains. We have developed a large animal model of mTBI with sheep brains by combining advanced MRI and finite element analysis as they mimic the human brain better. To the best of our knowledge, this study is the first mTBI neuroimaging study conducted on large animal brains to investigate the diffusional changes in the white matter tracts after mTBI. Our FE analysis revealed that such microstructural changes resulted in tissue softening as the extent of brain deformation increased on a subsequent head hit, indicating increased brain vulnerability after head impacts.

The long head of the biceps tendon undergoes multiaxial deformation during shoulder motion

The long head biceps tendon (LHBT) is presumed a common source of shoulder joint pain and injury. Despite common LHBT pathologies, diagnosis and preferred treatment remain frequently debated. This Short Communication reports the development of a subject-specific finite element model of the shoulder joint based on one subject’s 3D reconstructed anatomy and 3D in vivo kinematics recorded from bone-fixed electromagnetic sensors. The primary purpose of this study was to use the developed finite element model to investigate the LHBT mechanical environment during a typical shoulder motion of arm raising. Furthermore, this study aimed to assess the viability of material models derived from uniaxial tensile tests for accurate simulation of in vivo motion. The findings of our simulations indicate that the LHBT undergoes complex multidimensional deformations. As such, uniaxial material properties reported in the existing body of literature are not sufficient to simulate accurately the in vivo mechanical behavior of the LHBT. Further experimental tests on cadaveric specimens, such as biaxial tension and combinations of tension and torsion, are needed to describe fully the mechanical behavior of the LHBT and investigate its mechanisms of injury.

Impact of tricuspid annuloplasty device shape and size on valve mechanics—a computational study

BackgroundTricuspid valve disease significantly affects 1.6 million Americans. The gold standard treatment for tricuspid disease is the implantation of annuloplasty devices. These ring-like devices come in various shapes and sizes. Choices for both shape and size are most often made by surgical intuition rather than scientific rationale. MethodsTo understand the impact of shape and size on valve mechanics and to provide a rational basis for their selection, we used a subject-specific finite element model to conduct a virtual case study. That is, we implanted 4 different annuloplasty devices of 6 different sizes in our virtual patient. After each virtual surgery, we computed the coaptation area, leaflet end-systolic angles, leaflet stress, and chordal forces. ResultsWe found that contoured devices are better at normalizing end-systolic angles, whereas the one flat device, the Edwards Classic, maximized the coaptation area and minimized leaflet stress and chordal forces. We further found that reducing device size led to increased coaptation area but also negatively impacted end-systolic angles, stress, and chordal forces. ConclusionsBased on our analyses of the coaptation area, leaflet motion, leaflet stress, and chordal forces, we found that device shape and size have a significant impact on valve mechanics. Thereby, our study also demonstrates the value of simulation tools and device tests in “virtual patients.” Expanding our study to many more valves may, in the future, allow for universal recommendations.

Validation and evaluation of subject-specific finite element models of the pediatric knee

Finite element (FE) models have been widely used to investigate knee joint biomechanics. Most of these models have been developed to study adult knees, neglecting pediatric populations. In this study, an atlas-based approach was employed to develop subject-specific FE models of the knee for eight typically developing pediatric individuals. Initially, validation simulations were performed at four passive tibiofemoral joint (TFJ) flexion angles, and the resulting TFJ and patellofemoral joint (PFJ) kinematics were compared to corresponding patient-matched measurements derived from magnetic resonance imaging (MRI). A neuromusculoskeletal-(NMSK)-FE pipeline was then used to simulate knee biomechanics during stance phase of walking gait for each participant to evaluate model simulation of a common motor task. Validation simulations demonstrated minimal error and strong correlations between FE-predicted and MRI-measured TFJ and PFJ kinematics (ensemble average of root mean square errors < 5 mm for translations and < 4.1° for rotations). The FE-predicted kinematics were strongly correlated with published reports (ensemble average of Pearson's correlation coefficients (ρ) > 0.9 for translations and ρ > 0.8 for rotations), except for TFJ mediolateral translation and abduction/adduction rotation. For walking gait, NMSK-FE model-predicted knee kinematics, contact areas, and contact pressures were consistent with experimental reports from literature. The strong agreement between model predictions and experimental reports underscores the capability of sequentially linked NMSK-FE models to accurately predict pediatric knee kinematics, as well as complex contact pressure distributions across the TFJ articulations. These models hold promise as effective tools for parametric analyses, population-based clinical studies, and enhancing our understanding of various pediatric knee injury mechanisms. They also support intervention design and prediction of surgical outcomes in pediatric populations.

Automatic generation of subject-specific finite element models of the spine from magnetic resonance images.

The generation of subject-specific finite element models of the spine is generally a time-consuming process based on computed tomography (CT) images, where scanning exposes subjects to harmful radiation. In this study, a method is presented for the automatic generation of spine finite element models using images from a single magnetic resonance (MR) sequence. The thoracic and lumbar spine of eight adult volunteers was imaged using a 3D multi-echo-gradient-echo sagittal MR sequence. A deep-learning method was used to generate synthetic CT images from the MR images. A pre-trained deep-learning network was used for the automatic segmentation of vertebrae from the synthetic CT images. Another deep-learning network was trained for the automatic segmentation of intervertebral discs from the MR images. The automatic segmentations were validated against manual segmentations for two subjects, one with scoliosis, and another with a spine implant. A template mesh of the spine was registered to the segmentations in three steps using a Bayesian coherent point drift algorithm. First, rigid registration was applied on the complete spine. Second, non-rigid registration was used for the individual discs and vertebrae. Third, the complete spine was non-rigidly registered to the individually registered discs and vertebrae. Comparison of the automatic and manual segmentations led to dice-scores of 0.93-0.96 for all vertebrae and discs. The lowest dice-score was in the disc at the height of the implant where artifacts led to under-segmentation. The mean distance between the morphed meshes and the segmentations was below 1mm. In conclusion, the presented method can be used to automatically generate accurate subject-specific spine models.

Influence of stem design parameters on periprosthetic femoral fractures examined by subject specific finite element analyses.

Due to the increasing number of periprosthetic femoral fractures (PFF), the optimisation of implant design gains importance. For the presented research a validated, subject specific finite element model of a human femur with an inlying total hip stem was used to compare the influence of different geometrical implant parameters on the development of PFF. The heterogeneous bone tissue was modelled on the basis of computed tomography scans. A ductile damage model with element deletion was applied to simulate bone fracture in a load case re-enacting a stumbling scenario. The results were compared in terms of fracture load, subsidence and fracture pattern to analyse the influence of friction at the implant-bone interface, implant size and stem length. The results showed that higher friction coefficients lead to an increase of fracture load. Also, the usage of an oversized implant has a negligible effect while an undersized implant reduces the fracture load by 48.9% for the investigated femur. Lastly, a higher fracture load was reached with an elongated stem, but the bending and change in fracture path indicate a more distal force transmission and subsequent stress shielding in the proximal femur.

Cervical facet capsular ligament mechanics: Estimations based on subject-specific anatomy and kinematics.

To understand the facet capsular ligament's (FCL) role in cervical spine mechanics, the interactions between the FCL and other spinal components must be examined. One approach is to develop a subject-specific finite element (FE) model of the lower cervical spine, simulating the motion segments and their components' behaviors under physiological loading conditions. This approach can be particularly attractive when a patient's anatomical and kinematic data are available. We developed and demonstrated methodology to create 3D subject-specific models of the lower cervical spine, with a focus on facet capsular ligament biomechanics. Displacement-controlled boundary conditions were applied to the vertebrae using kinematics extracted from biplane videoradiography during planar head motions, including axial rotation, lateral bending, and flexion-extension. The FCL geometries were generated by fitting a surface over the estimated ligament-bone attachment regions. The fiber structure and material characteristics of the ligament tissue were extracted from available human cervical FCL data. The method was demonstrated by application to the cervical geometry and kinematics of a healthy 23-year-old female subject. FCL strain within the resulting subject-specific model were subsequently compared to models with generic: (1) geometry, (2) kinematics, and (3) material properties to assess the effect of model specificity. Asymmetry in both the kinematics and the anatomy led to asymmetry in strain fields, highlighting the importance of patient-specific models. We also found that the calculated strain field was largely independent of constitutive model and driven by vertebrae morphology and motion, but the stress field showed more constitutive-equation-dependence, as would be expected given the highly constrained motion of cervical FCLs. The current study provides a methodology to create a subject-specific model of the cervical spine that can be used to investigate various clinical questions by coupling experimental kinematics with multiscale computational models.

Stress-shielding resistant design of custom pelvic prostheses using lattice-based topology optimization

Endoprosthetic reconstruction of the pelvic bone using 3D-printed, custom-made implants has delivered early load-bearing ability and good functional outcomes in the short term to individuals with pelvic sarcoma. However, excessive stress-shielding and subsequent resorption of peri‑prosthetic bone can imperil the long-term stability of such implants. To evaluate the stress-shielding performance of pelvic prostheses, we developed a sequential modeling scheme using subject-specific finite element models of the pelvic bone-implant complex and personalized neuromusculoskeletal models for pre- and post-surgery walking. A new topology optimization approach is introduced for the stress-shielding resistant (SSR) design of custom pelvic prostheses, which uses 3D-printable porous lattice structures. The SSR optimization was applied to a typical pelvic prosthesis to reconstruct a type II+III bone resection. The stress-shielding performance of the optimized implant based on the SSR approach was compared against the conventional optimization. The volume of the peri‑prosthetic bone predicted to undergo resorption post-surgery decreased from 44 to 18%. This improvement in stress-shielding resistance was achieved without compromising the structural integrity of the prosthesis. The SSR design approach has the potential to improve the long-term stability of custom-made pelvic prostheses.

Human motion analysis and measurement techniques: current application and developing trend

Human motion analysis and measurement technology have been widely used in the fields of medical treatment, sports science, and rehabilitation. In clinical practice, motion analysis has been applied in the diagnosis and individualized treatment planning of various musculoskeletal diseases, and it is also an important objective scientific method to evaluate the therapeutic effect and the effectiveness of medical equipment. This study aimed to introduce the common modern motion capture measurement technology and equipment, the clinical application and limitations of motion analysis, and the possible development trend of motion analysis measuring techniques in the future. Motion analysis and measurement systems and medical image measurement and analysis technology have made landmark improvements over the past few decades in terms of orthopaedical biomechanics. Nevertheless, limitations still exist, both subjective and objective. All these drawbacks have promoted the exploration of the integrated methods that have now been widely used in motion analysis. The results of the case study about the subject-specific finite element modeling of the foot and sports shoe complex have also shown great consistency. Nevertheless, several possible future directions for motion analysis measuring techniques still exist. In the future, the progress of motion analysis and measurement methods will simultaneously drive the progress of orthopedics, rehabilitation, precision personalized medicine, and medical engineering.

The effect of thickness and elastic modulus of the anterior talofibular ligament on anterior ankle joint stiffness: A subject-specific finite element study.

Ankle sprain is a frequent type of sports injury leading to lateral ligament injury. The anterior talofibular ligament (ATFL) is a primary ligamentous stabilizer of the ankle joint and typically the most vulnerable ligament injured in a lateral ankle sprain (LAS). This study aimed to quantitively investigate the effect of the thickness and elastic modulus of ATFL on anterior ankle joint stiffness (AAJS) by developing nine subject-specific finite element (FE) models under acute injury, chronic injury, and control conditions of ATFL. A 120N forward force was applied at the posterior calcaneus leading to an anterior translation of the calcaneus and talus to simulate the anterior drawer test (ADT). In the results, the ratio of the forward force to the talar displacement was used to assess the AAJS, which increased by 5.85% in the acute group and decreased by 19.78% in the chronic group, compared to those of the control group. An empirical equation described the relationship between AAJS, thickness, and elastic modulus (R-square 0.98). The equation proposed in this study provided an approach to quantify AAJS and revealed the effect of the thickness and the elastic modulus of ATFL on ankle stability, which may shed light on the potential diagnosis of lateral ligament injury.

Computational comparison of medializing tibial tubercle osteotomy and trochleoplasty in patients with trochlear dysplasia.

Medial patellofemoral ligament reconstruction (MPFLR) has emerged as the procedure of choice for recurrent patellar dislocation. This addresses soft tissue injury but does not address underlying anatomic factors, including trochlear dysplasia, that are commonly present and increase risk of dislocation. Quantification of the stability offered by other surgical interventions, namely, medializing tibial tubercle osteotomy (mTTO) and trochleoplasty, with and without MPFLR, may provide insight for surgical choices in patients with trochlear dysplasia. We developed subject-specific finite element models based on magnetic resonance scans from a cohort of 20 patients with trochlear dysplasia and recurrent patellar dislocation. The objectives of this study were (1) to compare patella stability after mTTO and trochleoplasty procedures; (2) to evaluate whether it is necessary to perform an MPFLR in combination with the mTTO or trocheoplasty procedure; and (3) to quantify the robustness of patellar stability to variability in knee kinematics. Trochleoplasty performed better than mTTO at stabilizing the patella between 5° and 30° flexion. For both mTTO and trochleoplasty procedures, it was beneficial to also perform MPFLR-inclusion of MPFLR halved the magnitude of patellar laxity predicted in the simulations. Simulations that did not include any medial patellofemoral ligament restraintwere also more sensitive to variation in tibiofemoral internal-external kinematics.

COMPARTMENTAL FEMUR CORTICAL THICKNESS IN OLDER ADULTS DIFFERS BY DEMOGRAPHICS AND PHYSICAL FUNCTION

Abstract Previous studies have examined the relationship between cortical bone properties and fracture risk in older populations, but they have yet to examine how these properties change within the femur. Most hip fractures are classified by their location: femoral neck, intertrochanteric, or subtrochanteric fractures. Since cortical thickness can vary throughout the femur, a quantitative method of examining thickness by region was developed using a cortical mapping approach applied to computed tomography (CT) scans. Subject-specific finite element (FE) models of proximal femurs were created from baseline CT scans of 107 older adults (ages 60 – 85; 70% White; 72% Female) with obesity as classified by BMI (33.8 ± 4 kg/m2). An existing FE model was morphed to the segmented geometry of each subject’s proximal femur. A nearest neighbor search assigned cortical thickness values to the nearest finite element model node. Cortical thickness was grouped into four femoral compartments: femoral head, femoral neck, intertrochanteric and subtrochanteric regions. Pairwise paired t-tests indicated that cortical thickness differed between femoral compartments (p &amp;lt; 0.05). Multivariate regression models showed greater cortical thicknesses in femoral head, neck, and intertrochanteric regions in African Americans compared to Whites (p &amp;lt; 0.05). Additionally, these models suggest an association between cortical thickness and physical function assessments, such as the timed up-and-go test and leg muscle strength test. Since cortical thickness is not constant throughout the femur, future studies can use the framework developed here to assess each compartment individually and investigate the relationships between cortical thickness, demographics, and functional assessments.

MR-based quantitative measurement of human soft tissue internal strains for pressure ulcer prevention.

Pressure ulcers are a severe disease affecting patients that are bedridden or in a wheelchair bound for long periods of time. These wounds can develop in the deep layers of the skin of specific parts of the body, mostly on heels or sacrum, making them hard to detect in their early stages. Strain levels have been identified as a direct danger indicator for triggering pressure ulcers. Prevention could be possible with the implementation of subject-specific Finite Element (FE) models. However, generation and validation of such FE models is a complex task, and the current implemented techniques offer only a partial solution of the entire problem considering only external displacements and pressures, or cadaveric samples. In this paper, we propose an in vivo solution based on the 3D non-rigid registration between two Magnetic Resonance (MR) images, one in an unloaded configuration and the other deformed by means of a plate or an indenter. From the results of the image registration, the displacement field and subsequent strain maps for the soft tissues were computed. An extensive study, considering different cases (on heel pad and sacrum regions) was performed to evaluate the reproducibility and accuracy of the results obtained with this methodology. The implemented technique can give insight for several applications. It adds a useful tool for better understanding the propagation of deformations in the heel soft tissues that could generate pressure ulcers. This methodology can be used to obtain data on the material properties of the soft tissues to define constitutive laws for FE simulations and finally it offers a promising technique for validating FE models.