This thesis deals with developing a thermodynamically consistent model to simulate the electromechanical response of ionic polymer-metal composites based on Euler-Bernoulli beam theory. Constitutive assumptions are made for the Helmholtz free energy and the rate of dissipation. The governing equations involving small deformations are formulated using the conservation laws, the power theorem, and the maximum rate of dissipation hypothesis. The model is extended to solve large deformation cantilever beams involving pure bending which could be used in the characterization of the material parameters. A linear finite element solution along with a staggered time stepping algorithm is provided to numerically solve the governing equations of the small deformations problem under generalized electromechanical loading and boundary conditions. The results are in qualitative and quantitative agreement with the experiments performed on both Nafion and Flemion based Ionic Polymer-Metal Composite strips.