In order to exchange data between PDE(partial differential equation)modeling technique and the conventional CAD(computer aided design)systems,this paper presents a novel algorithm for approximating PDE surface by tensor product Bézier surface based on constrained optimization.The algorithm is also improved by the subdivision property of tensor product Bézier surface.The computational examples and the results of error analysis demonstrate the validity of the algorithm.