How graphene flexes and stretches under concomitant bending couples and tractions

We propose a geometrically and physically nonlinear discrete mechanical model of graphene that assigns an energetic cost to changes in bond lengths, bond angles, and dihedral angles. We formulate a variational equilibrium problem for a rectangular graphene sheet with assigned balanced forces and couples uniformly distributed over opposite side pairs. We show that the resulting combination of stretching and bending makes achiral graphene easier to bend and harder (easier) to stretch for small (large) traction loads. Our general developments hold for a wide class of REBO potentials; we illustrate them in detail by numerical calculations performed in the case of a widely used 2nd-generation Brenner potential.