We report numerical simulations of a trapped elastic vortex driven by a strong ac magnetic field H(t) = H sin ωt parallel to the surface of a superconducting film. The surface resistance and the power dissipated by an oscillating vortex perpendicular to the film surface were calculated as functions of H and ω for different spatial distributions, densities and strengths of pinning centers, including bulk pinning, surface pinning and cluster pinning. Our simulations were performed for both the Bardeen-Stephen viscous vortex drag and the Larkin-Ovchinnikov (LO) drag coefficient η(v) decreasing with the vortex velocity v. The local residual surface resistance Ri(H) calculated for different statistical realizations of the pinning potential exhibits strong mesoscopic fluctuations caused by local depinning jumps of a vortex segment as H increases, but the global surface resistance Ri(H) obtained by averaging Ri(H) over different pin configurations increases smoothly with the field amplitude at small H and levels off at higher fields. For strong pinning, the LO decrease of η(v) with v can result in a nonmonotonic field dependence of Ri(H) which decreases with H at higher fields, but cause a runaway instability of the vortex in a thick film for weak pinning. It is shown that overheating of a single moving vortex can produce the LO-like velocity dependence of η(v), but can mask the decrease of the surface resistance with H at a higher density of trapped vortices.