  Nanoindentation of bilayer graphene is studied using molecular dynamics
simulations. We compared our simulation results with those from elasticity
theory as based on the nonlinear F\"{o}ppl-Hencky equations with rigid boundary
condition. The force deflection values of bilayer graphene are compered to
those of monolayer graphene. Young's modulus of bilayer graphene is estimated
to be 0.8 TPa which is close to the value for graphite. Moreover, an almost
flat bilayer membrane at low temperature under central load has a 14$%$ smaller
Young's modulus as compared to the one at room temperature.
