  We extend the quark mean field model to the study of $\Lambda$ hypernuclei.
Without adjusting parameters, the properties of $\Lambda$ hypernuclei can be
described reasonably well. The small spin-orbit splittings for $\Lambda$ in
hypernuclei are achieved, while the $\Lambda$ single particle energies in the
present model are slightly underestimated as compared with the experimental
values. About 3% deviation from the quark model prediction for the
$\omega-\Lambda$ couplings is required in order to reproduce the experimental
single particle energies.
