在空間中生成剛性顆粒（注意是剛性顆粒）有下列幾種方法：

    1.修改關鍵字，構建粒子生成器模型生成隨機分布剛性顆粒

   2.使用python語言直接在ABAQUS中生成顆粒，并進行剛體綁定，使其成為剛性顆粒，或者直接生成解析剛體或離散剛體。

    方法1生成顆粒的隨機性較好，操作簡單。方法2直接在ABAQUS界面生成顆粒，當所需顆粒數量以萬為計量單位時，在前處理界面時就會卡死，對顯卡要求極高。因此，在僅考慮到這些弊端情況下，就已經使研究人員頭皮發麻，無從下手。

    在一些特定應用場合下，比如所需顆粒數量數以萬計，我們只能采用方法1生成顆粒，但我們不僅僅是需要顆粒，還需將這些顆粒與其它模型進行耦合求解計算，這個時候粒子生成器就會有局限性。此外，考慮到顆粒在空間中排布的多樣性，比如最典型的高斯分布，那么粒子生成器很難做到一步到位生成所需分布特征的顆粒。

    我們今天介紹的通過python編程操作ABAQUS輸入文件生成PD3D單元顆粒，其可操作性更強，我們可以不采用粒子生成器內部定義的隨機算法生成顆粒，用戶可以根據需求自定義顆粒分布算法，以契合實際工況。此外，可省去粒子生成顆粒的分析步，直接進行工況建模求解計算。

    本貼只是個人興趣，只提供思路，不提供源碼，用戶需了解ABAQUS的inp文件的書寫規則、python操作文件語法和生成顆粒的底層邏輯(分布模型)。感興趣的可以私信，提供編寫思路。

    下面我們采用這一方法生成直徑2mm、3mm、4mm和5mm的混合顆粒，數量為1000。具體生成結果如下圖所示。