本研究以「網格演化有限元素分析法」為基礎，主要發展高效能三維大地工程有限元素分析流程。傳統上，以有限元素分析在建構網格模型時，需考慮每階段分析上之幾何形狀，導致網格分割上的困難，並且於三維分析時可能產生大量的自由度量，因而需要龐大的求解計算能力。使用網格演化法分析，只需要針對現階段分析進行建模，使網格模型產生較為容易，並降低自由度量，進而提升分析效能。
在本研究中提出針對分析模型中不同材料性質之區域，分別進行其應力場傳遞之動作，減少不同材料區域間的互相影響。此外，並針對前一階段分析後之輪廓擷取方式，提出兩種演算法，使得在二維或三維分析下得以自動化執行網格演化有限元素分析法。
本文以大地工程分析為應用方向，並結合ABAQUS有限元素分析軟體，實作出本研究提出的分析流程，並與傳統有限元素分析之結果比較，以驗證此分析流程之可行性和正確性。本研究針對典型三維大地工程深開挖及隧道開挖分析，以網格演化法進行分析，所需計算時間分別為傳統有限元素法的12%和35%，可充分展現本分析方法之效能。

This study, based on “mesh evolving method”, develops a high-performance procedure for three-dimensional finite element analyses in geomechanics. Compared to conventional finite element analyses, the mesh-evolving method simplifies mesh-construction and reduces demand of computation power by meshing each excavation stage separately. Such meshing strategy reduces both meshing-difficulties and degrees-of-freedom, and results in improved computational efficiency.
This study improves mesh-evolving method in two aspects. 1) In mesh-to-mesh solution mapping, stresses in zones of different material properties are transferred independently to avoid interferences between different zones. 2) Furthermore, two algorithms are proposed in order to get geometrical profile from the calculation result of the previous step, and these algorithms enable 3D mesh-evolving method.
The mesh-evolving procedure was implemented with C++ programming language and commercial finite element software: ABAQUS. The proposed procedure is verified by comparing its calculation results with conventional finite element calculations. For typical 3D deep-excavation and tunneling analyses studied, the mesh-evolving method reduces the computation time by 88% and 65% compared with the conventional finite element analyses.

“ABAQUS實務入門引導”，愛發股份有限公司編著，全華科技圖書股份有限公司印行(2005)
宋俊宜(2007)，”有限元素網格對開挖分析結果影響之估計方法初探”，國立台灣科技大學營建工程所碩士論文
潘茂森(2006)，”網格演化有限元素分析於大地工程之應用” ，國立台灣科技大學營建工程所碩士論文
劉昌榮(2006)，”CPIM,改進之元素釋放法在工程結構分析之研究” ，國立中山大學海洋環境及工程學所碩士論文
“ABAQUS Analysis User’s Manual Version 6.5”, ABAQUS, Inc.
Agnantiaris, J. P. , Polyzos, D. , and Beskos, D. E. (1996) “Some studies on dual reciprocity BEM for elastodynamic analysis”, Computational Mechanics, 17:270-277
Christos Vrettos (1998). “THE BOUSSINESQ PROBLEM FOR SOILS WITH BOUNDED NON-HOMOGENEITY”, International Journal For Numerical And Analytical Methods In Geomechanics, 22:655-669
Hsieh, Y.-M. (2004). “Parallel Computation in Efficient Non-Linear Finite Element Analysis with Applications to Soft-Ground Tunneling Project” ,Massachusetts Institute of Technology, Cambridge.
Hardy, R. L. (1990). “Theory and applications of the multiquadrics – biharmonic method(20 years of discovery 1968-1988)” ,Computers and Mathematics with Applications, 19:163-208
Krysl, Petr. and Belytschko, Ted. (2001). “ESFLIB: A library to compute the element free Galerkin shape functions”, Computer Methods Applied Mechanics Engineering 190 (2001) 2181-2205.
“LAPACK Users' Guide Third Edition”, http://www.netlib.org/lapack/lug/index.html
Liu,G. R. , and Gu, Y.T. (1999). “A point Interpolation Method”, in Proc. 4th Asia-Pacific Conference on Computational Mechanics, December, Singapore, pp.1109~1014.
Liu, G. R. (2003). “Mesh free methods :moving beyond the finite element method” , CRC Press
Powell, M. J. D. (1994). “The uniform convergence of thin plate splines in two dimensions”, Numerische Mathematik, 68(1):107-128
Wang, J.G. , and Liu, G. R. (2002). “On the optimal shape paraeters of radial basis functions used for 2-D meshless methods”, Computer Methods In Applied Mechanics And Engineering, 191:2611-2630
Selvadurai, A. P. S. (2007).“The Analytical Method in Geomechanics”, Department of Civil Engineering and Applied Mechanics, McGil University.