Control Factors and Assessment Technique of Relative Density Using Pluviation Method for Saturated Model
-
摘要: 为满足动力离心液化试验相对密度(Dr)低的制样需求,实现砂雨法饱和模型制样相对密度的准确控制,建立了设备制模稳定性评价方法,通过自主研制一套适于饱和模型制样的鸭嘴式砂雨法装置,开展了三组干砂/饱和砂模型对比试验;通过分析出砂口尺寸、落距、移动速度等控制要素的影响,对新型装置制样性能进行了验证;采用微型动力触探仪测试饱和模型不同位置及深度的Dr空间分布,给出了模型均匀稳定性评价方法;建立描述砂雨法制样过程流速变化的数学模型和推导表达式,提出了控制稳态Dr的归一化标准. 研究结果表明:3 mm为低密实度制样的最佳出砂口尺寸;饱和模型Dr随水中落距的变化率为空中落距变化率的3.5倍,水中落距是饱和制样密实度的主导控制要素;出砂口移动速度最高达到颗粒落速的31%,对低密实度制样影响不可忽略;设备移速与落距对颗粒流速及试样Dr大小具有决定作用.Abstract: Accurate control of relative density (Dr) in saturated pluviation method is the fundamental requirement for low relative density sample preparation of dynamic centrifuge liquefaction tests. To achieve an evaluation method on pluviation stability, a set of duckbill pluviator which fits saturated modeling was developed. Three groups of dry and saturated sand model contrast tests were carried out to identify the impacts of nozzle size, drop distance and movement velocity, aiming to assess the capability of the developed sampling device. A mini dynamic penetrometer facility was applied to measure the spatial distribution of the saturated model’s Dr, giving a homogeneity evaluation of the models. A particle-flow-velocity model for pluviation was derived to give a normalize criteria for density stable control. The results show that the optimal size of nozzle is 3mm which satisfies the requirements of low compactness modeling. The change rate of Dr with thickness of water is 3.5 times of that with air drop height, indicating that thickness of water layer is the dominant factor on Dr of saturated model. The nozzle movement velocity is up to 31% of particle velocity which is non-trivial on the low density sample preparation. The travelling speed and drop height play a decisive role in particle flow velocity and Dr.
-
Key words:
- relative density (Dr) /
- model test /
- pluviation /
- saturated model /
- sampling factor
-
目前,无砟轨道是高速铁路和城市轨道交通的主要轨道结构形式[1]. 无砟轨道作为一种多层薄板体系,层间界面为结构薄弱环节,在列车荷载和环境共同作用下,结构层间易产生离缝[2-4]. 现场调研发现,在多雨地区或排水不良地段,水容易进入轨道结构层间,形成离缝冒浆等病害,加速无砟轨道结构劣化[5-8]. 因此,研究列车荷载与水耦合作用下无砟轨道水致损伤,首要的是明确结构层间水流特性.
列车荷载作用下无砟轨道层间离缝内水流特性问题,实质上是多层薄板体系内层间裂缝内水流特性问题. 调研发现,无砟轨道层间脱空病害主要发生在板边,实际无砟轨道层间脱空区域是不规则的形状,为方便分析计算,将层间脱空区域简化为矩形等规则形状. 基于此简化模型,相关学者设计了不同层间脱空状态的混凝土试件,结合现有的加载装置开展动水压力试验. 徐桂弘等[9]设计了与原型相似的几何模型试件,离缝位于试件板边中部,离缝尺寸为300 mm × 500 mm × 3 mm (深度 × 宽度 × 开口量). 试验过程中,施加的荷载为(60±35) kN,分析不同荷载频率下层间离缝内动水压力的变化. 发现荷载频率为15 Hz时,层间离缝动水压力最大,最大值为4.968 kPa. 该试验将混凝土层间脱空试件直接放置到水箱中,由于缺乏必要的防振动措施,试件在加载过程中会发生较为明显的振动甚至跳动,使试验测试结果与理论计算结果存在较大差异. 曹世豪等[10]在原试验装置上增设了必要的防振动措施,设计的裂缝尺寸为250 mm × 600 mm × 3 mm(深度 × 宽度 × 开口量),试验荷载为(25±10) kN,最大荷载频率为5 Hz,测得层间离缝最大动水压力为0.177 kPa,与理论计算结果较为接近. 但列车荷载具有频率高的特点,我国高速铁路设计时速在250~350 km/h,以时速250 km为例,一个转向架相邻轮载作用间隔时间为0.04 s,即荷载频率为25 Hz.
以上试验研究一方面受试验加载设备的限制,采用的荷载条件难以反映列车荷载频率高的特点,另一方面,由于设计的试件尺寸较大,在荷载作用下混凝土层间脱空区域的变形量很小,不能较好地模拟实际列车荷载作用下板边脱空区域的变形行为,也不能很好地反映实际列车荷载作用下无砟轨道层间动水压力变化. 因此,有必要基于无砟轨道层间脱空变形特性,设计能够有效模拟高频荷载作用下无砟轨道层间裂缝内水流运动的试验装置,为进一步分析列车荷载与水耦合作用下无砟轨道层间裂缝内水流特性提供试验基础.
1. 无砟轨道层间脱空力学行为
采用ANSYS Workbench建立无砟轨道层间离缝冒浆平面简化模型,网格划分如图1所示,对结构进行谐响应分析,荷载平均施加在承轨台上,荷载幅值取为钢轨支点力(扣件对轨道的作用力)的幅值[11],荷载频率取25 Hz;在不考虑弹性基础的影响下,支承层底部采取全约束. 各部件计算参数如表1所示. 在此基础上,在脱空区域的最内侧(尖端)、中点、边缘处分别设置测点1、2、3,监测不同测点处的垂向相对位移(脱空区域道床板底面与支承层表面垂向位移差值),分析不同脱空深度和开口量对无砟轨道层间脱空区域变形特性的影响.
表 1 材料参数Table 1. Material parameters部件 强度 弹性模量/GPa 泊松比 密度/(kg·m−3) 轨枕 C60 35.5 0.20 2500 道床板 C40 32.5 0.20 2450 支承层 C15 22.0 0.18 2400 1.1 脱空深度的影响
由于层间脱空区域扩展到承轨台下方时,会导致结构受力状态的改变,故首先分析不同脱空深度下脱空区域的变形行为. 当脱空深度L分别为0.10、0.15、0.20、0.25、0.30、0.35 m和0.40 m[10],开口量h=8 mm时,提取不同监测点垂向相对位移,负值表示向下运动,如图2所示. 由图可知:随着脱空深度的增大,脱空区域不同位置垂向相对位移近似线性增大,在测点1处变化幅度很小,垂向相对位移接近于0;在测点3处变化幅度最大,脱空深度增大4倍,垂向相对位移增大4倍~5倍;在不同脱空深度下,脱空区域垂向相对位移在最内侧最小,在边缘处最大,脱空区域垂向相对位移沿脱空边缘方向近似线性增大.
1.2 开口量的影响
取L = 0.30 m,分析开口量对脱空区域变形行为的影响,h = 2,4,6,8,10 mm[12],计算结果如图3所示. 由图3可知:随着开口量增大,脱空区域不同测点垂向相对位移近似相等,表明开口量的变化对脱空区域的变形行为影响很小;在不同开口量下,脱空区域垂向相对位移在最内侧最小,在边缘处最大,最大垂向相对位移均在0.220 mm内. 综合1.1、1.2节计算分析发现,当L = 0.30 m时,脱空区域最大垂向相对位移为0.201 mm,与文献[13]现场测试结果接近.
2. 层间脱空模拟装置设计
本文对轨道结构进行简化,将无砟轨道简化为道床板与支承层结构. 实际无砟轨道层间离缝形状是不规则的,假定离缝为规则的形状. 基于列车荷载作用下无砟轨道层间脱空区域的变形特性,结合课题组研制的高频加载装置,设计出如图4所示的模拟装置. 模拟装置由装配式钢板工装、预制水槽、高频加载装置、螺栓、混凝土板等构成,混凝土上、下板之间的空隙作为层间脱空区域. 模拟原理:装配式钢板工装的上板可视为悬臂结构,在高频荷载作用下,上板边缘垂向位移最大. 装配式钢板工装上下板通过螺栓与混凝土板联结,当钢板上板受力变形时,混凝土上板随之产生相同的位移,且保证混凝土边缘处的垂向位移最大,层间脱空区域垂向相对位移沿脱空边缘方向近似线性增大,符合列车荷载作用下脱空区域的变形特性.
根据第1节的计算分析,为模拟现场实际荷载作用下脱空区域的变形特征,设计层间脱空离缝深度为300 mm,宽度为200 mm,开口量在0~15 mm内调节. 混凝土上板尺寸为300 mm × 200 mm × 50 mm,混凝土下板尺寸为350 mm × 200 mm × 50 mm. 混凝土上、下板试件的强度等级为C40,试件内预埋螺栓与装配式钢板工装连接. 装配式钢板工装由10 mm的钢板焊接制作. 预制水槽尺寸为1600 mm × 800 mm × 150 mm. 水槽底面预埋6个螺栓,通过螺栓与装配式钢板工装连接,以固定结构,减少试验过程中结构振动对试验结果的影响. 高频加载装置放置到加载装置安置板上,荷载形式为正弦曲线,荷载频率可达50 Hz.
利用ANSYS Workbench对钢板工装进行谐响应分析,计算不同荷载作用下结构上板的垂向相对位移,与现场测试结果对比,确定试验实际加载值. 装置材料选用结构钢,弹性模量为 210.0 GPa,泊松比为 0.30. 装置底部设置为全约束,侧面约束纵向和横向位移. 当荷载频率为25 Hz,不同荷载幅值作用下结构上板边缘垂向相对位移如图5所示. 由图5可知,在荷载0.5~1.5 kN下,随着荷载增加,结构上板边缘位移近似线性增大,垂向相对位移值均在现场实测结果之内. 当荷载幅值为1.1 kN时,结构上板边缘垂向相对位移为0.220 mm,较为接近L = 0.30 m,h = 8 mm的工况下脱空区域边缘的垂向相对位移.
组装试验装置,在荷载幅值为1.1 kN、荷载频率为25 Hz的条件下,采用高精度接触式位移计测试板端位移,对比试验与理论计算结果,如图6所示. 在高频加载装置作用下,试验测试板端最大垂向相对位移约为0.20 mm,与理论计算结果0.22 mm接近,表明设计的层间脱空模拟装置可以用于模拟测试高频荷载作用下无砟轨道层间脱空区域的变形行为.
3. 无砟轨道层间动水压力试验
3.1 试验概况
结合前文设计的层间脱空装置,开展无砟轨道层间离缝动水压力试验. 向预制水槽中注入20 ℃的自来水,保证水位高于离缝上表面. 层间脱空区域后部及两侧利用防水胶黏接,保证离缝内的水仅沿着出口方向流动. 压力传感器测试端通过软管连接皮托管深入离缝内部,另一端与采集软件Smart Sensor连接,实时监测测点水压力. 各个测点等间距分布,相邻测点间距为50 mm,测点 ① 位于层间离缝最内侧,测点 ⑦ 位于离缝开口处,具体测量位置如图7所示. 由于脱空区域变形特征与脱空深度有关,而且为了方便测试,适应压力传感器尺寸,试验时开口量取为8、14 mm. 按照表2的4种设计工况进行试验.
表 2 试验工况Table 2. Test working conditions试验工况 开口量/mm 加载频率/Hz 目的 1 8 15 加载频率对动水压力的影响 2 8 25 3 8 25 开口量对动水压力的影响 4 14 25 3.2 试验结果及分析
3.2.1 动水压力时间分布特性
在加载频率25 Hz、层间离缝开口量8 mm的工况,提取加载后的试验数据,分析不同监测点的动水压力时程分布规律,如图8所示. 在高频荷载作用下,不同测点动水压力变化趋势总体呈现正弦形式,其值正负交替变化,随着荷载达到最大值后达到峰值. 加载时,由于道床板混凝土底面靠近支承层混凝土表面,使得层间脱空区域内水被排出,此时脱空区域内产生正压力;卸载时,由于道床板混凝土底面远离支承层混凝土表面,使得层间脱空区域外水被吸回脱空区域内,此时脱空区域内产生负压力. 以测点 ①为例,最大正压为15.794 kPa,最大负压为−16.766 kPa,压力差为32.560 kPa.
3.2.2 动水压力空间分布特性
在加载频率25 Hz、层间离缝开口量8 mm的条件下,提取不同测点动水压力的最大正压和最大负压,如图9所示. 动水压力沿离缝长度方向近似线性分布,在测点 ① 处动水压力最大,监测点越靠近离缝出口,动水压力越小.
3.2.3 加载频率对动水压力的影响
在离缝为8 mm的条件下改变加载频率大小,对比不同荷载频率下测点 ① 的动水压力时程曲线,如图10所示,当荷载频率为25 Hz时,该测点的最大动水压力约为15.794 kPa;当加载频率为15 Hz时,该处的最大动水压力约为1.646 kPa,频率增大10 Hz,最大动水压力增大接近10倍,说明荷载频率是无砟轨道层间离缝水压力的重要影响因素之一. 测得不同监测点动水压力沿离缝分布,如图11所示,随着加载频率增大,动水压力会大幅度增大. 在不同频率下,离缝内动水压力均在离缝深处最大,沿离缝出口方向减小.
3.2.4 离缝开口量对动水压力的影响
在加载频率为25 Hz时,对比不同离缝开口量下测点 ① 的动水压力时程曲线,如图12所示. 在开口量为8 mm时,测点 ① 的水压力为15.794 kPa;开口量为14 mm时,该处的水压力为8.320 kPa,说明离缝开口量是无砟轨道层间离缝水压力的重要影响因素之一. 测得不同监测点动水压力沿离缝分布如图13所示. 由图可知:随着离缝开口量的增大,动水压力会减小;在不同开口量下,离缝内动水压力均在离缝深处最大,沿离缝出口方向减小.
3.3 试验模型验证
根据图4试验模型,利用ANSYS Workbench与CFX建立双向瞬态流固耦合模型,将层间脱空模拟装置和混凝土试件设置为固体域;混凝土层间脱空区域设为流体域,并假定水流状态为层流;将混凝土试件层间离缝上下表面设置为流固耦合面,便于固体域与流体域进行数据传递. 根据表1设置材料参数,计算相同荷载条件下(25 Hz)层间裂纹不同测点处的动水压力. 不同测点动水压力理论计算结果与试验结果对比如图14所示. 试验结果与理论计算结果存在一定误差,可能原因是试验过程存在局部水压力损失,但二者整体趋势接近,表明试验模型具有一定的有效性.
4. 结 论
1) 随着脱空深度L的增大,脱空区域垂向相对位移显著增大,但开口量h对脱空区域垂向相对位移影响很小. 当L = 0.3 m时,脱空区域最大垂向相对位移为0.201 mm与现场测试结果接近. 脱空深度L是影响脱空区域垂向相对位移的主要因素.
2) 当荷载频率为25 Hz,幅值为1.1 kN时,设计的层间脱空模拟装置板端最大垂向相对位移为0.220 mm,与现场测试结果接近,表明可以利用该装置模拟实际列车荷载作用下无砟轨道层间脱空区域水流特性.
3) 在高频荷载作用下,无砟轨道混凝土层间离缝内动水压力整体变化趋势呈现正弦形式,当荷载值达到最大时,离缝内部水压力达到峰值;相同条件下,离缝内动水压力在离缝深处最大,沿离缝出口方向减小. 当荷载频率为25、15 Hz时,层间离缝内最大动水压力分别为15.794、1.646 kPa,随着荷载频率增大,动水压力增大. 当离缝开口量为8、14 mm时,层间离缝内最大动水压力分别为15.794、8.320 kPa,随着开口量增大,动水压力减小.
致谢:西南交通大学原创性科研仪器项目(XJ2021KJZK035)的资助.
-
HUANG A B, CHANG W J, HSU H H, et al. A mist pluviation method for reconstituting silty sand specimens[J]. Engineering Geology, 2015, 188(1): 1-9. 马险峰,孔令刚,方薇,等. 砂雨法试样制备平行试验研究[J]. 岩土工程学报,2014,36(10): 1791-1801. doi: 10.11779/CJGE201410005MA Xianfeng, KONG Linggang, FANG Wei, et al. Parallel tests on preparation of samples with sand pourer[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1791-1801. doi: 10.11779/CJGE201410005 李浩,罗强,张正,等. 砂雨法制备砂土地基模型控制要素试验研究[J]. 岩土工程学报,2014,36(10): 1872-1878. doi: 10.11779/CJGE201410015LI Hao, LUO Qiang, ZHANG Zheng, et al. Experimental study on control element of sand pourer preparation of sand foundation model[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1872-1878. doi: 10.11779/CJGE201410015 KHARI M, KASSIM K, ADNAN A. Sand samples’ preparation using mobile pluviator[J]. Arabian Journal for Science and Engineering, 2014, 39(10): 6825-6834. doi: 10.1007/s13369-014-1247-8 GADE V K, DASAKA S M. Development of a mechanized traveling pluviator to prepare reconstituted uniform sand specimens[J]. Journal of Materials in Civil Engineering, 2016, 28(2): 1-9. RAGHUNANDAN M, JUNEJA A, HSIUNG B. Preparation of reconstituted sand samples in the laboratory[J]. International Journal of Geotechnical Engineering, 2012, 6(1): 125-131. doi: 10.3328/IJGE.2012.06.01.125-131 IQBAL W. Sand pluviation technique[M]. [S.l.]: LAP Lambert Academic Publishing, 2012: 1-10 袁晓铭,曹振中,孙锐,等. 汶川8.0级地震液化特征初步研究[J]. 岩石力学与工程学报,2009,28(6): 1288-1296. doi: 10.3321/j.issn:1000-6915.2009.06.026YUAN Xiaoming, CAO Zhenzhong, SUN Rui, et al. Preliminary research on liquefaction characteristics of Wenchuan 8.0 earthquake[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(6): 1288-1296. doi: 10.3321/j.issn:1000-6915.2009.06.026 徐锡伟,陈桂华,于贵华,等. 芦山地震发震构造及其与汶川地震关系讨论[J]. 地学前缘,2013,20(3): 11-20.XU Xiwei, CHEN Guihua, YU Guihua, et al. Seismogenic structure of Lushan earthquake and its relationship with Wenchuan earthquake[J]. Earth Science Frontiers, 2013, 20(3): 11-20. 蒋关鲁,刘先峰,张建文,等. 高速铁路液化土地基加固的振动台试验研究[J]. 西南交通大学学报,2006,41(2): 190-196. doi: 10.3969/j.issn.0258-2724.2006.02.011JIANG Guanlu, LIU Xianfeng, ZHANG Jianwen, et al. Shaking table test of composite foundation reinforcement of saturated silty soil ground for high speed railway[J]. Journal of Southwest Jiaotong University, 2006, 41(2): 190-196. doi: 10.3969/j.issn.0258-2724.2006.02.011 刘汉龙. 土动力学与土工抗震研究进展综述[J]. 土木工程学报,2012,45(4): 148-164.LIU Hanlong. A review of recent advances in soil dynamics and geotechnical earthquake engineering[J]. China Civil Engineering Journal, 2012, 45(4): 148-164. BOULANGER R W, MONTGOMERY J. Nonlinear deformation analyses of an embankment dam on a spatially variable liquefiable deposit[J]. Soil Dynamics and Earthquake Engineering, 2016, 91: 222-233. doi: 10.1016/j.soildyn.2016.07.027 程新俊,景立平,崔杰,等. 不同场地沉管隧道振动台模型试验研究[J]. 西南交通大学学报,2017,52(6): 1113-1120. doi: 10.3969/j.issn.0258-2724.2017.06.011CHENG Xinjun, JING Liping, CUI Jie, et al. Research of shaking table model tests on immersed tunnels under different conitions[J]. Journal of Southwest Jiaotong University, 2017, 52(6): 1113-1120. doi: 10.3969/j.issn.0258-2724.2017.06.011 苏雷,凌贤长,唐亮,等. 可液化场地桥梁群桩基动力反应振动台试验研究[J]. 防灾减灾工程学报,2015,35(2): 186-191.SU Lei, LING Xianzhang, TANG Liang, et al. Shaking table tests on dynamic responses of pile group foundations for bridge in liquefiable ground[J]. Journal of Disaster Prevention and Mitigation Engineering, 2015, 35(2): 186-191. 冯研,蒋关鲁,陈伟志,等. 离心模型试验预测复合地基沉降的精度[J]. 西南交通大学学报,2014,49(1): 105-110. doi: 10.3969/j.issn.0258-2724.2014.01.017FENG Yan, JIANG Guanlu, CHEN Weizhi, et al. Accuracy of settlement prediction of composite foundation by centrifuge model tests[J]. Journal of Southwest Jiaotong University, 2014, 49(1): 105-110. doi: 10.3969/j.issn.0258-2724.2014.01.017 涂杰文,刘红帅,汤爱平,等. 基于离心振动台的堆积型滑坡加速度响应特征[J]. 岩石力学与工程学报,2015,34(7): 1361-1369.TU Jiewen, LIU Hongshuai, TANG Aiping, et al. Acceleration response of colluvial landslide based on centrifugal shaking table test[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(7): 1361-1369. 王维早,许强,郑光,等. 强降雨诱发缓倾堆积层边坡失稳离心模型试验研究[J]. 岩土力学,2016,37(1): 87-95.WANG Weizao, XU Qiang, ZHENG Guang, et al. Centrifugal model tests on sliding failure of gentle debris slope under rainfall[J]. Rock and Soil Mechanics, 2016, 37(1): 87-95. 王永志,WILSON D W,KHOSRAVI M,等. 动力离心模型试验循环剪应力-剪应变反演方法对比[J]. 岩土工程学报,2016,38(2): 271-277.WANG Yongzhi, WILSON D W, KHOSRAVI M, et al. Evaluation of cyclic shear stress-strain using inverse analysis techniques in dynamic centrifuge tests[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(2): 271-277. 徐超,贾斌,罗玉珊,等. 圬工与加筋土组合式挡墙离心模型试验[J]. 同济大学学报:自然科学版,2015,43(3): 379-385.XU Chao, JIA Bin, LUO Yushan, et al. Centrifuge model tests of behavior of masonry and reinforced soil composite retaining wall[J]. Journal of Tongji University:Natural Science, 2015, 43(3): 379-385. 蔡正银,张晨,黄英豪. 冻土离心模拟技术研究进展[J]. 水利学报,2017,48(4): 398-407.CAI Zhengyin, ZHANG Chen, HUANG Yinghao. A review on the development of geotechnical centrifuge modeling technique on frozen ground engineering[J]. Journal of Hydraulic Engineering, 2017, 48(4): 398-407. 郑健,李育超,陈云敏. 底泥固结对污染物运移影响的超重力离心试验模拟[J]. 浙江大学学报(工学版),2016,50(1): 8-15. doi: 10.3785/j.issn.1008-973X.2016.01.002ZHENG Jian, LI Yuchao, CHEN Yunmin. Centrifuge test modeling of impact of sediment consolidation on contaminant transportation[J]. Journal of Zhejiang University (Engineering Science), 2016, 50(1): 8-15. doi: 10.3785/j.issn.1008-973X.2016.01.002 安鹏,邢义川,张爱军,等. 基于离心模型试验的深厚湿陷性黄土自重湿陷性评价研究[J]. 四川大学学报(工程科学版),2016,48(6): 23-30.AN Peng, XING Yichuan, ZHANG Aijun, et al. Research on evaluation of self-weight collapsibility for large-thickness collapsible loess using centrifugal model test[J]. Journal of Sichuan University (Engineering Science Edition), 2016, 48(6): 23-30. TOBITA T, ASHINO T, REN J, et al. Kyoto University LEAP-GWU-2015 tests and the importance of curving the ground surface in centrifuge modelling[J]. Soil Dynamics and Earthquake Engineering, 2017, 113: 650-662. 周燕国,李永刚,丁海军,等. 砂土液化后再固结体变规律表征与离心模型试验验证[J]. 岩土工程学报,2014,36(10): 1838-1845. doi: 10.11779/CJGE201410011ZHOU Yanguo, LI Yonggang, DING Haijun, et al. Characterization of reconsolidation volumetric strain of liquefied sand and validation by centrifuge model tests[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1838-1845. doi: 10.11779/CJGE201410011 MOHAMMADI S D, NIKOUDEL M R, RAHIMI H, et al. Application of the dynamic cone penetrometer (DCP) for determination of the engineering parameters of sandy soils[J]. Engineering Geology, 2008, 101(3): 195-203. -