Efficient materially nonlinear [Formula: see text]FE solver for simulations of trabecular bone failure.

An efficient solver for large-scale linear [Formula: see text] simulations was extended for nonlinear material behavior. The material model included damage-based tissue degradation and fracture. The new framework was applied to 20 trabecular biopsies with a mesh resolution of [Formula: see text]. Suitable material parameters were identified based on two biopsies by comparison with axial tension and compression experiments. The good parallel performance and low memory footprint of the solver were preserved. Excellent correlation of the maximum apparent stress was found between simulations and experiments ([Formula: see text]). The development of local damage regions was observable due to the nonlinear nature of the simulations. A novel elasticity limit was proposed based on the local damage information. The elasticity limit was found to be lower than the 0.2% yield point. Systematic differences in the yield behavior of biopsies under apparent compression and tension loading were observed. This indicates that damage distributions could lead to more insight into the failure mechanisms of trabecular bone.