The two-dimensional numerical simulation of the flow over a cylinder forced to oscillate in the streamwise direction for Re = 200 is performed in CFX ANSYS. The controlled-vibration comprises of prescribed inline vibration from displacement amplitude-to-cylinder diameter A/D = 0.05 up to 0.5 with the excitation frequency ratios of 1, 1.5 and 2 including the harmonic and superharmonic excitation regions. The immersed boundary method is used to model the cylinder oscillation. Modal decomposition of the transverse velocity field via the proper orthogonal decomposition (POD) method is applied to uncover the interaction of symmetric and antisymmetric modes of the near wake. A model using the first two POD modes is developed based on symmetry group equivariance. The model predicts the mode interactions and bifurcated solution branches for all cases, and is shown to be in good agreement with numerical as well as previous experimental results. Lock-on is determined for a range of values of the oscillation amplitudes and frequency ratios. It is shown that the lock-on range widens with increasing nondimensional oscillation amplitude. The asymmetric 2S, P + S and symmetric pattern S with symbol S for a single vortex and P for a vortex pair shed per cycle, as well as a regime in which vortex formation is not synchronized with cylinder motion are observed in the cylinder wake depending on the combination of oscillation amplitude and frequency ratio. The frequency ratio variation from 1 to 2 leads to the switching from asymmetric to symmetric modes. The symmetric flow pattern corresponds to a near zero lift coefficient on the cylinder.