A method of patient-specific modeling of mechanical properties of bone comprising: obtaining a bending stiffness of a bone specimen non-invasively and non-destructively in a dynamic 3-point bending test using a device for measuring the bending stiffness of bone in vivo; determining a flexural rigidity of the bone specimen based on the obtained K in accordance with the equation EI=KL 3 /48 where L represents the length of the bone specimen; characterized by: imaging the bone specimen with at least one of computed tomography, dual-energy X-ray absorptiometry, and magnetic resonance imaging to create a geometric model of the bone specimen; and instructing a processor to perform the steps of: (i) creating a mathematical mechanical model of the bone specimen from the geometric model; (ii) assigning an elastic modulus to the bone specimen of the mathematical mechanical model; (iii) determining the flexural rigidity of the bone specimen from simulating the mathematical mechanical model; (iv) determining a discrepancy between the flexural rigidity of the bone specimen from the simulated mechanical model and the flexural rigidity of the bone specimen based on the obtained K; (v) adjusting the elastic modulus of the bone specimen of the simulated mechanical model to minimize the discrepancy and repeating steps (iii) and (iv); (vi) repeating step (v) until an optimized elastic modulus is determined where the discrepancy is reduced below or to a predetermined threshold; and (vii) applying the optimized elastic modulus to the simulated mechanical model to determine a strength of the bone, where the optimized elastic modulus represents a patient-specific E of the bone specimen. obtaining a bending stiffness of a bone specimen non-invasively and non-destructively in a dynamic 3-point bending test using a device for measuring the bending stiffness of bone in vivo; determining a flexural rigidity of the bone specimen based on the obtained K in accordance with the equation EI=KL 3 /48 where L represents the length of the bone specimen; imaging the bone specimen with at least one of computed tomography, dual-energy X-ray absorptiometry, and magnetic resonance imaging to create a geometric model of the bone specimen; and instructing a processor to perform the steps of: (i) creating a mathematical mechanical model of the bone specimen from the geometric model; (ii) assigning an elastic modulus to the bone specimen of the mathematical mechanical model; (iii) determining the flexural rigidity of the bone specimen from simulating the mathematical mechanical model; (iv) determining a discrepancy between the flexural rigidity of the bone specimen from the simulated mechanical model and the flexural rigidity of the bone specimen based on the obtained K; (v) adjusting the elastic modulus of the bone specimen of the simulated mechanical model to minimize the discrepancy and repeating steps (iii) and (iv); (vi) repeating step (v) until an optimized elastic modulus is determined where the discrepancy is reduced below or to a predetermined threshold; and (vii) applying the optimized elastic modulus to the simulated mechanical model to determine a strength of the bone, where the optimized elastic modulus represents a patient-specific E of the bone specimen. (i) creating a mathematical mechanical model of the bone specimen from the geometric model; (ii) assigning an elastic modulus to the bone specimen of the mathematical mechanical model; (iii) determining the flexural rigidity of the bone specimen from simulating the mathematical mechanical model; (iv) determining a discrepancy between the flexural rigidity of the bone specimen from the simulated mechanical model and the flexural rigidity of the bone specimen based on the obtained K; (v) adjusting the elastic modulus of the bone specimen of the simulated mechanical model to minimize the discrepancy and repeating steps (iii) and (iv); (vi) repeating step (v) until an optimized elastic modulus is determined where the discrepancy is reduced below or to a predetermined threshold; and (vii) applying the optimized elastic modulus to the simulated mechanical model to determine a strength of the bone, where the optimized elastic modulus represents a patient-specific E of the bone specimen.