This work presents the development of a unified gradient electromechanical theory for thin flexoelectric beams considering both direct and converse flexoelectric effects. The two-way coupled electromechanical theory is developed starting from 3D variational formulation by considering an electric field-strain based free energy function. The formulation incorporates mechanical as well as electrical size effects. The coupled 3D theory is specialized to isotropic materials and a 1D beam theory for composite flexoelectric curved beams is derived using the classical Kirchhoff assumptions. The beam theory is solved using a novel C2 continuous finite element framework for different loading and boundary conditions. Our finite element results are verified with analytical solutions for a simply-supported flexoelectric beam operating in both actuator and sensor modes. The results are also compared with existing literature for the special case of a passive micro-beam. Our computational framework is subsequently used to perform various parametric studies to analyze the effect of electrical and mechanical length scale parameters, geometric parameters like the radius of curvature, flexoelectric layer thickness etc., on the response of the beam. Also, contribution of converse flexoelectricity in the overall response of the flexoelectric beam is compared with that of the direct effect. Our simulation results predict that the converse effect is significant (≈ 10-25% of the direct effect) for a wide range of thickness and length scale parameter values. It is also observed that the effective electromechanical coupling coefficient, calculated in terms of the voltage developed across the flexoelectric layer thickness, is higher in flexoelectric materials compared to piezoelectric materials at smaller length scales (thickness of the order of a few microns). The trends observed in our simulation results also agree with recent experimental findings.