  The small mass and atomic-scale thickness of graphene membranes make them
highly suitable for nanoelectromechanical devices such as e.g. mass sensors,
high frequency resonators or memory elements. Although only atomically thick,
many of the mechanical properties of graphene membranes can be described by
classical continuum mechanics. An important parameter for predicting the
performance and linearity of graphene nanoelectromechanical devices as well as
for describing ripple formation and other properties such as electron
scattering mechanisms, is the bending rigidity, {\kappa}. In spite of the
importance of this parameter it has so far only been estimated indirectly for
monolayer graphene from the phonon spectrum of graphite, estimated from AFM
measurements or predicted from ab initio calculations or bond-order potential
models. Here, we employ a new approach to the experimental determination of
{\kappa} by exploiting the snap-through instability in pre-buckled graphene
membranes. We demonstrate the reproducible fabrication of convex buckled
graphene membranes by controlling the thermal stress during the fabrication
procedure and show the abrupt switching from convex to concave geometry that
occurs when electrostatic pressure is applied via an underlying gate electrode.
The bending rigidity of bilayer graphene membranes under ambient conditions was
determined to be $35.5^{+20}_{-15}$ eV. Monolayers have significantly lower
{\kappa} than bilayers.
