Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Civil-Comp Proceedings
ISSN 1759-3433
CCP: 102
PROCEEDINGS OF THE FOURTEENTH INTERNATIONAL CONFERENCE ON CIVIL, STRUCTURAL AND ENVIRONMENTAL ENGINEERING COMPUTING
Edited by:
Paper 105

Inverse Estimation of Rat Vertebrae Stiffness using Large-Scale Micro-Structural Finite Element Models

O. Jirousek1, P. Zlamal1, I. Jandejsek1, D. Kytyr1 and D. Schmidt2

1Department of Biomechanics, Institute of Theoretical and Applied Mechanics,
Academy of Sciences of the Czech Republic, Prague, Czech Republic
2Czech Technical University in Prague, Faculty of Transportation Sciences
Prague, Czech Republic

Full Bibliographic Reference for this paper
O. Jirousek, P. Zlamal, I. Jandejsek, D. Kytyr, D. Schmidt, "Inverse Estimation of Rat Vertebrae Stiffness using Large-Scale Micro-Structural Finite Element Models", in , (Editors), "Proceedings of the Fourteenth International Conference on Civil, Structural and Environmental Engineering Computing", Civil-Comp Press, Stirlingshire, UK, Paper 105, 2013. doi:10.4203/ccp.102.105
Keywords: large-scale finite element simulation, .

Summary

The possibility of using microstructural finite element models of whole bones to inversely calculate their mechanical properties is described in this paper. In this study, L2 vertebra of a Wistar rat were used to numerically determine its stiffness and the results are compared with the experimentally obtained value. The results demonstrate the possibility of using these large-scale finite element simulations to predict the overall bone stiffness. The micro-structural models used in this study are voxel-based, i.e. each finite element represents one spatial pixel (called a voxel). The micro-structural models are developed from a series of micro-CT images. The resolution of the input image data which is needed to capture the complex microstructure of trabecular bone in sufficient detail is discussed.

Tissue material properties are based on our previous nanoindentation study. This paper is focused on the parallel solution strategy employed to solve the large-scale finite element models utilizing existing open-source programs. Because of the high resolution of the microtomographic images the resulting finite element model of the vertebra is composed of approximately 9 million of hexahedral elements. To solve this large problem a parallel strategy must be employed. The main aim was to demonstrate the scalability of selected solvers (the preconditioned conjugate gradient and multifrontal massively parallel sparse direct solver) for these large voxel finite element models. Two architectures are tested: i) distributed memory system, and ii) shared memory system.

The results show, that the estimated overall stiffness of the vertebral body is not greatly influenced by the resolution of micro-CT images (albeit it must well capture the micro-structural characteristics of the trabecular bone) provided that the tissue material properties (elastic modulus and Poissons ratio at the level of individual trabeculae) are known.

purchase the full-text of this paper (price £20)

go to the previous paper
go to the next paper
return to the table of contents
return to the book description
purchase this book (price £65 +P&P)