Influence of the mineral staggering on the elastic properties of the mineralized collagen fibril in lamellar bone.

In this work, a three-dimensional finite element model of the staggered distribution of the mineral within the mineralized collagen fibril has been developed to characterize the lamellar bone elastic behavior at the sub-micro length scale. Minerals have been assumed to be embedded in a collagen matrix, and different degrees of mineralization have been considered allowing the growth of platelet-shaped minerals both in the axial and the transverse directions of the fibril, through the variation of the lateral space between platelets. We provide numerical values and trends for all the elastic constants of the mineralized collagen fibril as a function of the volume fraction of mineral. In our results, we verify the high influence of the mineral overlapping on the mechanical response of the fibril and we highlight that the lateral distance between crystals is relevant to the mechanical behavior of the fibril and not only the mineral overlapping in the axial direction.

Numerical modelling of the mechanical behaviour of an osteon with microcracks.

In this work, we present two strategies for the numerical modelling of microcracks and damage within an osteon. A numerical model of a single osteon under compressive diametral load is developed, including lamellae organized concentrically around the haversian canal and the presence of lacunae. Elastic properties have been estimated from micromechanical models that consider the mineralized collagen fibrils reinforced with hydroxyapatite crystals and the dominating orientation of the fibrils in each lamella. Microcracks are simulated through the node release technique, enabling propagation along the lamellae interfaces by application of failure criteria initially conceived for composite materials, in particular the Brewer and Lagacé criterion for delamination. A second approach is also presented, which is based on the progressive degradation of the stiffness at the element level as the damage increases. Both strategies are discussed, showing a good agreement with experimental evidence reported by other authors. It is concluded that interlaminar shear stresses are the main cause of failure of an osteon under compressive diametral load.

Homogenized stiffness matrices for mineralized collagen fibrils and lamellar bone using unit cell finite element models.

Mineralized collagen fibrils have been usually analyzed like a two-phase composite material where crystals are considered as platelets that constitute the reinforcement phase. Different models have been used to describe the elastic behavior of the material. In this work, it is shown that when Halpin-Tsai equations are applied to estimate elastic constants from typical constituent properties, not all crystal dimensions yield a model that satisfy thermodynamic restrictions. We provide the ranges of platelet dimensions that lead to positive definite stiffness matrices. On the other hand, a finite element model of a mineralized collagen fibril unit cell under periodic boundary conditions is analyzed. By applying six canonical load cases, homogenized stiffness matrices are numerically calculated. Results show a monoclinic behavior of the mineralized collagen fibril. In addition, a 5-layer lamellar structure is also considered where crystals rotate in adjacent layers of a lamella. The stiffness matrix of each layer is calculated applying Lekhnitskii transformations, and a new finite element model under periodic boundary conditions is analyzed to calculate the homogenized 3D anisotropic stiffness matrix of a unit cell of lamellar bone. Results are compared with the rule-of-mixtures showing in general good agreement.