在IC封裝領域，模流分析(Mold flow analysis)可預測成品生產前的潛在缺陷，而實體網格品質將影響分析結果的準確性。其中六面體網格(Hexahedron)的品質最佳，但自動化建構難度高。六面體網格由六片四邊形表面網格(Quadrilateral mesh)構成，隨表面網格逐層堆疊形成六面體網格，故四邊形網格的建構極為關鍵。因此，本實驗室以Paving 演算法為基礎，進行非結構化四邊形網格自動建構技術的開發，過程中，網格自交與共線為最常見且嚴重的問題，將導致網格建構失敗。本研究以14個IC CAD模型合計共進行超過1,400組網格尺寸的測試，並針對該技術的整體流程進行改善與優化，包含網格自交偵測，並著重於網格修正、網格局部建構以及函式合併等三個方面。另外，本研究亦開發多模型多尺寸自動測試模組，以提升測試與偵錯效率，並驗證其技術之穩定性與計算效率。測試結果顯示，所有模型皆可完成網格建構，且大多數模型的建構時間可縮短，以建構時間最久之模型為例，平均建構時間減少6.77%。故本研究可提升非結構化四邊形網格自動建構技術的穩定性與計算效率。

In the field of IC packaging, mold flow analysis is used to predict potential defects before actual production, and the accuracy of this analysis is highly dependent on the quality of the mesh. Among various mesh types, hexahedral meshes offer the highest quality but are challenging to generate automatically. Hexahedral meshes are formed by stacking quadrilateral surface meshes layer by layer, making the construction of quadrilateral meshes critically important. Therefore, this study focuses on developing an automatic unstructured quadrilateral mesh generation technique based on the Paving algorithm. During this process, mesh self-intersections and collinearity are the most common and severe issues, often leading to meshing failure. To address this, the study tested over 1,400 mesh size configurations across 14 IC CAD models and introduced a series of improvements to the overall meshing workflow. These include self-intersection detection and enhancements in three key areas: Mesh modification, local mesh construction, and function merging. Additionally, a multi-model, multi-size automated testing module was developed to improve testing and debugging efficiency, and to validate the stability and computational performance of the proposed technique. The test results demonstrated successful mesh generation for all models, with a reduction in meshing time for most cases. For the model with the longest original build time, the average meshing time was reduced 6.77%. This research successfully enhances the stability and computational efficiency of automatic unstructured quadrilateral mesh generation.

[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, 1991.
[2] R. Cass, “A general three dimensional all quadrilateral surface mesh generation algorithm”, MS. Thesis, Brigham Young University, 1992.
[3] R. J. Cass, S. E. Benzley, R. J. Meyers and T. D. Blacker, “Generalized 3-D paving: an automated quadrilateral surface mesh generation algorithm”, International Journal for Numerical Methods in Engineering, Vol. 39, No. 9, pp. 1475-1489, 1996.
[4] D. R. White and P. Kinney, “Redesign of the paving algorithm: robustness enhancements through element by element meshing”, Proceedings of 6th International Meshing Roundtable, Park City, United States, Oct. 13-15, pp. 323-335, 1997.
[5] P. Kinney, “CleanUp: improving quadrilateral finite element meshes”, Proceedings of 6th International Meshing Roundtable, Park City, United States, Oct. 13-15, pp. 449-461, 1997.
[6] T. J. Tautges, R. R. Lober, C. Vaughan, and J. Jung, “The design of a parallel adaptive paving all-quadrilateral meshing algorithm”, Proceedings of 6th International Symposium on Computational Fluid Dynamics, Lake Tahoe, USA, Sep. 4-8, 1995.
[7] S. J. Owen, M. L. Staten, S. A. Canann and S. Saigal, “Q-Morph: an indirect approach to advancing front quad meshing”, International Journal for Numerical Methods in Engineering, Vol. 44, No. 9, pp. 1317-1340, 1999.
[8] A. V. Skovpen, “Modified algorithm for unstructured quadrilateral meshing”, RFNC-VNIITF, 2004.
[9] R. Miyazaki and K. Harada, “Transformation of a closed 3D triangular mesh to a quadrilateral mesh based on feature edges”, International Journal of Computer Science and Network Security, Vol. 9, No. 5, pp. 30-37, 2009.
[10] M. L. Staten, S. J. Owen and T. D. Blacker, “Unconstrained paving and plastering: a new idea for all hexahedral mesh generation”, Proceedings of the 14th International Meshing Roundtable, pp. 399-416, 2005.
[11] J. Sarrate and A. Huerta, “Efficient unstructured quadrilateral mesh generation”, International Journal for Numerical Methods in Engineering, Vol. 49, No. 10, pp. 1327-1350, 2000.
[12] M. Xinwu and L. Sun, “An automatic approach to constrained quadrilateral mesh generation”, Engineering Computations, Vol. 37, No. 3, pp. 929-951, 2020.
[13] J. F. Remacle, J. Lambrechts, B. Seny, E. Marchandise, A. Johnen and C. Geuzaine, “Blossom-Quad: a non-uniform quadrilateral mesh generator using a minimum cost perfect matching algorithm”, International Journal for Numerical Methods in Engineering, Vol. 89, No. 9, pp. 1102-1119, 2012.
[14] D. Zint and R. Grosso, “A hybrid approach to fast indirect quadrilateral mesh generation”, Numerical Geometry, Grid Generation and Scientific Computing, pp. 281-294, 2021.
[15] D. Bommes, B. Lévy, N. Pietroni, E. Puppo, C Silva, M. Tarini and D. Zorin, “Quad-mesh generation and processing: a survey”, Computer Graphics Forum, Vol. 32, No. 6, pp. 51-76, 2013.
[16] A. A. Al-bana, S. Zahran, S. T. Shahen and M. M. Ahmed, “Performance analysis and validation of different meshing element size used for computational fluid dynamics simulation”, 2023 International Telecommunications Conference (ITC-Egypt), pp.165-168, Alexandria, Egypt, Jul. 18-20, 2023.
[17] 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, Tokyo, Japan, Apr. 21-25, 2023.
[18] J. Y. Lai, P. Putrayudanto, D. H. Chen, Q. Wang, J. Huang, P. P. Song and Y. C. Tsai, “Development of automatic quadratic meshing technique with unequal spreading intervals on outer and inner contours for IC CAD models”, System Innovation for an Artificial Intelligence Era, CRC Press, pp. 11-18, 2024.
[19] 陳定輝，「應用於非結構化四邊形網格建構之輪廓撒點與網格品質改善技術發展」，國立中央大學碩士論文，2023。
[20] 王齊，「四邊形網格自動建構之多尺寸輪廓撒點技術研究」，國立中央大學碩士論文，2024。
[21] 蔡敬崙，「四邊形網格自動建構之網格品質改善研究」，國立中央大學碩士論文，2024。
[22] 李奇勳，「非結構化四邊形網格自動建構研究」，國立中央大學碩士論文，2024。
[23] 戴宇辰，「混合結構化與非結構化四邊形網格之自動化區域劃分方法發展」，國立中央大學碩士論文，2023。
[24] 吳弈翰，「IC CAD 模型之混合結構化與非結構化四邊形網格自動建構技術發展」，國立中央大學碩士論文，2024。
[25] 梁秉傑，「不同網格類型對於IC CAD模型模流分析之影響探討」，國立中央大學碩士論文，2024。
[26] Rhinoceros. Avaliable: https://www.rhino3d.com. [Accessed May 20, 2025].
[27] OpenNURBS. Avaliable: https://www.rhino3d.com/tw/features/developer/ opennurbs. [Accessed May 20, 2025].
[28] B-rep structure. Available: https://developer.rhino3d.com/guides/cpp/brep-data-structure. [Accessed May 20, 2025].
[29] 曾思齊，「IC CAD模型混合結構化與非結構化四邊形網格自動建構技術優化」，國立中央大學碩士論文，2025。
[30] 蘇敬勛，「表面網格自動建構之網格品質優化方法發展」，國立中央大學碩士論文，2025。
[31] 劉弘堅，「四邊形網格自動建構之輪廓撒點方法與四邊形網格轉三角形網格之研究」，國立中央大學碩士論文，2025。
[32] 賴景義，翁文德，逆向工程理論與應用，全華科技圖書股份有限公司，臺北市，2004。