Optical second harmonic generation (SHG) from nanostructured graphene has been studied in the framework of classical electromagnetism using a surface integral equation method. Single disks and dimers are considered, demonstrating that the nonlinear conversion is enhanced when a localized surface plasmon resonance is excited at either the fundamental or second harmonic frequency. The proposed approach, beyond the electric dipole approximation used in the quantum description, reveals that SHG from graphene nanostructures with centrosymmetric shapes is possible when retardation effects and the excitation of high plasmonic modes at the second harmonic frequency are taken into account. Several SHG effects similar to those arising in metallic nanostructures, such as the silencing of the nonlinear emission and the design of double resonant nanostructures, are also reported. Finally, it is shown that the SHG from graphene disk dimers is very sensitive to a relative vertical displacement of the disks, opening new possibilities for the design of nonlinear plasmonic nanorulers.