Hyperon puzzle has been a big problem especially after the discovery of 2-solar-mass neutron star. Conflicts exist regarding the appearance of strange hadrons in the core parts of neutron stars. One promising way to solve this problem is to do the study as microscopically as possible using good enough hyperon interactions as input. This is not easy because there are not so much available hypernuclei for the performing of scattering experiments (only a few NY part, and no YY part), which means the present hyperon interactions are far from what we wish for. Alternatively one can refer to accurate measurements of hypernuclei properties, which can be compared with theoretical results based on many body models for the deriving of various features of hyperon interactions. For this, it is necessary to connect the realistic hyperon interactions with hypernuclei data. The present study is along this line. The hyperon interactions used are from the Nijmegen group (NSC97f and ESC08b), the model for infinite nuclear matter is the Brueckner-Hartree-Fock model, and the model employed for hypernuclei is the Skyrme-Hartree-Fock model.