這里的CPT用wall來模擬，wall可以記錄錐尖的法向力和側摩阻力。

為了方便觀察土體變形，這里對土體進行了染色。上圖的紅藍色便是染色的結果。

這個算例比較簡單，直接灌入便可以，注意灌入速度，我這里只是測試代碼，所以速度比較大，而且顆粒數不多，得到的力的結果不是很理想，這里就不放出來了。

下面給出位移動圖

力的監測代碼：

def jiance
    faxiang=math.sqrt(wall.force.contact.x(wpl)*wall.force.contact.x(wpl)+wall.force.contact.y(wpl)*wall.force.contact.y(wpl))
    cemozuli=wall.force.contact.y(wpu)
end

應該沒有錯，各位朋友可以參考一下。