堤防沉降计算——单元扭曲不知如何解决

waternj2008

如附件，做了一个堤防沉降的计算。由于几期填筑分层很细，加上有放坡，所以导致剖分很麻烦。虾米我做了半天处理，把所有部分划分成可以利用结构化网格进行剖分，但是，利用单元检查中的“分析检查”总是出现黄色的区域。（Number of elements :  86062,   Analysis errors:  0 (0%),  Analysis warnings:  22502 (26.1463%)）
提交计算后也老是出现单元扭曲的警告，然后就半天不动直至不收敛而跳出。
请教各位大侠，如何处置！

waternj2008

上传个细部放大图。
模型四周施加了侧向位移约束，底部施加了三向位移约束，模拟沉降为主。

zsq-w

模型中网格质量太差。在mesh步verify一下 你会发现很多呈黄色。

lz应该重新划分网格，画完verify下没有黄色和红色出现就表示这个网格质量已经很好了。

waternj2008

参照版主的建议，尝试对网格进一步做处理，最好的情况是能将Analysis warnings的比率缩小到3%左右，但是还是不是很好。出现的黄色区域仍然集中在需重点开展沉降分析的堤防顶部。 尝试用扫描网格，甚至会出现粉红色错误区域。 在论坛中搜索也未见如此这样多级放坡的模型。哎，无法了，只好将CAE文件上传，请诸位大侠指点迷津。多谢！ 另：作为新手，向在Aba里畅游的高手们致敬！

论坛附件大小限制不是256kb吗，怎么我153kb的文件就是传不上来？？？
原来是100Kb，哎只好拆成两个发了。

zsq-w

因为lz模型的图层很多，所以首先必须向lz确认两个问题才能看再如何划分这个网格：

1  这些分割（partition）是不是全都是必须的？ 因为这么分割导致有些尖角有16度到29度，对于这么小的角，画出来的网格必然有些四边形的最小角度是16度了，那种都是不好的网格。
  去掉一些分割可以改善这个问题。But，我对土木不熟悉，也不知道是不是允许删除一些partiton的
2  是不是只能划六面体？允许四面体存在么？ 有些区域划出四面体可能可以得到更好质量的网格。

一个建议：上部的几个土层，lz都用的是structure划分，可以尝试一下sweep。lz也可以感受下sweep和structure划分的区别。（一般的，画网格优先采用strcuture的，不过有时候也可以考虑下sweep和bottem-up）

waternj2008

感谢版主热心回答。
那些分割有的是不同土层的分层（下覆4层-分别是粉砂、粘土、粉砂、粗砂），有的是不同阶段人工填土形成的分界。特别是顶部的填土区域，由于是抛填，所以形成了很多个斜坡——每一级都是一个梯形故分割较多。
由于分析的时候用到了自己开发的材料模型，所以理论上最好用四边形单元（空间八节点六面体单元），尝试了扫描网格sweep、也尝试了缩减积分的20节点单元仍不收敛。
痛苦中！！

zsq-w

能理解用六面体单元。

对于这么对尖角的几何，画出的网格只能凑合用了。也可以考虑下六面体+四面体混合网格。

waternj2008

采用Explict，加沙漏控制（stiffness），计算通过。单元仍用的C3D8R，网格密度稍微增加了些，线性缩减积分。

waternj2008

还是上面的那个问题。
堤防分级堆载的过程中，后续堆载经常不收敛，故此，想采用了restart的方式，不断调整荷载步和增量步。
但是，在利用restart计算的时候，出现了如下错误提示：

***ERROR: in keyword *NSET, file "Taicang-geo-E-sandFill-rec-sand2.inp", line
           5: The input file for the original analysis does not contain parts
           and assemblies.  Upon restart, only assembly level sets may be
           defined.  The set name given by the NSET/ELSET parameter should not
           contain any periods which imply an assembly of part instances.  The
           unrecognizable name is: ASSEMBLY.PART-1-1._PICKEDSET22
***NOTE: DUE TO AN INPUT ERROR THE ANALYSIS PRE-PROCESSOR HAS BEEN UNABLE TO  INTERPRET SOME DATA.  SUBSEQUENT ERRORS MAY BE CAUSED BY THIS OMISSION

即在restart之前step里定义地节点、单元集合在restart里找不到对应的信息了-unrecognizable name is: ASSEMBLY.PART-1-1._PICKEDSET22。

有没有哪位大侠碰到类似的情况过?
问题产生的原因可能与如下过程有关：为了地平衡，用了应力导出后作为初始应力条件（initial conditions） 叠加的方式进行，故此，关闭了块体显示——mdb.models['模型名称].setValues(noPartsInputFile=ON)。但是不用restart就没有上述问题，可以一路计算下去？