The concept of acoustic black hole (ABH) designates a decrease in bending stiffness and a localized added damping within a thin structure aiming to localize and efficiently dissipate vibrations. The small thickness results in large amplitude vibrations, which suggests to take into account geometrical nonlinearities. The subject has been previously studied on a beam with a one dimensional ABH. Previous research studied a geometrically nonlinear plate embedding a one dimensional ABH. This study aims at solving the Von Karman equations for a variable thickness plate, via a modal projection using the eigenmodes obtained with a finite element model. The use of a finite element model in this context allows to take into consideration more complex geometries. The methodology is then applied to a rectangular plate embedding a large circular ABH with added damping. The convergence of the nonlinear coefficients is studied and time domain computations show the relative importance of the geometrical nonlinearity.