Abstract: The initial intersection point between surfaces is found by using the convexity property and subdivision method of triangular Bèzier patches，then it is refined by iteration．The whole intersection line is traced out from the initial point and cuts the original surface into two parts．Local area along the trimming border is retriangulated to improve the mesh quality．Tests show that this method is simple，robust and suitable for CAD modeling in reverse engineering．