Experimental study and numerical modelling on creep and creep recovery characteristics of geomembrane
-
摘要: 开展了3种不同厚度的HDPE土工膜6种荷载水平下的蠕变与蠕变恢复试验,对土工膜蠕变与蠕变恢复特性及其影响因素进行深入分析。研究表明:荷载施加与卸除的瞬间,土工膜产生瞬时变形,之后变形速率随时间逐渐变小。当荷载水平较低(如小于30%)时,土工膜蠕变随时间持续增加并较快地趋于稳定状态;当荷载水平较高(如大于40%~50%)时,土工膜将进入第二阶段蠕变,变形以恒定速率持续增长,不再稳定。当荷载水平超过20%时,完全卸载后土工膜蠕变恢复时产生残余变形,且残余变形随荷载水平增大近似呈线性增加。土工膜初始模量受荷载水平与蠕变过程的双重影响,当荷载水平上升至40%~50%时,其初始模量将显著下降超100 MPa,其中在加载初期初始模量下降约40%,经过蠕变后卸载初始模量进一步下降约30%。此外,改进建立了考虑材料黏弹性的土工膜蠕变联合元件模型及蠕变恢复模型,与试验数据对比验证了所建模型能较好地反映不同荷载水平下土工膜的蠕变及蠕变恢复过程。Abstract: The creep and creep recovery of HDPE geomembrane with three different thicknesses under six different load levels are investigated. The creep and creep recovery characteristics of the geomembrane and their influencing factors are deeply analyzed. The test results show that the geomembrane deforms instantly when the load is applied and removed, and then the deformation rate gradually slows down with time. When the load level is lower than 30%, the deformation of the geomembrane tends to be stable with time. When the load level is higher than 40%~50%, the strain of the geomembrane increases at a constant rate in the second stage creep. When the applied load level exceeds 20%, the geomembrane will produce residual deformation after creep recovery and complete unloading, and the residual deformation will increase approximately linearly with the increase of load level. The initial modulus of the geomembrane is influenced by both the load level and the creep process. When the load level of the geomembrane rises to 40%~50%, its initial modulus will drop significantly by more than 100 MPa. In addition, a composite element model for the creep and creep recovery of geomembrane considering the viscoelastic properties of materials is established. The numerical simulation shows that the proposed model can well reflect the creep and creep recovery process of the geomembrane under different load levels.
-
Keywords:
- HDPE geomembrane /
- creep /
- creep recovery /
- load level /
- composite element model
-
0. 引言
随着经济的不断发展和城市人口的增长,中国城市内河的污染问题日益凸显。据统计中国有80%以上的城市内河受到了不同程度污染[1]。河湖底泥的污染物主要分为有机质、氮磷营养盐、重金属三大类。
底泥污染物含量通常受底泥物理-化学性质的影响较大,大部分的有机质与黏土矿物结合在一起,并随细颗粒含量增加而增加,且有机质含量与底泥的比表面积线性相关[2];矿物表面的吸附过程对于有机质的保存起着重要作用,不同黏土矿物对有机质的吸附机理不同[3];底泥液限、塑限和塑性指数与黏粒中有机碳含量、黏粒含量与蒙脱石含量显著相关[4];不同矿物和不同有机质对底泥的物理性状影响显著不同[5]。
底泥中的有机质、营养盐和各种重金属影响底泥的物理性质和工程性质,进而影响底泥的处理技术与效果。因此,分析河湖底泥污染物与底泥物理–化学性质的相关性,对污染底泥的处置以及资源化利用具有重要的工程意义。
已有的研究成果表明了河湖底泥污染物与底泥物理–矿物成分密切相关,但是底泥物理–矿物成分–污染性状关联性的实例研究较少。本研究针对福州市晋安区河道的5处代表性污染底泥,进行了物性指标、矿物成分与污染物含量试验,同时搜集已有的国内外不同底泥污染物含量数据,分析了底泥物理–矿物成分与底泥污染物含量的关联性,并且探讨了底泥中不同污染物含量的相关关系,为河湖底泥的污染治理和处理技术选择提供科学依据。
1. 材料与方法
选取福州晋安区水系5个代表性点位,分别用A,B,C,D,E表示,如图 1所示。该水系有两条干流,分别是凤坂河和浦东河,浦东河有3条支流,分别是福兴河、新厝河、淌洋河,浦东河干流的最下游处为一个公园内的人工湖。B位于凤坂河干流的中段,A位于浦东河干流的上游点,C位于新厝河与浦东河干流的汇集处,D位于淌洋河与浦东河干流的汇集处,E位于浦东河下游的人工湖处。
采集A,B,C,D,E共5处0~10 cm深度的表层底泥,测定有机质(OM)、总氮(TN)、总磷(TP)和重金属(Cu,Zn,Ni,Pb),测定方法见表 1。底泥的颗粒组成、黏土矿物组成和界限含水率见表 2,采用筛分法和密度计法对底泥进行颗粒分析,分别采用Casagrande法和搓条法测定液限wL和塑限wP,采用X射线衍射法测定底泥黏粒中主要矿物成分,包括伊利石(I)、高岭石(K)、绿泥石(C)、蒙脱石(S)的含量。
表 1 本研究底泥污染物测定方法Table 1. Method for determination of pollutants in sediments序号 测试项目 测试方法 试验标准 1 OM 烧失量法 ASTM D2974 2 TN 凯氏法 HJ717—2014 3 TP 钼锑抗分光光度法 HJ 632—2011 4 重金属 ICP-MS法 US EPA 3050B 表 2 本研究底泥颗粒组成、黏土矿物组成及界限含水率Table 2. Particle composition, clay mineral composition and atterberg limits of sediments(%) 底泥 颗粒组成 矿物组成 wL wP Clay Silt Sand I K C S A 42.4 47.5 10.2 19 48 24 9 79.2 35.3 B 12.3 78.8 8.9 22 59 19 0 44.0 31.5 C 12.3 82.9 4.8 33 47 20 0 38.5 23.3 D 26.0 64.0 10.0 27 53 20 0 83.6 35.0 E 28.3 62.4 9.3 24 46 30 0 111.9 44.0 表 3显示了来源于文献的具有不同颗粒级配、界限含水率、黏土矿物组成,以及不同污染物含量的河道底泥数据,结合本文的试验数据,分析底泥物理–矿物成分–污染性状的关联性。
表 3 不同文献收集的底泥数据Table 3. Database of sediment pollutants compiled from literatures序号 颗粒级配 界限
含水率黏土矿物 污染物 主要污染来源 参考文献 OM TN TP Cu Zn Ni Pb 1 √ — — √ √ √ — — — — — 魏岚等[6] 2 √ — — √ √ √ — — — — — Xia等[7] 3 — — — √ √ √ — — — — 生活污水 孙广垠等[8] 4 √ — — — √ √ — — — — 废水、肥料 余成等[9] 5 √ — — — — — √ √ √ √ 废水 El-Sayed等[10] 6 √ — — — — — √ √ √ √ 养殖场 Wang等[11] 7 — — √ √ — — — — — — — Khim[12] 8 — — √ √ — — — — — — — Andrade等[13] 9 — — — √ — — √ √ √ √ 生活污水 Nguyen等 [14] 10 — — — √ — — √ √ √ √ 生活污水 牛红义等[15] 11 — — — — √ √ √ √ √ √ 废水 严玉林[16] 12 — √ — √ — — — — — — — 徐日庆等[17] 13 — √ — √ — — — — — — — Stanchi等 [18] 14 — √ — — — — √ — — — — Phanija等 [19] 15 — √ — — — — — √ — — — 储亚等[20] 16 — √ — — — — — — — √ — Ayodele等 [21] 17 — √ — — — — √ √ — √ — 吕伟豪[22] 2. 底泥试验结果与分析
2.1 污染物与底泥颗粒级配的关系
(1)有机质与底泥颗粒级配的关系
底泥有机质与细颗粒含量的关系绘制于图 2中,可以发现底泥的细颗粒与有机质之间具有较强的相关性,有机质含量随细颗粒含量的增加而增加。底泥有机质含量随细颗粒含量的关系曲线的斜率不同,斜率越大表明底泥中的细颗粒对有机质的吸附作用越强。细颗粒具有较大的比表面积,有利于对有机质的吸附和聚集。底泥中有机质不仅与颗粒级配有关,还有底泥附近的污染源和环境有关。本研究河道底泥位于城市居民区,周围有大量排污管道将居民生活废水排入河道中,使得底泥中含有较高的有机质,文献[6,7]的样品分别取自水库底泥和海湾底泥中,周围没有人为污染源,由于水库的流动性小于海湾,使得水库底泥的有机质含量>海湾底泥的有机质含量。
(2)总氮、总磷与底泥颗粒级配的关系
底泥细颗粒含量与总氮、总磷含量的关系绘制于图 3中,可以发现同一河道底泥的总氮、总磷含量随底泥细颗粒含量的增加而增加,这与有机氮、有机磷易于吸附在细颗粒上有关。底泥周边环境,黏土矿物成分的不同造成了总氮、总磷含量与细颗粒含量的关系曲线的斜率不同。
(3)重金属含量与底泥颗粒级配的关系
底泥细颗粒含量与重金属含量的关系绘制于图 4中。由图 4可见,底泥中重金属的积累受底泥颗粒级配的影响,底泥的Cu,Zn,Ni含量随底泥细颗粒含量的增加而增加,由于粒度影响底泥的比表面积、孔隙体积以及活性组分,使得底泥细颗粒具有强吸附能力,有利于重金属元素的汇集。同时,底泥粒径越细,所含有机质也越多,对重金属的吸附络合作用也越强。本研究的底泥重金属含量较高,与沿河汽车修理厂等工厂废水的长期污染有关,且本研究底泥中的有机质含量较高,使得重金属元素大量累积。
2.2 有机质与黏土矿物的关系
为了研究黏土矿物组成与底泥有机质含量的关系,选取文献[12,13]黏粒含量在40%±1%范围内的底泥,其有机质含量与各黏土矿物组成的关系绘制于图 5中,可以发现,对于相同黏粒含量的底泥,不同黏土矿物对有机质含量有不同的影响,蒙脱石与高岭土对有机质的吸附和储存能力较强,且与有机质含量呈现明显的正相关,伊利石与有机质含量呈弱负相关。虽然伊利石的比表面积大于高岭石,但是本研究对比发现高岭石含量高的底泥中有机质含量较伊利石多,其原因可能是黏土矿物对有机质存在选择性的吸附,不同的黏土矿物保存着不同的有机组分,高岭石易于吸附有机质中的—CH2基团,而在底泥中含量较多有机质是胡敏酸,—CH2是其主要官能团,易与高岭石吸附结合。这一现象有待今后积累更多的试验数据,开展进一步的探讨。
2.3 底泥各污染物之间的关系
(1)底泥有机质与总氮、总磷的关系
底泥有机质含量与总氮、总磷含量的关系绘制于图 6中。由图 6(a)可以发现底泥总氮含量随有机质含量的增加而增加,由图 6(b)可以看出底泥中总磷含量随有机质含量增加的规律不明显,有机质与总磷含量的相关性较有机质与总氮含量的相关性弱。底泥中的氮素有95%以上存在于有机物质中,因此总氮含量与有机质含量呈显著正相关。
(2)底泥有机质与重金属含量的关系
底泥有机质与重金属含量的关系绘制于图 7中。可以发现不同底泥中的重金属含量差异巨大,某些重金属元素的含量甚至相差100倍以上,本研究底泥的重金属含量明显远大于文献[14,15]底泥,这与重金属污染源有关,本研究底泥河道沿线有不锈钢加工厂、汽修厂等众多污染源,造成底泥中重金属污染严重,文献[14,15]底泥的主要污染源为生活污水,因此文献[14,15]底泥的重金属污染程度较本研究底泥轻。底泥重金属含量随有机质含量的增加而增加,不同文献底泥关系曲线的斜率不同,重金属污染源对曲线斜率的大小影响很大。
(3)底泥总氮、总磷含量与重金属含量的关系
底泥总氮、总磷含量与重金属含量的关系绘制于图 8中。可以发现底泥的重金属含量随底泥总氮、总磷含量的增加而增加。有机质与总氮总磷的同源性,以及有机质对重金属的吸附和络合作用,使得底泥重金属含量与底泥总氮总磷含量同样具有正相关的关系。
2.4 污染物与底泥界限含水率的关系
(1)有机质与底泥界限含水率的关系
底泥有机质含量与液限、塑限和塑性指数的关系绘制于图 9中,可以发现底泥中液限、塑限及塑性指数随底泥有机质含量的增加而增加,底泥的有机质含量与液塑限及塑性指数之间具有较强的相关性,液塑限、塑性指数与有机质的关系式列于图中。有机物对液限和塑限的影响是通过改变土颗粒结合水膜的厚度来实现的,有机物具有较高的比表面积和较强的持水能力,可吸附在黏土矿物表面,形成较厚的结合水膜,从而提高底泥的液塑限。
(2)重金属含量与底泥界限含水率的关系
底泥重金属含量与液限、塑限及塑性指数的关系绘制于图 10中。可以发现文献[19~22]底泥的液塑限随着重金属含量的增加而减小。重金属离子对底泥液塑限的影响主要是引起了黏土矿物的聚集和双电层厚度的改变。本研究底泥的液塑限及塑性指数则随着重金属含量的增加而增加,这是因为本研究底泥中有机质含量较高,而文献[19~22]底泥中几乎不含有机质,重金属含量会随着有机质含量的增加而增加,且有机质对液塑限的增加作用大于重金属对液塑限的减小作用。
3. 底泥污染物及物理性质相关性分析
对底泥中的污染物与物理性质数据进行相关性分析,为了减少底泥所处环境因素对相关性分析的影响,对底泥污染物和物理性质按区域进行分析,后取其相关系数平均值。底泥污染物与物理性质的相关系数,如表 4所示。有机质、总氮及Cu,Zn,Ni重金属含量与底泥中黏粒含量和粉粒含量的相关性较高,相关系数均大于0.5以上;污染物与黏土矿物相关性不强,相关系数均小于0.5,这是由于黏土矿物在整个底泥颗粒中所占比重较小,影响力有限;有机质含量与各污染物含量之间的相关系数均较高,表明底泥污染性状与有机质含量密切相关;有机质对底泥液限、塑限、塑性指数的相关系数分别为0.915,0.916,0.797,这表明有机质对底泥的物理性质有着重要的影响,其他污染物与底泥物理性质的相关系数均较低,对底泥物理性质的影响较小。
表 4 底泥污染物及物理性质的相关系数Table 4. Correlation coefficients of sediment pollutants and physical properties污染物与底泥颗粒级配 污染物与黏土矿物 污染物与污染物 污染物与底泥物理性质 Clay OM 0.231 Illite OM -0.099 OM TN 0.809 OM wL 0.915 Silt OM 0.524 Kaolinite OM -0.185 OM TP 0.456 OM wP 0.916 Sand OM -0.485 Chlorite OM -0.194 TN TP 0.623 OM IP 0.797 Clay+Silt OM 0.717 Smectite OM 0.020 OM Cu 0.636 TN wL 0.254 Clay TN 0.686 Illite TN -0.249 OM Zn 0.794 TN wP 0.242 Clay TP 0.439 Illite TP -0.356 OM Ni 0.490 TN IP 0.161 Silt TN -0.026 Kaolinite TN -0.216 OM Pb 0.777 TP wL -0.009 Silt TP -0.097 Kaolinite TP -0.269 TN Cu 0.452 TP wP -0.005 Sand TN -0.763 Chlorite TN 0.485 TN Zn 0.603 TP IP -0.033 Sand TP 0.321 Chlorite TP 0.264 TN Ni 0.511 Cu wL -0.342 Clay+Silt TN 0.763 Smectite TN 0.197 TN Pb 0.433 Cu wP 0.300 Clay+Silt TP 0.321 Smectite TP 0.054 TP Cu 0.335 Cu IP -0.350 Clay Cu 0.355 Illite Cu -0.216 TP Zn 0.577 Zn wL -0.331 Clay Zn 0.363 Illite Zn -0.382 TP Ni 0.203 Zn wP -0.317 Clay Ni 0.335 Illite Ni -0.270 TP Pb 0.501 Zn IP 0.275 Clay Pb 0.335 Illite Pb -0.286 Cu Zn 0.655 Pb wL 0.067 Silt Cu 0.572 Kaolinite Cu -0.229 Cu Ni 0.610 Pb wP -0.365 Silt Zn 0.558 Kaolinite Zn -0.329 Cu Pb 0.551 Pb IP 0.112 Silt Ni 0.639 Kaolinite Ni -0.245 Zn Ni 0.539 Silt Pb 0.006 Kaolinite Pb 0.047 Zn Pb 0.729 Sand Cu -0.555 Chlorite Cu 0.258 Ni Pb 0.406 Sand Zn -0.593 Chlorite Zn 0.482 Sand Ni -0.641 Chlorite Ni 0.445 Sand Pb -0.292 Chlorite Pb 0.362 Clay+Silt Cu 0.554 Smectite Cu -0.137 Clay+Silt Zn 0.591 Smectite Zn -0.214 Clay+Silt Ni 0.635 Smectite Ni -0.159 Clay+Silt Pb 0.295 Smectite Pb 0.103 4. 结论
基于福州晋安东区五处河道底泥系列试验研究结果,结合搜集的独立试验数据,进行了底泥物理-矿物成分–污染性状关联性分析,得出4点结论。
(1)在细颗粒含量较高的底泥中,有利于污染物的吸附积累,随着细颗粒含量的增加,污染物含量近似呈线性增加的趋势。而在砂粒含量较高的底泥中,则缺少这种吸附能力,底泥中污染物含量低。
(2)底泥中不同黏土矿物对有机质含量有不同的影响,蒙脱石与高岭土与有机质含量呈现明显的正相关。
(3)底泥中有机质与氮磷营养盐一般具有同源性,同时有机质对重金属具有络合作用,底泥中有机质含量与总氮总磷含量,有机质含量与重金属含量,总氮总磷含量与重金属含量,均具有良好的线性关系。
(4)底泥中的液限、塑限和塑性指数随着底泥有机质含量的增加而增加,相比于有机质,重金属对底泥界限含水率的影响较小。
-
表 1 HDPE土工膜拉伸力学特性
Table 1 Properties of geomembranes
厚度/mm 屈服强度/(N·mm-1) 屈服应变/% 断裂强度/(N·mm-1) 断裂应变/% 0.5 5.71 12.49 5.04 506.38 1.0 10.35 22.36 11.45 416.15 1.5 15.75 17.00 15.68 610.32 表 2 土工膜蠕变与恢复初始模量
Table 2 Comparison of initial moduli under different load levels
膜厚/mm 较低荷载水平/MPa 较高荷载水平/MPa Es'/Es Eu'/Eu Eu/Es Eu'/Es Es Es' Eu Eu' 0.5 432.43 387.49 274.58 164.08 0.90 0.60 0.63 0.38 1.0 387.53 381.91 194.90 155.17 0.99 0.80 0.50 0.40 1.5 287.17 276.45 194.40 145.80 0.96 0.75 0.68 0.51 平均值 — — — — 0.95 0.71 0.60 0.43 表 3 1.0 mmHDPE土工膜蠕变及蠕变恢复模型参数
Table 3 Parameters of creep deformation and recovery model of 1.0 mm-thick HDPE geomembrane
荷载水平/% 蠕变 蠕变恢复 E1/MPa E2/MPa η1/(MPa·h) η2/(MPa·h) E′1/MPa E′2/MPa η′1/(MPa·h) 10 406.62 711.11 531.86 — 406.62 720.38 1852.12 20 381.84 475.39 1656.88 — 385.75 478.60 1762.96 30 374.15 244.93 784.05 — 353.37 306.75 682.31 40 205.54 262.52 611.81 53866.31 190.53 331.86 627.10 50 192.61 164.97 381.58 28712.66 131.71 249.36 209.89 60 186.56 176.94 450.14 26936.52 143.26 228.82 387.86 -
[1] MÜLLER W W. HDPE Geomembranes in Geotechnics[M]. Heidelberg: Springer Berlin Heidelberg, 2007.
[2] BHOWMIK R, SHAHU J T, DATTA M. Failure analysis of a geomembrane lined reservoir embankment[J]. Geotextiles and Geomembranes, 2018, 46(1): 52–65. doi: 10.1016/j.geotexmem.2017.10.005
[3] NISHIGATA T, NISHIDA K, NISHIKAWA T, et al. Creep characteristic and estimation of long term strength of HDPE geomembrane liner in waste landfill[J]. Doboku Gakkai Ronbunshu, 2002, 61(722): 75–83. https://www.jstage.jst.go.jp/article/jscej1984/2002/722/2002_722_75/_article
[4] KOERNER R M, HSUAN Y G, KOERNER G R, et al. Ten year creep puncture study of HDPE geomembranes protected by needle-punched nonwoven geotextiles[J]. Geotextiles and Geomembranes, 2010, 28(6): 503–513. doi: 10.1016/j.geotexmem.2009.12.014
[5] 王广月, 李华銮, 李艳琴. 复合土工膜蠕变性能的试验研究[J]. 岩土力学, 2009, 30(6): 1599–1603. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200906012.htm WANG Guang-yue, LI Hua-luan, LI Yan-qin. Experimental research on creep properties of composite geomembrane[J]. Rock and Soil Mechanics, 2009, 30(6): 1599–1603. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200906012.htm
[6] 杨武, 侍克斌, 何建新, 等. 不同膜厚复合土工膜的蠕变特性及模型研究[J]. 岩土工程学报, 2021, 43(5): 955–961. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202105025.htm YANG Wu, SHI Ke-bin, HE Jian-xin, et al. Creep characteristics and model study of composite geomembrane with different film thicknesses[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(5): 955–961. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202105025.htm
[7] LEADERMAN H. Creep and creep recovery in plasticized polyvinyl chloride[J]. Industrial & Engineering Chemistry, 1943, 35(3): 374–378. doi: 10.1007/BF01530912
[8] HAO A Y, CHEN Y Z, CHEN J Y. Creep and recovery behavior of kenaf/polypropylene nonwoven composites[J]. Journal of Applied Polymer Science, 2014, 131(17): 40726. doi: 10.1002/app.40726
[9] 敬凌霄. 多轴向经编聚酯织物增强膜材力学性能研究[D]. 上海: 东华大学, 2017. JING Ling-xiao. Mechanical Property of Flexible Membrane Materials Strengthened by Polyester Multi-Axial Warp Knitted Fabrics[D]. Shanghai: Donghua University, 2017. (in Chinese)
[10] 汪泽幸, 李帅, 谭冬宜, 等. 循环加载处理对聚氯乙烯涂层膜材料蠕变性能的影响[J]. 纺织学报, 2021, 42(7): 101–107. https://www.cnki.com.cn/Article/CJFDTOTAL-FZXB202107016.htm WANG Ze-xing, LI Shuai, TAN Dong-yi, et al. Effect of cyclic loading treatment on creep behavior of polyvinyl chloride coated membrane[J]. Journal of Textile Research, 2021, 42(7): 101–107. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-FZXB202107016.htm
[11] 颜世铛. 半晶态聚合物拉伸变形行为的微观机理[D]. 重庆: 重庆大学, 2013. YAN Shi-dang. Microscopic Mechanism of Tensile Deformation of Semicrystalline Polymer[D]. Chongqing: Chongqing University, 2013. (in Chinese)
[12] MONASSE B, QUEYROY S, LHOST O. Molecular Dynamics prediction of elastic and plastic deformation of semi-crystalline polyethylene[J]. International Journal of Material Forming, 2008, 1(1): 1111–1114. doi: 10.1007/s12289-008-0174-1
[13] LIU X L, HUANG Y J, DENG C, et al. Study on the creep behavior of polypropylene[J]. Polymer Engineering & Science, 2009, 49(7): 1375–1382.
[14] CHENG J J, POLAK M A, PENLIDIS A. An alternative approach to estimating parameters in creep models of high-density polyethylene[J]. Polymer Engineering & Science, 2011, 51(7): 1227–1235.
[15] LIU H T, POLAK M A, PENLIDIS A. A practical approach to modeling time-dependent nonlinear creep behavior of polyethylene for structural applications[J]. Polymer Engineering & Science, 2008, 48(1): 159–167. doi: 10.1002/pen.20942
[16] ZHAO X, ZHANG Q H, CHEN D J, et al. Enhanced mechanical properties of graphene-based poly(vinyl alcohol) composites[J]. Macromolecules, 2010, 43(5): 2357–2363. https://www.sciencedirect.com/science/article/pii/S0167577X13004230
[17] LI Y Q, YANG T Y, YU T, et al. Synergistic effect of hybrid carbon nantube–graphene oxide as a nanofiller in enhancing the mechanical properties of PVA composites[J]. Journal of Materials Chemistry, 2011, 21(29): 10844–10851.
[18] SAWICKI A, KAZIMIEROWICZ-FRANKOWSKA K. Creep behaviour of geosynthetics[J]. Geotextiles and Geomembranes, 1998, 16(6): 365–382. https://www.sciencedirect.com/science/article/pii/S026611449800020X