首页/文章/ 详情

开裂还是脱粘?热障涂层失效路径如何转变

1天前浏览5

DOI:10.1016/j.compstruct.2026.120428

推荐理由

涂层加厚、基体更易塑性变形,是否就能提高热障涂层可靠性?Bie等的研究表明,这些改变可能把损伤转移到另一条路径。本文适合关注涂层失效、有限变形本构和相场实现的读者,尤其值得看横向开裂、界面脱粘与局部剥落之间的竞争。

作者研究思路

核心贡献是把既有PF-CZM(相场正则化内聚模型)拓展到有限变形弹黏塑性,使体内与界面裂纹在同一框架中竞争。裂纹场d描述损伤,界面场η平滑界面附近的强度和韧性;η并不是另一套损伤。损伤驱动力纳入塑性耗散,界面韧性则按耗散一致性修正。

UMAT和UEL的耦合流程

图 4:先用UEL求η,再由UMAT与裂纹UEL交错传递驱动力和损伤退化。

核心公式

原文式24蠕变流动

式(24):τ为等效剪应力,S为流动阻力,θ为给定温度。热激活幂律描述蠕变,不代表本文已求解温度场。

原文式20损伤驱动力

式(20):驱动力结合最大主应力项与累积塑性功,使流动历史参与损伤演化。

原文式38、39、41界面正则化

式(38)(39)(41):η连接界面与体材料的断裂能、强度;带波浪号的界面断裂能需经能量一致性修正,不能直接拿输入值代替。

核心结果图

压缩蠕变标定前半时段
压缩蠕变标定后半时段

图 9按0—35 h、35—70 h拆分:50、70、100 MPa压缩蠕变数据用于标定,并非独立验证。模型反映应力依赖趋势,但100 MPa后期应变仍偏低。

薄涂层实验
厚涂层实验

图 11(a)(b):前作实验中,100 μm涂层以横裂为主,500 μm涂层出现界面扩展和局部剥落。

三点弯曲模型

图 12:以1000℃三点弯曲模拟对照,基体厚1.2 mm,跨度16 mm。采用平面应变和半模型,基体不计损伤。

100微米涂层损伤场
500微米涂层损伤场

图 13(a)(c):同为400 μm位移,较厚涂层的界面损伤更明显。色条为损伤d,不能当成塑性应变或温度;这也是特定边界下的厚度趋势。

厚涂层损伤演化

图 15(b)(f):500 μm涂层,位移60→400 μm。横裂到达界面后,进一步发展为脱粘和局部剥落。

基体硬化对局部剥落的影响
基体硬化对界面扩展的影响

图 16(a—d):硬化模量H从0增至0.9E,局部剥落减轻,但部分界面裂纹走得更远。因此,“塑性抑制脱粘”和“低硬化加剧剥落”可以同时成立。

低界面韧性损伤场

图 20(a):左、右半域依次上下排布,便于看清局部裂纹。

高界面韧性损伤场

图 20(d):同样按左、右半域上下排布,保留原参数及指向。

不同界面韧性的反力曲线

图 20(a)(d)、21(a):涂层断裂能固定20 N/m,界面断裂能由10升至100 N/m,脱粘受抑、横裂增多,而反力曲线基本接近。承载相似不意味着局部失效相同。

作者核心结论

涂层厚度、基体塑性和界面/涂层韧性比共同改变失效路径。增强界面可以把脱粘转向涂层开裂,但并未消除损伤。现有证据限于二维、均匀各向同性近似;孔隙、三维、高温拉伸及热化学耦合仍待补充。

对我们的启示

涂层设计不能只比较峰值反力,还应同时检查界面扩展、横裂与剥落。建模时先用蠕变数据约束流动,再用裂纹基准和局部形貌检验损伤;长度尺度、网格及不可逆实现也需单独核查。


来源:Bie等,Composite Structures 389(2026),120428。原式及原图节选;图9实验来自Echsler等(2004),图11来自Long等(2019)。

来源:我的博士日记
断裂化学UM裂纹材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-10-28
最近编辑:1天前
此生君子意逍遥
博士 签名征集中
获赞 62粉丝 134文章 164课程 0
点赞
收藏
作者推荐

PMMA的屈服后软化,怎样同时考虑温度、速率与自热?

DOI:10.1016/j.polymer.2022.124524推荐理由PMMA的屈服峰可以拟合,峰后的下降幅度却随温度和速率明显变化。这篇文章值得借鉴的,是把“屈服—软化—硬化—塑性发热”分别处理,再放进同一有限变形框架。作者研究思路Lan等沿用Argon双扭折和全网络模型:温度、速率相关模量控制屈服,链取向产生硬化。主要新增贡献是将固定软化比例改为温度、速率的经验函数,同时引入温度相关塑性功转热比例β,并修正硬化参数的温度关系。 图 1|原图3:作者假想的分子链演化,非直接微观观测。参数主要由文献中的6组压缩曲线标定,再以其他温度、速率和加载方式检验;作者用ABAQUS/Explicit的VUMAT实现,未开展新的实验。核心公式 图 2|原式1、3:θ为温度,γ̇ᵖ为塑性剪切率。把E写成温度、速率的函数,使屈服阻力随工况变化。 图 3|原式8、11—13:s趋向饱和值sₛₛ,h控制演化快慢;A=sₛₛ/s₀越小,软化幅度越大。新A同时纳入加载速率和温度。 图 4|原式44—45:ρ、c为密度和比热,β是塑性功转热比例。β随温度线性近似,取值限定在0—1;它影响自热引起的额外软化。核心结果图 图 5|原图13a:298K、0.0003—0.01s⁻¹;实验Ames2009。0.01s⁻¹时,考虑自热比等温假设更接近实验。 图 6|原图13f:383K、相同三速率;实验Ames2009。最低速时软化很弱,温度效应不能用固定A表示。 图 7|原图14a:0.1s⁻¹,298/323/343K;实验Ames2009。 图 8|原图14b:0.1s⁻¹,363/373/383K;同源实验。298/383K用于热参数标定,其余温度是额外比较。 图 9|原图15a:298K、1350s⁻¹;实验Dar与Zhang2015。 图 10|原图15d:298K、4300s⁻¹;实验Richeton2006。两图的屈服后预测均偏高,说明高速软化仍有不足;4300s⁻¹数据也参与标定。 图 11|原图16:1500s⁻¹、378—413K;实验Garcia-Gonzalez2020。近转变温区的屈服和软化总体吻合。 图 12|原图17:0.1s⁻¹温升模拟;298K实验点来自Arruda1995。高温下温升较小,但多温度温升仍待直接实验验证。 图 13|原图18a:550s⁻¹/373K与2000s⁻¹/423K剪切;实验Fleck1990。 图 14|原图18b:0.0059s⁻¹、353/363/378K拉伸;实验VanLoock与Fleck2018。压缩标定后可描述这些载荷,但屈服前非线性不足,拉伸比较主要限于真应变0.5以内。作者核心结论率温相关软化函数扩大了PMMA的描述范围。β未考虑速率、应变依赖,硬化参数的线性温度近似仍粗略。对我们的启示本文分析:宏观吻合不能直接证明微观链断裂。先分清结构软化与热软化,再分层标定、跨工况检验。高速偏差应优先结合温升数据检查β,而非只继续调软化参数。来源:Lan, Jiang, Wu. A constitutive model for poly-methyl-methacrylate over a wide range of strain rates and temperatures. Polymer 240 (2022), 124524。图及公式节选自原文。来源:我的博士日记

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈