联系人简介: 郭宝华, 男, 教授, 主要从事高分子材料方面的研究. E-mail: bhguo@mail.tsinghua.edu.cn
以尼龙材料的应力松弛行为作为研究对象, 考察初始应变为1.0%, 2.8%和5.1%的尼龙1010样品在温度区间293353 K的松弛曲线, 采用时间-温度等效叠加方法得到了松弛模量主曲线, 计算出叠加过程中的表观活化能、 松弛过程中的活化体积和应力辅助功. 结果表明, 整个松弛过程中的表观活化能和应力辅助功表现出相同的变化趋势, 体现出松弛过程中克服运动单元位垒的过程. 当293323 K区间的松弛曲线叠加时, 随着初始应变的增加, 表观活化能和应力辅助功均逐渐降低, 有助于聚合物内部的运动单元越过能垒发生松弛, 与松弛过程中的应力辅助热活化理论相一致; 当333353 K区间的松弛曲线叠加时, 不同初始应变样品的表观活化能均为260 kJ/mol, 应力辅助功均为60 MPa·nm3, 说明松弛过程中克服运动单元的能垒与应力作用无关. 根据松弛主曲线, 计算出了尼龙1010在1.0%, 2.8%和5.1% 3种形变下, 长时间范围内应力衰减与时间的关系, 为预测实际使用过程中的应力松弛行为提供了依据.
Stress relaxation behavior of plastic has great influence on its service life and safety. The master curves were generated by horizontal and vertical shifts of the short term curves of stress relaxation for Nylon 1010 at temperatures ranging from 293 K to 353 K using proposed time-temperature superposition equivalence method. The apparent activation energy was calculated from the Arrhenius plots of shift factor, and the values of activation volume and stress-assisted work were calculated from the master curves. The time-temperature superposition principle was found to be available for the stress relaxation curves. During the superposition interval of the experimental temperature from 293 K to 323 K, the apparent activation energy and stress-assisted work show a tendency to decrease gradually, which is helpful for the relaxation of the moving units in the polymer over the energy barrier. During the superposition interval of the experimental temperature from 333 K to 353 K, the apparent activation energy of different initial strain samples is about the same value of 260 kJ/mol and the stress-aided work is basically unchanged at the value of 60 MPa·nm3, which indicates that the energy barrier of the moving unit is independent of the stress effect. According to the mater curves of Nylon 1010, the relationship between stress decay and relaxation time in a long time range can be obtained under the strain of 1.0%, 2.8% and 5.1%, which provides a reference for predicting stress relaxation behavior in practical use.
工程塑料以其质量轻、 耐腐蚀及成型加工性能好等优点作为金属材料的轻量化替代品, 广泛用于工业各个领域, 其中作为结构件和紧固件, 在长期使用中因应力衰减导致性能不能满足要求是一种典型的失效模式, 研究材料的应力衰减特性, 掌握长时间范围内的应力衰减规律是需要解决的重要理论和应用问题[1]. 目前, 对于长时间范围内应力衰减行为的研究, 主要通过对实验测试的应力松弛曲线进行数学归纳得到松弛过程的动力学方程[2, 3, 4, 5, 6, 7, 8], 或者建立模型对应力松弛曲线进行模拟[9, 16], 从而预测出应力衰减失效时间. 研究发现, 对于一些结晶性聚合物在长时间的应力松弛过程中存在着不同的松弛机理[2, 17, 18], 即单一的动力学方程或数学模型不能反映实际的时间-应力特性; 在实验室进行单个样品的应力松弛实验由于时间较短难以推导出长时间实际使用中的应力衰减规律. 根据聚合物在一定条件下温度和外力作用时间对其力学性能影响具有等效关系的特点[19], 通过系列环境温度下实验得到的应力松弛曲线, 采用时间-温度等效叠加方法, 推导出常温下长时间的应力衰减行为来预测应力失效时间、 确定使用寿命是可行的方法.
尼龙1010是重要的工程塑料[20, 35], 从20世纪70年代开始广泛应用于机械、 化工及航空航天等多个领域, 我国对该材料的合成路线[20]、 结构性能表征[21, 28]、 加工特性分析[29, 30]和使用性能评估[32, 34]等方面进行了系统研究, 但长期应力下的应力衰减行为研究还有待深入.
时间-温度等效叠加方法通过测试聚合物不同温度下的模量-时间松弛曲线, 以一个确定的温度为参考, 将不同温度下的松弛模量曲线沿对数时间坐标平行和垂直移动至叠合, 得到聚合物在参考温度下较长时间的应力衰减行为[19, 36]. 时间-温度等效叠加方法最早应用于无定型聚合物松弛行为的时间-温度转换研究[1]. 研究发现, 半晶型聚合物在满足一定条件后也可以使用该方法[37, 40]. 这是由于以下3方面原因: (1) 应力松弛的初始应变在应力-应变曲线的线性区域; (2) 松弛样品在实验之前需要进行预处理以减小松弛过程中聚集态结构的变化; (3) 松弛过程中聚合物不应发生化学结构的改变. 目前, 时间-温度等效叠加方法已经广泛应用于半晶型聚合物的应力松弛研究[41, 45].
本文针对尼龙材料单次应力松弛实验难以测试出长时间范围内应力衰减的问题, 探讨以时温等效叠加方法进行尼龙1010应力松弛行为研究的可行性以及长时间范围内的应力衰减规律, 预测材料工程应用中的使用寿命; 分析时温等效叠加计算出的表观活化能随松弛时间增加和应力衰减时的变化关系, 研究松弛过程中分子运动单元需克服能垒的变化, 进一步探索半晶型聚合物长时间范围内的应力松弛机理.
尼龙1010应力松弛初始样品为挤出棒材, 直径 30 mm, 原材料为上海市赛璐璐厂产品. 树脂的相对黏度1.61, 加工后应力松弛测试段直径 13 mm, 长 95 mm. 应力松弛采用河北承德金建仪器公司的万能材料试验机测试, 初始应变拉伸速率为5 mm/min; TMA-60型热机械分析(TMA)仪, 日本Shimadzu公司, 采用膨胀模式, 测试温度283393 K, 升温速率5 K/min, 样品直径2 mm, 长度8 mm; Q2000型示差扫描量热仪(DSC), 美国TA公司, 测试温度233493 K, 升温速度10 K/min; Rigaku R-Axis Spider X射线分析仪, 装配成像系统, 射线源为金属钼, 扫描区域为200 μ m× 200 μ m, 测试电压50 kV, 测试电流80 mA, 曝光时间600 s, 样品尺寸30 mm× 12 mm× 2.5 mm.
应力松弛曲线参考文献[34]中的数据. 为了增加不同温度下曲线的重合区域, 参考文献[34]实验方法在323 K时进行了初始应变为1.0%, 2.8%和5.1%的应力松弛实验, 得到尼龙1010在293353 K间隔为10 K的应力松弛曲线. 在此基础上, 按照式(1)将应力-时间曲线转换为模量-时间坐标(图1), 然后根据时间-温度叠加原理[19], 以293 K为参考温度, 将不同温度下的曲线沿对数时间坐标轴移动, 多条曲线叠加后得到样品在参考温度下长时间范围的模量-时间主曲线.
式中: σ 为松弛应力; ε 为应变; E为松弛模量.
曲线叠加时分为水平移动和垂直移动两步[19]. 水平移动参考式(2)(4)进行:
式中: ρ 为样品密度; t为松弛时间; T1和t1分别为参考温度和平移后的时间; T2为实验温度; α T为平移因子, 是保持曲线形状不变条件下的时间标尺位移量[19].
曲线叠加时的垂直移动参考下式进行:
式中: ρ (T1)T1/[ρ (T2)T2]可用β T表示, 为垂直修正因子, 是曲线在垂直方向的移动量, β T值与密度和温度2个参数相关. 式(6)中V0为松弛段样品初始体积; α V为体积膨胀系数, 表示随着环境温度的变化, 样品体积变化与初始体积之比; α l为线膨胀系数, 表示随着环境温度的变化时松弛段样品长度变化与初始长度之比; L0为松弛段样品初始长度. 对于膨胀系数较小的各向同性材料, 体积膨胀系数可按线膨胀系数的3倍计算[46], 实验中通过测试线膨胀系数来计算β T.
对于半晶型聚合物, 在采用时温等效原理进行松弛曲线的叠加时, 应尽量减小松弛前后样品取向程度和结晶度等聚集态结构的变化[37], 因此样品在松弛实验之前进行了退火处理, 并且实验时将初始应变控制在应力-应变曲线的线性区域[34].
图2为尼龙1010样品线膨胀性测试曲线. 可以看出, 在293 K时, 样品长度L293 K约为8.07 mm, 随着温度的升高, 曲线斜率基本不变, 在玻璃化转变温度附近未出现明显的转折, 线膨胀系数可以看作1个常数. 通过对曲线的线性拟合, 得到曲线的斜率b为0.93, 根据式(7)计算出线膨胀系数为1.16× 10-4 K-1. 将该值代入式(5)中, 能够计算出不同温度下的β T值.
图3为初始应变为5.1%时, 不同温度下尼龙1010样品应力松弛实验后, 沿拉伸方向制备的厚度为2.5 mm的薄片的二维X射线衍射图. 可以看到, 在应力松弛实验前后, 不同初始应变样品的衍射图均出现了3个衍射环, 由内到外分别对应于尼龙1010的(002), (100)和(110/010) 3个晶面方位角的衍射[23], 样品呈现典型的尼龙α 晶型[31]. 图4(A)为图3中(002)晶面衍射的Azimuthal图. 可以看出, 尼龙1010样品在松弛实验前后均未产生衍射峰, 样品的(002)取向为0. 图4(B)为图3中(100)晶面衍射的Azimuthal图. 曲线经高斯拟合后计算[47]得到的取向度列于表1. 可以看到, 松弛实验前后样品的(100)取向无明显改变(51.2%).
![]() | Fig.3 X-Ray diffraction patterns of Nylon 1010 samples (A) Initial; (B) 5.1% strain at 293 K; (C) 5.1% strain at 313 K; (D) 5.1% strain at 353 K. |
![]() | Table 1 Crystallinity(Xc) and orientation degree(Π ) of Nylon 1010 samples |
针对样品松弛前后可能产生的热物理参数变化, 对松弛前后的吸放热性能进行了DSC测试, 得到应力松弛前后样品在环境温度变化时的吸放热值, 根据式(8)对比松弛前后结晶度的差别. 其中Δ H为熔融焓, Δ
图5为初始应变为5.1%时, 尼龙1010在不同温度下进行应力松弛实验前后的DSC曲线. 由式(8)相应的熔融焓计算的结晶度见表1. 可以看出, 松弛实验前后各样品的DSC曲线形状和Δ H值均未出现明显的差别, 说明松弛实验前后样品没有明显的结晶度(38.4%)变化. 此外, 在温度为 318 K 附近的尼龙1010玻璃化转变温度区[34], 图5中曲线没有出现明显的吸热峰, 与图2中线膨胀曲线没有突然转折的现象相对应, 说明样品在松弛实验前的热处理提高了聚集态结构的稳定性.
以293 K为参考温度, 通过时间-温度叠加方法得到不同初始应变下的松弛模量主曲线, 结果如图6所示. 可以看出, 对于不同温度下的松弛曲线进行叠加时, 3条主曲线在大部分的叠加区均重合较好, 只是在松弛的初始阶段, 相邻2条曲线会出现偏差. 这种现象与Hejlova等[8]和Gupta等[48]的实验结果相同. 这是由于在松弛的初始阶段, 应变刚刚达到预定值, 样品处于由强制拉伸到突然静止的状态, 分子链在失去了强制的拉伸应力时, 产生很大的反作用力, 松弛速率快速增加至最大值[34]. 由于这一阶段的松弛并不完全是由分子自发的热运动造成的, 所以在曲线叠加时与相邻的温度曲线不重合. 参考Hejlova等[8]的工作, 曲线叠加时不考虑初始阶段不重合的部分, 同时为方便对松弛主曲线进行计算, 本文绘制的松弛模量主曲线也只考虑不同温度下松弛曲线的重合部分.
图7为平移因子(α T)与温度的关系. 可以看出, 在323333 K温度范围内, 不同初始应变曲线存在着转折区, 该温度区域正是尼龙1010的玻璃化转变区. 说明了在玻璃化转变温度前后应力松弛的机理并不相同, 依据阿累尼乌斯方程[49], 相应计算出的松弛过程中与温度相关的表观活化能也会有所差别.
![]() | Fig.7 Temperature dependence of logarithm of shift factor for the samples (A) lnα T and T; (B) lnα T and T-1. |
将图7(A)中的坐标转换为lnα T与绝对温度的倒数关系, 如图7(B)所示. 参考下式:
式中: Δ Hα 为松弛过程中的表观活化能, 表示松弛过程中的温度依赖性[49]; R 为气体常数. 将293323和333353 K的平移因子分别拟合计算出活化能, 结果见表2.
![]() | Table 2 Apparent activation energies(Δ Ha) in Arrhenius plots of Fig.6 |
从表2可以看出, 在实验温度范围为293323 K时, 随着样品松弛时初始形变的增加, 表观活化能逐渐降低. 原因是随着初始应变的增加, 样品受到的应力更大, 强迫聚合物内部的运动单元越过能垒发生松弛, 即松弛克服的能垒降低, 表观活化能呈现出逐渐降低的趋势, 与应力辅助的热活化理论相对应[34]. 当温度大于333 K时, 3个样品松弛过程中的表观活化能均为260 kJ/mol, 与样品初始应变无关. 原因是在此时较高的环境温度下, 1.0%应变的应力辅助作用已使链段可以自由运动(尼龙1010的Tg为343 K), 因此在更高的应力作用下促进松弛的效果与1.0%应变时应力促进松弛的效果基本相同, 活化单元都表现为链段运动.
研究[34]发现尼龙1010在应力松弛过程中会先后呈现2种机理: (1) 松弛的初始阶段, 应力有助于聚合物运动单元能垒的降低, 松弛行为满足应力辅助热活化理论; (2) 松弛过程中应力辅助功近似不变的阶段. 由于在单个温度下松弛实验的时间相对较短, 单次实验下松弛曲线不能很好地区分2种机理产生的条件. 将图6的松弛模量叠加主曲线去除相邻实验温度曲线中初始阶段不重合的部分, 得到3种初始应变下的松弛模量叠加主曲线, 其活化体积和应力辅助功也应该与293 K时间实验下的松弛曲线具有相似的变化规律, 因此可以利用叠加主曲线来分析随着等效松弛时间增加或应力衰减下的松弛机理.
根据叠加主曲线计算出对应的时间坐标、 松弛应力与应力辅助功的关系如图8所示. 其中应力辅助功为活化体积与松弛过程中真应力的乘积, 活化体积[50]的计算如下:
式中: σ 0和E0为初始应力和初始模量; σ * 为松弛过程中的真应力; σ ∞ 为残余应力, 对于初始应变为1.0%, 2.8%和5.1%的残余应力分别为1.8, 5.2和10.0 MPa[34]; Va为活化体积; σ * Va为应力辅助功, 指聚合物运动单元从当前位置跃迁到下一位置时应力所做的功[50]; k为波尔兹曼常数.
由图8(A)可以看出, 在松弛初始阶段, 相同时间段下随着初始应变的增加应力辅助功逐渐降低, 参考图6, 该时段对应的实验温度处于323 K 以下松弛曲线叠加段, 与表观活化能随初始应变的增加逐渐降低的趋势一致, 在此阶段初始应力越大, 样品松弛越容易, 即应力有助于聚合物内部的运动单元越过能垒发生松弛, 导致松弛活化能和松弛能垒的降低, 因此从表观活化能的角度解释了松弛过程中的应力辅助热活化理论. 随着松弛时间的增加, 对应于玻璃化转变温度以上的松弛曲线叠加段, 计算出的表观活化能近似恒定, 说明此阶段松弛过程中克服运动单元的能垒与应力作用无关, 相应图8(A)中3种初始应变的应力辅助功均为60 MPa· nm3. 图8(B)为不同初始应变下松弛应力的衰减与应力辅助功的关系. 可以看出, 不同的初始应变产生初始应力值不同, 但随着应力的衰减, 应力辅助功均逐渐减小至定值后基本保持不变; 对于不同初始应变下应力辅助功均逐渐减小, 不重合的曲线部分对应于表观活化能下降的阶段, 应力辅助功保持不变的部分对应于表观活化能近似恒定的阶段. 综合以上2种情况, 整个松弛过程中的表观活化能和应力辅助功均表现出先下降后为定值的相同变化趋势, 分别对应于应力辅助热活化松弛机理和活化单元协同运动松弛机理.
进行应力松弛研究最重要的工程应用是掌握零部件在实际使用中的应力下降情况, 预测松弛失效时间. 根据松弛模量主曲线, 以年度为横坐标, 分别以剩余应力和剩余应力与初始应力比值为纵坐标, 计算出了尼龙1010在 1.0%, 2.8%和5.1% 3种形变下实际应力的下降情况, 结果如图9所示. 由图9(D)可以看出, 随着初始形变的增加, 松弛过程中应力衰减的速度相应增加, 以应力衰减到初始应力的70%为例, 可以计算出初始应变为5.1%, 2.8%和1.0%样品的应力衰减时间分别为2 h, 15 d和16 a. 参考图9(A)(C)可以计算出, 其初始应力分别为56, 31和11 MPa, 应力下降为初始应力70%时残余的应力分别为39, 23和8 MPa.
![]() | Fig.9 Master curves of relaxation samples under the initial strains of 1.0%(A), 2.8%(B), 5.1%(C) and the comparison of different strains(D) |
针对尼龙1010应力松弛曲线特性, 由于该材料玻璃化转变温度较低, 因此常温下保持其固定形变时, 在应力辅助作用下因形变产生的应力容易使运动单元活化, 表观活化能降低, 进入应力辅助的热活化松弛阶段, 此时的松弛速率会快速增加, 从而应力松弛现象更为明显.
尼龙1010样品应力松弛前后结晶度和取向度均未出现明显的差别, 按照时间-温度叠加方法得到的松弛模量主曲线重合良好, 整个松弛过程中的表观活化能和应力辅助功表现出相同的变化趋势. 在松弛初始阶段, 相同时间范围内随着初始应变的增加表观活化能逐渐降低, 即应力有助于聚合物内部的运动单元越过能垒发生松弛, 符合应力辅助热活化理论; 在随后的松弛阶段, 不同初始应变样品的表观活化能保持定值, 说明此时松弛过程中克服运动单元的能垒与应力作用无关, 与长时间松弛后应力辅助功保持恒定的现象一致. 根据松弛模量主曲线, 计算出了尼龙1010在1.0%, 2.8%和5.1% 3种形变下, 长时间范围内应力衰减与时间的关系, 为预测尼龙1010材料实际使用过程中的应力松弛行为提供了依据.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
[38] |
|
[39] |
|
[40] |
|
[41] |
|
[42] |
|
[43] |
|
[44] |
|
[45] |
|
[46] |
|
[47] |
|
[48] |
|
[49] |
|
[50] |
|