关于数值计算中初始自平衡问题

qinjianshe

     Flac3D模拟隧道开挖，当考虑流固耦合时，其采用的有效应力在计算中是通过总应力Szz,Sxx,Syy与空隙水应力u求得，但是我们知道，ko值不能应用于总应力中侧向应力Sxx,Syy的计算，所以在总应力计算中，常规的ko值求出的侧向总应力有问题，从而导致侧向有效应力计算有误。
  请高手对此问题发表意见，我自己是反过来求总应力的。

cafe0519

现在就加分,欢迎两位多来这边聊聊,提提意见。谢谢。

fengjun4318

对于考虑地下水的情况：
由于符号显示问题，符号说明如下，s¢v-竖向有效应力;s¢h-横向有效应力;sv－竖向总应力；sh－横向总应力。
s¢v = (sv – u)  (竖向有效应力)
s¢h = (sv – u) Ko  (横向有效应力)
sh = s¢h  + u  = (sv – u) Ko＋u (横向总应力)，
如果认为s¢h = sv Ko－u ，sh =svKo就错了。
在数值计算中，可以按以下办法实现正确的横向总应力计算：
1．  单独计算地下水引起的空隙水压力。即暂时将水位以上的土体单元kill掉，水位以下的土体重度设为水的重度，侧压力系数设为1.0，计算出静水压力。
2．  计算总应力。将第一步计算出的应力作为初始应力，此时将水位以上的单元激活，重度使用天然容重，水位以下土体重度采用浮容重，再给定相应的伯松比，即可计算出总应力。
两次计算均采用相同网格。

qinjianshe

多谢cafe0519加分。
楼上老兄对于有效应力的计算没问题，关键是在flac3d计算中，其现给出总应力后软件自身求有效应力。关于单于死活是一个好建议，但不知flac3d中如何实现。请老兄指教