CLUSTERING A LARGE NUMBER OF FACES FOR 2-DIMENSIONAL MESH GENERATION

Transcription

1 CLUSTERING A LARGE NUMBER OF FACES FOR -DIMENSIONAL MESH GENERATION Keisuke Inoue, Takayuki Itoh, Atsushi Yamada 3, Tomotake Furuhata 4, Kenji Shimada 5,,3,4 IBM Japan Ltd., Yamato-shi, Kanagawa, Japan {inoue,itot,ayamada}@trl.ibm.co.jp, furuhata@jp.ibm.com 5 Carnegie Mellon University, Pittsburgh, PA., U.S.A. shimada@cmu.edu ABSTRACT This paper presents an automated face clustering method used as a preprocess of mesh generation for finite element analysis. Hundreds or thousands of faces are often contained in a CAD model designed in detail, and many of them are smaller than a single mesh element. This complexity of a CAD model makes a meshing process difficult both in quality and in speed. Our clustering method decomposes a CAD model into several regions, each of which is geometrically proper for a meshing process. In the algorithm, we start from the state where every single face makes its own region, and then repeat selecting a pair of adjacent regions and merging them into one region until there remains no mergeable pair of regions. The selection of the most suitable pair and the mergeability test are done based on several geometric indices about a pair of regions. The validity of the method is demonstrated with results of clustering and mesh generation on a real-scale CAD model. Keywords: mesh generation, clustering, simplification, finite element analysis. INTRODUCTION Increasingly, designers prefer to conduct finite element analysis (FEA) on original CAD models rather than to create specific models for analysis alone. However, a CAD model in the final design stage can consist of hundreds or even thousands of faces, each of which is represented as a trimmed parametric surface and is usually smaller than a preferred mesh element size. Although the mesh generation problem is considered to be a discretization process that divides a large surface into small elements, for actual CAD models, the situation is contrary. That is, one mesh element often covers a number of faces. Therefore, it is meaningless to generate a mesh for every face and combine the meshes into one. We have proposed a physically-based approach to generating a mesh for -dimensional objects [][][3][4]. However, when we apply the method to a whole model consisting of a large number of faces including tiny ones, the following problems arise: The whole parametric space is too complex since it consists of parametric spaces of a large number of faces. Iterative calculation for deciding node placement converges very slowly. The quality of the generated mesh is low both topologically and geometrically. These problems become more serious as a CAD model becomes more complex. Our solution to this problem takes the following three steps (See Figure ):. Decompose the entire region into a set of proper regions by face clustering.. Generate a mesh for each region using simplified geometry, i.e., a projected plane. 3. Merge all generated meshes into one mesh. By introducing simplified geometry instead of using the actual complicated geometry of a region, the meshing process

2 can be greatly accelerated and stabilized without worsening the quality of mesh. Practically, ignoring small details contributes much to the quality of the mesh as well since singularities such as CAD model errors or very small curvatures usually reside in those details. In that sense, the first step realizing proper regions is the key in our threestep solution in order to generate a well-shaped mesh in a reasonable time. However, it is very labor-intensitive to perform the first step manually. In this paper, we propose an automated method for the first step in which all faces in a model are clustered into several geometrically proper regions for mesh generation. In our approach, a pair of two adjacent faces are merged into a larger region iteratively. Through the process, each region grows to be a well-shaped region. The criterion for merger is determined based on a combination of geometric indices. This method can be categorized as a simplification method of geometric models for meshing [5][6], but it is different from other methods since it divides the entire region into wellsized sub-regions. This algorithm is similar to the polygon simplification methods [7][8][9][] because they also merge triangular polygons by prioritizing according to a combination of geometric indices. However, we apply very different geometric indices from theirs because our target is not graphics applications but mesh generation in which there are various mesher- or application-dependent requirements for output regions. Moreover, input faces in our problem are characterized by a large variety in: Size: some are tiny and others are so large as to cover a third of the CAD model. Complexity: some faces have very complicated boundaries and/or larger variations of normal vectors than others. These features make our problem very different from those dealt with in other papers. The remainder of the paper is organized as follows. In Section, after describing the problem, we outline our technical approach. Then in Section 3, we describe how to determine whether two regions should be merged or not and how to prioritize possible mergers. In Section 4, we show results of clustering and meshing on a real-scale CAD model. face input CAD model step : meshes of regions region step : clustering step 3: merged mesh Figure. Three-step solution. TECHNICAL APPROACH. Problem Description Our algorithm accepts the following as inputs: A set of faces with adjacency links A mark for each pair of adjacent faces Input faces are assumed to be connected without overlapping each other and to make a -manifold shell as a whole. As for the marks, the user can set either of two marks, INTERNAL or BOUNDARY, on any face pair. The mark INTERNAL represents that the edges shared by the pair of faces must not be on the boundary of a region, while the mark BOUNDARY represents that they must be on a boundary. Every unmarked face pair is listed as a candidate for iterative merger. Our algorithm clusters all input faces into several groups, each of which makes a connected region suitable for mesh generation. Due to a combinatorial explosion, this clustering problem can have a large number of solutions. Two trivial solutions are ) one region for the whole shape, and ) every face making its own region. However, by considering requirements for a well-shaped region, the proper clustering is found somewhere between the two extremes in most cases. A well-shaped region satisfies the following conditions: Area size: The area of a region should be large enough in comparison with that of a mesh element. A small region causes a poor-quality mesh. Boundary smoothness: A boundary of a region should be as smooth as possible. A sharp angle at a vertex on the boundary causes poor quality of mesh elements in that area. Region flatness: Each region should be as flat as possible. A large variation of normal vectors over a region causes failure of the meshing process on the projected plane. The definition of a well-shaped region can vary depending on the method of mesh generation and/or later applications that work on the mesh. For example, region flatness is specifically important only when a projected plane is adopted as a simplified geometry for meshing process, while the former two are required more generally. In any case, our clustering method is customizable to various requirements stated in terms of geometry.. Arc Contraction on an Adjacency Graph By mapping a face to a graph node, and adjacency of a pair of faces to a graph arc respectively, we define an adjacency graph from the set of input faces (See Figure ). Since input faces make a -manifold shell, the adjacency graph becomes planar. A face-merging operation corresponds to an arccontraction operation in this graph.

3 3 4 5 (a) a set of faces The arc-contraction operation consists of two steps (See Figure 3). In the first step, a target arc is deleted and its two terminal nodes are merged into one. At the same time, other arcs incident to the two nodes are connected to the merged node. In the second step, if multiple arcs are produced in the first step, they are unified into a new arc. In the unifying operation, if at least one of the multiple arcs is marked with BOUNDARY, the unified arc is marked with BOUNDARY. If at least one of the multiple arcs is marked with INTERNAL, the unified arc is marked with INTERNAL. There could be a case where one arc of the multiple arcs is BOUNDARY and the other INTERNAL. This means that the input marks are contradicting themselves and they must be corrected. Step (b) an adjacency graph Figure. Conversion from faces to a graph 3 4 Step 5 8. Re-calculate scores for all arcs connected to the merged node. 9. Jump to Step 4. All arcs marked with INTERNAL are resolved in the preprocess; therefore, they never appear in the list. On the other hand, arcs marked with BOUNDARY can appear in the list due to the unification of multiple arcs, even though they are not in the list initially. 3. CONDITIONS FOR MERGING REGIONS 3. Scoring and Mergeability The score is a scalar value calculated for each arc and it is used to prioritize arcs for contraction. The arc with the top score, which is considered to be the most preferable candidate for contraction, is always selected as a target arc. Without proper termination conditions, a sequential contraction based on scores always ends up with a situation where all faces are merged into one region. Therefore, in order to terminate the algorithm in a preferable state, we introduce a test of mergeability, which is a Boolean value determined for each arc. If the target arc selected based on scoring does not pass a mergeability test, the algorithm terminates assuming that there remains no proper arc for contraction. Figure 3. Arc contraction.3 Outline of the Algorithm Our algorithm is basically a sequential arc contraction on an adjacency graph. The process starts from the initial state where every graph node corresponds to one face. Then, a pair of adjacent regions is merged in an iterative manner. If n mergers take place, the total number of regions is reduced by n. The algorithm consists of the following steps:. As a preprocess, contract every arc marked with INTERNAL.. Make a list of all unmarked arcs. 3. Calculate a score for each arc in the list. 4. Take the top-scored arc from the list. 5. If the arc is marked with BOUNDARY, remove the arc from the list and jump to Step If the arc passes the mergeability test (this topic will be discussed in Section 3), contract the arc and remove the arc from the list. At the same time, remove arcs produced during the arc contraction from the list and add the unified arcs. 7. If the list is empty, terminate the algorithm. 3. Geometric Indices for Scoring and Mergeability In Section., we listed three items: () Area Size, () Boundary Smoothness, and (3) Region Flatness, to measure the preference of a region. In this section, we describe indices concretely to calculate a score and mergeability. 3.. Area Size We employ an area index for scoring a pair of regions: Area index ( A ). The area index is defined as follows: A = ( S S or S S ), ( otherwise). where S and S are the areas of the two regions, and S is the minimum allowable area. Infinity, for which we actually use a very large number, implies that the pair of regions has a high priority for merger. 3.. Boundary Smoothness As for boundary smoothness, we employ two kinds of indices:

4 Ratio of shared boundary ( R sh ) Average improvement of contact angle ( ) The ratio of shared boundary, R sh, contributes to the long boundary shared by the face pair and is defined as: l l R sh = max,, L L where l stands for the length of the shared boundary, while L and L stand for the perimeters of the two regions respectively. As shown in Figure 4, a larger value of Rsh implies that two regions share a relatively longer boundary. Figure 4. Ratio of shared boundary The average improvement of contact angle,, reflects the contribution to make the boundary simple. The value is calculated by the following equation: N = N ( φi α i ) + ( φ i β i ), where N is the number of endpoints of discrete segments of the shared boundary (usually equals ), while α, β, and φ are angles at an endpoint as shown in Figure 5. The ideal value of φ is π, which means a straight boundary curve around the contact point. A positive value of implies that the merger of the pair contributes to make the boundary simpler. Figure 5. Angles at a contact point 3..3 Region Flatness For region flatness, we employ three kinds of indices: Gauss map angle (θ ) Variance of normal vector ( σ ) Partial derivative ofσ ( λ ) The Gauss map angle θ is defined as the half angle of the cone circumscribing the Gauss map of normal vectors over a region as shown in Figure 6. A larger angle will contribute to a lower score. Figure 6. Gauss map angle and average normal Variance of normal vector ( σ ) is defined as: σ = { ang( m, n) } S ds, where ang( m, n) represents the angle of two vectors, vector n represents a normal vector at a point on the surface, vector m represents an average normal vector over the surface as shown in Figure 6. Since σ for a merged region can be approximated by: Sσ + S σ σ = S + S SS + ( S + S ) { ang(, )}, m m where S i, σ i, and m i ( i =, ) are respectively area, variance of normal vectors, and an average normal vector for each region. Using the above equation, we define the partial derivative of σ, λ as follows: σ λ = max S σ, S S = S = A large value of λ indicates that the variance increases rapidly when the two faces are merged. 3.3 Combining Indices The geometric indices described in the previous section reflect a certain characteristic of a region or its improvement caused by merger. We use them to define a score function that realizes the desired strategy of clustering. The most primitive way is to use just an index as a scoring function. To allow a more tailored algorithm, a weighted sum of indices is introduced, while the most general form of the score function can be written as: F r, r ) = F( g ( r, r ), g ( r, ),...), ( r where r and r are two regions and index. g i is a geometric The mergeability criterion can be, in the simplest form, an inequality about a geometric index or score function. But in the general form, it is a combination of inequalities about geometric indices. Since different indices reflect different aspects for a pair of regions, improving one index usually leads to worsening another. Therefore, the score function and mergeability

5 criterion should be carefully tuned, and the setting is mesheror application-dependent. 3.4 Multiple-phase implementation In addition to tailoring the scoring function and mergeability criterion, we can apply the clustering algorithm in multiple phases. In other words, after reducing the number of regions with a certain setting for the scoring and mergeability conditions, we can start the second phase from that state with a different setting, and iterate this. In order to allow further mergers, settings in later phases should be looser in some sense. From experiments, we found it effective to apply a couple of settings that focus on different geometric aspects. This hybrid implementation introduces a far larger variety in customizing the clustering algorithm, which should be heuristically explored. Some results will be shown in the next section. 4. RESULTS AND DISCUSSIONS The proposed method was implemented in C and executed on an IBM PowerStation RS/6 Model 43P (PowerPC 64e, 375MHz). Some results are presented in this section. In the following experiments, we used a CAD model of a mechanical part shown in Figure 7, which consists of 74 faces. In Figure 8, each face is colored in different color. There are many tiny or thin faces that are smaller than desired mesh element sizes. We implemented three simple settings for the score function and mergeability criterion, as shown in Table. With these settings, we demonstrate how the clustering method works. Table. Settings of scoring and mergeability Type A: Gauss map angle Type B: variance of normal vector Type C: ratio of shared edges scoring function mergeability criterion F = θ + A θ θ F = λ + A σ σ F = Rsh + A R sh R Figure shows the results of clustering faces by using the scoring pattern Type A shown in Table, and by setting the threshold θ to 5., 45., and 83. (deg). The total numbers of clusters (i.e., regions) are 34, 8, and 36, respectively. They show that the larger θ causes the smaller number of clusters. Figure shows the results by using the scoring pattern Type B shown in Table, and by setting the threshold σ to.5, 7., and 8. (deg). The total numbers of clusters are 7, 8, and 9, respectively. They show that the larger causes the smaller number of clusters. Figure shows the results by using the scoring pattern Type C shown in Table, and by setting the threshold R to.5,.4, and.3. The total numbers of clusters in these results are 33, 87, and 3, respectively. They show that the smaller R causes the smaller clusters. We compared the three scoring patterns in Table, by calculating the following two indices for regions: σ R l : Ratio of the length of the boundary to the circumference of an equivalent circle σ : Standard deviation of normal vector (i.e., square root of variance of normal vector) Table shows the values of indices for three results, which include similar numbers of clusters. It shows that clusters generated by Type C have the smaller R l values. Therefore, it seems that these clusters have simpler boundary shapes than those of other clustering results. It also shows that clusters generated by Types A and B have the smaller σ values than those generated by Type C. Table. Comparison of three types of clusterings pattern Type A Type B Type C mergeability threshold θ = 45. [deg] σ = 7. [deg] R =.4 N (#cluster) R l (average) R l (minimum).9.9. R l (maximum) σ (average) [deg] σ (minimum) [deg]... σ (maximum) [deg] We also generated meshes over the regions generated by the clustering processes. Our implementation first projects the boundaries of regions onto planes, and then generates quaddominant meshes on the planes by using the square packing method [4]. The meshes are finally projected back onto the surfaces of the regions, and merged into one single mesh (See Figure 9). Finally, we evaluated the quality of the meshes generated after the clustering processes. To measure the quality of the meshes, we calculated the ratio of undesirable elements to the total number of elements in a mesh. In this evaluation we measured the following six ratios: N : ratio of triangular elements N : ratio of elements that have small (under 3 degrees) angles

6 N 3 : ratio of elements that have wide (over 35 degrees) angles N 4 : ratio of elements that have short (under.4 times the given element pitch) edges N 5 : ratio of elements that have long (over.4 times the given element pitch) edges N 6 : ratio of folded (at least one angle between two normal vectors is over 3 degrees) elements Table 3 shows the results of evaluating the quality of the meshes. The results are also shown in Figure 4. In Table 3, Type denotes that a mesh is generated over one region consisting of all faces, while Type F denotes that the mesh is generated over regions that are clustered with a multiplephase implementation. For the experimental multiple-phase implementation (Type F), we employed three phases with different settings of the scoring function and mergeability criterion. Strategies in the three phases are outlined as: Phase. Merge regions using strict mergeability criterion on boundary smoothness and region flatness. Phase. Merge regions using mergeability criterion on area only (That is, make every region greater than a certain area). Phase 3. Merge regions using medium mergeability criterion focusing on boundary smoothness. The number of regions was initially 74, and was reduced to 898, 69, and 3 after each phase. Table 3. Comparison of mesh qualities Type Type A Type B Type C Type F N N N N N N regions generated by Type C have better boundaries, but the variance of normal vectors is not always as good. The mesh generated by Type F (i.e., multiple-phase implementation) seems to have no defective values in any aspect. Depending on the requirements of the application, multiple-phase implementation would be the best alternative, but it still needs more extensive work. The whole clustering result is shown in Figure CONCLUSION We have presented a method for clustering a large number of faces for the purpose of -dimensional mesh generation. When dealing with a large CAD model actually used in industry, decomposition of the whole model into well-sized, well-shaped regions contributes significantly to the improvement of the quality of the generated mesh. Our central idea is a sequential merger of adjacent regions. The most promising pair for making a well-shaped region is selected based on a score that is a combination of geometric indices. Experimental results show that well-shaped regions for meshing are obtained by using the multi-phase approach. By configuring the combination of geometric indices and tuning parameters, the method is customizable to various applications that have different requirements on the wellshapedness of a region. The potential limitation of our approach is that it always selects the best available pair at a moment in the clustering process. This approach does not always lead to an optimal solution. To obtain a better result, tactical merging might be necessary, such as accepting an ill-shaped region once and then merging it into a well-shaped region. As future work, instead of tactical merging, we are thinking of developing postprocess functions, such as exchanging faces between two regions locally to achieve better regions. Another topic would be a systematic way of tailoring the combination of geometric indices and of determining threshold values for mergeability tests. Since it is difficult to define a good projection plane for the whole shape, the mesh generated by Type contains many ill-shaped elements as we expected. This resulting mesh is not suitable for FEA. Meshes based on other clustering types have many fewer undesirable elements. It therefore proves the validity of our method. The meshes generated by Type A and Type B have particularly better N 5 and N 6 values, since the regions are generated considering the variance of normal vectors. However, their other characteristics are not as good as those of other types. This is because the regions generated by Type A or B do not always have well-shaped boundaries. Some features of the mesh generated by Type C are opposite to those of the meshes generated by Type A or B, since the

7 REFERENCES [] Kenji Shimada, Physically-Based Mesh Generation: Automated Triangulation of Surfaces and Volumes, PhD thesis, Massachusetts Institute of Technology (993). [] Kenji Shimada, and Takayuki Itoh, Automated Conversion of D Triangular Mesh into Quadrilateral Mesh, Proc. Intl. Conf. on Computational Engineering Science 95, pp (995). [3] Kenji Shimada, Atsushi Yamada, and Takayuki Itoh, Anisotropic triangularation of parametric surfaces via close packing of ellipsoids, Proc. 6th International Meshing Roundtable, pp (997). [4] Kenji Shimada, Jia-Huei Liao, and Takayuki Itoh, Quadrilateral Meshing with Directionality Control through the Packing of Square Cells, Proc. 7th International Meshing Roundtable, pp.6-75 (998). [6] Anton V. Mobley, Michael P. Carroll, and Scott A. Canann, An Object Oriented Approach to Geometry Deseaturing for Finite Element Meshing, Proc. 7th International Meshing Roundtable, pp (998). [7] William J. Schroeder, Jonathan A. Zarge, and William E. Lorensen, Decimation of Triangle Meshes, SIGGRAPH 9 Proc., pp.65-7 (99). [8] Hugues Hoppe, Tony DeRose, Tom Duchamp, John McDonald, and Werner Stuetzle, Mesh Optimization, SIGGRAPH 93 Proc., pp.9-6 (993). [9] Michael Garland, and Paul S. Heckbert, Surface Simplification Using Quadric Error Metrics, SIGGRAPH 97 Proc., pp. 9- (997). [] Andrew J. Willmott, Paul S. Heckbert, and Michael Garland, Face Cluster Radiosity, Eurographics Workshop on Rendering, (999). [5] Alla Sheffer, Ted Blacker, Jan Clements, and Michel Bercovier, Virtual Topology Operations for Meshing, Proc. 6th International Meshing Roundtable, pp (997). Figure 7. Input CAD model with 74 faces

8 Figure 8. Colored faces Figure 9. Merged mesh

9 (a) θ =5.(deg) (34 clusters) (b) θ =45.(deg) (8 clusters) (c) θ =83.(deg) (36 clusters) Figure. Type A clusterings with different thresholds

10 (a) σ =.5 (deg) (7 clusters) (b) σ =7. (deg) (8 clusters) (c) σ = 8. (deg) (9 clusters) Figure. Type B clusterings with different thresholds

11 (a) R =.5 (33 clusters) (b) R =.4 (87 clusters) (c) R =.3 (3 clusters) Figure. Type C clusterings with different thresholds

12 Figure 3. Type F (multiple-phase) clustering (3 clusters) (a) N (b) N,N 3, and N 4 (c) N 5, and N 6 Figure 4. Quality of meshes