By using the principle of stress release method, a numerical solution method of convergence-confinement characteristic curve was proposed for tunnels with various cross-sections based on the convergence-confinement analysis method. The validity of the proposed method was verified by comparison with the convergence-confinement analytical method, and compared with the usual numerical methods. Finally, the calculation error of the numerical solution method for convergence-confinement characteristic curve was analyzed. The results show that the numerical solution method of convergence-confinement characteristic curve is suitable for non-circular tunnel, and due to the stress release achieved by applying support reaction to nodes, the method is also suitable for the tunnel with complicated stress condition. Taking the Dianzhong diversion project as an example,the numerical solution method of convergence-confinement curve was compared with the calculation results of the convergence-confinement analytical method, and the validity of the numerical solution method of convergence-confinement curve was verified. Despite some error in the result of the numerical solution method of convergence-confinement curve, the overall relative error is less than 10%. The application of the numerical solution method for convergence-confinement characteristic curve is more extensive,which can provide reference for the optimization design of primary support and safety evaluation of tunnel engineering, and the numerical solution method has a certain value of engineering application.