利用軟件PFC3D建立砂石料的直剪模型，采用墻體的伺服功能對剪切盒的上壓板施加豎向壓力，模擬了在不同頂部壓力作用下的直剪過程。

100kPa頂部壓力時球顆粒接觸力及剪切力-位移曲線如下：

200kPa頂部壓力時球顆粒接觸力及剪切力-位移曲線如下：

300kPa頂部壓力時球顆粒接觸力及剪切力-位移曲線如下：

不同頂部壓力時剪切力-位移曲線對比結果：

剪切前后模型對比如下：

剪切前：

剪切后：

剪切后上下剪切盒銜接部位展示：

建模思路及代碼如下：