Abstract: An interpolating 3^1/2-subdivision method based on Bézier triangle is proposed. Assume every vertex of a triangular mesh is equipped with a unit normal vector, we subdivide the mesh in three main steps: construct a cubic Bézier triangle interpolating the vertices and vertex normals of every triangle on the control mesh; compute the center point and center normal vector for every Bézier triangle; re-triangulate the newly computed center points and old vertices. A method to subdivide boundary triangles or singular triangles with thin shapes is also proposed. Since every new face vertex and new vertex normal is only computed based on the vertices and vertex normal of a local triangle, this method is very easy to implement and it is also possible to accelerate the subdivision by parallel computing. Numerical examples show that the proposed method can produce smooth and fair-looking subdivision surfaces. Moreover, the shape of the limit surface can be edited effectively by adjusting the normal direction of the initial vertices.