two major kinematic hardening models are employed in this article to study the behavior of Euler-Bernoulli beams by a finite element formulation of the plasticity problem. The proposed finite element formulation uses a variable stiffness matrix in each incremental step reflecting the yield surface movement. Moreover, complete formulation is presented for both load and deformation controlled cases. Examples are worked out for the Ziegler-Prager and the Armstrong-Frederick theories, to show the stress-strain behavior under cyclic symmetric and asymmetric flexural loading. The results have been graphically illustrated in plots of the stress-strain curves and are compared to the published and experimental ones. It was observed that Ziegler-Prager theory for isotropic cases with symmetric and asymmetric loading conditions. The presented simulation results confirm that the anisotropic cases with symmetric loading exhibit a ratcheting response. While the results show agreement with published ones; it was also observed that the two kinematic hardening theories do not show similar responses of reverse plasticity or ratcheting for Euler-Bernoulli beams in all the example cases.