PFC模拟CPT（标准灌入试验）

这里的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

应该没有错，各位朋友可以参考一下。