粉煤灰-矿渣基地聚物混凝土的抗碳化性能

原元; 赵人达; 占玉林; 李福海; 成正清; 李健

doi:10.3969/j.issn.0258-2724.20191151

粉煤灰-矿渣基地聚物混凝土的抗碳化性能

doi: 10.3969/j.issn.0258-2724.20191151

原元^1,,
赵人达¹,
占玉林^{1, 2},
李福海^1, ,,
成正清¹,
李健¹

1.
西南交通大学土木工程学院，四川成都 610031
2.
西南交通大学土木工程材料研究所，四川成都 610031

基金项目: 国家自然科学基金（51778531）；四川省科研计划（2019YFG0001，2019YJ0219）；教育部产学合作协同育人项目（201801098032）

详细信息

作者简介:
原元（1992—），男，博士研究生，研究方向为高使用性混凝土材料性能及结构应用，E-mail：1181491115@qq.com

通讯作者:
李福海（1979—），男，高级工程师，博士，研究方向为新型材料及混凝土结构耐久性，E-mail：lifuhai2007@home.swjtu.edu.cn

中图分类号: TU528.41
计量
- 文章访问数: 667
- HTML全文浏览量: 263
- PDF下载量: 22
- 被引次数: 0
出版历程
- 收稿日期: 2019-12-02
- 修回日期: 2020-05-27
- 网络出版日期: 2020-09-21
- 刊出日期: 2020-09-21

Carbonation Resistance of Fly Ash-Slag Based Geopolymer Concrete

YUAN Yuan^1
,,
ZHAO Renda¹,
ZHAN Yulin^{1, 2},
LI Fuhai^{1
, ,},
CHENG Zhengqing¹,
LI Jian¹

1.
School of Civil Engineering，Southwest Jiaotong University，Chengdu 610031, China
2.
Institute of Civil Engineering Materials，Southwest Jiaotong University，Chengdu 610031, China

摘要

摘要:
F级粉煤灰-矿渣基地聚物混凝土，即GPC-10（矿渣掺量10%，80 °C高温养护）和GPC-50（矿渣掺量50%，标准养护）力学性能良好，为进一步研究其抗碳化性能, 首先，对这两种地聚物混凝土进行了快速碳化试验，并与作为对照组的普通水泥混凝土(OPCC)进行了比较，通过抗压强度和劈裂抗拉强度评价了碳化对混凝土的损伤；其次，为分析损伤原因，分别通过X射线能谱分析（EDS）和压汞测试（MIP），对碳化后的成分和孔结构进行了研究；最后，建立了两种地聚物混凝土的碳化模型. 研究结果表明：相比OPCC，地聚物混凝土的抗碳化能力薄弱，尤其是钙含量较高的GPC-50，其主要产物C—A—S—H会与CO₂反应而发生分解，导致孔隙率增大，进而加快了碳化速率，且碳化深度与时间呈线性关系；OPCC、GPC-10以及GPC-50的28 d碳化深度分别达到了2.0、9.2、18.8 mm.
- F级粉煤灰 /
- 矿渣 /
- 地聚物混凝土 /
- 快速碳化 /
- 孔结构 /
- 碳化模型
Abstract:
Two kinds of Class F fly ash-slag based geopolymer concrete, namely GPC-10 (10% slag content, 80 °C high-temperature curing) and GPC-50 (50% slag content, standard curing), are characterized by good mechanical properties. To understand better their carbonation resistance, rapid carbonation tests were carried out on the two kinds of geopolymer concrete in comparison with the ordinary Portland cement concrete (OPCC). The damage of concrete was evaluated in terms of compressive strength and splitting tensile strength. To analyze the damage cause, the composition and pore structure of the carbonized materials were investigated using X-ray energy spectroscopy (EDS) and mercury intrusion test (MIP), respectively. On this basis, carbonation models of the two kinds of concrete were established. The results show that compared with the OPCC, the geopolymer concrete has a weak carbonation resistance, especially for GPC-50, the type with high calcium content, in which the main product C—A—S—H was decomposed during carbonation, leading to an increase in porosity and thus accelerating the carbonation rate, and the carbonation depth has a linear relation with time. The 28 d carbonization depth of OPCC, GPC-10 and GPC-50 reached 2.0, 9.2 and 18.8 mm, respectively.
- class F fly ash /
- slag /
- geopolymer concrete /
- rapid carbonation test /
- pore structure /
- carbonation models

HTML全文

随着无人驾驶技术的发展，未来道路上的智能网联车（connected automated vehicle，CAV）将逐渐普及. 但最新研究表明^[1]，预计2045年道路上L4级别的CAV渗透率仅能达到24.8%. 因此，未来很长一段时间内，道路上将普遍呈现由常规车（regular vehicle，RV）、网联人工驾驶车（connected vehicle，CV）和CAV组成的混合交通流. 网联车（CV和CAV）能够提供海量的车辆轨迹数据，而这些数据蕴含了丰富的时空交通信息，是交通状态估计的基础数据之一. 但由于RV无法提供自身的轨迹数据，在智能网联环境下全样本车辆轨迹数据仍难以获得. 因此，如何利用已有的轨迹数据来重构所有车辆的轨迹具有重要的意义^[2].

近年来，有许多学者对车辆轨迹重构问题进行了探讨，主要为基于经典交通流理论和基于统计概率模型的方法. 这些方法按轨迹重构数据来源可分为三大类：定点检测器、移动检测器和多源数据. Coifman^[3]基于双线圈检测器数据和交通流理论，通过在时空图中估计路段行程时间重构车辆轨迹. Rao等^[4]利用车牌识别数据构建了一种车辆路径重构方法，进一步可基于路径对交通网络的历史OD （original-destination）进行估计. Hao等^[5]、Wan等^[6]、Shan等^[7]基于概率论的思想，分别采用随机模型、期望最大化算法（expectation-maximum，EM）及最大似然估计，从微观角度重构浮动车任意相邻采样点之间的运动轨迹. Li等^[8]提出了基于K最近邻回归算法利用不完整轨迹来扩充轨迹样本量. Xie等^[9]基于修正的交通信号配时方案、车辆通过时间和车辆行程时间构建微观仿真环境，实现了对干道车辆轨迹重构. Mehran等^[10]通过融合秒级浮动车数据、音视频数据（audio video interleaved，AVI）及信号控制参数，提出了一种基于变分理论城市干道车辆运行轨迹估计方法，并利用东京一条包含五个交叉口的单行道实测数据进行了方法验证和误差分析. 唐克双等^[11]通过融合视频和定点检测器流量数据，基于交通流理论和交通仿真思想，构建了一种不依赖高频浮动车数据，适用于多车道的轨迹重构方法. 综上分析可知，现有研究均能够对车辆轨迹数据进行有效重构，但均只考虑了人工驾驶环境，未考虑未来智能网联混合交通流的情况. 随着智能网联技术的发展，CAV除了能提供自身轨迹数据外，还能通过感知设备（摄像头和雷达）提供自身一定范围内的车辆轨迹数据. 因此，CAV采集的移动传感数据将会给轨迹重构提供一种新的思路. 如何充分考虑混合交通流环境下不同车辆类型提供的数据，构建全样本车辆轨迹重构模型，是未来CAV环境下车辆轨迹重构亟待解决的关键问题.

为解决该问题，本文拟研究混有RV、CV和CAV的混合交通流全样本车辆轨迹重构方法. 充分利用CAV和CV提供的轨迹数据信息，提出基于跟驰特性的混合交通流轨迹重构模型，并通过仿真试验验证模型在不同交通流密度和CAV渗透率条件下的合理性与有效性.

1. 跟驰特性分析

在本文研究中，不考虑车辆换道的影响，选取较为简单的场景来阐述轨迹重构思想. 在微观交通流分析中，一般采用跟驰模型描述车辆之间的跟驰行为. 该模型认为当前车辆是否加减速与本身速度、前后车速度差和位置差有关. 因此，在已知部分轨迹的情况下，采用该模型可分析相邻轨迹中存在其他车辆轨迹的可能性.

跟驰模型主要可分为五种类型：刺激反应模型、安全间距模型、社会力模型、优化速度模型和低阶线性模型^[12]，均可作为本文轨迹重构的基础模型. 而作为社会力模型的一种，智能驾驶员模型（intelligent driver model，IDM）能够很好地描述熟练驾驶员的驾驶习惯，应用范围更广. 因此，本文以IDM为例描述车辆之间的跟驰特性，如式（1）所示.

$\begin{split} & {a_p}\left( t \right) = {a_{\rm{M}}}\left[ {1 - {{\left( {\dfrac{{{v_p}\left( t \right)}}{{{v_{\rm{f}}}}}} \right)}^\delta } - } \right. \\&\quad \left. {{{\left( {\dfrac{{{s_0} + T{v_p}\left( t \right) + \dfrac{{{v_p}\left( t \right)\left[ {{v_p}\left( t \right) - {v_{p - 1}}\left( t \right)} \right]}}{{2\sqrt {{a_{\rm{M}}}b} }}}}{{{x_{p - 1}}\left( t \right) - {x_p}\left( t \right)}}} \right)}^2}} \right] \end{split},$

(1)

式中： ${a_p}(t)$ 为时刻t车辆 $p$ 的加速度； ${v_{\rm{f}}}$ 为期望速度； ${a_{\rm{M}}}$ 为最大加速度；T为安全车头时距； $b$ 为舒适减速度的绝对值； ${s_0}$ 为最小停车间距； ${v_p}(t)$ 为时刻t车辆 $p$ 的速度； ${x_p}(t)$ 为时刻t车辆 $p$ 的位置； $\delta$ 为加速度指数（一般取4）.

2. 基于IDM的车辆状态估计

2.1 数据环境分析

本文研究的混合交通流中有三种类型的车辆：RV、CV和CAV. 其中，CAV除提供自身的轨迹数据外，还可通过感知设备获取周围一定范围内的车辆位置和状态信息，因此，能提供周围一定范围内车辆的轨迹数据. 混合交通流中不同车辆组成如图1所示. CAV探测范围内的车辆信息均可采集，不论是RV还是CV. CV由于不具备感知设备，所以只能提供自身的轨迹数据. 而RV由于其不具备或未打开联网功能，所以既不能提供自身的轨迹数据，也不能提供周围车辆的轨迹数据.

图 1 混合交通流中不同车辆组成

Figure 1. Composition of different types of vehicles in mixed traffic flow

下载: 全尺寸图片幻灯片

从以上分析可知，在该混合交通流环境下，主要需要对CAV探测范围之外的RV进行轨迹重构. 混合交通流轨迹时空图如所示. CAV两辆（CAV1、CAV2），CV三辆（CV1、CV2、CV3）；将阴影区域A中前一辆CAV探测到的最后一辆车记为 ${{\rm{V}}_{n - 1}}$ ，后面一辆CV记为 ${{\rm{V}}_n}$ ；红色虚线表示CAV轨迹；蓝色点线表示CV轨迹；黑色实线表示RV轨迹；阴影区域A和B即为需要重构轨迹的区域. 在跟驰模型中，本车的跟驰行为主要与前车有关，即前车的轨迹直接影响本车的轨迹. 以阴影区域A为例，若区域A中无其他车辆， ${{\rm{V}}_n}$ 会跟随 ${{\rm{V}}_{n - 1}}$ ，本车 ${{\rm{V}}_n}$ 与前车 ${{\rm{V}}_{n - 1}}$ 构成跟驰行为，可采用跟驰模型进行描述. 但若区域A中还存在未被观测的RV，则本车 ${{\rm{V}}_n}$ 用跟驰模型计算的加速度会与实际加速度存在较大差异. 因此，构建轨迹重构模型试图寻找区域A内的最佳轨迹条数，来实现这种差异的最小化.

图 2 混合交通流轨迹时空图

Figure 2. Spatio-temporal diagram of mixed traffic flow trajectories

下载: 全尺寸图片幻灯片

CV提供的轨迹数据可分为是否在CAV的检测范围内. 当CV在CAV的检测范围内，如图2中的CV1和CV3，这两辆车的轨迹可直接由CAV观测，故该车提供的数据有无对整个模型的求解没有影响. 当CV不在CAV检测范围内时，如图2中的CV2，此时需要利用该轨迹和CAV提供的轨迹重构阴影区域A和B中的车辆轨迹.

2.2 估计插入车辆的速度

车辆时空轨迹示意如所示，图中： ${A_1}{\text{,}} {A_2}{\text{,}}\cdots{\text{,}} {A_k}$ 为时空区域； $k$ 为时间总步数； $\Delta t$ 为时间步长； ${x_n}\left( {\Delta t} \right)$ 为 $\Delta t$ 时 ${{\rm{V}}_n}$ 的位置. 考虑到交通流稳定时，前后车的速度差别较小. 因此，第k个时空区域 ${A_k}$ 中车辆 ${{\rm{V}}_n}$ 的速度 $v_n$ 的估计如式（2）所示^[13].

图 3 车辆轨迹示意

Figure 3. Schematic diagram of vehicle trajectories

下载: 全尺寸图片幻灯片

$\begin{split} & v_n = \frac{{d\left( {{A_k}} \right)}}{{T\left( {{A_k}} \right)}} =\\ &\quad \frac{{ {{x_{n - 1}}\left( {t + \Delta t} \right) - {x_{{{n}} - 1}}(t)} + {{x_{{n}}}\left( {t + \Delta t} \right) - {x_{{n}}}(t)} }}{{2\Delta t}} \end{split},$

(2)

式中： $d\left( {{A_k}} \right)$ 和 $T\left( {{A_k}} \right)$ 分别为区域 ${A_k}$ 的位置差和时间差.

由于插入的RV在区域 ${A_k}$ 中，因此，该车在时刻 $k\Delta t$ 的速度可以用区域 ${A_k}$ 的速度 $v\left( {{A_k}} \right)$ 进行估计.

2.3 估计插入车辆的位置

在任意时刻，插入车辆的位置都需满足与前后车之间的安全间距. 因此，将V_n前插入的第i辆车记为 ${{\rm{V}}_{n,i}}$ ， ${{\rm{V}}_{n,i}}$ 的位置x_n，i只能出现在一定的区间范围内，如式（3）所示.

${x_n}\left( t \right) + {s_0} + l \leqslant {x_{n{\text{,}}i}}\left( t \right) \leqslant {x_{n - 1}}\left( t \right) - {s_0} - l,$

(3)

式中：l为车身长度.

当获得插入车辆的速度和可能位置区间后，可将 ${{\rm{V}}_{n,i}}$ 和 ${{\rm{V}}_n}$ 分别作为前车和后车，而 ${{\rm{V}}_{n,i}}$ 插入位置 ${x_{n,i}}\left( t \right)$ 会对后车 ${{\rm{V}}_n}$ 产生影响，这种影响会反映到后车的加速度上. 因此，可以根据IDM代入式（1）计算 ${{\rm{V}}_n}$ 的加速度 ${{\hat a}_n}\left( {{x_{n,i}}\left( t \right)} \right)$ .

当 ${x_{n,i}}\left( t \right)$ 在可行域内取不同的值时， ${{\rm{V}}_n}$ 会得到不同的加速度. 而当估计加速度 ${\hat a_n}\left( {{x_{n,i}}(t)} \right)$ 与实际加速度 ${a_n}\left( t \right)$ 最接近时，对应的 ${x_{n,i}}\left( t \right)$ 即为要插入车辆 ${{\rm{V}}_{n,i}}$ 的位置.

因此，可构建目标函数求得最优的插入车辆位置，同时考虑到车辆不同时刻的位置即为车辆的轨迹，可选用 ${{\rm{V}}_n}$ 的实际加速度和估计加速度之间均方根误差最小为优化目标，如式（4）所示.

$\min \sqrt {\frac{1}{k}{{\sum\limits_{t = 0}^{k\Delta t} {\left[ {{{\hat a}_n}\left( {{x_{n,i}}(t)} \right) - {a_n}\left( t \right)} \right]^2} }}}.$

(4)

综上，该优化模型的目标函数为式（4），约束条件为 ${{x_{n,i}}}$ 和 ${{\hat a}_n}( {{x{_{n,i}(t)}}})$ ，为典型的非线性优化模型，可用MATLAB自带的函数直接进行求解.

3. 基于IDM的车辆轨迹重构

由第2节的分析可知，得出插入车辆后，可由模型计算其插入的位置和速度. 因此，车辆轨迹重构模型即为需要寻找已知轨迹中间未知轨迹的数量，插入车辆轨迹示意如图4所示. 当插入车辆数已知时，可采用第2节中的方法依次求得每辆插入车辆轨迹的具体位置和速度. 考虑到车辆之间的最小安全间距，插入车辆数m的最大可能取值为N，如式（5）所示.

图 4 插入车辆轨迹示意

Figure 4. Schematic diagram of inserted vehicle trajectories

下载: 全尺寸图片幻灯片

$\begin{split} & N = {\rm{ent}}\left( {\rm{min}} \left\{ {\dfrac{{{x_{n - 1}}\left( 0 \right) - {x_n}\left( 0 \right)}}{{{s_0} + l}},\dfrac{{x_{n - 1}}\left( {\Delta t} \right) - {x_n}\left( {\Delta t} \right)}{{{s_0} + l}},}\cdots, \right. \right. \\ &\quad \left. \left. { \dfrac{{{x_{n - 1}}\left( {k\Delta t} \right) - {x_n}\left( {k\Delta t} \right)}}{{{s_0} + l}}}\right\}\right) . \\[-20pt] \end{split}$

(5)

在车辆 ${{\rm{V}}_n}$ 和车辆 ${{\rm{V}}_{n - 1}}$ 之间依次插入的m辆车( ${{\rm{V}}_{n,1}},{{\rm{V}}_{n,2}},\cdots,{{\rm{V}}}_{n,m}$ ）， ${{\rm{V}}_{n,m}}$ 为距离 ${{\rm{V}}_n}$ 最近的车辆，即为 ${{\rm{V}}_n}$ 的前车. 插入m辆车后，类比式（2），可粗略估计 ${{\rm{V}}_{n,m}}$ 的速度和位置分别如式（6）、（7）所示.

${v_{n,m}}\left( t \right) = {v_n}\left( t \right) - \frac{{ {{v_n}\left( t \right) - {v_{n - 1}}\left( t \right)} }}{{ {m + 1} }},$

(6)

${x_{n,m}}\left( t \right) = {x_n}\left( t \right) - \frac{{{{x_n}\left( t \right) - {x_{n - 1}}\left( t \right)} }}{{{m + 1} }}.$

(7)

在估计插入的速度和位置后，可重构 ${{\rm{V}}_{n,m}}$ 的轨迹. 根据IDM由 ${{\rm{V}}_{n,m}}$ 的轨迹可求出 ${{\rm{V}}_n}$ 的加速度. 插入m辆RV后代入式（1）可计算 ${{\rm{V}}_n}$ 的理论加速度.

此处，选用 ${{\rm{V}}_n}$ 加速度的均方根误差作为衡量标准，则均方根误差最小时对应的 $m$ 即为最佳插入的车辆数，如式（8）所示.

$\mathop {{\rm{arg min}}}\limits_m \sqrt {\frac{1}{k}{{\sum\limits_{t = 0}^{k\Delta t} {\left[ {{{\hat a}_n}\left( {{x_{n,i}}(t)} \right) - {a_n}\left( t \right)} \right]^2} }}}.$

(8)

因此，该问题也为一个非线性优化问题，目标函数为式（8），约束条件为式（1）、（6）、（7）. 综上，当确定最佳插入车辆数量以后，可采用第2节中的方法对每一插入车辆轨迹进行精确地重构.

4. 数值仿真分析

4.1 仿真参数设置

为验证模型的有效性，基于IDM设计元胞自动机进化规则，模拟单车道高速公路交通流运行规律. 所有车辆之间跟驰行为均用IDM描述^[14-15]. 仿真路段长度1000 m，元胞长度0.1 m，仿真时间1000 s，仿真步长1 s. 为使仿真数据更趋向于真实数据，在仿真过程给车辆设置随机减速模拟真实道路上的车流运行状况.

本文共设计两种基准试验方案，第一种方案是将CAV的渗透率恒定为8%，CV渗透率恒定为20%，改变试验中的交通流密度，研究不同交通流密度对模型的影响. 第二种方案将交通流密度设定为70 辆/km，改变CAV和CV的渗透率，研究不同渗透率对模型的影响. 本文IDM的参数^[16]取值为： ${s_0}$ =2 m， ${a_{\rm{M}}}$ =1 m/s²， $T$ =1.5 s， $b$ =2 m/s²， $l$ =4.5 m.

4.2 交通流密度的影响

为研究道路上交通流密度对本文模型结果的影响，将CAV和CV的渗透率分别设定为8%和20%. 首先，可基于仿真试验获取路段的流量-密度-速度关系如图5所示.

图 5 流量-密度-速度关系

Figure 5. Relationship between volume，density and speed of traffic flow

下载: 全尺寸图片幻灯片

由图5可知：仿真路段中的最佳密度为40 辆/km，即当交通流密度大于40 辆/km时，路段上的车辆开始出现拥挤；当交通流密度大于70 辆/km时，车辆的平均速度低于20 km/h，而高速公路中很少出现车辆速度小于20 km/h的情况. 因此，为研究本文模型在不同交通状态下的重构效果，此处将交通流密度设为20、30、40、50、60、70 辆/km. 为避免随机因素的影响，同一交通流密度下采用不同的随机种子进行仿真和轨迹重构，取其平均值作为最终结果. 轨迹重构图像示例如图6所示. 不同交通流密度条件下求得的插入车辆数如表1所示.

图 6 交通流密度为60辆/km时轨迹时空图

Figure 6. Spatio-temporal diagram of the trajectories when traffic flow density is 60 veh/km

下载: 全尺寸图片幻灯片

由表1可知：当交通流密度小于等于40 辆/km时，会出现插入的车辆数小于实际车辆数；而当交通流密度大于40 辆/km时，则出现插入的车辆数大于实际车辆数. 从百分比误差来看，当出现拥堵前，随着交通流密度的增加，模型求得的插入车辆数的误差逐渐减小，插入辆数的百分比误差也逐渐减小；当出现拥堵后，随着交通流密度的增加，模型求得的插入车辆数的误差逐渐增加，插入辆数的百分比误差也逐渐增大. 分析可知，开始拥堵之前，大多数车辆自由行驶，且车辆之间相互干扰较少，所以误差较大. 开始拥堵后，当两辆车之间插入的车辆数较大时，模型求得的插入车辆数会产生一定的误差. 随着交通流密度的增大，车辆之间的间距逐渐减小，出现两辆车之间插入的车辆数较大的情况增多，故插入辆数的误差逐渐增大.

表 1 不同交通流密度条件下插入车辆数估计

Table 1. Estimation of the numbers of inserted vehicles under different traffic densities

交通流密度/ （辆•km⁻¹）	CAV/ 辆	插入车辆/辆	实际车辆/辆	插入误差/辆	插入误差百分比/%
20	2	5.80	6.40	0.60	9.38
30	2	14.05	13.80	0.35	2.54
40	3	15.80	16.05	0.25	1.56
50	4	21.75	21.38	0.31	1.50
60	5	27.85	25.75	2.10	8.16
70	6	31.40	27.25	4.30	15.78
注：表中插入车辆数为多次仿真试验的平均值，因此会出现不是整数的情况；插入误差百分比为插入辆数的误差与实际插入车辆数的百分比.

下载: 导出CSV

| 显示表格

当交通流密度小于60 辆/km时，本文模型求得的插入车辆数的误差小于1辆，误差小于10%. 综上分析可知，本文模型求得的插入车辆数较为准确，适用于高速公路的各种交通状态.

为进一步评估本文模型的轨迹插入效果，选用每一时刻插入车辆的实际和估计位置之间的平均绝对误差（MAE，e_MAE）、平均绝对百分比误差（MAPE，e_MAE）和均方根误差（RMSE，e_MAE）表征重构轨迹的准确性，分别如式（9）～（11）所示.

${e_{\rm{MAE}}} = \frac{1}{k}\sum\limits_{t = 0}^{k\vartriangle t} {\left| {\hat x(t) - x(t)} \right|},$

(9)

${e_{\rm{MAPE}}} = \frac{1}{k}\sum\limits_{t = 0}^{k\Delta t} {\left| {\frac{{\hat x(t) - x(t)}}{{x(t)}}} \right|}\; \times \; {100{\text{%}} },$

(10)

${e_{\rm{RMSE}}} = \sqrt {\frac{1}{k}\sum\limits_{t = 0}^{k\Delta t} {{{\left[ {\hat x(t) - x(t)} \right]}^2}} },$

(11)

式中： $\hat x(t)$ 为插入车辆的估计位置； $x(t)$ 为插入车辆的实际位置.

根据式（9）～（11），可得不同交通流密度下的误差值，如表2所示. 由表2可知：当交通流密度为20 辆/km时，由于需要插入的车辆数过少，所以重构轨迹车辆位置的平均绝对误差相较其他密度较高. 当交通流密度大于20 辆/km时，重构轨迹车辆位置的误差随着交通流密度的增大而增加，但变动幅度不大. 在不同的交通流密度条件下，重构轨迹车辆位置的平均绝对误差的平均值为12.97 m，平均绝对百分比误差的平均值为0.37%，均方根误差的平均值为19.55 m. 从整体上来说，本文模型在不同交通流密度下重构轨迹车辆位置的误差较为稳定.

表 2 不同交通流密度条件下轨迹重构误差

Table 2. Trajectory reconstruction error under different traffic density conditions

交通流密度/ （辆•km⁻¹）	CAV/ 辆	MAE/ m	MAPE/ %	RMSE/ m
20	2	15.02	0.36	18.28
30	2	11.90	0.29	14.30
40	3	12.14	0.31	18.68
50	4	12.30	0.35	20.39
60	5	13.09	0.42	22.19
70	6	13.39	0.52	23.47
平均值		12.97	0.37	19.55

下载: 导出CSV

| 显示表格

4.3 CAV和CV渗透率的影响

CAV和CV提供的数据作为模型的已知数据，研究其渗透率对模型结果的影响尤为重要.

4.3.1 CAV渗透率的影响

为研究不同CAV渗透率对模型结果的影响，将交通流密度设为70 辆/km，CV渗透率设为20%. CAV的检测范围设为100 m. 若车辆在道路上均匀分布，则CAV检测范围包含7辆车. 但考虑到车辆的随机分布和感知设备的误差，此处假设CAV能记录检测范围内车辆数为6. 将CAV渗透率分别设置为4%、6%、8%、10%、12%和14%，研究不同渗透率下的车辆轨迹重构效果，其中，插入车辆数估计误差如表3所示，轨迹重构结果示例如图7所示.

表 3 不同渗透率条件下插入车辆数估计

Table 3. Estimation of the numbers of inserted vehicles under different CAV penetration rates

CAV 渗透率/%	CAV/辆	插入车辆/辆	实际车辆/辆	插入误差/辆	插入误差百分比/%
4	3	45.50	38.70	6.80	17.57
6	4	41.50	35.20	5.80	16.48
8	6	31.40	27.25	4.30	15.78
10	7	29.44	25.44	4.00	15.72
12	8	26.80	23.20	3.60	15.52
14	10	19.60	17.00	2.60	15.29

下载: 导出CSV

| 显示表格

图 7 CAV渗透率为14%时轨迹时空图

Figure 7. Spatio-temporal diagram of the trajectories when CAV penetration rate is 14%

下载: 全尺寸图片幻灯片

由表3可知：当道路上的交通流密度为70 辆/km时，随着CAV渗透率的增加，需要插入的车辆数减少，轨迹重构的难度减小，求出的插入车辆数的误差和百分比误差均逐渐减小；当渗透率只有4%的情况下，模型也能够较好地估计插入轨迹的数量. 因此，本文模型在未来不同渗透率环境下均能很好应用. 进一步，可计算重构轨迹车辆位置误差如表4所示.

表 4 CAV不同渗透率条件下轨迹重构误差

Table 4. Trajectory reconstruction error under different CAV penetration rates

CAV渗透率/%	CAV/辆	MAE/m	MAPE/%	RMSE/m
4	3	16.76	0.68	33.77
6	4	15.69	0.63	28.95
8	6	13.39	0.52	23.47
10	7	12.98	0.51	21.86
12	8	12.23	0.49	19.90
14	10	10.45	0.42	16.45

下载: 导出CSV

| 显示表格

由表4可知：当道路上的交通流密度为70 辆/km时，随着CAV渗透率的增加，重构轨迹车辆位置的误差逐渐减小，即本文模型的准确性逐渐提高. 以MAE为例，当CAV的渗透率在4%～8%范围内时，平均每增加一辆CAV，平均绝对误差减小1.12 m. 当CAV的渗透率在8%～14%范围内时，平均每增加一辆CAV，平均绝对误差减小0.73 m. 因此，当CAV渗透率大于8%时，平均绝对误差趋于稳定.

4.3.2 CV渗透率的影响

将交通流密度设为70 辆/km，CAV渗透率设为8%，研究CV渗透率为16%、18%、20%、22%、24%时轨迹重构结果，如表5所示.

表 5 CV不同渗透率轨迹重构结果

Table 5. Reconstruction results of different CV penetration rates

CV 渗透率/%	CV/辆	插入误差/辆	插入误差百分比/%	MAE/m	MAPE/ %	RMSE/ m
16	11	5.40	0.18	17.46	0.69	32.43
18	13	4.33	0.16	13.98	0.55	24.27
20	14	4.30	0.16	13.39	0.52	23.47
22	15	4.20	0.15	12.93	0.52	21.50
24	17	3.80	0.14	11.50	0.48	19.11

下载: 导出CSV

| 显示表格

由表5可知：随着CV渗透率的提高，插入辆数的误差逐渐减小，插入数量的百分比误差也逐渐减小，即轨迹重构效果不断提升. 此外，车辆位置的误差也越小，即轨迹重构的准确性有所提高.

为研究CAV和CV渗透率对试验结果的影响，以CAV渗透率为横轴，CV渗透率为纵轴，绘制均方根误差热力图，如图8所示.

图 8 均方根误差热力图

Figure 8. Heat map of RMSE

下载: 全尺寸图片幻灯片

由图8可知：均方根误差趋势为随着CAV或CV渗透率的提高，轨迹重构效果逐渐变好. 但与CV渗透率相比，CAV渗透率对结果的影响更大. 分析可知，增加一辆CAV能够获得道路上多辆车的轨迹数据，而增加一辆CV，只能获取该车的轨迹数据. 因此，当CAV渗透率较高时，已知轨迹的车辆较多，此时CV对试验结果的影响较小.

5. 结　论

1）在出现拥堵前，随着交通流密度的增加，本文模型求得的插入车辆数的误差逐渐减小；而在出现拥堵后，随着交通流密度的增加，本文模型求得的插入车辆数的误差逐渐增大.

2）本文模型重构轨迹车辆位置的平均绝对误差在不同的交通流密度条件下均在可接受范围内，即本文模型在不同的交通流密度下均能较为准确地重构所有车辆的轨迹.

3）重构轨迹的准确性随着CAV和CV渗透率的增加而提高. 与CV渗透率相比，CAV渗透率对结果的影响更大. 且当CAV渗透率较低时，本文模型也能较好地重构所有车辆的轨迹.

本文证明了基于跟驰特性的车辆轨迹重构模型的有效性，但选取的场景简单，仅研究了单车道高速公路车辆轨迹重构模型，未来将结合换道模型，研究智能网联环境下多车道车辆轨迹重构模型.

致谢：西南交通大学综合交通大数据应用技术国家工程实验室交大数科创新中心项目（JDSKCXZX202003）和重庆交通大学重庆市交通运输工程重点实验室开放课题（2018TE01）的资助.

图 1 3种混凝土的pH值变化（0、3、14、28 d）

Figure 1. pH changes of three types of concrete （0、3、14、28 d）

下载: 全尺寸图片幻灯片

图 2 碳化对混凝土力学性能的影响

Figure 2. Effect of carbonation on mechanical properties of concrete

下载: 全尺寸图片幻灯片

图 3 地聚物凝胶EDS谱图

Figure 3. EDS spectra of geopolymer gel

下载: 全尺寸图片幻灯片

图 4 MIP分析结果

Figure 4. MIP analysis results

下载: 全尺寸图片幻灯片

图 5 GPC-50与GPC-10的碳化模型

Figure 5. Carbonation models of GPC-50 and GPC-10

下载: 全尺寸图片幻灯片

表 1 粉煤灰、矿渣及水泥化学成分

Table 1. Compositions ratios of fly ash，slag and cement %

类型	SiO₂	Al₂O₃	Fe₂O₃	CaO	P₂O₅	Na₂O	K₂O	MnO	MgO	SO₃	Tio₂	SrO
粉煤灰	62.04	25.5	4.28	3.96	0.31	0.46			1.27	0.73	1.33	0.12
矿渣	34.11	15.36	0.83	35.99		0.40	0.62	1.07	6.58	2.50	2.41	0.12
水泥	20.01	5.68	3.20	65.89	0.08	0.09	1.17	0.19	0.81	2.00

下载: 导出CSV

表 2 混凝土配合比设计

Table 2. Mix proportion design and numbering of the concrete specimens kg/m³

编号	粉煤灰	炉渣	水泥	沙	砾石	Na₂SiO₃	NaOH	水	减水剂	养护条件
GPC-10	346.7	38.5		601.7	1203.5	165.7	66.2			80 ℃ sealed
GPC-50	192.6	192.6		601.7	1203.5	165.7	66.2			标准
OPCC	56.0	56.0	448.0	626.0	1022.0			174.7	2.5	标准

下载: 导出CSV

表 3 碳化对混凝土力学性能的影响

Table 3. Effect of carbonation on mechanical properties of concrete

编号	3 d	7 d	14 d	21 d	28 d
GPC-50
GPC-10
OPCC

下载: 导出CSV

表 4 原子相对数量

Table 4. Relative number of atoms %

编号	O	C	Si	Al	Na	Ca	Ca/Na	Si/Al
GPC-10	59.80	18.10	11.80	3.30	3.90	3.10	0.79	3.60
GPC-50	54.10	29.40	5.80	6.50	0.70	3.40	4.80	0.89

下载: 导出CSV

参考文献(30)

[1]	阿列克谢耶夫. 钢筋混凝土结构中钢筋腐蚀与保护[M]. 黄可信, 吴兴祖, 蒋仁敏, 等译. 北京: 中国建筑工业出版社, 1983: 19-35.
[2]	金伟良, 赵羽习. 混凝土结构耐久性[M]. 北京: 科学出版社, 2002: 16-25.
[3]	MARQUES P F, CARLOS C, NUNES A. Carbonation service life modelling of RC structures for concrete with Portland and blended cements[J]. Cement and Concrete Composites, 2013, 37(3): 171-184.
[4]	JUENGER M C G, WINNEFELD F, PROVIS J L, et al. Advances in alternative cementitious binders[J]. Cement and Concrete Research, 2011, 41(12): 1232-1243. doi: 10.1016/j.cemconres.2010.11.012
[5]	SHI C J, FERNÁNDEZ-JIMÉNEZ A, PALOMO A. New cements for the 21st century:the pursuit of an alternative to Portland cement[J]. Cement and Concrete Research, 2011, 41(7): 750-763. doi: 10.1016/j.cemconres.2011.03.016
[6]	TENNAKOON C, SHAYAN A, SANJAYAN J G, et al. Chloride ingress and steel corrosion in geopolymer concrete based on long term tests[J]. Materials and Design, 2017, 116: 287-299. doi: 10.1016/j.matdes.2016.12.030
[7]	BAKHAREV T. Resistance of geopolymer materials to acid attack[J]. Cement and Concrete Research, 2005, 32(4): 658-670.
[8]	BAKHAREV T. Durability of geopolymer materials in sodium and magnesium sulfate solutions[J]. Cement and Concrete Research, 2005, 35(6): 1233-1246. doi: 10.1016/j.cemconres.2004.09.002
[9]	ZHUANG X Y, CHEN L, KOMARNENI S, et al. Fly ash-based geopolymer:clean production,properties and applications[J]. Journal of Cleaner Production, 2016, 125: 253-267. doi: 10.1016/j.jclepro.2016.03.019
[10]	ZHANG J, SHI C J, ZHANG Z H, et al. Durability of alkali-activated materials in aggressive environments:a review on recent studies[J]. Construction and Building Materials, 2017, 152(2): 598-613. doi: 10.1016/j.conbuildmat.2017.07.027
[11]	BERNAL S A, PROVIS J L, WALKLEY B, et al. Gel nanostructure in alkali-activated binders based on slag and fly ash,and effects of accelerated carbonation[J]. Cement and Concrete Research, 2013, 53: 127-144. doi: 10.1016/j.cemconres.2013.06.007
[12]	LAW D W, ADAM A A, MOLYNEAUX T K, et al. Long term durability properties of class F fly ash geopolymer concrete[J]. Materials and Structures, 2015, 48(3): 721-731. doi: 10.1617/s11527-014-0268-9
[13]	BAKHAREV T, SANJAYAN J G, CHENG Y B. Resistance of alkali-activated slag concrete to carbonation[J]. Cement and Concrete Research, 2001, 31(9): 1277-1283. doi: 10.1016/S0008-8846(01)00574-9
[14]	BADAR M S, KUPWADE-PATIL K, BERNAL S A, et al. Corrosion of steel bars induced by accelerated carbonation in low and high calcium fly ash geopolymer concretes[J]. Construction and Building Materials, 2014, 61: 79-89. doi: 10.1016/j.conbuildmat.2014.03.015
[15]	PASUPATHY K, BERNDT M, CASTEL A, et al. Carbonation of a blended slag-fly ash geopolymer concrete in field conditions after 8 years[J]. Construction and Building Materials, 2016, 125: 661-669. doi: 10.1016/j.conbuildmat.2016.08.078
[16]	CRIADO M, PALOMO A, FERNÁNDEZ-JIMÉNEZ A. Alkali activation of fly ashes. part 1:effect of curing conditions on the carbonation of the reaction products[J]. Fuel, 2005, 84(16): 2048-2054. doi: 10.1016/j.fuel.2005.03.030
[17]	黄琪,石宵爽,王清远,等. 粉煤灰基地聚物混凝土的碳化性能研究[J]. 中国农村水利水电,2015(7): 121-125,130. doi: 10.3969/j.issn.1007-2284.2015.07.031 HUANG Qi, SHI Xiaoshuang, WANG Qingyuan, et al. Research on carbonation of fly ash geopolymeric concrete[J]. China Rural Water and Hydropower, 2015(7): 121-125,130. doi: 10.3969/j.issn.1007-2284.2015.07.031
[18]	黄琪,石宵爽,王清远,等. 再生粗骨料对粉煤灰基地聚物混凝土碳化性能的影响[J]. 硅酸盐通报,2015,34(5): 1264-1269,1281. HUANG Qi, SHI Xiaoshuang, WANG Qingyuan, et al. Effect of recycled coarse aggregate on carbonation resistance of fly ash geopolymeric concrete[J]. Bulletin of the Chinese Ceramic Society, 2015, 34(5): 1264-1269,1281.
[19]	陈晓星,曹海琳,翁履谦,等. 碱激发水泥砂浆碳化行为研究[J]. 武汉理工大学学报,2014,36(3): 18-22. CHEN Xiaoxing, CAO Hailin, WENG Lvqian, et al. Research on carbonation process of alkali-activated cement mortars[J]. Journal of Wuhan University of Technology, 2014, 36(3): 18-22.
[20]	贺鹏飞. 混凝土碳化研究进展[C]//《工业建筑》2018年全国学术年会论文集（上册）. 北京: 工业建筑杂志社, 2018: 316-321.
[21]	KUMAR S S, VASUGI J, AMBILY P S, et al. Development and determination of mechanical properties of fly ash and slag blended geopolymer concrete[J]. International Journal of Scientific and Engineering Research, 2013, 4(8): 1-5.
[22]	NATH P, SARKER P K. Effect of GGBFS on setting,workability and early strength properties of fly ash geopolymer concrete cured in ambient condition[J]. Construction and Building Materials, 2014, 66: 163-171. doi: 10.1016/j.conbuildmat.2014.05.080
[23]	杨世玉,赵人达,靳贺松,等. 单组分地聚物砂浆的力学性能和微观结构分析[J]. 西南交通大学学报,2021,56(1): 101-107. YANG Shiyu, ZHAO Renda, JIN Hesong, et al. Mechanical performance and microstructure of single component geopolymer mortar[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 101-107.
[24]	ZHANG H E, SHI X S, WANG Q Y. Effect of curing condition on compressive strength of fly ash geopolymer concrete[J]. ACI Materials Journal, 2018, 115(2): 191-196.
[25]	中华人民共和国住房和城乡建设部. 普通混凝土长期性能和耐久性能试验方法标准: GB/T 50082—2009[S]. 北京: 中国建筑工业出版社, 2009.
[26]	中华人民共和国建设部. 普通混凝土力学性能试验方法标准: GB/T 50081—2002[S]. 北京: 中国建筑工业出版社, 2003.
[27]	BAKHAREV T. Geopolymeric materials prepared using class F fly ash and elevated temperature curing[J]. Cement and Concrete Research, 2004, 35(6): 1224-1232.
[28]	ZHAO R D, YUAN Y, CHENG Z Q, et al. Freeze-thaw resistance of class F fly ash-based geopolymer concrete[J]. Construction and Building Materials, 2019, 222(3): 474-483. doi: 10.1016/j.conbuildmat.2019.06.166
[29]	杨绿峰,成荻,刘才勇,等. 矿物掺合料混凝土碳化分析的实用预测模型[J]. 混凝土,2016(7): 79-83. doi: 10.3969/j.issn.1002-3550.2016.07.019 YANG Lufeng, CHENG Di, LIU Caiyong, et al. Practical prediction model of carbonation depth for concrete with mineral admixtures[J]. Concrete, 2016(7): 79-83. doi: 10.3969/j.issn.1002-3550.2016.07.019
[30]	SAETTA A V, VITALIANI R V. Experimental investigation and numerical modeling of carbonation process in reinforced concrete structures part I:theoretical formulation[J]. Cement and Concrete Research, 2004, 34(4): 571-579. doi: 10.1016/j.cemconres.2003.09.009

施引文献

附加材料(0)

访问统计

点击查看大图

图(5) / 表(4)

计量

文章访问数: 667
HTML全文浏览量: 263
PDF下载量: 22
被引次数: 0

1. 跟驰特性分析
2. 基于IDM的车辆状态估计
2.1 数据环境分析
2.2 估计插入车辆的速度
2.3 估计插入车辆的位置
3. 基于IDM的车辆轨迹重构
4. 数值仿真分析
4.1 仿真参数设置
4.2 交通流密度的影响
4.3 CAV和CV渗透率的影响
5. 结　论

粉煤灰-矿渣基地聚物混凝土的抗碳化性能

doi: 10.3969/j.issn.0258-2724.20191151

作者简介: 原元（1992—），男，博士研究生，研究方向为高使用性混凝土材料性能及结构应用，E-mail：1181491115@qq.com

通讯作者: 李福海（1979—），男，高级工程师，博士，研究方向为新型材料及混凝土结构耐久性，E-mail：lifuhai2007@home.swjtu.edu.cn

计量

出版历程

Carbonation Resistance of Fly Ash-Slag Based Geopolymer Concrete

1. 跟驰特性分析

2. 基于IDM的车辆状态估计

2.1 数据环境分析

2.2 估计插入车辆的速度

2.3 估计插入车辆的位置

3. 基于IDM的车辆轨迹重构

4. 数值仿真分析

4.1 仿真参数设置

4.2 交通流密度的影响

4.3 CAV和CV渗透率的影响

4.3.1 CAV渗透率的影响

4.3.2 CV渗透率的影响

5. 结 论

计量

出版历程

目录

1. 跟驰特性分析

2. 基于IDM的车辆状态估计

2.1 数据环境分析

2.2 估计插入车辆的速度

2.3 估计插入车辆的位置

3. 基于IDM的车辆轨迹重构

4. 数值仿真分析

4.1 仿真参数设置

4.2 交通流密度的影响

4.3 CAV和CV渗透率的影响

5. 结 论

作者简介:
原元（1992—），男，博士研究生，研究方向为高使用性混凝土材料性能及结构应用，E-mail：1181491115@qq.com

通讯作者:
李福海（1979—），男，高级工程师，博士，研究方向为新型材料及混凝土结构耐久性，E-mail：lifuhai2007@home.swjtu.edu.cn

5. 结　论

5. 结　论