在UDEC计算中,经常有小伙伴会弹出以下的错误窗口:

udec Contact Overlap Too Great

提示“块体嵌入过大”的错误。块体嵌入的示意图如下所示:

udec Contact Overlap Too Great

遇到这种错误不要慌张,在错误提示窗口中其实已经说明了可能的最常见的原因,就是法向刚度过低。如果还不能解决问题,在帮助文档的 USER’S GUIDE 中专门有一节来说明计算过程中可能遇到的各种问题的解决办法,里面就有嵌入过大的常见可能原因及解决办法。

udec Contact Overlap Too Great

在该节 3 PROBLEM SOLVING WITH UDEC3.9 Tips and Advice 中的 11. Contact Overlap 小结中这样写道:

The error message “Contact Overlap Too Great” can occur during cycling if one block penetrates too far into another. The maximum amount of overlap allowed by the code is one-half the rounding length. If such an error occurs, it is usually necessary to restart the problem from an earlier state. However, before restarting, it is important to identify the cause of the error and correct it.
Useful information concerning the contact location(s) where overlap occurred is given preceding the contact overlap message. Also, use the PLOT overlap command to identify the blocks involved. The following possible causes for contact overlap can be identified.
Joint Normal Stiffness Too Low – If the joint normal stiffness is unrealistically low for the loads applied, the blocks will penetrate too far into each other. This cause can often be identified by (a) plotting a close-up of the affected area; or (b) printing contacts in the same area. See note 4, above, for a remedy to this problem.
Numerical Instability – Numerical instabilities, characterized by increasing amplitude of oscillations, result from timesteps that are too large. History plots that show wild fluctuations are indicative of a numerical instability. The only way to correct a numerical instability without changing other problem parameters is to reduce the timestep by using the FRACTION command. Increasing mass damping parameters can often hide instability but will not likely eliminate it. UDEC automatically determines a timestep that is stable for most cases. However, situations may arise when this timestep is too large. Several situations have been identified as causing numerical instability:

如果一个块体嵌入到另一个块体中太深,则在循环计算期间可能会弹出错误消息:

Contact Overlap Too Great

程序默认允许的最大嵌入量是 rounding 长度的一半。 如果发生此类错误,通常必须从较早的状态重新开始计算。 但是,在重新计算之前,确定错误原因并进行纠正很重要。可以采用 PLOT overlap 命令来识别所涉及的块体。下面可能是块体嵌入过大的原因:

1.节理法向刚度太低

如果节理法向刚度相对于施加的荷载是不切实际的低,则块体将插入到其他块体太深。这个原因常常通过如下识别:(a)绘制影响区域闭合图或者(b)打印在同一区域的接触面。对于修正该问题,见上述 note 4

2.数值不稳定性

数值不稳定性的特征是振荡幅度增加,这是由于时间步长太大所致。History 记录图上显示出的剧烈波动就是数值不稳定的标志。在不更改其他问题参数的情况下纠正数值不稳定性的唯一方法是使用 FRACTION 命令减少时间步长。 质量阻尼参数的增加通常可以掩盖不稳定性,但不可能消除它。 UDEC自动确定在大多数情况下稳定的时间步长。然而,当这个时间步太大时就可能出现不稳定性了。已经确定了几种导致数值不稳定的情况:

(1) use of large values of stiffness proportional damping at high frequencies (see Section 4.2 in Special Features); (2) use of high values of joint dilation; (3) use of problem geometries in which one block contacts many (more than 3) other blocks on one side; and (4) use of nonreflecting (i.e., viscous) boundaries in which the bounding material is significantly stiffer than the material in the problem domain.
If the cause of the contact overlap cannot be identified, it may be necessary to use the SET cscan command. This command causes the location of the center of rounded corners to be updated more frequently, resulting in more accurate calculation of contacts and overlap.

  1. 在高频下采用了较大的刚度阻尼成分(与质量阻尼对应的阻尼) <见理论与背景的第4.2节>;
  2. 采用了较高的节理剪胀角;
  3. 釆用的几何模型中有一个块体在一个边与很多(多于3个)其他块体接触;
  4. 采用了非反射 (即粘滞)边界,且边界材料的刚度远远高于模型内部材料的刚度。

如果错误原因还不能找到,则可能有必要采用 SET cScan 命令。此命令使圆角中心的位置更加频繁地更新,从而可以更精确地计算接触和重叠。

Another procedure for eliminating this error is to increase the overlap tolerance at the start of the problem. Use the command SET ovtol. This procedure is useful if the original rounding length was very small and/or the problem geometry involved blocks with very acute angles. However, this approach can also be very dangerous if not used thoughtfully, because erroneous results can occur if ovtol is too large.

消除此错误的另一种方法是在问题开始时增大嵌入“容许值”,即使用命令 SET ovtol 进行设置。 如果初始设置的rounding 长度非常小,和(或)问题几何形状涉及具有非常尖锐锐角的块体,该方法很有用。 但是,如果不慎重使用该方法,也会非常危险,因为如果自行设置的 ovtol 太大,可能会产生错误的计算结果。

当然,上面的方法只是一些常见原因,很可能出现以上所有方法都没用的情况,这时就需要具体分析所写的命令流及检查模型的正确性了,不能一概而论。


长风破浪会有时,直挂云帆济沧海。