A variational technique has been developed to evaluate the static stress intensity factors of mixed mode problems with mesh free method in this paper. The stiffness is evaluated by regular domain integrals and shape functions are determined by both radial basis function (RBF) interpolation and moving least-square (MLS) method. The stress intensity factors are obtained by two boundary integrals with variation of crack length. The applications of proposed technique to two-dimensional fracture mechanics have been presented with several examples. Comparisons are made with benchmark solutions.