AbstractA practical finite element method is presented for geometrically non‐linear (von Kármán) plate problems. Symmetric strain energy tensors provide an efficient formulation and solution, and allow the effects of initial imperfections to be considered. An unbalanced force iteration is employed, with the advantage that the time consuming Gauss decomposition needs to be performed only once. The convergence is controlled by an automatically computed relaxation factor, and this makes the method applicable in advanced stage of non‐linearity.Several different triangular elements are used in explicitly the required tensors. Buckling and post‐buckling solutions are outlined and demonstrated by numerical examples.