Voronoi and Delaunay tessellations are applied to pattern recognition of atomic environments and to investigation of the nonlocal order in molecular-dynamics ͑MD͒-simulated materials. The method is applicable also to materials generated using other computer techniques such as Monte Carlo. The pattern recognition is based on an analysis of the shapes of the Voronoi polyhedron ͑VP͒. A procedure for contraction of short edges and small faces of the polyhedron is presented. It involves contraction to vertices of all edges shorter than a certain fraction x of the average edge length, with concomitant contraction of the associated faces. Thus, effects of fluctuations are eliminated, providing ''true'' values of the geometric coordination numbers f , both local and averaged over the material. Nonlocal order analysis involves geometric relations between Delaunay simplexes. The methods proposed are used to analyze the structure of MD-simulated solid lead ͓J.