等效线性黏弹性本构的应用

有个困惑就是，做试验可以获得不同围压下土体的最大剪切模量，以及剪切模量比与剪应变，阻尼比和剪应变的曲线，然后根据不同围压与最大剪切模量的关系就知道了公式中的k和n，做模拟的时候，在材料属性输入k,n,v,w，关键字中输入各土体单元的震前围压，剪切模量比，阻尼比，最大剪应变。如果假定震前围压为100，那在迭代过程用所用到的剪切模量比与剪应变及阻尼比的曲线就是100kpa所对应的曲线。但是如果考虑震前围压，就是先做静力分析求出各单元的有效应力作为关键字输入中的第一列，那这样的话基本一层土是一个应力，也就是一层土一种围压，一种围压对应一个最大剪切模量和关系曲线，迭代的时候不可能取每层土对应围压的下土体的关系曲线，那么要用哪个围压下的关系曲线？对于正常固结土，最大的围压是密度*g*h，做土力学实验获得以上关系曲线是根据土体深度来加的围压，那如果现在土体密度是2，模拟土体厚度为60m，最大围压就是1200kpa，迭代的时候要用围压1200kpa对应的关系曲线吗？可是考虑静应力之后，每层土一个围压，只有最下层土体围压才是1200