In this paper, we propose a new method to construct energy-minimizing cubic B-spline interpolation curves by minimizing the approximated stretch energy. The construction of a B-spline interpolation curve with a minimal approximated stretch energy can be addressed by solving a sparse linear system. The proof of both the existence and uniqueness of the solution for the linear system is provided. In addition, we analyze the computational cost of cubic B-spline curves with an approximated stretch energy, which is close to that of the ordinary interpolation method with cubic B-splines without the requirement of stretch energy.