Mechanical Performance and Microstructure of Single Component Geopolymer Mortar
-
摘要: 为了研究不同NaOH浓度、胶砂比及溶胶比在多种固化温度下粉煤灰地聚物砂浆的强度发展规律及其微观机理,进行了力学性能试验,并用扫描电镜(SEM)和压汞试验(MIP)分析了其微观形貌、孔径分布. 分析结果表明: 对浓度为10%的NaOH溶液制备的地聚物砂浆试件,即使在很高的温度下固化也没有观察到明显的强度发展;随着NaOH浓度增加或固化温度上升,单组分地聚物砂浆的抗压和抗弯强度均可获得最佳值,该值出现位置由热固化温度和NaOH浓度的共同决定;地聚物的孔径分布微分曲线为单峰分布,控制NaOH的浓度可以大幅减小孔径微分曲线的峰值,显著降低地聚物的孔隙率;粉煤灰颗粒在NaOH的作用下逐渐溶解,并在其表面形成胶凝物质,当Na+ 浓度较低时,地聚物较少,通过控制NaOH浓度可以使得地聚物变得密实,提高其抗压强度;基于热力学关系的分形模型描述地聚物孔结构形态的效果最好,其次为孔轴线模型,空间填充模型和海绵体模型只能较好地描述胶凝孔隙和过渡孔隙的孔结构分形维数;基于热力学关系与基于海绵模型分形维数计算值均在2.0~3.0之间,与一般水泥基材料的结果相近;适当调整NaOH的浓度可以改善地聚物的孔隙结构.Abstract: In order to investigate the strength development rule and microscopic mechanism of fly ash based geopolymer mortar with different NaOH concentrations, cement-sand ratios and solution-cement ratios at various curing temperatures, the mechanical properties test were carried out. The microscopic morphology and pore size distribution were analyzed by scanning electron microscope (SEM) and mercury intrusion porosimetry (MIP), respectively. The results show as follows: For geopolymer mortar specimens prepared using NaOH solution with a concentration of 10%, no significant strength development was observed even at very high curing temperatures. With the increase of NaOH concentration or curing temperature, the compressive strength and flexural strength of single component geopolymer mortar can obtain their best value, and the location of the value is determined by both curing temperature and NaOH concentration. The pore size distribution differential curve of the geopolymer is a unimodal distribution. Controlling the concentration of can greatly reduce the peak value of the pore diameter differential curve and significantly reduce the porosity of the geopolymer. Fly ash particles dissolve gradually under the action of NaOH and form cementitious material on its surface. When the concentration of Na+ is low, the geopolymer polymerization is less; controlling the concentration of NaOH can make the geopolymer become dense and improve their compressive strength. The fractal model based on thermodynamic relation is the best one to describe the structure of the geopolymer pore, followed by the pore axis model, while the spatial filling model and the sponge model can only describe the fractal dimension of the pore structure of gel pores and transition pores. The calculated values of fractal dimension based on thermodynamic relationship and sponge model are between 2.0 and 3.0, which is similar to the results of general cement-based materials. Proper adjustment of the concentration of NaOH can improve the pore structure of the geopolymer.
-
道岔是机车车辆从一股轨道转入或越过另一股轨道的线路设备,是铁路轨道的重要组成部分,也是线路上的薄弱环节[1]. 轨距和轨底坡是高速道岔重要的设计技术条件,其参数的设置直接影响行车安全和行车品质. 60N钢轨在区间线路上的使用取得了成功,为提升列车过岔的平稳性和轮轨接触力学特性,在高速道岔区同样需要应用60N钢轨. 列车运行过程中,车轮始终处于动态磨损状态,对轮轨接触行为影响较大.
致谢:京沪高速铁路股份有限公司科技研究项目(京沪科研-2020-11).
针对线路上的轨距和轨底坡参数,国内外学者做了大量的研究. 杜星等[2]建立了LMD车轮和CHN60钢轨匹配的动力学模型,分析同一轨道在不同轨底坡条件下的动力学行为发现,轨底坡变化对列车直线运行时的平稳性、舒适性影响很大. 钱瑶等[3]对比分析了不同轨底坡下60N钢轨和高速车轮LMA、XP55、S1002G匹配时的轮轨接触行为,匹配较优的轨底坡是1/20和1/30. 陈嵘等[4-5]研究了我国地铁线路常用的LM型面与CHN60钢轨在不同轨距和非对称轨底坡下的轮轨接触特性,发现增大轨距和改变轨底坡可改善轮轨匹配关系. Cui等[6]提出一种优化轮轨廓形的“正向求解法”,车辆动力学行为结果表明,轨底坡1/30、轨距1435 mm是优化后廓形的最优轨道参数. Sánchez等[7]开发了一个严谨的测量轨距和轨底坡的程序确保测量过程的准确性. 李超等[8]分析了转辙器区采用动态轨距加宽技术的轮轨接触作用,该技术可减轻尖轨磨耗和滚动接触疲劳. Ye等[9]通过建立多体动力学模型,分析轨道参数对车轮磨耗的影响,发现轨道参数对车轮磨耗和脱轨安全性有较大影响. 闫正等[10]分析了高速动车组车轮踏面和高速60N钢轨道岔断面的静态接触特性,发现适当地增大轨距和轨底坡有利于改善轮轨接触状态. 上述研究表明,改变轨距和轨底坡参数对轮轨接触行为影响较大,且现行的1/40轨底坡和1435 mm轨距往往不是最优的轨道参数,而国内的研究主要集中在区间线路或是针对CHN60钢轨道岔,因此,有必要对新型350 km/h 60N钢轨18号高速道岔的合理轨距和轨底坡展开研究.
本文根据实测LMA磨耗车轮型面和60N钢轨高速道岔关键断面,基于迹线法原理和三维非赫兹滚动接触理论,建立道岔区轮轨滚动接触模型,计算不同轨距和轨底坡下的接触几何参数和静力学指标,并与CHN60钢轨高速道岔进行对比.
1. 道岔关键断面和磨耗车轮
1.1 60N钢轨关键断面
利用CAD导出各关键断面离散数据,将尖轨和基本轨分开,并线性插值,得到如图1所示60N钢轨18号高速道岔转辙器区钢轨模型,X为沿钢轨的纵向里程坐标,Y为钢轨横坐标,Z为钢轨竖向坐标,尖轨顶宽35 mm钢轨断面位于轮载过渡段,其接触行为较为复杂,轮轨相互作用剧烈,对研究道岔区的轮轨接触行为具有代表性,因此,本文选取尖轨顶宽35 mm关键断面进行计算. 60N与CHN60钢轨高速道岔在尖轨顶宽35 mm处钢轨廓形见图2.
1.2 LMA磨耗车轮
跟踪记录某线路上运行的CRH2型动车组上LMA车轮型面演变[11],从标准LMA磨耗车轮到运营里程15万km,再到25万km,磨耗量增大较为明显,本文取这三种磨耗车轮进行分析. 车轮凹型磨耗对轮对接触具有较大影响[12],随着列车运营里程的增加,凹型磨耗车轮的磨耗量大致呈线性增大,如图3所示. 对磨耗车轮廓形通过三次样条函数进行插值、拟合.
2. 计算原理
2.1 轮轨接触几何计算原理
利用二维迹线法原理,不考虑车轮摇头角,由最小距离法搜索轮轨接触点. 接触点位置的求解有两个等价几何条件:1) 轮轨接触点处轮轨垂直距离为0,非接触点轮轨垂直距离大于0;2) 轮轨接触点处轮轨的轮廓线具有相同的斜率. 本文采用条件1求解轮轨接触点,再用条件2对结果加以验证.
2.2 UIC 519标准计算等效锥度
EN 15302标准[13]中,等效锥度的计算采用轮对周期运动的假设,这个方法称之为UIC 519标准[14]等效锥度. 自由轮对在轨道上的运动用微分方程表示为
¨y+v2er0Δr=0, (1) 式中:
y 为轮对在钢轨上的横向位移;v 为车辆前进速度;e 为接触点跨距;r0 为名义滚动圆半径;Δr 为右侧和左侧车轮滚动半径差.最后应用Klingel公式计算等效锥度
γe :γe=2er0(π λ)2 (2) 式中:
λ 为轮对在钢轨上做正弦运动的波长,其求解见文献[14].2.3 轮轨滚动接触力学原理
本文在接触力学部分采用Kalker的三维非赫兹滚动接触理论及其数值程序CONTACT,这是目前为止最为完善的滚动接触理论[15]. 该理论将轮轨接触问题转化为数学规划问题,利用Bossinesq-Cerruti公式可得轮轨滚动接触离散模型为:
{minCpJj=12pIiAIiJjpJj+[(g0J−q)pJz+(WJτ−uJτ)pJτ]A0,pJz⩾0,|pJτ|⩽bJ, ∀x∈Ac, (3) 式中:
C 为接触问题的余能;I、J 为接触斑上离散单元的编号;i、j、z 为X、Y、Z 坐标轴的方向;τ 为切向,对应XOY 平面;pIi 为单元I 沿方向i 的接触应力;pJj、pJz、pJτ 分别为单元J上沿Y轴、Z轴及方向τ 的接触应力;AIiJj 为影响函数,表示单元J方向j上单位力引起的单元I方向i上位移;g0J 为单元J中心处轮轨变形前接触面间的法向间隙;q为渗透量;WJτ 为单元J中心处从前时刻t′ 到当前时刻t总的刚性滑动量;uJτ 为轮轨滚动接触过程中时刻t时单元J中心处的弹性位移差;A0 为单元J的面积;bJ 为是单元J中心处Coulomb极限摩擦力;x为轮轨在X轴上坐标;Ac 为接触斑内坐标集合;车轮和尖轨的接触为异型接触,容易产生接触疲劳现象. 本文采用基于安定图模型的表面滚动接触疲劳因子预测不同工况下的轮轨表面滚动接触疲劳伤损. 轮轨接触斑中任一点的轮轨表面滚动接触疲劳系数
fI(•) 定义为{fI(x,y)=ft−kpz(x,y),ft=√(px(x,y))2+(py(x,y))2pz(x,y), (4) 式中:
ft 为牵引系数;y、z 分别为接触斑在Y、Z 轴上坐标值;px(•) 、py(•) 、pz(•) 分别为轮轨纵向、横向和法向接触应力;k 为材料自身的剪切屈服强度,取350 MPa.式(4)是根据赫兹接触理论得到,对于非赫兹接触问题,依据条带法,认为条带中间位置符合赫兹假设条件,从而将式(4)的应用扩展至非赫兹接触的范围,求得接触斑内任一单元的轮轨表面滚动接触疲劳系数,轮轨表面滚动接触疲劳因子定义为接触斑内滚动接触疲劳系数的最大值.
3. 计算结果与分析
3.1 轮轨接触点
轮轨几何接触是解释轮轨接触关系的基础[4]. 利用迹线法原理计算轮轨接触点,计算参数有:轮背距1353 mm,名义滚动圆半径460 mm,轮背到名义滚动圆处水平距离为70 mm,以Y轴正向为正,轮对横移量取−12 ~ 12 mm,横移量步长取0.5 mm. 以往研究表明,加宽轨距和改变轨底坡往往能够改善轮轨匹配关系[3-5],因此,对轨距1433、1435、1437、1439 mm,轨底坡1/10、1/20、1/30、1/40、1/50进行计算分析.
不同轨距和轨底坡参数下的接触点分布如图4和图5所示. 由图4可见:随着轨距的变化,轮轨接触点分布存在明显的差异;随着轨距的增加,轮轨接触点更多地分布在尖轨顶部,不易发生轮缘接触;随着轨距的增大,轮轨接触点由基本轨转移到尖轨所需要的横移量增大,也即轮载过渡延后,有利于减小尖轨受力,但会增大轮载过渡时轮轨接触点的跳跃,横向不平顺增大;在相同轨距下,随着车轮的磨耗,轮轨接触点更多地分布在基本轨,轮载过渡位置延后,横向不平顺增大.
由图5可见:不同轨底坡条件下,轮轨接触点的分布存在较大差异. 当轨底坡为1/10、1/20时,发生轮缘接触需要的轮对横移量最大,但和其余轨底坡相差较小;当轨底坡为1/10、1/20时,轮载过渡延后,横向不平顺增大;不同轨底坡下,发生轮载过渡时所需要的横移量最大为1/10,其次为1/20和1/30,最后为1/40和1/50.
3.2 等效锥度
等效锥度作为轮轨接触线性化指标,被广泛用于表征轮轨接触几何特征. 计算等效锥度的方法有简化法、谐波法和UIC 519法,UIC 519法采用轮对周期运动的假设,计算更准确[16],本文通过UIC 519标准[14]计算等效锥度.
在不同轨距下,磨耗车轮和60N关键断面匹配时的等效锥度见图6. 增大轨距有利于减小车轮踏面的等效锥度,从而提升列车过岔的平稳性;当轨距为1439 mm时,等效锥度基本上小于0.05,横移量相同条件下,其等效锥度是轨距1435 mm时的20%左右,极大减小了车轮等效锥度,提升了列车过岔平稳性;轨距变化对运营里程为25万km车轮的等效锥度影响较小.
图7为不同轨底坡和不同磨耗车轮条件下的等效锥度. 对于标准车轮和运营里程25万km车轮,1/30、1/40、1/50轨底坡条件下的等效锥度相差较小;当轨底坡为1/10、1/20,在横移量小于6 mm时,和轨底坡为1/30、1/40、1/50相比,等效锥度普遍较大,车辆过岔平稳性较差,横移量大于8 mm时,结果相反;在车轮运营里程为15万km时,1/30轨底坡条件下的等效锥度较小,列车过岔平稳性较好.
3.3 接触斑面积
相同法向轮轨力作用下,轮轨接触斑面积越大,其接触应力越小. 轮轨接触应力是影响轮轨磨耗和接触疲劳的重要因素. 利用非赫兹接触理论,计算轮轨接触斑面积和滚动接触疲劳因子,分析轨道参数取值对轮轨静力学接触行为的影响. 不同轨道参数工况下的参数取值:轴重14 t,单侧车轮轮心施加一半轴重,剪切模量82 GPa,泊松比0.28,摩擦系数0.3,划分网格单元0.2 mm × 0.2 mm,轮对横移量取0 ~ 12 mm. 选取右轮轨作为分析对象.
不同轨距和不同磨耗车轮条件下的接触斑面积如图8所示. 当横移量小于9 mm时,轨距和接触斑面积参数呈负相关,而当轮对横移量大于9 mm时,轨距和接触斑面积参数大致呈正相关;轮对横移量小于8 mm时,不同轨距条件下的轮轨接触斑面积相差较小,但当轮对横移量大于8 mm时,轨距越大,轮轨接触斑面积普遍越大;在轮对横移量大于8 mm时,轨距加宽有利于增大轮轨接触斑面积,减小轮轨接触应力.
不同轨底坡和不同磨耗车轮条件下的接触斑面积如图9所示. 1/10轨底坡下的轮轨接触斑面积普遍较小,且随车轮磨耗量的增大,其接触斑面积远小于其余轨底坡条件下的,说明1/10轨底坡下的轮轨接触力学性能较差,且会随着车轮磨耗变得更差. 由标准轮轨条件下,轨底坡和接触斑面积大致呈负相关. 随车轮磨耗量增大,1/30 ~ 1/50轨底坡的接触斑面积相差较小.
3.4 滚动接触疲劳因子
利用式(4)计算轮轨滚动接触疲劳因子,分析轮轨在不同轨道参数下的接触疲劳现象. 不同轨道参数下的表面滚动接触疲劳因子分布如图10、11所示. 由图10可知:在轮对横移量小于7 mm时,不同轨距下的表面滚动接触疲劳因子相差较小,在横移量大于 7 mm时,增大轨距可延缓轮轨表面进入滚动接触疲劳区;车轮运营里程达到15万km,在横移量为3 ~ 7 mm时,易发生轮轨间的两点接触,导致轮轨表面材料易进入疲劳区,从而产生疲劳破坏;增大轨距有利于减少轮轨材料出现接触疲劳现象,延长轮轨服役寿命.
由图11可知:标准轮轨条件下,轨底坡和滚动接触疲劳因子大致呈正相关;车轮磨耗导致轮轨滚动接触疲劳因子减小,原因是轮轨过渡位置延后,有利于减少尖轨磨耗,在横移量大于8 mm时,1/10和1/20轨底坡能延缓轮轨材料进入滚动接触疲劳区,但在横移量较小时,其滚动接触疲劳因子明显较大;1/10和1/20轨底坡下的轮轨滚动接触疲劳因子普遍较大,易引起轮轨材料进入滚动接触疲劳区,降低轮轨材料的使用寿命,轨底坡为1/30、1/40、1/50时,滚动接触疲劳因子相差较小.
随着车轮磨耗加深,道岔区的轮轨接触行为变得更为复杂. 运营里程15万km车轮和60N钢轨高速道岔在轨距1439 mm条件下发生轮轨两点接触时的接触斑分布如图12所示,接触斑中箭头指向表示切向应力合力的方向,箭头长短代表合力大小,滑动区轮轨表面间发生了相对滑移,黏着区轮轨间存在滑动趋势,但没有相对滑移. 由图12(a)可见,轨距1439 mm,轮对横移量5 ~ 7 mm时,发生了轮载过渡行为,随着横移量的增大,轮载更多由尖轨承载. 由图12(b)可见:在轮对横移量5 ~ 6 mm时,尖轨上存在较大的滑动区,在轮对横移量7 mm时,基本轨上存在较大的滑动区,两点接触导致了较大滑动区的存在,易导致车轮在钢轨上空转,引起道岔钢轨的磨损,从而缩短道岔区钢轨服役寿命. 需要指出的是,法向接触应力较小时同样可能存在较大滑动区,如图12中轮对横移量为5、7 mm时所示.
3.5 60N和CHN60钢轨高速道岔轮轨接触几何对比
为对比60N和CHN60钢轨高速道岔的轮轨接触几何行为,在轨距1439 mm和轨底坡1/30、1/40、1/50条件下,计算60N、CHN60钢轨高速道岔35 mm顶宽关键断面与三种磨耗车轮的等效锥度,如图13.
由图13可知:相同轨底坡和轮对横移量条件下,相比于CHN60钢轨,60N钢轨高速道岔区的等效锥度更小,列车过岔平稳性更优;车轮运营里程达到25万km时,当车轮横移量小于5 mm时,60N钢轨高速道岔区等效锥度普遍更小,具有较好的过岔平稳性,车轮横移量大于5 mm时,60N钢轨高速道岔区等效锥度普遍更大,轮对的对中性能更好;随车轮磨耗量增加,等效锥度整体呈增加趋势.
4. 结 论
1) 轨距对轮轨接触行为影响较大. 轨距加宽有利于减少轮缘接触,较大程度减小等效锥度,提升列车过岔的平稳性;轨距加宽可减小轮对横移量大于8 mm时的轮轨接触应力和表面滚动接触疲劳因子,减少轮轨材料发生接触疲劳,延长尖轨使用寿命. 轮对横移量小于8 mm时,轨距加宽对运营里程25万km车轮的接触性能影响较小.
2) 轨底坡对轮轨接触行为影响较大. 标准轮轨条件下,轨底坡和接触斑面积大致呈负相关,与滚动接触疲劳因子大致呈正相关. 轨底坡为1/10和1/20时,轮载过渡位置延后,横向不平顺增大,车轮横移量小于6 mm时,等效锥度普遍较大;1/10和1/20轨底坡下的接触斑面积普遍较小,轮轨滚动接触疲劳因子普遍较大,较易引起轮轨材料进入滚动接触疲劳区,降低轮轨材料的使用寿命,且1/10轨底坡对车轮磨耗的适应性较差. 轨底坡为1/30、1/40、1/50时,轮轨接触参数相差较小,匹配性能较优.
3) 和CHN60钢轨高速道岔相比,60N钢轨的等效锥度普遍更小,列车过岔平稳性更优;车轮运营里程为25万km时,当轮对横移量小于5 mm时,60N钢轨的等效锥度普遍更小,当轮对横移量大于5 mm时,结果相反.
4) 车轮磨耗易引起道岔区轮轨间的两点接触,在较小轮轨法向接触应力下,接触斑上也易出现较大滑动区,导致车轮空转,引起钢轨伤损. 随车轮磨耗量增加,轮轨间等效锥度整体呈增加趋势.
致谢:京沪高速铁路股份有限公司科技研究项目(京沪科研-2020-11)的支持.
-
-
XIE J, WANG J, ZHANG B, et al. Physicochemical properties of alkali activated GGBS and fly ash geopolymeric recycled concrete[J]. Construction and Building Materials, 2019, 204: 384-398. doi: 10.1016/j.conbuildmat.2019.01.191 CRETESCU I, HARJA M, TEODOSIU C, et al. Synthesis and characterisation of a binder cement replacement based on alkali activation of fly ash waste[J]. Process Safety and Environmental Protection, 2018, 119: 23-35. doi: 10.1016/j.psep.2018.07.011 WIANGLOR K, SINTHUPINYO S, PIYAWORAPAIBOON M, et al. Effect of alkali-activated metakaolin cement on compressive strength of mortars[J]. Applied Clay Science, 2017, 141: 272-279. doi: 10.1016/j.clay.2017.01.025 PART W K, RAMLI M, CHEAH C B. An overview on the influence of various factors on the properties of geopolymer concrete derived from industrial by-products[J]. Construction and Building Materials, 2015, 77: 370-395. doi: 10.1016/j.conbuildmat.2014.12.065 BAKHAREV T. Geopolymeric materials prepared using class F fly ash and elevated temperature curing[J]. Cement and Concrete Research, 2005, 35(6): 1224-1232. doi: 10.1016/j.cemconres.2004.06.031 RATTANASAK U, CHINDAPRASIRT P. Influence of NaOH solution on the synthesis of fly ash geopolymer[J]. Minerals Engineering, 2009, 22(12): 1073-1078. doi: 10.1016/j.mineng.2009.03.022 PANIAS D, GIANNOPOULOU I P, PERRAKI T. Effect of synthesis parameters on the mechanical properties of fly ash-based geopolymers[J]. Colloids and Surfaces A:Physicochemical and Engineering Aspects, 2007, 301(1): 246-254. 龙涛, 石宵爽, 王清远, 等. 粉煤灰基地聚物再生混凝土的力学性能和微观结构[J]. 四川大学学报(工程科学版), 2013, 45(增刊1): 43-47.LONG Tao, SHI Xiaoshuang, WANG Qingyuan, et al. Mechanical properties and microstructure of fly ash based geopolymeric polymer recycled concrete[J]. Journal of Sichuan University (Engineering Science Edition), 2013, 45(S1): 43-47. PAN Z, TAO Z, CAO Y F, et al. Compressive strength and microstructure of alkali-activated fly ash/slag binders at high temperature[J]. Cement and Concrete Composites, 2018, 86: 9-18. doi: 10.1016/j.cemconcomp.2017.09.011 HAMIDI R M, MAN Z, AZIZLI K A. Concentration of NaOH and the Effect on the properties of fly ash based Geopolymer[J]. Procedia Engineering, 2016, 148: 189-193. doi: 10.1016/j.proeng.2016.06.568 ATIŞ C D, GÖRÜR E B, Karahan O, et al. Very high strength (120 MPa)class F fly ash geopolymer mortar activated at different NaOH amount,heat curing temperature and heat curing duration[J]. Construction and Building Materials, 2015, 96: 673-678. doi: 10.1016/j.conbuildmat.2015.08.089 SOMNA K, JATURAPITAKKUL C, KAJITVICHYANUKUL P, et al. NaOH-activated ground fly ash geopolymer cured at ambient temperature[J]. Fuel, 2011, 90(6): 2118-2124. doi: 10.1016/j.fuel.2011.01.018 李星烨. 碱激发粉煤灰砂浆性能研究[D]. 成都: 西南交通大学, 2018 CHINDAPRASIRT P, JATURAPITAKKUL C, CHALEE W, et al. Comparative study on the characteristics of fly ash and bottom ash geopolymers[J]. Waste Management, 2009, 29(2): 539-543. doi: 10.1016/j.wasman.2008.06.023 LEE W K W, DEVENTER J S. The effects of inorganic salt contamination on the strength and durability of geopolymers[J]. Colloids and Surfaces A:Physicochemical and Engineering Aspects, 2002, 211(2): 115-126. YANG W C, GE Y, YUAN J, et al. Pore fractal characteristic of cement pastes with inorganic salts[C]//International Conference of Concrete Pavement. Haikou: [s.n.], 2009: 101-115. JI X, CHAN S Y, FENG N. Fractal model for simulating the space-filling process of cement hydrates and fractal dimensions of pore structure of cement-based materials[J]. Cement and Concrete Research, 1997, 27(11): 1691-1699. doi: 10.1016/S0008-8846(97)00157-9 金珊珊,张金喜,李爽. 混凝土孔结构分形特征的研究现状与进展[J]. 混凝土,2009(10): 34-37,42.JIN Shanshan, ZHANG Jinxi, LI Shuang. Current situation and development of fractal characteristic of pore structure of concrete[J]. Concrete, 2009(10): 34-37,42. ZHANG B, LIU W, LIU X. Scale-dependent nature of the surface fractal dimension for bi-and multi-disperse porous solids by mercury porosimetry[J]. Applied Surface Science, 2006, 253(3): 1349-1355. doi: 10.1016/j.apsusc.2006.02.009 -