Submerged flexible aquatic vegetation exists widely in nature and achieves multiple functions mainly through fluid–structure interactions (FSIs). In this paper, the evolution of large-scale vortices above the vegetation canopy and its effect on flow and vegetation dynamics in a two-dimensional (2-D) laminar flow are investigated using numerical simulations under different bending rigidity γ and gap distance d. According to the variation of large-scale vortex size and intensity, the evolution process is divided into four distinct zones in the streamwise direction, namely the ‘developing’ zone, ‘transition’ zone, ‘dissipation’ zone and ‘interaction’ zone, and different evolution sequences are further classified. In the ‘developing’ zone, the size and intensity of the large-scale vortex gradually increase along the array, while they decrease in the ‘dissipation’ zone. The supplement of vegetation oscillating vortices to large-scale vortices is the key to the enhancement of the latter. The most obvious dissipation of large-scale vortices occurs in the ‘transition’ zone, where the position of the large-scale vortex is significantly uplifted. The effects of γ and d on the evolution of the large-scale vortex are discussed. In general, the features of vegetation swaying vary synchronously with those of large-scale vortices. The flow above the canopy is dominated by large-scale vortices, and the development of flow characteristics such as time-averaged velocity profile and Reynolds stress are closely related to the evolution of large-scale vortices. The flow inside the canopy, however, is mainly affected by the vortex shed by the vegetation oscillation, which leads to the emergence of negative time-averaged velocity and negative Reynolds stress.