Silk fibers have good biocompatibility and mechanical properties, which make them attractive in biomaterial applications as well as textile industries. It is believed that the superior mechanical property is associated with the crystalline β-sheet structure in the fiber; but a deeper understanding of the structure-property relationship is still needed for full exploitation of its physical properties. Especially, accurate information on hydrogen-bonding interactions within β-sheet domains at the nanoscale and their spatial distributions at the mesoscale are critically needed. In this study, we demonstrate the selective detection of crystalline β-sheet domains in Bombyx mori silk fiber using sum frequency generation (SFG) spectroscopy and its use to determine the angular distribution of the β-sheet crystallites with respect to the fiber axis. Numerical simulations of the SFG signal of the amide-I band were carried out using tensors based on the B2 symmetry of the D2 point group and compared with experimental data. This comparison found that the crystalline β-sheet domains are aligned along the fiber axis with a standard deviation of ∼27° and parallel to the fiber surface with a standard deviation of ∼5°. It was also found that the amide bands in the SFG spectra cannot be fully explained with the assumption that the crystalline β-sheet vibrations can be described with the D2 point group. Being able to monitor the amide group vibrations sensitive to both interchain hydrogen bonding and crystallite orientations, SFG analysis has a potential to unveil the structure-mechanical property relationship that may not be readily assessable with other characterization techniques.