Homogenized finite element analysis of distal tibia sections: Achievements and limitations.

High-resolution peripheral quantitative computed tomography (HR-pQCT) based micro-finite element (µFE) analysis allows accurate prediction of stiffness and ultimate load of standardised (∼1 cm) distal radius and tibia sections. An alternative homogenized finite element method (hFE) was recently validated to compute the ultimate load of larger (∼2 cm) distal radius sections that include Colles' fracture sites. Since the mechanical integrity of the weight-bearing distal tibia is gaining clinical interest, it has been shown that the same properties can be used to predict the strength of both distal segments of the radius and the tibia. Despite the capacity of hFE to predict structural properties of distal segments of the radius and the tibia, the limitations of such homogenization scheme remain unclear. Therefore, the objective of this study is to build a complete mechanical data set of the compressive behavior of distal segments of the tibia and to compare quantitatively the structural properties with the hFE predictions. As a further aim, it is intended to verify whether hFE is also able to capture the post-yield strain localisation or fracture zones in such a bone section, despite the absence of strain softening in the constitutive model. Twenty-five fresh-frozen distal parts of tibias of human donors were used in this study. Sections were cut corresponding to an in-house triple-stack protocol HR-pQCT scan, lapped, and scanned using micro computed tomography (µCT). The sections were tested in compression until failure, unloaded and scanned again in µCT. Volumetric bone mineral density (vBMD) and bone mineral content (BMC) were correlated to compression test results. hFE analysis was performed in order to compare computational predictions (stiffness, yield load and plastic deformation field pattern) with the compressive experiment. Namely, strain localization was assessed based on digital volume correlation (DVC) results and qualitatively compared to hFE predictions by comparing mid-slices patterns. Bone mineral content (BMC) showed a good correlation with stiffness (R2 = 0.92) and yield (R2 = 0.88). Structural parameters also showed good agreement between the experiment and hFE for both stiffness (R2 = 0.96, slope = 1.05 with 95 % CI [0.97, 1.14]) and yield (R2 = 0.95, slope = 1.04 [0.94, 1.13]). The qualitative comparison between hFE and DVC strain localization patterns allowed the classification of the samples into 3 categories: bad (15 sections), semi (8), and good agreement (2). The good correlations between BMC or hFE and experiment for structural parameters were similar to those obtained previously for the distal part of the radius. The failure zones determined by hFE corresponded to registration only in 8 % of the cases. We attribute these discrepancies to local elastic/plastic buckling effects that are not captured by the continuum-based FE approach exempt from strain softening. A way to improve strain localization hFE prediction would be to use longer distal segments with intact cortical shells, as done for the radius. To conclude, the used hFE scheme captures the elastic and yield response of the tibia sections reliably but not the subsequent failure process.

Validation of a fall rate prediction model for community-dwelling older adults: a combined analysis of three cohorts with 1850 participants.

BACKGROUND: Fragility fractures in older adults are often caused by fall events. The estimation of an expected fall rate might improve the identification of individuals at risk of fragility fractures and improve fracture prediction. METHODS: A combined analysis of three previously developed fall rate models using individual participant data (n = 1850) was conducted using the methodology of a two-stage meta-analysis to derive an overall model. These previously developed models included the fall history as a predictor recorded as the number of experienced falls within 12 months, treated as a factor variable with the levels 0, 1, 2, 3, 4 and ≥ 5 falls. In the first stage, negative binomial regression models for every cohort were fit. In the second stage, the coefficients were compared and used to derive overall coefficients with a random effect meta-analysis. Additionally, external validation was performed by applying the three data sets to the models derived in the first stage. RESULTS: The coefficient estimates for the prior number of falls were consistent among the three studies. Higgin's I2 as heterogeneity measure ranged from 0 to 55.39%. The overall coefficient estimates indicated that the expected fall rate increases with an increasing number of previous falls. External model validation revealed that the prediction errors for the data sets were independent of the model to which they were applied. CONCLUSION: This analysis suggests that the fall history treated as a factor variable is a robust predictor of estimating future falls among different cohorts.

Does tissue fixation change the mechanical properties of dry ovine bone extracellular matrix?

Tissue fixation is a prevalent method for bone conservation. Bone biopsies are typically fixed in formalin, dehydrated in ethanol, and infiltrated with polymethyl methacrylate (PMMA) Since some experiments can only be performed on fixed bone samples, it is essential to understand how fixation affects the measured material properties. The aim of this study was to quantify the influence of tissue fixation on the mechanical properties of cortical ovine bone at the extracellular matrix (ECM) level with state-of-the-art micromechanical techniques. A small section from the middle of the diaphysis of two ovine tibias (3.5 and 5.5 years old) was cut in the middle and polished on each side, resulting in a pair of mirrored surfaces. For each pair, one specimen underwent a fixation protocol involving immersion in formalin, dehydration with ethanol, and infiltration with PMMA. The other specimen (mirrored) was air-dried. Six osteons were selected in both pairs, which could be identified in both specimens. The influence of fixation on the mechanical properties was first analyzed using micropillar compression tests and nanoindentation in dry condition. Additionally, changes in the degree of mineralization were evaluated with Raman spectroscopy in both fixed and native bone ECM. Finally, micro tensile experiments were conducted in the 3.5-year fixed ovine bone ECM and compared to reported properties of unfixed dry ovine bone ECM. Interestingly, we found that tissue fixation does not alter the mechanical properties of ovine cortical bone ECM compared to experiments in dry state. However, animal age increases the degree of mineralization (p = 0.0159) and compressive yield stress (p = 0.041). Tissue fixation appears therefore as a valid conservation technique for investigating the mechanical properties of dehydrated bone ECM.

Tensile Mechanical Properties of Dry Cortical Bone Extracellular Matrix: A Comparison Among Two Osteogenesis Imperfecta and One Healthy Control Iliac Crest Biopsies.

Osteogenesis imperfecta (OI) is a genetic, collagen-related bone disease that increases the incidence of bone fractures. Still, the origin of this brittle mechanical behavior remains unclear. The extracellular matrix (ECM) of OI bone exhibits a higher degree of bone mineralization (DBM), whereas compressive mechanical properties at the ECM level do not appear to be inferior to healthy bone. However, it is unknown if collagen defects alter ECM tensile properties. This study aims to quantify the tensile properties of healthy and OI bone ECM. In three transiliac biopsies (healthy n = 1, OI type I n = 1, OI type III n = 1), 23 microtensile specimens (gauge dimensions 10 × 5 × 2 µm3) were manufactured and loaded quasi-statically under tension in vacuum condition. The resulting loading modulus and ultimate strength were extracted. Interestingly, tensile properties in OI bone ECM were not inferior compared to controls. All specimens revealed a brittle failure behavior. Fracture surfaces were graded according to their mineralized collagen fibers (MCF) orientation into axial, mixed, and transversal fracture surface types (FST). Furthermore, tissue mineral density (TMD) of the biopsy cortices was extracted from micro-computed tomogra[hy (µCT) images. Both FST and TMD are significant factors to predict loading modulus and ultimate strength with an adjusted R 2 of 0.556 (p = 2.65e-05) and 0.46 (p = 2.2e-04), respectively. The influence of MCF orientation and DBM on the mechanical properties of the neighboring ECM was further verified with quantitative polarized Raman spectroscopy (qPRS) and site-matched nanoindentation. MCF orientation and DBM were extracted from the qPRS spectrum, and a second mechanical model was developed to predict the indentation modulus with MCF orientation and DBM (R 2 = 67.4%, p = 7.73e-07). The tensile mechanical properties of the cortical bone ECM of two OI iliac crest biopsies are not lower than the one from a healthy and are primarily dependent on MCF orientation and DBM. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Mechanical testing and comparison of porcine tissue, silicones and 3D-printed materials for cardiovascular phantoms.

Background: Cardiovascular phantoms for patient education, pre-operative planning, surgical training, haemodynamic simulation, and device testing may help improve patient care. However, currently used materials may have different mechanical properties compared to biological tissue. Methods/Aim: The aim of this study was to investigate the mechanical properties of 3D-printing and silicone materials in comparison to biological cardiovascular tissues. Uniaxial cyclic tension testing was performed using dumbbell samples from porcine tissue (aorta, pulmonary artery, right and left ventricle). Flexible testing materials included 15 silicone (mixtures) and three 3D-printing materials. The modulus of elasticity was calculated for different deformation ranges. Results: The modulus of elasticity (0%-60%) for the aorta ranged from 0.16 to 0.18 N/mm2, for the pulmonary artery from 0.07 to 0.09 N/mm2, and for the right ventricle as well as the left ventricle short-axis from 0.1 to 0.16 N/mm2. For silicones the range of modulus of elasticity was 0.02-1.16 N/mm2, and for the 3D-printed materials from 0.85 to 1.02 N/mm2. The stress-strain curves of all tissues showed a non-linear behaviour in the cyclic tensile testing, with a distinct toe region, followed by exponential strain hardening behaviour towards the peak elongation. The vessel samples showed a more linear behaviour comparted to myocardial samples. The silicones and 3D printing materials exhibited near-linearity at higher strain ranges, with a decrease in stiffness following the initial deformation. All samples showed a deviation between the loading and unloading curves (hysteresis), and a reduction in peak force over the first few cycles (adaptation effect) at constant deformation. Conclusion: The modulus of elasticity of silicone mixtures is more in agreement to porcine cardiovascular tissues than 3D-printed materials. All synthetic materials showed an almost linear behaviour in the mechanical testing compared to the non-linear behaviour of the biological tissues, probably due to fibre recruitment mechanism in the latter.

Composition and micromechanical properties of the femoral neck compact bone in relation to patient age, sex and hip fracture occurrence.

Current clinical methods of bone health assessment depend to a great extent on bone mineral density (BMD) measurements. However, these methods only act as a proxy for bone strength and are often only carried out after the fracture occurs. Besides BMD, composition and tissue-level mechanical properties are expected to affect the whole bone's strength and toughness. While the elastic properties of the bone extracellular matrix (ECM) have been extensively investigated over the past two decades, there is still limited knowledge of the yield properties and their relationship to composition and architecture. In the present study, morphological, compositional and micropillar compression bone data was collected from patients who underwent hip arthroplasty. Femoral neck samples from 42 patients were collected together with anonymous clinical information about age, sex and primary diagnosis (coxarthrosis or hip fracture). The femoral neck cortex from the inferomedial region was analyzed in a site-matched manner using a combination of micromechanical testing (nanoindentation, micropillar compression) together with micro-CT and quantitative polarized Raman spectroscopy for both morphological and compositional characterization. Mechanical properties, as well as the sample-level mineral density, were constant over age. Only compositional properties demonstrate weak dependence on patient age: decreasing mineral to matrix ratio (p = 0.02, R2 = 0.13, 2.6 % per decade) and increasing amide I sub-peak ratio I∼1660/I∼1683 (p = 0.04, R2 = 0.11, 1.5 % per decade). The patient's sex and diagnosis did not seem to influence investigated bone properties. A clear zonal dependence between interstitial and osteonal cortical zones was observed for compositional and elastic bone properties (p < 0.0001). Site-matched microscale analysis confirmed that all investigated mechanical properties except yield strain demonstrate a positive correlation with the mineral fraction of bone. The output database is the first to integrate the experimentally assessed microscale yield properties, local tissue composition and morphology with the available patient clinical information. The final dataset was used for bone fracture risk prediction in-silico through the principal component analysis and the Naïve Bayes classification algorithm. The analysis showed that the mineral to matrix ratio, indentation hardness and micropillar yield stress are the most relevant parameters for bone fracture risk prediction at 70 % model accuracy (0.71 AUC). Due to the low number of samples, further studies to build a universal fracture prediction algorithm are anticipated with the higher number of patients (N > 200). The proposed classification algorithm together with the output dataset of bone tissue properties can be used for the future comparison of existing methods to evaluate bone quality as well as to form a better understanding of the mechanisms through which bone tissue is affected by aging or disease.

A multi-stack registration technique to improve measurement accuracy and precision across longitudinal HR-pQCT scans.

BACKGROUND: Recent applications of high-resolution peripheral quantitative computed tomography (HR-pQCT) have demonstrated that changes in local bone remodelling can be quantified in vivo using longitudinal three-dimensional image registration. However, certain emerging applications, such as fracture healing and joint analysis, require larger multi-stack scan regions that can result in stack shift image artifacts. These artifacts can be detrimental to the accurate alignment of the bone structure across multiple timepoints. The purpose of this study was to establish a multi-stack registration protocol for evaluating longitudinal HR-pQCT images and to assess the accuracy and precision error in comparison with measures obtained using previously established three-dimensional longitudinal registration. METHODS: Three same day multi-stack HR-pQCT scans of the radius (2 stacks in length) and tibia (3 stacks in length) were obtained from 39 healthy individuals who participated in a previous reproducibility study. A fully automated multi-stack registration algorithm was developed to re-align stacks within a scan by leveraging slight offsets between longitudinal scans. Stack shift severity before and after registration was quantified using a newly proposed stack-shift severity score. The false discovery rate for bone remodelling events and precision error of bone morphology and micro-finite element analysis parameters were compared between longitudinally registered scans with and without the addition of multi-stack registration. RESULTS: Most scans (82 %) improved in stack alignment or maintained the lowest stack shift severity score when multi-stack registration was implemented. The false discovery rate of bone remodelling events significantly decreased after multi-stack registration, resulting in median false detection of bone formation and resorption fractions between 3.2 to 7.5 % at the radius and 3.4 to 5.3 % at the tibia. Further, precision error was significantly reduced or remained unchanged in all standard bone morphology and micro-finite element analysis parameters, except for total and trabecular cross-sectional areas. CONCLUSION: Multi-stack registration is an effective strategy for accurately aligning multi-stack HR-pQCT scans without modification of the image acquisition protocol. The algorithm presented here is a viable approach for performing accurate morphological analysis on multi-stack HR-pQCT scans, particularly for advanced application investigating local bone remodelling in vivo.

Corrigendum: History of falls and fear of falling are predictive of future falls: Outcome of a fall rate model applied to the Swiss CHEF trial cohort.

[This corrects the article DOI: 10.3389/fragi.2022.1056779.].

In situ synchrotron radiation µCT indentation of cortical bone: Anisotropic crack propagation, local deformation, and fracture.

The development of treatment strategies for skeletal diseases relies on the understanding of bone mechanical properties in relation to its structure at different length scales. At the microscale, indention techniques can be used to evaluate the elastic, plastic, and fracture behaviour of bone tissue. Here, we combined in situ high-resolution SRµCT indentation testing and digital volume correlation to elucidate the anisotropic crack propagation, deformation, and fracture of ovine cortical bone under Berkovich and spherical tips. Independently of the indenter type we observed significant dependence of the crack development due to the anisotropy ahead of the tip, with lower strains and smaller crack systems developing in samples indented in the transverse material direction, where the fibrillar bone ultrastructure is largely aligned perpendicular to the indentation direction. Such alignment allows to accommodate the strain energy, inhibiting crack propagation. Higher tensile hoop strains generally correlated with regions that display significant cracking radial to the indenter, indicating a predominant Mode I fracture. This was confirmed by the three-dimensional analysis of crack opening displacements and stress intensity factors along the crack front obtained for the first time from full displacement fields in bone tissue. The X-ray beam significantly influenced the relaxation behaviour independent of the tip. Raman analyses did not show significant changes in specimen composition after irradiation compared to non-irradiated tissue, suggesting an embrittlement process that may be linked to damage of the non-fibrillar organic matrix. This study highlights the importance of three-dimensional investigation of bone deformation and fracture behaviour to explore the mechanisms of bone failure in relation to structural changes due to ageing or disease. STATEMENT OF SIGNIFICANCE: Characterising the three-dimensional deformation and fracture behaviour of bone remains essential to decipher the interplay between structure, function, and composition with the aim to improve fracture prevention strategies. The experimental methodology presented here, combining high-resolution imaging, indentation testing and digital volume correlation, allows us to quantify the local deformation, crack propagation, and fracture modes of cortical bone tissue. Our results highlight the anisotropic behaviour of osteonal bone and the complex crack propagation patterns and fracture modes initiating by the intricate stress states beneath the indenter tip. This is of wide interest not only for the understanding of bone fracture but also to understand other architectured (bio)structures providing an effective way to quantify their toughening mechanisms in relation to their main mechanical function.

Precision of bone mechanoregulation assessment in humans using longitudinal high-resolution peripheral quantitative computed tomography in vivo.

Local mechanical stimuli in the bone microenvironment are essential for the homeostasis and adaptation of the skeleton, with evidence suggesting that disruption of the mechanically-driven bone remodelling process may lead to bone loss. Longitudinal clinical studies have shown the combined use of high-resolution peripheral quantitative computed tomography (HR-pQCT) and micro-finite element analysis can be used to measure load-driven bone remodelling in vivo; however, quantitative markers of bone mechanoregulation and the precision of these analyses methods have not been validated in human subjects. Therefore, this study utilised participants from two cohorts. A same-day cohort (n = 33) was used to develop a filtering strategy to minimise false detections of bone remodelling sites caused by noise and motion artefacts present in HR-pQCT scans. A longitudinal cohort (n = 19) was used to develop bone imaging markers of trabecular bone mechanoregulation and characterise the precision for detecting longitudinal changes in subjects. Specifically, we described local load-driven formation and resorption sites independently using patient-specific odds ratios (OR) and 99 % confidence intervals. Conditional probability curves were computed to link the mechanical environment to the remodelling events detected on the bone surface. To quantify overall mechanoregulation, we calculated a correct classification rate measuring the fraction of remodelling events correctly identified by the mechanical signal. Precision was calculated as root-mean-squared averages of the coefficient of variation (RMS-SD) of repeated measurements using scan-rescan pairs at baseline combined with a one-year follow-up scan. We found no significant mean difference (p < 0.01) between scan-rescan conditional probabilities. RMS-SD was 10.5 % for resorption odds, 6.3 % for formation odds, and 1.3 % for correct classification rates. Bone was most likely to be formed in high-strain and resorbed in low-strain regions for all participants, indicating a consistent, regulated response to mechanical stimuli. For each percent increase in strain, the likelihood of bone resorption decreased by 2.0 ± 0.2 %, and the likelihood of bone formation increased by 1.9 ± 0.2 %, totalling 38.3 ± 1.1 % of strain-driven remodelling events across the entire trabecular compartment. This work provides novel robust bone mechanoregulation markers and their precision for designing future clinical studies.

The elasto-plastic nano- and microscale compressive behaviour of rehydrated mineralised collagen fibres.

The hierarchical design of bio-based nanostructured materials such as bone enables them to combine unique structure-mechanical properties. As one of its main components, water plays an important role in bone's material multiscale mechanical interplay. However, its influence has not been quantified at the length-scale of a mineralised collagen fibre. Here, we couple in situ micropillar compression, and simultaneous synchrotron small angle X-ray scattering (SAXS) and X-ray diffraction (XRD) with a statistical constitutive model. Since the synchrotron data contain statistical information on the nanostructure, we establish a direct connection between experiment and model to identify the rehydrated elasto-plastic micro- and nanomechanical fibre behaviour. Rehydration led to a decrease of 65%-75% in fibre yield stress and compressive strength, and 70% in stiffness with a 3x higher effect on stresses than strains. While in agreement with bone extracellular matrix, the decrease is 1.5-3x higher compared to micro-indentation and macro-compression. Hydration influences mineral more than fibril strain with the highest difference to the macroscale when comparing mineral and tissue levels. The effect of hydration seems to be strongly mediated by ultrastructural interfaces while results provide insights towards mechanical consequences of reported water-mediated structuring of bone apatite. The missing reinforcing capacity of surrounding tissue for an excised fibril array is more pronounced in wet than dry conditions, mainly related to fibril swelling. Differences leading to higher compressive strength between mineralised tissues seem not to depend on rehydration while the lack of kink bands supports the role of water as an elastic embedding influencing energy-absorption mechanisms. STATEMENT OF SIGNIFICANCE: Characterising structure-property-function relationships in hierarchical biological materials helps us to elucidate mechanisms that enable their unique properties. Experimental and computational methods can advance our understanding of their complex behaviour with the potential to inform bio-inspired material development. In this study, we close a gap for bone's fundamental mechanical building block at micro- and nanometre length scales. We establish a direct connection between experiments and simulations by coupling in situ synchrotron tests with a statistical model and quantify the behaviour of rehydrated single mineralised collagen fibres. Results suggest a high influence of hydration on structural interfaces, and the role of water as an elastic embedding by outlining important differences between wet and dry elasto-plastic properties of mineral nanocrystals, fibrils and fibres.

Textile Design of an Intervertebral Disc Replacement Device from Silk Yarn.

Low back pain is often due to degeneration of the intervertebral discs (IVD). It is one of the most common age- and work-related problems in today's society. Current treatments are not able to efficiently restore the full function of the IVD. Therefore, the aim of the present work was to reconstruct the two parts of the intervertebral disc-the annulus fibrosus (AF) and the nucleus pulposus (NP)-in such a way that the natural structural features were mimicked by a textile design. Silk was selected as the biomaterial for realization of a textile IVD because of its cytocompatibility, biodegradability, high strength, stiffness, and toughness, both in tension and compression. Therefore, an embroidered structure made of silk yarn was developed that reproduces the alternating fiber structure of +30° and -30° fiber orientation found in the AF and mimics its lamellar structure. The developed embroidered ribbons showed a tensile strength that corresponded to that of the natural AF. Fiber additive manufacturing with 1 mm silk staple fibers was used to replicate the fiber network of the NP and generate an open porous textile 3D structure that may serve as a reinforcement structure for the gel-like NP.

2D-3D reconstruction of the proximal femur from DXA scans: Evaluation of the 3D-Shaper software.

Introduction: Osteoporosis is currently diagnosed based on areal bone mineral density (aBMD) computed from 2D DXA scans. However, aBMD is a limited surrogate for femoral strength since it does not account for 3D bone geometry and density distribution. QCT scans combined with finite element (FE) analysis can deliver improved femoral strength predictions. However, non-negligible radiation dose and high costs prevent a systematic usage of this technique for screening purposes. As an alternative, the 3D-Shaper software (3D-Shaper Medical, Spain) reconstructs the 3D shape and density distribution of the femur from 2D DXA scans. This approach could deliver a more accurate estimation of femoral strength than aBMD by using FE analysis on the reconstructed 3D DXA. Methods: Here we present the first independent evaluation of the software, using a dataset of 77 ex vivo femora. We extend a prior evaluation by including the density distribution differences, the spatial correlation of density values and an FE analysis. Yet, cortical thickness is left out of this evaluation, since the cortex is not resolved in our FE models. Results: We found an average surface distance of 1.16 mm between 3D DXA and QCT images, which shows a good reconstruction of the bone geometry. Although BMD values obtained from 3D DXA and QCT correlated well (r 2 = 0.92), the 3D DXA BMD were systematically lower. The average BMD difference amounted to 64 mg/cm3, more than one-third of the 3D DXA BMD. Furthermore, the low correlation (r 2 = 0.48) between density values of both images indicates a limited reconstruction of the 3D density distribution. FE results were in good agreement between QCT and 3D DXA images, with a high coefficient of determination (r 2 = 0.88). However, this correlation was not statistically different from a direct prediction by aBMD. Moreover, we found differences in the fracture patterns between the two image types. QCT-based FE analysis resulted mostly in femoral neck fractures and 3D DXA-based FE in subcapital or pertrochanteric fractures. Discussion: In conclusion, 3D-Shaper generates an altered BMD distribution compared to QCT but, after careful density calibration, shows an interesting potential for deriving a standardized femoral strength from a DXA scan.

Development of a personalized fall rate prediction model in community-dwelling older adults: a negative binomial regression modelling approach.

BACKGROUND: Around a third of adults aged 65 and older fall every year, resulting in unintentional injuries in 30% of the cases. Fractures are a frequent consequence of falls, primarily caused in individuals with decreased bone strength who are unable to cushion their falls. Accordingly, an individual's number of experienced falls has a direct influence on fracture risk. The aim of this study was the development of a statistical model to predict future fall rates using personalized risk predictors. METHODS: In the prospective cohort GERICO, several fall risk factor variables were collected in community-dwelling older adults at two time-points four years apart (T1 and T2). Participants were asked how many falls they experienced during 12 months prior to the examinations. Rate ratios for the number of reported falls at T2 were computed for age, sex, reported fall number at T1, physical performance tests, physical activity level, comorbidity and medication number with negative binomial regression models. RESULTS: The analysis included 604 participants (male: 122, female: 482) with a median age of 67.90 years at T1. The mean number of falls per person was 1.04 and 0.70 at T1 and T2. The number of reported falls at T1 as a factor variable was the strongest risk factor with an unadjusted rate ratio [RR] of 2.60 for 3 falls (95% confidence interval [CI] 1.54 to 4.37), RR of 2.63 (95% CI 1.06 to 6.54) for 4 falls, and RR of 10.19 (95% CI 6.25 to 16.60) for 5 and more falls, when compared to 0 falls. The cross-validated prediction error was comparable for the global model including all candidate variables and the univariable model including prior fall numbers at T1 as the only predictor. CONCLUSION: In the GERICO cohort, the prior fall number as single predictor information for a personalized fall rate is as good as when including further available fall risk factors. Specifically, individuals who have experienced three and more falls are expected to fall multiple times again. TRIAL REGISTRATION: ISRCTN11865958, 13/07/2016, retrospectively registered.

Investigating the post-yield behavior of mineralized bone fibril arrays using a 3D non-linear finite element unit-cell model.

In this study, we propose a 3D non-linear finite element (FE) unit-cell model to investigate the post-yield behavior of mineralized collagen fibril arrays (FAY). We then compare the predictions of the model with recent micro-tensile and micropillar compression tests in both axial and transverse directions. The unit cell consists of mineralized collagen fibrils (MCFs) embedded in an extrafibrillar matrix (EFM), and the FE mesh is equipped with cohesive interactions and a custom plasticity model. The simulation results confirm that MCF plays a dominant role in load bearing prior to yielding under axial tensile loading. Damage was initiated via debonding in shear and progressive sliding at the MCF/EFM interface, and resulted in MCF pull-out until brittle failure. In transverse tensile loading, EFM carried most of the load in pre-yield deformation, and then mixed normal/shear debonding between MCF and EFM began to form, which eventually produced brittle delamination of the two phases. The loading/unloading FE analysis in compression along both axial and transverse directions demonstrated perfect plasticity without any reduction in elastic modulus, i.e., damage due to the interfaces as seen in micropillar compression. Beyond the brittle and ductile nature of the stress-strain curves, in tensile and compressive loading, the simulated post-yield behavior and failure mechanism are in good quantitative agreement with the experimental observations. Our rather simple but efficient unit-cell FE model can reproduce qualitatively and quantitatively the mechanical behavior of bone ECM under tensile and compressive loading along the two main orientations. The model's integration into higher length scales may be useful in describing the macroscopic post-yield and failure behavior of trabecular and cortical bone in greater detail.

Personalized loading conditions for homogenized finite element analysis of the distal sections of the radius.

The microstructure of trabecular bone is known to adapt its morphology in response to mechanical loads for achieving a biomechanical homeostasis. Based on this form-function relationship, previous investigators either simulated the remodeling of bone to predict the resulting density and architecture for a specific loading or retraced physiological loading conditions from local density and architecture. The latter inverse approach includes quantifying bone morphology using computed tomography and calculating the relative importance of selected load cases by minimizing the fluctuation of a tissue loading level metric. Along this concept, the present study aims at identifying an optimal, personalized, multiaxial load case at the distal section of the human radius using in vivo HR-pQCT-based isotropic, homogenized finite element (hFE) analysis. The dataset consisted of HR-pQCT reconstructions of the 20 mm most distal section of 21 human fresh-frozen radii. We simulated six different unit canonical load cases (FX palmar-dorsal force, FY ulnar-radial force, FZ distal-proximal force, MX moment about palmar-dorsal, MY moment about ulnar-radial, MZ moment about distal-proximal) using a simplified and efficient hFE method based on a single isotropic bone phase. Once we used a homogeneous mean density (shape model) and once the original heterogeneous density distribution (shape + density model). Using an analytical formulation, we minimized the deviation of the resulting strain tensors ε(x) to a hydrostatic compressive reference strain ε0, once for the 6 degrees of freedom (DOF) optimal (OPT) load case and for all individual 1 DOF load cases (FX, FY, FZ, MX, MY, MZ). All seven load cases were then extended in the nonlinear regime using the scaled displacements of the linear load cases as loading boundary conditions (MAX). We then compared the load cases and models for their objective function (OF) values, the stored energies and their ultimate strength using a specific torsor norm. Both shape and shape + density linear-optimized OPT models were dominated by a positive force in the z-direction (FZ). Transversal force DOFs were close to zero and mean moment DOFs were different depending on the model type. The inclusion of density distribution increased the influence and changed direction of MX and MY, while MZ was small in both models. The OPT load case had 12-15% lower objective function (OF) values than the FZ load case, depending on the model. Stored energies at the optimum were consistently 142-178% higher for the OPT load case than for the FZ load case. Differences in the nonlinear response maximum torsor norm ‖t‖ were heterogeneous, but consistently higher for OPT_MAX than FZ_MAX. We presented the proof of concept of an optimization procedure to estimate patient-specific loading conditions for hFE methods. In contrast to similar models, we included canonical load cases in all six DOFs and used a strain metric that favors hydrostatic compression. Based on a biomechanical analysis of the distal joint surfaces at the radius, the estimated load directions are plausible. For our dataset, the resulting OPT load case is close to the standard axial compression boundary conditions, usually used in HR-pQCT-based FE analysis today. But even using the present simplified hFE model, the optimized linear six DOF load case achieves a more homogeneous tissue loading and can absorb more than twice the energy than the standard uniaxial load case. The ultimate strength calculated with a torsor norm was consistently higher for the 6-DOF nonlinear model (OPT_MAX) than for the 1-DOF nonlinear uniaxial model (FZ_MAX). Defining patient-specific boundary conditions may decrease angulation errors during CT measurements and improve repeatability as well as reproducibility of bone stiffness and strength estimated by HR-pQCT-based hFE analysis. These results encourage the extension of the present method to anisotropic hFE models and their application to repeatability data sets to test the hypothesis of reduced angulation errors during measurement.

Influence of aging on mechanical properties of the femoral neck using an inverse method.

Today, we are facing rapid aging of the world population, which increases the incidence of hip fractures. The gold standard of bone strength assessment in the laboratory is micro-computed finite element analysis (µFEA) based on micro-computed tomography (µCT) images. In clinics, the standard method to assess bone fracture risk is based on areal bone mineral density (aBMD), measured by dual-energy X-ray absorptiometry (DXA). In addition, homogenized finite element analysis (hFEA) constructed from quantitative computed tomography reconstructions (QCT) predicts clinical bone strength more accurately than DXA. Despite considerable evidence of degradation of bone material properties with age, in the past fifty years of finite element analysis to predict bone strength, bone material parameters remained independent of age. This study aims to assess the influence of age on apparent modulus, yield stress, and strength predictions of the human femoral neck made by laboratory-available bone volume fraction (BV/TV) and µFEA; and by clinically available DXA and hFEA. Using an inverse method, we test the hypothesis that FEA material parameters are independent of age. Eighty-six human femora were scanned with DXA (aBMD) and with QCT. The femoral necks were extracted and scanned at 16 µm resolution with µCT. The grayscale images were downscaled to 32 µm and 65 µm for linear and non-linear analyses, respectively, and segmented. The µFE solver ParOSolNL (non-linear) and a standard hFEA method were applied to the neck sections with the same material properties for all samples to compute apparent modulus, yield stress, and strength. Laboratory-available BV/TV was a good predictor of apparent modulus (R2 = 0.76), almost as good as µFEA (R2 = 0.79). However, yield stress and strength were better predicted by µFEA (R2 = 0.92, R2 = 0.86, resp.) than BV/TV (R2 = 0.76, R2 = 0.76, resp.). For clinically available variables, prediction of apparent modulus was better with hFEA than aBMD (R2 = 0.67, R2 = 0.58, resp.). hFEA outperformed aBMD for predictions of yield stress (R2 = 0.63 vs R2 = 0.34 for female and R2 = 0.55 for male) and strength (R2 = 0.48 vs R2 = 0.33 for female and R2 = 0.15 for male). The inclusion of age did not improve the multiple linear models for apparent modulus, yield stress, and strength. The resolution of the µFE meshes seems to account for most morphological changes induced by aging. The errors between the simulation and the experiment for apparent modulus, yield stress, and strength were age-independent, suggesting no rationale for correcting tissue material parameters in the current FE analysis of the aging femoral neck.

Assessing minipig compact jawbone quality at the microscale.

Preclinical studies often require animal models for in vivo experiments. Particularly in dental research, pig species are extensively used due to their anatomical similarity to humans. However, there is a considerable knowledge gap on the multiscale morphological and mechanical properties of the miniature pigs' jawbones, which is crucial for implant studies and a direct comparison to human tissue. In the present work, we demonstrate a multimodal framework to assess the jawbone quantity and quality for a minipig animal model that could be further extended to humans. Three minipig genotypes, commonly used in dental research, were examined: Yucatan, Göttingen, and Sinclair. Three animals per genotype were tested. Cortical bone samples were extracted from the premolar region of the mandible, opposite to the teeth growth. Global morphological, compositional, and mechanical properties were assessed using micro-computed tomography (micro-CT) together with Raman spectroscopy and nanoindentation measurements, averaged over the sample area. Local mineral-mechanical relationships were investigated with the site-matched Raman spectroscopy and micropillar compression tests. For this, a novel femtosecond laser ablation protocol was developed, allowing high-throughput micropillar fabrication and testing without exposure to high vacuum. At the global averaged sample level, bone relative mineralization demonstrated a significant difference between the genotypes, which was not observed from the complementary micro-CT measurements. Moreover, bone hardness measured by nanoindentation showed a positive trend with the relative mineralization. For all genotypes, significant differences between the relative mineralization and elastic properties were more pronounced within the osteonal regions of cortical bone. Site-matched micropillar compression and Raman spectroscopy highlighted the differences between the genotypes' yield stress and mineral to matrix ratios. The methods used at the global level (averaged over sample area) could be potentially correlated to the medical tools used to assess jawbone toughness and morphology in clinics. On the other hand, the local analysis methods can be applied to quantify compressive bone mechanical properties and their relationship to bone mineralization.

Evaluation of Load-To-Strength Ratios in Metastatic Vertebrae and Comparison With Age- and Sex-Matched Healthy Individuals.

Vertebrae containing osteolytic and osteosclerotic bone metastases undergo pathologic vertebral fracture (PVF) when the lesioned vertebrae fail to carry daily loads. We hypothesize that task-specific spinal loading patterns amplify the risk of PVF, with a higher degree of risk in osteolytic than in osteosclerotic vertebrae. To test this hypothesis, we obtained clinical CT images of 11 cadaveric spines with bone metastases, estimated the individual vertebral strength from the CT data, and created spine-specific musculoskeletal models from the CT data. We established a musculoskeletal model for each spine to compute vertebral loading for natural standing, natural standing + weights, forward flexion + weights, and lateral bending + weights and derived the individual vertebral load-to-strength ratio (LSR). For each activity, we compared the metastatic spines' predicted LSRs with the normative LSRs generated from a population-based sample of 250 men and women of comparable ages. Bone metastases classification significantly affected the CT-estimated vertebral strength (Kruskal-Wallis, p < 0.0001). Post-test analysis showed that the estimated vertebral strength of osteosclerotic and mixed metastases vertebrae was significantly higher than that of osteolytic vertebrae (p = 0.0016 and p = 0.0003) or vertebrae without radiographic evidence of bone metastasis (p = 0.0010 and p = 0.0003). Compared with the median (50%) LSRs of the normative dataset, osteolytic vertebrae had higher median (50%) LSRs under natural standing (p = 0.0375), natural standing + weights (p = 0.0118), and lateral bending + weights (p = 0.0111). Surprisingly, vertebrae showing minimal radiographic evidence of bone metastasis presented significantly higher median (50%) LSRs under natural standing (p < 0.0001) and lateral bending + weights (p = 0.0009) than the normative dataset. Osteosclerotic vertebrae had lower median (50%) LSRs under natural standing (p < 0.0001), natural standing + weights (p = 0.0005), forward flexion + weights (p < 0.0001), and lateral bending + weights (p = 0.0002), a trend shared by vertebrae with mixed lesions. This study is the first to apply musculoskeletal modeling to estimate individual vertebral loading in pathologic spines and highlights the role of task-specific loading in augmenting PVF risk associated with specific bone metastatic types. Our finding of high LSRs in vertebrae without radiologically observed bone metastasis highlights that patients with metastatic spine disease could be at an increased risk of vertebral fractures even at levels where lesions have not been identified radiologically.

Unified validation of a refined second-generation HR-pQCT based homogenized finite element method to predict strength of the distal segments in radius and tibia.

INTRODUCTION: HR-pQCT based micro finite element (µFE) analyses are considered as "gold standard" for virtual biomechanical analyses of peripheral bone sites such as the distal segment of radius and tibia. An attractive alternative for clinical use is a homogenized finite element method (hFE) based on constitutive models, because of its much shorter evaluation times and modest computational resource requirements. Such hFE models have been experimentally validated for the distal segment of the radius, but neither for the distal segments of the tibia nor for both measurement sites together. Accordingly, the aim of the present study was to refine and experimentally validate an hFE processing pipeline for in vivo prediction of bone strength and stiffness at the distal segments of the radius and the tibia, using only one unified set of material properties. MATERIAL AND METHODS: An existing hFE analysis procedure was refined in several aspects: 1) to include a faster evaluation of material orientation based on the mean surface length (MSL) method, 2) to distinguish cortical and trabecular bone compartments with distinct material properties and 3) to directly superimpose material properties in mixed phase elements instead of densities. Based on an existing dataset of the distal segment of fresh-frozen radii (double sections 20.4 mm, n = 21) and a newly established dataset of the distal segment of fresh-frozen tibiae (triple sections, 30.6 mm, n = 25), a single set of material properties was calibrated on the radius dataset and validated on the tibia dataset by comparing hFE stiffness and ultimate load with respective experimental results, obtained by compressing the samples on a servo-hydraulic testing machine at a monotonic and quasi-static displacement rate up to failure. RESULTS: Using the identified set of material properties, the hFE-predicted stiffness and failure load were in excellent agreement with respective experimental results at both measurement sites (radius stiffness R2 = 0.93, slope = 1.00, intercept = 479 N/mm2/radius ultimate load: R2 = 0.97, slope = 1.00, intercept = 679 N; tibia stiffness R2 = 0.96, slope = 1.01, intercept = -1027 N/mm2/tibia ultimate load: R2 = 0.97, slope = 1.04, intercept = 394 N; combined dataset stiffness R2 = 0.95, slope = 1.01, intercept = -230 N/mm2/combined dataset ultimate load: R2 = 0.97, slope = 1.03, intercept = 495 N). DISCUSSION AND CONCLUSION: In conjunction with unified BV/TV calibration, the established hFE pipeline accurately predicts experimental stiffness and ultimate load of distal multi-sections at the radius and tibia. Processing time for non-linear analysis was substantially reduced compared to previous µFE and hFE methods but could be further minimized by estimating bone strength based on a fast and linear analysis like as is currently done with µ FE.