《Materials & Design》：Nonlinear hyper-viscoelastic modeling of first harmonic response in CB-reinforced SSBRs using LAOS: Study of Payne effect and molecular weight distribution
编辑推荐：
专业摘要翻译：
炭黑（CB）增强溶聚丁苯橡胶（SSBR）的动态力学行为对其在高性能轮胎等应用中的性能表现起着决定性作用。填充弹性体的一个关键特征是Payne效应，即在振荡剪切过程中，储能模量随应变幅值增大而显著降低的非线性现象。该效应源于填料网络及聚合物-填料
专业摘要翻译：
炭黑（CB）增强溶聚丁苯橡胶（SSBR）的动态力学行为对其在高性能轮胎等应用中的性能表现起着决定性作用。填充弹性体的一个关键特征是Payne效应，即在振荡剪切过程中，储能模量随应变幅值增大而显著降低的非线性现象。该效应源于填料网络及聚合物-填料相互作用的破坏与重建，强烈影响材料在动态载荷下的刚度和能量耗散能力。因此，准确理解和建模Payne效应对优化材料配方和预测服役性能至关重要。与此同时，橡胶分子量及其分布显著影响粘弹性响应、链缠结密度和填料分散性，从而控制整体力学行为。精确表征这些分子参数，尤其是在结合填料诱导的非线性时，对于将微观结构特征与宏观性能联系起来至关重要。大振幅振荡剪切（LAOS）测试利用橡胶加工分析仪（RPA）等仪器，为研究CB增强SSBR在真实大变形条件下的非线性粘弹性行为提供了强有力的实验途径。根据文献调研，大多数关于橡胶LAOS分析的研究以实验为主。将实验数据与先进的数学模型相结合，可增强对数据的解读，提供更深入的见解并具备准确预测材料行为的能力。
本研究应用非线性超粘弹性本构模型来描述CB增强SSBR在LAOS条件下的动态响应。通过整合分子量分布效应和填料网络动力学，该模型旨在定量捕捉Payne效应，并解释分子架构与填料增强之间的相互作用。本工作的创新之处在于将传统上为时间相关松弛和蠕变现象开发并验证的非线性超粘弹性模型创新性地应用于大振幅条件下的动态振荡剪切行为分析。尽管此类模型已被广泛用于表征时域中的粘弹性响应，但其向频域和LAOS数据的扩展仍然有限。通过利用该模型的数学框架，本研究弥合了CB增强SSBR中时间相关本构行为与非线性动态响应之间的差距。该方法能够在统一的本构描述中同时考虑非线性弹性、粘塑性流动和填料网络动力学。此外，在该模型中纳入分子量分布效应，增强了其对复杂填充弹性体的预测准确性。然而，由于现有仪器的限制，本研究仅限于第一谐波响应的分析，无法捕捉应力随时间变化的详细非正弦周期变化。尽管存在此局限性，对第一谐波的关注仍为橡胶化合物的非线性粘弹性行为提供了有价值且有意义的见解。
论文主体解读文章：
1. 研究背景与问题提出
炭黑（CB）增强溶聚丁苯橡胶（SSBR）作为高性能轮胎等工业应用中的关键材料，其动态力学性能直接决定了制品的服役表现。填充橡胶体系中一个核心的非线性现象是Payne效应，表现为在振荡剪切过程中，储能模量随应变幅值的增大而显著下降。这一现象源于填料网络以及聚合物-填料相互作用的可逆破坏与重建，深刻影响着材料在动态载荷下的刚度与能量耗散能力。因此，准确建模Payne效应对于优化材料配方和预测服役性能具有重要的工程意义。
与此同时，橡胶的分子量及其分布（MWD）显著影响其粘弹性响应、链缠结密度以及填料分散状态，从而在宏观层面控制整体力学行为。然而，将分子结构参数与填料诱导的非线性行为耦合起来进行定量描述，一直是本领域的难点。大振幅振荡剪切（LAOS）测试利用橡胶加工分析仪（RPA）等设备，为研究材料在大变形条件下的非线性粘弹性行为提供了强有力的实验手段。文献调研表明，现有研究多集中于实验表征，而将实验数据与先进数学模型结合、以增强物理解读和预测能力的工作仍相对有限。因此，本研究旨在通过非线性超粘弹性本构建模，定量捕捉CB增强SSBR的Payne效应，并阐明分子量分布与填料网络动力学之间的耦合机制。
2. 研究内容与主要结论
研究人员采用基于Zener型框架的非线性超粘弹性模型，将平衡态超弹性分支（branch A，用于描述长期、率无关响应并包含Mullins效应）与非平衡态粘弹性分支（branch B，由超弹性元件与非线性粘壶串联构成）相结合，以统一描述材料的非线性弹性、粘塑性流动和填料网络动力学。通过运动学与应力分解，建立了总变形梯度与总柯西应力的本构关系。该模型的创新之处在于将传统上用于时域松弛和蠕变现象验证的本构框架，拓展至频域LAOS数据的分析，从而弥合了时间相关本构行为与非线性动态响应之间的理论鸿沟。研究结果表明，该模型能够定量描述CB增强SSBR在LAOS条件下的第一谐波响应，并有效捕捉Payne效应随应变幅值的变化规律。分子量分布、填料含量和体系粘度被证实是控制Payne效应的关键因素。模型预测的储能模量（G′）和损耗模量（G″）与橡胶加工分析仪实验数据高度吻合，验证了模型的有效性。填料分散分析和结合橡胶含量测试进一步为模型参数提供了微观结构层面的支持。这项研究的意义在于，为复杂填充弹性体提供了一种能够在统一框架内同时考虑非线性弹性、粘塑性流动和填料网络动力学的预测工具，有助于从分子层面优化橡胶配方设计。
3. 关键技术方法
本研究主要采用以下关键技术方法：其一，建立基于Zener型框架的非线性超粘弹性本构模型，结合平衡态超弹性分支与非平衡态粘弹性分支，并引入非线性流动法则描述粘性流动；其二，通过运动学分解（将非平衡分支的变形梯度分解为弹性与粘性部分）和应力分解（总应力为平衡与非平衡应力之和）推导本构方程；其三，利用橡胶加工分析仪（RPA）开展大振幅振荡剪切（LAOS）实验，获取储能模量与损耗模量的第一谐波响应数据；其四，采用填料分散性分析和结合橡胶含量测试，从微观结构层面验证模型参数；其五，通过参数拟合程序确定材料模型参数，实现模型预测与实验数据的定量对比。
4. 研究结果
4.1 非线性超粘弹性本构方程的建立
研究人员构建了Zener型非线性超粘弹性模型，其中平衡分支A采用超弹性本构描述长期率无关响应，并考虑Mullins应力软化效应；非平衡分支B由超弹性元件与非线性粘壶串联组成，用于描述时间相关响应。通过运动学分解，非平衡分支的变形梯度张量被分解为弹性部分和粘性部分；通过应力分解，总柯西应力表示为平衡应力与非平衡应力之和。这一框架为后续储能模量和损耗模量的计算奠定了理论基础。
4.2 储能模量与损耗模量的计算方法
基于所建立的本构方程，研究人员推导了在振荡剪切条件下储能模量（G′）和损耗模量（G″）的计算方法。通过将模型预测的应力响应与LAOS实验获取的第一谐波数据相对比，实现了模型参数的系统识别。该方法使得模型能够定量描述储能模量随应变幅值增大而下降的Payne效应特征。
4.3 材料模型参数的确定
通过拟合不同填料含量和不同分子量分布SSBR化合物的实验数据，研究人员确定了模型参数。结果表明，模型参数与填料网络强度和聚合物基体粘弹性特性具有明确的对应关系。填料含量越高，Payne效应越显著；分子量分布越宽，粘弹性损耗行为越复杂。模型参数的系统标定为预测不同配方材料的动态响应提供了可靠基础。
4.4 模型验证与微观结构支持
研究人员将模型预测结果与橡胶加工分析仪（RPA）实测的储能模量和损耗模量数据进行对比，验证了模型在不同应变幅值和不同配方条件下的适用性。同时，填料分散性分析和结合橡胶含量测试结果表明，填料分散状态和聚合物-填料相互作用强度与模型参数高度相关，从微观结构层面为模型的有效性提供了独立佐证。
5. 讨论与结论
讨论部分指出，Payne效应并非单一机制所能解释，而是填料网络破坏重建、聚合物-填料界面相互作用以及分子链缠结动力学共同作用的结果。分子量分布通过影响链缠结密度和填料分散性，间接调控了非线性响应的强度。模型将时域本构框架成功拓展至频域LAOS分析，为复杂填充弹性体的非线性粘弹性行为提供了统一的数学描述。研究结论强调，所建立的非线性超粘弹性模型能够定量捕捉CB增强SSBR在LAOS条件下的第一谐波响应，分子量分布、填料含量和体系粘度是控制Payne效应的核心因素。尽管受限于仪器能力，研究仅分析了第一谐波响应，但所得结果已为理解填充橡胶的非线性粘弹性机制和优化材料设计提供了重要的理论工具与实验依据。未来工作可进一步拓展至高阶谐波分析，以获取更完整的非线性响应信息。
