在進行模流分析之前，需要將CAD模型轉化為實體網格，以便在分析軟件中進行後續計算。目前，四面體網格發展較為成熟，應用較為普遍。但是，為了減少整體的網格數量，以達到減少計算時間的目的，需要構建六面體網格進行模流分析。而在生成六面體的實體網格之前，必須先建構四邊形的表面網格，然而生成高品質的四邊形網格仍然存在著許多的問題。本研究提出了一種考慮網格頂點相鄰關係的方法，來改善四邊形網格品質。首先，本研究使用B-rep結構的資料，對面上的各個輪廓邊進行分類，並透過分類結果進行均勻輪廓節點或是非均勻輪廓節點的建構。接著，應用本實驗室開發的網格自動建構技術進行網格生成，並在最後透過本研究開發的品質修復演算法，來改善網格的最終品質，其中包含平滑化以及網格合併等程序。最後，透過展弦比以及最小角度來評估生成的四邊形網格，並根據其結果進行分析，其中，所有測試案例中，99.8%的點可以正確輔助網格建構，而經過品質修復後，品質不好的網格占總體網格數量的不到1%，故證明本研究以正確的輔助網格的建構，並改善其網格品質使整體的網格品質優良。

Before performing mold flow analysis, it is necessary to convert the CAD model into a solid mesh for subsequent calculation in the analysis software. At present, the development of tetrahedral mesh is relatively mature, and its application is relatively common. However, in order to reduce the overall number of grids and achieve the purpose of reducing calculation time, it is necessary to construct a hexahedral grid for mold flow analysis. Before generating the hexahedral solid mesh, the quadrilateral surface mesh must be constructed, but there are still many problems in generating high-quality quadrilateral mesh. This study proposes a method that considers the neighbor relationship of mesh vertices to improve the quality of quadrilateral meshes. First, this research uses the B-rep structure data to classify each silhouette edge on the surface, and constructs uniform or non-uniform contour nodes through the classification results. Then, apply the mesh automatic construction technology developed in this laboratory for mesh generation, and finally use the quality restoration algorithm developed in this research to improve the final quality of the mesh, including smoothing and mesh merging procedures. Finally, the generated quadrilateral mesh is evaluated by aspect ratio and minimum angle, and analyzed according to the results. In all test cases, 99.8% of the points can correctly assist mesh construction, and after quality repair, the quality is not good. Good grids account for less than 1% of the total number of grids, so it is proved that this research uses correct auxiliary grid construction and improves its grid quality to make the overall grid quality excellent.

[1] T. D. Blacker and M. B. Stephenson, "Paving: a new approach to automated quadrilateral mesh generation," International Journal for Numerical Methods in Engineering, Vol. 32, No. 4, pp. 811-847, 1990.
[2] X. Roca, J. Sarrate and A. Huerta, "Mesh projection between parametric surfaces," Communications in Numerical Methods in Engineering, Vol. 22, No. 6, pp. 591-603, 2006.
[3] W. J. Gordon and C. A. Hall, "Construction of curvilinear co-ordinate systems and applications to mesh generation," International Journal for Numerical Methods in Engineering, Vol. 7, No. 4, pp. 461-477, 1973.
[4] S. H. Lo, " A new mesh generation scheme for arbitrary planar domains," International Journal for Numerical Methods in Engineering, Vol. 21, No. 8, pp. 1403-1426, 1985.
[5] W. R. Buell and B. A. Bush, "Mesh generation-a survey," Journal of Manufacturing Science and Engineering, Vol. 95, No.1, pp. 332-338 , 1973.
[6] L. F. Diachin, P. Knupp, T. Munson and S. Shontz, "A comparison of two optimization methods for mesh quality improvement," Engineering with Computers, Vol. 22, No.2, pp. 61-74 , 2006.
[7] E. Amezua, M.V. Hormaza, A. Hernández and M. B. G. Ajuria ," A method for the improvement of 3D solid finite-element meshes," Advances in Engineering Software, Vol. 22, No. 1, pp. 45-53 , 1995.
[8] L. A. Freitag, " On combining Laplacian and optimization-based mesh smoothing techniques," Trend in Unstructured Mesh Generation, Vol. 220, pp. 37-43, 1997.
[9] M. S. Joun and M. C. Lee, " Quadrilateral finite-element generation and mesh quality control for metal forming simulation," International Journal for Numerical Methods in Engineering, Vol. 40, No. 21, pp. 4059-4075, 1997.
[10] J. Kim, T. Panitanarak and S. M. Shontz, " A multiobjective mesh optimization framework for mesh quality improvement and mesh untangling," International Journal for Numerical Methods in Engineering, Vol. 94, No. 1, pp. 20¬¬-42, 2013.
[11] H. Yagou, Y. Ohtake and A. Belyaev, "Mesh smoothing via mean and median filtering applied to face normals," Proc. of Geometric Modeling and Processing , pp. 124-131, 2002.
[12] R. E. Jones, "A self-organizing mesh generation," Journal of Pressure Vessel Technology, Vol. 96, No. 3, pp. 193-199,1974.
[13] Rhinoceros. Available: https://www.rhino3d.com. [Accessed 18 June 2023].
[14] OpenNURBS. Available:https://www.rhino3d.com/tw/features/developer/ opennurbs. [Accessed 18 June 2023].
[15] B-rep structure. Available: https://developer.rhino3d.com/guides/cpp/brep-data-structure. [Accessed 18 June 2023].
[16] 宋培溥，「應用於模流分析之薄殼元件CAD模型特徵辨識與分解技術發展」，國立中央大學博士論文，2021.
[17] 吳添瑜，「B-rep模型凸起特徵辨識及分解技術發展」，國立中央大學碩士論文，2022.
[18] 方聲翔，「自動化分解薄殼CAD模型殼本體之增強演算法發展」，國立中央大學碩士論文，2022.
[19] 余晉歡，「CAD模型體積分解技術之優化及改進」，國立中央大學碩士論文，2023.
[20] 鄭旭佑，「自動化多邊形分解與結構化網格生成方法研究」，國立中央大學碩士論文，2023.
[21] J. Y. Lai, P. Putrayudanto, D. H. Chen,J. H. Huang, P. P. Song and Y. C. Tsai , "An enhanced paving algorithm for automatic quadratic generation of IC CAD models,"The 9th IEEE International Conference on Applied System Innovation,[Online], 2023.
[22] 戴宇辰，「混合結構化與非結構化四邊形網格之自動區域劃分方法發展」，國立中央大學碩士論文，2023.