双网络水凝胶的机械性能研究
- 格式：pdf
- 大小：910.20 KB
- 文档页数：6
Open Journal of Natural Science 自然科学, 2020, 8(5), 398-403Published Online September 2020 in Hans. /journal/ojnshttps:///10.12677/ojns.2020.85049双网络水凝胶的机械性能研究王超，弓莹，于蓉蓉*，惠小健西京学院理学院，陕西西安收稿日期：2020年8月13日；录用日期：2020年8月25日；发布日期：2020年9月1日摘要本文在分子动力学和蒙特卡洛模拟的分子模拟基础上，利用PEO SN和PAA SN两种单网络水凝胶聚合成PEO-PAA DN水凝胶进行水凝胶的性能研究。
通过给出机械性能评估模型，利用单轴模拟拉伸的方法将DN 水凝胶沿x轴拉伸，绘制了DN水凝胶的应力应变曲线，发现应变大于100%时，DN水凝胶的应力开始发生突变，同时模拟得出两种SN水凝胶网络呈现协同效应。
为了进一步探究这一现象，从能量的角度分析，发现PEO-PAA DN水凝胶应力变化的主要原因在于键拉伸能和角弯曲能这两个因素。
此外，分子量较小的PEO网络在DN水凝胶中以100%应变就可以达到了完全伸展，即分子量值较小的组分对DN水凝胶的机械强度起主导作用。
因此，综上所述可以通过控制每层网络的分子量来提高双网络水凝胶的机械强度。
关键词双网络水凝胶，机械性能，分子模拟Research of Study on Mechanical Properties of Dual Network HydrogelsChao Wang, Ying Gong, Rongrong Yu*, Xiaojian HuiSchool of Science, Xijing University, Xi’an ShaanxiReceived: Aug. 13th, 2020; accepted: Aug. 25th, 2020; published: Sep. 1st, 2020AbstractIn this paper, the PEO-PAA DN hydrogels polymerized by the two single network hydrogels PEO SN and PAA SN are discussed for hydrogel performance, based on molecular dynamics and monte *通讯作者。
王超 等carlo simulated molecular simulation. First of all, the mechanical performance evaluation model is provided. According to uniaxial simulated tensile method, the stress-strain curves of DN hydrogels are drawn by stretching the DN hydrogel along the x axis. In addition, it is found that the stress of DN hydrogel begins to change when the strain is greater than 100%. In the same time, a syner-gistic effect of two SN hydrogel networks is obtained. Furthermore, to further explore this pheno-menon, it is researched that the main reasons for the change of stress in PEO-PAA DN hydrogels are bond stretching energy and angular bending energy from the point of view of energy. And the PEO network with smaller molecular weight can reach full extension at 100% strain in DN hydro-gels. In other words, the components with small molecular values dominate the mechanical strength of DN hydrogels. Thus, in summary, the mechanical strength of the double network hy-drogel can be improved by controlling the molecular weight of each layer network. KeywordsDual Network Hydrogel, Mechanical Properties, Molecular SimulationCopyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0)./licenses/by/4.0/1. 引言水凝胶是一种具有亲水性的三维网络交联结构的高分子网络体系[1]。
水凝胶性质柔软，能保持一定的形状，在一定条件下吸收或者释放大量水分，同时自身对环境友好，不会造成大气环境、水环境、土壤环境的污染，因此，近年来它作为一种环保型的功能材料，越来越多的被应用于工农生活等各个领域。
研究水凝胶这个领域是一个多学科交叉课题，将涉及力学，物理，化学，材料，生命科学等，是当今最具有挑战的高技术前沿研究领域之一。
水凝胶是水膨胀的亲水聚合物网络，含有大量的水。
然而传统的水凝胶太脆，在比较低的应力或应变下会发生断裂，在一定程度上限制了水凝胶的应用。
因此，研究者们开展了大量研究工作，开发出了由两层单网络凝胶合成的双网络水凝胶[2]，极大的提升了水凝胶的机械性能。
龚剑萍和他的同事等[3] [4]就报告了其构成的双DN 水凝胶比单个DN 凝胶的断裂能力大20倍。
双网络水凝胶优异的机械性能和韧性是用于农用薄膜的好材料，同时超强的吸水能力在土壤固水方面也有较好前景。
在单网络凝胶中，聚环氧乙烷(PEO)和聚丙烯酸(PAA)水凝胶都具有较低的毒性和超强的亲水性，正是农业中稀缺的材料，但其机械脆性阻碍了其应用。
本文利用分子动力学模拟软件LAMMPS 构建了PEO 和PAA 两种SN 水凝胶，并且利用这两种单网络水凝胶聚合了PEO-PAA DN 水凝胶，在原子级别模拟的基础上，针对这两种单网络水凝胶和双网络水凝胶(PEO-PAA DN)的机械性能做了探究，为我们进一步研究水凝胶材料应用于农业生产打下了坚实的基础。
2. 机械性能评估模型[5]为了评估DN 水凝胶的机械性能，我们对水凝胶进行拉伸变形。
在模拟中，沿x 轴均匀施加应变，利用蒙特卡洛方法使每个原子随机运动，同时根据水凝胶的形变，选择性的确定每个原子新的坐标位置。
因此我们可以通过维里公式计算应力：1,,111N N xxx ij x ij i j i r f V σ−==+=∑∑王超 等其中V ，,x ij r 和,x ij f 分别表示水凝胶体积，原子i 和j 之间的矢量位置以及在x 分量上原子i 受到的力。
考虑到水凝胶是交联的网络结构，我们在模拟拉伸过程中降低了y 轴与z 轴的维数，使得泊松比为0.5，在同一个方向上保持了恒定体积。
通过对每个原子的力学计算，进行求和计算整个模拟盒的应力，我们绘制了水凝胶在拉伸过程中的应力变化曲线如图1。
Figure 1. Stress variation curve图1. 应力变化曲线通过应力应变曲线可得，在小于100%应变时，PEO-PAA DN 水凝胶的应力近似于PEO 水凝胶与PAA 水凝胶应力之和。
因此，我们认为此时DN 凝胶中的两个单网络水凝胶是相互独立的，没有产生交联。
当应变程度大于100%时，DN 水凝胶所展现的机械性能并不是两个单网络凝胶的简单相加，而是协同效应，表明交联在此时有效。
但应变达到250%时，PEO SN 的应力也开始急剧增加，为了解释这一现象，我们将从能量学的角度研究水凝胶的变化。
通过对水凝胶的分析，给定以下力场形式[6]：total vdW Q bond angle torsion inversion E E E E E E E =+++++其中,,,,,,total vdW Q bond angle torsion inversion E E E E E E E 分别表示总能量，范德华分量，静电分量，健拉伸分量，角弯曲分量，扭转分量和反转分量。
模拟过程中，原子的运动方程采用velocity-Verlet 算法[7]进行积分，在一定的精度上获得了原子的位置及速度量，并且每步积分中只需要存储一个时刻的状态变量，保证了在模拟过程中的稳定性。
聚合物的原子电荷采用电荷平衡[8]方法计算，同时利用PPPM 方法[9]计算静电的相互作用。
使用LAMMPS 模拟[10] [11]绘制出总能量与应变的变化曲线，如下图2。
从图2中可以看出，随着应变逐渐增加，两种SN 水凝胶的势能没有发生明显变化，而DN 水凝胶在100%应变后，其伸长率逐渐增大，且分子量较小的PEO 网络比PAA 网络更早达到完全伸展。
依据橡胶弹性处理网络结构的理论[12] [13] [14]，我们认为在水凝胶变形产生的应力是由构象熵引起的，即交联网络中能量不是应变的函数。
我们对其他分量的能量分析发现，随着应变逐渐增加，PEO-PAA DN 水凝胶的角弯能和键拉伸能表现出显著增加，其他能量均无显著变化，因此可以认为这两个能量是引起DN 水凝胶应力变化的主要原因。
在LAMMPS 中均可利用胡克定律描述键角弯曲势能函数和键拉伸势能函数，绘制其角弯能和健拉伸能的变化曲线，如下图3和图4。
王超等Figure 2. Total energy variation curve图2. 总能量变化曲线Figure 3. Angular bending energy variation curve图3. 角弯曲分量变化曲线Figure 4. Tensile energy change curve图4. 键拉伸分量变化曲线王超等从图3中显示，当DN水凝胶的角弯能随着应变的增加，也表现出增加的趋势，尤其是150%应变之后，变化趋势非常明显，而两种SN水凝胶没有显著变化。
