Methods of patient-specific modeling of the mechanical properties of bone and related systems. The methods include obtaining a bending stiffness (K) of a bone specimen non-invasively and non-destructively in a dynamic 3-point bending test, creating a mathematical mechanical model of the bone specimen, assigning an elastic modulus (E) to the bone specimen of the mathematical mechanical model, determining the flexural rigidity of the bone specimen from simulating the mathematical mechanical model, determining a discrepancy between the flexural rigidity of the bone specimen from the simulated mechanical model and based on the obtained K, adjusting the elastic modulus of the simulated mechanical model to minimize the discrepancy, adjusting the elastic modulus until an optimized elastic modulus is determined where the discrepancy is reduced below a predetermined threshold, and applying the optimized elastic modulus to the simulated mechanical model to determine a strength of the bone.
