Risk of fracture modelling in healthy and osteoporotic tibia

Magnus Gislason, Sylvie Coupaud, Keisuke Sasagawa, Yuji Tanabe, Mariel Purcell, David Allan, Elizabeth Tanner

Research output: Contribution to conferenceProceeding

Abstract

Background: Osteoporosis has detrimental effects on the structural integrity of the skeletal system. The degeneration of bone material can be seen in spinal cord injury (SCI) patients due to lack of mechanical stimulation following extensive muscle paralysis. Aim: The aim of the study was to simulate and compare the mechanical response of the tibia of both able bodied (AB) subjects and SCI patients under four different types of loading, compression, A-P bending, M-L bending and torsion, and map how the risk of fracture (RF) correlated with bone mineral density (BMD) and type of loading. The stress distribution was calculated using the finite element method. Material and methods: Eight male SCI subjects with varying degrees of bone loss and four age-matched AB male control subjects were recruited for the study. The bones were scanned in vivo at the Queen Elizabeth National Spinal Injuries Unit (Glasgow) using a multi-slice peripheral Quantitative Computed Tomography, and individualised three dimensional finite element models were created of each subject's tibia from the scans using Mimics (Materialise). The material assignment was based on the greyscale pixel values of the scan images from which BMD was calculated. An empirical power law was used to estimate the Young’s modulus and the ultimate tensile strength for each element from the BMD. Results: Analysis of the material distribution in the tibia showed large numbers of elements with low BMD for the SCI group whereas a larger proportion of high density bone was seen in the AB group. The average BMD was lower at the epiphyses for the SCI group, but overall little difference was seen around the diaphysis between the groups. The finite element modelling was carried out using Abaqus (v.6.11, Simulia). The analysis was performed using 4 node tetrahedral elements and the element density was on average 2.9 elements/mm^3. The von Mises stress values of each element were exported and the percentage of elements exceeding an estimated RF threshold calculated. The results showed an increase in percentage of elements exceeding the RF threshold for the SCI group in bending and torsion compared to compression which is in agreement with clinical experience. On average, the SCI group showed greater numbers of elements exceeding the RF threshold than the AB group for all four loading conditions. The majority of the fractured elements were located around the epiphyses for all subjects. Conclusions and clinical implications: Using the analysis methods described can provide clinicians with useful patient-specific information about the amount of loading that can be put safely on the osteoporotic bone. That is of importance in SCI patient rehabilitation where musculoskeletal rehabilitation techniques include stimulating the bone via mechanical methods such as vibration or electrical stimulation of the muscles. Using such subject-specific analysis it is possible to determine the location along the bone that is most likely to fail along with the risk associated.

Conference

ConferenceEuropean Orthopaedic Research Society's 22nd Annual Meeting
CountryFrance
CityNantes
Period2/07/144/07/14

Fingerprint

Tibia
Spinal Cord Injuries
Bone
Bone Density
Bone and Bones
Minerals
Epiphyses
Rehabilitation
Diaphyses
Spinal Injuries
Patient rehabilitation
Muscles
Torsional stress
Muscle
Tensile Strength
Elastic Modulus
Compaction
Vibration
Paralysis
Electric Stimulation

Keywords

  • osteoporosis
  • skeletal system
  • spinal chord injury
  • SCI
  • patient rehabilitation
  • musculoskeletal rehabilitation

Cite this

Gislason, M., Coupaud, S., Sasagawa, K., Tanabe, Y., Purcell, M., Allan, D., & Tanner, E. (2014). Risk of fracture modelling in healthy and osteoporotic tibia. European Orthopaedic Research Society's 22nd Annual Meeting, Nantes, France.
Gislason, Magnus ; Coupaud, Sylvie ; Sasagawa, Keisuke ; Tanabe, Yuji ; Purcell, Mariel ; Allan, David ; Tanner, Elizabeth. / Risk of fracture modelling in healthy and osteoporotic tibia. European Orthopaedic Research Society's 22nd Annual Meeting, Nantes, France.
@conference{918779b2f00a43469eff1401d5f10279,
title = "Risk of fracture modelling in healthy and osteoporotic tibia",
abstract = "Background: Osteoporosis has detrimental effects on the structural integrity of the skeletal system. The degeneration of bone material can be seen in spinal cord injury (SCI) patients due to lack of mechanical stimulation following extensive muscle paralysis. Aim: The aim of the study was to simulate and compare the mechanical response of the tibia of both able bodied (AB) subjects and SCI patients under four different types of loading, compression, A-P bending, M-L bending and torsion, and map how the risk of fracture (RF) correlated with bone mineral density (BMD) and type of loading. The stress distribution was calculated using the finite element method. Material and methods: Eight male SCI subjects with varying degrees of bone loss and four age-matched AB male control subjects were recruited for the study. The bones were scanned in vivo at the Queen Elizabeth National Spinal Injuries Unit (Glasgow) using a multi-slice peripheral Quantitative Computed Tomography, and individualised three dimensional finite element models were created of each subject's tibia from the scans using Mimics (Materialise). The material assignment was based on the greyscale pixel values of the scan images from which BMD was calculated. An empirical power law was used to estimate the Young’s modulus and the ultimate tensile strength for each element from the BMD. Results: Analysis of the material distribution in the tibia showed large numbers of elements with low BMD for the SCI group whereas a larger proportion of high density bone was seen in the AB group. The average BMD was lower at the epiphyses for the SCI group, but overall little difference was seen around the diaphysis between the groups. The finite element modelling was carried out using Abaqus (v.6.11, Simulia). The analysis was performed using 4 node tetrahedral elements and the element density was on average 2.9 elements/mm^3. The von Mises stress values of each element were exported and the percentage of elements exceeding an estimated RF threshold calculated. The results showed an increase in percentage of elements exceeding the RF threshold for the SCI group in bending and torsion compared to compression which is in agreement with clinical experience. On average, the SCI group showed greater numbers of elements exceeding the RF threshold than the AB group for all four loading conditions. The majority of the fractured elements were located around the epiphyses for all subjects. Conclusions and clinical implications: Using the analysis methods described can provide clinicians with useful patient-specific information about the amount of loading that can be put safely on the osteoporotic bone. That is of importance in SCI patient rehabilitation where musculoskeletal rehabilitation techniques include stimulating the bone via mechanical methods such as vibration or electrical stimulation of the muscles. Using such subject-specific analysis it is possible to determine the location along the bone that is most likely to fail along with the risk associated.",
keywords = "osteoporosis , skeletal system, spinal chord injury, SCI, patient rehabilitation, musculoskeletal rehabilitation",
author = "Magnus Gislason and Sylvie Coupaud and Keisuke Sasagawa and Yuji Tanabe and Mariel Purcell and David Allan and Elizabeth Tanner",
year = "2014",
month = "7",
language = "English",
note = "European Orthopaedic Research Society's 22nd Annual Meeting ; Conference date: 02-07-2014 Through 04-07-2014",

}

Gislason, M, Coupaud, S, Sasagawa, K, Tanabe, Y, Purcell, M, Allan, D & Tanner, E 2014, 'Risk of fracture modelling in healthy and osteoporotic tibia' European Orthopaedic Research Society's 22nd Annual Meeting, Nantes, France, 2/07/14 - 4/07/14, .

Risk of fracture modelling in healthy and osteoporotic tibia. / Gislason, Magnus; Coupaud, Sylvie; Sasagawa, Keisuke; Tanabe, Yuji; Purcell, Mariel; Allan, David; Tanner, Elizabeth.

2014. European Orthopaedic Research Society's 22nd Annual Meeting, Nantes, France.

Research output: Contribution to conferenceProceeding

TY - CONF

T1 - Risk of fracture modelling in healthy and osteoporotic tibia

AU - Gislason, Magnus

AU - Coupaud, Sylvie

AU - Sasagawa, Keisuke

AU - Tanabe, Yuji

AU - Purcell, Mariel

AU - Allan, David

AU - Tanner, Elizabeth

PY - 2014/7

Y1 - 2014/7

N2 - Background: Osteoporosis has detrimental effects on the structural integrity of the skeletal system. The degeneration of bone material can be seen in spinal cord injury (SCI) patients due to lack of mechanical stimulation following extensive muscle paralysis. Aim: The aim of the study was to simulate and compare the mechanical response of the tibia of both able bodied (AB) subjects and SCI patients under four different types of loading, compression, A-P bending, M-L bending and torsion, and map how the risk of fracture (RF) correlated with bone mineral density (BMD) and type of loading. The stress distribution was calculated using the finite element method. Material and methods: Eight male SCI subjects with varying degrees of bone loss and four age-matched AB male control subjects were recruited for the study. The bones were scanned in vivo at the Queen Elizabeth National Spinal Injuries Unit (Glasgow) using a multi-slice peripheral Quantitative Computed Tomography, and individualised three dimensional finite element models were created of each subject's tibia from the scans using Mimics (Materialise). The material assignment was based on the greyscale pixel values of the scan images from which BMD was calculated. An empirical power law was used to estimate the Young’s modulus and the ultimate tensile strength for each element from the BMD. Results: Analysis of the material distribution in the tibia showed large numbers of elements with low BMD for the SCI group whereas a larger proportion of high density bone was seen in the AB group. The average BMD was lower at the epiphyses for the SCI group, but overall little difference was seen around the diaphysis between the groups. The finite element modelling was carried out using Abaqus (v.6.11, Simulia). The analysis was performed using 4 node tetrahedral elements and the element density was on average 2.9 elements/mm^3. The von Mises stress values of each element were exported and the percentage of elements exceeding an estimated RF threshold calculated. The results showed an increase in percentage of elements exceeding the RF threshold for the SCI group in bending and torsion compared to compression which is in agreement with clinical experience. On average, the SCI group showed greater numbers of elements exceeding the RF threshold than the AB group for all four loading conditions. The majority of the fractured elements were located around the epiphyses for all subjects. Conclusions and clinical implications: Using the analysis methods described can provide clinicians with useful patient-specific information about the amount of loading that can be put safely on the osteoporotic bone. That is of importance in SCI patient rehabilitation where musculoskeletal rehabilitation techniques include stimulating the bone via mechanical methods such as vibration or electrical stimulation of the muscles. Using such subject-specific analysis it is possible to determine the location along the bone that is most likely to fail along with the risk associated.

AB - Background: Osteoporosis has detrimental effects on the structural integrity of the skeletal system. The degeneration of bone material can be seen in spinal cord injury (SCI) patients due to lack of mechanical stimulation following extensive muscle paralysis. Aim: The aim of the study was to simulate and compare the mechanical response of the tibia of both able bodied (AB) subjects and SCI patients under four different types of loading, compression, A-P bending, M-L bending and torsion, and map how the risk of fracture (RF) correlated with bone mineral density (BMD) and type of loading. The stress distribution was calculated using the finite element method. Material and methods: Eight male SCI subjects with varying degrees of bone loss and four age-matched AB male control subjects were recruited for the study. The bones were scanned in vivo at the Queen Elizabeth National Spinal Injuries Unit (Glasgow) using a multi-slice peripheral Quantitative Computed Tomography, and individualised three dimensional finite element models were created of each subject's tibia from the scans using Mimics (Materialise). The material assignment was based on the greyscale pixel values of the scan images from which BMD was calculated. An empirical power law was used to estimate the Young’s modulus and the ultimate tensile strength for each element from the BMD. Results: Analysis of the material distribution in the tibia showed large numbers of elements with low BMD for the SCI group whereas a larger proportion of high density bone was seen in the AB group. The average BMD was lower at the epiphyses for the SCI group, but overall little difference was seen around the diaphysis between the groups. The finite element modelling was carried out using Abaqus (v.6.11, Simulia). The analysis was performed using 4 node tetrahedral elements and the element density was on average 2.9 elements/mm^3. The von Mises stress values of each element were exported and the percentage of elements exceeding an estimated RF threshold calculated. The results showed an increase in percentage of elements exceeding the RF threshold for the SCI group in bending and torsion compared to compression which is in agreement with clinical experience. On average, the SCI group showed greater numbers of elements exceeding the RF threshold than the AB group for all four loading conditions. The majority of the fractured elements were located around the epiphyses for all subjects. Conclusions and clinical implications: Using the analysis methods described can provide clinicians with useful patient-specific information about the amount of loading that can be put safely on the osteoporotic bone. That is of importance in SCI patient rehabilitation where musculoskeletal rehabilitation techniques include stimulating the bone via mechanical methods such as vibration or electrical stimulation of the muscles. Using such subject-specific analysis it is possible to determine the location along the bone that is most likely to fail along with the risk associated.

KW - osteoporosis

KW - skeletal system

KW - spinal chord injury

KW - SCI

KW - patient rehabilitation

KW - musculoskeletal rehabilitation

UR - http://mymarketingteam.conference-services.net/resources/1789/3989/pdf/EORS2014_0288.pdf

UR - http://eors2014.org/

UR - http://www.eors.eu/

M3 - Proceeding

ER -

Gislason M, Coupaud S, Sasagawa K, Tanabe Y, Purcell M, Allan D et al. Risk of fracture modelling in healthy and osteoporotic tibia. 2014. European Orthopaedic Research Society's 22nd Annual Meeting, Nantes, France.