Elastic response of trabecular bone under compression calculated using the firm and floppy boundary lattice element method.

Micro-Finite Element analysis (µFEA) has become widely used in biomechanical research as a reliable tool for the prediction of bone mechanical properties within its microstructure such as apparent elastic modulus and strength. However, this method requires substantial computational resources and processing time. Here, we propose a computationally efficient alternative to FEA that can provide an accurate estimation of bone trabecular mechanical properties in a fast and quantitative way. A lattice element method (LEM) framework based on the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) open-source software package is employed to calculate the elastic response of trabecular bone cores. A novel procedure to handle pore-material boundaries is presented, referred to as the Firm and Floppy Boundary LEM (FFB-LEM). Our FFB-LEM calculations are compared to voxel- and geometry-based FEA benchmarks incorporating bovine and human trabecular bone cores imaged by micro Computed Tomography (µCT). Using 14 computer cores, the apparent elastic modulus calculation of a trabecular bone core from a µCT-based input with FFB-LEM required about 15 min, including conversion of the µCT data into a LAMMPS input file. In contrast, the FEA calculations on the same system including the mesh generation, required approximately 30 and 50 min for voxel- and geometry-based FEA, respectively. There were no statistically significant differences between FFB-LEM and voxel- or geometry-based FEA apparent elastic moduli (+24.3% or +7.41%, and +0.630% or -5.29% differences for bovine and human samples, respectively).

Mechanical loading of ex vivo bovine trabecular bone in 3D printed bioreactor chambers.

Previous ex vivo bone culture methods have successfully implemented polycarbonate (PC) bioreactors to investigate bone adaptation to mechanical load; however, they are difficult to fabricate and have been limited to a 5 mm maximum specimen height. The objective of this study was to validate a custom-made 3D printed MED610TM bioreactor system that addresses the limitations of the PC bioreactor and assess its efficacy in ex vivo bone culture. Twenty-three viable trabecular bone cores (10 mm height by 10 mm diameter) from an 18-month-old bovine sternum were cultured in MED610TM bioreactors with culture medium at 37 °C and 5% CO2 for 21-days. Bone cores were ranked based on their day 0 apparent elastic modulus (Eapp) and evenly separated into a "Load" group (n = 12) and a control group (n = 11). The Load group was loaded five times per week with a sinusoidal strain waveform between -1000 and -5000 µÎµ for 120 cycles at 2 Hz. Eapp was assessed on day 0, 8, and 21 using quasi-static tests with a -4000 µÎµ applied strain. Over 21-days, the Eapp of Load group samples tended to increase by more than double the control group (53.4% versus 20.9%) and no visual culture contamination was observed. This study demonstrated that bone organ culture in 3D printed MED610TM bioreactors replicated Eapp trends found in previous studies with PC bioreactors. However, further studies are warranted with a larger sample size to increase statistical power and histology to assess cell viability and bone mineral apposition rate.