A new three dimensional simulation method is introduced to study the workspace of a 6 PSS (P denotes a prismatic kinematic pair, S denotes a spherical kinematic pair) parallel machine tool. This algorithm adopts the method of numerical analysis to investigate the boundary points in a series of sections which form the surface of the workspace. That is, to study such points that have the largest polar radius on a certain section in a system of polar coordinates according to conditions of constraint. The constraint conditions considered in the article include the maximum and minimum displacements of each dieblock, the maximum and minimum angles of oscillation in each hinge. By converting the constraint inequalities into constraint equations, the largest polar radius corresponding to every constraint condition can be evaluated and the minimum one is used to decide the boundary point. This algorithm greatly simplifies the computational process and can be used to analyze any section of the workspace. It provides a theoretical basis for the structural design of such a machine tool.