CN114441436B - 一种基于路面纹理评价抗滑性能的分析方法 - Google Patents
一种基于路面纹理评价抗滑性能的分析方法 Download PDFInfo
- Publication number
- CN114441436B CN114441436B CN202111564062.1A CN202111564062A CN114441436B CN 114441436 B CN114441436 B CN 114441436B CN 202111564062 A CN202111564062 A CN 202111564062A CN 114441436 B CN114441436 B CN 114441436B
- Authority
- CN
- China
- Prior art keywords
- texture
- test piece
- elevation
- section
- piece plate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 16
- 238000012360 testing method Methods 0.000 claims abstract description 132
- 238000005096 rolling process Methods 0.000 claims abstract description 47
- 239000010426 asphalt Substances 0.000 claims abstract description 31
- 239000000203 mixture Substances 0.000 claims abstract description 29
- 238000004321 preservation Methods 0.000 claims abstract description 11
- 239000004568 cement Substances 0.000 claims abstract description 4
- 238000005520 cutting process Methods 0.000 claims abstract description 4
- 238000004519 manufacturing process Methods 0.000 claims abstract description 4
- 238000010276 construction Methods 0.000 claims description 45
- 238000000034 method Methods 0.000 claims description 35
- 238000011156 evaluation Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 13
- 230000003068 static effect Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 6
- 238000005056 compaction Methods 0.000 claims description 5
- 230000001419 dependent effect Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000011056 performance test Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 10
- 239000000463 material Substances 0.000 description 7
- 238000011160 research Methods 0.000 description 5
- 238000010998 test method Methods 0.000 description 5
- 206010039203 Road traffic accident Diseases 0.000 description 4
- 239000011384 asphalt concrete Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012417 linear regression Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000007906 compression Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000465 moulding Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000004567 concrete Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000013521 mastic Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000010972 statistical evaluation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N19/00—Investigating materials by mechanical methods
- G01N19/02—Measuring coefficient of friction between materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
- G01B11/0608—Height gauges
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/30—Adapting or protecting infrastructure or their operation in transportation, e.g. on roads, waterways or railways
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Automation & Control Theory (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Road Repair (AREA)
Abstract
本发明涉及一种基于路面纹理评价抗滑性能的分析方法，包括步骤一，制作试件板，步骤二，根据T0702‑2011选择沥青混合料的拌和温度和碾压温度进行拌和与碾压；步骤三，在水泥地面上按照试件板的试模尺寸切割成型槽，以放置试模，将保温后的沥青混合料放入试模内成型试件板；步骤四，利用激光纹理扫描仪对试件板上的测点表面进行扫描，以获取路面表面纹理的高程数据。本发明建立了基于纹理的抗滑性能预测模型，可以实现基于路面纹理的抗滑性能精准预测，并且通过激光技术方便快捷获取路面纹理，避免因抗滑性能测试而阻断交通或威胁测试人员安全。
Description
技术领域
本发明涉及路面抗滑性能分析技术领域，特别涉及一种基于路面纹理评价抗滑性能的分析方法。
背景技术
路面有充足的抗滑性能对公路的行车安全至关重要，但随着时间推移，沥青路面的沥青薄膜脱落和磨损，导致集料暴露在外，在行车荷载和环境等的影响下，集料发生磨光导致路面抗滑性能下降。道路交通事故的发生，是由驾驶员、环境和道路条件共同决定的，研究表明路面抗滑性能下降与交通事故的发生有一定的联系。对路面抗滑性能及时监测，掌握抗滑性能的衰减状况，并对抗滑性能不足的路段及时进行养护，可以降低因路面抗滑性能下降而导致的交通事故率。但相关研究的前提和基础是对路面抗滑性能现状的评价和评判。
目前抗滑性能常用的评价方法可分为两种：直接测试和间接测试。常用的直接测试方法主要包括摆值法、动态旋转式摩擦系数测试法、横向力系数测试车法等，其中摆值法和动态旋转式摩擦系数法均可适用于室内外测试，但前者由于设备性价比高，应用更为广泛，后者价格昂贵，仅部分单位有购置；常用间接测试法主要包含传统的铺沙法和激光纹理扫描法，前者从宏观纹理角度间接地对路面抗滑性能进行评价，在国内外被广泛应用，但测量结果受操作人员的影响较大，测试方法比较粗糙，测量结果有一定的离散性，且测得的构造深度也未能全面、精确地表征路面的纹理的信息，后者具有测试速度快、精度高、操作便捷等优点，能够实现路面表面纹理的高效、准确采集，且随着激光扫描技术的迅速发展，开始逐渐利用更为精密的试验设备对路面表面纹理进行识别和分析。
路面纹理可以分为微观纹理(波长可达0.5mm)、宏观纹理(波长0.5mm～50mm)、大构造(波长50mm～500mm)和不平整构造(波长大于500mm)，其中宏观和微观纹理是抗滑性能的关注对象。目前针对路面纹理评价的研究，主要侧重宏观纹理，对路面微观纹理关注相对较少。实际上，作为轮胎与路面直接接触的部位，路面的宏观纹理和微观纹理对轮胎-路面间的抗滑性能均有影响。此外，宏观纹理的波段范围很广，深入研究其对抗滑性能影响时，有待继续细分。
针对上述描述内容，本专利采用国内外广泛应用的抗滑性能性能评价方法摆值法进行抗滑性能测试，采用方便、快速、准确的激光技术并结合滤波法获得宏观纹理分波段评价指标、微观纹理和常规纹理评价指标，建立抗滑性能和纹理的关联，为快速评价抗滑性能提供了一种方法，避免了抗滑性能测试对现行交通的影响及对测试人员安全的威胁。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的一个目的在于提出一种基于路面纹理评价抗滑性能的分析方法，以解决背景技术中所提到的问题，克服现有技术中存在的不足。
为了实现上述目的，本发明提供一种基于路面纹理评价抗滑性能的分析方法，包括以下步骤，
步骤一，按照T0702-2011成型标准马歇尔试件，以计算出标准试件对应的标准密度，将标准密度作为试件板的密度，选取试件板的尺寸并计算其体积，将体积和标准密度相乘并乘以系数1.03最终求得试件板所需的沥青混合料的质量，按照试件板的尺寸制作试模；
步骤二，根据T0702-2011选择沥青混合料的拌和温度，根据拌和温度将沥青混合料进行拌和，将拌和后的沥青混合料放置在恒温箱中并按照保温温度和保温时间进行保温；
步骤三，在水泥地面上按照试件板的试模尺寸切割成型槽，以放置试模，将保温后的沥青混合料放入试模内，并通过手推式小型压路机按照设定碾压参数对沥青混合料进行碾压，最终获得试件板；
步骤四，利用激光纹理扫描仪对试件板上的测点表面进行扫描，以获取路面表面纹理的高程数据；
步骤五，对步骤四中获取的高程数据按照公式一进行坏点处理，以将坏点去除；按照公式二对高程数据进行倾斜度和偏置误差处理；
Zi1＝min(zn,zm) (1)
其中，i——为路面断面构造线中坏点的编号；
m——断面构造线中距离坏点最近的前一个点的编号；
n——断面构造线中距离坏点最近的后一个点的编号；
zm——为编号为m的断面构造线的高程；
zn——为编号为n的断面构造线的高程；
Zi1——为坏点处理后的断面高程；
Zi＝zi-b1i-b0 (2)
其中，b1、b0——分别为倾斜度和偏置误差的消除系数；
i——为断面构造线中采集点的水平编号；
N——为断面构造线在评价长度范围内的采集点总数；
zi——为倾斜度和偏置误差消除前i编号处的断面高程；
Zi——为倾斜度和偏置误差消除后i编号处的断面高程值；
步骤六，将步骤五得到的高程数据进行平均断面深度MPD、均方根Rms、偏斜度Rsk以及不同波段的高程方差的计算，其中不同波段的高程方差包括宏观纹理长波段方差L、宏观纹理短波段方差S、微观纹理波段方差M；
步骤七，以试件板上测点的中心为控制点，进行BPN测试，得到摆值BPN；
步骤八，以BPN为因变量，以MPD、Rms、Rsk、L、S和M共六个参数为自变量，选用筛选参数较佳的Stepwise逐步回归法，以0.05为进入概率，以0.10为剔除标准，进行参数的选择，通过上述的多元线性回归方法建立了路面抗滑性能BPN与路面纹理的关系式如下所示，
BPN＝21.828*S-8.041*MPD+59.461 (3)
其中，BPN——摆值；
S——宏观纹理短波段方差；
MPD——测点的平均断面深度。
由上述任一方案优选的是，步骤二中，拌和温度为155～175℃，保温时间为1h，保温温度比碾压温度高5℃。
由上述任一方案优选的是，步骤三中的碾压参数包括碾压往返次数、碾压压力值和碾压温度，所述碾压往返次数为24次，按照碾压顺序依次为静压8次，振动压实8次、静压8次，碾压温度为150～170℃，碾压压力值为15kN。
由上述任一方案优选的是，在碾压之前，在沥青混合料的表面铺设厚度为5mm的橡胶块。
由上述任一方案优选的是，步骤四中，采用激光纹理扫描仪依次对试件板中间轮迹带上沿行车方向依次扫描的三个位置形成三个测点，每个测点的扫描面积为104mm×72.01mm，扫描每个测点均能够得到163条长度为104mm的均匀分布的断面构造线，相邻测点间距为0.009525mm，相邻两条断面构造线的间距为0.4445mm。
由上述任一方案优选的是，MPD的计算参照ISO-13473-1标准，每个测点的MPD是该测点163条断面构造线的断面构造深度的均值，试件板的MPD是三个测点的MPD值的平均值。
由上述任一方案优选的是，采用Butterworth高通和低通滤波器将三个测点的断面构造线滤波成三个不同的波段，分别为波段为5mm～50mm的宏观纹理长波段，波段为0.5mm～5mm的宏观纹理短波段，波段为0.024mm～0.5mm的微观纹理波段。
由上述任一方案优选的是，所述均方根Rms的计算公式如下，
其中，N——断面构造线采集点总数；
Zi——经过倾斜度和偏置误差处理的断面高程值。
由上述任一方案优选的是，所述偏斜度Rsk的计算公式如下，
其中，Rms——均方根；
N——断面构造线采集点总数；
Zi——经过倾斜度和偏置误差处理的断面高程值。
由上述任一方案优选的是，所述不同波段的高程方差计算公式如下，
其中，Zi——经特定频段滤波之后的断面高程；
——经过倾斜度和偏置误差处理的断面高程值均值。
与现有技术相比，本发明所具有的优点和有益效果为：
1、已有研究对路面宏观纹理关注较多，微观纹理关注相对较少，且路面宏观纹理包含的纹理波段范围较广(0.5mm～50mm)，不便于研究路面纹理对路面抗滑性能的影响。本申请采用Butterworth高通和低通滤波器将路面断面构造线滤波成三个不同的波段，5mm～50mm的宏观纹理长波段，0.5mm～5mm的宏观纹理短波段，0.024mm～0.5mm的微观纹理波段，从宏观纹理分段和微观纹理对路面纹理进行评价，提升了路面纹理评价的精确性。
2、摆式仪虽然性价比比较高，但测试过程中需要清洁路面、仪器调平、指针调零、校核滑块长度、洒水等细节，操作过程繁琐，测试效率较低，且测试结果受到操作人员影响较大。本申请通过建立抗滑性能BPN与路面表面纹理的关联，基于路面纹理可以通过激光技术获得，且具有测试速度快、精度高、操作便捷等优点，能够实现路面表面纹理的高效、准确采集，且随着激光扫描技术的迅速发展，更为精密的激光试验设备逐渐应用道路领域，为实现路面抗滑性能评价提供一种新方法。
3、本申请基于36组抗滑性能BPN数据和纹理数据，通过Stepwise方法筛选进入模型的自变量参数，最终建立的回归方程是一个基于纹理参数的最优自变量组合，实现了从路面纹理角度预测抗滑性能，促进了抗滑性能监测和衰变、路面及时养护等工作的开展，有利于降低抗滑性能不足导致的交通事故率。
4、本申请基于路面纹理建立的路面抗滑性能BPN预测模型，可以实现从路面纹理角度预测抗滑性能。在高精度车载激光的辅助下，能实现路面纹理的高效、快速、精确采集，避免了抗滑性能测试封闭路段对已有交通的影响，同时也有效避免了车流对抗滑性能测试人员安全的潜在威胁。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为试件板碾压方向和纹理测点示意图；
图2a为坏点处理前不同路面断面构造线示意图；
图2b为坏点处理后不同路面断面构造线示意图；
图3a为倾斜度和偏置误差处理前的不同路面断面构造线示意图；
图3b为倾斜度和偏置误差处理后的不同路面断面构造线示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
如图1-3所示，本发明实施例的一种基于路面纹理评价抗滑性能的分析方法，包括以下步骤，
步骤一，按照T0702-2011成型标准马歇尔试件，以计算出标准试件对应的标准密度，将标准密度作为试件板的密度，选取试件板的尺寸并计算其体积，将体积和标准密度相乘并乘以系数1.03最终求得试件板所需的沥青混合料的质量，按照试件板的尺寸制作试模；
步骤二，根据T0702-2011选择沥青混合料的拌和温度，根据拌和温度将沥青混合料进行拌和，将拌和后的沥青混合料放置在恒温箱中并按照保温温度和保温时间进行保温；考虑试件板成型采用的沥青种类不同，根据T0702-2011选择沥青混合料板拌和，拌和温度为155～175℃，保温时间为1h，保温温度比碾压温度高5℃；
步骤三，在水泥地面上按照试件板的试模尺寸切割成型槽，以放置试模，将保温后的沥青混合料放入试模内，装料完成后，在沥青混合料的表面铺设厚度为5mm的橡胶块，并通过手推式小型压路机按照设定碾压参数对沥青混合料进行碾压，并记录碾压方向，按照碾压参数碾压完成后，将试模连同试件从成型槽中取出，常温放置48h后脱模使试模和试件板分离，最终获得试件板；
其中，手推式压路机的型号为YCYL-600C，碾压参数包括碾压往返次数、碾压压力值和碾压温度，所述碾压往返次数为24次，按照碾压顺序依次为静压8次，振动压实8次、静压8次，碾压温度为150～170℃，碾压压力值为15kN。
具体操作时，为模拟实际路面碾压过程(非大孔隙路面)，试件板碾压往返24次，其中依次静压8次、振动压实8次、静压8次，依次模拟实际路面施工过程中的初压、复压和终压(大孔隙路面除外，与此碾压工艺不同)；在移动试件板和拆模之前，用记号笔记录压路机在试件板的碾压方向。
在沥青混合料的表面铺设厚度为5mm的橡胶块，使沥青混合料表面在碾压过程中受力更加均衡，形成的纹理更加均匀，并尽量减少压路机造成的油渍污染。
步骤四，利用激光纹理扫描仪对试件板上的测点表面进行扫描，以获取路面表面纹理的高程数据；
步骤五，对步骤四中获取的高程数据按照公式一进行坏点处理，以将坏点去除；
Zi1＝min(zn,zm) (1)
其中，i——为路面断面构造线中坏点的编号；
m——断面构造线中距离坏点最近的前一个点的编号；
n——断面构造线中距离坏点最近的后一个点的编号；
zm——为编号为m的断面构造线的高程；
zn——为编号为n的断面构造线的高程；
Zi1——为坏点处理后的断面高程；
为避免沥青路面颜色暗以及集料颗粒遮挡对纹理仪扫描结果的影响，对采集到的路面纹理数据进行坏点处理，采用同一扫描断面该坏点前后的最近的较低高程替换坏点高程。经坏点处理后的路面表面更接近真实的路面形貌。坏点处理前后不同路面断面构造线分别如图2a和2b所示，不同断面构造线坏点率均在10％以内，符合ISO要求(ISO.No.13473-1,1997)。
为避免试件板不平整或扫描过程中试件板放置不平整对纹理结果的影响，按照公式二对高程数据进行倾斜度和偏置误差处理，处理前后示意图如图3a和3b所示；
Zi＝zi-b1i-b0 (2)
其中，b1、b0——分别为倾斜度和偏置误差的消除系数；
i——为断面构造线中采集点的水平编号；
N——为断面构造线在评价长度范围内的采集点总数；
zi——为倾斜度和偏置误差消除前i编号处的断面高程；
Zi——为倾斜度和偏置误差消除后i编号处的断面高程值；
步骤六，将步骤五得到的高程数据进行平均断面深度MPD、均方根Rms、偏斜度Rsk以及不同波段的高程方差的计算，其中不同波段的高程方差包括宏观纹理长波段方差L、宏观纹理短波段方差S、微观纹理波段方差M；
同时，以常用的MPD作为试件板三个测点纹理深度评价指标，以该三个测点的均值作为该试件板的纹理深度代表值，并与标准马歇尔试件的纹理深度进行对比，检验成型的试件板是否与标准马歇尔试件一致。检验参考规范《公路路基路面现场测试规程》(JTG3450-2019)附录B进行，当马歇尔试件的MPD位于所成型试件板三个测点的均值±1.15*标准差范围时，认为试件板成型符合标准要求。
步骤七，根据公路路基路面现场测试规程JTG 3450-2019中T0964-2008摆式仪测试方法，以试件板上测点的中心为控制点，进行BPN测试，得到摆值BPN；
步骤八，采用SPSS软件，以BPN为因变量，以MPD、Rms、Rsk、L、S和M共六个参数为自变量，选用筛选参数较佳的Stepwise逐步回归法，以0.05为进入概率，以0.10为剔除标准，进行参数的选择，就是在回归过程中，每引入一个变量，都会对已在方程中的变量进行检验，本次分析是以0.05为进入概率，以0.10为剔除标准，每次增加一个因变量后，检验模型的自变量在回归方程中是否依旧显著，即若显著性概率小于0.05，则自变量被选入回归方程；若显著性概率大于0.1，则从回归方程中剔除，依次类推，最终建立的回归方程是一个最优的自变量组合。从多个参数中选择出均符和标准的参数，方程内的自变量参数均有统计学意义，而方程外的自变量参数均无统计学意义；
通过上述的多元线性回归方法建立了路面抗滑性能BPN与路面纹理的关系式如下所示，
BPN＝21.828*S-8.041*MPD+59.461 (3)
其中，BPN——摆值；
S——宏观纹理短波段方差；
MPD——测点的平均断面深度。
步骤四中，采用激光纹理扫描仪依次对试件板中间轮迹带上沿行车方向依次扫描的三个位置形成三个测点，每个测点的扫描面积为104mm×72.01mm，扫描每个测点均能够得到163条长度为104mm的均匀分布的断面构造线，相邻测点间距为0.009525mm，相邻两条断面构造线的间距为0.4445mm。测点间距和相邻两条断面构造线间距的选取参照表1所示的数据进行选取。
表1已有激光纹理扫描设备常用的采集参数
进一步的，MPD的计算参照ISO-13473-1标准，每个测点的MPD是该测点163条断面构造线的断面构造深度的均值，试件板的MPD是三个测点的MPD值的平均值。MPD是常用的纹理深度评价指标，Rms和Rsk常用的路面纹理统计评价指标。不同波段的高程方差L、S、M是经过滤波后计算的纹理的方差，作为不同波段纹理的表征。
进一步，采用Butterworth高通和低通滤波器将三个测点的断面构造线滤波成三个不同的波段，分别为波段为5mm～50mm的宏观纹理长波段，波段为0.5mm～5mm的宏观纹理短波段，波段为0.024mm～0.5mm的微观纹理波段。将断面构造线滤波成不同的波段并分别进行考虑，提升了路面纹理评价的精确性。
所述均方根Rms的计算公式如下，
其中，N——断面构造线采集点总数；
Zi——经过倾斜度和偏置误差处理的断面高程值。
所述偏斜度Rsk的计算公式如下，
其中，Rms——均方根；
N——断面构造线采集点总数；
Zi——经过倾斜度和偏置误差处理的断面高程值。
所述不同波段的高程方差计算公式如下，
其中，Variance——不同波段范围纹理的高程方差；
Zi——经特定频段滤波之后的断面高程；
——过倾斜度和偏置误差处理的断面高程值均值。
基于路面纹理评价抗滑性能通过以下实现数据进行具体说明。
选取AC(Asphalt Concrete)和SMA(Stone Mastic Asphalt)两种混合料分别进行试件板的成型，每种混合料类型选取粗、中、细三个级配，每个级配成型2块试件板，共计成型6种级配合计12块试件板。两种混合料的三种级配及油石比如表2所示。
表2AC和SMA两种混合料的三种级配和油石比
将混合料分别按照拌和温度165℃进行拌和，拌和后的混合料放置在恒温箱中按照保温温度165℃进行保温1h，将保温后的混合料装入对应的试模内，并在其表面铺设厚度为5mm的橡胶块，然后通过手推式小型压路机对混合料进行碾压，碾压温度为160℃。为模拟实际路面碾压过程(非大孔隙路面)，试件板碾压往返24次，其中依次静压8次、振动压实8次、静压8次，依次模拟实际路面施工过程中的初压、复压和终压(大孔隙路面除外，与此碾压工艺不同)，碾压完成，最终成型12块试件板。
采用激光纹理扫描仪依次对每块试件板中间轮迹带上沿行车方向依次扫描的三个位置形成三个测点，每个测点的扫描面积为104mm×72.01mm，扫描每个测点均能够得到163条长度为104mm的均匀分布的断面构造线，相邻测点间距为0.009525mm，相邻两条断面构造线的间距为0.4445mm，最终获得12块试件板合计36个测点的路面纹理高程数据。
对12块试件板对应的36组路面纹理高程数据按照公式一进行坏点处理，按照公式二对其进行倾斜度和偏置误差处理，最后得到符合要求的高程数据。
然后根据处理后的高程数据计算路面表面纹理数据，主要包括平均断面深度MPD、均方根Rms、偏斜度Rsk以及不同波段的高程方差的计算。
其中，MPD的计算参照ISO-13473-1标准，每个测点的MPD是该测点163条断面构造线的断面构造深度的均值，试件板的MPD是三个测点的MPD值的平均值；均方根Rms按照公式四计算、偏斜度Rsk按照公式五计算，不同波段的高程方差按照公式六计算。
以常用的MPD作为试件板三个测点纹理深度评价指标，以该三个测点MPD的均值作为该试件板的纹理深度代表值。采用相同的纹理扫描参数和纹理仪，获取标准马歇尔试件的纹理深度进行，并与之对比，检验成型的试件是否与标准马歇尔试件一致。检验参考规范《公路路基路面现场测试规程》(JTG 3450-2019)附录B进行，当马歇尔试件的MPD位于所成型试件板三个测点MPD的均值±1.15*标准差范围时，认为试件板成型符合标准要求。12块板纹理数据分析表明，上述成型工艺所成型的试件板纹理与标准马歇尔试件基本一致，如表3所示。
表3 12块试件板纹理与标准马歇尔试件对比结果
对12块试件板合计36个测点按照公路路基路面现场测试规程JTG 3450-2019中T0964-2008摆式仪测试方法，以试件板上测点的中心为控制点，进行BPN测试，得到摆值BPN。
通过以上步骤共获得36组纹理数据和抗滑性能数据，相关数据见表4。
表4 36组纹理和抗滑详细数据
采用SPSS软件对上述36组数据进行处理，以BPN为因变量，以MPD、Rms、Rsk、L、S和M共六个参数为自变量，选用筛选参数较佳的Stepwise逐步回归法，以0.05为进入概率，以0.10为剔除标准，通过多元线性回归方法建立路面抗滑性能BPN与路面纹理的联系，并剔除相关性差的指标，得到公式三，相关系数为0.50，实现基于纹理数据的路面抗滑性能BPN的预测。
对上述公式三进行验证，采用与上述12块试件板相同材料、相同成型工艺，成型不同于上述级配的试件板2块，如表5所示，按照上述步骤得到6组纹理数据和抗滑数据，如表6所示，把纹理数据代入公式三中，计算得出BPN预估值，该2块试件板不同测点BPN预测值和实测值结果如表7所示，数据偏差率在1.1％～5.7％，偏差率可以接受，因此采用公式三建立的BPN与纹理数据的关联，从路面纹理角度预测抗滑性能有效。
表5不同于12块试件板的2块试件板级配和油石比
表6回归模型验证使用的6组原始数据
表7不同测点BPN预测值和实测值结果
本方法建立了基于纹理的抗滑性能预测模型，可以实现基于路面纹理的抗滑性能精准预测，避免抗滑性能测试导致交通中断或测试人员安全区受威胁。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
本领域技术人员不难理解，本发明包括上述说明书的发明内容和具体实施方式部分以及附图所示出的各部分的任意组合，限于篇幅并为使说明书简明而没有将这些组合构成的各方案一一描述。凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (10)
1.一种基于路面纹理评价抗滑性能的分析方法，其特征在于，包括以下步骤：
步骤一，按照T0702-2011成型标准马歇尔试件，以计算出标准试件对应的标准密度，将标准密度作为试件板的密度，选取试件板的尺寸并计算其体积，将体积和标准密度相乘并乘以系数1.03，最终求得试件板所需的沥青混合料的质量，按照试件板的尺寸制作试模；
步骤二，根据T0702-2011选择沥青混合料的拌和温度，根据拌和温度将沥青混合料进行拌和，将拌和后的沥青混合料放置在恒温箱中并按照保温温度和保温时间进行保温；
步骤三，在水泥地面上按照试件板的试模尺寸切割成型槽，以放置试模，将保温后的沥青混合料放入试模内，并通过手推式小型压路机按照设定的碾压参数对沥青混合料进行碾压，最终获得试件板；
步骤四，采用激光纹理扫描仪对试件板中间轮迹带上沿行车方向依次扫描三个位置形成三个测点，并对三个测点的表面进行扫描，以获取路面表面纹理的高程数据；采用Butterworth高通和低通滤波器将三个测点的断面构造线滤波成三个不同的波段；
步骤五，对步骤四中获取的高程数据按照公式一进行坏点处理，以将坏点去除；按照公式二对高程数据进行倾斜度和偏置误差处理；
公式一为：Zi1＝min(zn,zm)
其中，i——断面构造线中坏点的编号；
m——断面构造线中距离坏点最近的前一个点的编号；
n——断面构造线中距离坏点最近的后一个点的编号；
zm——编号为m的断面构造线的高程，mm；
zn——编号为n的断面构造线的高程，mm；
Zi1——坏点处理后的断面高程，mm；
公式二为：Zi＝zi-b1i-b0
其中，b1、b0——分别为倾斜度和偏置误差的消除系数；
i——断面构造线中采集点的水平编号；
N——断面构造线在评价长度范围内的采集点总数；
zi——倾斜度和偏置误差消除前i编号处的断面高程，mm；
Zi——倾斜度和偏置误差消除后i编号处的断面高程，mm；
步骤六，将步骤五得到的高程数据进行平均断面深度MPD、均方根Rms、偏斜度Rsk以及不同波段的高程方差的计算，其中不同波段的高程方差包括宏观纹理长波段方差L、宏观纹理短波段方差S和微观纹理波段方差M；
步骤七，以试件板上测点的中心为控制点，进行BPN测试，得到摆值BPN；
步骤八，以BPN为因变量，以MPD、Rms、Rsk、L、S和M六个参数为自变量，选用逐步回归法，以0.05为进入概率，以0.10为剔除标准，进行参数的选择，通过上述的逐步回归法建立路面抗滑性能BPN与路面纹理的关系式；
该关系式为：BPN＝21.828*S-8.041*MPD+59.461
其中，BPN——摆值，无量纲；
S——宏观纹理短波段方差，mm2；
MPD——测点的平均断面深度，mm。
2.如权利要求1所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，步骤二中，拌和温度为155-175℃，保温时间为1h，保温温度比碾压温度高5℃。
3.如权利要求2所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，步骤三中，碾压参数包括碾压往返次数、碾压压力值和碾压温度，所述碾压往返次数为24次，按照碾压顺序依次为静压8次、振动压实8次、静压8次，所述碾压温度为150-170℃，所述碾压压力值为15kN。
4.如权利要求3所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，在碾压之前，在沥青混合料的表面铺设厚度为5mm的橡胶块。
5.如权利要求4所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，步骤四中，每个测点的扫描面积均为104mm×72.01mm，扫描每个测点均能够得到163条长度为104mm的均匀分布的断面构造线，相邻测点的间距为0.009525mm，相邻两条断面构造线的间距为0.4445mm。
6.如权利要求5所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，所述MPD的计算参照ISO-13473-1标准，每个测点的MPD是指该测点163条断面构造线的断面构造深度的平均值，试件板的MPD是指三个测点的断面构造深度的平均值。
7.如权利要求6所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，三个不同的波段分别为：波段为5-50mm的宏观纹理长波段、波段为0.5-5mm的宏观纹理短波段、波段为0.024-0.5mm的微观纹理波段。
8.如权利要求7所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，所述均方根Rms的计算公式为
其中，Rms——均方根，mm；
N——断面构造线采集点总数；
Zi——倾斜度和偏置误差消除后i编号处的断面高程，mm。
9.如权利要求8所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，所述偏斜度Rsk的计算公式为
其中，Rsk——偏斜度，无量纲；
N——断面构造线采集点总数；
Zi——倾斜度和偏置误差消除后i编号处的断面高程，mm。
10.如权利要求9所述的一种基于路面纹理评价抗滑性能的分析方法，其特征在于，所述不同波段的高程方差的计算公式为
其中，Variance——不同波段的高程方差，mm2；
Hi——Zi经特定频段滤波之后的断面高程，mm；
——Zi经特定频段滤波之后的断面高程的平均值，mm。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111564062.1A CN114441436B (zh) | 2021-12-20 | 2021-12-20 | 一种基于路面纹理评价抗滑性能的分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111564062.1A CN114441436B (zh) | 2021-12-20 | 2021-12-20 | 一种基于路面纹理评价抗滑性能的分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114441436A CN114441436A (zh) | 2022-05-06 |
| CN114441436B true CN114441436B (zh) | 2024-04-16 |
Family
ID=81364650
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111564062.1A Active CN114441436B (zh) | 2021-12-20 | 2021-12-20 | 一种基于路面纹理评价抗滑性能的分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114441436B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114994054B (zh) * | 2022-05-26 | 2023-06-20 | 哈尔滨工业大学 | 路表纹理平均构造深度与平均断面深度关系的确定方法 |
| CN116973302B (zh) * | 2023-09-20 | 2024-01-23 | 四川省交通勘察设计研究院有限公司 | 一种路面抗滑性能检测方法、装置、设备和介质 |
| CN117074291B (zh) * | 2023-10-17 | 2024-01-02 | 西南交通大学 | 一种非接触式纹理摩擦预测方法 |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928334A (zh) * | 2012-08-14 | 2013-02-13 | 北京建筑工程学院 | 沥青混合料抗滑耐久性能的评价方法 |
| KR20140046358A (ko) * | 2012-10-10 | 2014-04-18 | 이명진 | 산사태 위험성 분석 시스템 및 방법 |
| CN104155201A (zh) * | 2014-08-12 | 2014-11-19 | 北京建筑大学 | 一种融雪路面评价方法 |
| CN104504721A (zh) * | 2015-01-08 | 2015-04-08 | 中国科学院合肥物质科学研究院 | 一种基于Gabor小波变换纹理描述的非结构化道路检测方法 |
| CN105716969A (zh) * | 2016-04-12 | 2016-06-29 | 长安大学 | 一种基于三维激光检测技术的沥青路面打毛效能评价方法 |
| CN105957064A (zh) * | 2016-04-24 | 2016-09-21 | 长安大学 | 一种沥青路面表面构造二维测试评价系统及方法 |
| CN106918553A (zh) * | 2017-03-26 | 2017-07-04 | 长安大学 | 基于aims集料成像系统的粗集料抗滑性能评价方法 |
| EP3201396A1 (en) * | 2014-10-03 | 2017-08-09 | Flint Trading, Inc. | Preformed thermoplastic pavement marking and method for high skid resistance with maintained high retroreflectivity |
| CN107886721A (zh) * | 2017-11-09 | 2018-04-06 | 西华大学 | 一种交通事故道路状况数据采集方法 |
| CN108562727A (zh) * | 2018-04-24 | 2018-09-21 | 武汉科技大学 | 一种沥青混合料表面纹理磨光行为的评价方法 |
| CN110055867A (zh) * | 2018-10-29 | 2019-07-26 | 交通运输部公路科学研究所 | 一种沥青铺面抗滑安全风险实时感知预警方法 |
| CN110738636A (zh) * | 2019-09-12 | 2020-01-31 | 北京建筑大学 | 基于数字图像技术测试沥青混合料抗滑性能的方法及系统 |
| CN112270677A (zh) * | 2020-11-16 | 2021-01-26 | 武汉大学 | 一种路面宏观纹理获取及典型指标计算方法 |
| CN113218332A (zh) * | 2021-04-30 | 2021-08-06 | 长安大学 | 一种基于降雨条件下积沙沥青路面宏观纹理的评价方法 |
| CN113418861A (zh) * | 2021-06-17 | 2021-09-21 | 东南大学 | 一种用于室内研究路面材料抗滑性能衰变规律的装置及方法 |
| CN113686858A (zh) * | 2021-08-31 | 2021-11-23 | 山东省交通科学研究院 | 一种基于表面纹理构造的路面抗滑性能评价方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6599057B2 (en) * | 2001-06-14 | 2003-07-29 | Kmc Enterprises, Inc. | Cold in-place recycling of bituminous material |
| US20150140294A1 (en) * | 2004-04-02 | 2015-05-21 | Flint Trading, Inc. | Preformed Thermoplastic Pavement Marking and Method for High Skid Resistance with Maintained High Retroreflectivity |
| US8085987B2 (en) * | 2007-10-26 | 2011-12-27 | Ahmed Shalaby | Method and tool for surface texture evaluation |
| EP3007906A2 (en) * | 2013-06-09 | 2016-04-20 | Apple Inc. | Laser-formed features |
| US20190362510A1 (en) * | 2018-05-24 | 2019-11-28 | Lu Sun | Method and system for evaluating friction coefficient and skid resistence of a surface |
-
2021
- 2021-12-20 CN CN202111564062.1A patent/CN114441436B/zh active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928334A (zh) * | 2012-08-14 | 2013-02-13 | 北京建筑工程学院 | 沥青混合料抗滑耐久性能的评价方法 |
| KR20140046358A (ko) * | 2012-10-10 | 2014-04-18 | 이명진 | 산사태 위험성 분석 시스템 및 방법 |
| CN104155201A (zh) * | 2014-08-12 | 2014-11-19 | 北京建筑大学 | 一种融雪路面评价方法 |
| EP3201396A1 (en) * | 2014-10-03 | 2017-08-09 | Flint Trading, Inc. | Preformed thermoplastic pavement marking and method for high skid resistance with maintained high retroreflectivity |
| CN104504721A (zh) * | 2015-01-08 | 2015-04-08 | 中国科学院合肥物质科学研究院 | 一种基于Gabor小波变换纹理描述的非结构化道路检测方法 |
| CN105716969A (zh) * | 2016-04-12 | 2016-06-29 | 长安大学 | 一种基于三维激光检测技术的沥青路面打毛效能评价方法 |
| CN105957064A (zh) * | 2016-04-24 | 2016-09-21 | 长安大学 | 一种沥青路面表面构造二维测试评价系统及方法 |
| CN106918553A (zh) * | 2017-03-26 | 2017-07-04 | 长安大学 | 基于aims集料成像系统的粗集料抗滑性能评价方法 |
| CN107886721A (zh) * | 2017-11-09 | 2018-04-06 | 西华大学 | 一种交通事故道路状况数据采集方法 |
| CN108562727A (zh) * | 2018-04-24 | 2018-09-21 | 武汉科技大学 | 一种沥青混合料表面纹理磨光行为的评价方法 |
| CN110055867A (zh) * | 2018-10-29 | 2019-07-26 | 交通运输部公路科学研究所 | 一种沥青铺面抗滑安全风险实时感知预警方法 |
| CN110738636A (zh) * | 2019-09-12 | 2020-01-31 | 北京建筑大学 | 基于数字图像技术测试沥青混合料抗滑性能的方法及系统 |
| CN112270677A (zh) * | 2020-11-16 | 2021-01-26 | 武汉大学 | 一种路面宏观纹理获取及典型指标计算方法 |
| CN113218332A (zh) * | 2021-04-30 | 2021-08-06 | 长安大学 | 一种基于降雨条件下积沙沥青路面宏观纹理的评价方法 |
| CN113418861A (zh) * | 2021-06-17 | 2021-09-21 | 东南大学 | 一种用于室内研究路面材料抗滑性能衰变规律的装置及方法 |
| CN113686858A (zh) * | 2021-08-31 | 2021-11-23 | 山东省交通科学研究院 | 一种基于表面纹理构造的路面抗滑性能评价方法 |
Non-Patent Citations (7)
| Title |
|---|
| Development and testing of a multivariable accelerated abrasion machine to characterize the polishing wear of pavement by tires;Ren WY, et al;《Surface topography metrology and properities》;20190802;第7卷(第3期);第035006页 * |
| Pavement skid resistance as a function of pavement surface and aggregate texture properties;Li Q, et al;《International journal of pavement engineering》;20200915;第21卷(第10期);第1159-1169页 * |
| 基于分段变维的沥青路表纹理磨光行为分析;周兴林，等;《中国公路学报》;20190430;第32卷(第4期);第187-195页 * |
| 基于微观形貌特征的沥青路面抗滑性能评价研究进展;初秀民，等;《交通与计算机》;20070228(第01期);第61-65页 * |
| 基于高精度激光断面高程的沥青路面MTD测量;丁世海，等;《东南大学学报(自然科学版)》;20200120;第50卷(第01期);第137-142页 * |
| 沥青路面三维纹理分形维数及其抗滑性能;钱振东，等;《中南大学学报》;20161031;第47卷(第10期);第3590-3596页 * |
| 沥青路面三维纹理构造特征研究;刘林，等;《公路交通科技(应用技术版)》;20190415;第15卷(第04期);第46-50页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114441436A (zh) | 2022-05-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114441436B (zh) | 一种基于路面纹理评价抗滑性能的分析方法 | |
| US11492766B2 (en) | Method of hot recycling repairing by optimizing proportion of asphalt mixture on pavement alignment variation section | |
| Fwa et al. | Comparison of permeability and clogging characteristics of porous asphalt and pervious concrete pavement materials | |
| Wang et al. | Evaluation of the polishing resistance characteristics of fine and coarse aggregate for asphalt pavement using Wehner/Schulze test | |
| Kassem et al. | Development of predictive model for skid loss of asphalt pavements | |
| Cross et al. | Selection of aggregate properties to minimize rutting of heavy duty pavements | |
| Prowell et al. | Evaluation of circular texture meter for measuring surface texture of pavements | |
| Bendtsen et al. | Acoustic aging of asphalt pavements: A Californian/Danish comparison | |
| Doyle et al. | Laboratory assessment of skid resistance for high RAP content warm mixed asphalt | |
| CN105865871A (zh) | 一种沥青混合料标准试件模拟施工态碾压成型方法及设备 | |
| Kutay et al. | Development of an acceptance test for chip seal projects | |
| Nicholls | Asphalt mixture specification and testing | |
| Kutay et al. | An acceptance test for chip seal projects based on image analysis. | |
| Miller et al. | Using Georgia loaded-wheel tester to predict rutting | |
| Kowalski et al. | Identification of laboratory technique to optimize Superpave HMA surface friction characteristics | |
| Obulareddy | Fundamental characterization of Louisiana HMA mixtures for the 2002 mechanistic-empirical design guide | |
| Kowalski et al. | Development of a laboratory procedure to evaluate the influence of aggregate type and mixture proportions on the frictional characteristics of flexible pavements | |
| CN115081813A (zh) | 一种钢渣沥青混合料路面磨耗层抗滑性能评价方法 | |
| Jacobson | Study of the wear resistance of bituminous mixes to studded tyres: Tests with slabs of bituminous mixes inserted in roads and in the VTIs road simulator | |
| Kassem et al. | Effect of compaction on skid resistance of asphalt pavements | |
| Pranjić et al. | Pavement surface macrotexture analysis | |
| Ozdemir | An acceptance test for chip seal projects based on image analysis | |
| CN112612994B (zh) | 沥青稳定碎石2倍标准差确定级配方法 | |
| Kassem et al. | Laboratory evaluation of friction loss and compactability of asphalt mixtures. | |
| Woodward et al. | Use of 3D modelling techniques to better understand road surface textures |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |