Patient-specific biomechanical modeling using Finite Element Analysis (FEA) is pivotal for understanding the structural health of bones, optimizing surgical procedures, assessing outcomes, and validating medical devices, aligning with guidance issued by standards and regulatory bodies. Accurate mapping of image-to-mesh-material is crucial given bone's heterogeneous composition. This study aims to rigorously assess mesh convergence and evaluate the sensitivity of material grouping strategies in quantifying bone strength. Subject-specific geometry and nonlinear material properties were derived from computed tomography (CT) scan data of one cadaveric human vertebral body. Linear tetrahedral elements with varying edge lengths between 2mm and 0.9mm were then generated to study the mesh convergence. To compare the effectiveness of different grouping strategies, three approaches were used: Modulus Gaping (a user-defined absolute threshold of Young's modulus ranging from 500 MPa to 1 MPa), Percentual Thresholding (relative parameter thresholds ranging from 50% to 1%), and Adaptive clustering (unsupervised k-means-based clustering ranging from 10 to 200 clusters). Adaptive clustering enables a constant number of unique material properties in cross-specimen studies, improving the validity of results.Introduction
Method
Functional Spine Units (FSUs) play a vital role in understanding biomechanical characteristics of the spine, particularly bone fracture risk assessment. While established models focus on simulating axial compression of individual bones to assess fracture load, recent models underscore the importance of understanding fracture load within FSUs, offering a better representation of physiological conditions. Despite the limited number of FSU fracture studies, they predominantly rely on a linear material model with an annulus fibrosus Young's modulus set at 500 MPa, significantly higher than stiffness values (ca. 4 MPa) utilized in other FSU and spine section biomechanical models. Thus, this study aims to study the effect of varying annulus fibrosus stiffness on FSU fracture load, aiming to identify physiologically relevant biomechanical parameters. Subject-specific geometry and material properties of bones were derived from computed tomography (CT) image data of five human cadaveric FSU specimens. The annulus fibrosus and nucleus pulposus were manually recreated and assigned linear elastic material properties. By subjecting the model to axial compression, the fracture load of the FSU was deduced from the peak of the force-displacement graph. To explore the effect of stiffness of the annulus fibrosus on simulated fracture load, we conducted a parameter study, varying stiffness values from the high 500 MPa to a more physiologically relevant 25 MPa, aiming to approximate values applied in FSU kinematic models while achieving bone fracture.Introduction
Method
