Development of miniature dynamic triaxial apparatus for microct scanning and its experimental verification
-
摘要: 土体的微观结构对其宏观力学性质有着很大的影响,为了能更深入研究土体在循环荷载作用下的微观结构演化与其宏观力学特性之间的联系,研制了一种适用于工业显微CT(computerized tomography)扫描的微型动三轴仪。该仪器主要由控制采集箱系统、主机系统、气水转换器组成。其具有体积较小、试样小型化(直径×高=10 mm×20 mm)、便于操作、兼容性强等优点,可在无需改装已有CT设备前提下对试样进行旋转扫描。以丰浦砂干砂样为材料,采用微型和常规动三轴仪进行了不同围压(50,100,200,300 kPa)和循环荷载作用下的不固结不排水试验,并对测得的初始动弹性模量E0,动剪切模量Gd,阻尼比
λ 进行了对比分析。结果表明:微型动三轴仪测得的初始动弹性模量E0,动剪切模量Gd,阻尼比λ与常规动三轴仪测得的结果基本相近。试验结果误差较小,验证了该仪器的可靠性。Abstract: The microstructure of soil has a great influence on the macroscopic mechanical properties. Therefore, to study the relationship between the macroscopic mechanical properties and the microstructure evolution of soil under the cyclic loading, a miniature dynamic triaxial apparatus is developed, which is suitable for industrial micro CT (computerized tomography) scanning. The apparatus is mainly composed of the acquisition control system, loading device system, and air-liquid converter. It has the advantages of small volume, miniaturization of the specimen (diameter × height=10 mm×20 mm), convenient to operate, strong compatibility, and the specimen can be rotating-scanned without refitting the existing CT equipment. By using the miniature dynamic triaxial apparatus and conventional dynamic triaxial apparatus, a series of unconsolidated-undrained triaxial tests on Toyoura sand under cyclic loading and different confining pressures (50, 100, 200 and 300 kPa) are carried out, and the initial dynamic elastic modulus (E0), dynamic shear modulus (Gd) and damping ratio (λ ) are comparatively analyzed. The results show that the experimental curves of the initial dynamic elastic modulus (E0), dynamic shear modulus (Gd) and damping ratio (λ ) exhibit a similar evolution. The difference of experimental data is small, thus, the reliability of the miniature dynamic triaxial apparatus is verified. -
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 日本丰浦砂物理性能参数
Table 1 Physical parameters of Toyoura sand
试样 土粒相对密度Gs 平均粒径D50/mm 不均匀系数Cu 最大孔隙比emax 最小孔隙比emin 丰浦砂 2.65 0.2 1.32 0.977 0.597 表 2 试验方案
Table 2 Experimental schemes
序号 试验仪器 试样 干密度ρd/(g·cm-3) 频率f/Hz 围压 σ3 /kPa1 微型 丰浦砂 1.483 0.2 50,100,200,300 2 常规 丰浦砂 1.483 0.2 50,100,200,300 表 3 动剪切模量与动剪应变
Table 3 Dynamic shear modulus and dynamic shear strains
微型动三轴仪 常规三轴仪 相对误差/% 动剪应变γd/% 动剪切模量Gd/kPa 动剪应变γd/% 动剪切模量Gd/kPa 0.012 508 0.012 534.4 4.9 0.031 408 0.031 387.2 5.4 0.062 264.8 0.062 286.4 7.5 0.12 217.6 0.12 230 5.4 0.25 170 0.25 177.49 4.2 0.5 133.8 0.5 140.69 4.9 2 69.43 2 70.31 1.3 表 4 初始动弹性模量
Table 4 Initial dynamic elastic moduli
围压 初始动弹性模量E0/MPa 相对误差/% 微型动三轴仪 常规三轴仪 50 1.43 1.37 4.2 100 2.50 2.38 4.8 200 14.29 13.70 4.1 300 17.86 17.83 0.2 表 5 不同围压条件下阻尼比
Table 5 Damping ratios under different confining pressures
围压 阻尼比 λ 相对误差/% 微型动三轴仪 常规动三轴仪 50 0.2366 0.2349 0.6 100 0.2429 0.2527 3.8 200 0.2572 0.2745 6.3 300 0.2717 0.2831 4.0 -
[1] 蒋明镜. 现代土力学研究的新视野—宏微观土力学[J]. 岩土工程学报, 2019, 41(2): 195-254. JIANG Ming-jing. New paradigm for modern soil mechanics: Geomechanics from micro to macro[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(2): 195-254. (in Chinese)
[2] AMOROSI A, RAMPELLO S. An experimental investigationinto the mechanical behaviour of a structured stiff clay[J]. Géotechnique, 2007, 57(2): 153-166. doi: 10.1680/geot.2007.57.2.153
[3] JIANG M J, YU H S, HARRIS D. A novel discrete model for granular material incorporating rolling resistance[J]. Computers and Geotechnics, 2005, 32(5): 340-357. doi: 10.1016/j.compgeo.2005.05.001
[4] COOP M R, SORENSEN K K, FREITAS T B, et al. Particle breakage during shearing of a carbonate sand[J]. Géotechnique, 2004, 54(3): 157-163. doi: 10.1680/geot.2004.54.3.157
[5] JIANG M J, SUN Y G, LI L Q, et al. Contact behavior of idealized granules bonded in two different interparticledistances: an experimental investigation[J]. Mechanics of Materials, 2012, 55(14): 1-15.
[6] 王艳丽, 程展林. CT扫描技术在我国土工试验中的应用研究进展[J]. 地震工程学报, 2015, 37(增刊1): 35-39. WANG Yan-li, CHENG Zhan-lin. Progress in the application of CT scanning technology in Chinese soil tests[J]. China Earthquake Engineering Journal, 2015, 37(S1): 35-39. (in Chinese)
[7] 胡波, 龚壁卫, 童军, 等. 单向冻结冻融循环三轴仪的研制及其应用[J]. 长江科学院院报, 2015, 32(2): 128-132. HU Bo, GONG Bi-wei, TONG Jun, et al. Development and application of unilateral freezing triaxial apparatus for freeze-thaw cycle[J]. Journal of Yangtze River Scientific Research Institute, 2015, 32(2): 128-132. (in Chinese)
[8] 郑剑锋, 马巍, 赵淑萍, 等. 与CT机配套的便携式控温压力舱的研制与应用[J]. 冰川冻土, 2014, 36(4): 928-933. ZHENG Jian-feng, MA Wei, ZHAO Shu-ping, et al. The portable three-terminal controlling temperature apparatus matched to computerized tomography: development and application[J]. Journal of Glaciology and Geocryology, 2014, 36(4): 928-933. (in Chinese)
[9] 程壮, 王剑锋. 用于颗粒土微观力学行为试验的微型三轴试验仪[J]. 岩土力学, 2018, 39(3): 1123-1129. CHENG Zhuang, WANG Jian-feng. A mini-triaxial apparatus for testing of micro-scale mechanical behavior of granular soils[J]. Rock and Soil Mechanics. 2018, 39(3): 1123-1129. (in Chinese)
[10] 蒋明镜, 李光帅, 曹培, 等. 用于土体宏微观力学特性测试的微型三轴仪研制[C]//第29届全国土工测试学术研讨会, 2020, 呼和浩特. JIANG Ming-jing, LI Guang-shuai, CAO Pei, et al. Developmentof the miniature triaxial apparatus for testing of macro and micromechanical behavior of soils[C]//The 29th National Conference on Geotechnical Testing, 2020, Huhhot. (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)