PubMed Health⌕ Search

Biomedical subjects

B van Rietbergen

Publications and source records attributed to B van Rietbergen.

At least 19 recordsLinked to original sources

Comparison of micro-level and continuum-level voxel models of the proximal femur.

Continuum-level finite element (FE) models became standard computational tools for the evaluation of bone mechanical behavior from in vivo computed tomography scans. Such scans do not account for the anisotropy of the bone. Instead, local mechanical properties in the continuum-level FE models are assumed isotropic and are derived from bone density, using statistical relationships. Micro-FE models, on the other hand, incorporate the anisotropic structure in detail. This study aimed to quantify the effects of assumed isotropy, by comparing continuum-level voxel models of a healthy and a severely osteoporotic proximal femur with recently analyzed micro-FE models of the same bones. The micro-model element size was coarsened to generate continuum FE models with two different element sizes (0.64 and 3.04 mm) and two different density-modulus relationships found in the literature for wet and ash density. All FE models were subjected to the same boundary conditions that simulated a fall to the side, and the stress and strain distributions, model stiffness and yield load were compared. The results indicated that the stress and strain distributions could be reproduced well with the continuum models. The smallest differences between the continuum-level model and micro-level model predictions of the stiffness and yield load were obtained with the coarsest element size. Better results were obtained for both continuum-element sizes when isotropic moduli were based on ash density rather than wet density.

Anisotropy↗

A theoretical framework for strain-related trabecular bone maintenance and adaptation.

It is assumed that density and morphology of trabecular bone is partially controlled by mechanical forces. How these effects are expressed in the local metabolic functions of osteoclast resorption and osteoblast formation is not known. In order to investigate possible mechano-biological pathways for these mechanisms we have proposed a mathematical theory (Nature 405 (2000) 704). This theory is based on hypothetical osteocyte stimulation of osteoblast bone formation, as an effect of elevated strain in the bone matrix, and a role for microcracks and disuse in promoting osteoclast resorption. Applied in a 2-D Finite Element Analysis model, the theory explained the formation of trabecular patterns. In this article we present a 3-D FEA model based on the same theory and investigated its potential morphological predictability of metabolic reactions to mechanical loads. The computations simulated the development of trabecular morphological details during growth, relative to measurements in growing pigs, reasonably realistic. They confirmed that the proposed mechanisms also inherently lead to optimal stress transfer. Alternative loading directions produced new trabecular orientations. Reduction of load reduced trabecular thickness, connectivity and mass in the simulation, as is seen in disuse osteoporosis. Simulating the effects of estrogen deficiency through increased osteoclast resorption frequencies produced osteoporotic morphologies as well, as seen in post-menopausal osteoporosis. We conclude that the theory provides a suitable computational framework to investigate hypothetical relationships between bone loading and metabolic expressions.

Adaptation, Physiological↗

A fibril-reinforced poroviscoelastic swelling model for articular cartilage.

From a mechanical point of view, the most relevant components of articular cartilage are the tight and highly organized collagen network together with the charged proteoglycans. Due to the fixed charges of the proteoglycans, the cation concentration inside the tissue is higher than in the surrounding synovial fluid. This excess of ion particles leads to an osmotic pressure difference, which causes swelling of the tissue. The fibrillar collagen network resists straining and swelling pressures. This combination makes cartilage a unique, highly hydrated and pressurized tissue, enforced with a strained collagen network. Many theories to explain articular cartilage behavior under loading, expressed in computational models that either include the swelling behavior or the properties of the anisotropic collagen structure, can be found in the literature. The most common tests used to determine the mechanical quality of articular cartilage are those of confined compression, unconfined compression, indentation and swelling. All theories currently available in the literature can explain the cartilage response occurring in some of the above tests, but none of them can explain these for all of the tests. We hypothesized that a model including simultaneous mathematical descriptions of (1) the swelling properties due to the fixed-change densities of the proteoglycans and (2) the anisotropic viscoelastic collagen structure, can explain all these test simultaneously. To study this hypothesis we extended our fibril-reinforced poroviscoelastic finite element model with our biphasic swelling model. We have shown that the newly developed fibril-reinforced poroviscoelastic swelling (FPVES) model for articular cartilage can simultaneously account for the reaction force during swelling, confined compression, indentation and unconfined compression as well as the lateral deformation during unconfined compression. Using this theory it is possible to analyze the link between the collagen network and the swelling properties of articular cartilage.

Animals↗

A three-dimensional digital image correlation technique for strain measurements in microstructures.

A three-dimensional digital image correlation technique is presented for strain measurements in open-cell structures such as trabecular bone. The technique uses high-resolution computed tomography images for displacement measurements in the solid structure. In order to determine the local strain-state within single trabeculae, a tetrahedronization method is used to fill the solid structure with tetrahedrae. Displacements are calculated at the nodes of the tetrahedrae. The displacement data is subsequently converted to a deformation tensor in each of the tetrahedral element centers with a least-squares estimation method. Because the trabeculae are represented by a mesh, it is possible to deform this mesh according to the deformation tensor and, at the same time, visualize the calculated local strain in the deformed mesh with a finite element post-processing tool. In this way, the deformation of a single trabecula from an aluminum foam sample was determined and validated with rendered images of the three-dimensional sample. A precision analysis showed that a rigid translation or rotation does not affect the accuracy. Typical values for the standard deviation in the displacement and strain components are 2.0 microm and 0.01, respectively. Presently, the precision limits the technique to strain measurements beyond the yield strain.

Algorithms↗

Thermal analysis of bone cement polymerisation at the cement-bone interface.

The two major problems that have been reported with the use of polymethylmethacrylate (PMMA) cement are thermal necrosis of surrounding bone due to the high heat generation during polymerisation and chemical necrosis due to unreacted monomer release. Computer models have been used to study the temperature and monomer distribution after cementation. In most of these models, however, polymerisation is modelled as temperature independent and cancellous bone is modelled as a continuum. Such models thus cannot account for the expected important role of the trabecular bone micro-structure. The aim of this study is to investigate the distribution of temperature and monomer leftover at the cancellous bone-cement interface during polymerisation for a realistic trabecular bone-cement micro-structure and realistic temperature-dependent polymerisation kinetics behaviour. A 3-D computer model of a piece of bovine cancellous bone that underwent pressurization with bone-cement was generated using a micro-computed tomography scanner. This geometry was used as the basis for a finite element model and a temperature-dependent problem for bone cement polymerisation kinetics was solved to simulate the bone cement polymerisation process in the vicinity of the interface. The transient temperature field throughout the interface was calculated, along with the polymerisation fraction distribution in the cement domain. The calculations revealed that the tips of the bone trabeculae that are embedded in the cement attain temperatures much higher than the average temperature of the bone volume. A small fraction of the bone (10%) is exposed to temperatures exceeding 70 degrees C, but the exposure time to these high temperatures is limited to 50s. In the region near the bone, the cement polymerisation fraction (about 84%) is less than that in the centre (where it is reaching values of over 96%). An important finding of this study thus is the fact that the bone tissue that is subjected to the highest temperatures is also subjected to high leftover monomer concentration. Furthermore the maximum bone temperature is reached relatively early, when monomer content in the neighbouring cement is still quite high.

Animals↗

Stresses in the local collagen network of articular cartilage: a poroviscoelastic fibril-reinforced finite element study.

Osteoarthritis (OA) is a multifactorial disease, resulting in diarthrodial joint wear and eventually destruction. Swelling of cartilage, which is proportional to the amount of collagen damage, is an initial event of cartilage degeneration, so damage to the collagen fibril network is likely to be one of the earliest signs of OA cartilage degeneration. We propose that the local stresses and strains in the collagen fibrils, which cause the damage, cannot be determined dependably without taking the local arcade-like collagen-fibril structure into account. We investigate this using a poroviscoelastic fibril-reinforced FEA model. The constitutive fibril properties were determined by fitting numerical data to experimental results of unconfined compression and indentation tests on samples of bovine patellar articular cartilage. It was demonstrated that with this model the stresses and strains in the collagen fibrils can be calculated. It was also exhibited that fibrils with different orientations at the same location can be loaded differently, depending on the local architecture of the collagen network. To the best of our knowledge, the present model is the first that can account for these features. We conclude that the local stresses and strains in the articular cartilage are highly influenced by the local morphology of the collagen-fibril network.

Animals↗

Image-based micro-finite-element modeling for improved distal radius strength diagnosis: moving from bench to bedside.

Although osteoporosis is characterized by quantitative (mass) and qualitative (structural) changes, standard clinical techniques (dual-energy X-ray absorptiometry, DXA) only measure the former. Three-dimensional micro-finite-element (micro-FE) models based on high-resolution images can account for structural aspects as well, and it has recently been shown that an improved prediction of distal radius strength is possible with micro-FE analysis. A clinical application of this technique, however, is limited by its high imaging and computational demands. The objective of this study is to investigate if an improved prediction of bone strength can be obtained as well when only a small part of the radius is used for micro-FE modeling. Images of a 1-cm region of the metaphysis of the distal radius of 54 cadaver arms (mean age: 82 +/- 9 SD) made with a three-dimensional peripheral quantitative computed tomography (pQCT) device at 165- micro m resolution formed the basis for micro-FE models that were used to predict the bone failure load. Following imaging, specimens were experimentally compressed to failure to produce a Colles'-type fracture. Failure loads predicted from micro-FE analyses agreed well with those measured experimentally (R2 = 0.66, p < 0.001). Lower correlations were observed with bone mass (R2 = 0.48, p < 0.001) and microstructural parameters (R2 = 0.47, p < 0.001). Hence, even when only a small region is modeled, micro-FE analysis provides an improved prediction of radius strength.

Absorptiometry, Photon↗

Mechanical consequences of different scenarios for simulated bone atrophy and recovery in the distal radius.

Metabolic bone diseases such as osteoporosis usually cause a decrease in bone mass and a deterioration of bone microarchitecture leading to a decline in bone strength. Methods to predict bone strength in patients are currently based on bone mass only. It has been suggested that an improved prediction of bone strength might be possible if structural changes are taken into account as well. In this study we evaluated which structural parameters (other than bone mass) are the best predictors for changes in bone mechanical properties of the human radius after different bone atrophy scenarios and whether the original strength of the affected bone can be recovered if bone loss is restored by thickening of the remaining structures. To answer these questions, a human radius was measured with a microcomputer tomography scanner to extract the full three-dimensional architecture of the distal radius at an isotropic resolution of 80 microm. Eight models with modified bone architecture were created and the mechanical variations due to these modifications were studied using microfinite element (micro-FE) simulations. In four models mass was lowered by 20%, either by reducing cortical thickness, trabecular thickness, or number of trabeculae or by overall thinning of structures. In the other four models bone mass was restored to the original value using a trabecular bone thickening procedure. The micro-FE analyses revealed that most load was carried by the cortical bone. For this reason, bone strength was affected most in the reduced cortical thickness model. For the same reason, the trabecular bone atrophy scenarios, all of which affected bone strength in a very similar way, resulted in less dramatic bone strength reduction. The restoration of bone mass did not recover the original bone strength. These findings demonstrate that the importance of different parameters for the prediction of bone strength also depends on the mechanical loading. This could explain why results of earlier studies on the importance of structural parameters can be inconsistent and site-dependent.

Adult↗

Pathways of load-induced cartilage damage causing cartilage degeneration in the knee after meniscectomy.

Results of both clinical and animal studies show that meniscectomy often leads to osteoarthritic degenerative changes in articular cartilage. It is generally assumed that this process of cartilage degeneration is due to changes in mechanical loading after meniscectomy. It is, however, not known why and where this cartilage degeneration starts. Load induced cartilage damage is characterized as either type (1)--damage without disruption of the underlying bone or calcified cartilage layer--or type (2), subchondral fracture with or without damage to the overlying cartilage. We asked the question whether cartilage degeneration after meniscectomy is likely to be initiated by type (1) and/or type (2) cartilage damage. To investigate that we applied an axisymmetric biphasic finite element analysis model of the knee joint. In this model the articular cartilage layers of the tibial and the femoral condyles, the meniscus and the bone underlying the articular cartilage of the tibia plateau were included. The model was validated with data from clinical studies, in which the effects of meniscectomy on contact areas and pressures were measured. It was found that both the maximal values and the distributions of the shear stress in the articular cartilage changed after meniscectomy, and that these changes could lead to both type (1) and type (2) cartilage damage. Hence it likely that the cartilage degeneration seen after meniscectomy is initiated by both type (1) and type (2) cartilage damage.

Anisotropy↗

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.

The authors have developed a system for the characterization of trabecular bone structure from high-resolution MR images. It features largely automated coil inhomogeneity correction, trabecular bone region segmentation, serial image registration, bone/marrow binarization, and structural calculation steps. The system addresses problems of efficiency and inter- and intraoperator variability inherent in previous analyses. The system is evaluated on repetitive scans of 8 volunteers for both two-dimensional (2D) apparent structure calculations and three-dimensional (3D) mechanical calculations using micro-finite element analysis. Coil correction methods based on a priori knowledge of the coil sensitivity and on low-pass filtering of the high-resolution mages are compared and found to perform similarly. Image alignment is found to cause small but significant changes in some structural parameters. Overall the automated system provides on the order of a 3-fold decrease in trained operator time over previous manual methods. Reproducibility is found to be dependent on image quality for most parameters. For 7 subjects with good image quality, reproducibility of 2-4% is found for 2D structural parameters, while 3D mechanical parameters vary by 4-9%, with percent standardized coefficients of variation in the ranges of 15-34% and 20-38% respectively.

Adult↗

In vivo assessment of architecture and micro-finite element analysis derived indices of mechanical properties of trabecular bone in the radius.

Measurement of microstructural parameters of trabecular bone noninvasively in vivo is possible with high-resolution magnetic resonance (MR) imaging. These measurements may prove useful in the determination of bone strength and fracture risk, but must be related to other measures of bone properties. In this study in vivo MR imaging was used to derive trabecular bone structure measures and combined with micro-finite element analysis (microFE) to determine the effects of trabecular bone microarchitecture on bone mechanical properties in the distal radius. The subjects were studied in two groups: (I) postmenopausal women with normal bone mineral density (BMD) (n = 22, mean age 58 +/- 7 years) and (II) postmenopausal women with spine or femur BMD -1 SD to -2.5 SD below young normal (n = 37, mean age 62 +/- 11 years). MR images of the distal radius were obtained at 1.5 T, and measures such as apparent trabecular bone volume fraction (App BV/TV), spacing, number and thickness (App TbSp, TbN, TbTh) were derived in regions of interest extending from the joint line to the radial shaft. The high-resolution images were also used in a micro-finite element model to derive the directional Young's moduli (E1, E2 and E3), shear moduli (G12, G23 and G13) and anisotropy ratios such as E1/E3. BMD at the distal radius, lumbar spine and hip were assessed using dual-energy X-ray absorptiometry (DXA). Bone formation was assessed by serum osteocalcin and bone resorption by serum type I collagen C-terminal telopeptide breakdown products (serum CTX) and urinary CTX biochemical markers. The trabecular architecture displayed considerable anisotropy. Measures of BMD such as the ultradistal radial BMD were lower in the osteopenic group (p<0.01). Biochemical markers between the two groups were comparable in value and showed no significant difference between the two groups. App BV/TV, TbTh and TbN were higher, and App TbSp lower, in the normal group than the osteopenic group. All three directional measures of elastic and shear moduli were lower in the osteopenic group compared with the normal group. Anisotropy of trabecular bone microarchitecture, as measured by the ratios of the mean intercept length (MIL) values (MIL1/MIL3, etc.), and the anisotropy in elastic modulus (E1/E3, etc.), were greater in the osteopenic group compared with the normal group. The correlations between the measures of architecture and moduli are higher than those between elastic moduli and BMD. Stepwise multiple regression analysis showed that while App BV/TV is highly correlated with the mechanical properties, additional structural measures do contribute to the improved prediction of the mechanical measures. This study demonstrates the feasibility and potential of using MR imaging with microFE modeling in vivo in the study of osteoporosis.

Aged↗

High-resolution MRI and micro-FE for the evaluation of changes in bone mechanical properties during longitudinal clinical trials: application to calcaneal bone in postmenopausal women after one year of idoxifene treatment.

OBJECTIVE: To investigate whether recently developed in vivo high-resolution magnetic resonance-imaging and micro-finite element techniques can monitor changes in bone mechanical properties during long-term clinical trials aiming at evaluating the efficacy of new drugs for the treatment of osteoporosis. DESIGN: Comparison of baseline and follow-up mechanical parameters calculated using micro-finite element analysis of the calcaneus for subjects participating in a study investigating the effect of idoxifene. BACKGROUND: Contemporary measurements for the evaluation of bone mechanical properties, based on dual-energy X-ray absorptiometry measurements, are not very accurate and require large trial populations. METHODS: A total of 56 postmenopausal subjects received either a placebo, 5 mg or 10 mg per day of idoxifene. Magnetic resonance-images of the calcaneus were made at baseline and after one year. Mechanical parameters of a trabecular volume of interest in the calcaneus were calculated using micro-finite element analysis. RESULTS: Although there were no significant differences between the mean changes in the treated groups and the placebo group, there were significant changes from baseline within groups after one year of treatment. Significant changes, however, were found only for mechanical parameters and only in the treated groups. CONCLUSIONS: The present study is the first demonstration that longitudinal changes in bone mechanical properties due to trabecular micro-architectural changes may be quantified in long-term clinical studies. Since significant changes in mechanical parameters were obtained for the treated groups whereas no significant change in bone mass was found we conclude that the application of these techniques may increase the clinical significance of these trials. RELEVANCE: A precise diagnosis of in vivo bone mechanical properties that accounts for (changes in) trabecular bone architecture is of particular importance for longitudinal clinical trials aiming at evaluating the efficacy of new drugs since it can lead to clinically relevant results from shorter follow-up intervals and may enable a reduction of the number of patients involved in the trial.

Bone Density↗

Estimation of distal radius failure load with micro-finite element analysis models based on three-dimensional peripheral quantitative computed tomography images.

There is increasing evidence that, in addition to bone mass, bone microarchitecture and its mechanical load distribution are important factors for the determination of bone strength. Recently, it has been shown that new high-resolution imaging techniques in combination with new modeling algorithms based on the finite element (FE) method can account for these additional factors. Such models thus could provide more relevant information for the estimation of bone failure load. The purpose of the present study was to determine whether results of whole-bone micro-FE (microFE) analyses with models based on three-dimensional peripheral quantitative computer tomography (3D-pQCT) images (isotropic voxel resolution of 165 microm) could predict the failure load of the human radius more accurately than results with dual-energy X-ray absorptiometry (DXA) or bone morphology measurements. For this purpose, microFE models were created using 54 embalmed cadaver arms. It was assumed that bone failure would be initiated if a certain percentage of the bone tissue (varied from 1% to 7%) would be strained beyond the tissue yield strain. The external force that produced this tissue strain was calculated from the FE analyses. These predictions were correlated with results of real compression testing on the same cadaver arms. The results of these compression tests were also correlated with results of DXA and structural measurements of these arms. The compression tests produced Colles-type fractures in the distal 4 cm of the radius. The predicted failure loads calculated from the FE analysis agreed well with those measured in the experiments (R(2) = 0.75 p < 0.001). Lower correlations were found with bone mass (R(2) = 0.48, p < 0.001) and bone structural parameters (R(2) = 0.57 p < 0.001). We conclude that application of the techniques investigated here can lead to a better prediction of the bone failure load for bone in vivo than is possible from DXA measurements, structural parameters, or a combination thereof.

Aged↗

Micro-FE analyses of bone: state of the art.

The ability to provide a complete characterization of elastic properties of bone has vastly improved our understanding of trabecular bone mechanical properties. Based on this information, it was possible to validate several mechanical concepts related to the elastic behavior of trabecular bone that could not be validated earlier. With recently developed micro-CT scanners and the availability of large parallel computer systems, this technique has also enabled the determination of physiological bone tissue loading conditions from very large microFE models that can represent whole human bones in detail. Such analyses can provide the data needed for a better understanding of bone failure processes or cell mediated load adaptive remodeling processes. Computational demands for whole bone analyses, however, are still excessive. Unlike linear stress and strain analyses, the application of PFE to study non-linear processes, in particular bone failure mechanisms, is still in an early phase Results of recent studies, however, are promising and indicate that an accurate prediction of bone failure with these techniques is possible. Compelling features of such analyses are that they enable multi-axial failure criteria at the apparent level to be developed using primarily computational methods as well as that they can provide a basis for detailed analysis of micro-mechanics associated with trabecular failure at the apparent level. The application of microFE techniques to analyze bone in vivo is in an early stage as well. First results have indicated that, although the resolution of presently available in vivo imaging techniques (i.e. pQCT and MR) is much less than that of images used so far for uFE analyses, the technique can provide meaningful elastic properties of trabecular bone in vivo in most cases. It is expected that the remaining uncertainties in the microFE results can be eliminated as soon as the resolution of in vivo images is improved. With the fast developments in pQCT(47) and MR scanning, this will probably be in the near future.

Biomechanical Phenomena↗

Load transfer and stress shielding of the hydroxyapatite-ABG hip: a study of stem length and proximal fixation.

We questioned whether the excellent clinical results reported for the hydroxyapatite-coated ABG hip are consistent with results of preclinical computer simulation methods for the prediction of strain-adaptive bone remodeling patterns around implants. We also investigated whether a further improvement of the results can be expected if complete proximal load transfer is enforced by reducing stem length to the metaphyseal area only. Resorption patterns predicted were in good agreement with radiographic clinical findings. Mechanical failure of the bone-implant interface was judged based on a Hoffman interface stress criterion, which was low. Reducing stem length hardly increased interface failure probability, but it did not reduce stress shielding, either. From a mechanical point of view, reducing stem length to the metaphyseal area only is not advantageous.

Arthroplasty, Replacement, Hip↗

High-resolution three-dimensional-pQCT images can be an adequate basis for in-vivo microFE analysis of bone.

Micro-finite element (microFE) models based on high-resolution images have enabled the calculation of elastic properties of trabecular bone in vitro. Recently, techniques have been developed to image trabecular bone structure in vivo, albeit at a lesser resolution. The present work studies the usefulness of such in-vivo images for microFE analyses, by comparing their microFE results to those of models based on high-resolution micro-CT (microCT) images. Fifteen specimens obtained from human femoral heads were imaged first with a 3D-pQCT scanner at 165 microns resolution and a second time with a microCT scanner at 56 microns resolution. A third set of images with a resolution of 165 microns was created by downscaling the microCT measurements. The microFE models were created directly from these images. Orthotropic elastic properties and the average tissue von Mises stress of the specimens were calculated from six FE-analyses per specimen. The results of the 165 microns models were compared to those of the 56 microns model, which was taken as the reference model. The results calculated from the pQCT-based models, correlated excellent with those calculated from the reference model for both moduli (R2 > 0.95) and for the average tissue von Mises stress (R2 > 0.83). Results calculated from the downscaled micro-CT models correlated even better with those of the reference models (R2 > 0.99 for the moduli and R2 > 0.96 for the average von Mises stress). In the case of the 3D-pQCT based models, however, the slopes of the regression lines were less than one and had to be corrected. The prediction of the Poisson's ratios was less accurate (R2 > 0.45 and R2 > 0.67) for the models based on 3D-pQCT and downscaled microCT images respectively). The fact that the results from the downscaled and original microCT images were nearly identical indicates that the need for a correction in the case of the 3D-pQCT measurements was not due to the voxel size of the images but due to a higher noise level and a lower contrast in these images, in combination with the application of a filtering procedure at 165 micron images. In summary: the results of microFE models based on in-vivo images of the 3D-pQCT can closely resemble those obtained from microFE models based on higher resolution microCT system.

Adult↗

Correlation between pre-operative periprosthetic bone density and post-operative bone loss in THA can be explained by strain-adaptive remodelling.

Periprosthetic adaptive bone remodelling after total hip arthroplasty can be simulated in computer models, combining bone remodelling theory with finite element analysis. Patient specific three-dimensional finite element models of retrieved bone specimens from an earlier bone densitometry (DEXA) study were constructed and bone remodelling simulations performed. Results of the simulations were analysed both qualitatively and quantitatively. Patterns of predicted bone loss corresponded very well with the DEXA measurements on the retrievals. The amount of predicted bone loss, measured quantitatively by simulating DEXA on finite element models, was found to be inversely correlated with the initial bone mineral content. It was concluded that the same clinically observed correlation can therefore be explained by mechanically induced remodelling. This finding extends the applicability of numerical pre-clinical testing to the analysis of interaction between implant design and initial state of the bone.

Absorptiometry, Photon↗

The role of an effective isotropic tissue modulus in the elastic properties of cancellous bone.

Conceptually, the elastic characteristics of cancellous bone could be predicted directly from the trabecular morphology--or architecture--and by the elastic properties of the tissue itself. Although hardly any experimental evidence exists, it is often implicitly assumed that tissue anisotropy has a negligible effect on the apparent elastic properties of cancellous bone. The question addressed in this paper is whether this is actually true. If it is, then micromechanical finite element analysis (micro-FEA) models, representing trabecular architecture, using an 'effective isotropic tissue modulus' should be able to predict apparent elastic properties of cancellous bone. To test this, accurate multi-axial compressive mechanical tests of 29 whale bone specimens were simulated with specimen-specific micro-FEA computer models built from true three-dimensional reconstructions. By scaling the micro-FEA predictions by a constant tissue modulus, 92% of the variation of Young's moduli determined experimentally could be explained. The correlation even increased to 95% when the micro-FEA moduli were scaled to the isotropic tissue moduli of individual specimens. Excellent agreement was also found in the elastic symmetry axes and anisotropy ratios. The prediction of Poisson's ratios was somewhat less precise at 85% correlation. The results support the hypothesis; for practical purposes, the concept of an 'effective isotropic tissue modulus' concept is a viable one. They also suggest that the value of such a modulus for individual cases might be inferred from the average tissue density, hence the degree of mineralization. Future studies must clarify how specific the tissue modulus should be for different types of bone if adequate predictions of elastic behavior are to be made in this way.

Animals↗