Despite the importance of multiaxial failure of trabecular bone in many biomechanical applications, to date no complete multiaxial failure criterion for human trabecular bone has been developed. By using experimentally validated nonlinear high-resolution, micro-mechanical finite-element models as a surrogate for multiaxial loading experiments, we determined the three-dimensional normal strain yield surface and all combinations of the two-dimensional normal-shear strain yield envelope. High-resolution finite-element models of three human femoral neck trabecular bone specimens obtained through micro-computed tomography were used. In total, 889 multiaxial-loading cases were analyzed, requiring over 41,000 CPU hours on parallel supercomputers. Our results indicated that the multiaxial yield behavior of trabecular bone in strain space was homogeneous across the specimens and nearly isotropic. Analysis of stress-strain curves along each axis in the 3-D normal strain space indicated uncoupled yield behavior, whereas substantial coupling was seen for normal-shear loading. A modified super-ellipsoid surface with only four parameters fit the normal strain yield data very well with an arithmetic error±SD less than –0.04±5.1%. Furthermore, the principal strains associated with normal-shear loading showed excellent agreement with the yield surface obtained for normal strain loading (arithmetic error±SD<2.5±6.5%). We conclude that the four-parameter “Modified Super-Ellipsoid” yield surface presented here describes the multiaxial failure behavior of human femoral neck trabecular bone very well