27.12.2012 Views

ARUP; ISBN: 978-0-9562121-5-3 - CMBBE 2012 - Cardiff University

ARUP; ISBN: 978-0-9562121-5-3 - CMBBE 2012 - Cardiff University

ARUP; ISBN: 978-0-9562121-5-3 - CMBBE 2012 - Cardiff University

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

osteoporotic (OP). The sample volume fractions (BV/TV) range from 14.6 to 51.1%.<br />

The scans are 170x170x170 voxels at a resolution of 30µm.<br />

Figure 1: Locations of 5x5x5mm cubic samples within the femoral head. The<br />

nomenclature is as follows: MP – medial posterior, MSC – medial superior central,<br />

MC – medial central, MIC – medial inferior central, MA – medial anterior, CP –<br />

central posterior, SC – superior central, CC – central central, IC – inferior central,<br />

CA – central anterior.<br />

4.2 Procedure<br />

The 5x5x5mm volumes of interest (VOIs) were directly considered as representative<br />

volume elements as supported by an examination of the continuum assumption for<br />

trabecular bone [23, 24, 25].<br />

The µ-CT scans were first binarised into solid and void using an identical procedure for<br />

all samples comprising contrast enhancement and automatic thresholding based on an<br />

analysis of the image. A morphological analysis was then performed to determine<br />

indices of micro-architecture including bone volume fraction (BV/TV) and trabecular<br />

thickness (Tb.Th).<br />

In order to perform mechanical tests on the samples, a finite element mesh was<br />

generated from the binarised images using 10-noded tetrahedral elements. The tissuelevel<br />

elastic properties were E = 18GPa and v = 0.3, which were applied to all elements.<br />

Three uniaxial and three pure shear strain tests were applied independently to each<br />

sample [26, 27]. The resulting stresses were averaged volumetrically, and from this the<br />

full anisotropic elastic stiffness tensor was computed as<br />

E= 1<br />

V ∫ V<br />

E t M dV (1)<br />

where V is the sample volume, Et is the tissue stiffness tensor, and M is the local<br />

structure tensor. The stiffness tensor was then optimised to find the most orthotropic<br />

orientation, and the error in the orthotropy assumption was calculated [26]. The direct<br />

stiffness (E1111) was also maximised separately in order to find its largest value and<br />

direction.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!