E. Sales

Micro-CT based finite element models of cancellous bone predict accurately displacement once the boundary condition is well replicated: A validation study

Non-destructive 3D micro-computed tomography (microCT) based finite element (microFE) models are used to estimate bone mechanical properties at tissue level. However, their validation remains challenging. Recent improvements in the quantification of displacements in bone tissue biopsies subjected to staged compression, using refined Digital Volume Correlation (DVC) techniques, now provide a full field displacement information accurate enough to be used for microFE validation. In this study, three specimens (two humans and one bovine) were tested with two different experimental set-ups, and the resulting data processed with the same DVC algorithm. The resulting displacement vector field was compared to that predicted by microFE models solved with three different boundary conditions (BC): nominal force resultant, nominal displacement resultant, distributed displacement. The first two conditions were obtained directly from the measurements provided by the experimental jigs, whereas in the third case the displacement field measured by the DVC in the top and bottom layer of the specimen was applied. Results show excellent relationship between the numerical predictions (x) and the experiments (y) when using BC derived from the DVC measurements (UX: y=1.07x-0.002, RMSE: 0.001mm; UY: y=1.03x-0.001, RMSE: 0.001mm; UZ: y=x+0.0002, RMSE: 0.001 mm for bovine specimen), whereas only poor correlation was found using BCs according to experiment set-ups. In conclusion, microFE models were found to predict accurately the vectorial displacement field using interpolated displacement boundary condition from DVC measurement.

Evaluating the macroscopic yield behaviour of trabecular bone using a nonlinear homogenisation approach

Computational homogenisation approaches using high resolution images and finite element (FE) modelling have been extensively employed to evaluate the anisotropic elastic properties of trabecular bone. The aim of this study was to extend its application to characterise the macroscopic yield behaviour of trabecular bone. Twenty trabecular bone samples were scanned using a micro-computed tomography device, converted to voxelised FE meshes and subjected to 160 load cases each (to define a homogenised multiaxial yield surface which represents several possible strain combinations). Simulations were carried out using a parallel code developed in-house. The nonlinear algorithms included both geometrical and material nonlinearities. The study found that for tension-tension and compression-compression regimes in normal strain space, the yield strains have an isotropic behaviour. However, in the tension-compression quadrants, pure shear and combined normal-shear planes, the macroscopic strain norms at yield have a relatively large variation. Also, our treatment of clockwise and counter-clockwise shears as separate loading cases showed that the differences in these two directions cannot be ignored. A quadric yield surface, used to evaluate the goodness of fit, showed that an isotropic criterion adequately represents yield in strain space though errors with orthotropic and anisotropic criteria are slightly smaller. Consequently, although the isotropic yield surface presents itself as the most suitable assumption, it may not work well for all load cases. This work provides a comprehensive assessment of material symmetries of trabecular bone at the macroscale and describes in detail its macroscopic yield and its underlying microscopic mechanics.