COMSOL Multiphysics® 5.2a 版本中发布了一项针对结构力学接触建模的新功能，它可以协助您模拟那些相互接触后就粘在一起的物体（粘附），以及受力后分开的物体（剥离），同时还包括全内聚力的模拟。让我们一起学习如何使用 COMSOL Multiphysics 的新功能来处理上述情况。
使物体相互粘结：如何模拟粘附
当对分离的固体施加压缩力，将其紧压在一起时，边界上的机械接触会使固体产生形变，以使接触边界互相契合。而如果用拉力将域分开，接触便会随即消失。这一效应可使用
COMSOL Multiphysics 中传统的接触建模进行模拟。如果物体没有分离，而是保持粘结状态，说明它们可维持拉伸力或粘附力。
事实表明，在模拟与接触和粘附现象有关的力时，我们需要格外注意切向上的力。当物体接触时，可能出现三种“相切状态”：无摩擦滑动、有摩擦滑动和摩擦粘着。除此之外还有其他的复杂因素，在许多接触过程中，只有满足某些物理条件时，两个边界才会开始粘着。例如，某种粘附材料只有在超过特定温度的环境中进行加工处理，才能有效发挥粘附作用。不过，借助
COMSOL 5.2a 版本“结构力学模块”中全新的粘附与剥离功能，便能实现对上述所有现象的准确表征。
现在我们来探讨由胶水（真实的胶水或者具有相似功能的物质）粘结起来的两个固体零件。COMSOL Multiphysics 中新的“粘附”节点是连接两个边界的关键，您可以在“模型开发器”模型树中的“接触”节点下找到这个子节点。
想要在接触建模中使用“粘附”子节点，必须先勾选“罚函数”选项。罚函数的作用可被视为使用了一个刚性的单向弹簧对来模拟接触。当两个边界相互挤压时，它们之间会形成一个虚拟的弹性薄层。在激活“粘附”模式后，弹簧切换为双向，并具有了切向刚度。如果两个边界之间存在真正的粘附层，您可以参考真实的材料数据来确定刚度。如果和实际情况相反，您可以采用高刚度值，这样就能将两个边界虚拟焊接在一起。
针对描述两层结构之间粘附现象，您可以基于以下四个准则作为依据：
-
当超过一定的接触压力时
-
当两个边界彼此在一定距离之内时
-
从分析的最初阶段开始
-
当执行了用户自定义的布尔表达式时
发生粘着行为的圆柱体的应力和形变。
让我们首先仔细观察一下上方的动画：一个圆柱被下压到一个柔性支撑面后回复。动画中的接触是利用由接触压力触发的粘附行为来模拟的。在圆柱向上运动的过程中，支撑基础被向上拉起。由于柔性基础必须保证圆柱的曲率恒定不变，所以圆柱向上运动时下方的弯曲应力几乎保持不变。另外，您能观察到边界上的应力主要集中在由粘附转换为非粘附状态过渡的部分。在现实中，这个效应会导致粘附层发生部分剥离，我们将在下文中解释这一现象。
软件允许用户对粘附的激活准则进行自定义，这为接触建模提供了巨大的灵活性。比如说，您可以模拟需要一定的温度或时间才能发挥作用的胶水。我们将通过一个新案例来展现 COMSOL 的灵活性优势。
在下一个案例中，在支撑物上以一个恒定速度拉动滑块，与此同时，从初始室温开始对整个装配体进行加热。根据粘附的激活准则，接触边界的温度需高于 365 K。下方的一组图显示了示例的几何结构及其边界条件，并应用了用户自定义的粘附激活准则。
左图：机械边界条件（三个边界上的辊支承，滑块一侧的限定水平位移，以及滑块与基础之间的接触）。右图：热边界条件（滑块顶部与基础底部的对流热通量）。图像对两个零间的热接触进行了模拟。
输入用户自定义的粘附激活准则。
执行瞬态分析时，只要接触边界的温度低于 365 K，移动滑块就几乎不受到任何阻力。一旦温度达到 365 K，就会产生弹性粘结，因此继续移动滑块就需要施加更大的力。
图像绘制了剪切应力与界面温度的变化。
另外有一点需要注意：该仿真的模型树中包含“热接触对”节点，此节点使两个零件之间的热阻随接触压力的变化而变化。由于垂直方向的热膨胀受到阻碍，接触压力会随着时间的增加而变大。因此，仿真开始时两个零件之间的传热相对较小。
使物体分离：如何模拟剥离
到目前为止，我们已经讨论了如何将两个边界粘结在一起。此外，您还能模拟此过程的相反情况：用足够大的力将两个边界扯开。借助“粘附”子节点中的“剥离”栏，我们便能够对这类现象进行控制。
剥离设置。
您可以使用新发布的剥离功能来对两层结构之间的剥离进行模拟，还可以描述连续性材料的裂纹扩展现象。对于后者，我们必须先了解裂纹路径才能够执行分析。内聚力模型（cohesive zone model，简称 CZM）便是基于以下应用条件的材料模型：
-
粘结层的应力随着边界分离距离的增加而增大。
-
在两个边界的分离距离到达一定程度之前，粘结层表现出线弹性行为。
-
在物体的弹性变形达到峰值后，应力随着进一步变形而减小。
-
消耗完一定的能量后，两层之间的粘接将彻底断裂。
-
如果在粘结完全断裂之前卸载负载，层结构将被视为受到了损坏，其弹性刚度则随之减小。
下图显示了在纯拉伸情况下法向应力与边界分离之间的关系。所需输入数据包括最大应力
与蓝色曲线下的面积，后者可以解释为能量释放率
。红色斜线表示弹性路径，部分受损结构在卸载应力过程会遵循这一路径。
