The equilibrium position of an Euler-Bernoulli beam is investigated. Using the discrete Euler-Lagrange and Lagrange-d’Alembert principles, the behavior of the beam is simulated using variational integrators and, in particular, AVIs (Asynchronous Variational Integrators). Special emphasis is placed on the discrete mechanics and geometric structure underlying stress resultant beam models.