为什么初始应力图会变成这样？

geostability

这种情况遇到的很多了,你单独考虑重力的情况下,计算一下,就不会有你的现象了.因为这个施加的szzsxxsyy有一定的偏差的,

who_is

首先初始化应力的时候应该是 grav*dens，你用了10 ，可能计算会多余一点。从你的图中感觉好像两个地方密度不一样。这是solve以后的？

qinjianshe

最好把程序传全点，不然很难帮你调试指出问题的

chunhuo420

我给你提两个和主题无关的建议：
1）set grav 0 0 -9.81
ini dens 3950
ini szz -7.9e6 grad 0,0,3.95e4
ini sxx -2.63e6 grad 0,0,1.32e4
ini syy -2.63e6 grad 0,0,1.32e4

用到ini初始平衡的时候要和重力计算的相同才有用，
可以免去重力平衡的时间，你这里明显是不相等的

2）group tun4 range group cave3 not group tun3 not group cave2 not group tun2 not group side1 not group tun1 not group cave1 not group lower_cave not

分组的时候不需要这么做的，只要在你要group 的块体建模后面 直接加上就可以了  （注意：不是在建模，并回车后group，而是直接group ）

希望对你有点帮助！

geostability

本想帮你计算一下，你的模型太大了，107000个单元，我的电脑显示内存不足，哈哈。加不加ini，影响重力平衡的时间。可以直接用grav了，这个一般计算要准确采用的,同意chunhuo420 的观点。

freddie_1999

网格不一致，请regen或attach

who_is

ini 只是减少达到初始平衡的时间。我说得9.8时因为你把重力加速度设置为9.8，而grad 0 0 3.95e4种3.95e4明显是根据加速度为10来算的，所以并不平衡。
而且无论如何设置，模型表面处的szz应该为0。

qjbscc

我觉得是网格化分密度的问题

feiyang52013

我不知道什么问题，