A HOMOGENIZATION SAMPLING PROCEDURE FOR CALCULATING TRABECULAR BONE EFFECTIVE STIFFNESS AND TISSUE-LEVEL STRESS

被引:256
作者
HOLLISTER, SJ [1 ]
BRENNAN, JM [1 ]
KIKUCHI, N [1 ]
机构
[1] UNIV MICHIGAN,DEPT MECH ENGN & APPL MECH,COMPUTAT MECH LAB,ANN ARBOR,MI 48109
关键词
D O I
10.1016/0021-9290(94)90019-1
中图分类号
Q6 [生物物理学];
学科分类号
071011 ;
摘要
A homogenization sampling procedure is introduced which allows computation of effective trabecular bone stiffness and individual trabecula level stress based on precise models of trabecular bone architecture. Three-dimensional digitized images of 53 trabecular bone specimens with a resolution of 50 mum per voxel were directly converted into three-dimensional finite element meshes by making each voxel an 8-node isoparametric brick element. Owing to the large mesh of 8000 elements, an element-by-element preconditioned conjugate gradient (EBEPCG) program was written to solve the local homogenization finite element equations. Predicted effective stiffness measures correlated well with experimental results (R2 > 0.73). The predicted effective stiffnesses tended to under estimate the experimental values. Average absolute errors in effective stiffness estimates ranged between 31 and 38% for the sampling procedure compared to a range 49-150% for a regression fit to volume fraction squared. Trabecula level stress ranged between - 200 and + 300 times that predicted by analyzing trabecular bone as a continuum. Both tensile and compressive tissue stresses were engendered by a continuum compressive stress. Trabecula level strain energy density (SED) ranged between 0 and 100 times the continuum SED value for two trabecular specimens. In conclusion, the homogenization sampling procedure consistently predicted the influence of trabecular bone architecture on effective stiffness. It can also provide trabecular tissue stress and strain estimates for arbitrary global loading of whole bones. Tissue stresses and strains showed large variations compared to corresponding continuum level quantities.
引用
收藏
页码:433 / 444
页数:12
相关论文
共 49 条
  • [1] Bakhvalov N. S., 1989, HOMOGENISATION AVERA
  • [2] AN APPROACH FOR TIME-DEPENDENT BONE MODELING AND REMODELING - THEORETICAL DEVELOPMENT
    BEAUPRE, GS
    ORR, TE
    CARTER, DR
    [J]. JOURNAL OF ORTHOPAEDIC RESEARCH, 1990, 8 (05) : 651 - 661
  • [3] FINITE-ELEMENT ANALYSIS OF A 3-DIMENSIONAL OPEN-CELLED MODEL FOR TRABECULAR BONE
    BEAUPRE, GS
    HAYES, WC
    [J]. JOURNAL OF BIOMECHANICAL ENGINEERING-TRANSACTIONS OF THE ASME, 1985, 107 (03): : 249 - 256
  • [4] ELEMENT-BY-ELEMENT LINEAR AND NONLINEAR SOLUTION-SCHEMES
    CAREY, GF
    JIANG, BN
    [J]. COMMUNICATIONS IN APPLIED NUMERICAL METHODS, 1986, 2 (02): : 145 - 153
  • [5] TRABECULAR BONE-DENSITY AND LOADING HISTORY - REGULATION OF CONNECTIVE-TISSUE BIOLOGY BY MECHANICAL ENERGY
    CARTER, DR
    FYHRIE, DP
    WHALEN, RT
    [J]. JOURNAL OF BIOMECHANICS, 1987, 20 (08) : 785 - +
  • [6] THE ELASTIC-MODULI OF HUMAN SUBCHONDRAL, TRABECULAR, AND CORTICAL BONE TISSUE AND THE SIZE-DEPENDENCY OF CORTICAL BONE MODULUS
    CHOI, K
    KUHN, JL
    CIARELLI, MJ
    GOLDSTEIN, SA
    [J]. JOURNAL OF BIOMECHANICS, 1990, 23 (11) : 1103 - 1113
  • [7] EVALUATION OF ORTHOGONAL MECHANICAL-PROPERTIES AND DENSITY OF HUMAN TRABECULAR BONE FROM THE MAJOR METAPHYSEAL REGIONS WITH MATERIALS TESTING AND COMPUTED-TOMOGRAPHY
    CIARELLI, MJ
    GOLDSTEIN, SA
    KUHN, JL
    CODY, DD
    BROWN, MB
    [J]. JOURNAL OF ORTHOPAEDIC RESEARCH, 1991, 9 (05) : 674 - 682
  • [8] THE RELATIONSHIP BETWEEN THE ELASTICITY TENSOR AND THE FABRIC TENSOR
    COWIN, SC
    [J]. MECHANICS OF MATERIALS, 1985, 4 (02) : 137 - 147
  • [9] COMPACT-BONE - NUMERICAL-SIMULATION OF MECHANICAL CHARACTERISTICS
    CROLET, JM
    AOUBIZA, B
    MEUNIER, A
    [J]. JOURNAL OF BIOMECHANICS, 1993, 26 (06) : 677 - 687
  • [10] CROLET JM, 1988, 12TH P ANN M AM SOC, P177