This work develops a procedure to find classes of Lagrangian densities that describe generalizations of the Abelian Maxwell-Higgs, the Chern-Simons-Higgs and the Maxwell-Chern-Simons-Higgs models. The investigation focuses on the construction of models that support vortices that obey the stressless condition and lead to first order differential equations which are compatible with the equations of motion. The results induce the appearance of constraints that restrict the choice of the Lagrangian densities, but help us to introduce an auxiliary function that allows to calculate the energy without knowing the explicit form of the solutions. * Electronic address: bazeia@fisica.ufpb.br † Electronic address:In the current work, we concentrate on vortices in relativistic models described in (2, 1) spacetime dimensions. As it is known, vortices were firstly studied by Helmholtz in Ref. [4] and are commonly found in fluid mechanics [5]. They are also present in condensed matter when one studies superconductors. As it is well-known, when superconductors are below a critical temperature they expel the magnetic field, a phenomenom known as the Meissner effect [6]. However, working with the Ginzburg-Landau theory of superconductivity [7], in Ref.[8] Abrikosov noticed that vortices also appear in type II superconductors when exposed to an external electromagnetic field in a specific range of values.The Ginzburg-Landau theory is nonrelativistic. Nevertheless, it is possible to find relativistic field theories that support vortex solutions. The first model was proposed by Nielsen and Olesen in Ref. [9]; it consists of a complex scalar field minimally coupled with a Maxwell gauge field under under the action of the local U (1) gauge symmetry. The equations of the fields that describe the problem, however, are of second order involving the two aforementioned fields that interact in a nontrivial way. Even so, by setting an additional condition to the stress tensor, a classical solution was found in Ref.[10] by de Vega and Schaposnik.Moreover, in Ref.[11] Bogomol'nyi found first order differential equations compatible with the equations of motion by minimizing the energy of the field configurations. The procedure developed by Bogomol'nyi works for kinks, vortices and monopoles, and joined the work of Prasad and Sommerfield [12] on monopoles, to make what is now called BPS states, which represent solutions of first order differential equations that solve the equations of motion and minimize the energy of the non-trivial static field configurations that describe the topological structures.