Investigation of the mechanical interaction of the trabecular core with an external shell using rapid prototype and finite element models

Are we crying Wolff? 3D printed replicas of trabecular bone structure demonstrate higher stiffness and strength during off-axis loading

Roux's principle of bone functional adaptation postulates that bone tissue, and particularly trabecular bone tissue, responds to mechanical stimuli by adjusting (modeling) its architecture accordingly. Hence, it predicts that the new modeled trabecular structure is mechanically improved (stiffer and stronger) in line with the habitual in vivo loading direction. While previous studies found indirect evidence to support this theory, direct support was so far unattainable. This is attributed to the fact that each trabecular bone is unique, and that trabecular bone tissue tends to be damaged during mechanical testing. Consequently, a unique modeled trabecular structure can be mechanically tested only along one direction and a comparison to other directions for that specific structure is impossible. To address this issue, we have 3D printed 10 replicas of a trabecular structure from a sheep talus cropped along the 3 principal axes of the bone and in line with the principal direction of loading (denoted on-axis model). Next, we have rotated the same cropped trabecular structure in increments of 10° up to 90° to the bone principal direction of loading (denoted off-axis models) and printed 10 replicas of each off-axis model. Finally, all on-axis and off-axis 3D printed replicas were loaded in compression until failure and trabecular structure stiffness and strength were calculated. Contrary to our prediction, and conflicting with Roux's principle of bone functional adaptation, we found that a trabecular structure loaded off-axis tended to have higher stiffness and strength values when compared to the same trabecular structure loaded on-axis. These unexpected results may not disprove Roux's principle of bone functional adaptation, but they do imply that trabecular bone adaptation may serve additional purposes than simply optimizing bone structure to one principal loading scenario and this suggests that we still don't fully understand bone modeling in its entirety.

Modeling the Mechanical Consequences of Age-Related Trabecular Bone Loss by XFEM Simulation

The elderly are more likely to suffer from fracture because of age-related trabecular bone loss. Different bone loss locations and patterns have different effects on bone mechanical properties. Extended finite element method (XFEM) can simulate fracture process and was suited to investigate the effects of bone loss on trabecular bone. Age-related bone loss is indicated by trabecular thinning and loss and may occur at low-strain locations or other random sites. Accordingly, several ideal normal and aged trabecular bone models were created based on different bone loss locations and patterns; then, fracture processes from crack initiation to complete failure of these models were observed by XFEM; finally, the effects of different locations and patterns on trabecular bone were compared. Results indicated that bone loss occurring at low-strain locations was more detrimental to trabecular bone than that occurring at other random sites; meanwhile, the decrease in bone strength caused by trabecular loss was higher than that caused by trabecular thinning, and the effects of vertical trabecular loss on mechanical properties were more severe than horizontal trabecular loss. This study provided a numerical method to simulate trabecular bone fracture and distinguished different effects of the possible occurrence of bone loss locations and patterns on trabecular bone.

A Computational Approximation to Model Variation in Cancellous Bone Screw Pull-Out

Cancellous bone screws are used to achieve good pull-out characteristics when connected to cancellous bone. This study examines some screw characteristics, including pitch and inner diameter, using a model of cancellous bone with a range of bone apparent densities (ADs). This was achieved using bone geometry based on microCT-scanned cancellous bone and converted into a geometric model using mimics® software. The finite element (FE) models were produced in ansys®. The calculated reaction force for pull-out of 0.2 mm shows the influence of design parameters. Change in the proximal half angle increased the stiffness by about 15% in line with the experimental findings of others. An increase in pull-out reaction force with an increase in bone AD was also observed. However, when a particular screw geometry in lower AD bone was modeled and then rotated through 180 deg on plan, a significant reduction in reaction force was noted. Further models with screws of similar geometry in the same location showed similar reductions in reaction force and hence pull-out stiffness. Examination of the geometry of the bone/screw interface indicates that in certain positions there is little cancellous bone to support the implant—leading to low pull-out reaction forces, which is very difficult to predict. The study also examined the effect of increasing the bone stiffness adjacent to the implant, and concluded that, even in bone of low AD, increases in pull-out stiffness might be achieved.

A parametric finite element analysis of the compacted bone–cement interface following balloon kyphoplasty

Treating fractures of the spine is a major challenge for the medical community with an estimated 1.4 million fractures per annum worldwide. While a considerable volume of study exists on the biomechanical implications of balloon kyphoplasty, which is used to treat these fractures, the influence of the compacted bone–cement region properties on stress distribution within the vertebral body remains unknown. The following article describes a novel method for modelling this compacted bone–cement region using a geometry-based approach in conjunction with the knowledge of the bone volume fractions for the native and compacted bone regions. Three variables for the compacted region were examined, as follows: (1) compacted thickness, (2) compacted region Young’s modulus and (3) friction coefficient. Results from the model indicate that the properties of the compacted bone–cement region can affect stresses in the cortical bone and cement by up to +28% and −40%, respectively. These findings demonstrate the need for further investigation into the effects of the compacted bone–cement interface using computational and experimental methods on multi-segment models.

Finite element micro-modelling of a human ankle bone reveals the importance of the trabecular network to mechanical performance: new methods for the generation and comparison of 3D models

Most modelling of whole bones does not incorporate trabecular geometry and treats bone as a solid non-porous structure. Some studies have modelled trabecular networks in isolation. One study has modelled the performance of whole human bones incorporating trabeculae, although this required considerable computer resources and purpose-written code. The difference between mechanical behaviour in models that incorporate trabecular geometry and non-porous models has not been explored. The ability to easily model trabecular networks may shed light on the mechanical consequences of bone loss in osteoporosis and remodelling after implant insertion. Here we present a Finite Element Analysis (FEA) of a human ankle bone that includes trabecular network geometry. We compare results from this model with results from non-porous models and introduce protocols achievable on desktop computers using widely available softwares. Our findings show that models including trabecular geometry are considerably stiffer than non-porous whole bone models wherein the non-cortical component has the same mass as the trabecular network, suggesting inclusion of trabecular geometry is desirable. We further present new methods for the construction and analysis of 3D models permitting: (1) construction of multi-property, non-porous models wherein cortical layer thickness can be manipulated; (2) maintenance of the same triangle network for the outer cortical bone surface in both 3D reconstruction and non-porous models allowing exact replication of load and restraint cases; and (3) creation of an internal landmark point grid allowing direct comparison between 3D FE Models (FEMs).

Finite element analysis of a micromechanical model of bone and a new 3D approach to validation

Finite Element Analysis (FEA) is now widely used to analyse the mechanical behaviour of bone structures. Ideally, simulations are validated against experimental data. To date, validation of Finite Element Models (FEMs) has been 2 Dimensional (2D) only, being based on comparison with surface-mounted strain gauge readings. In this study we present a novel 3-Dimensional (3D) approach to validation that allows comparison of modelled with experimental results between any two points in 3D space throughout the structure, providing magnitude and direction data for comparison, internally and externally. Specifically, we validate a FEM of a rat tibia, including trabecular network geometry, using a material testing stage housed within a microCT scanner. We further apply novel landmark based morphometric approaches to more effectively compare modelled and experimental results. 542 landmark points on the cortical and trabecular bone surfaces of the model were selected and validated in 3D against experimental data. This approach may hold considerable potential in fields wherein a better understanding of the mechanical behaviour of trabecular networks is important, e.g., the studies of osteoporosis and trabecular loss after orthopaedic implant insertion.

An approximate model for cancellous bone screw fixation

This paper presents a finite element (FE) model to identify parameters that affect the performance of an improved cancellous bone screw fixation technique, and hence potentially improve fracture treatment. In cancellous bone of low apparent density, it can be difficult to achieve adequate screw fixation and hence provide stable fracture fixation that enables bone healing. Data from predictive FE models indicate that cements can have a significant potential to improve screw holding power in cancellous bone. These FE models are used to demonstrate the key parameters that determine pull-out strength in a variety of screw, bone and cement set-ups, and to compare the effectiveness of different configurations. The paper concludes that significant advantages, up to an order of magnitude, in screw pull-out strength in cancellous bone might be gained by the appropriate use of a currently approved calcium phosphate cement.

Comparison of trabecular bone behavior in core and whole bone samples using high-resolution modeling of a vertebral body

Computational analysis of trabecular bone normally involves the modeling of (experimental tests of) cored samples. However, the lack of constraint on the sides of the extracted trabecular bone samples limits the information that can be inferred regarding true in situ behavior. Here, the element-by-element voxel-based finite element method was applied via, a custom-written software suite (FEEBE), to a 72 microm resolution model of an ovine vertebra. The difference between the apparent modulus of eight concentric core cylinders when modeled as part of the whole bone (containing 84 x 10(6) degrees of freedom) and independent of the whole bone was investigated. The results showed that cored trabecular bone apparent modulus depended significantly on the core diameter when modeled as an extracted core (r (2) = 0.975) and as part of a whole bone (r (2) = 0.986). The cause of this result was separated into the side-artifact effect and bone volume fraction (BV/TV) effect. For the independently modeled cores, the apparent modulus of an inner core region of interest varied with increasing thickness of the outer annulus. This was attributed to the side-artifact effect, given that the BV/TV of the core region was constant. Within the whole trabecular structure, the side artifact was eliminated as the entire bone structure was modeled. However, a BV/TV effect influenced the apparent modulus depending on the size of the core selected for determining apparent modulus. Changing the size of the core varied the overall BV/TV of the core, and this significantly (r (2) = 0.999) influences the apparent modulus. Therefore, determining a 'true' apparent modulus for trabecular bone was not achievable. The independently modeled cores consistently under-predict the in vivo apparent modulus. It is recommended that if a 'true' apparent modulus is required, the BV/TV at which it is required needs to be first determined. Apparent modeling of entire bones at microscale resolution allowed regions of low and high tissue strains to be identified, consistent with patterns of trabecular bone remodeling and resorption reported in literature. The basivertebral vein cavity underwent the highest strains within the entire vertebral body, suggesting that failure might initiate here, despite containing visibly thicker struts and plate trabeculae. Although computationally expensive, analysis of the entire vertebral body provided a full picture of in situ trabecular bone deformation.