The nanostructure of bone consists of two hierarchical levels called 1) Mineralized Collagen Fibril (MCF) and 2) Microfibril. At both levels, the nanostructure contains different volume fractions and arrangements of highly tough collagen and highly stiff mineral phases, which yields unique properties like lightweight, high strength, stiffness, and toughness. A finite element (FE) model has been developed, where the inclusions are considered as one-dimensional embedded elements, and the matrix is modeled with two-dimensional plane stress continuum elements. Two-scale asymptotic homogenization approach has been used to derive the mechanical properties of MCF. In addition to this, the effect of uncertainties associated with geometric and material parameters of mineral and collagen phases has also been considered directly in the FE mesh. Different realizations of MCF are generated using in-house MATLAB code, and Monte-Carlo type simulations are performed to obtain the characteristic elastic modulus and strength of MCF. The results are found to be in good agreement with experiments and molecular dynamics simulations. The present model has closely predicted the local distribution of stresses and strains with a significant reduction in the computational cost. The work has been published in the Journal of Biomechanics (click here) and the Journal of Biomechanical Engineering (click here).