光滑粒子流體動力學(SPH)及DYNA相關問題的案例
SPH光滑粒子流體動力學方法入門
一、SPH法簡介
SPH方法的理論思想源于粒子法，在SPH方法中，所有單元都是由占有獨立空間且具有獨立質量的有限個粒子近似表示的。它本質上基于拉格朗日方法，但又不同于邊界元法和有限元法，無需借助網格。因此，SPH函數f(X)的積分表達式可表示為支持域內所有粒子疊加求和的離散化形式。
SPH算法現已發展成為比較成熟的計算三維物理問題的有效方法并被推廣到流體動力學、空氣動力學、切削，高速碰撞大變形等力學分析的各個方面。
二、SPH粒子的生成方式
1.直接生成
在ls-prepost中用SPHGEN命令直接生成，生成方式選擇box，輸入SPH粒子所占的空間坐標。
XYZ 3個軸的起始坐標點，坐標終點分別填入，該坐標可以用全局坐標或自定義的坐標。一般來說在lspp中的坐標操作均以全局坐標為主，局部坐標運用的很少，離子的密度填在numXYZ處，就是XYZ軸方向有幾個粒子。這個是生成的在100mm3的空間內分布的例子。用這種方式建模不方便，在模型復雜情況下不容易進行操作。一般不采用這種方法。
2.通過網格生成
通過網格生成的方式就是將已劃分網格的物體轉化為SPH粒子，通過solid nodes 或solid centry方式將網格的節點轉化為粒子。
這里需要注意，SPH粒子如果通過網格生成，網格一定要是六面體網格才行，普通的四面體網格雖然能生成粒子，但是在計算的時候由于粒子分布不均勻會導致計算結果不準確。
一般我在workbench中做前處理，將要轉化成SPH粒子的單元用sweep方式劃分網格。
展開 ABAQUS中的光滑粒子流體動力學 ( SPH ) 方法
有限單元轉化成SPH粒子，可參考本人帖子：
https://www.yqgqt.org.cn/post/435476
概述
光滑粒子流體動力學（SPH）方法是一種無網格數值方法。通常的有限元分析中需要定義節點和單元，而該方法用點的集合來描述給定的部件，無須定義單元。在SPH法中這些點通常被稱為粒子或擬顆粒。
圖1中對比了兩種方法。兩個離散模型描述的都是瓶子里裝的液體。左邊的模型是由流體占據的傳統四面體網格；在右邊，同樣的流體體積是由離散點的集合表示的。注意，后者情況下沒有網格連接這些點（粒子），它們無需像左邊傳統的有限元定義多節點單元從而保持連通性。在ABAQUS中除了直接定義SPH粒子外，還支持先定義傳統的連續單元，然后在分析開始時或在分析過程中將單元網格自動轉換成粒子。
圖1 有限單元和SPH顆粒的分布
光滑質點流體動力學（SPH）是一種純拉格朗日方法，它允許通過插值性質直接離散化一個給定的連續性方程組而無需定義空間網格。SPH的主要優勢是無固定網格，對于流體流動、結構大變形和自由表面等難題，該方法處理得相對自然恰當。
SPH的核心并非基于在壓縮中彼此碰撞或在張力作用下表現出粘性行為的離散顆粒(球)。相反，它是將連續偏微分方程組巧妙離散化的一種方法，這一點與有限元法非常相似。SPH利用插值來近似域中任意點的場變量值，粒子某個變量值通過對相鄰粒子對應的值疊加求和來近似，這些粒子以下角標j來區分，其核函數為W（非零）如下。
SPH的核心是核函數，它可以被理解為一種在一定光滑長度h范圍內其他臨近粒子對研究粒子影響程度的權函數，如圖2所示。其中，光滑長度h決定了對某個點的插值產生影響的粒子數目。
展開 SPH(光滑粒子流體動力學)-模擬水蝕 ￥10
</p><p>目前采用SPH方法實現單個水平沖擊金屬涂層基體的過程，具體詳細步驟大家可以自行去研究cae和inp文件，如果有不明白的地方，可</p><p>在此感謝Usim大佬的支持，大家可以搜索會員名字 Usim ，去他的主頁看看，不是一般的NB，動力顯示分析的大手。</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/201911/f83382c467a74f39a6ad8326a928ae9c.gif" title="SPH.gif" alt="SPH.gif" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/201911/f83382c467a74f39a6ad8326a928ae9c.gif?image_process=/format,webp/quality,q_40/resize,w_400" data-pc-src="https://img.jishulink.com/upload/201911/f83382c467a74f39a6ad8326a928ae9c.gif?
展開 