聚合算法——内聚力测试

前处理不对，得不到合理的内聚力曲线。（卡在这3天了）

合理的曲线应是下图（图1）所示，没有数值震荡，且随着开裂的增加而减小。

1.1.1 参数设置：如果得到震荡的曲线，（图2）。可能是刚度太小，单元变形不可忽略。可以调整刚度或改变外荷载作用位置。

1.1.2 建模。

如果只插入一个聚合单元，外荷载、约束如下图（图3）所示。理论上，蓝色部分应向上竖直运动，红色部分保持不动。聚合单元上下两个平面保持平行。如果出现下图（图4）所示，上一层节点对不平行的位移，则是单元不对称导致。上一层的节点，分别共同属于1,2,3个单元。

这种不对称，即使是施加同样的外荷载（速度荷载），也会导致节点位移不一致，从而发生断裂的时间不一致，位移也不一致，测试的峰值也不一致。会出现如下图的曲线，如果强行使面上所有节点对一致，再进行cohesive 计算，峰值也对不上，因为同时失效的时刻，三对的法向、切向应力，p2>p3>p1，大概是3:2:1的样子。开裂位移也不一致。

所以建模时，应使用等边三角形单元，且保证施力面上各点是对称的（至少划分情况一致）。