In this work, a numerical approach is presented for solving problems of finitely deformed membrane structures made of hyperelastic material and subjected to external pressure loadings. The problem is geometrically non-linear since it deals with finite displacements and rotations as well as finite strains. Furthermore, the material nonlinear behaviour is described by use of a hyperelastic constitutive relationship. Usually, the classical finite element method based on satisfaction of the equilibrium equations requires that geometric stiffness matrix and large displacement matrix are taken into account. Moreover, structure instabilities or bifurcations may occur under large displacements and material non-linearity, which must be predicted by a specific treatment. To correctly handle all the previous numerical problems, the classical finite element method usually leads to highly CPU time consuming programs.In this paper, the membrane structures are solved by directly minimizing the total potential energy. Although theoretically equivalent to the traditional finite element method, the proposed approach proves to be a more robust and efficient technique. All the numerical computations are performed using a three-dimensional flexible triangular finite element, which is able to handle the abovementioned non-linearities. The finite element has no bending stiffness and the tangent matrix is singular for the rotation degrees of freedom. The numerical solution is carried out by means of an iterative method like the conjugate gradient or Newton one.The strain energy is given by a quasi-incompressible hyperelastic potential density in the reference configuration written in the following form:$$\Psi=\frac{\mu}{2}(\mathrm{tr}(C)-3)-\mu\ln(J)+\frac{\lambda}{2}\ln^2(J)$$First, the proposed finite element is validated through the well-known Hencky problem [1] which consists of the inflation of a circular membrane clamped at the rim. For this problem we recall the analytical solution in the case of isotropic linear elastic behavior [2,3], before the comparison between the analytical solution and the finite element results. Next, an extension to a nonlinear hyperelastic behavior is treated, which allows us to exhibit the pressure value which leads to the material instability of the circular membrane.The second application is related to problem of an orthotropic rectangular membrane. We focus on two points: first we present a semi-analytical solution in the case of orthotropic linear behavior for the rectangular membrane. Then, this solution is used to validate once again the finite element in the case of isotropic behavior.