我国高速铁路事业近年来发展迅猛，在高速列车的长期高频冲击下，铁路基础设施的服役性能在不断发生变化^[1-4]. 因此，需要对铁路系统的主要结构组件（包括钢轨、扣件、轨枕、无砟道床等）的结构健康状态进行长期持续监测^[5]. 轨枕作为将车轮荷载从钢轨传递到道床的关键结构，若能嵌入传感功能，则可提供列车荷载作用下的枕上压力等关键信息，其在一定程度上反映了高速列车载荷下的轮轨力状态.

目前，对轨道结构的监测主要分为车载监测和地面监测^[6]. 金鑫等^[7]分析了轮轨作用力对轮轴结构受力特性的影响，提出了在车辆轮轴上粘贴传感器的轮轨力测试方法. 宫雪^[8]结合轮辐测力与车轴测力的特点，提出了辐轴结合法，并在模拟滚动台上进行了验证. 上述车载测试方法，尽管能够达到较高的测量精度，但均需制备特殊的测力轮对，测试成本较高，且只能用于周期性的轮轨作用关系检测，难以满足全天候安全服役状态监测的需要. 因此，进行地面监测也十分必要.

Zhang等^[9]在轨下胶垫中植入了压电压力传感器，在实验室内完成标定后，将其布置到运营线路中进行现场测试，并获取了列车通过时轨枕承受的压力大小. Sadeghi^[10]利用压电压力传感器研究了钢轨和轨枕、轨枕和道床之间的相互作用力及其线性关系. 宋颖等^[11]提出了基于压电薄膜的轮轨力测试方法，并通过有限元法计算验证了监测原理的可行性. 由于电信号传感元件在防水、抗电磁干扰、耐高温、耐腐蚀等方面性能的不足，在长期监测过程中会不可避免地产生损坏，其监测稳定性和可靠性大大降低.

光纤光栅传感器是以光波为载体，光纤为媒介，通过捕捉反射光波的变化特征，实现对其相关物理量的测试^[12]. 近年来，光纤光栅在大型工程安全监测领域得到了迅速发展和广泛应用^[13-16]. 孟汇^[17]研究了光纤光栅传感器在复杂环境中的使用年限，预测寿命在20年左右. 在高速铁路运营线路的复杂工况下，光纤光栅传感器监测的稳定性以及耐久性具有较大优势. 除此之外，光纤光栅传感器尺寸较小、有良好的动态特性，其测试精度与测试量程满足高速铁路监测需求.

国内外学者对光纤光栅传感器在铁路领域的应用展开了大量研究. Wei等^[18]在钢轨垂向、纵向安装布拉格光栅（FBG）光纤传感器，根据实验列车以不同速度经过时钢轨的应变响应，找到了最适合的布置位置，然后在运营线路钢轨上安装光纤传感器测量了钢轨的应变变化. 张兆亭等^[19]使用光纤光栅传感器进行了动态载荷下的钢轨应变监测实验，结果表明，光纤光栅应变传感器适用于钢轨应变的监测需要，具有良好的工作性能. 高亮等^[20]提出了一种基于光纤光栅的轮轨力长期监测方法，通过在钢轨跨中断面的3处不同部位粘贴光纤光栅传感器，采集列车通过时的钢轨修正应变，结合现场标定的轮轨力-钢轨应变线性关系，可长期实时获取轮轨相互作用力. 但上述研究均基于表面接触式监测，由于传感器及其导线外露，在长期监测中有脱落的风险.

因此，可通过在轨道结构的混凝土部件内植入传感器作为监测方法的补充手段，Butler等 ^[21-23]研究了有砟铁路道床、轨枕状态的监测方法，在预应力混凝土轨枕的制造阶段将光纤光栅准分布式传感系统埋入其中，并初步讨论了自感知轨枕用于铁路道床荷载传递监测的可行性. 上述研究是基于有砟轨道的预应力轨枕，目前对高速铁路无砟轨道轨枕、承轨台位置荷载的传递监测还处于空白.

鉴于现有研究的不足，本文以CRTS Ⅰ型双块式轨枕为研究对象，首次提出了一种双块式无砟轨道枕上压力监测标定方法. 将光纤光栅传感器埋入双块式轨枕中，通过反力架静载试验对传感器进行标定，并通过有限元仿真进行修正. 采用该技术可实现对轨枕枕上压力的长期监测，为轮轨力监测提供一种间接监测方法，也可为研究高速铁路轮轨荷载传递、完善无砟轨道结构的强度计算理论和方法提供重要依据.

1.   双块式轨枕及其传感器布置

1.1   光纤光栅传感器及其测试原理

本文试验选取如图1所示的埋入式光纤光栅FRP-OF （fiber reinforced polymer-optical fiber）应变传感器，其尺寸小，埋入混凝土后不影响结构的应力状态与应力分布.

图  1  FRP-OF应变传感器

Figure  1.  FRP-OF strain sensor

下载: 全尺寸图片 幻灯片

该传感器采用无胶封装，克服了传统胶封光纤光栅传感器的耐久性问题，并且在实际工程中布设简单，量程大、精度高，动态特性良好. 具体参数见表1.

表  1  FRP-OF应变传感器技术参数

Table  1.  Technical parameters of FRP-OF strain sensor

传感器技术参数	数值
量程/με	−1500～3 000
分辨率/με	±0.1
测量精度/%	0.5
应变系数/ （pm·με）	0.8～1.2
线性度/%	99.9
重复性误差/%	≤0.5
中心波长/nm	1528～1568
反射率/%	≥90
标距尺寸/mm	60
工作温度范围/℃	−30～80

下载: 导出CSV 

| 显示表格

FRP-OF应变传感器的结构示意及其工作原理如图2所示. 传感器端部采用锚头限位，使其与被测混凝土协同变形，布拉格光栅（FBG）被封装在纤维增强复合材料（FRP）中，通过光纤（OF）传导光信号. 当入射光经过布拉格光栅时，满足相位匹配条件的光会被反射回来，不满足布拉格相位匹配的则会被透射出去. 利用光纤光栅对压力的敏感特性，制成了FRP-OF应变传感器.

图  2  FRP-OF应变传感器结构及其工作原理

Figure  2.  Structure of FRP-OF strain sensor and its working principle

下载: 全尺寸图片 幻灯片

当满足相位匹配条件时，光纤的谐波长如式（1）所示，其中任何一个参数发生变化都会引起光纤谐振波长发生偏移，由此可得布拉格光栅谐振波长偏移量如式（2）所示，当外界应变或者温度发生变化时，引起光纤光栅传播有效折射率n_e和光栅周期T的变化，通过光纤光栅解调仪可以检测出谐振波长的偏移量，得到输出波长，如式（3）所示，通过计算波长变化量与出厂时传感器的应变标定系数K_b的比值可以计算出相应的应变变化.

式中：ε为测试应变值；λ为输出波长；λ₀为初始波长；$\Delta {{{n}}_{\rm{e}}}$为光纤光栅受到外界应力时的有效折射率变化；$\Delta T$为光纤光栅受到外界应力时的光栅周期变化.

1.2   传感器埋设

本文试验所用的双块式轨枕由某轨枕厂制造. 图3为传感器植入位置示意，其中传感器中心距离承轨台表面约为80 mm，距离螺栓孔116.5 mm，位于承轨台正下方.

图  3  双块式轨枕结构以及传感器埋设位置示意

Figure  3.  Bi-block sleeper structure and sensor embedment position

下载: 全尺寸图片 幻灯片

如图4所示，在轨枕模具中通过轧带垂向固定FRP-OF应变传感器，用于感知当轨枕表面受到压力时轨枕内部的垂向应变，传感器导线沿着桁架筋从轨枕底部边缘引出.

图  4  双块式轨枕生产阶段安装传感器示意

Figure  4.  Schematic diagram of installing sensor in production stage of bi-block sleeper

下载: 全尺寸图片 幻灯片

在传感器安装完毕后向模具中浇筑C60混凝土，在振捣台上将混凝土振捣均匀后送入养护室进行蒸汽养护，待其蒸养24 h完毕后进行脱模，将制作好的4根双块式轨枕的8块轨枕块编号为1#～8#，随后进行标定试验.

2.   标定试验

2.1   试验方法

由于SK-Ⅱ型双块式轨枕下部有桁架筋结构，无法将其直接放在压力试验机上进行静载标定. 为测试轨枕内传感器的感知能力以及找到枕上压力-枕内应变的线性关系，本文设计了一套反力架工装，如图5所示. 安装轨下胶垫、铁垫板、侧向挡块以及反力架，将千斤顶放置于铁垫板和反力架钢板之间，并且在千斤顶上安装一个压电压力传感器用于记录加载数值，构成了一套自反力加载系统. 通过千斤顶对轨枕进行加载，而反力架钢板作为反力点，铁垫板和轨下胶垫均匀地将荷载均匀传递到轨枕表面.

图  5  反力架工装结构示意

Figure  5.  Schematic diagram of structure of reaction frame

下载: 全尺寸图片 幻灯片

如图6所示，使用DH5983便携式动态信号分析仪采集枕上压力数据，使用光纤光栅解调仪采集埋入轨枕的FRP-OF应变传感器的波长变化，采样频率均为3 Hz. 进行预加载，使结构进入正常的工作状态，使反力架和轨枕之间贴合紧密，预加载为10 kN. 加载过程中观察光纤光栅应变计波长变化情况以及卸载后波长是否恢复到加载前的状态，并且观察放置于轨枕上方的压电传感器数值是否稳定，准备工作完成后开始正式加载.

图  6  通过反力架-千斤顶进行标定试验

Figure  6.  Calibration test with reaction frame and jack

下载: 全尺寸图片 幻灯片

在正式加载中使用千斤顶将载荷缓慢增加到20 kN，待其稳定约30 s后卸载至0，每个轨枕块经过3次加载-卸载循环. 高速列车轴重一般为150 kN，轮载为75 kN，根据文献[24]轮对正下方轨枕的轮载分配率约为33%，枕上压力约为24.75 kN. 混凝土单轴受压时的本构模型在前1000 με的应力-应变曲线为线性^[25]，故轨枕在服役过程中实际受力超过20 kN时的应力应变也为线性变化的.

2.2   试验结果

通过千斤顶对轨枕施加荷载，如图7（a）所示，轨枕内的FRP-OF应变传感器对枕上压力具有良好的感知性能，所测得的应变与载荷同步变化. 当载荷增加到20 kN时，轨枕在3个加载周期的应变分别达到了3.8、3.5 με和3.7 με，由于载荷迅速减小到0，轨枕内应变传感器的应变也急剧减小到0. 图7（b）为3次加载试验枕内应变-枕上压力的曲线，3次加载周期的曲线重合程度较高，线性度良好. 图7（c）为第1次加载时1#～8# 号轨枕块的枕内应变-枕上压力曲线，不同轨枕块之间的线性变化趋势一致，但斜率略有不同. 如图7（d）所示，为了得到枕上压力和枕内应变的线性关系斜率，对采样点进行线性拟合，得到线性回归曲线方程，如式（4）所示. 线性回归曲线的损失函数J（K, b）如式（5）所示（K为斜率，b为标定截距），用损失函数来衡量回归模型的误差，这个损失函数的值越小，直线能更好地拟合实际数据. 因此，对J（K, b）的K求导，使其导数为0即可得到K的最优解. 线性回归曲线的斜率K可由式（6）计算. 1#～8# 号轨枕块在反力架加载时的枕内应变-枕上压力斜率约为4.90～5.28，如表2所示.

图  7  标定试验数据

Figure  7.  Calibration test data

下载: 全尺寸图片 幻灯片

基于反力架的静载试验初步表明，植入FRP-OF应变传感器在轨枕中具有良好的传感性能，并初步得到每块轨枕的标定曲线及标定系数，通过标定系数($\bar K$)，为斜率的平均值和截距可以得到轨枕枕上压力的计算如式（7）.

式中：F为枕上压力；ε₁₁为标定试验枕内传感器测试应变；n为采样点数量；F_i为枕上压力的第i个数据样本点；ɛ_i为枕内应变的第i个数据样本点；$\bar{F}$为所有样本的平均枕上压力；$\bar{\mathrm{\varepsilon }}$为所有样本的平均枕内应变.

表  2  标定试验的荷载值为20 kN时的数据

Table  2.  Data obtained at load of 20 kN in calibration test

轨枕块 编号	枕内应变 平均值/με	枕内应变-枕 上压力斜率 平均值/ （kN·με−1）	枕内应变-枕 上压力截距 平均值/kN
1#	3.7	5.03	0.61
2#	3.6	4.98	1.11
3#	3.9	4.90	0.73
4#	3.9	4.93	0.76
5#	3.9	4.78	0.50
6#	3.7	5.10	0.69
7#	3.6	5.28	0.10
8#	3.7	5.21	1.14

下载: 导出CSV 

| 显示表格

3.   自感知性能测试试验仿真分析

通过标定试验初步验证了枕内应变-枕上压力具有良好的线性关系. 为进一步验证试验数据的正确性，通过有限元仿真模拟双块式轨枕在实际工况受力时的枕内应变，并对标定试验的枕上压力计算公式进行修正.

3.1   模型建立

如图8所示，采用有限元软件建立双块式轨枕、道床板的有限元模型，在轨枕的内部预留一个和光纤光栅传感器尺寸相当的孔，随后将传感器装配进预留的孔里与轨枕块进行约束.

图  8  轨枕块和道床板有限元模型

Figure  8.  Finite element model of sleeper block and track slab

下载: 全尺寸图片 幻灯片

轨枕、道床板和传感器的模型材料参数如表3所示.

表  3  有限元模型材料参数

Table  3.  Material parameters of finite element model

部件	材料	弹性 模量/GPa	泊松比	密度/ (kg·m−3)
轨枕块	C60 混凝土	36.5	0.20	2640
道床板	C40 混凝土	32.5	0.24	2440
传感器	FRP	72.0	0.20

下载: 导出CSV 

| 显示表格

反力架静力加载标定试验的边界条件如图9（a）所示，将轨枕底部的反力架约束区域进行完全固定，并且对其轨枕表面施加分布荷载；轨枕在道床板中实际的边界条件如图9（b）所示，将轨枕与道床板进行绑定约束，并且将道床板固定在底座板上，对轨枕块表面施加均布荷载.

图  9  不同工况的边界条件

Figure  9.  Boundary conditions of different working conditions

下载: 全尺寸图片 幻灯片

3.2   仿真计算分析

如图10所示，标定试验中1# 轨枕块的实际测量枕内应变-枕上压力线性拟合斜率为5.03，仿真模拟的斜率为4.99，误差为0.80%，且1#～8# 轨枕块与仿真数据的误差均在5.00%以内. 仿真数据和实测数据基本一致，由此说明模型建立正确并且现场试验所得数据可靠.

图  10  反力架标定试验数据与仿真数据对比

Figure  10.  Comparison of reaction frame calibration test data and simulation data

下载: 全尺寸图片 幻灯片

图11为轨枕在不同工况下施加20 kN枕上压力时的应变云图，其中，图11（a）为反力架加载工况下内部传感器埋入点应变大小约为3.77 με；图11（b）为轨枕在道床板实际工况下，内部传感器埋入点的应变大小约为3.40 με.

图  11  仿真模拟不同工况下的应变云图

Figure  11.  Strain nephogram at different working conditions

下载: 全尺寸图片 幻灯片

由于利用反力架进行标定时的边界条件与轨枕在道床板中的实际边界条件不同，通过有限元仿真模拟两种工况的传力区别，找到在两种工况下枕内应变的线性拟合曲线方程，如式（9）. 得到了两种工况下枕内应变的换算关系如图12所示.

图  12  反力架约束工况和实际工况枕内应变数值线性拟合

Figure  12.  Linear fitting of strain values under reaction frame constraint condition and actual condition

下载: 全尺寸图片 幻灯片

根据式（9）对反力架约束工况的标定试验进行修正，修正后的枕上压力换算公式如式（10）所示.

式中：ɛ₂₂为实际工况下的轨枕应变；α为标定工况与实际工况的修正系数，取1.35；b取0.1；$\beta$为标定工况与实际工况线性拟合曲线截距.

4.   现场测试分析

将埋入FRP-OF应变传感器的轨枕安装在我国某设计速度为350 km/h的高速铁路运营线路上，设置光纤光栅解调仪的采集频率为1 000 Hz，使用低通滤波器过滤高频干扰信号后得到如图13（a）所示的应变曲线，该数据来源于图4中的1# 号轨枕块的传感器. 当8车厢编组的CRH380AJ列车以350 km/h速度通过，枕内应变约为4～6 με. 如图13（b）所示，根据式（9）的标定系数将枕内应变换算为枕上压力，幅值约为30～42 kN，根据文献[24]所使用的测力垫板进行枕上压力测试，幅值范围为32～38 kN，说明本文提出的枕上压力测试方法是可信的，证明了该监测标定方法的使用价值.

图  13  高速铁路运营线路现场测试数据时程

Figure  13.  Time history diagram of field test data from high-speed railway operating line

下载: 全尺寸图片 幻灯片

5.   结　论

本文提出了一种利用埋入式应变传感器监测双块式轨枕枕上压力的新方法. 首次将FRP-OF应变传感器植入双块式轨枕中，进行了轨枕块反力架静力加载标定试验，通过有限元仿真计算进行验证，并修正了枕上压力计算公式，最后在现场试验证明了其可行性，具体得出了以下结论.

1） 提出了一种双块式轨枕枕上压力的监测新方法，并通过现场试验和理论计算验证了此方法的可行性. 为实现轨枕支点反力的长期监测、高速铁路轮轨荷载传递研究、无砟轨道结构的强度计算理论和方法完善、列车轮对状态监测提供了试验基础.

2） 通过将反力架安装在轨枕上进行静力加载试验，枕上压力和枕内应变之间的线性关系良好，其标定系数在4.90～5.28，说明了FRP-OF光纤光栅传感器植入轨枕内部具有较好的感知能力且重复性良好.

3） 通过仿真计算验证了反力架静力加载试验的测试数据误差在5.00%以内. 轨枕在道床板中实际工况的应变值略低于反力架加载工况，最后通过仿真模拟对反力架静载试验的标定试验进行修正，得出了枕上压力的计算方法.

4） 在350 km/h高速列车动载作用下测试到轨枕应变约为4～6 με，根据本文的标定方法计算出枕上压力幅值为30～42 kN，证明了该方法具有较高的使用价值.