We evaluate the cross section for the pp → pp f 1 (1285) and p p → p p f 1 (1285) reactions at near threshold energies relevant for the HADES and PANDA experiments at GSI-FAIR. We assume that at energies close to the threshold the ωω → f 1 (1285) and ρ 0 ρ 0 → f 1 (1285) fusion processes are the dominant production mechanisms. The vertex for the VV → f 1 coupling is derived from an effective coupling Lagrangian. The g ρρ f 1 coupling constant is extracted from the decay rate of f 1 (1285) → ρ 0 γ using the vector-meson-dominance ansatz. We assume g ωω f 1 = g ρρ f 1 , equality of these two coupling constants, based on arguments from the naive quark model and vectormeson dominance. The amplitude for the VV → f 1 fusion, supplemented by phenomenological vertex form factors for the process, is given. The differential cross sections at energies close to the threshold are calculated. In order to determine the parameters of the model the γp → f 1 (1285)p reaction is discussed in addition and results are compared with the CLAS data. The possibility of a measurement by HADES@GSI is presented and discussed. We performed a Monte Carlo feasibility simulations of the pp → pp f 1 reaction for √ s = 3.46 GeV in the π + π − π + π − (not shown explicitly) and π + π − π + π − π 0 final states using the PLUTO generator. The latter one is especially promising as a peak in the π + π − η should be observable by HADES.