首先感謝中南王同學提供的滾刀圖形，是用犀牛軟件調整網格，具體的我也不太明白，這里講一下導入進來后的一些操作。

對于三維滾刀，有前進速度，也就是盾構機轉動的線速度。還有切入速度，也就是盾構機頂入的速度。還有滾動速度，是滾刀轉動的角速度。轉速6.28rad/s，前進速度0.04，切入速度是前進速度的0.2。

速度比較簡單，難點在于根據滾刀的位置調整滾刀的轉動中心。

我寫了一個dat用于處理導入的geometry,里面主要的功能是得到geo的坐標，還有擴大縮小geo，移動geo的基本操作。這里滾刀是7寸,直徑0.423m。

將geo移動到合適的位置后，生成滾刀進行加載，具體看代碼體會，下面放上處理的結果。

總的運動圖上面可以看到。

給出俯視圖，顆粒按豎向位置染色。

紅色是表面的，再下面是橘黃色和黃色。滾刀沒有顯示，可以看到滾刀的侵入過程。

裂紋圖

這里可以看到滾刀侵入過程中產生的裂紋。

這里將試樣切片了，可以看到滾刀侵入過程中力的變化。豎向力會略大于橫向力。

綠色是豎向力

扭矩這里寫了一下，沒有驗證過，大家可以和實際對照一下。