Turbulence in tokamaks generates radially sheared zonal flows (ZFs). Their oscillatory counterparts, geodesic acoustic modes (GAMs), appear due to the action of the magnetic field curvature. The GAMs can also be driven unstable by an anisotropic energetic particle (EP) population leading to the formation of global radial structures, called EGAMs. The EGAMs might play the role of an intermediate agent between the EPs and thermal plasma, by redistributing EP energy to the bulk plasma through collisionless wave-particle interaction. In such a way, the EGAMs might contribute to the plasma heating. Thus, investigation of EGAM properties, especially in the velocity space, is necessary for precise understanding of the transport phenomena in tokamak plasmas. In this work, the nonlinear dynamics of EGAMs is investigated with the help of a Mode-Particle-Resonance (MPR) diagnostic recently implemented in the global gyrokinetic (GK) particle-in-cell code ORB5. This enables to investigate the relative importance and the evolution of the resonances responsible for the ion and electron Landau damping, and of the EP drive. An ASDEX Upgrade discharge is chosen as a reference case for this investigation due to its rich EP nonlinear dynamics.