Short-baseline neutrino anomalies suggest the existence of low-mass ( m sim O(1)~eV) sterile neutrinos u_s. These would be efficiently produced in the early universe by oscillations with active neutrino species, leading to a thermal population of the sterile states seemingly incompatible with cosmological observations. In order to relieve this tension it has been recently speculated that new secret interactions among sterile neutrinos, mediated by a massive gauge boson X (with M_X << M_W), can inhibit or suppress the sterile neutrino thermalization, due to the production of a large matter potential term. We note however, that they also generate strong collisional terms in the sterile neutrino sector that induce an efficient sterile neutrino production after a resonance in matter is encountered, increasing their contribution to the number of relativistic particle species N_ eff. Moreover, for values of the parameters of the u_s- u_s interaction for which the resonance takes place at temperature Tlesssim few MeV, significant distortions are produced in the electron (anti)neutrino spectra, altering the abundance of light element in Big Bang Nucleosynthesis (BBN). Using the present determination of $^4$He and deuterium primordial abundances we determine the BBN constraints on the model parameters. We find that $^2$H/H density ratio exclude much of the parameter space if one assume a baryon density at the best fit value of Planck experiment, Omega_B h^2= 0.02207, while bounds become weaker for a higher Omega_B h^2=0.02261, the 95 % C.L. upper bound of Planck. Due to the large error on its experimental determination, the helium mass fraction Y_p gives no significant bounds.