Computationally efficient multiscale finite element model is employed for the nonlinear dynamic response of graphene sheet.

Equations of motion are solved using Newmark’s time integration and shooting technique for steady state frequency response.

Effects of material and geometric nonlinearities, size, boundary conditions, and damping on the response are investigated.

Hardening response with the dominant effect of geometric nonlinearity compared to material nonlinearity.