Distribution and basic characteristics of new transparent clay
-
摘要: 人工合成透明土是研究土体内部的变形及渗流情况的良好工具,传统透明土多针对砂土,且在量化透明度方面尚未形成统一的评价系统。针对上述问题,提出一种新型透明黏土——AVC透明土的配制方法,和一种透明土的透明度量化评价方法,并用该方法对比分析AVC透明土与传统透明土的透明度;开展固结试验与微型十字板剪切试验,得到该透明土的压缩固结特性与强度影响因素。试验结果表明:AVC透明黏土具有250 mm可视厚度,透明度优于传统透明土;压缩特性与天然黏土十分相似,固结系数与淤泥相似;2.5%浓度AVC透明土比4.5%浓度LAPONITERD透明土的峰值强度更高。Abstract: The synthetic transparent soil is a good tool to study deformation and seepage inside the soil. The traditional transparent soil is mostly aimed at sand, and a unified evaluation system has been not yet formed in terms of quantitative transparency. In view of the above problems, a new type of transparent clay, AVC transparent soil formulation method, and a transparent soil transparency quantitative evaluation method are propesed to compare and analyze the transparency of AVC transparent soil and traditional transparent soil. The compression and solidification characteristics and strength factors of the transparent soil are obtained. The test results show that the AVC transparent clay has visual thickness of 250 mm, and the transparency is better than that of the traditional transparent soil. The compression characteristic is very similar to that of the natural clay, and the solidification coefficient is similar to that of the silt. The peak strength of AVC transparent soil of 2.5% is higher than that of the concentration of LAPONITERD transparent soil of 4.5%.
-
Keywords:
- transparent clay /
- preparation /
- transparency /
- strength /
- compression characteristics
-
0. 引言
膨胀土具有较强的亲水性,对于处于季节冻土区的渠道工程、边坡工程等,反复冻融循环对膨胀土的劣化作用会极大影响工程稳定性[1-2]。已有的研究表明,冻融作用对土体的性质改变相对复杂,部分特征并无统一认识[3-4],典型的如内摩擦角,变大、变小、基本不变的结论均有。由于冻融循环作用下,不同土体响应机制的复杂性,针对膨胀土特定土体的冻融分析在工程实践中具有更好的指导意义。
目前针对膨胀土的冻融循环研究已有了较大进展,蔡正银等[2]以高寒渠道膨胀土为研究对象,在干湿冻融耦合的试验条件下系统讨论了冻融作用下膨胀土的裂隙演化规律。Lu等[5]、Tang等[6]则讨论了不同状态下的膨胀土与循环冻融次数相关的劣化规律,研究表明,随着冻融循环次数的增加,膨胀土的破坏强度、弹性模量、抗剪强度等力学指标劣化加剧,且第一次冻融循环下其劣化最为明显。
目前有关冻融循环下的损伤模型,更多是基于连续介质理论展开。损伤力学中的损伤因子D部分研究需要借助其他测试手段进行确定,如卢再华等[7]通过CT扫描技术定义损伤因子,其构建的损伤模型包含23个参数,均需测试确定。目前主流的损伤因子定义是在已有的应力应变测试成果上进行,其中假设土体微元满足weibull分布,基于应变等价原理开展损伤模型构建获得了广泛应用。杨明辉等[8]、李向东等[9]进行了基质吸力和外部荷载耦合状态下的非饱和土损伤模型构建。Tang等[6]基于多种技术手段综合研究了损伤程度与冻融循环次数之间的关系,并构建了冻融循环作用下的统一损伤模型。总结目前的土的损伤模型,其损伤定义时,主要考虑了弹性模量的退化,核心均集中到应力应变曲线的模拟,但是对于土体本身的性质如压硬性等考虑不足。
传统弹性力学中的应变是微应变,但是土工试验中的土体变形量较高,因此完全按照传统应变去研究其土体性质无法充分考虑土体的内在特质。基于此,通过引入自然应变的概念、weibull统计分布规律、应变等价原理等,推导构建适用于膨胀土的损伤模型,并以损伤模型为基础进行冻融损伤分析等工作。
1. 基于自然应变的冻融损伤模型
1.1 自然应变适用性
以图1中所示的岩土体为分析模型,其满足均质各向同性的基本假定,受荷时,假定无新的损伤发生,其应力应变关系用HOOK定律进行描述,
dσ=Kdε, (1) 式中,
σ 为岩土材料截面应力(压为正),K为岩土体的刚度,ε 为当前应力状态下的应变量。根据弹性力学理论,岩土体的应变可为
dε=−dH/H, (2) 式中,H为当前应力状态下的试件高度。
式(2)计算的应变,Freed将其定义为自然应变或真实应变[10]。具体到应用阶段通过下式计算岩土体应变量更为广泛:
dε=−dH/H0, (3) 式中,H0为未受荷时的试件高度。对于高刚度材料,通常变形量较原来的形状较为微小,数量级相差在2个以上,此时式(2),(3)基本等效。式(3)由于更简洁的表达形式,在岩土工程建设中具有广泛的应用场景。但对于部分软岩、松散土体、断裂结构、多孔介质等,全部采用式(3)去分析其力学特征,则可能会阻碍对岩土本质特征的正确认识。
对于土体而言,土颗粒是高刚度材料,由于空气和水的存在,其形成的多介质土体则不完全满足式(3)的适用条件。尤其是对于亲水性矿物含量高的膨胀土而言,饱和时,会极大的削弱土体的刚度,干湿循环、冻融循环的条件下,会进一步导致其刚度的降低。因此在分析其力学特征时,式(2)具有更好的适用性。
1.2 基于应变等价原理的冻融损伤变量
岩土体的损伤产生主要有受荷和环境影响2部分。针对产生损伤岩土体,目前常用的损伤变量定义方法是基于应变等价原理实现,即可将损伤状态的材料在有效应力作用下的应变与无损状态下的应变进行等效,其无损状态的材料所受的应力称为名义应力。原理图如图1所示,表达式为
ε=σE′=σ′E, (4) σA=σ′A′, (5) 式中,
σ 为损伤状态下的有效应力,A为损伤状态下有效承载面积,E'为损伤状态下的岩土体模量,σ'等效无损状态下的名义应力,A'等效无损状态下有效承载面积,E为等效无损状态下的岩土体模量。损伤变量可用岩土体模量来表示,
D=1−E'/E。 (6) 更一般的情况下,岩土体的损伤可用岩土体变形模量的损伤来定义,目前多数文献均采用此定义[6, 8-9]。
基于式(2)的自然应变概念,可进一步建立岩土体刚度和传统方法变形模量的关系,
E=σ1−exp(−σ/K)。 (7) 将式(7)代入式(6),可得自然应变下的岩土体损伤变量为
D=1−K'/K, (8) 式中,K'为损伤状态下的土体刚度。式(6),(8)的形式基本一致,但是代表的是两种不同应变概念下的损伤定义方式。对于高刚度材料,E=K,传统的应变定义是自然应变定义下的一种特殊情况。
对于环境损伤Df(冻融循环、干湿循环)和受荷损伤Ds耦合的情况,总的损伤变量D可表示为
D=1−(1−Df)(1−Ds)。 (9) 1.3 受荷损伤模型构建
假设土体内部微元服从概率分布,可用Weibull分布定量描述微元强度的分布特征[11],
φ(F)=mF0(FF0)m−1e−(FF0)m=dDsdF, (10) 式中,φ(F)为土体微元强度概率分布函数,F0,m为Weibull分布的参数,F为微元的强度。
对式(10)进行积分可得,受荷损伤的损伤变量Ds与土体微元强度F的关系,
Ds=1−exp[−(F/F0)m]。 (11) 土体微元的强度准则采用D-P准则进行计算[12]。土体微元的强度F为
F=αI′1+√J′2, (12) 式中,
α 为材料参数,I′1 为有效应力状态下的应力第一不变量,J2'为有效应力状态下的偏应力第二不变量,α=2sinφ√3(3−sinφ) ,I′1=σ′1+2σ′3 ,J′2=13(σ′1−σ′3)2 ,} (13) 其中,φ土体微元内摩擦角,
σ′1 ,σ′3 为常规三轴条件下的竖向有效应力和有效围压。通过自然应变的概念,有效应力作用下的最大主应力方向的应变
ε′1 为ε′1=1−exp(−σ′1−σ′3K')(σ′1−σ′3)(σ′1−2μσ′3), (14) 式中,K'为土体的刚度,μ为泊松比。常规三轴试验时的固结过程,测得应力应变曲线主要是偏应力
(σ1−σ3) 和除去初始应变之后的应变ε1 ,基于应变等价原理,将有效应力进行转换,可将式(14)改写为σ1−σ3=−K0(1−Ds)ln(1−ε1), (15) 式中,K0为偏应力为0时的初始刚度。设孔隙水压力为PW,则联立式(12)~(15),土体基于自然应变的损伤本构模型如下:
σ1−σ3=−K0ln(1−ε1)exp[−(FF0)m], (16) 其中,
F={2sinφ√3(3−sinφ)[(σ1−σ3)+3(σ3−PW)]+1√3(σ1−σ3)}⋅−K0ln(1−ε1)σ1−σ3。 (17) 2. 基于遗传算法的模型参数确定
从式(16),(17)中可以发现,对于损伤参量无法直接测量的参数主要有初始刚度K0、土体微元内摩擦角φ、Weibull分布的参数F0和m。由于土体结构的复杂性,为避免误差累积,降低模型本身的精度,具体参数确定时采用智能算法遗传算法实现。遗传算法本质上属于生物进化算法,通过模拟生物进化中的优胜劣汰,对于具体的问题给出最优解[13]。对于给定的应力应变曲线,限定参数的取值范围,通过数值计算软件MATLAB编程,进行参数自动搜索,全局寻优,具体原理如图2所示。
3. 损伤模型验证
3.1 冻融循环试验数据
为进一步研究冻融循环后,膨胀土的损伤特性,利用Tang等[6]。的试验数据进行进一步的分析。试验土样充分饱水后进行冻融循环,-10℃保持12 h后,10℃保持12 h为1次冻融循环。试验采用常规三轴压缩试验,试验条件为固结排水,因此可不考虑超静孔隙水压力影响,试验围压为100 kPa。文献[6]采用的循环冻融次数NFT分别为0,1,3,7,11。试验结果如图3所示。
3.2 损伤模型验证
通过MATLAB编程,遗传算法可实现损伤模型参数的自动获取。通过实测的应变应变数据模拟,得到到实测值和模拟值的对比如图4所示,模拟参数如表1所示。从图4中可明显的发现,基于自然应变推导的损伤模型可以很好的模拟膨胀土在不同冻融循环下的应力应变关系,模拟精度较高。
表 1 损伤模型拟合参数Table 1. Fitting parameters for damage modelNFT K0/MPa φ/(°) F0 m 0 15.96 35.31 154.05 0.3174 1 14.89 35.83 153.06 0.3205 3 12.96 34.18 146.31 0.3200 7 12.42 35.30 145.52 0.3202 11 11.90 34.86 143.63 0.3209 3.3 损伤模型参数分析
从表1可以明显发现,4个参数中初始刚度K0、参数F0呈现明显的劣化特征,土体微元内摩擦角φ,m相对变化不大。Weibull分布中的参数F0和m通常代表着宏观平均强度和均布系数。由于土体结构本身的劣化,对土体形成一定的损伤,导致F0的减小。由于冻融过程中,土颗粒的团聚效应,均布系数m有略微增加的趋势,验证了论文中本构模型对于分析冻融损伤的合理性。冻融损伤产生的一个重要原因就是孔隙的增加,结构劣化,因此膨胀土的冻融损伤同样表现在初始刚度K0的明显降低。
选取受冻融影响明显的初始刚度K0和参数F0作为分析对象,绘制其与冻融循环次数的关系曲线如图5所示。从图5中可以明显看出,初始刚度K0和参数F0的劣化规律均满足一定的负指数关系,随着冻融次数的增加,劣化速度减慢,逐渐趋于一稳定值。
4. 讨论
通过建立的自然应变损伤模型,高精度的模拟了膨胀土冻融循环后的应力应变曲线,以及分析了冻融损伤对模型参数的影响。为进一步研究冻融损伤的损伤演化过程,利用式(9)绘制损伤因子D的变化曲线如图6所示。通过图6可以看出,在应力应变曲线的后期,不同冻融循环下的损伤因子基本上趋于一致。在应力应变曲线前期,随着冻融循环次数的增加,冻融损伤速度也进一步加快。但是这种加快的趋势随循环次数的增加开始减弱。
事实上,基于自然应变的膨胀土损伤模型优势不仅在于模拟精度上的提升,而且可以解决传统应变模式下对于土体本身性质认识的局限性。通过自然应变损伤模型,将非损伤变形和损伤变形进行一定区分,绘制图7的应力变形曲线。
通过图7可以发现,非损伤变形部分呈现逐渐密实的状态,即对应加载过程中不存在土体损伤时的变形状态,变形曲线特征基本上与土体的侧限压缩试验对应。三轴压缩时,由于损伤变形贡献了主要的变形量,因此曲线特征基本上与损伤变形曲线一致。
5. 结论
通过引入自然应变的概念,基于Weibull统计理论和应变等价原理,构建了损伤模型。在模型验证的基础上进一步分析膨胀土的冻融损伤特征。通过研究得出以下结论。
(1)构建的自然应变损伤模型可高精度的模拟膨胀土不同冻融循环下的应力应变曲线特征。遗传算法在模型验证中自动化程度高,可实现参数的自动寻优。
(2)膨胀土初始刚度K0和参数F0的冻融劣化明显,与冻融循环次数呈负指数关系。内摩擦角呈现出一定的波动下降趋势,均布系数m有略微增加的趋势。
(3)从不同冻融循环下的损伤因子演化曲线可以看出,冻融损伤主要发生在应力应变曲线前期,且随着冻融循环次数增加,应力应变曲线前期的损伤速度有加快趋势。应力应变后期损伤过程则基本一致。
(4)绘制非损伤变形和损伤变形,非损伤变形呈现逐渐压密的状态,损伤变形贡献了更多的变形量。其变形特征与侧限压缩试验曲线和三轴压缩曲线具有较好的对应性。
-
表 1 Aristoflex AVC物理力学性质
Table 1 Physical and mechanical of properties Aristoflex AVC
外观(20℃) 含水率/% Ph 值(1%浓度) 黏度/MPa(1%水溶液) 气味 白色粉末 7.53 4.0~6.0 48000~80000 特征性气味 表 2 透明度量化试验方案
Table 2 Schemes of transparency quantitative pilot
编号 浓度/% 厚度/mm 时间/d TQ1-1 1.00 100 7 TQ1-2 1.50 100 7 TQ1-3 2.00 100 7 TQ1-4 2.50 100 7 TQ2-2 2.50 100 14 TQ2-3 2.50 100 21 TQ2-4 2.50 100 28 TQ3-1 2.50 50 7 TQ3-2 2.50 150 7 TQ3-3 2.50 200 7 TQ3-4 2.50 250 7 TQ3-5 2.50 300 7 表 3 压缩试验加载路径
Table 3 Load paths of compression tests
土样编号 浓度/% 加载等级/ kPa CT-1 1.0 6.25→12.50→18.75→25.00→31.25→37.50→43.75→50.00→56.25 kPa CT-2 1.5 CT-3 2.0 25.00 kPa 下预压8h,卸载至6.25 kPa 持续8 h,稳定后加载 CT-4 2.5 表 4 不同浓度AVC透明黏土压缩试验参数
Table 4 Parameters of compression tests of different concentrations of AVC transparent soil
参数 e0 Cc Cs Cy/(cm2·s-1) CT1-1 52.5 16.3 0.7970 1.13×10-3 CT1-2 34.8 7.6 0.5480 1.09×10-3 CT1-3 26.0 4.4 0.4820 9.90×10-4 CT1-4 20.7 2.6 0.2998 8.04×10-4 表 5 AVC透明黏土十字板剪切试验方案
Table 5 Schemes of crossboard shear tests on AVC transparent clay
试验编号 浓度/% 附加荷载/ kPa 静置时间/d VST1-1 1.5 0 7 VST1-2 2.0 0 7 VST1-3 2.5 0 7 VST2-1 1.5 1 7 VST2-2 2.0 1 7 VST2-3 2.5 1 7 VST3-1 1.5 2 7 VST3-2 2.0 2 7 VST3-3 2.5 2 7 VST4-1 2.5 2 7 VST4-2 2.5 2 1 VST4-3 2.5 2 10 VST4-4 2.5 2 14 -
[1] ALLERSMA H. Photo-elastic stress analysis and strain in simple shear: IUTAM Conference on Deformation and Failure of Granular Materials[C]//Goteborg, Swed, 1982.
[2] SADEK S, ISKANDER M G, et al. Geotechnical properties of transparent silica[J]. Canadian Geotechnical Journal, 2002, 39(1): 111-124. doi: 10.1139/t01-075
[3] 孔纲强, 刘璐, 刘汉龙, 等. 玻璃砂透明土变形特性三轴试验研究[J]. 岩土工程学报, 2013, 35(6): 1140-1146. KONG Gang-qiang, LIU Wei, LIU Han-long, et al. Three-axis test study of transparent soil deformation characteristics of glass sand[J]. Journal of Geotechnical Engineering, 2013, 35(6): 1140-1146. (in Chinese)
[4] 孔纲强, 孙学谨, 肖扬, 等. 透明土与标准砂压缩变形特性对比试验研究[J]. 岩土工程学报, 2016, 38(10): 1895-1903. doi: 10.11779/CJGE201610020 KONG Gang-qiang, SUN Xue-ming, XIAO Yang, et al. A comparative experimental study of the deformation characteristics of transparent soil and standard sand compression[J]. Journal of Geotechnical Engineering, 2016, 38(10): 1895-1903. (in Chinese) doi: 10.11779/CJGE201610020
[5] 孔纲强, 孙学谨, 李辉, 等. 孔隙液体对玻璃砂透明土强度特性影响研究[J]. 岩土工程学报, 2016, 38(2): 377-384. KONG Gang-qiang, SUN Xue-cheng, LI Hui, et al. Study on the effect of pore liquid on the strength characteristics of transparent soil in glass sand[J]. Journal of Geotechnical Engineering, 2016, 38(2): 377-384. (in Chinese)
[6] 李亮. 透明土合成及物理力学特性研究[D]. 杭州: 浙江大学, 2014. LI Liang. Study of Transparent Soil Synthesis and Physical and Physical and Physical Properties[D]. Hangzhou: Zhejiang University, 2014. (in Chinese)
[7] WALLACE J F, RUTHERFORD C J. Geotechnical properties of laponite rd®[J]. Geotechnical Testing Journal, 2015, 38(5): 20140211. doi: 10.1520/GTJ20140211
[8] HAKHAMANESHI M, BLACK J A. Shear strength of transparent gelita-effect of mixture ratio, displacement rate and over-consolidation ratio[C]//Chicago Illinois, US, 2016.
[9] HAKHAMANESHI M, BLACK J A, TATARI A. Optical characterization of transparent gelita using the modulation transfer function[C]//Chicago Illinois, US, 2016.
[10] STANIER S A. Modelling the Behaviour of Helical Screw Piles[D]. University of Sheffield, Department of Civil and Structural Engineering, 2011.
[11] BLACK J A, TAKE W A. Quantification of optical clarity of transparent soil using the modulation transfer function[J]. Geotechnical Testing Journal, 2015(5): 588-602.
[12] 房后国. 深圳湾结构性淤泥土固结机理及模型研究[D]. 吉林: 吉林大学, 2005. FANG Hou-guo. Study on the Solidification of Structural Silt in Shenzhen Bay[D]. Jilin: Jilin University, 2005. (in Chinese)
[13] 吴雪婷. 温州浅滩淤泥固结系数与固结应力关系研究[J]. 岩土力学, 2013, 34(6): 1675-1680. WU Xue-xuan. Study on the relationship between shoal sludge solidification coefficient and solidification stress in Wenzhou[J]. Rock and Soil Machanics, 2013, 34(6): 1675-1680. (in Chinese)
[14] 张长生, 高明显, 強小俊. 深圳后海湾海相淤泥固结系数变化规律研究[J]. 岩土工程学报, 2013, 35(增刊1): 247-252. ZHANG Chang-sheng, GAO Ming-xian, QIANG Xiao-jun. Study on the law of change of the solidification coefficient of sea silt in the back bay of Shenzhen[J]. Journal of Geotechnical Engineering, 2013, 35(S1): 247-252. (in Chinese)
-
期刊类型引用(5)
1. 谢柏涵,孔元元,张学飞,孙东彦,王智慧. 冻融循环作用下盐渍土损伤模型的建立及验证. 人民长江. 2025(02): 183-189+200 . 百度学术
2. 孙超,谷鸿雁,郭浩天. 冻融循环次数对粉质黏土强度劣化影响的试验研究. 吉林建筑大学学报. 2024(05): 66-71 . 百度学术
3. 吕思清,朱杰兵,汪斌,祝永锁. 冻融荷载耦合作用下含开口裂隙砂岩宏细观损伤模型研究. 岩石力学与工程学报. 2023(05): 1124-1135 . 百度学术
4. 李新宇,凌贤长,曲娜. 考虑温度效应的冻结膨胀土统计损伤模型. 吉林大学学报(工学版). 2023(08): 2339-2349 . 百度学术
5. 张雁,薛晓雨,杨晓蕴. 低温状态下纤维加筋膨胀土数值模拟试验. 中国科技论文. 2022(12): 1351-1357 . 百度学术
其他类型引用(5)