The study of the mechanical behaviour of trabecular bone has extensively employed micro-level finite element (μFE) models generated from images of real bone samples. It is now recognized that the key determinants of the mechanical behaviour of bone are related to its micro-architecture. The key indices of micro-architecture, in turn, depend on factors such as age, anatomical site, sex, and degree of osteoporosis. In practice, it is difficult to acquire sufficient samples that encompass these variations. In this preliminary study, a method of generating virtual finite element (FE) samples of trabecular bone is considered. Virtual samples, calibrated to satisfy some of the key micro-architectural characteristics, are generated computationally. The apparent level elastic and post-elastic mechanical behaviour of the generated samples is examined: the elastic mechanical response of these samples is found to compare well with natural trabecular bone studies conducted by previous investigators; the post-elastic response of virtual samples shows that material non-linearities have a much greater effect in comparison with geometrical non-linearity for the bone densities considered. Similar behaviour has been reported by previous studies conducted on real trabecular bone. It is concluded that virtual modelling presents a potentially valuable tool in the study of the mechanical behaviour of trabecular bone and the role of its micro-architecture.
SimonU.AbelJ.ClaesL.Minimal size of cubic μFE-models of trabecular bone to predict the apparent stiffness. In Proceedings of the Seventh International Symposium on Computer methods in biomechanics and biomedical engineering (CMBBE2006), Antibes, France, 2006, pp. 661–665 (First Numerics, Cardiff).
2.
van RietbergenB.MajumdarS.PistoiaW.Assessment of cancellous bone mechanical properties from micro-FE models based on micro-CT, PQCT and MT images. Technol. Health Care, 1998, 6(5–6), 413–420.
3.
FaganM. J.DobsonC. A.OngE. K.BMU-Sim3D, a bone modelling environment. In Proceedings of the Sixth International Symposium on Computer methods in biomechanics and biomedical engineering (CMBBE 2006), Madrid, Spain, 2004, Vol. 19, Suppl. 1, pp. S365–S365 (First Numerics, Cardiff).
4.
HuiskesR.RuimermanR.van LentheG. H.Effects of mechanical forces on maintenance and adaptation of form in trabecular bone. Nature, 2000, 405, 704–706.
5.
UlrichD.van RietbergenB.WeinansH.Finite element analysis of trabecular bone structure: A comparison of image-based meshing techniques. J. Biomechanics, 1998, 31, 1187–1192.
6.
GibsonL. J.The mechanical behaviour of cancellous bone. J. Biomechanics, 1985, 18(5), 317–328.
7.
GuoX. E.ZyssetP. K.GoldsteinS. A.Study of post-yield behaviour of trabecular bone using a 3-D microstructural model. In Proceedings of the ASME Conference on Advances in bioengineering, San Francisco, Caifornia, USA, 12–17 November, 1995, vol. 31, pp. 165–166 (American Society of Mechanical Engineers, New York).
8.
KowalczykP.Orthotropic properties of cancellous bone modelled as parameterised cellular material. Computer Meth. Biomechanics Biomed. Engng, 2006, 9(3), 135–147.
9.
UlrichD.van RietbergenB.LaibA.The ability of three-dimensional structural indices to reflect mechanical aspects of trabecular bone. Bone, 1999, 25(1), 55–60.
10.
van RietbergenB.Micro-FE analysis of bone, state of the art. In Proceedings of the ASME Conference on Bioengineering, BED-vol 50. 2001, pp. 107–108. American Society of Mechanical Engineers, New York).
11.
CarterD. R.HayesW. C.The compressive behaviour of bone as a two-phase porous structure. J. Bone Jt Surg., 1977, 59, 954–962.
12.
RiceJ. C.CowinS. C.BowmanJ. A.On the dependence of the elasticity and strength of cancellous bone on apparent density. J. Biomechanics, 1988, 21(2), 155–168.
13.
WirtzD. C.SchiffersN.PandorfT.Critical evaluation of known bone material properties to realize anisotropic FE-simulation of the proximal femur. J. Biomechanics, 2000, 33, 1325–1330.
14.
WhitehouseW. J.The quantitative morphology of anisotropic trabecular bone. J. Microsc., 1974, 101, 153–168.
15.
OdgaardA.GundersonH. J. G.Quantification of connectivity in cancellous bone, with special emphasis on 3-D reconstructions. Bone, 1993, 14, 173–182.
16.
CowinS. C.LabordeA. J.The relationship between the elasticity tensor and the fabric tensor. Mechanics Mater., 1985, 4(2), 137–147.
17.
SahaP. K.ChaudhuriB. B.Detection of 3-D simple points for topology preserving transformations with application to thinning. IEEE Trans. Pattern Analysis Mach. Intell., 1994, 16, 1028–1032.
18.
StauberM.MüllerR.Volumetric spatial decomposition of trabecular bone into rods and plates - a new method for local bone morphometry. Bone, 2006, 38, 475–484.
19.
PankajP.PerroneC.BeesonF. A.Virtual modelling of normal and osteoporotic bone. In Proceedings of the seventh International Symposium on Computer methods in biomechanics and biomedical engineering, 2006, pp. 717–721 (First Numerics, Cardiff).
20.
HildebrandT.LaibA.MüllerR.Direct three-dimensional morphometric analysis of human cancellous bone: Microstructural data from spine, femur, iliac crest, and calcaneus. J. Bone Mineral Res., 1999, 14(7), 1167–1174.
21.
LaibA.NewittD. C.LuY.New model-independent measures of trabecular bone strength applied to in vivo high-resolution MR images. Osteoporosis Int., 2002, 13, 130–136.
22.
NewittD. C.van RietbergenB.MajumdarS.Processing and analysis of in vivo high-resolution MR images of trabecular bone for longitudinal studies: Reproducibility of structural measures and micro-finite element analysis derived mechanical properties. Osteoporosis Int., 2002, 13, 278–287.
23.
TurnerC. H.RhoJ.TakanoY.The elastic properties of trabecular and cortical bone tissues are similar: Results from two microscopic measurement techniques. J. Biomechanics, 1999, 32, 437–441.
24.
NieburG. L.A computational investigation of multiaxial failure in trabecular bone. PhD Thesis, Mechanical Engineering, Graduate Division, University of California, Berkley, California, USA, 2000.
25.
WilliamsJ. L.LewisJ. L.Properties and an anisotropic model of cancellous bone from the proximal tibial epiphysis. J. Biomed. Engng, 1982, 104, 50–56.
26.
KabelJ.van RietbergenB.DalstraM.The role of an effective isotropic tissue modulus in the elastic properties of cancellous bone. J. Biomechanics, 1999, 32, 673–680.
27.
van RietbergenB.OdgaardA.KabelJ.Relationships between bone morphology and bone elastic properties can be accurately quantified using high resolution computer reconstructions. J. Orthop. Res., 1998, 16, 23–28.
28.
KabelJ.van RietbergenB.OdgaardA.Constitutive relationships of fabric, density, and elastic properties in cancellous bone architecture. Bone, 1999, 25(4), 481–486.
29.
KeavenyT. M.WachtelE. F.FordC. M.Differences between the tensile and compressive strengths of bovine tibial trabecular bone depends on modulus. J. Biomechanics, 1994, 27(9), 1137–1146.
30.
PhillipsA. T. M.PankajP.UsmaniA. S.The effect of acetabular cup size on the short-term stability of revision hip arthroplasty: A finite element investigation. Proc. Instn Mech. Engrs Part H: J. Engineering in Medicine, 2004, 218, 239–249.
PankajP.LinE.Material modelling in the seismic response analysis for the design of RC framed structures. Engng Struct., 2005, 27(7), 1014–1023.
33.
KeyakJ. H.RossiS. A.Prediction of femoral fracture load using finite element models: An examination of stress- and strain-based failure theories. J. Biomechanics, 2000, 33(2), 209–214.
34.
TaiK.UlmF. J.OrtizC.Nanogranular origins of the strength of bone. Nono Lett., 2006, 6(11), 2520–2525.
35.
de BorstR.PankajP.BicanicN.A note on singularity indicators for Mohr Coulomb type yield criteria. Computers Struct., 1991, 39, 219–220.
36.
NieburG. L.KeavenyT. M.Biaxial yield strains in bovine tibial trabecular bone do not depend of volume fraction. In Proceedings of the ASME Conference on Bioengineering, BED-vol. 50, 2001, pp. 23–24 (American Society of Mechanical Engineers, New York).
37.
DingM.DalstraM.DanielsenC. C.Age variations in the properties of human tibial trabecular bone. J. Bone Jt Surg. Br., 1997, 79, 995–1002.
38.
EvansG. F.LebowM.Regional differences in some of the physical properties of the human femur. J. Appl. Physiology, 1951, 3, 563–572.
39.
KopperdahlD. L.KeavenyT. M.Yield strain behaviour of trabecular bone. J. Biomechanics, 1998, 31, 601–608.
40.
PistoiaW.van RietbergenB.LochmllerE.-M.Estimation of distal radius failure load with micro-finite element analysis models based on three dimensional peripheral quantitative computed tomography images. Bone, 2002, 30(6), 842–848.
41.
van RietbergenB.PistoiaW.UlrichD.Prediction of trabecular bone failure parameters using a tissue failure and μFE analysis. J. Computer Modelling Simulation Medicine, 2000, 1(2), 98–101.
42.
HarriganT. P.JastyM.MannR. W.Limitations of the continuum assumption in cancellous bone. 1988, 21(4), J. Biomechanics, 269–275.
CiarelliM. J.GoldsteinS. A.KuhnJ. L.Evaluation of orthogonal mechanical properties and density of human trabecular bone from the major metaphyseal regions with materials testing and computed tomography. J. Orthop. Res., 1991, 9(5), 674–682.
45.
KnaussP.Material properties and strength behaviour of spongy bone tissue at the coxal human femur. Biomed. Technqs, 1981, 26, 200–210.
46.
LotzJ. C.GerhartT. N.HayesW. C.Mechanical properties of trabecular bone from the proximal femur: A quantitative CT study. J. Computer Assisted Tomography, 1990, 14, 107–114.
47.
GuoX. E.KimC. H.Mechanical consequence of trabecular bone loss and its treatment: A three-dimensional model simulation. Bone, 2002, 30(2), 404–411.
48.
ParfittA. M.Age-related structural changes in trabecular and cortical bone: Cellular mechanisms and biomechanical consequences. Calcified Tissue Int., 1984, 36, 123–128.
49.
BevillG.EswaranS. K.GuptaA.Influence of bone volume fraction and architecture on computed large-deformation failure mechanisms in human trabecular bone. Bone, 2006, 39, 1218–1225.
50.
SölkenJ. S.KinneyJ. H.On the importance of geometric nonlinearity in finite-element simulations of trabecular bone failure. Bone, 2003, 33, 494–504.
