In this manuscript, classical molecular dynamics simulations (MD) have been applied to study the short and medium range order of very complex vanadium containing glasses with the aim of improving the first microscopic picture of such materials. A rigid ionic force-field has been extended to include the V 5+ -O, V 4+ -O and Cu 2+ -O interatomic pair parameters and tested to reproduce structural properties of known crystal phases with quite good accuracy. Then the structure of Na 2 O-SiO 2 , CaO-MgO-Al 2 O 3 -SiO 2 and Na 2 O-P 2 O 5 glass compositions in which vanadium is present in the range 1-72 wt.% (0.3-60 mol.%) have been fully described in terms of vanadium local structure and Qn distributions. A fairly good agreement was found with experimental data further validating our computational models and providing a computational approach that could be used and extend to investigate in detail the structural information (V-V distances, V-O-V linkages and BO/NBO) directly correlated to macroscopic properties of application interest.