abaqus接触区域选择不恰当导致 domain decomposition不均匀及计算时间长的问题
新建了一个仿真模型,跑的时候异常缓慢,查看sta文件发现domain的划分很不均匀,domain 16 占比八十二,其他剩下的个加起来只有百分之十几。具体如下。
DOMAIN DECOMPOSITION INFORMATION
DOMAIN 1 HAS 1772 ELEMENTS AND 2696 NODES with WEIGHT 1.1510
DOMAIN 2 HAS 2215 ELEMENTS AND 3114 NODES with WEIGHT 1.1512
DOMAIN 3 HAS 3521 ELEMENTS AND 4516 NODES with WEIGHT 1.1513
DOMAIN 4 HAS 2383 ELEMENTS AND 3316 NODES with WEIGHT 1.1512
DOMAIN 5 HAS 3520 ELEMENTS AND 4548 NODES with WEIGHT 1.1510
DOMAIN 6 HAS 1963 ELEMENTS AND 2783 NODES with WEIGHT 1.1510
DOMAIN 7 HAS 2559 ELEMENTS AND 3497 NODES with WEIGHT 1.1514
DOMAIN 8 HAS 2358 ELEMENTS AND 3289 NODES with WEIGHT 1.1512
DOMAIN 9 HAS 3521 ELEMENTS AND 4477 NODES with WEIGHT 1.1512
DOMAIN 10 HAS 2014 ELEMENTS AND 2862 NODES with WEIGHT 1.1512
DOMAIN 11 HAS 2518 ELEMENTS AND 3408 NODES with WEIGHT 1.1510
DOMAIN 12 HAS 3521 ELEMENTS AND 4455 NODES with WEIGHT 1.1512
DOMAIN 13 HAS 2360 ELEMENTS AND 3240 NODES with WEIGHT 1.1510
DOMAIN 14 HAS 2516 ELEMENTS AND 3397 NODES with WEIGHT 1.1513
DOMAIN 15 HAS 2015 ELEMENTS AND 2804 NODES with WEIGHT 1.1511
DOMAIN 16 HAS 100248 ELEMENTS AND 134172 NODES with WEIGHT 82.733
查阅帮助手册Parallel Execution in Abaqus/Explicit这一节,里面提到:
A contact pair cannot be split across parallel domains, but separate contact pairs are not restricted to be in the same parallel domain. A contact pair that uses the kinematic contact algorithm requires that all the nodes associated with the involved surfaces be within a single parallel domain and not be shared with any other parallel domains. A contact pair that uses the penalty contact algorithm requires that the associated nodes be part of a single parallel domain, but these nodes may also be part of other parallel domains if they are not defined in an adaptive mesh domain. Analyses in which a large percentage of nodes are involved in contact might not scale well if contact pairs are used, especially with kinematic enforcement of contact constraints. General contact does not limit the domain decomposition boundaries.
接触对不能跨平行域拆分,但单独的接触对不限于在同一平行域中。使用运动学接触算法的接触对要求与所涉及的曲面相关联的所有节点都在单个平行域内,并且不与任何其他平行域共享。使用惩罚接触算法的接触对要求关联节点是单个并行域的一部分,但是如果这些节点没有在自适应网格域中定义,则它们也可以是其他并行域的一部分。如果使用接触对,尤其是在接触约束的运动学强制下,大比例节点参与接触的分析可能无法很好地扩展。一般接触不限制区域分解边界。
猜测是某一个接触里包含了太多的节点。将接触的节点范围减小,重新提交检查,domain 16占比果然减小。
DOMAIN 1 HAS 5725 ELEMENTS AND 6974 NODES with WEIGHT 4.4425
DOMAIN 2 HAS 5791 ELEMENTS AND 7193 NODES with WEIGHT 4.4426
DOMAIN 3 HAS 6988 ELEMENTS AND 9097 NODES with WEIGHT 4.4422
DOMAIN 4 HAS 6980 ELEMENTS AND 8503 NODES with WEIGHT 4.4428
DOMAIN 5 HAS 6837 ELEMENTS AND 8377 NODES with WEIGHT 4.4425
DOMAIN 6 HAS 5847 ELEMENTS AND 7055 NODES with WEIGHT 4.4423
DOMAIN 7 HAS 5963 ELEMENTS AND 7232 NODES with WEIGHT 4.4422
DOMAIN 8 HAS 5865 ELEMENTS AND 7234 NODES with WEIGHT 4.4421
DOMAIN 9 HAS 7991 ELEMENTS AND 10520 NODES with WEIGHT 4.4430
DOMAIN 10 HAS 5762 ELEMENTS AND 6977 NODES with WEIGHT 4.4421
DOMAIN 11 HAS 5876 ELEMENTS AND 7189 NODES with WEIGHT 4.4419
DOMAIN 12 HAS 8101 ELEMENTS AND 10610 NODES with WEIGHT 4.4424
DOMAIN 13 HAS 5861 ELEMENTS AND 7256 NODES with WEIGHT 4.4427
DOMAIN 14 HAS 5878 ELEMENTS AND 7078 NODES with WEIGHT 4.4420
DOMAIN 15 HAS 5806 ELEMENTS AND 6909 NODES with WEIGHT 4.4423
DOMAIN 16 HAS 39537 ELEMENTS AND 55286 NODES with WEIGHT 33.364