Accuracy of beam theory for estimating bone tissue modulus and yield stress from 3-point bending tests on rat femora

Elucidation of the radius and ulna fracture mechanisms in toy poodle dogs using finite element analysis

Fractures occurring in the distal radius and ulna of toy breed dogs pose distinctive challenges for veterinary practitioners, requiring specialized treatment approaches primarily based on anatomical features. Finite Element Analysis (FEA) was applied to conduct numerical experiments to determine stress distribution across the bone. This methodology offers an alternative substitute for directly investigating these phenomena in living dog experiments, which could present ethical obstacles. A three-dimensional bone model of the metacarpal, carpal, radius, ulna, and humerus was reconstructed from Computed Tomography (CT) images of the toy poodle and dachshund forelimb. The model was designed to simulate the jumping and landing conditions from a vertical distance of 40 cm to the ground within a limited timeframe. The investigation revealed considerable variations in stress distribution patterns between the radius and ulna of toy poodles and dachshunds, indicating notably elevated stress levels in toy poodles compared to dachshunds. In static and dynamic stress analysis, toy poodles exhibit peak stress levels at the distal radius and ulna. The Von Mises stresses for toy poodles reach 90.07 MPa (static) and 1,090.75 MPa (dynamic) at the radius and 1,677.97 MPa (static) and 1,047.98 MPa (dynamic) at the ulna. Conversely, dachshunds demonstrate lower stress levels for 5.39 MPa (static) and 231.79 MPa (dynamic) at the radius and 390.56 MPa (static) and 513.28 MPa (dynamic) at the ulna. The findings offer valuable insights for modified treatment approaches in managing fractures in toy breed dogs, optimizing care and outcomes.

Method for Evaluating Cortical Bone Young's Modulus: Numerical Twin Reconstruction, Finite Element Calculation, and Microstructure Analysis

The determination of bone mechanical properties remains crucial, especially to feed up numerical models. An original methodology of inverse analysis has been developed to determine the longitudinal elastic modulus of femoral cortical bone. The method is based on a numerical twin of a specific three-point bending test. It has been designed to be reproducible on each test result. In addition, the biofidelity of the geometric acquisition method has been quantified. As the assessment is performed at the scale of a bone shaft segment, the Young's modulus values obtained (between 9518.29 MPa and 14181.15 MPa) are considered average values for the whole tissue, highlighting some intersubject variability. The material microstructure has also been studied through histological analysis, and bone-to-bone comparisons highlighted discrepancies in quadrants microstructures. Furthermore, significant intrasubject variability exists since differences between the bone's medial-lateral and anterior-posterior quadrants have been observed. Thus, the study of microstructures can largely explain the differences between the elastic modulus values obtained. However, a more in-depth study of bone mineral density would also be necessary and would provide some additional information. This study is currently being setup, alongside an investigation of the local variations of the elastic modulus.

Contribution to the 3R Principle: Description of a Specimen-Specific Finite Element Model Simulating 3-Point-Bending Tests in Mouse Tibiae

Bone mechanical properties are classically determined by biomechanical tests, which normally destroy the bones and disable further histological or molecular analyses. Thus, obtaining biomechanical data from bone usually requires an additional group of animals within the experimental setup. Finite element models (FEMs) may non-invasively and non-destructively simulate mechanical characteristics based on material properties. The present study aimed to establish and validate an FEM to predict the mechanical properties of mice tibiae. The FEM was established based on µCT (micro-Computed Tomography) data of 16 mouse tibiae. For validating the FEM, simulated parameters were compared to biomechanical data obtained from 3-point bending tests of the identical bones. The simulated and the measured parameters correlated well for bending stiffness (R2 = 0.9104, p < 0.0001) and yield displacement (R2 = 0.9003, p < 0.0001). The FEM has the advantage that it preserves the bones’ integrity, which can then be used for other analytical methods. By eliminating the need for an additional group of animals for biomechanical tests, the established FEM can contribute to reducing the number of research animals in studies focusing on bone biomechanics. This is especially true when in vivo µCT data can be utilized where multiple bone scans can be performed with the same animal at different time points. Thus, by partially replacing biomechanical experiments, FEM simulations may reduce the overall number of animals required for an experimental setup investigating bone biomechanics, which supports the 3R (replace, reduce, and refine) principle.

Finite element analysis of the Union Plate in treating elderly acetabular fracture patients

BACKGROUND

Elderly acetabular fractures are one of the more difficult types of fractures to resolve. For patients at this age, the more common type of fracture is comminuted. How to better fix this type of fracture has always been an issue of concern. This study was performed to observe the mechanical properties of different internal fixation methods used in treating elderly acetabular fracture patients.

METHODS

A model of a comminuted acetabular fracture in osteoporosis was established, consisting of an anterior column-posterior hemitransverse fracture with disruption of the quadrilateral surface. Fixation of the acetabular fracture model using a reconstruction plate and Union Plates was simulated. For the different internal fixation methods, static and transient modal analyses were performed under different loads, with an action time of 0.21 s and an analysis time of 0.7 s. The stress of the model was observed in the static analysis, and the displacement of the nodes and the entire model in the U1 direction was observed in the transient modal analysis.

RESULTS

In the static analysis, the stress of the osteoporosis model, the suprapectineal pelvic reconstruction plate model, the infrapectineal quadrilateral surface buttress plate model, and the suprapectineal quadrilateral surface buttress plate model were 42.62 MPa, 37.49 MPa, 44.39 MPa, and 46.15 MPa, respectively. The stress was mainly distributed near the suprapubic branch. The corresponding displacement in the U1 direction was 0.1500 mm, 0.1020 mm, 0.0836 mm, and 0.0990 mm, respectively. In the transient modal analysis, there was a significant difference in displacement between the different models (P < 0.05). When different loads were applied with the same fixation method, there was no significant difference in model displacement (P > 0.05).

CONCLUSION

Static and transient modal analyses show that the infrapectineal quadrilateral surface buttress plate or the suprapectineal quadrilateral surface buttress plate has an advantage in maintaining the stability of fracture fragments when fixing comminuted acetabular fractures in elderly individuals. The infrapectineal quadrilateral surface buttress plate also presents better biomechanical results.

Evaluation of experimental, analytical, and computational methods to determine long-bone bending stiffness

Methods used to evaluate bone mechanical properties vary widely depending on the motivation and environment of individual researchers, clinicians, and industries. Further, the innate complexity of bone makes validation of each method difficult. Thus, the purpose of the present research was to quantify methodological error of the most common methods used to predict long-bone bending stiffness, more specifically, flexural rigidity (EI). Functional testing of a bi-material porcine bone surrogate, developed in a previous study, was conducted under four-point bending test conditions. The bone surrogate was imaged using computed tomography (CT) with an isotropic voxel resolution of 0.625 mm. Digital image correlation (DIC) of the bone surrogate was used to quantify the methodological error between experimental, analytical, and computational methods used to calculate EI. These methods include the application of Euler Bernoulli beam theory to mechanical testing and DIC data; the product of the bone surrogate composite bending modulus and second area moment of inertia; and finite element analysis (FEA) using computer-aided design (CAD) and CT-based geometric models. The methodological errors of each method were then compared. The results of this study determined that CAD-based FEA was the most accurate determinant of bone EI, with less than five percent difference in EI to that of the DIC and consistent reproducibility of the measured displacements for each load increment. CT-based FEA was most accurate for axial strains. Analytical calculations overestimated EI and mechanical testing was the least accurate, grossly underestimating flexural rigidity of long-bones.