Loading [MathJax]/extensions/TeX/boldsymbol.js
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

路面典型局部变形激励下的振动舒适性评价

郭雄峰 陈长 冯奇

郭雄峰, 陈长, 冯奇. 路面典型局部变形激励下的振动舒适性评价[J]. 西南交通大学学报, 2022, 57(2): 376-383. doi: 10.3969/j.issn.0258-2724.20200041
引用本文: 郭雄峰, 陈长, 冯奇. 路面典型局部变形激励下的振动舒适性评价[J]. 西南交通大学学报, 2022, 57(2): 376-383. doi: 10.3969/j.issn.0258-2724.20200041
GUO Xiongfeng, CHEN Zhang, FENG Qi. Vibration Serviceability Evaluation Under Typical Local Deformations of Pavement[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 376-383. doi: 10.3969/j.issn.0258-2724.20200041
Citation: GUO Xiongfeng, CHEN Zhang, FENG Qi. Vibration Serviceability Evaluation Under Typical Local Deformations of Pavement[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 376-383. doi: 10.3969/j.issn.0258-2724.20200041

路面典型局部变形激励下的振动舒适性评价

doi: 10.3969/j.issn.0258-2724.20200041
基金项目: 国家自然科学基金(71471134)
详细信息
    作者简介:

    郭雄峰(1992—),男,工程师,硕士,研究方向为路面养护与管理,E-mail:tjrommel@163.com

    通讯作者:

    陈长(1977—),男,副教授,博士,研究方向为交通基础设施养护与管理,E-mail:gxf@tongji.edu.cn

  • 中图分类号: U418

Vibration Serviceability Evaluation Under Typical Local Deformations of Pavement

  • 摘要:

    为研究沥青路面局部变形对行驶质量的影响,对路面养护管理中常见的错台型、双断面型和曲面型3类局部变形激励下的振动舒适性进行了评价. 首先,建立了局部变形的纵断面简化模型,并基于1/4车理论推导出路面局部变形激励下的人体竖向加速度时域表达式;其次,设计频率加权滤波器对竖向加速度进行滤波,得到瞬时频率加权加速度,从而求解出不同工况下均方根加速度,实现了振动舒适性的时域加权评价;最后,提出了局部变形的评价标准. 分析结果表明:均方根加速度与局部变形高差呈正比关系;同等变形程度下,双断面型变形的行驶感受最差;当错台型、双断面型及曲面型3类局部变形分别达到4.76、2.78、3.26 cm时,行驶质量均达到极不舒适状态.

     

  • 行驶舒适性与路面的平整程度、车辆状态、运行速度等方面密切相关,被广泛认为是路面服务性能评价的重要参考指标. 路面的不平整激励是引起汽车振动并导致舒适性下降的核心因素,道路表面横向、 纵向、水平向的各种波长及幅值的复杂变形组合构成了路面的不平度[1]. 作为激励源,不平度的连续变化会引起车辆的随机振动,局部突变则会引起车辆的瞬态振动,二者共同作用,通过车身悬架系统传递至人体,最终引起其心理和生理上的舒适性感受变化.

    国际上用于评估行驶舒适性的指标众多,例如舒适度指数(ride number,RN)、国际平整度指数(international roughness index,IRI)、行驶质量指数(riding quality index,RQI)以及加权加速度均方根(weighed root-mean-square acceleration,WRMSA)等. 利用这些评价指标,现有文献研究集中在宏观层面的路段行驶质量和人车振动特性进行分析[2-4],但对于路面局部变形(如桥头差异沉降、坑槽、窨井沉陷)导致的冲击振动问题却鲜有涉及.

    已有研究表明路面局部变形引起的冲击振动对行驶舒适性具有显著影响[5-6],但尚未明确其中的定量关系,部分文献也提出过路面中桥头差异沉降及窨井沉陷等局部变形的控制阈值[7-9],但是缺乏有效的理论支撑. 另有文献提出以较短长度作为行驶质量评价单元能够更加准确地描述路面的服务性能[10],但是较短的单元长度可能会造成路面波长截断现象[2],导致评价结果失真. 由于路面局部变形随机存在于路面中,且其导致的振动具有瞬时性特点,因此,难以通过现场试验获取变形激励与振动感受之间的定量关系.

    鉴于此,本文从理论上研究局部变形对舒适性的影响. 首先,建立典型路面局部变形简化模型作为振动激励源,基于两自由度车辆强迫振动模型,从理论上推导出振动模型中簧上质量位移和加速度时间历程函数,再根据ISO 2631-1[11]中的频率加权系数分布特点设计频率加权滤波器对竖向加速度进行滤波,求解不同工况下的均方根加速度,评估人体的振动感受,最后提出不同类型局部变形的评价标准,以期为路面精细化养护管理提供依据.

    文献[2, 12]对上海城市道路进行详细调研发现,引起道路纵断面高程局部异常波动的不平整因素主要是窨井变形、桥头差异沉降、路面修补等,各类因素的统计结果如图1所示.

    图  1  路面局部变形因素统计
    Figure  1.  Statistical factors of local deformation in pavement

    图1可知,窨井损坏和修补(包括路面开挖回填)是最主要的变形形式,其次是桥头差异沉降. 但是,实际路面中路面修补一般面积较大,且对行车影响十分有限;路面开挖则一般为临时性养护管理;另外,由于轮胎存在包络特性,伸缩缝问题也不具显著影响. 为进一步研究变形的影响,本文根据其余常见局部变形的形态特征,将其分为3类:

    1) 错台型,描述同一断面的高程突变,如桥头差异沉降;

    2) 双断面型,描述路面规则的凹凸变形,行进中的轮胎会在极短时间内经过变形区域的两端,如窨井沉陷或拱起;

    3) 曲面型,描述路面较均匀的凸起,如路面拥包等.

    建立激励模型需首先假定各类局部变形的函数形式,由于路面局部变形处纵断面高程一般具有突变特性,传统简谐函数难以描述其断面特点,因此对3类变形进行简化(如图2所示),以高差Δg作为变形衡量指标,纵断面高程y(t)如式(1)~(3)所示.

    1) 错台型

    y(t)={0,t=0,Δg,0<tl/v, (1)

    式中:l 为变形部分长度; t 为时间; v 为车速.

    2) 双断面型

    y(t)={0,t=0,Δg,0<td/v,0,t>d/v, (2)

    式中:d为凹凸部分直径或长度.

    3) 曲面型

    y(t)={0,t=0,4Δgvl2(vt2lt),0<tl/v,0,t>l/v. (3)
    图  2  典型路面局部变形纵断面示意
    Figure  2.  Profile diagrams of typical local deformation

    振动舒适性评价一般以加速度为基础指标,已有研究中加速度多为现场实测值,不同车辆之间难以对比. 本文通过两自由度车辆强迫振动模型的求解,从理论上推导人体的竖向位移和加速度响应.

    为统一各国的平整度度量标准,世界银行曾在巴西进行了大量平整度试验(international road roughness experiment),试验中首次构想出两自由度1/4理想车模型[13],由于理想车的运动方程恒定,固定车辆参数后,相当于利用同一辆车在任意地方行驶,因此,其测试方法具有很好的移植性. 1/4车模型如图3所示,图中:ms为簧上质量,mt为簧下质量,Ks为悬架刚度系数,Kt为轮胎刚度系数,Cs为悬架阻尼系数,ZsZt分别为簧上质量和簧下质量的绝对位移.

    图  3  1/4车模型
    Figure  3.  Quarter car model

    运动方程为

    msZs+Cs(ZsZt)+Ks(ZsZt)=0 (4)
    mtZt+Cs(ZtZs)+Ks(ZtZs)+Kt(Zty)=0. (5)

    将式(4)、(5)等式左右两侧同时除以ms,得

    Zs+c(ZsZt)+k2(ZsZt)=0 (6)
    uZt+c(ZtZs)+k2(ZtZs)+k1(Zty)=0 (7)

    式中:系数 c=Cs/ms=6.00s1 系数u=mt/ms=0.15;系数k1=Kt/ms=653s2系数k2=Ks/ms=63.3s2.

    路面局部变形对两自由度车辆振动系统的激励不具有简谐特性和周期特性,在这种激振作用下,系统通常仅有瞬态振动,故为求得加速度的时域函数,可以采用卷积积分法[14]求解簧上质量的加速度响应,具体过程如下:

    1) 求解1/4车振动系统对于单位脉冲激励δ(t)的响应函数;

    2) 利用卷积积分法求解路面局部变形激励作用下的系统响应.

    簧上质量的竖向位移可近似认为是人体的竖向位移. 由于运动方程存在耦合,故需要对方程进行解耦,运动方程如式(8)所示.

    {\boldsymbol{M}}\mathop {\boldsymbol{Z}}\limits^{ \bullet \bullet } + {\boldsymbol{C}}\mathop {\boldsymbol{Z}}\limits^ \bullet + {\boldsymbol{KZ}} = {\boldsymbol{Y}}, (8)

    式中:MCKZY分别为质量矩阵、阻尼矩阵、刚度矩阵、位移矩阵、输入矩阵,

    \begin{split}& {\boldsymbol{M}} = \left[ {\begin{array}{*{20}{c}} {{m_{\text{s}}}}&0 \\ 0&{{m_{\text{t}}}} \end{array}} \right],\;\;\; \\ & {\boldsymbol{C}} = \left[ {\begin{array}{*{20}{c}} {{C_{\text{s}}}}&{ - {C_{\text{s}}}} \\ { - {C_{\text{s}}}}&{{C_{\text{s}}}} \end{array}} \right], \\& {\boldsymbol{K}} = \left[ {\begin{array}{*{20}{c}} {{K_{\text{s}}}}&{ - {K_{\text{s}}}} \\ { - {K_{\text{s}}}}&{{K_{\text{s}}} + {K_{\text{t}}}} \end{array}} \right],\;\; \\& {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {{Z_{\text{s}}}} & {{Z_{\text{t}}}} \end{array}} \right]^{\text{T}},\; \\& \;{\boldsymbol{Y}} = \left[ {\begin{array}{*{20}{c}} 0 & {{K_t}y} \end{array}} \right]^{\text{T}}. \\ \end{split}

    运动方程中阻尼矩阵无法通过实模态变换实现对角化,因此采用复模态变换求解[15]. 首先,令式(8) {\boldsymbol{Y}} =\bf{ 0} 得到自由振动方程,设特解为 {\boldsymbol{\mu }}{{\rm{e}}^{{\boldsymbol{\lambda }}t}} {\boldsymbol{\mu }} λ分别为特征向量和特征值,可求解振动系统的特征方程为

    \left| {{{\boldsymbol{\lambda }}^2}{\boldsymbol{M}} + {\boldsymbol{\lambda C}} + {\boldsymbol{K}}} \right| = {\boldsymbol{0}}. (9)

    二自由度系统存在4个特征值,且为2对共轭复数. 将1/4车模型参数代入特征方程可得4个特征值分别为

    \left\{ {\begin{array}{*{20}{l}} {{\lambda _1} = - 20.4 + 64.4{\rm{i}},\;}\\ {{\lambda _2} = - 2.6 + 7.3{\rm{i}},}\\ {{\lambda _3} = - 20.4 - 64.4{\rm{i}},\;}\\ {{\lambda _4} = - 2.6 - 7.3{\rm{i}}.} \end{array}} \right. (10)

    代入自由振动方程,对应的特征向量构成复模态矩阵{\boldsymbol{\mu }} = \left[ {\begin{array}{*{1}{c}}{{{\boldsymbol{\mu }}_1}}\;\;{{{\boldsymbol{\mu }}_2}}\;\;{{{\boldsymbol{\mu }}_3}}\;\;{{{\boldsymbol{\mu }}_4}}\end{array}} \right],其中:

    \left\{ \begin{array}{*{20}{l}} {{\boldsymbol{\mu }}_1} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{11}}}\;\;{{\mu _{12}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{ - 4.2 + 10.5{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ {{\boldsymbol{\mu }}_2} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{21}}}\;\;{{\mu _{22}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{0.074 + 0.054{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ {{\boldsymbol{\mu }}_3} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{31}}}\;\;{{\mu _{32}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{ - 4.2 - 10.5{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ {{\boldsymbol{\mu }}_4} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{41}}}\;\;{{\mu _{42}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{0.074 - 0.054{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ \end{array} \right. (11)

    由于上述复模态矩阵无法直接变换,因此这里引入状态变量 {\boldsymbol{x}} 对方程解耦[15],补充恒等式{\boldsymbol{M}}\mathop {\boldsymbol{Z}}\limits^ \bullet - {\boldsymbol{M}}\mathop {\boldsymbol{Z}}\limits^ \bullet = {\bf{0}},将运动方程转换为状态方程:

    \mathop {\widetilde{M}}\limits \mathop {\boldsymbol{x}}\limits^ \bullet + \mathop {\widetilde{K}}\limits {\boldsymbol{x}} = \mathop {\widetilde{Y}}\limits \text{,} (12)

    式中:

    {\boldsymbol{x}} = {\left[ {\begin{array}{*{20}{c}} {\mathop {\boldsymbol{Z}}\limits^ \bullet }&{\boldsymbol{Z}} \end{array}} \right]^{\rm{T}}};\;\;\mathop {\boldsymbol{\widetilde{M}}}\limits = \left[ {\begin{array}{*{20}{c}} {\bf{0}}&{\boldsymbol{M}} \\ {\boldsymbol{M}}&{\boldsymbol{C}} \end{array}} \right];
    \mathop {\boldsymbol{\widetilde{K}}}\limits = \left[ {\begin{array}{*{20}{c}} { - {\boldsymbol{M}}}&{\bf{0}} \\ {\bf{0}}&{\boldsymbol{K}} \end{array}} \right];\;\;\mathop {\boldsymbol{\widetilde{Y}}}\limits = {\left[ {\begin{array}{*{20}{c}} {\bf{0}}&{\boldsymbol{Y}} \end{array}} \right]^{\rm{T}}}.

    特征向量构成复模态矩阵:

    {\boldsymbol{V}} = \left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {{\lambda _1}{{\boldsymbol{\mu }}_1}} \\ {{{\boldsymbol{\mu }}_1}} \end{array}}&{\begin{array}{*{20}{c}} {{\lambda _2}{{\boldsymbol{\mu }}_2}} \\ {{{\boldsymbol{\mu }}_2}} \end{array}}&{\begin{array}{*{20}{c}} {{\lambda _3}{{\boldsymbol{\mu }}_3}} \\ {{{\boldsymbol{\mu }}_3}} \end{array}}&{\begin{array}{*{20}{c}} {{\lambda _4}{{\boldsymbol{\mu }}_4}} \\ {{{\boldsymbol{\mu }}_4}} \end{array}} \end{array}} \right] \text{.} (13)

    {\boldsymbol{V}} 具有关于{\boldsymbol{\widetilde{M}}}\mathop {\boldsymbol{\widetilde{K}}}\limits的加权正交性,满足

    {\lambda _i}\mathop {\boldsymbol{\widetilde{M}}}\limits {{\boldsymbol{V}}_i} + \mathop {\boldsymbol{\widetilde{K}}}\limits {{\boldsymbol{V}}_i} = 0,\;\;\;\;i = 1,2,3,4, (14)

    式中: {{\boldsymbol{V}}_i} 为矩阵V的第i个列向量.

    代入\mathop {\boldsymbol{\widetilde{M}}}\limits\mathop {\boldsymbol{\widetilde{K}}}\limits,并作复模态变换 {\boldsymbol{x}} = {\boldsymbol{Vp}} ,其中p为模态变换响应列阵,代入式(12)可对振动方程解耦:

    \begin{split}& {{\rm{diag}}\left( {{m_1},{m_2},{m_3},{m_4}} \right)\mathop {p}\limits^ \bullet + } \;{\rm{diag}} ( - {\lambda _1}{m_1}, - {\lambda _2}{m_2},\\&\quad - {\lambda _3}{m_3}, - {\lambda _4}{m_4} ) {\boldsymbol{p}} = {{\boldsymbol{V}}^{\rm{T}}}{\tilde{\boldsymbol Y}} \text{,}\end{split} (15)

    式中:{m_i} = {{\boldsymbol{V}}_i}^{\rm{T}}\mathop {\boldsymbol{\widetilde{M}}}\limits {{\boldsymbol{V}}_i}

    将式(10)、(11)代入式(15),求得

    \left\{ {\begin{array}{*{20}{l}} {{m_1} = \left( {1\;730 - 1\;774{\rm{i}}} \right){m_{\rm{s}}}},\\ {{m_2} = {m_{\rm{s}}}14{\rm{i}}},\\ {{m_3} = \left( {1\;730 + 1\;774{\rm{i}}} \right){m_{\rm{s}}}},\\ {{m_4} = - {m_{\rm{s}}}14{\rm{i}}}. \end{array}} \right. (16)

    当初始输入为单位脉冲信号,式(15)化简为

    \mathop {{p_i}}\limits^ \bullet - {\lambda _i}{p_i} = \frac{{{\mu _{i2}}}}{{{m_i}}}\delta \left( t \right)\; \text{.} (17)

    式中: {p_i} 为响应列阵中第i个值.

    式(17)等式两边取跨零积分,得

    \int_{0 - \varepsilon }^{0 + \varepsilon } {\left( {\mathop {{p_i}}\limits^ \bullet - {\lambda _i}{p_i}} \right)} {\rm{d}}t = \frac{{{\mu _{i2}}}}{{{m_i}}}\int_{0 - \varepsilon }^{0 + \varepsilon } {\delta \left( t \right){\rm{d}}t} \text{,} (18)

    式中: \varepsilon 为极小量.

    得到系统的单位脉冲激励响应为

    {h_i}(t) = \frac{{{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}} \text{.} (19)

    利用卷积积分方法,可求得在任意激励{K_{\text{t}}}y(t)下的响应为

    {p_i} = \int_0^t {h\left( {t - \tau } \right)} {K_{\text{t}}}y\left( \tau \right){\rm{d}}\tau = \dfrac{{{K_{\text{t}}}{\mu _{i2}}}}{{{m_i}}}\int_0^t {{{\rm{e}}^{{\lambda _i}\left( {t - \tau } \right)}}y\left( \tau \right){\rm{d}}\tau } , (20)

    代入 {\boldsymbol{x}} = {\boldsymbol{Vp}} 得到

    \begin{split}& \left[ {\begin{array}{*{20}{c}} {\mathop {{\rm{ }}{\boldsymbol{Z}}}\limits^ \bullet }\\ {\boldsymbol{Z}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {{\lambda _1}{{\boldsymbol{\mu }}_1}}\\ {{{\boldsymbol{\mu }}_1}} \end{array}\;\;\begin{array}{*{20}{c}} {{\lambda _2}{{\boldsymbol{\mu }}_2}}\\ {{{\boldsymbol{\mu }}_2}} \end{array}\;\;\begin{array}{*{20}{c}} {{\lambda _3}{{\boldsymbol{\mu }}_3}}\\ {{{\boldsymbol{\mu }}_3}} \end{array}\;\;\begin{array}{*{20}{c}} {{\lambda _4}{{\boldsymbol{\mu }}_4}} \\ {{{\boldsymbol{\mu }}_4}} \end{array}} \end{array}} \right]\;\times\\&\quad {\left[ {\begin{array}{*{20}{c}} {{p_1}\;\;{p_2}\;\;{p_3}\;\;{p_4}} \end{array}} \right]^{\rm{T}}}. \end{split} (21)

    可得

    {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {\displaystyle \sum\limits_{i = 1}^4 {{\mu _{i1}}{p_i}} }\\ {\displaystyle \sum\limits_{i = 1}^4 {{\mu _{i2}}{p_i}} } \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{K_{\rm{t}}}{\mu _{i1}}{\mu _{i2}}}}{{{m_i}}}\int\nolimits_0^t {{{\rm{e}}^{{\lambda _i}\left( {t - \tau } \right)}}y\left( \tau \right){\rm{d}}\tau } } }\\ {\displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{K_{\rm{t}}}{\mu _{i2}}^2}}{{{m_i}}}\int\nolimits_0^t {{{\rm{e}}^{{\lambda _i}\left( {t - \tau } \right)}}y\left( \tau \right){\rm{d}}\tau } } } \end{array}} \right], (22)

    {Z_{\text{s}}} 即为 y\left( t \right) 激励下的人体竖向位移.

    由式(22)可知:人体竖向位移的变化振动与系统的参数及路面激励相关,1/4车模型参数固定后,位移则是路面激励的函数. 对位移求二阶导数可以得到加速度函数[16],下面分别推导3类典型局部变形激励下的加速度响应.

    2.3.1   错台型

    车辆在经过错台型变形会经历两个阶段,以桥头跳车为例,第一阶段为驶上但是未离开桥面,第二阶段为驶离桥面. 车辆经过两处错台,引起两次跳车. 但是,实际桥梁具有一定长度,汽车经过第一处错台导致的加速度响应会在短时间内衰减为0,而此时车辆仍未驶下桥面,故每次错台视作为一处变形. 考虑驶上错台的加速度时间历程,即0 \leqslant t \leqslant l/ v,联立式(1)、(22)求解得

    {Z_{\text{s}}} = \left\{ \begin{gathered} 0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g\sum\limits_{i = 1}^4 {\frac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - 1} \right)} ,\;\;\;0 \lt t \leqslant l /v, \hfill \\ \end{gathered} \right. (23)

    对其求二阶导数,得到竖向加速度为

    {a_{\text{s}}} = \left\{ \begin{gathered} 0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}}} ,\;\;\;0 \lt t \leqslant l/ v. \hfill \\ \end{gathered} \right. (24)
    2.3.2   双断面型

    双断面型变形的长度一般较小,例如窨井的直径一般仅为0.5 m左右,导致车辆驶过瞬间会连续经过变形两端,且间隔时间极短,两次振动响应会发生叠加,所以车轮驶上驶下变形的两次冲击都需要考虑在内,加速度响应为分段函数,联立式(2)、(22)求解,结果如式(25)所示.

    {Z_{\text{s}}} = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} 0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - 1} \right),\;\;\;0 \lt t \lt d/v,} \hfill \\ \end{gathered} \\ {{K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - {{\rm{e}}^{{\lambda _i}\left( {t - d/v} \right)}}} \right),\;\;\;t \geqslant d/v.} } \end{array}} \right. (25)

    对式(25)求二阶导,得到竖向加速度为

    {a_{\text{s}}} = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} \;0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}},\;\;\;0 \lt t \lt d/v,} \hfill \\ \end{gathered} \\ {{K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - {{\rm{e}}^{{\lambda _i}\left( {t - d/v} \right)}}} \right),\;\;\;t \geqslant d/v.} } \end{array}} \right. (26)
    2.3.3   曲面型

    此类变形断面为平滑曲线,在变形长度较大时,振动响应会减弱. 联立式(3)、(22)求解,结果如式(27)所示. 同理可得竖向加速度结果如式(28)所示.

    \qquad \qquad \begin{split} {Z_{\rm{s}}} = \left\{ \begin{array}{l} \;0\;,\;\;\;t = 0,\;\;\\ \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{4\Delta gv{K_{\rm{t}}}{\mu _{i2}}}}{{{\lambda _i}^3{m_i}{l^2}}}\left( \begin{array}{l} v{\lambda _i}^2{t^2} - 2v{{\rm{e}}^{{\lambda _i}t}} + {\lambda _i}l{{\rm{e}}^{{\lambda _i}t}} - {\lambda _i}^2lt + 2v{\lambda _i}t - {\lambda _i}l + 2v \end{array} \right),} \;\;\;0 < t < l/v,\\ \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{4\Delta gv{K_{\rm{t}}}{\mu _{i2}}}}{{{\lambda _i}^3{m_i}{l^2}}}{{\rm{e}}^{{\lambda _i}\left( {t - l/v} \right)}}\left( \begin{array}{l} {\lambda _i}l{{\rm{e}}^{{\lambda _i}l/v}} - 2v{{\rm{e}}^{{\lambda _i}l/v}} + {\lambda _i}l + 2v \end{array} \right),} \;\;\;t \geqslant l/v. \end{array} \right. \end{split} (27)
    {a_{\text{s}}} = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} \;0,\;\;\;t = 0,\;\; \hfill \\ \sum\limits_{i = 1}^4 {\frac{{4\Delta gv{K_{\text{t}}}{\mu _{i2}}}}{{{\lambda _i}{m_i}{l^2}}}\left( {\left( {{\lambda _i}l - 2v} \right){{\rm{e}}^{{\lambda _i}t}} + 2v} \right),\;\;\;0 \lt t \lt l/v,} \hfill \\ \end{gathered} \\ \begin{gathered} \sum\limits_{i = 1}^4 {\frac{{4\Delta gv{K_{\text{t}}}{\mu _{i2}}}}{{{\lambda _i}{m_i}{l^2}}}{{\rm{e}}^{{\lambda _i}\left( {t - l/v} \right)}}\left( {{\lambda _i}l{{\rm{e}}^{{\lambda _i}l/v}} + {\lambda _i}l - 2v{{\rm{e}}^{{\lambda _i}l/v}} + 2v} \right),} \hfill \\ \;\;\;\;\;\;\;t \geqslant l/v.\;\; \hfill \\ \end{gathered} \end{array}} \right. (28)

    局部变形激励导致的人车振动为瞬时冲击,不具周期性特征,ISO 2631-1[11]推荐采用运行均方根法来评价这种路面偶然激励导致的振动感受,其形式为

    {a_w}\left( {{t_0}} \right) = \sqrt {\frac{1}{{{\tau _0}}}\int_{{t_0} - {\tau _0}}^{{t_0}} {a_{\rm{w}}^2\left( t \right){\rm{d}}t} } \text{,} (29)

    式中: {t_0} 为瞬时激励发生时刻; {\tau _0} 为持续时间; {a_{\rm{w}}}\left( t \right) 为瞬时频率加权加速度.

    欲求 {a_{\rm{w}}}\left( t \right) ,需要对竖向加速度的时域值进行频率加权滤波. 由第2节可知,竖向加速度时域值与行车速度、变形长度、变形高差相关,因本文主要研究不同类型局部变形引起的人体振动,故先统一车速和变形长度两个变量. 根据现行城市道路规范[17],道路设计速度为20~60 km/h,本文车速取40 km/h;实际路面中双断面型变形多表现为窨井的凸起或沉陷损坏,窨井的直径一般为0.6~0.7 m[2],本文取0.6 m.

    人体对不同频率的振动敏感程度不同,一般认为0.5~80.0 Hz是振动舒适性评价的考虑范围,竖向振动频率加权函数[18]

    {W_{\text{k}}}\left( f \right) = \left\{ \begin{gathered} 0.5,\;\;\;0.5 \lt f \leqslant 2.0, \hfill \\ {f/4},\;\;\;2.0 \lt f \leqslant 4.0, \hfill \\ 1.0,\;\;\;4.0 \lt f \leqslant 12.5, \hfill \\ {12.5/f},\;\;\;12.5 \lt f \leqslant 80.0, \hfill \\ \end{gathered} \right. (30)

    式中:f 为竖向振动频率.

    {W_{\text{k}}} 在低频段和高频段均呈现抑制性,而在中间频段则具有较好的通过性,所以可以设计带通滤波器进行频率加权滤波. 因运行均方根计算对加速度相位要求不高,故采用IIR (infinite impulse response)数字滤波器设计,并借助MATLAB信号处理工具箱完成. 巴特沃斯滤波器原型具有通带内最大平坦和随频率升高而单调递减特性,与标准中的 {W_{\text{k}}} 具有较高契合度,因此本文采用巴特沃斯带通滤波器.

    数字滤波器通过间接法设计,步骤如下:

    步骤1 确定数字滤波器技术指标;

    步骤2 将数字滤波器技术指标转换为模拟低通滤波器技术指标;

    步骤3 设计模拟低通滤波器(滤波器原型);

    步骤4 通过频率转换,将模拟低通滤波器转变为模拟带通等滤波器;

    步骤5 模拟滤波器数字化.

    由MATLAB中buttord函数先确定滤波器阶数,再由buttap函数设计出滤波器原型,并借助lp2bp函数实现模拟低通至模拟带通的转换,最后由bilinear函数完成数字化.

    根据 {W_{\text{k}}} 渐近线形式,频率加权函数在4.0~12.5 Hz,权重保持为1,考虑巴特沃斯滤波器的过渡特性,试算后取通带范围为2.5~18.0 Hz,通过上述过程可最终设计出频率加权函数滤波器,本文同时设计了一个对比滤波器,两者的区别在于阻带最小衰减量不同. 设计滤波器的幅频特性与ISO 2631-1[11]中的频率加权函数对比如图4所示.

    图  4  竖向加速度的频率加权曲线
    Figure  4.  Frequency weighting curves of vertical acceleration

    图4可见,设计滤波器在高频段具有一定误差,但在25.0 Hz以下人体最敏感的低频段吻合度高,且在80.0 Hz之后迅速衰减为0,而对比滤波器则与之相反. 考虑到路面局部变形的激振频率主要为低频振动(如图5),高频段的权重实际意义不大,即保证通带频率范围的特性更为关键,因此,本文设计的带通滤波器适用性更高.

    图  5  竖向加速度频谱
    Figure  5.  Spectrum of vertical acceleration

    因局部变形高差与时间变化无关,可视为固定系数,因此,将代表性车速和变形长度代入式(24)、(26)和式(28)可得竖向加速度的时域函数,通过带通滤波器后可得到瞬时频率加权加速度,频率记权前后的加速度变化情况如图6所示.

    根据运行均方根的定义,通过数值积分方法计算得到3类局部变形作用下的均方根加速度,如式(31)~(33)所示. 其中,加速度持续时间取为1 s[11].

    图  6  滤波前后的竖向加速度对比
    Figure  6.  Vertical acceleration comparison before and after filtering

    1) 错台型

    {a_{{\rm{w}}}}= {K_{\text{t}}}\Delta g{\left[ {{{\int_0^1 {\left( {\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}}} } \right)} }^2}{\text{d}}t} \right]^{0.5}} = 42.04\Delta g \text{.} (31)

    2) 双断面型

    \begin{split}& {a_{{\rm{w}}}} = {K_{\text{t}}}\Delta g \;\times \hfill \\&\;\; {\left[ {\int_0^{\tfrac{l}{v}} {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}}} } \right)}^2}} {\rm{d}}t + \int_{\tfrac{l}{v}}^1 {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - {{\rm{e}}^{{\lambda _i}\left( {t - \tfrac{l}{v}} \right)}}} \right)} } \right)}^2} {\rm{d}}t} } \right]^{0.5}} \hfill \\&\;\; =72.07\Delta g\text{.} \\ \end{split} (32)

    3) 曲面型

    {a_{{\rm{w}}}} = \frac{{4\Delta gv{K_{\text{t}}}}}{l}{\left( {A + B} \right)^{0.5}}= 61.28\Delta g \text{,} (33)

    式中:

    \begin{gathered} A = \int_0^{\tfrac{l}{v}} {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {\left( {{\lambda _i}l - 2v} \right){{\rm{e}}^{{\lambda _i}t}} + 2v} \right)} } \right)}^2}} {\rm{d}}t ;\hfill \\ B = \int_{\tfrac{l}{v}}^1 {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}{{\rm{e}}^{{\lambda _i}\left( {t - l/v} \right)}}\left( {{\lambda _i}l{{\rm{e}}^{{\lambda _i}l/v}} + {\lambda _i}l - 2v{{\rm{e}}^{{\lambda _i}l/v}} + 2v} \right)\;} } \right)}^2} {\rm{d}}t}. \hfill \\ \end{gathered}

    变形高差与均方根加速度之间的理论关系如图7所示. 结合ISO 2631-1[11]给出的舒适度评价标准,本文计算出的局部变形评价标准如表1所示. 当错台型、双断面型及曲面型局部变形分别达到 4.76、2.78、3.26 cm时,行驶质量均达到极不舒适.

    综合上述分析,3类局部变形激励作用时的均方根加速度均和变形高差呈正比关系;在典型行驶车速40 km/h和变形长度0.6 m情况下,同样的变形高差作用,双断面型对人体舒适性影响最强烈,曲面型次之,错台型影响相对最小. 值得注意的是,实际路面中,曲面型变形多表现为路面拥包等损坏,长度不一. 当式(33)中局部变形长度l较大时,均方根加速度aw 的比例系数会减小,如图8所示. 本文所取的0.6 m变形长度几乎为曲面型变形的最不利长度,当长度约为0.3 ~ 1.5 m,同等高差时曲面型影响大于错台型,路面养护管理中应当重点考虑该长度范围.

    图  7  均方根加速度与局部变形高差的关系
    Figure  7.  Relationship between root-mean-square acceleration and height difference of local deformation
    表  1  局部变形高差评价标准
    Table  1.  Evaluation criteria for height difference of local deformation
    均方根加速度/
    (m • s−2)
    振动感受g 范围/cm
    错台型双断面型曲面型
    < 0.315 没有不
    舒适
    < 0.75 < 0.45 < 0.51
    0.315~0.630 有些不
    舒适
    0.75~
    1.50
    0.45~
    0.87
    0.51~
    1.03
    0.500~1.000 比较不
    舒适
    1.19~
    2.38
    0.69~
    1.39
    0.82~
    1.63
    0.800~1.600 不舒适 1.90~
    3.81
    1.11~
    2.22
    1.31~
    2.61
    1.250~2.500 很不舒适 2.97~ 5.95 1.73~ 3.47 2.04~ 4.08
    > 2.000 极不舒适 > 4.76 > 2.78 > 3.26
    下载: 导出CSV 
    | 显示表格
    图  8  变形长度与比例系数的关系
    Figure  8.  Relationship between length of deformation and evaluation coefficient

    1) 对路面中常见的局部变形形式进行分析总结,建立了错台型、双断面型和曲面型3类典型局部变形的纵断面模型.

    2) 求解两自由度汽车强迫振动模型,推导出任意路面变形激励作用时人体竖向位移的时域变化函数;根据3类局部变形的纵断面特征从理论上得出人体竖向加速度的时域解,并对函数式中变形长度和行驶车速两个参数作出合理假定.

    3) 依据ISO 2631-1评价标准中竖向加速度的频率加权系数分布,设计了频率加权函数数字带通滤波器;对3类局部变形激励下的竖向加速度信号进行滤波,得到瞬时频率加权加速度并求解出均方根加速度;最后根据标准中舒适度划分方法确定了基于舒适性的3类局部变形评价标准. 结果表明:3类局部变形激励作用下,均方根加速度均与局部变形高差呈正比关系,同等变形程度下,双断面型变形的行驶不适度最高.

    4) 路面养护管理中,可结合高精度检测设备及时获取路面严重错台型及双断面型局部变形的位置,做好养护和修复工作,提高路面服务性能和寿命.

    5) 文中研究结果适用于城市道路沥青路面,城市道路水泥路面可以参考;对于公路路面中的类似问题,需要做进一步的研究.

  • 图 1  路面局部变形因素统计

    Figure 1.  Statistical factors of local deformation in pavement

    图 2  典型路面局部变形纵断面示意

    Figure 2.  Profile diagrams of typical local deformation

    图 3  1/4车模型

    Figure 3.  Quarter car model

    图 4  竖向加速度的频率加权曲线

    Figure 4.  Frequency weighting curves of vertical acceleration

    图 5  竖向加速度频谱

    Figure 5.  Spectrum of vertical acceleration

    图 6  滤波前后的竖向加速度对比

    Figure 6.  Vertical acceleration comparison before and after filtering

    图 7  均方根加速度与局部变形高差的关系

    Figure 7.  Relationship between root-mean-square acceleration and height difference of local deformation

    图 8  变形长度与比例系数的关系

    Figure 8.  Relationship between length of deformation and evaluation coefficient

    表  1  局部变形高差评价标准

    Table  1.   Evaluation criteria for height difference of local deformation

    均方根加速度/
    (m • s−2)
    振动感受g 范围/cm
    错台型双断面型曲面型
    < 0.315 没有不
    舒适
    < 0.75 < 0.45 < 0.51
    0.315~0.630 有些不
    舒适
    0.75~
    1.50
    0.45~
    0.87
    0.51~
    1.03
    0.500~1.000 比较不
    舒适
    1.19~
    2.38
    0.69~
    1.39
    0.82~
    1.63
    0.800~1.600 不舒适 1.90~
    3.81
    1.11~
    2.22
    1.31~
    2.61
    1.250~2.500 很不舒适 2.97~ 5.95 1.73~ 3.47 2.04~ 4.08
    > 2.000 极不舒适 > 4.76 > 2.78 > 3.26
    下载: 导出CSV
  • [1] 孙立军. 道路与机场设施管理学[M]. 北京: 人民交通出版社, 2009.
    [2] 刘婷. 城市道路平整度评价方法研究[D]. 上海: 同济大学, 2013.
    [3] CANTISANI G, LOPRENCIPE G. Road roughness and whole body vibration:evaluation tools and comfort limits[J]. Journal of Transportation Engineering, 2010, 136(9): 818-826. doi: 10.1061/(ASCE)TE.1943-5436.0000143
    [4] 王锋,张金喜. 基于驾乘舒适性的沥青混凝土路面评价[J]. 西南交通大学学报,2014,49(4): 687-692. doi: 10.3969/j.issn.0258-2724.2014.04.020

    WANG Feng, ZAHNG Jinxi. Evaluation of asphalt concrete pavement based on passenger ride comfort level[J]. Journal of Southwest Jiaotong University, 2014, 49(4): 687-692. doi: 10.3969/j.issn.0258-2724.2014.04.020
    [5] AHN S J. Discomfort of vertical whole-body shock-type vibration in the frequency range of 0.5 to 16 Hz[J]. International Journal of Automotive Technology, 2010, 11(6): 909-916. doi: 10.1007/s12239-010-0108-z
    [6] HUSTON D R, ZHAO X, JOHNSON C C. Whole-body shock and viabration:frequency and amplitude dependence of comfort[J]. Journal of Sound and Vibration, 2000, 230(4): 964-970.
    [7] 陶向华. 路桥过渡段差异沉降控制标准与人车相互作用[D]. 南京: 东南大学, 2005.
    [8] JAME H L, SCOTT M O, TIMOTHY D S. Differential movement at emankment bridge structure interface in Illinois[J]. Transportation Research Record, 1998, 1633: 53-60. doi: 10.3141/1633-07
    [9] 何智龙. 上海道路窨井病害调研及机理分析[J]. 城市道桥与防洪,2010(8): 51-53,57. doi: 10.3969/j.issn.1009-7716.2010.08.016

    HE Zhilong. Investigation and mechanism analysis of damage of road basket well in Shanghai[J]. Urban Roads Bridges & Flood Control, 2010(8): 51-53,57. doi: 10.3969/j.issn.1009-7716.2010.08.016
    [10] 王雪如. 基于乘车舒适性的高速公路路面平整度评价方法研究[D]. 北京: 北京工业大学, 2014.
    [11] International Organization for Standardization. Mechanical vibration and shock—evaluation of human exposure to whole-body vibration—part 1, general requirements: ISO 2631-1[S]. Geneva: [s.n.], 1997.
    [12] 黄舫. 考虑路面局部不平整的城市道路行驶质量评价方法研究[D]. 上海: 同济大学, 2015.
    [13] SAYERS M W, GILLESPIE T D, PATERSON W. Guidelines for conducting and calibrating road roughness measurements[M]. Washington D C: International Bank for Reconstruction and Development, 1986.
    [14] 靳晓雄, 张立军, 江浩. 汽车振动分析[M]. 上海: 同济大学出版社, 2002.
    [15] 李惠彬. 振动理论与工程应用[M]. 北京: 北京理工大学出版社, 2006.
    [16] 彭华. 交通资产管理系统优化框架的研究[D]. 上海: 同济大学, 2010.
    [17] 中华人民共和国住房和城乡建设部. 城市道路工程技术规范: GB 51286—2018[S]. 北京: 中国建筑工业出版社, 2018.
    [18] 余志生. 汽车理论[M]. 北京: 机械工业出版社, 2004.
  • 加载中
图(8) / 表(1)
计量
  • 文章访问数:  291
  • HTML全文浏览量:  186
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-15
  • 修回日期:  2020-07-28
  • 刊出日期:  2020-08-24

目录

/

返回文章
返回