The multiscale meso-mechanics model of viscoelastic cortical bone

Biomechanics and Modeling in Mechanobiology - Cortical bone is a complex hierarchical structure consisting of biological fiber composites with transversely isotropic constituents, whose...     

Specifications for machining the bovine cortical bone in relation to its microstructure

Until date, many devices have been developed for cutting human bones during orthopedic surgeries. However, bones are anisotropic material, and their machining characteristics depend on the tool feed direction. In this study, microcutting of the bovine cortical bone is performed and its structure observed under a microscope. Furthermore, the formation of cutting chips and measurement of the cutting force during bone machining are dynamically observed while considering the anisotropy of bone tissue. In particular, the fracture of secondary osteons and crack propagation in bones are observed and analyzed. The results indicate that when the cut depth exceeds 20 μm and is greater than the interval of concentric lamellae, cracks are formed together with chips. A new method for bone machining is proposed. This method is based on the characteristics of crack propagation in bones and is expected to produce low mechanical stress and realize highly efficient and precise machining of living tissues such as bones.     

Application of the multiscale FEM to the modeling of cancellous bone

This paper considers the application of multiscale finite element method (FEM) to the modeling of cancellous bone as an alternative for Biot’s model, the main intention of which is to decrease the extent of the necessary laboratory tests. At the beginning, the paper gives a brief explanation of the multiscale concept and thereafter focuses on the modeling of the representative volume element and on the calculation of the effective material parameters, including an analysis of their change with respect to increasing porosity. The latter part of the paper concentrates on the macroscopic calculations, which is illustrated by the simulation of ultrasonic testing and a study of the attenuation dependency on material parameters and excitation frequency. The results endorse conclusions drawn from the experiments: increasing excitation frequency and material density cause increasing attenuation.     

Micromechanics modeling of Haversian cortical bone properties.

A finite-element micromechanics model for Haversian cortical bone tissue has been developed and studied. The model is an extension of two-dimensional micromechanics techniques for fiber-reinforced composite materials. Haversian systems, or secondary osteons, are considered to be the fiber component, and interstitial lamellar bone the matrix material. The cement line is included as an 'interphase' component along the fiber/matrix interface. The model assumes a regular repeatable spacing of the longitudinally aligned continuous fibers and is, therefore, restricted to approximating Haversian cortical bone in its present form. Haversian porosity is modeled explicitly by incorporating a hollow fiber to represent the Haversian canal. Solutions have been obtained by applying uniform macroscopic stresses to the boundaries of the repeating unit cell model. Macroscopic mechanical property predictions correspond reasonably well with the experimental data for cortical bone, but are necessarily dependent on the input properties for each constituent, which are not well established. The predicted variation in the elastic modulus with porosity is not as sensitive as that observed experimentally. Stresses within the constituents can also be modeled with this method and are demonstrated to deviate from the macroscopic applied stress levels.     

Dependence of ultrasonic attenuation on bone mass and microstructure in bovine cortical bone

The development of the axial transmission technique now enables in vivo evaluation of cortical bone quality, which plays an important role in bone fragility. Cortical bone is a complex multiscale material, which may be made of different types of microstructure. The interaction between ultrasound and cortical bone remains unclear and most studies have been confined to wave speed analysis. The first aim of this study is to investigate the dependence of the frequency-dependent attenuation on the type of bone microstructure. The second goal is to determine whether broadband ultrasonic attenuation (BUA) is related to volumetric bone mineral density (vBMD) and mass density. Parallelepipedic samples of bovine cortical bone were cut from three specimens and tested in the axial, radial and tangential directions using an ultrasonic transmission device. BUA was evaluated over a 1-MHz wide bandwidth around 4MHz. In addition, the microstructure of each sample was determined using an optical microscope. BUA values measured in porotic microstructure are significantly higher than in Haversian microstructure. The lowest BUA values are obtained for plexiform microstructure. For all structures, BUA in the axial direction is significantly smaller than in the radial and tangential directions. Moreover, BUA is correlated with both vBMD and density (determination coefficient (R2) equal to 0.44 and 0.65, respectively, in the axial direction). BUA variations can be explained by scattering and viscoelastic mechanisms. This study suggests that BUA measurements have the potential to discriminate among different cortical bone microstructures in addition to providing material properties.     

Nonlinear multiscale analysis of coronary atherosclerotic vulnerable plaque artery: fluid-structural modeling with micromechanics

Biomechanics and Modeling in Mechanobiology - A unique three-dimensional (3D) computational multiscale modeling approach is proposed to investigate the influence of presence of microcalcification...     

Numerical simulation of streaming potentials due to deformation-induced hierarchical flows in cortical bone

Bone is a very dynamic tissue capable of modiA,fing its composition, microstructure, and overall geometry in response to the changing biomechanical needs. Streaming potential has been hypothesized as a mechanotransduction mechanism that may allow osteocytes to sense their biomechanical environment. A correct understanding of the mechanism for streaming potential will illuminate our understanding of bone remodeling, such as the remodeling associated with exercise hypertrophy, disuse atrophy, and the bone remodeling arounid implants. In the current research, a numerical model based on the finite element discretization is proposed to simulate the fluid flows through the complicated hierarchical flow system and to calculate the concomitant stress generated potential (SGP) as a result of applied mechanical loading. The lacunae-canaliculi and the matrix microporosity are modeled together as discrete one-dimensional flow channels superposed in a biphasic poroelastic matrix. The cusplike electric potential distribution surrounding the Haversian canal that was experimentallv observed and reported in the literature earlier was successfully reproduced by the current numerical calculation.     

A new numerical concept for modeling hydroxyapatite in human cortical bone

This research presents a new modelling procedure which allows the computation of the physical properties of the human cortical bone, considered as a strongly heterogeneous medium consisting of bony architecture and the physical properties of the two basic components: the collagen and the hydroxyapatite (Hap). The numerical simulations are based on the homogenisation theory, however, since the size of the Hap crystals are small compared to the size of a collagen stick, a new entity (the elementary volume of mineral content (EVMC)) is defined at the nanoscopic scale. This model permits the testing of all the possible structural configurations that may be present and suggests that the anisotropy of the bone is not only induced by the haversian structure but by the properties of the Hap crystals and their special organisation.     

A new numerical concept for modeling hydroxyapatite in human cortical bone

This research presents a new modelling procedure which allows the computation of the physical properties of the human cortical bone, considered as a strongly heterogeneous medium consisting of bony architecture and the physical properties of the two basic components: the collagen and the hydroxyapatite (Hap). The numerical simulations are based on the homogenisation theory, however, since the size of the Hap crystals are small compared to the size of a collagen stick, a new entity (the elementary volume of mineral content (EVMC)) is defined at the nanoscopic scale. This model permits the testing of all the possible structural configurations that may be present and suggests that the anisotropy of the bone is not only induced by the haversian structure but by the properties of the Hap crystals and their special organisation.     

Tensile damage and its effects on cortical bone

Plexiform bovine bone samples are repeatedly loaded in tension along their longitudinal axis. In order to induce damage in the bone tissue, bone samples are loaded past their yield point. Half of the bone samples from the damaged group were stored in saline to allow for viscoelastic recovery while the others were decalcified. Tensile tests were conducted on these samples to characterize the effects of damage on the mechanical behavior of the organic matrix (decalcified samples) as well as on bone tissue (stored in saline). The ultimate strain of the damaged decalcified bone is 29% higher compared to that of non-damaged decalcified (control) bone. The ultimate stresses as well as the elastic moduli are similar in both decalcified groups. This phenomenon is also observed in other collagenous tissue (tendon and ligament). This may suggest that damage in bone is caused by shear failure of the organic matrix; transverse separation of the collagen molecules or microfibrils from each other. In contrast, there is a trend towards lowered ultimate strains in damaged bone, which is soaked in saline, with respect to control bone samples (not damaged). The damaged bone tissue exhibits a bi-linear behavior in contrast to the mechanical behavior of non-damaged bone. The initial elastic modulus (below 55 MPa) and ultimate strength of damaged bone are similar to that in non-damaged bone.     

Subject-specific multiscale modeling of aortic valve biomechanics

A Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range (\(0.36-0.87\)) with a strong correlation with the organ length-scale radial strain (\(R^{2}= 0.95\)); conversely, circumferential cell strains spanned a very narrow range (\(0.75-0.88\)) showing no correlation with the circumferential strain at the organ level (\(R^{2}= 0.02\)). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.

     

Cohesive finite element modeling of age-related toughness loss in human cortical bone

Although the age-related loss of bone quality has been implicated in bone fragility, a mechanistic understanding of the relationship is necessary for developing diagnostic and treatment modalities in the elderly population at risk of fracture. In this study, a finite element based cohesive zone model is developed and applied to human cortical bone in order to capture the experimentally shown rising crack growth behavior and age-related loss of bone toughness. The cohesive model developed here is based on a traction–crack opening displacement relationship representing the fracture processes in the vicinity of a propagating crack. The traction–displacement curve, defining the cohesive model, is composed of ascending and descending branches that incorporate material softening and nonlinearity. The results obtained indicate that, in contrast to initiation toughness, the finite element simulations of crack growth in compact tension (CT) specimens successfully capture the rising R-curve (propagation toughness) behavior and the age-related loss of bone toughness. In close correspondence with the experimentally observed decrease of 14–15% per decade, the finite element simulation results show a decrease of 13% in the R-curve slope per decade. The success of the simulations is a result of the ability of cohesive models to capture and predict the parameters related to bone fracture by representing the physical processes occurring in the vicinity of a propagating crack. These results illustrate that fracture mechanisms in the process zone control bone toughness and any modification to these would cause age-related toughness loss.     

Impact of the porous microstructure on the overall elastic properties of the osteonal cortical bone

Mechanical properties of bones are largely determined by their microstructure. The latter comprises a large number of diverse pores. The present paper analyzes a connection between structure of the porous space of the osteonal cortical bone and bone's overall anisotropic elastic moduli. The analysis is based on recent developments in the theory of porous materials that predict the anisotropic effective moduli of porous solids in terms of pores' shapes, orientations and densities. Bone's microstructure is modeled using available micrographs. The calculated anisotropic elastic constants for porous cortical bone are, mostly, in agreement with available experimental data. The influence of each of the pore types on the overall moduli is examined. The results of the analysis can also be used to estimate the extent of mineralization (hydroxyapatite content) if the overall porosity and the effective moduli are known and, vice versa, to estimate porosity from the measured moduli and the extent of mineralization.     

Impact loading intensifies cortical bone (re)modeling and alters longitudinal bone growth of pubertal rats

Biomechanics and Modeling in Mechanobiology - Physical exercise is important for musculoskeletal development during puberty, which builds bone mass foundation for later in life. However, strenuous...     

Micromechanical modeling of elastic properties of cortical bone accounting for anisotropy of dense tissue

The paper analyzes the connection between microstructure of the osteonal cortical bone and its overall elastic properties. The existing models either neglect anisotropy of the dense tissue or simplify cortical bone microstructure (accounting for Haversian canals only). These simplifications (related mostly to insufficient mathematical apparatus) complicate quantitative analysis of the effect of microstructural changes – produced by age, microgravity, or some diseases – on the overall mechanical performance of cortical bone. The present analysis fills this gap; it accounts for anisotropy of the dense tissue and uses realistic model of the porous microstructure. The approach is based on recent results of Sevostianov et al. (2005) and Saadat et al. (2012) on inhomogeneities in a transversely-isotropic material. Bone?s microstructure is modeled according to books of Martin and Burr (1989), Currey (2002), and Fung (1993) and includes four main families of pores. The calculated elastic constants for porous cortical bone are in agreement with available experimental data. The influence of each of the pore types on the overall moduli is examined.     

Multiscale modeling of cancellous bone considering full coupling of mechanical,electric and magnetic effects

Modeling of cancellous bone has important applications in the detection and treatment of fatigue fractures and diseases like osteoporosis. In this paper, we present a fully coupled multiscale approach considering mechanical, electric and magnetic effects by using the multiscale finite element method and a two-phase material model on the microscale. We show numerical results for both scales, including calculations for a femur bone, comparing a healthy bone to ones affected by different stages of osteoporosis. Here, the magnetic field strength resulting from a small mechanical impact decreases drastically for later stages of the disease, confirming experimental research.