This paper presents an efficient adaptive analysis procedure being able to operate in the framework of the node-based smoothed point interpolation method (NS-PIM). The NS-PIM uses three-node triangular cells and is very easy to be implemented, which make it an ideal candidate for adaptive analysis. In the present adaptive procedure, a new error indicator is devised for NS-PIM settings; two ways are proposed to calculate the local critical value; a simple h-type local refinement scheme is adopted and Delaunay technology is used for regenerating optimal new mesh. A number of typical numerical examples involving stress concentration and solution singularities have been tested. The results demonstrate that the present procedure achieves much higher convergence rate results compared to the uniform refinement, and can obtain upper bound solution in strain energy. (C) 2011 Elsevier Inc. All rights reserved.