Variation of the mechanical properties of PMMA to suit osteoporotic cancellous bone

Poly(methyl methacrylate) (PMMA) is by far the most frequently used bone substitute material for vertebroplasty. However, there are serious complications, such as cement leakage and an increased fracture rate of the adjacent vertebral bodies. The latter may be related to the mechanical properties of the augmented segment within the osteoporotic spine. A possible counter-measure is prophylactic augmentation at additional levels, but this aggravates the risk for the patient. Introduction of pores is a possible method to reduce the inherent high stiffness of PMMA. This study investigates the effect of porosity on the mechanical properties of PMMA bone cement. Different fractions of a highly viscous liquid were mixed into the PMMA during preparation. An open-porous material with adjustable mechanical properties resulted after removal of the aqueous phase. Different radiopacifiers were admixed to investigate their suitability for vertebroplasty. The final material was characterized mechanically by compressive testing, microscopically and radiologically. In addition, the monomer release subsequent to hardening was measured by means of gas chromatography. The Young's modulus in compression could be varied between 2800 +/- 70 MPa and 120 +/- 150 MPa, and the compression ultimate strength between 170 +/- 5 MPa and 8 +/- 9 MPa for aqueous fractions ranging between 0 and 50% of volume. Only a slight decrease of the Young's modulus and small changes of ultimate strength were found when the mixing time was increased. An organic hydrophilic and lipophilic radiopacifier led to a higher Young's modulus of the porous material; however, the ultimate strength was not significantly affected by adding different radiopacifiers to the porous cement. The radiopacity was lost after washing the aqueous phase out of the pores. No separation occurred between the aqueous and the PMMA phase during injection into an open porous ceramic material. The monomer released was found to increase for increasing aqueous fractions, but remained comparable in magnitude to standard PMMA. This study demonstrates that a conventional PMMA can be modified to obtain a range of mechanical properties, including those of osteoporotic bone.

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.

Augmentation of mechanical properties in osteoporotic vertebral bones--a biomechanical investigation of vertebroplasty efficacy with different bone cements

Recent clinical trials have reported favorable early results for transpedicular vertebral cement reinforcement of osteoporotic vertebral insufficiencies. There is, however, a lack of basic data on the application, safety and biomechanical efficacy of materials such as polymethyl-methacrylate (PMMA) and calciumphospate (CaP) cements. The present study analyzed 33 vertebral pairs from five human cadaver spines. Thirty-nine vertebrae were osteoporotic (bone mineral density < 0.75 g/cm2), 27 showed nearly normal values. The cranial vertebra of each pair was augmented with either PMMA (Palacos E-Flow) or experimental brushite cement (EBC), with the caudal vertebra as a control. PMMA and EBC were easy to inject, and vertebral fillings of 20-50% were achieved. The maximal possible filling was inversely correlated to the bone mineral density (BMD) values. Cement extrusion into the spinal canal was observed in 12% of cases. All specimens were subjected to axial compression tests in a displacement-controlled mode. From load-displacement curves, the stiffness, S, and the maximal force before failure, Fmax, were determined. Compared with the native control vertebrae, a statistically significant increase in vertebral stiffness and Fmax was observed by the augmentation. With PMMA the stiffness increased by 174% (P = 0.018) and Fmax by 195% (P = 0.001); the corresponding augmentation with EBC was 120% (P = 0.03) and 113% (P = 0.002). The lower the initial BMD, the more pronounced was the augmentation effect. Both PMMA and EBC augmentation reliably and significantly raised the stiffness and maximal tolerable force until failure in osteoporotic vertebral bodies. In non-porotic specimens, no significant increase was achieved.

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.

In situ micropillar compression reveals superior strength and ductility but an absence of damage in lamellar bone

Ageing societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the hierarchical organization of bone. Here, we investigate the mechanical response under monotonic and cyclic compression of both single osteonal lamellae and macroscopic samples containing numerous osteons. Micropillar compression tests in a scanning electron microscope, microindentation and macroscopic compression tests were performed on dry ovine bone to identify the elastic modulus, yield stress, plastic deformation, damage accumulation and failure mechanisms. We found that isolated lamellae exhibit a plastic behaviour, with higher yield stress and ductility but no damage. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behaviour of bone at the microscale to a quasi-brittle response driven by the growth of cracks along interfaces or in the vicinity of pores at the macroscale.

Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They 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 were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

Biomechanical properties of the atlantoaxial joint with naturally-occurring instability in a toy breed dog. A comparative descriptive cadaveric study

The biomechanical properties of the atlanto-axial joint in a young Yorkshire Terrier dog with spontaneous atlantoaxial instability were compared to those of another young toy breed dog with a healthy atlantoaxial joint. The range-of-motion was increased in flexion and lateral bending in the unstable joint. In addition, lateral bending led to torsion and dorsal dislocation of the axis within the atlas. On gross examination, the dens ligaments were absent and a longitudinal tear of the tectorial membrane was observed. These findings suggest that both ventral and lateral flexion may lead to severe spinal cord compression, and that the tectorial membrane may play a protective role in some cases of atlantoaxial instability.

Properties of granular analogue model materials: A community wide survey

We report the material properties of 26 granular analogue materials used in 14 analogue modelling laboratories. We determined physical characteristics such as bulk density, grain size distribution, and grain shape, and performed ring shear tests to determine friction angles and cohesion, and uniaxial compression tests to evaluate the compaction behaviour. Mean grain size of the materials varied between c. 100 and 400 μm. Analysis of grain shape factors shows that the four different classes of granular materials (14 quartz sands, 5 dyed quartz sands, 4 heavy mineral sands and 3 size fractions of glass beads) can be broadly divided into two groups consisting of 12 angular and 14 rounded materials. Grain shape has an influence on friction angles, with most angular materials having higher internal friction angles (between c. 35° and 40°) than rounded materials, whereas well-rounded glass beads have the lowest internal friction angles (between c. 25° and 30°). We interpret this as an effect of intergranular sliding versus rolling. Most angular materials have also higher basal friction angles (tested for a specific foil) than more rounded materials, suggesting that angular grains scratch and wear the foil. Most materials have an internal cohesion in the order of 20–100 Pa except for well-rounded glass beads, which show a trend towards a quasi-cohesionless (C < 20 Pa) Coulomb-type material. The uniaxial confined compression tests reveal that rounded grains generally show less compaction than angular grains. We interpret this to be related to the initial packing density after sifting, which is higher for rounded grains than for angular grains. Ring-shear test data show that angular grains undergo a longer strain-hardening phase than more rounded materials. This might explain why analogue models consisting of angular grains accommodate deformation in a more distributed manner prior to strain localisation than models consisting of rounded grains.