Second-Generation Orthopaedic Measurement Device for in-vivo Assessment of Bone Quality

There are two main types of bone in the human body, trabecular and cortical bone. Cortical bone is primarily found on the outer surface of most bones in the body while trabecular bone is found in vertebrae and at the end of long bones (Ross 2007). Osteoporosis is a condition that compromises the structural integrity of trabecular bone, greatly reducing the ability of the bone to absorb energy from falls. The current method for diagnosing osteoporosis and predicting fracture risk is measurement of bone mineral density. Limitations of this method include dependence on the bone density measurement device and dependence on type of test and measurement location (Rubin 2005). Each year there are approximately 250,000 hip fractures in the United States due to osteoporosis (Kleerekoper 2006). Currently, the most common method for repairing a hip fracture is a hip fixation surgery. During surgery, a temporary guide wire is inserted to guide the permanent screw into place and then removed. It is believed that directly measuring this screw pullout force may result in a better assessment of bone quality than current indirect measurement techniques (T. Bowen 2008-2010, pers. comm.). The objective of this project is to design a device that can measure the force required to extract this guide wire. It is believed that this would give the surgeon a direct, quantitative measurement of bone quality at the site of the fixation. A first generation device was designed by a Bucknell Biomedical Engineering Senior Design team during the 2008- 2009 Academic Year. The first step of this project was to examine the device, conduct a thorough design analysis, and brainstorm new concepts. The concept selected uses a translational screw to extract the guide wire. The device was fabricated and underwent validation testing to ensure that the device was functional and met the required engineering specifications. Two tests were conducted, one to test the functionality of the device by testing if the device gave repeatable results, and the other to test the sensitivity of the device to misalignment. Guide wires were extracted from 3 materials, low density polyethylene, ultra high molecular weight polyethylene, and polypropylene and the force of extraction was measured. During testing, it was discovered that the spring in the device did not have a high enough spring constant to reach the high forces necessary for extracting the wires without excessive deflection of the spring. The test procedure was modified slightly so the wires were not fully threaded into the material. The testing results indicate that there is significant variation in the screw pullout force, up to 30% of the average value. This significant variation was attributed to problems in the testing and data collection, and a revised set of tests was proposed to better evaluate the performance of the device. The fabricated device is a fully-functioning prototype and further refinements and testing of the device may lead to a 3rd generation version capable of measuring the screw pullout force during hip fixation surgery.

A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads

Due to the inherent limitations of DXA, assessment of the biomechanical properties of vertebral bodies relies increasingly on CT-based finite element (FE) models, but these often use simplistic material behaviour and/or single loading cases. In this study, we applied a novel constitutive law for bone elasticity, plasticity and damage to FE models created from coarsened pQCT images of human vertebrae, and compared vertebral stiffness, strength and damage accumulation for axial compression, anterior flexion and a combination of these two cases. FE axial stiffness and strength correlated with experiments and were linearly related to flexion properties. In all loading modes, damage localised preferentially in the trabecular compartment. Damage for the combined loading was higher than cumulated damage produced by individual compression and flexion. In conclusion, this FE method predicts stiffness and strength of vertebral bodies from CT images with clinical resolution and provides insight into damage accumulation in various loading modes.

Biomechanics of meniscal horn attachments

Menisci are anchored to the tibia by means of ligament-like structures called meniscal attachments. Failure material properties of bovine meniscal attachments were obtained. There were no significant differences in the structural properties or ultimate stress between the meniscal attachments (p>0.05). Furthermore, Glycosaminoglycan (GAG) fraction and crimping frequency was obtained for each attachment using histology and differential interference contrast (DIC) respectively. Results showed that the anterior attachment’s insertion had the greatest GAG fraction when compared to the posterior attachment’s insertion. Crimp frequency of the collagen fibrils was homogeneous along the length. Moreover, Scanning Electron Microscopy (SEM) technique was used to reveal the morphology of collagen in human meniscal attachments. Its midsubstance was composed of collagen fascicles running parallel to the longitudinal axis, with a few fibrils running obliquely, and others transversely. There were no differences between attachments for crimping angle or length. Since ligamentous-type tissues are comprised mainly of water, the fluid pressure within meniscal horn attachments was measured using a Fiber Optic Microsensor (FOM). Four cadaveric human joints were subjected to 2BW compressive load (ramp) at 0-, 15-, and 30-degrees of flexion for a minute and then the load was hold for 20 minutes (equilibrium). There were significant differences between 0- and 15- (p1– c5) were obtained. Significant differences were found on the straightened collagen fibers coefficient (c5) between MP and LA attachments (p

The influence of yield surface shape and damage in the depth-dependent response of bone tissue to nanoindentation using spherical and Berkovich indenters

Prevention and treatment of osteoporosis rely on understanding of the micromechanical behaviour of bone and its influence on fracture toughness and cell-mediated adaptation processes. Postyield properties may be assessed by nonlinear finite element simulations of nanoindentation using elastoplastic and damage models. This computational study aims at determining the influence of yield surface shape and damage on the depth-dependent response of bone to nanoindentation using spherical and conical tips. Yield surface shape and damage were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic-to-total work ratio is well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not statistically significant (p<0.0001). For spherical tips, damage was not a significant parameter (p<0.0001). The gained knowledge can be used for developing an inverse method for identification of postelastic properties of bone from nanoindentation.

Embedding of human vertebral bodies leads to higher ultimate load and altered damage localisation under axial compression

Computer tomography (CT)-based finite element (FE) models of vertebral bodies assess fracture load in vitro better than dual energy X-ray absorptiometry, but boundary conditions affect stress distribution under the endplates that may influence ultimate load and damage localisation under post-yield strains. Therefore, HRpQCT-based homogenised FE models of 12 vertebral bodies were subjected to axial compression with two distinct boundary conditions: embedding in polymethylmethalcrylate (PMMA) and bonding to a healthy intervertebral disc (IVD) with distinct hyperelastic properties for nucleus and annulus. Bone volume fraction and fabric assessed from HRpQCT data were used to determine the elastic, plastic and damage behaviour of bone. Ultimate forces obtained with PMMA were 22% higher than with IVD but correlated highly (R2 = 0.99). At ultimate force, distinct fractions of damage were computed in the endplates (PMMA: 6%, IVD: 70%), cortex and trabecular sub-regions, which confirms previous observations that in contrast to PMMA embedding, failure initiated underneath the nuclei in healthy IVDs. In conclusion, axial loading of vertebral bodies via PMMA embedding versus healthy IVD overestimates ultimate load and leads to distinct damage localisation and failure pattern.

Electrospinning Auricular Shaped Scaffolds for Tissue Engineering

Poly(ɛ)caprolactone scaffolds have been electrospun directly into an auricular shaped conductive mould. Bovine chondrocytes were harvested from articular cartilage and seeded onto 16 of the produced scaffolds, which received either an ethanol (group A) or a plasma treatment (group B) for sterilisation before seeding. The seeded scaffolds were cultured for 3 weeks in vitro and analysed with regard to total DNA and GAG content as well as the expression of AGG, COL1, COL2, MMP3 and MMP13. Rapid cell proliferation and GAG accumulation was observed until week 2. However, total DNA and GAG content decreased again in week 3. qPCR data shows a slight increase in the expression of anabolic genes and a slight decrease for the catabolic genes, with a significant difference between the groups A and B only for COL2 and MMP13.

NEW TOOLS FOR MONITORING GAMMA CAMERA UNIFORMITY

Detector uniformity is a fundamental performance characteristic of all modern gamma camera systems, and ensuring a stable, uniform detector response is critical for maintaining clinical images that are free of artifact. For these reasons, the assessment of detector uniformity is one of the most common activities associated with a successful clinical quality assurance program in gamma camera imaging. The evaluation of this parameter, however, is often unclear because it is highly dependent upon acquisition conditions, reviewer expertise, and the application of somewhat arbitrary limits that do not characterize the spatial location of the non-uniformities. Furthermore, as the goal of any robust quality control program is the determination of significant deviations from standard or baseline conditions, clinicians and vendors often neglect the temporal nature of detector degradation (1). This thesis describes the development and testing of new methods for monitoring detector uniformity. These techniques provide more quantitative, sensitive, and specific feedback to the reviewer so that he or she may be better equipped to identify performance degradation prior to its manifestation in clinical images. The methods exploit the temporal nature of detector degradation and spatially segment distinct regions-of-non-uniformity using multi-resolution decomposition. These techniques were tested on synthetic phantom data using different degradation functions, as well as on experimentally acquired time series floods with induced, progressively worsening defects present within the field-of-view. The sensitivity of conventional, global figures-of-merit for detecting changes in uniformity was evaluated and compared to these new image-space techniques. The image-space algorithms provide a reproducible means of detecting regions-of-non-uniformity prior to any single flood image’s having a NEMA uniformity value in excess of 5%. The sensitivity of these image-space algorithms was found to depend on the size and magnitude of the non-uniformities, as well as on the nature of the cause of the non-uniform region. A trend analysis of the conventional figures-of-merit demonstrated their sensitivity to shifts in detector uniformity. The image-space algorithms are computationally efficient. Therefore, the image-space algorithms should be used concomitantly with the trending of the global figures-of-merit in order to provide the reviewer with a richer assessment of gamma camera detector uniformity characteristics.

IN-VIVO CT DOSIMETRY DURING VIRTUAL COLONOSCOPY

Virtual colonoscopy (VC) is a minimally invasive means for identifying colorectal polyps and colorectal lesions by insufflating a patient’s bowel, applying contrast agent via rectal catheter, and performing multi-detector computed tomography (MDCT) scans. The technique is recommended for colonic health screening by the American Cancer Society but not funded by the Centers for Medicare and Medicaid Services (CMS) partially because of potential risks from radiation exposure. To date, no in‐vivo organ dose measurements have been performed for MDCT scans; thus, the accuracy of any current dose estimates is currently unknown. In this study, two TLDs were affixed to the inner lumen of standard rectal catheters used in VC, and in-vivo rectal dose measurements were obtained within 6 VC patients. In order to calculate rectal dose, TLD-100 powder response was characterized at diagnostic doses such that appropriate correction factors could be determined for VC. A third-order polynomial regression with a goodness of fit factor of R2=0.992 was constructed from this data. Rectal dose measurements were acquired with TLDs during simulated VC within a modified anthropomorphic phantom configured to represent three sizes of patients undergoing VC. The measured rectal doses decreased in an exponential manner with increasing phantom effective diameter, with R2=0.993 for the exponential regression model and a maximum percent coefficient of variation (%CoV) of 4.33%. In-vivo measurements yielded rectal doses ranged from that decreased exponentially with increasing patient effective diameter, in a manner that was also favorably predicted by the size specific dose estimate (SSDE) model for all VC patients that were of similar age, body composition, and TLD placement. The measured rectal dose within a younger patient was favorably predicted by the anthropomorphic phantom dose regression model due to similarities in the percentages of highly attenuating material at the respective measurement locations and in the placement of the TLDs. The in-vivo TLD response did not increase in %CoV with decreasing dose, and the largest %CoV was 10.0%.

Image-based biomechanical assessment of vertebral body and intervertebral disc in the human lumbar spine

Life expectancy continuously increases but our society faces age-related conditions. Among musculoskeletal diseases, osteoporosis associated with risk of vertebral fracture and degenerative intervertebral disc (IVD) are painful pathologies responsible for tremendous healthcare costs. Hence, reliable diagnostic tools are necessary to plan a treatment or follow up its efficacy. Yet, radiographic and MRI techniques, respectively clinical standards for evaluation of bone strength and IVD degeneration, are unspecific and not objective. Increasingly used in biomedical engineering, CT-based finite element (FE) models constitute the state-of-art for vertebral strength prediction. However, as non-invasive biomechanical evaluation and personalised FE models of the IVD are not available, rigid boundary conditions (BCs) are applied on the FE models to avoid uncertainties of disc degeneration that might bias the predictions. Moreover, considering the impact of low back pain, the biomechanical status of the IVD is needed as a criterion for early disc degeneration. Thus, the first FE study focuses on two rigid BCs applied on the vertebral bodies during compression test of cadaver vertebral bodies, vertebral sections and PMMA embedding. The second FE study highlights the large influence of the intervertebral disc’s compliance on the vertebral strength, damage distribution and its initiation. The third study introduces a new protocol for normalisation of the IVD stiffness in compression, torsion and bending using MRI-based data to account for its morphology. In the last study, a new criterion (Otsu threshold) for disc degeneration based on quantitative MRI data (axial T2 map) is proposed. The results show that vertebral strength and damage distribution computed with rigid BCs are identical. Yet, large discrepancies in strength and damage localisation were observed when the vertebral bodies were loaded via IVDs. The normalisation protocol attenuated the effect of geometry on the IVD stiffnesses without complete suppression. Finally, the Otsu threshold computed in the posterior part of annulus fibrosus was related to the disc biomechanics and meet objectivity and simplicity required for a clinical application. In conclusion, the stiffness normalisation protocol necessary for consistent IVD comparisons and the relation found between degeneration, mechanical response of the IVD and Otsu threshold lead the way for non-invasive evaluation biomechanical status of the IVD. As the FE prediction of vertebral strength is largely influenced by the IVD conditions, this data could also improve the future FE models of osteoporotic vertebra.

Nitinol Stent Oversizing in Patients Undergoing Popliteal Artery Revascularization: A Finite Element Study

Nitinol stent oversizing is frequently performed in peripheral arteries to ensure a desirable lumen gain. However, the clinical effect of mis-sizing remains controversial. The goal of this study was to provide a better understanding of the structural and hemodynamic effects of Nitinol stent oversizing. Five patient-specific numerical models of non-calcified popliteal arteries were developed to simulate the deployment of Nitinol stents with oversizing ratios ranging from 1.1 to 1.8. In addition to arterial biomechanics, computational fluid dynamics methods were adopted to simulate the physiological blood flow inside the stented arteries. Results showed that stent oversizing led to a limited increase in the acute lumen gain, albeit at the cost of a significant increase in arterial wall stresses. Furthermore, localized areas affected by low Wall Shear Stress increased with higher oversizing ratios. Stents were also negatively impacted by the procedure as their fatigue safety factors gradually decreased with oversizing. These adverse effects to both the artery walls and stents may create circumstances for restenosis. Although the ideal oversizing ratio is stent-specific, this study showed that Nitinol stent oversizing has a very small impact on the immediate lumen gain, which contradicts the clinical motivations of the procedure.