首页/文章/ 详情

两类缠结,如何连接聚合物软化、硬化与剪切带?

35分钟前浏览1

DOI:10.1016/j.euromechsol.2023.105015

推荐理由

聚合物压缩后为何先软化、再硬化?相似的加载曲线,又为何对应不同卸载和剪切带形貌?这篇文章将两类缠结映射到本构内变量,适合想把微观解释、宏观曲线与有限元云图连起来的同学。

作者研究思路

作者简化其2022年的热耦合模型,去除温度影响,聚焦室温低应变率。永久缠结PE形成背应力,动态缠结DE的解离与再生分别参与软化和额外硬化,橡胶网络提供另一部分硬化。

贡献在于保留这套物理分工,并以PC、PS压缩及PC锻造检验简化模型、分析局部化;两类缠结概念继承前作。采用ABAQUS/Explicit与VUMAT实现。

作者并联本构结构

图 1|分子间阻力与橡胶网络并联。原图4示意。

核心公式

原文动态缠结演化式

图 2|原式(33)–(36)。ξ₁、ξ₂为内变量;κ₁、κ₂为对应阻力。解离阈值随累计塑性应变γᵖ衰减;再生项随塑性伸长λᵖ增长。

原文永久缠结背应力

图 3|原式(30)、(31)。A为PE内张量,M_b为背应力;b控制饱和量级,ς影响演化速度。

原文塑性流动率

图 4|原式(39)。τ已扣除PE背应力,τ_eq再扣除DE阻力与压力项αₚp,决定塑性流动。

核心结果图

PC与PS单调压缩比较

图 5|室温10⁻⁴s⁻¹,上PC、下PS。PC大应变硬化更明显。原图6a/7a;实验分别来自Mirkhalaf2017、Wendlandt2005,PS与标定数据同源。

PC循环压缩比较

图 6|PC、10⁻⁴s⁻¹。软化后中等应变卸载较好,峰附近及大应变仍有偏差。原图8,实验Wang2019。

再生动态缠结参数对照

图 7|PC、10⁻³s⁻¹,改变ξ₂饱和值。再生DE可缓解过度弹性恢复;这是参数分析。原图B5b。

PS曲线与剪切带演化

图 8|PS、10⁻³s⁻¹。薄带交叉、边界反射,后期展宽并派生新细带。原图10a/e/i,云图为累计γᵖ。

PC剪切带演化

图 9|PC、10⁻³s⁻¹。带更快展宽融合,后期较均匀。原图15a/i;PC与PS的试样尺寸、网格不同。

PS软化更强、硬化较弱,容易维持尖锐带;PC硬化较强,更容易带动邻近材料流动。

网格敏感性对照

图 10|PS,三网格27×20、54×40、108×80;云图ε=0.22。宏观曲线接近,细网格带更窄、细节更多;各色标不同。原图11。

PC平面应变宏观比较

图 11|PC、10⁻³s⁻¹平面应变压缩。软化、平台和硬化基本吻合,屈服前非线性仍有偏差。原图14d,实验Mirkhalaf2017。

PC锻造形状和力位移

图 12|室温PC锻造:模速0.02mm/s,位移至4.6mm后卸载。形状及力位移基本一致。原图19,实验Srivastava2010。

PC锻造局部塑性演化

图 13|原图20d/i:贯穿剪切带到后期塑性区融合。模拟采用刚模、无摩擦,忽略温升;云图不是实测场。

作者核心结论

简化模型可描述室温低速下的主要压缩与成形响应。局部化同时受软化、硬化和几何影响;峰附近与大应变卸载还需引入黏弹性。

对我们的启示

本文分析:不要仅凭单调曲线选机制。应结合循环卸载、不同应力状态与局部场,并检查网格;曲线吻合不能直接证明缠结演化或带宽预测正确。


来源:Lan等,A physically-based constitutive model for amorphous glassy polymers in large deformations,European Journal of Mechanics / A Solids 104 (2024) 105015。公式、配图与封面节选自原论文。

来源:我的博士日记
AbaqusDeform非线性UM材料控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-10-23
最近编辑:35分钟前
此生君子意逍遥
博士 签名征集中
获赞 62粉丝 133文章 160课程 0
点赞
收藏
作者推荐

精品晶体塑性文章推荐

文章题目:A benchmarking CPFEM study on dimensionality reduction: quantifying grain-scale strain deviations across crystal structures, morphologies, sizes, and textures文章 DOI:10.1016/j.ijsolstr.2026.114237推荐理由:这篇文章主要研究了一个在 DIC+CPFEM 分析中经常遇到的问题:二维晶体塑性模型到底能够在多大程度上代表真实的三维晶粒变形?目前很多研究会根据表面 EBSD 数据直接建立二维或准二维 CPFEM 模型,再与 DIC 测得的局部应变场进行比较。但是二维模型缺少厚度方向的晶粒结构以及三维邻域约束,因此即使表面的晶粒形貌和取向完全一致,其局部应变响应也可能与三维模型存在明显差异。作者针对这一问题,系统比较了 FCC、BCC 和 HCP 三种晶体结构,并进一步考虑了晶粒形貌、晶粒尺寸、初始织构以及边界条件对二维—三维应变偏差的影响。相比单纯比较应变云图,这篇文章还采用 Bland–Altman 方法对对应晶粒进行逐一统计分析,因此对于做 DIC、EBSD、CPFEM以及局部应变场分析的研究比较有参考价值。作者的理论框架:作者采用经典的有限变形晶体塑性模型,包括变形梯度乘法分解、基于 Schmid 定律的滑移系分切应力、幂律滑移模型以及晶体塑性硬化模型,并通过 UMAT 在 Abaqus/Standard 中实现。FCC 考虑 12 个 滑移系;BCC 考虑 ({110}\langle111\rangle) 和 24 个滑移系;HCP 则考虑基面、柱面以及锥面 (c+a) 滑移。作者分别建立了三维 RVE 和准二维 RVE。三维模型尺寸为 100×100×100 μm,准二维模型尺寸为 100×100×5 μm。二维模型由对应三维模型表面直接截取得到,因此二者具有相同的表面晶粒形貌。模型采用 C3D8 单元,网格尺寸约为 1 μm,三维模型约包含 100 万个单元,而准二维模型约为 5 万个单元。作者考虑了三种典型晶粒形貌:等轴晶、层片晶和柱状晶。同时设置不同晶粒尺寸和不同初始织构,从而系统分析这些微观因素对于二维模型准确性的影响。边界条件方面,三维模型采用四周期边界条件(4P),二维模型分别采用四周期边界和普通单轴拉伸边界。因此:3D-4P 与 2D-4P 的比较主要反映维度降低本身的影响;而 3D-4P 与 2D-UT 的比较则同时包含维度降低和边界条件简化的影响。这种设置可以将二维模型误差的两个主要来源进行一定程度的分离。作者没有只根据应变云图判断二维和三维模型“像不像”,而是将二维和三维模型中的对应晶粒逐个匹配,并计算每个晶粒的应变差。随后采用 Bland–Altman 分析方法,通过中位数以及 2.5% 和 97.5% 分位点定义 95% 一致性区间,从而定量评价二维模型相对于三维模型产生的系统偏差以及随机离散。作者的主要结果:首先,对于等轴晶随机取向模型,FCC、BCC 和 HCP 在二维化以后都会产生明显的晶粒尺度应变偏差。其中 FCC 的偏差最小,BCC 次之,HCP 最大。作者给出的结果分别约为:FCC:30%BCC:50%HCP:63%说明晶体塑性各向异性越强,材料对于三维约束缺失通常越敏感。另外,二维模型表现出一个比较明显的规律:低应变晶粒容易被低估,而高应变晶粒容易被高估。所以即使宏观平均应变差别不大,二维模型仍可能明显放大局部应变不均匀性。晶粒形貌同样会明显影响二维—三维偏差。对于 HCP,在周期边界条件下,等轴晶模型的偏差范围约为 1.20,而层片晶下降至 0.92,降低约 23%。这说明几何各向异性与晶体各向异性之间存在明显耦合,不能简单认为某一种晶体结构的二维模型误差始终固定不变。晶粒尺寸的影响则更加复杂。作者发现,晶粒尺寸对整体系统偏差影响不大,但会明显改变二维模型误差的离散程度,而且不同晶体结构、不同晶粒形貌得到的规律甚至可能完全相反。例如等轴晶条件下,随着晶粒尺寸减小,BCC 和 FCC 的偏差有所增加,而 HCP 的偏差反而减小;对于层片晶,BCC 和 FCC 又表现出了相反的尺寸效应。最后作者研究了初始织构。结果表明,织构既可能降低二维化误差,也可能放大误差。例如 FCC 的 Cube 织构可以明显降低二维—三维偏差,而 Copper 和 Brass 织构反而会增加偏差;HCP 的 Basal 织构则能够明显降低二维模型的应变离散。可以看到,这篇文章虽然没有提出新的晶体塑性本构模型,但是比较系统地回答了一个很实际的问题:二维 CPFEM 到底在什么情况下能够比较可靠地代表三维模型?作者的结果表明,二维模型的准确性不仅与晶体结构有关,还受到晶粒形貌、晶粒尺寸、织构以及边界条件的共同影响。尤其对于局部应变集中、疲劳裂纹萌生和损伤预测等问题,由于真正关心的是少数高应变晶粒,因此直接使用二维 CPFEM 时需要更加谨慎。来源:我的博士日记

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