Material properties assignment to finite element models of bone structures: a new method.

Finite element analysis (FEA) is widely adopted to investigate the mechanical behaviour of bone structures. Computed tomography (CT) data are frequently used to generate FE models of bone. If properly calibrated, CT images are capable of providing accurate information about the bone morphology and tissue density. The aim of this work was to develop a special program able to read a CT data set as well as the FEA mesh generated from it, and to assign to each element of the mesh the material properties derived from the bone tissue density at the element location. The program was tested on phantom data sets and was adopted to evaluate the effects of the discrete description of the bone material properties. A three-dimensional FE model was generated automatically from a 16 bit CT data set of a distal femur acquired in vivo. The strain energy density (SED) was evaluated for each model element for increasing model complexity (number of different material cards assigned to the model). The computed SED were strongly dependent on the material mapping strategy.

[1]  J. O´Rourke,et al.  Computational Geometry in C: Arrangements , 1998 .

[2]  S A Marom,et al.  Computer aided stress analysis of long bones utilizing computed tomography. , 1990, Journal of biomechanics.

[3]  P Ducheyne,et al.  The mechanical behaviour of intracondylar cancellous bone of the femur at different loading rates. , 1977, Journal of biomechanics.

[4]  H. Skinner,et al.  Correlations between orthogonal mechanical properties and density of trabecular bone: use of different densitometric measures. , 1994, Journal of biomedical materials research.

[5]  J. Currey,et al.  Young's modulus, density and material properties in cancellous bone over a large density range , 1992 .

[6]  T. Keller,et al.  Young's modulus, bending strength, and tissue physical properties of human compact bone , 1990, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[7]  L Cristofolini,et al.  A comparative study on different methods of automatic mesh generation of human femurs. , 1998, Medical engineering & physics.

[8]  W. Hayes,et al.  The compressive behavior of bone as a two-phase porous structure. , 1977, The Journal of bone and joint surgery. American volume.

[9]  E Y Chao,et al.  A survey of finite element analysis in orthopedic biomechanics: the first decade. , 1983, Journal of biomechanics.

[10]  J H Keyak,et al.  Automated three-dimensional finite element modelling of bone: a new method. , 1990, Journal of biomedical engineering.

[11]  R Van Audekercke,et al.  The mechanical characteristics of cancellous bone at the upper femoral region. , 1983, Journal of biomechanics.

[12]  R. B. Ashman,et al.  Relations of mechanical properties to density and CT numbers in human bone. , 1995, Medical engineering & physics.

[13]  Joseph O'Rourke,et al.  Computational Geometry in C. , 1995 .

[14]  W. Hayes,et al.  Prediction of vertebral body compressive fracture using quantitative computed tomography. , 1985, The Journal of bone and joint surgery. American volume.

[15]  S. Goldstein,et al.  Evaluation of orthogonal mechanical properties and density of human trabecular bone from the major metaphyseal regions with materials testing and computed tomography , 1991, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[16]  E. Schneider,et al.  Estimation of mechanical properties of cortical bone by computed tomography , 1991, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[17]  S. Cowin,et al.  On the dependence of the elasticity and strength of cancellous bone on apparent density. , 1988, Journal of biomechanics.

[18]  William E. Lorensen,et al.  Marching cubes: a high resolution 3D surface construction algorithm , 1996 .