Unless the targeted mesh is developable, metric distortion is inevitable during the process of surface mesh parameterization, thus one important objective of all involved parametric studies is to reduce the metric distortion. In order to further reduce area and angle distortion, a novel method of boundary-free mesh parameterization is presented in the paper. Firstly, the initial boundary-fixed conformal parameterization from 3D surface mesh patch to a plane is performed in the method. Then, based on the initial parameterization, the iterations of boundary-free quasi-harmonic parameterization are developed, where the tensor field is updated in each iterative step and the principal curvature direction is utilized to terminate the iteration. The solution of the novel method is convenient to calculate since it involves a series of linear systems. In our novel parameterization method, lower metric distortion and considerable efficiency have been obtained in experiments.