MARC接触问题及其解决方法.pdfVIP

接触问题及其解决方法 综述 • 本技术手册提供了在Marc和Mentat中使用接触经常遇到的问题的解决方法。 • MARC用直接约束法处理接触问题，在整个过程中跟踪运动接触体的轨迹，当发生接触时， 使用边界条件直接约束运动体——两者的运动约束转化成了节点自由度的约束和节点 力的约束。如果程序能够预测出什么时候发生约束，这种方法非常准确。在Marc中通过 CONTACT选项执行这个程序。因为不需要预先知道那里发生接触，所以在本程序中不用 提前指定参考单元，可以模拟复杂的接触。 • Marc中接触的详细用法见《程序输入手册》“接触”一章。 • Marc中接触的详细求解过程见《理论手册》（卷A）中的附录。 • 下面对接触的用法分别做一下解释： 接触分析的设置 • 在接触分析中最常见的问题怎样解决？下面的介绍可能对你有所帮助： a. 增加LOADCASE中的载荷步数。 b. 确保在JOBS的INITIAL LOADS中包括固定边界条件。 c. 接触分析允许两个或多个接触体之间发生任意大的相对变形。在一些分析中的 变形不可避免地会超出几何线性分析中小变形的设定范围（最大旋转角度的正 切值约定于角度本身，使用弧度值）。如果是这种情况，激活LARGE DISP参数。 d. 线-线之间的接触会出现数值问题，导致收敛困难。如果他们落入接触区域的容 限之内，Marc的接触算法将会设置接触的节点。容限缺省值设置为定义的所有 接触体的最小单元边长的1/20。在有些情况下，这个值太大了，尤其是线-线接 触。线-线网格不是拓扑相等的，可能一个或另外一个网格上的多数节点落入接 触区域容限内，人为设置接触应力。在接触区域容限设置时推荐使用偏离系数。 而且下降的接触区域容限推荐使用系数5或10。接触体使用解析选项，接触的线 或面使用nurbs曲线或coons曲面，很可能消除这个问题。一种特殊情况是接触 表面上的网格是连续的。即一个网格上的节点和其他网格上的节点一致。在这 种情况下，因为所有的节点都是一致的，不应该设置人为应力。但是，实践证 明如果接触体定义为离散的而不是解析的可能生成人为应力。 e. 通过运行一个虚构的一个增量步（没载荷）的载荷工况，看在只有接触的条件 下是否有应力产生，调试接触模型。这将提供任何几何问题或解析接触定义（如 果使用）的信息。 f. 在有边界条件的单元节点，当发生接触时，边界条件自动消失，代之接触条件。 如果单元后来不再接触，Marc不会恢复消失了的边界条件，这可能会导致刚体 模态。 g. 确保在RUN窗口下的“Extended Precision Input File”选项被激活——这有 很多不同：如果单元的坐标系已经输入单位系统，这些单位导致很多明显数值， 趋向于产生机器四舍五入误差。改变单位使尺寸更“大”，即使用mm而不是m。 h. 将“iterative increment splitting”（第一个接触片的域7）设置为3（在JOBS “iterative penetration”检查）。这是2001版本中的新方法，在很多情况下， 优于先前的增量分离程序。在2003版本中这很可能变成缺省值。推荐使用这个 设置。在旧的增量分离程序中当接触即将发生时，增量被分离，不幸地是这会 导致输出文件中“多重”变量的混乱。使用新的功能在同一个增量中考虑了每 个方面。 注意在2001版本中使用自动增量加载（弧长）时，因为存在一个缺陷，不能使 用迭代穿透探测。 i. 缩放比例：因为当使用米做单位时，一些穿透的数值量级非常小。为了避免圆 整误差也许必须将模型量纲（和模型所有相关的数量，如杨氏模量等）由米改