SPH方法在大變形情況下的計算效果還是可以的，Abaqus集成SPH算法也有幾個版本了，每個版本或多或少都會對SPH進行改進。用Abaqus的SPH算法做了水動力中經典的潰壩算例。

SPH算例和其他有限元算例很類似，需要注意的點有：

1.   在Abaqus中，粒子并不是在初始時刻就存在的，只有在滿足一定條件的情況下才會由網格轉化為粒子。轉化為粒子的條件有時間條件、應變條件等，比如選擇時間條件，可以指定開始轉化的物理時間，方法是：在指定物體的網格類型時，選擇Explicit，與Standard不同的是Explicit單元庫中的單元可以選擇轉化為粒子（“Conversion to particles”），轉化規則有Time、Strain和Stress，如果轉化門檻選擇為0，則在開始時刻就轉化為粒子。

2. 接觸問題：直接采用General Contact即可，其他接觸形式包括面面接觸、自接觸等，試了一下效果并沒有更好。接觸特性（“Contact Property”）選擇一個法向“Hard Contact”和切向“Frictionless”即可。

3. 邊界條件：由于粒子是在提交計算后由網格轉化來的，所以邊界條件貌似無法直接施加到粒子上。潰壩的算例是采用了一個矩形盒子作為邊界，矩形盒子材料是鋁，由于水域的尺寸太大（100m），盒子的厚度指定的太小（0.01m），還需要將盒子的四周固支，否則盒子受力很大就會發生很大變形。

注意：本算例基于Abaqus 2018版本，較低版本可能不一定適用，以下內容為cae文件和inp文件，給需要的朋友，無其他技術討論，其實上面1、2、3點就是主要的技術點了。