如何在 COMSOL Multiphysics 中模拟黏附与剥离

2016年 7月 28日

自 5.2a 版本起,COMSOL Multiphysics® 引入了许多新功能来提高结构力学接触建模能力,它可以帮助您模拟那些相互接触后就黏在一起的物体(黏附),以及相互分离的物体(剥离),包括完整的内聚力模拟。这篇博客,让我们一起学习如何使用 COMSOL Multiphysics 的新功能来处理上述情况。

使物体相互黏结:如何模拟附着力

当对物理上分离的固体施加压缩力,使它们紧紧地压在一起时,边界上的机械接触就会使固体发生变形,从而使接触边界吻合。如果施加的是将域分开的拉力,就不会产生接触。这种效应可以通过 COMSOL Multiphysics 中传统的接触建模功能来模拟。如果物体黏在一起,说明它们支持拉伸力或附着力。

事实表明,在模拟与接触和黏附现象有关的力时,需要格外注意切向上的力。接触的物体可能存在三种“相切状态”:无摩擦滑动、有摩擦滑动和摩擦黏接。除此之外,另外一个复杂的问题是,在一些接触过程中,只有某些物理条件得到满足时,两个边界才会开始黏附。例如,黏附材料只有在超过特定温度时才会固化,进而有效黏附。不过,借助 COMSOL 5.2a 版本结构力学模块中全新的黏附与剥离功能,就能准确表征上述所有现象。

我们以用胶水(真实的胶水或者具有相似功能的物体)黏接起来的两个固体零件为例来说明。COMSOL Multiphysics 中新的 黏附 节点是连接两个边界的关键,您可以在模型开发器中的模型树接触 节点下找到这个子节点。

要在接触建模中使用 黏附 子节点,必须先勾选 罚函数 选项。罚函数的作用可以看成使用一个刚性的单向弹簧对来模拟接触。当两个边界相互挤压时,它们之间会形成一个虚拟的弹性薄层。激活“黏附”模式后,弹簧被切换为双向,并具有切向刚度。如果两个边界之间存在真实的黏附层,可以参考真实的材料数据来确定刚度。否则就采用一个高刚度值将两个边界虚拟焊接在一起。

当涉及描述两个层之间的黏附现象时,您可以在以下4个准则中进行选择:

  • 超过一定的接触压力时
  • 两个边界彼此在一定距离之内时
  • 从分析的最初阶段开始
  • 满足用户自定义的布尔表达式时

 

发生黏附行为时,圆柱体的应力和变形。

首先,让我们仔细观察一下上方的动画:一个圆柱被压向一个柔性支撑后再向上回复。动画中的接触是通过由接触压力触发的黏附行为来模拟的。在圆柱体向上回复运动的过程中,支撑基础被向上拉起。由于柔性基础必须保证圆柱的曲率恒定不变,所以圆柱下方的弯曲应力在向上运动时几乎保持不变。另外,还可以观察到,边界上的应力主要集中在由黏附转向非黏附状态过渡的部分。在真实应用中,此效应会导致黏附层发生部分剥离,这一点我们将在下文中解释。

软件允许用户自定义黏附准则,这为接触建模提供了非常大的灵活性。比如说,您可以模拟需要特定的温度或时间才能固化的胶水。让我们来看另一个案例。

在下列案例中,以恒定的速度沿着一个较大的基础拉动滑块,同时从室温开始对整个装配体进行加热。选择黏附开始的激活准则,使接触边界的温度需高于 365 K。下面的一组图中突出显示了示例中的几何结构及其边界条件,并应用了用户自定义的黏附激活准则。

突出显示了带机械边界条件的滑块几何结构。
图像显示了滑块模型的热边界条件。

左图:机械边界条件(三个边界上的辊支承,滑块一侧的限定水平位移,以及滑块与基础之间的接触)。右图:热边界条件(滑块顶部与基础底部的对流热通量)。绘图模拟了两个零件之间的热接触。

屏幕截图展示了如何在 COMSOL Multiphysics 中输入用户自定义的黏附准则。
输入一个用户自定义的黏附激活准则。

当进行瞬态分析时,只要接触边界的温度低于 365 K,移动滑块几乎不受到任何阻力。一旦温度达到 365 K,就会产生弹性黏接,因此继续移动滑块就需要施加更大的力。

 

使用彩色表示域温度的瞬态分析。绿色等值线表示温度为 365 K 的黏附温度。箭头表示移动滑块所需的力。

绘图展示了移动滑块所需的力与边界温度的关系。
剪切应力与界面温度绘图。

还有一点需要注意,这个模拟通过两个零件之间的热阻与接触压力的函数,来包括模型树中的热接触对节点。由于垂直方向的热膨胀受到阻碍,接触压力会随着时间增加。因此,仿真开始时两个零件之间的传热相对较小。

使物体分离:如何模拟剥离

到目前为止,我们已经讨论了如何将两个边界黏接在一起。此外,您还能模拟此过程的相反情况:用足够大的力将两个边界拉开。我们可以通过 黏附 子节点中的 剥离 选项,对这类现象进行控制。

屏幕截图展示了 COMSOL Multiphysics 中的剥离设置。
剥离设置。

您可以使用新的剥离功能模拟两层结构之间的剥离,还可以描述连续性材料的裂纹扩展现象。对于后者,我们必须先了解裂纹路径才能够进行分析。内聚力模型(CZM)就是基于以下应用条件建立的材料模型:

  • 黏接层的应力随着边界分离距离的增加而增大。
  • 在两个边界的分离距离到达一定程度之前,黏接层表现出线弹性行为。
  • 在物体的弹性变形达到峰值后,应力随着进一步变形而减小。
  • 消耗完一定的能量后,两层之间的黏接将彻底断裂。
  • 如果在黏接完全断裂之前卸载载荷,层结构将被视为受到损坏,其弹性刚度将随之减小。

下图显示了在纯拉伸情况下法向应力与边界分离之间的关系。需要输入的数据包括最大应力 \sigma_{\textrm{ic}} 和蓝色曲线下方的面积,后者可以解释为能量释放率 G_{\textrm{iC}}。红色斜线表示弹性路径,出现在部分受损结构卸载应力后。

图表展示了线性分离定律中应力与边界分离的关系。
线性分离定律中应力与边界分离的关系。

剪切应力与剪切变形的关系曲线与上图类似。因此,对于纯剪切和纯拉伸情况,剥离行为被单独定义。对于更为普遍的情况,则采用混合模式的分离定律。从本质上来讲,该定律提供了两类基本情况的权重,类似情况还包括如何将有效应力用于塑性的多轴应力状态描述。

针对不同的本构定律,您可以从三个选项中做出相应的选择。除了上图绘制的线性分离定律之外,也可采用多项式定律与多轴分离定律。

两张图像分别展示了多项式定律与多线性分离定律。
分别显示了多项式定律与多线性分离定律。

为了演示如何模拟剥离行为,我们使用COMSOL 案例库结构力学模块分支下的层压复合材料的混合模式剥离教程模型为例来说明。这个案例模拟了混合模式弯曲( MMB)测试 的实验装置,目的是测试研究层压复合梁的剥离行为。在该几何结构中,我们能够创建定义清晰的拉伸与剪切混合状态。

用于 MMB 测试的试样的几何结构。
用于 MMB 测试的样本的几何结构。

如下图所示,在测试过程中,当载荷超过临界值,初始裂纹开始扩展。结构随着裂纹的增长变得更加柔韧,从而使力减小。

图表描绘了梁外缘上的力与位移的关系。
梁外缘上的力与位移的比较。

在这里,必须在位移条件受到控制的前提下进行物理实验与仿真,否则一旦载荷达到峰值,裂纹增长就会变得不稳定。

 

上层结构内的有效应力与施加于梁外缘的力。随着载荷的减小,剥离行为清晰可见。变形具有实际的尺度。

在这个示例中,我们从开始仿真时便施加了黏附作用带来的影响。需要注意的是,仿真虽然不包含实际的黏接过程,但可以将两种效应结合在一起。

在本文的第一个示例中,黏附层的一端明显承受了极大的压力。如果我添加一条剥离规则,情况会怎样呢?下图显示了由此得到的最终仿真结果。分离行为导致层结构的黏附程度变小。在转折点上以及在回复运动过程中,两个研究的结果一致。不过,之后黏附层的压力超过了剥离极限并导致接触消失。

仅施加黏附的仿真结果,以及同时施加黏附和剥离的仿真结果。
仅施加黏附作用后的最终仿真状态(上图),以及黏附与剥离作用并存时的最终仿真状态(下图)。

其他的接触建模场景:始终保持无摩擦接触

想象这样一个场景:结构中发生了黏附行为,但是黏附层的剪切刚度被设为 0。这表示两个边界在法向上相互连接,但能在切向方向上自由滑动。例如模拟轴承一类的装置便是利用了这一特性。

下面的动画展示了两个半柱面体的连接方式。较大的柱面体包含了一个阻碍水平位移的弹性弹簧,即具有辊支承条件,故无法进行垂直和旋转运动。较小的柱面体被指定进行平移与旋转。在仿真开始时,两个物体之间的力主要表现为张力。仿真快结束时,较小的物体推动较大的物体。在整个仿真中,两个连接边界上的切向力为 0。

 

半柱面体仅在法向上相互连接。

使用COMSOL Multiphysics® 提升结构力学接触建模能力

COMSOL Multiphysics 5.2a 版本中引入的全新黏附与剥离建模功能,为创建高保真的结构力学仿真提供了许多新的可能性。这些方法尤其适用于分析包含相互黏接零件的制造工艺。模拟剥离对于研究结构的最大载荷承受力具有重要意义。借助全新的接触建模功能,可以快速获取精确的仿真结果,进而可帮助我们为一系列工业领域开发出更加高效、可靠的制造工艺流程。

下一步


评论 (22)

正在加载...
yue jiang
yue jiang
2020-08-09

请问对一个局部层合板特征频率进行仿真,如何设置其接触呢

陈勇勇 陈
陈勇勇 陈
2021-07-20

特征频率分析是属于线性范畴,不能考虑接触等非线性特征;即使定义考虑了,软件可能也仅考虑初始状态作为整个结构连接状态。个人理解哈

yue jiang
yue jiang
2020-08-09

How to set up the contact of a local laminated plate by simulating its characteristic frequency?

hao huang
hao huang
2021-08-11 COMSOL 员工

Hi yue jiang

Setup the contact as normal: assemble the parts ,define contact pair in the definition , and setup contact conditions in the physical interface.

Hao Tang
Hao Tang
2021-10-31

您好,请问一下可以模拟液体在金属表面上的粘附力吗

越 赵
越 赵
2021-11-24 COMSOL 员工

您好,您所述的液体在金属表面的粘附力,其实属于双相流中润湿壁中的接触角的设置,与本文涉及的粘附不是一个概念。

启明 郭
启明 郭
2023-09-04

您好,请问您解决液滴和固体之间黏附力的仿真了吗

阳 欧
阳 欧
2022-02-28

请问,可以用于黏弹性薄膜和弹簧的接触吗?

hao huang
hao huang
2022-03-01 COMSOL 员工

可以,您可以尝试用固体力学接口建模

藩 徐
藩 徐
2022-03-25

请问,如何保证被剥离的两个界面间是按要求粘合的

Kaixi Tang
Kaixi Tang
2022-03-28 COMSOL 员工

您好,力学中是有多种模型去定义粘附和剥离的。本例中以及固体力学下 ’脱层‘ 是采用内聚力模型(Cohesive Zone Model)去定义何时会发生剥离。

瑜 鲸
瑜 鲸
2022-04-04

你好请问锚杆与围岩的接触,需要怎么设置呢

连凯 许
连凯 许
2022-07-03

请问这位朋友,您的问题解决了吗,我近期也想探究下能否用comsol模拟锚杆围岩的传力机理。

连凯 许
连凯 许
2022-07-03

您好,请问接触能否用于模拟锚杆与灌浆体之间的界面粘结呢?

Kaixi Tang
Kaixi Tang
2022-07-08 COMSOL 员工

你好,这种界面的粘接可以考虑使用博客中提到的黏附与剥离,也可以对粘接层进行损伤力学的仿真分析。

鑫 邱
鑫 邱
2023-06-20

但是我在模拟锚杆的拉拔过程中,对于“剥离”下设置的剪切强度是6MPa,为什么出来的应力应变曲线,峰值应力大约在12MPa左右。同样的,我设置剪切强度是3MPa,最后的峰值应力基本在6MPa左右

鑫 邱
鑫 邱
2023-06-20

鑫 邱2023-06-20
我在模拟锚杆的拉拔过程中,对于“剥离”下设置的剪切强度是6MPa,为什么出来的应力应变曲线,峰值应力大约在12MPa左右。同样的,我设置剪切强度是3MPa,最后的峰值应力基本在6MPa左右,这个2倍的关系时怎么出现的呢。在这个过程中我认为是只有剪切位移的,拉伸的比重是0

Kaixi Tang
Kaixi Tang
2023-06-25 COMSOL 员工

你好,使用CZM模型对黏附和剥离去进行仿真时,峰值应力的大小主要是和“牵引分离定律”、“混合模式准则”以及“混合模式度指数”的选择和设置有关的。所以说你提到的2倍关系主要是和模型中这些方面设置有关的,当改变相关设置后也不一定会一直存在2倍的关系。

HC L
HC L
2023-12-26

您好!
”牵引分离定律“、”混合模式准则”、“混合模式度指数”,这些部分涉及的参数都是什么意义呢?特别是“混合模式度指数”alpha代表什么?如何设置可以模拟仅存在法向拉伸破坏的情况?

Kaixi Tang
Kaixi Tang
2023-12-28 COMSOL 员工

您好,这些参数对应的详细信息以及物理含义,您可以在COMSOL文件夹下的帮助文档中查看:C:\Program Files\COMSOL\COMSOL62\Multiphysics\doc\pdf\Structural_Mechanics_Module

温 锦锋
温 锦锋
2024-02-01

请问方便发一下这篇博客里面“移动滑块传热”的源文件吗?

温 锦锋
温 锦锋
2024-02-01

请问方便分享这篇博客里面“移动滑块传热黏附”案例的源文件吗?

浏览 COMSOL 博客