In the presence of relatively high freestream turbulence, laminar-turbulent transition in boundary layers follows a so-called bypass route. This process, in general, is described as follows: low-frequency freestream vortical disturbances (FSVDs) with large amplitude excite streaky disturbances that grow transiently in boundary layers; high-frequency disturbances amplify exponentially due to the secondary-instability regime; streaks break down and the intermittent turbulence spots emerge when the secondary modes reach finite amplitude. However, one process still remains unclear, namely, the mechanism that the secondary instability modes are excited by in the early stage of bypass transition. This paper first verifies the solutions of the boundary-region and edge-layer equations on describing the boundary-layer response to low- and high-frequency FSVDs, respectively. Then, using the above models as the inflow perturbation, the generation of the secondary instability modes is studied by direct numerical simulations. The transformation of Wigner-Ville distribution is employed on analyzing the numerical data, and the formation of the secondary instability modes in the vicinity of the neutral position and time is exhibited in detail. The receptivity process is described as the interaction between the high-frequency FSVDs and the non-parallel streaky disturbance that generates the secondary instability modes when the latter first becomes unstable. A simplified model is established by selecting a steady base flow according to the group velocity of the excited secondary modes, which is confirmed to be an efficient means to determine the receptivity coefficient. Moreover, the nonlinear evolution of the secondary instability modes is simulated, exhibiting the whole process of bypass transition.