关于三维有接触有其他部件的地应力平衡，求大家指教

abaqus的新手

模型情况：地层中有其他部件，与土体有接触，莫尔库伦本构，版本6.14。已知不能采用自动平衡地应力的方法平衡有接触的情况。
方法1:采用predefined field，也就是之前的initial conditions type=stress，geostatic。step设置为geostatic，type选择fixed，非对称求解器。结果不收敛，更改了土体参数仍不收敛，怀疑该方法是否能处理又部件有接触的情况。
方法2:自动平衡，题主做的管土接触，平衡时不加管道，在接触面上设置土体的位移边界条件，在后续分析中激活管道，step设置成automatic，非对称求解。结果是可以收敛，但结果不近人意，怀疑该结论的真实性。
方法3:采用生死单元，平衡时，杀死部件，部件所在位置由土体代替，采用自动平衡方法，在第二步remove该部分土体，激活部件进行后续分析。这个方法成功率相当高，但是总感觉离实际工况有差距。
方法4:将面面接触改为tie约束，再采用自动平衡后导入odb的方法。只成功过一次，感觉收敛程度基于参数。
以上就是这段时间以来试过的方法，刚刚接触abaqus以来两个月，感觉身体被掏空。总的来说，仍没有一个令我满意的方法，在这里，欢迎大家一起来讨论一下，交流一下方法

早上好888

楼主，我们做的应当是同类型的计算，我的是三维的，有接触的，地下是基岩，上覆堆积土，另外还考虑孔隙水压力，做流固耦合，但是计算初始地应力时很难收敛，把弹性模量调的很高了也不行

沙漠胡杨806

怎么样了，做好了么？问下边界条件怎么设置

yqhuodian

Abaqus土木仿真模拟工作室 版主qq443941211 欢迎加入学习

abaqus的新手

早上好888 发表于 2016-9-1 11:02
楼主，我们做的应当是同类型的计算，我的是三维的，有接触的，地下是基岩，上覆堆积土，另外还考虑孔隙水压 ...

接触改成tie试试

abaqus的新手

沙漠胡杨806 发表于 2016-9-3 15:12
怎么样了，做好了么？问下边界条件怎么设置

边界条件跟着文献资料里面来啊

沙漠胡杨806

这是什么原因，求解答
Error in job Job-1: in keyword *INITIALCONDITIONS, file "Job-1.inp", line 1970374: Unknown assembly node set ALLNODES
Error in job Job-1: in keyword *INITIALCONDITIONS, file "Job-1.inp", line 1970380: Unknown assembly element set ALLELEMENTS
Error in job Job-1: AN INITIAL CONDITION HAS BEEN SPECIFIED ON NODE SET ASSEMBLY_ALLNODES.  THIS NODE SET IS NOT ACTIVE IN THE MODEL
Error in job Job-1: AN INITIAL CONDITION HAS BEEN SPECIFIED ON NODE SET ASSEMBLY_ALLNODES.  THIS NODE SET IS NOT ACTIVE IN THE MODEL
Error in job Job-1: AN INITIAL CONDITION HAS BEEN SPECIFIED ON ELEMENT SET ASSEMBLY_ALLELEMENTS - THIS ELEMENT SET IS NOT ACTIVE IN THE MODEL
Error in job Job-1: 965820 elements have missing property definitions. The elements have been identified in element set ErrElemMissingSection.
Error in job Job-1: ELEMENT SET _Whole_Model_Gravity_Elset HAS NOT BEEN DEFINED
Job Job-1: Analysis Input File Processor aborted due to errors.
Error in job Job-1: Analysis Input File Processor exited with an error.
Job Job-1 aborted due to errors.

abaqus的新手

沙漠胡杨806 发表于 2016-9-8 23:42
这是什么原因，求解答
Error in job Job-1: in keyword *INITIALCONDITIONS, file "Job-1.inp", line 19703 ...

你这是导入地应力文件吗？感觉是文件节点不对，用的什么方法？

沙漠胡杨806

谢谢，问题已解决

yuzhenshuai

早上好888 发表于 2016-9-1 11:02
楼主，我们做的应当是同类型的计算，我的是三维的，有接触的，地下是基岩，上覆堆积土，另外还考虑孔隙水压 ...

你的问题解决了吗

willmutter

个人比较赞同方法3，土体和桩基的相互作用毕竟没那么强，主要把握土平衡