The classical radial point interpolation method (RPIM) is a powerful meshfree numerical technique for engineering computation. In the original RPIM, the moving support domain for the quadrature point is usually employed for the field function approximation, but the local supports of the nodal shape functions are always not in alignment with the integration cells constructed for numerical integration. This misalignment can result in additional numerical integration error and lead to a loss in computation accuracy. In this work, a modified RPIM (M-RPIM) is proposed to address this issue. In the present M-RPIM, the misalignment between the constructed integration cells and the nodal shape function supports is successfully overcome by using a fixed support domain that can be easily constructed by the geometrical center of the integration cell. Several numerical examples of free vibration analysis are conducted to evaluate the abilities of the present M-RPIM and it is found that the computation accuracy of the original RPIM can be markedly improved by the present M-RPIM.