The force method is a simple and accurate technique to obtain the stress intensity factors (SIF) for both modes I, II and also mixed I+II modes of fracture. The method uses the summation of internal nodal forces in the vicinity of the crack tip to compute SIFs. Recently, de Morais 1 showed that the force method is able to yield accurate SIF values from FE models constructed with regular meshes of linear elements. In this paper, the force method is applied successfully to general finite element meshes, in such a way that it can be used on crack propagation algorithms with arbitrary crack paths.