计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究.pdf

应用数学和力学编委会，：计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究范志庚，万 强，牛红攀，靳 凡（中国工程物理研究院 总体工程研究所，四川 绵阳）摘要：基于黏弹性基本理论，引入材料非线性特征，考虑了材料加载、保载应力松弛历史、老化效应以及黏弹性模型各运动单元退化的差异性，并从两种老化机制出发，获得了老化硅泡沫垫层力学模型以及长时应力松弛硅泡沫垫层接续加载力学模型 模型机理清晰，能够反映材料服役历史信息及其对力学效应的影响 关 键 词：硅泡沫垫层；非线性黏弹性；应力松弛；老化；力学模型中图分类号：文献标志码：，（，）：，：；引 言多层预紧结构是复杂装备中常用的一类结构形式，其采用软垫层材料通过过盈预紧装配的形式将众多物理功能件进行集成 其中，软垫层材料力学性能的稳定性通常是制约多层预紧结构长时可靠性的敏感因素 先进发泡工艺制备的硅泡沫软垫层是一种常备方案 为研究多层预紧结构全生命周期力学状态的演化变化规律，众多学者围绕硅泡沫软垫层的力学性能，如短时压缩、长时应力松弛和老化等开展了大量的理论应用数学和力学 卷 期 年 月 ，收稿日期：；修订日期：基金项目：国家自然科学基金委员会中国工程物理研究院 联合基金（）作者简介：范志庚（），男，高级工程师，博士（通讯作者：）引用格式：范志庚，万强，牛红攀，靳凡 计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究 应用数学和力学，（）：和实验工作 等提出了弹性可压缩固体模型 文献从超弹性材料一般性建模方法，即将应变能函数分为等容变形部分和体积变形部分，对常见的体积应变能函数表达式进行了概述 刘占芳等考虑到单轴压缩条件下孔隙率对硅橡胶泡沫应力应变关系的影响，提出了反映体积变形与孔隙度关系的体积应变能密度函数 文献介绍了 国家实验室（）建立的基于孔隙率、超弹性、各向同性的可压缩泡沫模型（），该模型具有简单的微观力学基础，并证明了其拟合现有数据的能力 在应力松弛研究方面 文献对硅泡沫试样进行了不同温度下的短时压缩应力松弛实验、测试和实验后样品的微观形貌观察，得到了试样的载荷保持率 剩余应力与初始应力之比随时间的变化曲线，研究了温度对硅泡沫材料短时应力松弛性能的影响 开展了 和 的恒定压缩量，以及在室温和规定高温下，不同密度和厚度的硅泡沫垫层 年的应力松弛实验 史平安等对硅泡沫垫层材料的应力松弛行为进行了数值模拟，获得了反应应力、温湿度效应的硅泡沫材料黏弹性本构模型，数值模拟结果与实验数据吻合较好 老化力学性能研究方面 张凯等开展了 压缩量、高温下硅泡沫材料 的贮存实验，获得了永久变形与温度以及贮存时间的关系 文献开展了辐照和热氧老化的实验研究，建立了老化材料力学性能表征模型，研究了辐照和热氧老化条件对材料压缩力学性能的影响 等基于力、热及其耦合作用对垫层力学性能退化影响规律的认识，实现了二阶 超弹性本构模型参数与力、热、老化时间的定量关联，建立了描述长时力热应力耦合影响的时变本构模型 目前研究尚存在以下问题：高聚物类材料在一个很宽的时间范围内均呈现有黏弹性，其松弛时间是一个覆盖范围很宽的连续谱 从线性黏弹性理论的观点来看，硅泡沫垫层由众多运动单元组成，这些运动单元具有不同的松弛时间，应力松弛状态下每个运动单元的松弛程度是不一样的 考虑到硅泡沫垫层材料非线性力学特性，应力松弛状态下每个运动单元的即时刚度与其初始状态以及回复后再加载状态是不一样的 显然相关研究是缺乏的，对长时服役复杂装备环境适应性研究支撑不够 现有模型多为唯象模型，缺乏对材料服役历史信息的描述，对复杂装备全寿命、全剖面、全流程评估需求构成挑战 本文基于黏弹性基本理论，考虑材料非线性特征，从加载历史出发，依次给出了静态压缩、应力松弛、长时应力松弛后接续加载 种力学状态的理论表征形式，建立了时变非线性弹簧和任意流体组合描述的，包含初始加载、应力松弛以及不可逆变形等信息的硅泡沫垫层非线性黏弹性模型 基于广义 模型的硅泡沫力学性能表征硅泡沫是由硅橡胶经发泡制成，黏弹性是其最典型的力学特征，常采用广义 模型形式进行表征 图 广义 模型 广义 模型表现为任意多组 模型的并联形式，机械模型如图 所示 对于硅泡沫垫层，其静态力学性能表现出典型的超弹特性 因此，图 中的弹簧为非线性弹簧，且有 （），（）式中，和 分别为第 个 单元的弹簧刚度和弹簧应变，（）为第 个 单元的弹簧应力函数 不妨假设单个 单元的弹簧应力具备与整体模型相一致的函数形式，即 （）（），（）式中，（）为第 个 单元的弹簧应力函数，（）为整体模型的应力函数，为第 个 单元的弹簧相对于模型的刚度比，则可以得到广义 模型形式的硅泡沫垫层应力应变关系为应 用 数 学 和 力 学 年 第 卷 ，（）（），（，），（）其中，分别为第 个 单元的应变、弹簧应变和黏壶应变，分别为第 个 单元的应力、弹簧应力和黏壶应力，为第 个 单元的黏壶应变率，为第 个 单元的黏壶黏性系数 当 时，基于此，可以得到如下依次变化的 种力学状态描述 初始静态压缩力学状态 ，（）（）（），（）（）应力松弛状态由式（）可以得到 （，）（）（，）（）（）由于（）表现为强非线性形式，通常由上式很难直接获得理论解 不妨假设 ，（）其中，为第 个 单元的载荷保持率，为整体模型的载荷保持率 当 时，；时，当 单元退化为弹簧单元时，进而可以得到 ，（）（）（）（）长时应力松弛后接续加载力学状态 （）（），（）（）（）可见，从广义 模型出发，欲获得硅泡沫软垫层长时应力松弛后叠加静态加载状态响应的预测，需要根据静态压缩和应力松弛过程，掌握材料静态压缩模型、初始压缩量以及所有 单元的载荷保持率 其中，（）由新鲜状态硅泡沫软垫层的压缩实验获得，和 可以由应力松弛实验数据获得 将图 模型分别退化为 模型和 参数标准线性固体模型，可得到两种模型描述的长时应力松弛后接续加载力学状态 模型 （），（）（）第 期 范志庚，等：计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究可见，从 模型出发，根据材料静态压缩模型、初始压缩量以及应力松弛率 载荷保持率，即可获得硅泡沫软垫层长时应力松弛后叠加静态加载状态响应的预测 参数标准线性固体模型 ，（），（）其中，分别为弹簧单元的弹簧以及 单元的弹簧相对于模型的刚度比，为整体模型的载荷保持率 可见，从 参数标准线性固体模型出发，根据材料静态压缩模型、初始压缩量、弹簧刚度比以及应力松弛率 载荷保持率，即可获得硅泡沫软垫层长时应力松弛后叠加静态加载状态响应的预测 基于广义黏弹塑性模型的硅泡沫力学性能表征实际上，在力热辐照等多因素长期作用下，硅泡沫材料不可避免发生老化等物理现象，进而出现材料硬化和压缩永久变形等力学效应，见图、图 图 材料老化不可逆变形 图 材料老化前后的压缩力学性能曲线 图 广义黏弹塑性模型 一方面，广义 模型及其退化模型均为完全可回复模型，难以描述硅泡沫材料老化带来的不可逆变形行为 同时，硅泡沫材料老化将带来材料硬化效应 亦即模型中各运动单元的弹簧刚度随时间变化，而现有模型很少描述弹簧刚度的时变效应 为此，拟采用图 所示的广义黏弹塑性模型开展长时贮存硅泡沫材料即时压缩力学性能表征，并考虑各运动单元的弹簧刚度随时间变化 可以得到图 广义黏弹性模型形式的硅泡沫垫层应力应变关系为 ，（，）（，），（，）（）（）应 用 数 学 和 力 学 年 第 卷基于此，同样可以得到如下依次变化的 种力学状态描述 初始静态压缩力学状态 ，（）（），（）其中，分别为图 中 弹簧、黏壶的应变 应力松弛状态 ，（，）（，），（，），（）其中，为第 个 单元的载荷保持率，为整体模型的载荷保持率，为不可逆应变 长时应力松弛后接续加载力学状态 （），（，），（）（，），（）其中，（，），由老化实验获得，和 可以由应力松弛实验数据获得 图 模型 将图 模型退化为图 所示的 参数 模型，可以得到 模型描述的长时应力松弛后接续加载力学状态 初始静态压缩力学状态 ，（）（）（），（）其中，分别为图 中 弹簧、弹簧、黏壶以及 黏壶的应变 应力松弛状态 （，）（，），（，）（，），（，），（）其中，为整体模型的载荷保持率，为不可逆应变 长时应力松弛后接续加载力学状态 （），（，），（，）（，），），（）其中，（，），由老化实验获得，可以由应力松弛实验数据获得 硅泡沫垫层应用研究 硅泡沫垫层静态压缩力学性能表征静态压缩状态下，硅泡沫垫层表现出超弹性行为特征 通常采用体积可压缩 模型描述：第 期 范志庚，等：计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究 （）（），（）式中，为等容变形应变能密度，为体积变形应变能密度，为材料参数，为主伸长率 采用二阶 形式，并基于体积完全可压缩假设，可以得到硅泡沫垫层单轴应力状态下的工程应力工程应变关系：（）（）（）（）采用式（）对某硅泡沫垫层压缩实验进行曲线拟合，如图 所示 图中，实心圆点为实测数据，曲线为拟合结果 可见，采用二阶 形式的体积可压缩泡沫模型可以精准再现硅泡沫垫层压缩力学性能 图 某硅泡沫垫层压缩实验曲线拟合 老化硅泡沫垫层力学性能表征和预测如前文所述，高聚物类材料在一个很宽的时间范围内均呈现有黏弹性，其松弛时间是一个覆盖范围很宽的连续谱 硅泡沫材料不外如是 也就是说，图 模型应用过程中，采用的 单元阶数越多，预测效果越佳 根据工程使用经验，综合精度和计算规模限制，本文选用 阶模型，即图 所示的 参数模型，开展长时应力松弛后接续加载状态下硅泡沫垫层的力学行为表征，同时对比分析了与新鲜状态、老化状态的差异 图 参数黏弹塑性模型 由前面的分析可以看出，材料的硬化与其变形及老化时间相关，其压缩永久变形与初始应力及时间相关，即 （，），（，）（）文献等关于硅泡沫材料老化性能的研究均表明，硅泡沫材料的老化性能变化主要来源于基体材料自身的老化以及泡孔结构的变化 基于此，建立唯象的硅泡沫材料老化模型，即 （，）（，）（）（）（）（）实际计算过程中，（，），由老化实验获得，和可以由应力松弛实验数据获得 此处暂不对硅泡沫材料的硬化模型和压缩永久变形的具体表征过程和表征形式做深入讨论，直接假定：，（）应 用 数 学 和 力 学 年 第 卷将式（）代入式（）可以得到给定条件下的老化垫层的压缩应力应变关系：（，）（）（）（）（）（）将式（）代入式（）和（），通过计算可以得到，长时应力松弛后接续加载力学状态的函数描述：（）（）（）（）（）（）由式（）、（）和（）得到新鲜状态、老化状态以及长时应力松弛后接续加载状态的应力应变关系如图（）所示 显然，取 ，意味着式（）和（）模型重点考虑了硅泡沫材料老化过程中的泡孔结构的变化机制，忽略了基体材料自身的老化效应 为此，进一步假设 ，并令（，），其他参数保持不变，即考虑硅泡沫材料老化过程中基体材料自身的老化，忽略泡孔结构的变化机制，可以得到老化垫层的压缩应力应变关系以及长时应力松弛后接续加载力学状态的函数描述 老化垫层的压缩应力应变关系如下：（，）（）（）长时应力松弛后接续加载力学状态的函数描述如下：（）（）（）（）（）由式（）、（）和（）得到新鲜状态、老化状态以及长时应力松弛后接续加载状态的应力应变关系如图（）所示（）基于泡孔结构变化老化机制（）基于基体材料老化机制（）（）图 种状态力学响应对比 对比图 和图 所示材料老化前后的硅泡沫压缩力学性能实验曲线和理论预测曲线可知，图（）预测趋势以及压缩曲线的演化和发展规律与图 实验结果较为一致 因此，从老化硅泡沫垫层压缩曲线的形态上看，基于泡孔结构变化机制的预测模型更接近于实测曲线 小 结本文通过考虑材料加载历史、黏弹性材料各运动单元应力松弛差异性以及两种老化机制，建立了老化硅泡沫垫层以及应力松弛状态硅泡沫垫层接续加载力学响应预测模型，获得以下几点认识：）老化硅泡沫垫层压缩力学性能与其老化机制直接相关，从老化硅泡沫垫层压缩曲线的形态上看，基于泡孔结构变化机制的预测模型更接近于实测曲线）由于硅泡沫垫层压缩力学性能数据表现出强非线性特征，硅泡沫垫层应力松弛和老化状态直接影响到各运动单元的刚度和剩余变形，进而影响到其后续压缩力学响应 因此，不能简单地采用新鲜或老化硅泡沫垫层性能来表征或预测此刻的硅泡沫垫层力学性能，必须综合考虑载荷历史以及材料退化效应）经历复杂载荷历程的硅泡沫垫层力学性能的量化和表征与材料的加载历程以及材料退化效应紧密第 期 范志庚，等：计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究关联，文中对其退化模型参数的量化和表征尚缺失，后续需开展系列的实验实现模型验证和参数量化 参考文献（）：，：，：刘占芳，励凌峰，胡文军 多孔硅橡胶有限变形的弹性行为 应用力学学报，（）：（，（）：（），（）：王翕，石耀刚，李明，等 温度对硅泡沫材料短时应力松弛性能的影响 化工新型材料，（）：（，（）：（）：史平安，尹益辉，王军，等 柔性硅泡沫薄片的短时松弛行为研究 固体力学学报，（）：（，（）：（）张凯，范敬辉，吴菊英，等 硅橡胶泡沫材料的热氧老化机理研究 合成材料老化与应用，（）：（，（）：（），：，（）：，：，：，：，：，：，：，：，（）：，（）：应 用 数 学 和 力 学 年 第 卷