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.

Intraoperative Mechanical Measurement of Bone Quality With the DensiProbe.

Reduced bone stock can result in fractures that mostly occur in the spine, distal radius, and proximal femur. In case of operative treatment, osteoporosis is associated with an increased failure rate. To estimate implant anchorage, mechanical methods seem to be promising to measure bone strength intraoperatively. It has been shown that the mechanical peak torque correlates with the local bone mineral density and screw failure load in hip, hindfoot, humerus, and spine in vitro. One device to measure mechanical peak torque is the DensiProbe (AO Research Institute, Davos, Switzerland). The device has shown its effectiveness in mechanical peak torque measurement in mechanical testing setups for the use in hip, hindfoot, and spine. In all studies, the correlation of mechanical torque measurement and local bone mineral density and screw failure load could be shown. It allows the surgeon to judge local bone strength intraoperatively directly at the region of interest and gives valuable information if additional augmentation is needed. We summarize methods of this new technique, its advantages and limitations, and give an overview of actual and possible future applications.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p ≤ 0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Continuum damage interactions between tension and compression in osteonal bone

Skeletal diseases such as osteoporosis impose a severe socio-economic burden to ageing societies. Decreasing mechanical competence causes a rise in bone fracture incidence and mortality especially after the age of 65 y. The mechanisms of how bone damage is accumulated under different loading modes and its impact on bone strength are unclear. We hypothesise that damage accumulated in one loading mode increases the fracture risk in another. This study aimed at identifying continuum damage interactions between tensile and compressive loading modes. We propose and identify the material constants of a novel piecewise 1D constitutive model capable of describing the mechanical response of bone in combined tensile and compressive loading histories. We performed several sets of loading–reloading experiments to compute stiffness, plastic strains, and stress-strain curves. For tensile overloading, a stiffness reduction (damage) of 60% at 0.65% accumulated plastic strain was detectable as stiffness reduction of 20% under compression. For compressive overloading, 60% damage at 0.75% plastic strain was detectable as a stiffness reduction of 50% in tension. Plastic strain at ultimate stress was the same in tension and compression. Compression showed softening and tension exponential hardening in the post-yield regime. The hardening behaviour in compression is unaffected by a previous overload in tension but the hardening behaviour in tension is affected by a previous overload in compression as tensile reloading strength is significantly reduced. This paper demonstrates how damage accumulated under one loading mode affects the mechanical behaviour in another loading mode. To explain this and to illustrate a possible implementation we proposed a theoretical model. Including such loading mode dependent damage and plasticity behaviour in finite element models will help to improve fracture risk analysis of whole bones and bone implant structures.

Biomechanical Role of Bone Anisotropy Estimated on Clinical CT Scans by Image Registration

Image-based modeling is a popular approach to perform patient-specific biomechanical simulations. Accurate modeling is critical for orthopedic application to evaluate implant design and surgical planning. It has been shown that bone strength can be estimated from the bone mineral density (BMD) and trabecular bone architecture. However, these findings cannot be directly and fully transferred to patient-specific modeling since only BMD can be derived from clinical CT. Therefore, the objective of this study was to propose a method to predict the trabecular bone structure using a µCT atlas and an image registration technique. The approach has been evaluated on femurs and patellae under physiological loading. The displacement and ultimate force for femurs loaded in stance position were predicted with an error of 2.5% and 3.7%, respectively, while predictions obtained with an isotropic material resulted in errors of 7.3% and 6.9%. Similar results were obtained for the patella, where the strain predicted using the registration approach resulted in an improved mean squared error compared to the isotropic model. We conclude that the registration of anisotropic information from of a single template bone enables more accurate patient-specific simulations from clinical image datasets than isotropic model.

The effects of long (C20/22) and short (C18) chain omega-3 fatty acids on keel bone fractures, bone biomechanics, behavior, and egg production in free-range laying hens.

Keel fractures in the laying hen are the most critical animal welfare issue facing the egg production industry, particularly with the increased use of extensive systems in response to the 2012 EU directive banning conventional battery cages. The current study is aimed at assessing the effects of 2 omega-3 (n3) enhanced diets on bone health, production endpoints, and behavior in free-range laying hens. Data was collected from 2 experiments over 2 laying cycles, each of which compared a (n3) supplemented diet with a control diet. Experiment 1 employed a diet supplemented with a 60:40 fish oil-linseed mixture (n3:n6 to 1.35) compared with a control diet (n3:n6 to 0.11), whereas the n3 diet in Experiment 2 was supplemented with a 40:60 fish oil-linseed (n3:n6 to 0.77) compared to the control diet (n3:n6 to 0.11). The n3 enhanced diet of Experiment 1 had a higher n3:n6 ratio, and a greater proportion of n3 in the long chain (C20/22) form (0.41 LC:SC) than that of Experiment 2 (0.12 LC:SC). Although dietary treatment was successful in reducing the frequency of fractures by approximately 27% in Experiment 2, data from Experiment 1 indicated the diet actually induced a greater likelihood of fracture (odds ratio: 1.2) and had substantial production detriment. Reduced keel breakage during Experiment 2 could be related to changes in bone health as n3-supplemented birds demonstrated greater load at failure of the keel, and tibiae and humeri that were more flexible. These results support previous findings that n3-supplemented diets can reduce fracture likely by increasing bone strength, and that this can be achieved without detriment to production. However, our findings suggest diets with excessive quantities of n3, or very high levels of C20/22, may experience health and production detriments. Further research is needed to optimize the quantity and type of n3 in terms of bone health and production variables and investigate the potential associated mechanisms.

Strength indices of the proximal femur and shaft in prepubertal female gymnasts

Introduction/Purpose: The role of impact loading activity on bone mass is well established; however, there are little data on the effects of exercise on bone geometry and indices of bone strength. The primary purpose of this study was to compare indices of bone strength at the proximal femur (PF) between elite premenarcheal gymnasts (N = 30) and age-matched controls (N = 30). Methods: Structural properties of the proximal femur were derived from the hip analyses program and included measurement of subperiosteal width, endosteal diameter, cross-sectional area, bone mineral density, cross-section moment of inertia (CSMI), and section modulus (Z). These parameters were measured for two regions of the PF: the narrow neck (NN), and the shaft (S). In addition, a strength index (S-SI) was calculated at the shaft by dividing the Z at the shaft by the femur length. A secondary purpose was to compare bone mineral content (BMC) values at the total body, lumbar spine, and three sites at the PF (neck, trochanter, and total) between the groups. All dependent values were compared adjusting for height and weight using an ANCOVA procedure and for relative lean body mass post hoc. Results: The gymnasts had significantly greater size-adjusted strength indices (CSMI, Z, and SI) at the NN and S. Gymnasts also had significantly greater size-adjusted BMC at all sites investigated. However, these differences disappeared when adjusted for relative lean body mass. Conclusion: When adjusted for body size, gymnasts had significantly greater indices of both axial strength and bending strength at the NN region of the PF and S, as well as a greater bone SI at the femoral shaft. These differences may be related to greater relative lean body mass attained in gymnastics training.

Effect of aestivation on long bone mechanical properties in the green-striped burrowing frog, Cyclorana alboguttata

The green-striped burrowing frog, Cyclorana alboguttata, survives extended drought periods by burrowing underground and aestivating. These frogs remain immobile within cocoons of shed skin and mucus during aestivation and emerge from their burrows upon heavy rains to feed and reproduce. Extended periods of immobilisation in mammals typically result in bone remodelling and a decrease in bone strength. We examined the effect of aestivation and, hence, prolonged immobilisation on cross-sectional area, histology and bending strength in the femur and tibiolibula of C alboguttata. Frogs were aestivated in soil for three and nine months and were compared with control animals that remained active, were fed and had a continual supply of water. Compared with the controls, long bone size, anatomy and bending strength remained unchanged, indicating an absence of disuse osteoporosis. This preservation of bone tissue properties enables C. alboguttata to compress the active portions of their life history into unpredictable windows of opportunity, whenever heavy rains occur.

The 'muscle-bone unit' during the pubertal growth spurt

Mechanostat theory postulates that developmental changes in bone strength are secondary to the increasing loads imposed by larger muscle forces. Therefore, the increase in muscle strength should precede the increase in bone strength. We tested this prediction using densitometric surrogate measures of muscle force (lean body mass, LBM) and bone strength (bone mineral content, BMC) in a study on 70 boys and 68 girls who were longitudinally examined during pubertal development. On the level of the total body, the peak in LBM accrual preceded the peak in BMC accretion by an average of 0.51 years in girls and by 0.36 years in boys. In the arms, the maximal increase in LBM was followed by arm peak BMC accrual after an interval of 0.71 years in girls and 0.63 years in boys. In the lower extremities, the maximal increase in LBM was followed by peak BMC accrual after an interval of 0.22 years in girls and 0.48 years in boys. A multiple regression model revealed that total body peak LBM velocity, but not peak height velocity and sex, was independently associated with total body peak BMC velocity (r(2) = 0.50; P < 0.001). Similarly, arm and leg peak LBM velocity, but not peak height velocity and sex, were independently associated with arm and leg peak BMC velocity, respectively (r(2) = 0.61 for arms, r(2) = 0.41 for legs; P < 0.001 in both cases). These results are compatible with the view that bone development is driven by muscle development, although the data do not exclude the hypothesis that the two processes are independently determined by genetic mechanisms. (C) 2004 Elsevier Inc. All rights reserved.

Physical activity and strength of the femoral neck during the adolescent growth spurt: A longitudinal analysis

Loading of the femoral neck (FN) is dominated by bending and compressive stresses. We hypothesize that adaptation of the FN to physical activity would be manifested in the cross-sectional area (CSA) and section modulus (Z) of bone, indices of axial and bending strength, respectively. We investigated the influence of physical activity on bone strength during adolescence using 7 years of longitudinal data from 109 boys and 121 girls from the Saskatchewan Paediatric Bone and Mineral Accrual Study (PBMAS). Physical activity data (PAC-Q physical activity inventory) and anthropometric measurements were taken every 6 months and DXA bone scans were measured annually (Hologic QDR2000, array mode). We applied hip structural analysis to derive strength and geometric indices of the femoral neck using DXA scans. To control for maturation, we determined a biological maturity age defined as years from age at peak height velocity (APHV). To account for the repeated measures within individual nature of longitudinal data, multilevel random effects regression analyses were used to analyze the data. When biological maturity age and body size (height and weight) were controlled, in both boys and girls, physical activity was a significant positive independent predictor of CSA and Z of the narrow region of the femoral neck (P < 0.05). There was no independent effect of physical activity on the subperiosteal width of the femoral neck. When leg length and leg lean mass were introduced into the random effects models to control for size and muscle mass of the leg (instead of height and weight), all significant effects of physical activity disappeared. Even among adolescents engaged in normal levels of physical activity, the statistically significant relationship between physical activity and indices of bone strength demonstrate that modifiable lifestyle factors like exercise play an important role in optimizing bone strength during the growing years. Physical activity differences were explained by the interdependence between activity and lean mass considerations. Physical activity is important for optimal development of bone strength. (c) 2005 Elsevier Inc. All rights reserved.

Stable Incretin Mimetics Counter Rapid Deterioration of Bone Quality in Type 1 Diabetes Mellitus.

International audience

Comparison of 3D finite element analysis derived stiffness and BMD to determine the failure load of the excised proximal femur

Introduction: Bone mineral density (BMD) is currently the preferred surrogate for bone strength in clinical practice. Finite element analysis (FEA) is a computer simulation technique that can predict the deformation of a structure when a load is applied, providing a measure of stiffness (Nmm−1). Finite element analysis of X-ray images (3D-FEXI) is a FEA technique whose analysis is derived froma single 2D radiographic image. Methods: 18 excised human femora had previously been quantitative computed tomography scanned, from which 2D BMD-equivalent radiographic images were derived, and mechanically tested to failure in a stance-loading configuration. A 3D proximal femur shape was generated from each 2D radiographic image and used to construct 3D-FEA models. Results: The coefficient of determination (R2%) to predict failure load was 54.5% for BMD and 80.4% for 3D-FEXI. Conclusions: This ex vivo study demonstrates that 3D-FEXI derived from a conventional 2D radiographic image has the potential to significantly increase the accuracy of failure load assessment of the proximal femur compared with that currently achieved with BMD. This approach may be readily extended to routine clinical BMD images derived by dual energy X-ray absorptiometry. Crown Copyright © 2009 Published by Elsevier Ltd on behalf of IPEM. All rights reserved

The role of periosteum : a neglected aspect of osteoporosis

Most current studies on the pathogenesis of osteoporosis emphasize the bone metabolic activities occurring on endosteal surfaces, whereas the periosteal aspect is somewhat neglected. In terms of bone physiology, periosteum plays a determining role in de novo cortical bone formation and cortical bone expansion through periosteum is the most efficient way of increasing bone strength against fractures. Despite the important role of periosteum in the pathogenesis and treatment of osteoporosis, little is known about the structural and cellular features of periosteum in osteoporosis. This chapter will focus on the major changes occurring in the periosteum of osteoporosis and possible implications of these changes in the pathogenesis of osteoporosis. The changes identified in the periosteum of osteoporosis are mainly located in the metaphyseal compartment, which include: (a) much thicker and more cellular cambial layer; (b) increased number of TRAP (tartrate resistant acid phosphatase), VEGF (vascular endothelial growth factor) cells and the degree of vascularization; and (c) enhanced expression of sympathetic nerve fibers. The structural and cellular changes of osteoporotic periosteum indicate that periosteum plays an important role in the cortical bone resorption in metaphyseal areas and this pathological process may be regulated by the sympathetic nervous system.

Comparison of 3D finite element analysis derived stiffness and BMD to determine the failure load of the excised proximal femur

Bone mineral density (BMD) is currently the preferred surrogate for bone strength in clinical practice. Finite element analysis (FEA) is a computer simulation technique that can predict the deformation of a structure when a load is applied, providing a measure of stiffness (N mm− 1). Finite element analysis of X-ray images (3D-FEXI) is a FEA technique whose analysis is derived from a single 2D radiographic image. This ex-vivo study demonstrates that 3D-FEXI derived from a conventional 2D radiographic image has the potential to significantly increase the accuracy of failure load assessment of the proximal femur compared with that currently achieved with BMD.

An experimental and finite element investigation of the biomechanics of vertebral compression fractures

Osteoporosis is a disease characterized by low bone mass and micro-architectural deterioration of bone tissue, with a consequent increase in bone fragility and susceptibility to fracture. Osteoporosis affects over 200 million people worldwide, with an estimated 1.5 million fractures annually in the United States alone, and with attendant costs exceeding $10 billion dollars per annum. Osteoporosis reduces bone density through a series of structural changes to the honeycomb-like trabecular bone structure (micro-structure). The reduced bone density, coupled with the microstructural changes, results in significant loss of bone strength and increased fracture risk. Vertebral compression fractures are the most common type of osteoporotic fracture and are associated with pain, increased thoracic curvature, reduced mobility, and difficulty with self care. Surgical interventions, such as kyphoplasty or vertebroplasty, are used to treat osteoporotic vertebral fractures by restoring vertebral stability and alleviating pain. These minimally invasive procedures involve injecting bone cement into the fractured vertebrae. The techniques are still relatively new and while initial results are promising, with the procedures relieving pain in 70-95% of cases, medium-term investigations are now indicating an increased risk of adjacent level fracture following the procedure. With the aging population, understanding and treatment of osteoporosis is an increasingly important public health issue in developed Western countries. The aim of this study was to investigate the biomechanics of spinal osteoporosis and osteoporotic vertebral compression fractures by developing multi-scale computational, Finite Element (FE) models of both healthy and osteoporotic vertebral bodies. The multi-scale approach included the overall vertebral body anatomy, as well as a detailed representation of the internal trabecular microstructure. This novel, multi-scale approach overcame limitations of previous investigations by allowing simultaneous investigation of the mechanics of the trabecular micro-structure as well as overall vertebral body mechanics. The models were used to simulate the progression of osteoporosis, the effect of different loading conditions on vertebral strength and stiffness, and the effects of vertebroplasty on vertebral and trabecular mechanics. The model development process began with the development of an individual trabecular strut model using 3D beam elements, which was used as the building block for lattice-type, structural trabecular bone models, which were in turn incorporated into the vertebral body models. At each stage of model development, model predictions were compared to analytical solutions and in-vitro data from existing literature. The incremental process provided confidence in the predictions of each model before incorporation into the overall vertebral body model. The trabecular bone model, vertebral body model and vertebroplasty models were validated against in-vitro data from a series of compression tests performed using human cadaveric vertebral bodies. Firstly, trabecular bone samples were acquired and morphological parameters for each sample were measured using high resolution micro-computed tomography (CT). Apparent mechanical properties for each sample were then determined using uni-axial compression tests. Bone tissue properties were inversely determined using voxel-based FE models based on the micro-CT data. Specimen specific trabecular bone models were developed and the predicted apparent stiffness and strength were compared to the experimentally measured apparent stiffness and strength of the corresponding specimen. Following the trabecular specimen tests, a series of 12 whole cadaveric vertebrae were then divided into treated and non-treated groups and vertebroplasty performed on the specimens of the treated group. The vertebrae in both groups underwent clinical-CT scanning and destructive uniaxial compression testing. Specimen specific FE vertebral body models were developed and the predicted mechanical response compared to the experimentally measured responses. The validation process demonstrated that the multi-scale FE models comprising a lattice network of beam elements were able to accurately capture the failure mechanics of trabecular bone; and a trabecular core represented with beam elements enclosed in a layer of shell elements to represent the cortical shell was able to adequately represent the failure mechanics of intact vertebral bodies with varying degrees of osteoporosis. Following model development and validation, the models were used to investigate the effects of progressive osteoporosis on vertebral body mechanics and trabecular bone mechanics. These simulations showed that overall failure of the osteoporotic vertebral body is initiated by failure of the trabecular core, and the failure mechanism of the trabeculae varies with the progression of osteoporosis; from tissue yield in healthy trabecular bone, to failure due to instability (buckling) in osteoporotic bone with its thinner trabecular struts. The mechanical response of the vertebral body under load is highly dependent on the ability of the endplates to deform to transmit the load to the underlying trabecular bone. The ability of the endplate to evenly transfer the load through the core diminishes with osteoporosis. Investigation into the effect of different loading conditions on the vertebral body found that, because the trabecular bone structural changes which occur in osteoporosis result in a structure that is highly aligned with the loading direction, the vertebral body is consequently less able to withstand non-uniform loading states such as occurs in forward flexion. Changes in vertebral body loading due to disc degeneration were simulated, but proved to have little effect on osteoporotic vertebra mechanics. Conversely, differences in vertebral body loading between simulated invivo (uniform endplate pressure) and in-vitro conditions (where the vertebral endplates are rigidly cemented) had a dramatic effect on the predicted vertebral mechanics. This investigation suggested that in-vitro loading using bone cement potting of both endplates has major limitations in its ability to represent vertebral body mechanics in-vivo. And lastly, FE investigation into the biomechanical effect of vertebroplasty was performed. The results of this investigation demonstrated that the effect of vertebroplasty on overall vertebra mechanics is strongly governed by the cement distribution achieved within the trabecular core. In agreement with a recent study, the models predicted that vertebroplasty cement distributions which do not form one continuous mass which contacts both endplates have little effect on vertebral body stiffness or strength. In summary, this work presents the development of a novel, multi-scale Finite Element model of the osteoporotic vertebral body, which provides a powerful new tool for investigating the mechanics of osteoporotic vertebral compression fractures at the trabecular bone micro-structural level, and at the vertebral body level.