多孔介质的地应力平衡

根据单向固结理论，有效应力=总应力-孔隙压应力。在非饱和问题中，ABAQUS采用的是：有效应力=总应力-饱和度*孔隙水压力。又由于ABAQUS对孔压正负号的规定和土力学中一致，而对应力正负号的规定和土力学中刚好相反，故上面公式中的“-”号要改为“+”号。对《 abaqus在土木工程中的应用》一书中的这个例子，土体顶面的总应力是为零的，但孔压为-10kpa，查-10kpa对应的饱和度得0.9789，再用上面的公式计算就得到了-9.789kpa。个人感觉不要死套公式，对每一点都分别算出总应力和孔隙水压力，再按上式计算有效应力。最后一点：《abaqus在土木工程中的应用》一书中给出的初始有效应力的计算公式有点问题，手册上面说的很详细，手册里面是把土体分为完全干燥区，半饱和区和完全饱和区三部分来考虑的。把手册上相关的部分传上来，希望对大家有所帮助。