PFC利用体积不变法模拟不排水双轴

    这个算例可以看出离散元在模拟散体颗粒力学性质方面的优越性。因为我自己硕士大论 文 做的是有限元循环塑性本构，用的是下负荷面模型实现的循环特性，主要的模型参数多达12个。而在离散元中，我实际上的材料参数只有emod、kratio、fric这三个，当然如果孔隙率、级配也算的话就比较多了。下面讲一下实现不排水循环双轴的实现思路。

    首先，不排水采用的方法为体积不变法。

    dx=(-A/y^2)*dy

用这个可以根据竖向速度计算横向速度，速度的方向如下图。

如下图为计算过程中面积的变化，有少许变化，但是这个量级可以看一下，很小。

如下图为计算过程中的应力路径变化。

这里已经算了一天了，基本上可以看出比较经典的滞回圈和液化现象。当p为0基本上可以认为是液化。

应力应变曲线也是比较经典的。这个可以自行和砂土实验对比一下。

这里记录了墙上力的变化，可以根据正应力的损失去计算静水压力。

注意：此处循环实现用的是匀速加载，监测应力到达一定值后反向。

代码分为：成样-加围压-循环加载