pts = Table[{ Cos[2 Pi u/6] Cos[v], Sin[2 Pi u/6] Cos[v], v}, {u, 6}, {v, -1, 1, 1/2}];
f = BSplineFunction[pts];
graph = ParametricPlot3D[f[u, v], {u, 0, 1}, {v, 0, 1}, Mesh -> None];
data = FirstCase[graph, _GraphicsComplex];
polygons = Cases[data, _Polygon, Infinity];