Multi-phonon excitations in atomic nuclei were observed very rarely although collective motions in quantum many-body systems are described as bosonic excitations. In particular, the first two-phonon $gamma$ vibrational ($2gamma$) excitation in odd-$A$ nuclei was reported in 2006 and only a few have been known. Quite recently, conspicuously enhanced $B(E2)$s feeding $2gamma$ states were observed in $^{105}$Nb and conjectured that their parent states are candidates of $3gamma$ states. In the present work, the model space is enlarged from the present authors previous calculation for $^{103}$Nb. The purpose is twofold: One is to see how the description of $2gamma$ states is improved, and the other is to examine the existence of collective $3gamma$ states, and when they exist, study their collectivity through calculating interband $B(E2)$s. The particle-vibration coupling model based on the cranking model and the random-phase approximation is used to calculate the vibrational states in rotating odd-$A$ nuclei. Interband $B(E2)$s are calculated by adopting the method of the generalized intensity relation. The present calculation reproduces the observed spectra of $0gamma$ - $2gamma$ states well and gives collective $3gamma$ states with enhanced $B(E2)$s to $2gamma$ states in $^{103}$Nb and $^{105}$Nb. The most collective $3gamma$ state with the highest $K$ at zero rotation is thought to be the main component of the observed band.