Cortical bone is a damageable elastic orthotropic material due to the orientation of
its components at different scales. Thus, to consistently describe its mechanical
behavior, it is necessary first to experimentally identify the anisotropy directions and
second to be able to implement them in a numerical model. In our previous paper, we
presented an original approach which enables one to get the numerical directions of
orthotropy throughout a whole femur. In the present paper, we propose a significant
improvement of such a computational approach in order to depict specific
orientations in the femoral diaphysis. To do so, a correction of the main
orientation of the Haversian canals has been applied based on experimental
measurements. The influence of such a correction on the global distribution of the
stresses has been quantified via two finite element applications which simulate
pure torsion and flexion-compression. Our main objective is to provide the
numerical tools to describe the orthotropic behavior of the cortical bone for any
femur.