沥青混合料变形特性

王 岚，陈 刚，邢永明，胡江三

（内蒙古工业大学土木工程学院，内蒙古呼和浩特010051）

沥青混合料变形特性

王 岚，陈 刚，邢永明，胡江三

（内蒙古工业大学土木工程学院，内蒙古呼和浩特010051）

采用数字散斑相关技术对橡胶粉（CR）改性沥青混合料和苯乙烯-丁二烯-苯乙烯三嵌段共聚物（SBS）改性沥青混合料进行变形场分析，利用Vic-3D图像处理软件对散斑图像进行处理，提取试件表面随着荷载变化的位移场和应变场.通过预切口半圆弯拉实验研究沥青混合料的断裂韧性.结果表明：水平位移变化速率和水平应变适用于沥青混合料裂纹性质及断裂特性的分析.在沥青混合料裂纹开裂过程中，开裂时的荷载、最大荷载与最大水平应变对应荷载并非同一荷载值.与SBS改性沥青混合料相比，长期温度老化和冻融循环前后，CR改性沥青混合料都具有更强的抗开裂能力，并且长期温度老化对2种改性沥青混合料影响较大.通过分析断裂韧性指标得到与DSCM实验结果相吻合的开裂特性，验证了DSCM对沥青混合料裂纹研究的适用性.

沥青混合料；数字散斑相关法；半圆弯拉实验；变形场；断裂韧度

沥青混合料开裂是沥青路面破坏的主要形式之一［1-3］.目前，各国学者对沥青混合料开裂行为的研究主要集中于单纯的数值模拟分析和宏观实验方法［4-5］.Mahmoud等［6］通过建立内聚力模型对沥青混合料裂纹的产生及裂缝扩展规律进行预测，并证实了内聚力模型对沥青混合料和沥青路面的适用性；通过扩展有限元对沥青混合料半圆试件建立内聚力模型，验证了模型对沥青混合料裂缝扩展规律预测的准确性，并得出拉应力是裂纹扩展主要原因的结论.由于数值模拟时存在将实验条件及材料参数理想化的假设，不能真实地反应沥青混合料破坏情况.Kim等［7］通过水煮法、水浸法等宏观实验对沥青混合料由于界面黏结性退化而产生的开裂破坏进行研究.但是，一方面，此方法包含的主观因素多，导致实验结果误差较大；另一方面，实验结果只能对最终破坏情况作出判断分析，不能对实验过程进行描述，因此，单纯宏观实验也不能对沥青混合料开裂进行全面分析.

Yue等［8-10］研究发现，从细观角度对沥青混合料的变形进行精确测试可以准确评价沥青混合料的开裂行为.数字散斑相关法（digital speckle correlation method，DSCM）作为一种全场、非接触的光学测量技术非常适用于具有各向异性的复合材料变形特性分析，因此，近年来，该方法被引入并应用于沥青混合料变形特性实验研究.DSCM在20世纪80年代由Peters等［11-12］提出，经过几十年的发展，现已被广泛地应用于固体力学变形研究.2002年，Kim等［13］将DSCM应用于沥青混合料.之后，Bjorn等［14-15］应用DSCM对沥青混合料应变场分析计算，验证了DSCM的优越性，并利用DSCM对沥青混合料裂纹扩展过程进行应变场分析，证明了DSCM对于沥青混合料这种具有异构性质的黏弹性材料有较高的适用性.

鉴于DSCM用于分析沥青混合料的优越性，本文在宏观半圆试件弯拉实验的基础上，结合DSCM对不同沥青混合料的断裂特性进行研究.通过预切口半圆试件弯拉实验对沥青混合料断裂韧性进行计算，并验证DSCM在沥青混合料研究中的适用性.

1 实验材料与方案

1.1 实验材料

1.1.1 沥青 本实验采用苯乙烯-丁二烯-苯乙烯三嵌段共聚物（styrene-butadiene-styrene，SBS）改性沥青和橡胶粉（crumb rubber，CR）改性沥青作为胶结材料.CR改性沥青是在90＃基质沥青中掺入质量分数为20%的30目胶粉颗粒经湿法加工制成，SBS改性沥青是在90＃基质沥青中掺入4%的SBS改性剂制成，相应技术指标如表1所示.PN为针入度（实验温度为5℃），tS为软化点，D为延度，tF为闪点.

表1 CR改性沥青和SBS改性沥青技术指标Tab.1 Properties of CR modified asphalt and SBS modified asphalt

1.1.2 沥青混合料 2种改性沥青混合料试件均采用AC-16级配，其中粗、细集料均采用玄武岩，矿料采用石灰岩矿粉.

1.1.3 半圆试件制作及实验过程 采用旋转压实仪成型沥青混合料圆柱形试件，并将圆柱形试件切割成一定尺寸的半圆试件.对半圆试件进行长期温度老化、水溶液以及盐溶液冻融循环处理.将半圆试件表面喷涂哑光白漆，与沥青混合料自身表面的黑色形成反差，从而在试件表面形成散斑（为了提高DSCM算法中图像子区域相关搜索的精确性，要求散斑试件的灰度值随机分布），并对试件进行半圆弯拉实验，利用Vic-3D软件对采集到的试件破坏过程图像进行分析计算.

1.2 实验方法及原理

1.2.1 半圆试件弯拉实验 半圆弯拉实验（semicircular bendin，SCB）试件的加载方式如图1所示，试件尺寸如下：直径2R＝100 mm，厚度B＝5 mm，预切口深度a分别为1、2和3 cm，支撑距离2D＝80 mm，加载速率v＝1 mm／min，实验温度为10℃.

图1 半圆弯拉试件加载示意图Fig.1 Schematic diagram of SCB test

1.2.2 DSCM的构成 DSCM是一种基于光学技术的图像处理方法，其测试系统如图2所示，主要由光源、CCD相机、图像采集系统及Vic-3D图像处理软件组成.同时启动试件加载系统与图像采集系统，通过CCD相机对加载过程中的试件散斑图像进行实时采集并储存于计算机中，将采集的散斑图像输入Vic-3D处理软件，对输入图像进行分析.

图2 DSCM测试系统Fig.2 Test system of DSCM

1.2.3 DSCM基本原理 DSCM是一种从物体表面随机分布的散斑直接提取位移场和应变场的测量方法.DSCM基本原理如下：给定物体变形前、后的2个散斑图像，在变形前、后的散斑场中识别出对应于变形前的散斑区域，判断物体变形后，散斑场经变换恢复到变形前对应的散斑场相似程度的指标，即散斑场之间的相关性，把变形测量问题转化为一个相关计算问题，如图3所示.图中，V为竖向位移，U为横向位移.将图中窗口区域的像素中心点P（x0，y0）定义为亚像素区域，通过预先设定的搜索方法及计算方法在变形图像中找到一个与参考图形中相似度最高的子区域P′（x0′，y 0′）.

考虑到图像采集时可能出现的曝光过度、光照不均等情况，采用零均值归一化的最小平方距离相关函数：

图3 图像变形前后计算区域的变形示意图Fig.3 Schematic figure of reference square subset before and after deformation

式中：Cf，g（¯P）为一个用来描述相似度的方程，f（x，y）为图3（a）中点P的灰度值，g（x′，y′）为图3（b）中点P′的灰度值，fm、gm分别为参考图像与变形图像中窗口区域灰度值，M为变形场内各个数据点的坐标值.

2 实验结果分析

2.1 沥青混合料变形与变形速率分析

图4 半圆试件开裂前后示意图Fig.4 Schematic diagram of semi-circle specimen before and after crack

如图4所示为试件开裂前后的示意图，其中图4（b）右侧标尺为图4（b）云图中水平应变的具体数值.图5为SBS改性沥青混合料半圆试件底边各点（如图4（a）所示）在开裂时位移、位移变化速率、应变以及应变率的关系曲线，图中x为位移，d x／d t为位移变化速率，dε／d t为应变，dε／d t为应变率，Pn为试件上所对应点.由图5（a）可知，试件开裂时底边各点的水平位移（U）在点P1～P21为负值（分析软件设定变形方向向右时为正方向），点P22～P28为正值，且点P22的水平位移有较大变化；竖向位移（V）由点P1～P21沿负方向（分析软件设定变形方向向上时为正方向）增大，由点P22～P28沿负向减小；离面位移（W）几乎相等.由图4（b）可知，试件的开裂点在点P22的位置，在竖向荷载作用下试件底部裂纹的主要型式为张开型（Ⅰ型裂纹），Ⅰ型裂纹在试件受水平拉力时产生，因此，水平位移对试件开裂影响较大.试件开裂后开裂处左、右位置点的水平位移会向不同的方向变化.竖向位移是由于试件在竖向荷载作用下产生竖向挠度.

图5 试件底边不同位置点开裂时变形及变形速率曲线Fig.5 Curves of deformation and deformation rate for different position points of specimen bottom when cracking

由图5（b）可知，水平位移变化速率d U／d t在开裂点（点P22）位置产生最大的突变，而竖向位移变化速率（d V／d t）和离面位移变化速率（d W／d t）相对d U／d t变化较小，且在开裂点位置没有突出变化.由图5（c）可知，水平应变（εU）较大值出现在点P11、P22处.由图4可知，点P22为试件开裂点，且点P11、P22均处于集料与沥青胶浆界面处.这是由于界面处的自身缺陷和较小模量而产生应力集中，从而产生较大应变，开裂处由于裂点的应力奇异场而具有更大应变值.竖直应变（εV）和离面应变（εW）在不同位置点几乎没有变化.由图5（d）可知，水平应变率（d x／d t）在点P22具有最大值，这是由于在开裂点位置具有最大水平应变的原因，而竖向应变率（dεV／d t）和离面应变率（dεW／d t）在开裂点与水平应变率相比变化较小.因此，水平变形以及变形速率对试件开裂最为敏感，这与试件主要受水平拉应力作用产生开裂相吻合.由于水平位移变化速率和水平应变在开裂点变化最为明显，本文主要研究试件的开裂破坏，选择水平位移变化速率和水平应变作为沥青混合料断裂性能指标进行研究.

2.2 荷载对水平应变的影响

如图6所示为SBS改性沥青混合料半圆试件底边开裂时全场水平应变云图.为了更好地区分沥青胶浆、集料与界面，取点P′1、P′2、P′3、P′4、P′5、P′6、P′7、P′8对应图4（a）中的点P1、P6、P11、P14、P17、P22、P24、P26进行研究.由图6可知，P′6具有全场最大水平应变，且位于P′3界面处的水平应变值也较其他各处大，这与图5（c）的分析结果吻合.如图7所示为试件底边点P′1～P′8处的水平应变随荷载的变化曲线，图中1.000、2.000、3.000、3.632以及4.078 k N为加载阶段对应的荷载.其中，3.632 k N为微裂缝产生时的荷载，4.078 k N为荷载峰值.为了将卸载阶段荷载与加载阶段荷载加以区分，用括号内数字表示卸载阶段荷载，如用（3.758）k N、（2.000）k N、（1.000）k N表示卸载阶段荷载，其中（3.758）k N为水平应变最大值对应的荷载.由图7可以看出，当荷载为1.000、2.000以及3.000 k N时，底边最大位移位于点P′3和P′5处（点P′3与P′5为沥青胶浆与粗集料界面）；当荷载增加至3.632 k N（微裂缝产生荷载）时，最大应变位于点P′6处，与图6应变云图相吻合，说明沥青混合料开裂时裂尖处具有全场最大εU.由图6可知，点P′3、P′5、P′6均为沥青胶浆与粗集料界面点，点P′3、P′5处集料尺寸大于点P′6，因此，点P′3、P′5处的集料与胶浆的接触面积大于点P′6处的接触面积.在复合材料力学理论中，受到外荷载作用时损伤总倾向于发生在界面结合部或其附近.这是由于结合材料界面附近不仅容易存在各种各样的缺陷，导致结合能力的降低，而且容易引发应力集中并产生残余应力，使界面处产生较高的应力水平，从而产生较大应变.点P′3、P′5处由于界面面积大，界面黏结力也较大，即对应的抵抗破坏能力比小界面要大，点P′6处小界面先被破坏而产生开裂.当裂纹形成时，由于裂纹尖端产生的应力奇异场，裂尖处具有最大应变.

图6 开裂时试件底边位置点及开裂云图Fig.6 Bottorn location points and cracking nephogram when cracking

图7 试件底边各点水平应变随荷载的变化曲线Fig.7 Changing curves of horizontal strain of specimen bottom’s points with loading

开裂点处水平应变的最大值并未出现在荷载峰值（4.078 k N）处，而其对应卸载阶段的荷载为（3.758）k N.因此，在沥青混合料开裂过程中，最大荷载与最大应变均未出现在开裂时刻，最大荷载出现在开裂后，而最大应变出现在最大荷载之后.这是由于沥青混合料开裂是特指微裂缝在底边形成时，而沥青混合料开裂过程分为微裂缝形成阶段、微裂缝汇集阶段与形成宏观裂缝阶段.在微裂缝形成时宏观裂缝并未产生，还需要继续施加荷载使微裂缝汇集贯通；当微裂缝贯通形成长裂缝时，试件达到破坏，此时对应的破坏荷载降低；宏观裂纹未形成时所有微裂纹尖端都处于应力集中状态，由于微裂缝多，应力较分散，应变也较分散，从而导致微裂缝形成时对应应变并非应变最大值；当宏观裂纹形成后，全场水平拉应力集中于裂纹尖端，此时裂尖处具有最大应力，因此具有最大拉应变.综上所述，微裂缝形成时对应荷载在荷载峰值之前，荷载峰值（4.078 k N）对应微裂缝汇集成长裂缝，最大水平应变对应宏观裂纹形成.

如图8所示为SBS改性沥青混合料半圆试件竖向裂缝形成时全场应变云图.为了研究试件竖向裂缝的性质，图8另取一裂缝方向接近竖直的试件为研究对象，沿裂缝取7个点对竖向裂缝进行研究.

图8 竖向裂缝位置点及应变云图Fig.8 Points and cloud picture of specimen for vertical crack

如图9所示为竖直方向7个点水平应变（εU）随荷载（L）变化关系曲线.由图9可知，开裂时点P″1（开裂点）的水平应变最大，这与图5分析结果一致，证明水平应变对研究试件裂缝形成适用性.图9中，1.000～6.000 k N为加载阶段荷载，（0.500）～（5.000）k N为卸载阶段荷载.由图9可知，εU值在5.000 k N时开始增长，最大应变对应荷载并非荷载峰值6.000 k N，而是点P″1、P″2、P″3、P″4在卸载阶段（5.000）k N时，εU达到峰值，点P″5、P″6、P″7在卸载阶段（4.000）k N时，εU达到峰值.εU值在荷载为5.000 k N时开始增长，此荷载为微裂缝产生荷载，这是由于微裂缝产生时的应力集中使应变开始增长；6.000 k N是微裂缝汇集为长裂纹荷载即加载峰值；（5.000）k N为点P″1、P″2、P″3、P″4处宏观裂缝形成荷载，（4.000）k N为点P″5、P″6、P″7处宏观裂缝形成荷载，这是由于点P″1、P″2、P″3、P″4处的微观裂缝先在同一荷载形成宏观裂缝，随裂纹扩展而使点P″5、P″6、P″7在另一荷载形成宏观裂缝，由此也可证明宏观裂纹是由于微观裂纹同时汇集而形成.以上分析与图7分析结果相吻合，可知在裂纹开裂过程中开裂时荷载、最大荷载与最大水平应变对应荷载并非同一荷载值，开裂时荷载在最大荷载之前，即试件微裂缝在最大荷载之前形成，最大水平应变对应荷载在最大荷载之后，即宏观裂缝在最大荷载之后形成.

图9 竖向裂缝各点水平应变随荷载的变化曲线Fig.9 Curves of horizontal strain at different points for vertical crack with loading

2.3 沥青混合料断裂性能分析

2.3.1 水平位移变化速率 如图10所示为在不同实验条件下CR与SBS改性沥青混合料对应的水平位移变化率d U／d t随加载时间变化曲线.由图5可知，d U／d t发生突变时间点为试件裂缝产生时刻.由图10可知，CR改性沥青混合料的裂缝产生d U／d t发生突变）时间点大约为280 s，经水溶液和盐溶液冻融循环的裂缝产生时间点大约都是210 s，经长期温度老化的裂缝产生的时间点大约为120 s，说明CR改性沥青混合料的开裂性能受长期温度老化的影响较受水溶液与盐溶液冻融循环的影响严重.由图10（b）可知，SBS改性沥青混合料的裂缝产生时间点大约为120 s，经水溶液冻融循环的裂缝产生时间点大约都是85 s，经盐溶液冻融循环的裂缝产生时间点大约都是75 s，经长期温度老化的裂缝产生时间点大约为65 s.由此可知，水溶液、盐溶液冻融循环以及温度老化都对SBS改性沥青混合料有较大影响，其中长期温度老化产生的影响最大.

对比分析图10（a）、（b）可知，CR改性沥青混合料及经冻融和温度老化的沥青混合料开裂时间均较对应的SBS改性沥青混合料晚，说明CR改性沥青混合料具有较好的抗开裂能力.由图10还可以看出，CR改性沥青混合料及经冻融和温度老化的沥青混合料开裂时d U／d t呈缓慢下降趋势，而SBS改性沥青混合料则呈迅速下降趋势，说明CR改性沥青混合料及经冻融和温度老化的沥青混合料在开裂时裂缝形成及扩展速率慢，基本呈塑性破坏，SBS改性沥青混合料及经冻融和温度老化的沥青混合料在开裂时裂缝形成及扩展速率快，基本呈脆性破坏.

2.3.2 不同实验条件下沥青混合料的断裂韧性分析 由于沥青混合料具有明显的弹塑性能，线弹性力学行为不适合描述沥青混合料的破坏行为，此时的力学行为应采用弹塑性断裂力学进行研究.弹塑性断裂力学的J积分是一个与线路无关的能量积分，是裂缝尖端附近弹塑性场强弱的一个平均值，无论是线弹性还是非线性的荷载-位移关系，都可采用J积分对沥青及改性沥青混合料的断裂性能进行描述.Little等［16］采用半圆弯拉实验研究沥青混合料的断裂性能，认为用断裂韧性JC作为沥青混合料的断裂参数是合适的，且断裂韧性越大其抗裂性能越好.

图10 不同沥青混合料水平位移变化速率随加载时间的变化曲线Fig.10 Curves of horizontal displacement rate with loading time for different asphalt mixtures

图11 不同实验条件下沥青混合料断裂韧性Fig.11 Fracture toughness of different asphalt mixtures in different experimental conditions

为了确定J积分值，采用不同预切口深度的半圆试件进行半圆弯拉（semi-circular bending，SCB）实验.SCB实验中断裂韧性的取得基于下式：

式中：E1、E2为2种不同锯缝深度试件的断裂应变能；a1、a2为2种锯缝深度；B1、B2为2种不同锯缝深度的试件厚度.

如图11所示为不同沥青混合料的断裂韧性关系图.由图可知，CR改性沥青混合料及经长期温度老化和冻融循环的改性沥青混合料的断裂韧性均较对应的SBS改性沥青混合料大，说明CR改性沥青混合料的抗裂性能优于SBS改性沥青混合料；并且，2种改性沥青混合料经长期温度老化后的断裂韧性值均为最小，由此可知长期温度老化对沥青混合料的抗裂性能的影响较水溶液与盐溶液冻融循环破坏的影响大.这与图10结果均吻合，也可验证DSCM在沥青混合料变形场研究中的适用性.

3 结 论

（1）通过水平变形及变形速率分析可知，采用水平位移变化速率和水平应变对沥青混合料裂纹性质及断裂特性进行分析研究.

（2）沥青混合料裂纹开裂过程中开裂时荷载、荷载峰值与最大水平应变对应荷载并非同一荷载值.

（3）CR改性沥青混合料及经长期温度老化和冻融循环的改性沥青混合料与其对应条件的SBS改性沥青混合料相比，具有较强的抗开裂能力.

（References）：

［1］DAS P K，JELIGEN D，BIRGISSON B.Evaluation of the low temperature cracking performance of asphalt mixtures utilizing HMA fracture mechanics［J］.Construction and Building Materials，2013，47：594- 600.

［2］MUNIANDY R，AKHIRN A B C M，HASSIM S，et al.Laboratory fatigue evaluation of modified and unmodified asphalt binders in stone mastic asphalt mixtures using a newly developed crack meander technique［J］.Journal of Materials in Civil Engineering，2014，59（3）：1- 8.

［3］ARAGO F T S，KIM Y R，LEE J，et al.Micromechanical model for heterogeneousasphalt concrete mixtures subjected to fracture failure［J］.Journal of Materials in Civil Engineering，2010，23（1）：30- 38.

［4］KIM Y R，BAEK C，UNDERWOOD B S，et al.Application of viscoelastic continuum damage model based finite element analysis to predict the fatigue performance of asphalt pavements［J］.KSCE Journal of Civil Engineering，2008，12（2）：109- 120.

［5］ZHU X，CHEN L.Numerical prediction of elastic modulus of asphalt concrete with imperfect bonding［J］.Construction and Building Materials，2012，35：45- 51.

［6］MAHMOUD E，SAADEH S，HAKIMELAHI H，et al.Extended finite-element modelling of asphalt mixtures fracture properties using the semi-circular bending test［J］.Road Materials and Pavement Design，2014，15（1）：153- 166.

［7］KIM Y R，PINTO I，PARK S W.Experimental evaluation of anti-stripping additives in bituminous mixtures through multiple scale laboratory test results［J］.Construction and Building Materials，2011，29（4）：386-393.

［8］YUE Z Q，CHEN S，THAM L G.Finite-element modeling of geomaterials using digital image processing［J］.Computers and Geotechnics，2003，30（5）：375- 397.

［9］SUTTON M A，WOLTERS W J，PETERS W H，et al.Determination of displacements using an improved digital correlation method［J］.Image and Vision Computing，1983，1（83）：133- 139.

［10］ZHANG X N，YUAN M M.Study on fatigue deformation of asphalt mixtures in digital speckle correlation method［J］.I nternational Journal of Digital Content Technology and its Applications，2012，6（17）：162- 170.

［11］PETERS W H，RANSON W F.Digital imaging techniques in experimental stress analysis［J］.Optical Engineering，1982，21（3）：427- 431.

［12］SUTTON M A，MCNEILL S R，JANG J.Effects of subpixel image restoration on digitalcorrelation error［J］.Journal of Optical Engineering，1988，27（10）：870- 877.

［13］KIM Y R，WEN H.Fracture energy from indirect tension test［J］.Journal of the Association of Asphalt Paving Technologists，2002，71：779- 793.

［14］BJORN B，ANTONIO M，ELENA R，et al.Characterisation of asphalt mixture cracking behaviour using the threepoint bending beam test［J］.International Journal of Pavement Engineering，2011，6（12）：569- 578.

［15］TAN Y Q，HOU M H，ZHANG L，et al.Studying the strain field distribution of asphalt mixture with the digitalspeckle correlation method［J］.Road Materials and Pavement Design，2014，15（1）：90- 101.

［16］LITTLE D N，MAHMOUB K.Engineering properties of first generation plasticized sulfur binders and low temperature fracture evaluation of plasticized sulfur paving mixtures［J］.Transportation Research Board，1985（1034）：103- 111.

Deformation characteristics of asphalt mixture

WANG Lan，CHEN Gang，XING Yong-ming，HU Jiang-san
（College of Civil Engineering，Inner Mongolia University of Technology，Hohhot 010051，China）

The deformation fields of crumb rubber（CR）modified asphalt mixture and styrene-butadienestyrene（SBS）modified asphalt mixture were analyzed and compared by using digital speckle correlation method（DSCM）.Vic-3D image manipulation software was applied to analyze speckle images.Deformation field and strain field of the samples’surface that changed along with load were extracted.Fracture toughness of the asphalt mixture was studied through pre-cut semi-circle flexural tensile test.Results show that the change rates of horizontal displacement and the horizontal strain are suitable for the analysis of the crack properties and fracture characteristics of asphalt mixture.The micro-crack initiation load，the maximum load and the maximum horizontal strain are not the same load value.Compared with the SBS modified asphalt mixture，the CR modified asphalt mixture has stronger ability to resist cracking before and after long-term aging and freeze-thaw cycle.The long-term temperature aging of asphalt mixture has greater impact on the two kinds of modified asphalt mixture.The cracking characteristics from analysis of fracture toughness index agree with the experimental results of DSCM，indicating that DSCM is suitable for research on crack in asphalt mixture.

asphalt mixture；digital speckle correlation method（DSCM）；semi-circular bending test；deformation field；fracture toughness

10.3785／j.issn.1008-973X.2015.09.026

U 443.33

A

1008- 973X（2015）09- 1805- 07

2014- 06- 05. 浙江大学学报（工学版）网址：www.journals.zju.edu.cn／eng

国家自然科学基金资助项目（11162012）.

王岚（1966－），女，教授，博导，从事道路工程材料研究.ORCID：0000-0002-8130-3702.E-mail：wanglan661018＠163.com

邢永明，男，教授，博导.ORCID：0000-0001-6726-4041.E-mail：xym＠imut.edu.cn