求助砂土液化的模拟

benbenmoon

最近一直在模拟一个震动台试验。主要是参考刘洋老师的博士论文（砂土液化破坏的细观力学机制与数值模拟）里的方法，用周期边界来模拟左右边界，通过底部的墙施加动荷载，然后想试试能不能得到液化状态。。。

但做出来的结果很不理想。。。发个帖子散尽家财求助达人。。。
遇到的主要问题如下：

有效应力/总应力怎么描述？
在模型里不同深度设置了测量圆，可以用measure 的s22作为竖直向有效应力吗？这个记录的是有效应力还是总应力？。。。还是应该循环contact_head再计算，这样算出来的又是总应力还是有效应力？。。。

另外，我设置用fc_pre(m,n)记录孔压，在加上动荷载之后孔压变化很乱，这个是不是因为流体的设置有问题。



流体部分设置如下：
    fluid model 0 _lx 0 _ly size 3 5
    fluid prop den 1e3 vis 1e-3
    fluid bou slip xl
    fluid bou slip xu
    fluid bou pres yu
    fluid bou pres yl
    fluid bou vel 0 0 yl
    fluid set grterm on
    fluid set buoterm on
    fluid set visterm on


命令流如下：

wzx2005

楼主问题解决了吗 向您请教 谢谢