A parametric study of hard tissue injury prediction using finite elements : consideration of geometric complexity, failure theory, sub-failure material properties, thresholding, and element characteristics.

This study quantifies the sensitivity of finite element-predicted fracture in the clavicle to several parameters. Clavicles were acquired from 14 donors (age range 14-56 years). Quasistatic axial compression tests were performed using a custom rig designed to replicate in situ boundary conditions. Prior to testing, high-resolution computed tomography (CT) scans were taken of each clavicle. From those images, finite element (FE) models were constructed. Factors varied parametrically included the density used to threshold cortical bone in the CT scans, the presence of trabecular bone, the mesh density. Young's modulus, the maximum stress, and the element type(shell vs. solid. triangular vs. quadrilateral surface elements). The experiments revealed significant variability in the peak force (2.41±0.72 kN)and displacement to peak force (4.9±1.1 mm), with age (p<0.05) and with some geometrical traits of the specimens. In the FE models, the failure force and location were moderately dependent upon the Young's modulus. The fracture force was highly sensitive to the maximum stress (80-110 MPa). Neither fracture location nor force was strongly dependent on mesh density. Both the fracture location and force were strongly dependent upon the threshold density used to define the thickness of the cortical shell. For the covering abstract see ITRD E144229.