無網格法

有限元（finite element method 簡稱 FEM）在學術和工業已經有了很多的成就 ^[4]  。但是，有限元也存在很多短板。其中最嚴重的之一就是有限元的計算結果嚴重依賴於網格的劃分質量 ^[4]  。高質量的有限元網格需要相關專家長時間的參與，尤其是對複雜的集合圖形 ^[5-6]  。即使使用高質量的網格，有限元模型也經常會在大形變的模擬中失敗 ^[7]  。

無網格法（meshless methods）的誕生，就是來消除有限元（finite element method）的短板 ^[4]  。最早的有關無網格法的文獻是由劍橋大學的L.B.Lucy ^[8]  和劍橋大學的R.A.Gingold, J.J.Monaghan ^[9]  分別於1977年提出。他們介紹了光滑粒子流體動力學(Smoothed particle hydrodynamics，簡稱SPH) ^[8-9]  。這兩篇著作的初衷是解決天體物理學（astrophysics）的問題，後來，SPH多用於流體動力學（fluid dynamics） ^[10]  。SPH法的計算是基於強形式（strong form），其餘的無網格方法都是在1990年代提出的，他們是基於弱形式（weak form）計算的。最開始主要是用於計算材料力學的 ^{[2]  [4]}  。Galerkin法（Element free Galerkin，簡稱EFG）是第一個基於弱形式的無網格法，於1994年被T. Belytschko等人提出。 ^[11] 

今天，無網格法廣泛的引用在流體 ^[12]  ，生物力學 ^{[5]  [13-14]}  等研究中，並取得諸多成果。

該法大致可分成兩類：一類是以Lagrange方法為基礎的粒子法（Particle method），如光滑粒子流體動力學（Smoothed particle hydrodynamics，簡稱SPH）法，和在其基礎上發展的運動粒子半隱式（Moving－particle semi－implicit，簡稱MPS）法等；另一類是以Euler方法為基礎的無格子法（Gridless methods），如無格子Euler/N—S算法（Gridless Euler/Navier-Stokes solution algorithm）和無單元Galerkin法（Element free Galerkin，簡稱EFG）等。

無網格方法可以方便地利用座標點計算模擬複雜形狀流場計算，但不足之處是在高雷諾數流動時提高數值計算精度較困難。

無網格方法中比較常見的還有徑向基函數方法（Radious Basis Function），主要使用某徑向基函數（如（MQ）f(r)=r^5）的組合，來逼近原函數。吳忠敏院士在這方面有比較突出的工作。以上方法中，無網格伽遼金法成為目前影響最大，應用最廣的無網格計算方法，現有的 LS-dyna，Abaqus，Radioss等商業軟件都加入了該方法的計算模塊。