Reference values and clinical predictors of bone strength for HR-pQCT-based distal radius and tibia strength assessments in women and men.

Reference values for radius and tibia strength using multiple-stack high-resolution peripheral quantitative computed tomography (HR-pQCT) with homogenized finite element analysis are presented in order to derive critical values improving risk prediction models of osteoporosis. Gender and femoral neck areal bone mineral density (aBMD) were independent predictors of bone strength. INTRODUCTION: The purpose was to obtain reference values for radius and tibia bone strength computed by using the homogenized finite element analysis (hFE) using multiple stacks with a HR-pQCT. METHODS: Male and female healthy participants aged 20-39 years were recruited at the University Hospital of Bern. They underwent interview and clinical examination including hand grip, gait speed and DXA of the hip. The nondominant forearm and tibia were scanned with a double and a triple-stack protocol, respectively, using HR-pQCT (XCT II, SCANCO Medical AG). Bone strength was estimated by using the hFE analysis, and reference values were calculated using quantile regression. Multivariable analyses were performed to identify clinical predictors of bone strength. RESULTS: Overall, 46 women and 41 men were recruited with mean ages of 25.1 (sd 5.0) and 26.2 (sd 5.2) years. Sex-specific reference values for bone strength were established. Men had significantly higher strength for radius (mean (sd) 6640 (1800) N vs. 4110 (1200) N; p < 0.001) and tibia (18,200 (4220) N vs. 11,970 (3150) N; p < 0.001) than women. In the two multivariable regression models with and without total hip aBMD, the addition of neck hip aBMD significantly improved the model (p < 0.001). No clinical predictors of bone strength other than gender and aBMD were identified. CONCLUSION: Reference values for radius and tibia strength using multiple HR-pQCT stacks with hFE analysis are presented and provide the basis to help refining accurate risk prediction models. Femoral neck aBMD and gender were significant predictors of bone strength.

Perspectives on the non-invasive evaluation of femoral strength in the assessment of hip fracture risk.

We reviewed the experimental and clinical evidence that hip bone strength estimated by BMD and/or finite element analysis (FEA) reflects the actual strength of the proximal femur and is associated with hip fracture risk and its changes upon treatment. INTRODUCTION: The risk of hip fractures increases exponentially with age due to a progressive loss of bone mass, deterioration of bone structure, and increased incidence of falls. Areal bone mineral density (aBMD), measured by dual-energy X-ray absorptiometry (DXA), is the most used surrogate marker of bone strength. However, age-related declines in bone strength exceed those of aBMD, and the majority of fractures occur in those who are not identified as osteoporotic by BMD testing. With hip fracture incidence increasing worldwide, the development of accurate methods to estimate bone strength in vivo would be very useful to predict the risk of hip fracture and to monitor the effects of osteoporosis therapies. METHODS: We reviewed experimental and clinical evidence regarding the association between aBMD and/orCT-finite element analysis (FEA) estimated femoral strength and hip fracture risk as well as their changes with treatment. RESULTS: Femoral aBMD and bone strength estimates by CT-FEA explain a large proportion of femoral strength ex vivo and predict hip fracture risk in vivo. Changes in femoral aBMD are strongly associated with anti-fracture efficacy of osteoporosis treatments, though comparable data for FEA are currently not available. CONCLUSIONS: Hip aBMD and estimated femoral strength are good predictors of fracture risk and could potentially be used as surrogate endpoints for fracture in clinical trials. Further improvements of FEA may be achieved by incorporating trabecular orientations, enhanced cortical modeling, effects of aging on bone tissue ductility, and multiple sideway fall loading conditions.

Validation of distal radius failure load predictions by homogenized- and micro-finite element analyses based on second-generation high-resolution peripheral quantitative CT images.

This study developed a well-standardized and reproducible approach for micro-finite element (mFE) and homogenized-FE (hFE) analyses that can accurately predict the distal radius failure load using either mFE or hFE models when using the approaches and parameters developed in this study. INTRODUCTION: Micro-FE analyses based on high-resolution peripheral quantitative CT (HR-pQCT) images are frequently used to predict distal radius failure load. With the introduction of a second-generation HR-pQCT device, however, the default modelling approach no longer provides accurate results. The aim of this study was to develop a well-standardized and reproducible approach for mFE and hFE analyses that can provide precise and accurate results for distal radius failure load predictions based on second-generation HR-pQCT images. METHODS: Second-generation HR-pQCT was used to scan the distal 20-mm section of 22 cadaver radii. The sections were excised and mechanically tested afterwards. For these sections, mFE and hFE models were made that were used to identify required material parameters by comparing predicted and measured results. Using these parameters, the models were cropped to represent the 10-mm region recommended for clinical studies to test their performance for failure load prediction. RESULTS: After identification of material parameters, the measured failure load of the 20-mm segments was in good agreement with the results of mFE models (R2 = 0.969, slope = 1.035) and hFE models (R2 = 0.966, slope = 0.890). When the models were restricted to the clinical region, mFE still accurately predicted the measured failure load (R2 = 0.955, slope = 1.021), while hFE predictions were precise but tended to overpredict the failure load (R2 = 0.952, slope = 0.780). CONCLUSIONS: It was concluded that it is possible to accurately predict the distal radius failure load using either mFE or hFE models when using the approaches and parameters developed in this study.

Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations.

Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r2=0.81, femur-specific fabric: r2=0.88, mean fabric: r2=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.

The thermal conductivity of cortical and cancellous bone.

Surgical interventions close to vulnerable structures, such as nerves, require precise handling of surgical instruments and tools. These tools not only pose the risk of mechanical damage to soft tissues, but they also generate heat, which can lead to thermal necrosis of bone or soft tissues. Researchers and engineers are trying to improve those tools through experimentation and simulations. To simulate temperature distributions in anatomical structures, reliable material constants are needed. Therefore, this study aimed at investigating the thermal conductivity of cortical and cancellous bone. Accordingly, a custom-made steady-state experimental setup was designed and validated. 6 bovine and 3 human cortical bone samples, as well as 32 bovine cancellous bone samples, with variable bone volume fraction were tested. The cancellous bone samples were scanned by micro-computed tomography (µCT) and micro-finite element (µFE) voxel models were created to calculate iteratively the thermal conductivity of the bone marrow. The experimental results provided 0.64 ± 0.04 W/mK for bovine cortical bone and 0.68 ± 0.01 W/mK for human cortical bone. A linear dependency of thermal conductivity on bone volume fraction was found for cancellous bone [R-square (R2) = 0.8096, standard error of the estimates (SEE) = 0.0355 W/mK]. The thermal conductivity of the bone marrow was estimated to be 0.42 ± 0.05 W/mK. These results will help to improve thermal finite element simulations of the human skeleton and aid the development of new surgical tools or procedures.

Characterizing microcrack orientation distribution functions in osteonal bone samples.

Prefailure microdamage in bone tissue is considered to be the most detrimental factor in defining its strength and toughness with respect to age and disease. To understand the influence of microcracks on bone mechanics it is necessary to assess their morphology and three-dimensional distribution. This requirement reaches beyond classic histology and stereology, and methods to obtain such information are currently missing. Therefore, the aim of the study was to develop a methodology that allows to characterize three-dimensional microcrack distributions in bulk bone samples. Four dumbbell-shaped specimens of human cortical bone of a 77-year-old female donor were loaded beyond yield in either tension, compression or torsion (one control). Subsequently, synchrotron radiation micro-computed tomography (SRµCT) was used to obtain phase-contrast images of the damaged samples. A microcrack segmentation algorithm was developed and used to segment microcrack families for which microcrack orientation distribution functions were determined. Distinct microcrack families were observed for each load case that resulted in distinct orientation distribution functions. Microcracks had median areas of approximately 4.7 µm2 , 33.3 µm2 and 64.0 µm2 for tension, compression and torsion. Verifying the segmentation algorithm against a manually segmented ground truth showed good results when comparing the microcrack orientation distribution functions. A size dependence was noted when investigating the orientation distribution functions with respect to the size of the volume of interest used for their determination. Furthermore, a scale separation between tensile, compressive and torsional microcracks was noticeable. Visual comparison to classic histology indicated that microcrack families were successfully distinguished. We propose a methodology to analyse three-dimensional microcrack distributions in overloaded cortical bone. Such information could improve our understanding of bone microdamage and its impact on bone failure in relation to tissue age and disease.

Effect of boundary conditions on yield properties of human femoral trabecular bone.

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7-72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient [Formula: see text]. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.

Rebound-associated vertebral fractures after discontinuation of denosumab-from clinic and biomechanics.

UNLABELLED: Rebound-associated vertebral fractures may follow treatment discontinuation of highly potent reversible bone antiresorptives, resulting from the synergy of rapid bone resorption and accelerated microdamage accumulation in trabecular bone. INTRODUCTION: The purposes of this study are to characterize rebound-associated vertebral fractures following the discontinuation of a highly potent reversible antiresorptive therapy based on clinical observation and propose a pathophysiological rationale. METHODS: This study is a case report of multiple vertebral fractures early after discontinuation of denosumab therapy in a patient with hormone receptor-positive non-metastatic breast cancer treated with an aromatase inhibitor. RESULTS: Discontinuation of highly potent reversible bone antiresorptives such as denosumab may expose patients to an increased fracture risk due to the joined effects of absent microdamage repair during therapy followed by synchronous excess activation of multiple bone remodelling units at the time of loss-of-effect. We suggest the term rebound-associated vertebral fractures (RVF) for this phenomenon characterized by the presence of multiple new clinical vertebral fractures, associated with either no or low trauma, in a context consistent with the presence of high bone turnover and rapid loss of lumbar spine bone mineral density (BMD) occurring within 3 to 12 months after discontinuation (loss-of-effect) of a reversible antiresorptive therapy in the absence of secondary causes of bone loss or fractures. Unlike atypical femoral fractures that emerge from failure of microdamage repair in cortical bone with long-term antiresorptive treatment, RVF originate from the synergy of rapid bone resorption and accelerated microdamage accumulation in trabecular bone triggered by the discontinuation of highly potent reversible antiresorptives. CONCLUSIONS: Studies are urgently needed to i) prove the underlying pathophysiological processes suggested above, ii) establish the predictive criteria exposing patients to an increased risk of RVF, and iii) determine appropriate treatment regimens to be applied in such patients.

Experimental validation of a nonlinear µFE model based on cohesive-frictional plasticity for trabecular bone.

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of µFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Quantitative analysis of imprint shape and its relation to mechanical properties measured by microindentation in bone.

Microindentation in bone is a micromechanical testing technique routinely used to extract material properties related to bone quality. As the analysis of microindentation data is based on assumptions about the contact between sample and surface, the aim of this study was 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 topology was measured by atomic force microscopy. Statistical shape modeling of the residual imprint allowed to define a mean shape and to describe the variability in terms of 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was found to be highly consistent and free of any pile up while differing mostly by depth between species and direction. A few of the topological parameters, in particular depth, showed significant but rather weak and inconsistent correlations to variations in mechanical properties. The mechanical response of bone as well as the residual imprint shape was highly consistent within each category. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small deviations from an ideally flat surface.

Cathepsin S controls adipocytic and osteoblastic differentiation, bone turnover, and bone microarchitecture.

Cathepsin S is a cysteine protease that controls adipocyte differentiation and has been implicated in vascular and metabolic complications of obesity. Considering the inverse relation of osteoblasts and adipocytes and their mutual precursor cell, we hypothesized that cathepsin S may also affect osteoblast differentiation and bone remodeling. Thus, the fat and bone phenotypes of young (3 months old) and aged (12 or 18 months old) cathepsin S knock-out (KO) and wild-type (WT) mice were determined. Cathepsin S KO mice had a normal body weight at both ages investigated, even though the amount of subscapular and gonadal fat pads was reduced by 20%. Further, cathepsin S deficiency impaired adipocyte formation (-38%, p<0.001), which was accompanied by a lower expression of adipocyte-related genes and a reduction in serum leptin, IL-6 and CCL2 (p<0.001). Micro-CT analysis revealed an unchanged trabecular bone volume fraction and density, while tissue mineral density was significantly lower in cathepsin S KO mice at both ages. Aged KO mice further had a lower cortical bone mass (-2.3%, p<0.05). At the microarchitectural level, cathepsin S KO mice had thinner trabeculae (-8.3%), but a better connected trabecular network (+24%). Serum levels of the bone formation marker type 1 procollagen amino-terminal-propeptide and osteocalcin were both 2-3-fold higher in cathepsin S KO mice as was the mineralized surface. Consistently, osteogenic differentiation was increased 2-fold along with an increased expression of osteoblast-specific genes. Interestingly, serum levels of C-terminal telopeptide of type I collagen were also higher (+43%) in cathepsin S KO mice as were histological osteoclast parameters and ex vivo osteoclast differentiation. Thus, cathepsin S deficiency alters the balance between adipocyte and osteoblast differentiation, increases bone turnover, and changes bone microarchitecture. Therefore, bone and fat metabolisms should be monitored when using cathepsin S inhibitors clinically.

Tissue properties of the human vertebral body sub-structures evaluated by means of microindentation.

PURPOSE: The better understanding of vertebral mechanical properties can help to improve the diagnosis of vertebral fractures. As the bone mechanical competence depends not only from bone mineral density (BMD) but also from bone quality, the goal of the present study was to investigate the anisotropic indentation moduli of the different sub-structures of the healthy human vertebral body and spondylophytes by means of microindentation. METHODS: Six human vertebral bodies and five osteophytes (spondylophytes) were collected and prepared for microindentation test. In particular, indentations were performed on bone structural units of the cortical shell (along axial, circumferential and radial directions), of the endplates (along the anterio-posterior and lateral directions), of the trabecular bone (along the axial and transverse directions) and of the spondylophytes (along the axial direction). A total of 3164 indentations down to a maximum depth of 2.5 µm were performed and the indentation modulus was computed for each measurement. RESULTS: The cortical shell showed an orthotropic behavior (indentation modulus, Ei, higher if measured along the axial direction, 14.6±2.8 GPa, compared to the circumferential one, 12.3±3.5 GPa, and radial one, 8.3±3.1 GPa). Moreover, the cortical endplates (similar Ei along the antero-posterior, 13.0±2.9 GPa, and along the lateral, 12.0±3.0 GPa, directions) and the trabecular bone (Ei= 13.7±3.4 GPa along the axial direction versus Ei=10.9±3.7 GPa along the transverse one) showed transversal isotropy behavior. Furthermore, the spondylophytes showed the lower mechanical properties measured along the axial direction (Ei=10.5±3.3 GPa). CONCLUSIONS: The original results presented in this study improve our understanding of vertebral biomechanics and can be helpful to define the material properties of the vertebral substructures in computational models such as FE analysis.

Early changes in biochemical markers of bone formation during teriparatide therapy correlate with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis.

UNLABELLED: Changes of the bone formation marker PINP correlated positively with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis (GIO) who received 18-month treatment with teriparatide, but not with risedronate. These results support the use of PINP as a surrogate marker of bone strength in GIO patients treated with teriparatide. INTRODUCTION: To investigate the correlations between biochemical markers of bone turnover and vertebral strength estimated by finite element analysis (FEA) in men with GIO. METHODS: A total of 92 men with GIO were included in an 18-month, randomized, open-label trial of teriparatide (20 µg/day, n = 45) and risedronate (35 mg/week, n = 47). High-resolution quantitative computed tomography images of the 12th thoracic vertebra obtained at baseline, 6 and 18 months were converted into digital nonlinear FE models and subjected to anterior bending, axial compression and torsion. Stiffness and strength were computed for each model and loading mode. Serum biochemical markers of bone formation (amino-terminal-propeptide of type I collagen [PINP]) and bone resorption (type I collagen cross-linked C-telopeptide degradation fragments [CTx]) were measured at baseline, 3 months, 6 months and 18 months. A mixed-model of repeated measures analysed changes from baseline and between-group differences. Spearman correlations assessed the relationship between changes from baseline of bone markers with FEA variables. RESULTS: PINP and CTx levels increased in the teriparatide group and decreased in the risedronate group. FEA-derived parameters increased in both groups, but were significantly higher at 18 months in the teriparatide group. Significant positive correlations were found between changes from baseline of PINP at 3, 6 and 18 months with changes in FE strength in the teriparatide-treated group, but not in the risedronate group. CONCLUSIONS: Positive correlations between changes in a biochemical marker of bone formation and improvement of biomechanical properties support the use of PINP as a surrogate marker of bone strength in teriparatide-treated GIO patients.

Advanced CT based in vivo methods for the assessment of bone density, structure, and strength.

Based on spiral 3D tomography a large variety of applications have been developed during the last decade to asses bone mineral density, bone macro and micro structure, and bone strength. Quantitative computed tomography (QCT) using clinical whole body scanners provides separate assessment of trabecular, cortical, and subcortical bone mineral density (BMD) and content (BMC) principally in the spine and hip, although the distal forearm can also be assessed. Further bone macrostructure, for example bone geometry or cortical thickness can be quantified. Special high resolution peripheral CT (hr-pQCT) devices have been introduced to measure bone microstructure for example the trabecular architecture or cortical porosity at the distal forearm or tibia. 3D CT is also the basis for finite element analysis (FEA) to determine bone strength. QCT, hr-pQCT, and FEM are increasingly used in research as well as in clinical trials to complement areal BMD measurements obtained by the standard densitometric technique of dual x-ray absorptiometry (DXA). This review explains technical developments and demonstrates how QCT based techniques advanced our understanding of bone biology.

DXA predictions of human femoral mechanical properties depend on the load configuration.

The aim of this study was to evaluate the ability of dual energy X-rays absorptiometry (DXA) areal bone mineral density (aBMD) measured in different regions of the proximal part of the human femur for predicting the mechanical properties of matched proximal femora tested in two different loading configurations. 36 pairs of fresh frozen femora were DXA scanned and tested until failure in two loading configurations: a fall on the side or a one-legged standing. The ability of the DXA output from four different regions of the proximal femur in predicting the femoral mechanical properties was measured and compared for the two loading scenarios. The femoral neck DXA BMD was best correlated to the femoral ultimate force for both configurations and predicted significantly better femoral failure load (R(2)=0.80 vs. R(2)=0.66, P<0.05) when simulating a side than when simulating a standing configuration. Conversely, the work to failure was predicted similarly for both loading configurations (R(2)=0.54 vs. R(2)=0.53, P>0.05). Therefore, neck BMD should be considered as one of the key factors for discriminating femoral fracture risk in vivo. Moreover, the better predictive ability of neck BMD for femoral strength if tested in a fall compared to a one-legged stance configuration suggests that DXA's clinical relevance may not be as high for spontaneous femoral fractures than for fractures associated to a fall.

A generalized anisotropic quadric yield criterion and its application to bone tissue at multiple length scales.

Nonlinear computational analysis of materials showing elasto-plasticity or damage relies on knowledge of their yield behavior and strengths under complex stress states. In this work, a generalized anisotropic quadric yield criterion is proposed that is homogeneous of degree one and takes a convex quadric shape with a smooth transition from ellipsoidal to cylindrical or conical surfaces. If in the case of material identification, the shape of the yield function is not known a priori, a minimization using the quadric criterion will result in the optimal shape among the convex quadrics. The convexity limits of the criterion and the transition points between the different shapes are identified. Several special cases of the criterion for distinct material symmetries such as isotropy, cubic symmetry, fabric-based orthotropy and general orthotropy are presented and discussed. The generality of the formulation is demonstrated by showing its degeneration to several classical yield surfaces like the von Mises, Drucker-Prager, Tsai-Wu, Liu, generalized Hill and classical Hill criteria under appropriate conditions. Applicability of the formulation for micromechanical analyses was shown by transformation of a criterion for porous cohesive-frictional materials by Maghous et al. In order to demonstrate the advantages of the generalized formulation, bone is chosen as an example material, since it features yield envelopes with different shapes depending on the considered length scale. A fabric- and density-based quadric criterion for the description of homogenized material behavior of trabecular bone is identified from uniaxial, multiaxial and torsional experimental data. Also, a fabric- and density-based Tsai-Wu yield criterion for homogenized trabecular bone from in silico data is converted to an equivalent quadric criterion by introduction of a transformation of the interaction parameters. Finally, a quadric yield criterion for lamellar bone at the microscale is identified from a nanoindentation study reported in the literature, thus demonstrating the applicability of the generalized formulation to the description of the yield envelope of bone at multiple length scales.

A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro.

PURPOSE: Femoral fracture is a common medical problem in osteoporotic individuals. Bone mineral density (BMD) is the gold standard measure to evaluate fracture risk in vivo. Quantitative computed tomography (QCT)-based homogenized voxel finite element (hvFE) models have been proved to be more accurate predictors of femoral strength than BMD by adding geometrical and material properties. The aim of this study was to evaluate the ability of hvFE models in predicting femoral stiffness, strength and failure location for a large number of pairs of human femora tested in two different loading scenarios. METHODS: Thirty-six pairs of femora were scanned with QCT and total proximal BMD and BMC were evaluated. For each pair, one femur was positioned in one-legged stance configuration (STANCE) and the other in a sideways configuration (SIDE). Nonlinear hvFE models were generated from QCT images by reproducing the same loading configurations imposed in the experiments. For experiments and models, the structural properties (stiffness and ultimate load), the failure location and the motion of the femoral head were computed and compared. RESULTS: In both configurations, hvFE models predicted both stiffness (R(2)=0.82 for STANCE and R(2)=0.74 for SIDE) and femoral ultimate load (R(2)=0.80 for STANCE and R(2)=0.85 for SIDE) better than BMD and BMC. Moreover, the models predicted qualitatively well the failure location (66% of cases) and the motion of the femoral head. CONCLUSIONS: The subject specific QCT-based nonlinear hvFE model cannot only predict femoral apparent mechanical properties better than densitometric measures, but can additionally provide useful qualitative information about failure location.

An anisotropic elastic-viscoplastic damage model for bone tissue.

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain [Formula: see text] , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

The mechanical behavior of PMMA/bone specimens extracted from augmented vertebrae: a numerical study of interface properties, PMMA shrinkage and trabecular bone damage.

Recently published compression tests on PMMA/bone specimens extracted after vertebral bone augmentation indicated that PMMA/bone composites were not reinforced by the trabecular bone at all. In this study, the reasons for this unexpected behavior should be investigated by using non-linear micro-FE models. Six human vertebral bodies were augmented with either standard or low-modulus PMMA cement and scanned with a HR-pQCT system before and after augmentation. Six cylindrical PMMA/bone specimens were extracted from the augmented region, scanned with a micro-CT system and tested in compression. Four different micro-FE models were generated from these images which showed different bone tissue material behavior (with/without damage), interface behavior (perfect bonding, frictionless contact) and PMMA shrinkage due to polymerization. The non-linear stress-strain curves were compared between the different micro-FE models as well as to the compression tests of the PMMA/bone specimens. Micro-FE models with contact between bone and cement were 20% more compliant compared to those with perfect bonding. PMMA shrinkage damaged the trabecular bone already before mechanical loading, which further reduced the initial stiffness by 24%. Progressing bone damage during compression dominated the non-linear part of the stress-strain curves. The micro-FE models including bone damage and PMMA shrinkage were in good agreement with the compression tests. The results were similar with both cements. In conclusion, the PMMA/bone interface properties as well as the initial bone damage due to PMMA polymerization shrinkage clearly affected the stress-strain behavior of the composite and explained why trabecular bone did not contribute to the stiffness and strength of augmented bone.

Microindentation can discriminate between damaged and intact human bone tissue.

Bone mineral density and microarchitecture was found to predict 70-95% of bone strength. Microdamage, as factor of bone quality, might help to explain the remaining uncertainties. The goal of this study was to investigate whether microindentation can discriminate between intact and severely damaged human vertebral bone tissue in vitro. One portion from each human vertebral slice (N=35) tested in compression in a previous study was embedded, polished and tested in wet conditions by means of microindentation. The indentation moduli and hardness (HV) of trabecular, osteonal and interstitial bone structural units were computed along the cranio-caudal direction. Each indented region was defined as damaged or intact as seen under a light microscope. A total of 1190 indentations were performed. While both hardness and indentation modulus were independent from gender, both mechanical properties were affected by damage and microstructure. The damaged regions showed 50% lower stiffness and hardness compared to undamaged ones. Interstitial bone was stiffer and harder (13.2±4.4 GPa and 44.7±20.3 HV) than osteonal bone (10.9±3.8 GPa and 37.8±17.3 HV), which was stiffer and harder than trabecular bone (8.1±3.0 GPa and 28.8±11.2 HV) indented in the transverse direction. Moreover, along the axial direction intact trabecular bone (11.4±4.3 GPa) was 16% less stiff than the intact interstitial bone and as stiff as intact osteonal bone. In conclusion microindentation was found to discriminate between highly damaged and intact tissue in both trabecular and cortical bone tested in vitro. It remains to be investigated whether this technique would be able to detect also the damage, which is induced by physiological load in vivo.