Orthotropic properties of cancellous bone modelled as parameterized cellular material

Analysis of stresses and strains in bone tissues and simulation of their adaptive remodelling require exhaustive information about distribution of constitutive properties of cancellous bone and their relationships to microstructural parameters. Homogenization of "equivalent" trabecular microstructures appears to be an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. The models, based on a space-filling dodecahedron, are fully three-dimensional and are parameterized with four shape parameters. Variation of the parameters allows to imitate most of typical microstructure patterns observed in real bones, along with a variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Coefficients of elastic orthotropic stiffness matrix are determined as tabularized functions of elastic constants versus the shape parameters. Additionally, they are correlated to apparent density and principal fabric tensor values. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement.     

Orthotropic properties of cancellous bone modelled as parameterized cellular material

Analysis of stresses and strains in bone tissues and simulation of their adaptive remodelling require exhaustive information about distribution of constitutive properties of cancellous bone and their relationships to microstructural parameters. Homogenization of “equivalent” trabecular microstructures appears to be an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. The models, based on a space-filling dodecahedron, are fully three-dimensional and are parameterized with four shape parameters. Variation of the parameters allows to imitate most of typical microstructure patterns observed in real bones, along with a variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Coefficients of elastic orthotropic stiffness matrix are determined as tabularized functions of elastic constants versus the shape parameters. Additionally, they are correlated to apparent density and principal fabric tensor values. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement.     

Trabecular surface remodeling simulation for cancellous bone using microstructural voxel finite element models

A computational simulation method for three-dimensional trabecular surface remodeling was proposed, using voxel finite element models of cancellous bone, and was applied to the experimental data. In the simulation, the trabecular microstructure was modeled based on digital images, and its morphological changes due to surface movement at the trabecular level were directly expressed by removing/adding the voxel elements from/to the trabecular surface. A remodeling simulation at the single trabecular level under uniaxial compressive loading demonstrated smooth morphological changes even though the trabeculae were modeled with discrete voxel elements. Moreover, the trabecular axis rotated toward the loading direction with increasing stiffness, simulating functional adaptation to the applied load. In the remodeling simulation at the trabecular structural level, a cancellous bone cube was modeled using a digital image obtained by microcomputed tomography (microCT), and was uniaxially compressed. As a result, the apparent stiffness against the applied load increased by remodeling, in which the trabeculae reoriented to the loading direction. In addition, changes in the structural indices of the trabecular architecture coincided qualitatively with previously published experimental observations. Through these studies, it was demonstrated that the newly proposed voxel simulation technique enables us to simulate the trabecular surface remodeling and to compare the results obtained using this technique with the in vivo experimental data in the investigation of the adaptive bone remodeling phenomenon.     

Mimetization of the elastic properties of cancellous bone via a parameterized cellular material

Bone tissue mechanical properties and trabecular microarchitecture are the main factors that determine the biomechanical properties of cancellous bone. Artificial cancellous microstructures, typically described by a reduced number of geometrical parameters, can be designed to obtain a mechanical behavior mimicking that of natural bone. In this work, we assess the ability of the parameterized microstructure introduced by Kowalczyk (Comput Methods Biomech Biomed Eng 9:135–147, 2006. doi: 10.1080/10255840600751473) to mimic the elastic response of cancellous bone. Artificial microstructures are compared with actual bone samples in terms of elasticity matrices and their symmetry classes. The capability of the parameterized microstructure to combine the dominant isotropic, hexagonal, tetragonal and orthorhombic symmetry classes in the proportions present in the cancellous bone is shown. Based on this finding, two optimization approaches are devised to find the geometrical parameters of the artificial microstructure that better mimics the elastic response of a target natural bone specimen: a Sequential Quadratic Programming algorithm that minimizes the norm of the difference between the elasticity matrices, and a Pattern Search algorithm that minimizes the difference between the symmetry class decompositions. The pattern search approach is found to produce the best results. The performance of the method is demonstrated via analyses for 146 bone samples.     

Elastic properties of cancellous bone: measurement by an ultrasonic technique

The high degree of porosity of cancellous bone makes elastic property measurement difficult by traditional mechanical testing methods. An ultrasonic technique is described with which mechanical properties of anisotropic, rigid, porous materials, such as cancellous bone, can be measured. The technique utilizes unique piezoelectric transducers operated in a continuous wave mode at a frequency of approximately 50 kHz. Both longitudinal and shear waves can be propagated and received with the transducers allowing both Young's moduli and shear moduli to be determined with the technique. A comparison between moduli measured with the ultrasonic technique and moduli measured with traditional mechanical testing shows the new method to be quite accurate in elastic property determination, (r2 = 0.935, Emech = 1.00E1dt + 23.3 MPa) (r2 = 0.656, Gmech = 1.08 Gult--3.3MPa).     

A new approach for assigning bone material properties from CT images into finite element models

Generation of subject-specific finite element (FE) models from computed tomography (CT) datasets is of significance for application of the FE analysis to bone structures. A great challenge that remains is the automatic assignment of bone material properties from CT Hounsfield Units into finite element models. This paper proposes a new assignment approach, in which material properties are directly assigned to each integration point. Instead of modifying the dataset of FE models, the proposed approach divides the assignment procedure into two steps: generating the data file of the image intensity of a bone in a MATLAB program and reading the file into ABAQUS via user subroutines. Its accuracy has been validated by assigning the density of a bone phantom into a FE model. The proposed approach has been applied to the FE model of a sheep tibia and its applicability tested on a variety of element types. The proposed assignment approach is simple and illustrative. It can be easily modified to fit users’ situations.     

Convergence behavior of high-resolution finite element models of trabecular bone

The convergence behavior of finite element models depends on the size of elements used, the element polynomial order, and on the complexity of the applied loads. For high-resolution models of trabecular bone, changes in architecture and density may also be important. The goal of this study was to investigate the influence of these factors on the convergence behavior of high-resolution models of trabecular bone. Two human vertebral and two bovine tibial trabecular bone specimens were modeled at four resolutions ranging from 20 to 80 microns and subjected to both compressive and shear loading. Results indicated that convergence behavior depended on both loading mode (axial versus shear) and volume fraction of the specimen. Compared to the 20 microns resolution, the differences in apparent Young's modulus at 40 microns resolution were less than 5 percent for all specimens, and for apparent shear modulus were less than 7 percent. By contrast, differences at 80 microns resolution in apparent modulus were up to 41 percent, depending on the specimen tested and loading mode. Overall, differences in apparent properties were always less than 10 percent when the ratio of mean trabecular thickness to element size was greater than four. Use of higher order elements did not improve the results. Tissue level parameters such as maximum principal strain did not converge. Tissue level strains converged when considered relative to a threshold value, but only if the strains were evaluated at Gauss points rather than element centroids. These findings indicate that good convergence can be obtained with this modeling technique, although element size should be chosen based on factors such as loading mode, mean trabecular thickness, and the particular output parameter of interest.     

Elastic constants of composites formed from PMMA bone cement and anisotropic bovine tibial cancellous bone

An ultrasonic pulse-transit time technique is used to determine the nine orthotropic engineering constants of 32 cement-cancellous bone composites as a function of volume fractions of bone ranging from 0.0 to 0.4. The composites are manufactured using well-aligned bovine cancellous bone from the proximal end of the tibia and low viscosity bone cement. Selected composites are also subjected to mechanical compression tests to compare with the ultrasonic results. There is excellent correlation between the dynamic or ultrasonically determined moduli and the static or mechanically determined moduli; the dynamic moduli are approximately twice the static moduli and this difference is thought to be due to the effect of strain rate. An orthotropic model is assumed requiring nine independent elastic constants to be determined. The dynamic Young's modulus in the direction of major trabecular alignment, E1, increases linearly from 4.9 to 10.4 GPa as bone volume fraction increases from 0 to 0.4; dynamic E2 and E3 values increase from 4.9 to 7 GPa as bone volume fractions increase from 0 to 0.4, with E2 being slightly higher than E3. The dynamic shear modulus, G12, increases from 1.8 to 3.0 GPa, and G31 and G23 increase slightly from 1.8 to 2.2 GPa as bone volume fractions increase from 0 to 0.4. The Poisson's ratios are more sensitive than the Young's moduli and shear moduli to experimental error in the velocity measurements. The mechanically tested modulus (static modulus) in the direction of major trabecular alignment, E1, increases with volume fraction of bone from 2.4 to 4.4 GPa as the bone volume fraction increases from 0 to 0.25; static E2 and E3 values are either equal to or lower than that of pure PMMA.     

Comparative analysis of bone remodelling models with respect to computerised tomography-based finite element models of bone

Subject-specific finite element models are an extensively used tool for the numerical analysis of the biomechanical behaviour of human bones. However, bone modelling is not an easy task due to the complex behaviour of bone tissue, involving non-homogeneous and anisotropic mechanical properties. Moreover, bone is a living tissue and therefore its microstructure and mechanical properties evolve with time in a known process called bone remodelling. This phenomenon has been widely studied, many being the numerical models that have been formulated to predict density distribution and its evolution in several bones. The aim of the present study is to assess the capability of a bone remodelling model to predict the bone density distribution of different types of human bone (femur, tibia and mandible) comparing the obtained results with the bone density estimated by means of computerised tomography. Good accuracy was observed for the bone remodelling predictions including the thickness of the cortical layer.     

Effect of microcomputed tomography voxel size on the finite element model accuracy for human cancellous bone

The level of structural detail that can be acquired and incorporated in a finite element (FE) analysis might greatly influence the results of microcomputed tomography (microCT)-based FE simulations, especially when relatively large bones, such as whole vertebrae, are of concern. We evaluated the effect of scanning and reconstruction voxel size on the microCT-based FE analyses of human cancellous tissue samples for fixed- and free-end boundary conditions using different combinations of scan/reconstruction voxel size. We found that the bone volume fraction (BV/TV) did not differ considerably between images scanned at 21 and 50 microm and reconstructed at 21, 50, or 110 microm (-0.5% to 7.8% change from the 21/21 microm case). For the images scanned and reconstructed at 110 microm, however, there was a large increase in BV/TV compared to the 21/21 microm case (58.7%). Fixed-end boundary conditions resulted in 1.8% [coefficient of variation (COV)] to 14.6% (E) difference from the free-end case. Dependence of model output parameters on scanning and reconstruction voxel size was similar between free- and fixed-end simulations. Up to 26%, 30%, 17.8%, and 32.3% difference in modulus (E), and average (VMExp), standard deviation (VMSD) and coefficient of variation (COV) of von Mises stresses, respectively, was observed between the 21/21 microm case and other scan/reconstruction combinations within the same (free or fixed) simulation group. Observed differences were largely attributable to scanning resolution, although reconstruction resolution also contributed significantly at the largest voxel sizes. All 21/21 microm results (taken as the gold standard) could be predicted from the 21/50 (r2adj= 0.91-0.99;p<0.001), 21/110 (r2adj =0.58-0.99;p<0.02) and 50/50 results (r2adj=0.61-0.97;p<0.02). While BV/TV, VMSD, and VMExp/sigma(z) from the 21/21 could be predicted by those from the 50/110 (r2adj =0.63-0.93;p<0.02) and 110/110 (r2adj =0.41-0.77;p<0.05) simulations as well, prediction of E, VMExp, and COV became marginally significant (0.04相似文献   

Biomechanical behaviour of cancellous bone on patellofemoral arthroplasty with Journey prosthesis: a finite element study

Isolated patellofemoral (PF) arthritis of the knee is a common cause of anterior knee pain and disability. Patellofemoral arthroplasty (PFA) is a bone conserving solution for patients with PF degeneration. Failure mechanisms of PFA include growing tibiofemoral arthritis and loosening of components. The implant loosening can be associated with bone resorption or fatigue-failure of bone by overload. This research work aims at determining the structural effects of the implantation of PF prosthesis Journey PFJ (Smith & Nephew, Inc., Memphis, TN, USA) on femoral cancellous bone. For this purpose, the finite element method is considered to perform computational simulations for different conditions, such as well-fixed and loosening scenarios. From the global results obtained, in the well-fixed scenario, a decrease in strain on cancellous bone was noticed, which can be related to bone resorption. In the loosening scenario, when the cement layer becomes inefficient, a significant increase in cancellous bone strain was observed, which can be associated with bone fatigue-failure.These strain changes suggest a weakness of the femur after PFA.     

Effects of thresholding techniques on microCT-based finite element models of trabecular bone

Microimaging based finite element analysis is widely used to predict the mechanical properties of trabecular bone. The choice of thresholding technique, a necessary step in converting grayscale images to finite element models, can significantly influence the predicted bone volume fraction and mechanical properties. Therefore, we investigated the effects of thresholding techniques on microcomputed tomography (micro-CT) based finite element models of trabecular bone. Three types of thresholding techniques were applied to 16-bit micro-CT images of trabecular bone to create three different models per specimen. Bone volume fractions and apparent moduli were predicted and compared to experimental results. In addition, trabecular tissue mechanical parameters and morphological parameters were compared among different models. Our findings suggest that predictions of apparent mechanical properties and structural properties agree well with experimental measurements regardless of the choice of thresholding methods or the format of micro-CT images.     

Finite element modelling of polymethylmethacrylate flow through cancellous bone

A mathematical model based on the Finite Element Method is developed to simulate the non-linear flow of acrylic bone cement through cancellous bone. The cancellous bone bed is modelled as a bed of parallel capillaries filled with equal spaced toroidal trabeculae. By manipulating the relative size of the torus and the capillary, the flow within bone of varying porosity is simulated. An apparent permeability based on the volume weighted average viscosity and Darcy's law is developed to describe the flow of the acrylic through the cancellous bone bed. The model predicts a cancellous bone permeability of 5.6 x 10(-9)-8.3 x 10(-9) m2 for linear flow. The non-linear behavior of the acrylic cement results in an increase of apparent permeability when compared to the permeability computed for linear flow. Estimates of penetration are achieved by running the model in a quasi-steady state fashion with pressure applied over a fixed time increment. Close agreement is shown between model predictions of penetration depth and experimental results available in the literature.     

Using 3D finite element models verified the importance of callus material and microstructure in biomechanics restoration during bone defect repair

Background: There is lack of further observations on the microstructure and material property of callus during bone defect healing and the relationships between callus properties and the mechanical strength. Methods: Femur bone defect model was created in rabbits and harvested CT data to reconstruct finite element models at 1 and 2 months. Three types of assumed finite element models were compared to study the callus properties, which assumed the material elastic property as heterogeneous (R-model), homogenous (H-model) or did not change from 1 to 2 months (U-model). Results: The apparent elastic moduli increased at 2 months (from 355.58 ± 132.67 to 1139.30 ± 967.43 MPa) in R-models. But there was no significant difference in apparent elastic moduli between R-models (355.58 ± 132.67 and 1139.30 ± 967.43 MPa) and H-models (344.79 ± 138.73 and 1001.52 ± 692.12 MPa) in 1 and 2 months. A significant difference of apparent elastic moduli was found between the R-model (1139.30 ± 967.43 MPa) and U-model group (207.15 ± 64.60 MPa) in 2 months. Conclusions: This study showed that the callus structure stability remodeled overtime to achieve a more effective structure, while the material quality of callus tissue is a very important factor for callus strength. At the meantime, this study showed an evidence that the material heterogeneity maybe not as important as it is in bone fracture model.     

Tensile and compressive properties of cancellous bone

The relationship between the mechanical properties of trabecular bone in tension and compression was investigated by non-destructive testing of the same specimens in tension and compression, followed by random allocation to a destructive test in either tension or compression. There was no difference between Young's modulus in tension and compression, and there was a strong positive correlation between the values (R = 0.97). Strength, ultimate strain and work to failure was significantly higher in tensile testing than in compressive testing.     

Scaffolds derived from cancellous bovine bone support mesenchymal stem cells' maintenance and growth

Since bone defects can lead to various disabilities, in recent years, many increasing attempts have been made in bone tissue engineering. In this regard, scaffolds have attracted a lot of attention as three dimensional substrates for cell attachment which improve successful tissue engineering. The aim of the present study was to provide an interconnected porous scaffold to facilitate cell infiltration. To do so, cancellous bone from bovine femur was dissected in fragments and decellularized by physicochemical methods, including snap freeze/thaw, rinsing in hot water and treatment with different solutions of sodium dodecyl sulfate (SDS). Histological analysis and 4′,6-diamidino-2-phenylindole staining revealed that the best results were obtained after treatment with 2.5%, 5%, and 8% SDS for 8, 3, or 1 h respectively, which significantly removed bone cells with intact trabeculae geometry. Further characterization of decellularized scaffolds by the compression tests also revealed no significant difference between elastic modulus values of the three different SDS treatments. Moreover, studying the ratio of bone trabeculae to bone surfaces (BT/BS) as assessed by Clemex vision software 3.5 showed that treatment with 2.5% SDS for 8 h resulted in a BT/BS score in the range of native bone and therefore this treatment was used for further experiments. Histological studies and scanning electron microscopy revealed rat mesenchymal stem cells integration, adhesion, and maintenance during the 2 and 7 d of culture in vitro. In conclusion, the present results support the effective role of SDS in cancellous bovine bone decellularization and also propensity of treated samples in providing a suitable three-dimentional environment to support the maintenance and growth of mesenchymal stem cells.     

Finite element prediction of fatigue damage growth in cancellous bone

Cyclic stresses applied to bones generate fatigue damage that affects the bone stiffness and its elastic modulus. This paper proposes a finite element model for the prediction of fatigue damage accumulation and failure in cancellous bone at continuum scale. The model is based on continuum damage mechanics and incorporates crack closure effects in compression. The propagation of the cracks is completely simulated throughout the damaged area. In this case, the stiffness of the broken element is reduced by 98% to ensure no stress-carrying capacities of completely damaged elements. Once a crack is initiated, the propagation direction is simulated by the propagation of the broken elements of the mesh. The proposed model suggests that damage evolves over a real physical time variable (cycles). In order to reduce the computation time, the integration of the damage growth rate is based on the cycle blocks approach. In this approach, the real number of cycles is reduced (divided) into equivalent blocks of cycles. Damage accumulation is computed over the cycle blocks and then extrapolated over the corresponding real cycles. The results show a clear difference between local tensile and compressive stresses on damage accumulation. Incorporating stiffness reduction also produces a redistribution of the peak stresses in the damaged region, which results in a delay in damage fracture.     

Biologic properties of mesenchymal stem cells derived from bone marrow and adipose tissue

The biologic characteristics of mesenchymal stem cells (MSCs) isolated from two distinct tissues, bone marrow and adipose tissue were evaluated in these studies. MSCs derived from human and non-human primate (rhesus monkey) tissue sources were compared. The data indicate that MSCs isolated from rhesus bone marrow (rBMSCs) and human adipose tissue (hASCs) had more similar biologic properties than MSCs of rhesus adipose tissue (rASCs) and human bone marrow MSCs (hBMSCs). Analyses of in vitro growth kinetics revealed shorter doubling time for rBMSCs and hASCs. rBMSCs and hASCs underwent significantly more population doublings than the other MSCs. MSCs from all sources showed a marked decrease in telomerase activity over extended culture; however, they maintained their mean telomere length. All of the MSCs expressed embryonic stem cell markers, Oct-4, Rex-1, and Sox-2 for at least 10 passages. Early populations of MSCs types showed similar multilineage differentiation capability. However, only the rBMSCs and hASCs retain greater differentiation efficiency at higher passages. Overall in vitro characterization of MSCs from these two species and tissue sources revealed a high level of common biologic properties. However, the results demonstrate clear biologic distinctions, as well.