疲劳预警分级提示下的驾驶行为响应特征

胥川; 郭启明; 王雪松

doi:10.3969/j.issn.0258-2724.20180254

疲劳预警分级提示下的驾驶行为响应特征

doi: 10.3969/j.issn.0258-2724.20180254

胥川^1,2,,
郭启明^3,4,
王雪松^3,4,5, ,

基金项目: 国家自然科学基金资助项目（61703352）；中央高校基本科研业务费专项资金资助项目（2682017CX021）；上海市科学技术委员会资助课题（15DZ1207502）

详细信息

作者简介:
胥川（1987—），男，讲师，博士，研究方向为交通安全、数据科学等，E-mail： xuchuan@swjtu.edu.cn

通讯作者:
王雪松（1977—），男，教授，研究方向为交通安全、交通统计分析、交通规划、驾驶模拟器应用等，E-mail： wangxs@tongji.edu.cn

中图分类号: U491.25
计量
- 文章访问数: 589
- HTML全文浏览量: 288
- PDF下载量: 33
- 被引次数: 17
出版历程
- 收稿日期: 2018-04-04
- 修回日期: 2018-07-02
- 网络出版日期: 2018-07-08
- 刊出日期: 2019-02-01

Driver Behavior Response to Drowsiness Alarming at Different Levels

XU Chuan^{1,2
,},
GUO Qiming^3,4,
WANG Xuesong^{3,4,5
, ,}

摘要

摘要: 确定疲劳预警提示时机是车载疲劳预警提示技术中的关键和难点，通过驾驶员对预警提示的驾驶行为响应特征来判断提示时机的合理性是解决该问题的新思路. 通过开展在疲劳分级预警提示环境下的驾驶模拟试验，采集驾驶员眼动（眼闭合时间比例）和车道偏移（车道偏移标准差、车道偏移均值、车道越线面积）指标，并通过成对Wilcoxon signed-rank检验，分析不同等级预警提示前后15 s时间内关键驾驶行为指标的差异. 研究结果表明：在设定的疲劳分级规则下，一般疲劳预警提示后，车道偏移标准差均值显著下降0.129 1，车道越线面积均值显著下降8.574 4；在严重疲劳预警提示后，眼闭合时间比例显著下降0.044 9，但车道偏移标准差和车道越线面积均未发生显著改变，驾驶员应该尽快停车休息.
- 驾驶模拟器 /
- 驾驶行为 /
- 提示时机 /
- 提示响应 /
- 疲劳驾驶
Abstract: To identify the timing of drowsiness driving warning is the key issue and a bottleneck of onboard drowsiness driving warning technology. Finding a rationale for warning timing using driver’s driving behavior response feature is an innovation. Therefore, after conducting a driving simulator experiment under the influence of drowsiness alarming, eye movement index, the percentage of eyelid closure (PERCLOS), and vehicle lateral position indexes (standard deviation of lateral position, average of lateral position, area of line exceeding) were recorded. Then, the differences in driving behavior between the pre-warning and post-warning in 15 seconds for each warning were compared using the paired Wilcoxon signed-rank test. The results demonstrated that under the drowsiness classification criterion, after the normal level of drowsiness warning, the mean of the standard deviation of lateral position and the mean of the area of line exceeding significantly dropped down by 0.129 1 and 8.574 4 respectively; after the serious level of drowsiness warning, although the mean of PERCLOS decreased by 0.044 9, the standard deviation of the lateral position and area of line exceeding did not significantly change and that the driver should stop and rest immediately.
- driving simulator /
- driving behavior /
- warning time /
- warning response /
- drowsiness driving

HTML全文

高承载能力和低固有频率之间的矛盾是低频隔振技术发展的瓶颈^[1-3]. 具有高静-低动刚度特性的非线性隔振系统可以很好地解决上述问题. Carrella等^[4-5]将一个垂直弹簧与2个倾斜弹簧并联，当承载质量偏离平衡位置时，倾斜弹簧产生的负刚度与垂直弹簧的正刚度抵消，构建了具有高静-低动刚度特性的隔振系统；周加喜等^[6]采用滚轮-凸轮结构设计了一种非线性隔振系统，当凸轮离开平衡位置，带有滚轮的水平弹簧与凸轮组成的机构在竖直方向产生负刚度，从而抵消竖直弹簧的正刚度. 上述方式往往存在布局、结构相对较大，一旦结构尺寸固定，起始隔振频率也就固定，无法对不同振源做出适应性调整的问题^[7-10]. 除此之外，目前还有采用欧拉梁或气动弹簧产生负刚度来构建非线性隔振系统^[7-8].

近年来，随着永磁材料性能的提升，采用永磁材料间的非线性力学特性构建具有高静-低动刚度特性的隔振系统，受到越来越多的重视. 李爽等^[11]利用双环永磁体结构提出一种具有高静-低动刚度特性的隔振器，双环永磁体提供负刚度，与机械弹簧产生的正刚度抵消. Zhou等^[12-13]使用2块电磁铁与一块永磁铁设计了一种磁力弹簧机构，通过调节电磁铁的电流大小与方向实现了磁力弹簧负刚度可调节.

借鉴上述方式的优点，本文提出一种利用永磁体结合空心电磁线圈，构建具有高静-低动刚度特性的非线性隔振系统，采用永磁体直接构建正负刚度，且提供轴向的承载力，通过改变电磁线圈中的电流大小、方向来改变结构中的电磁磁通，从而调节系统轴向刚度，以达到改变起始隔振频率的目的，且该系统结构更紧凑. 采用增量谐波平衡法求解该系统的位移传递率，考虑了传统方法分析二次非线性项较困难的问题^[14]，并分析了线圈电流变化对系统位移传递率的影响.

1. 隔振系统结构与工作原理

1.1 结　构

提出了一种具有高静-低动刚度特性的非线性隔振系统，由3个环形永磁体和2个空心线圈构成，如图1（a）. 3个永磁体同轴排列，两端电磁线圈分别与对应永磁体同中心同轴布置，且两端永磁体和线圈通过连接件刚性连接在框架的特定位置，框架刚性连接在基座上. 中间永磁体固定在轴上，随轴进行轴向的往复移动，中间永磁体通过轴与质量块连接. 所有永磁体都是轴向磁化的，相邻的永磁体磁化方向相反. 2个线圈中的电流极性始终相同，永磁体的磁化方向和线圈中的电流极性如图1（b）所示，图中绿色箭头为永磁体磁化方向，红色箭头为电流极性.

图 1 高静-低动刚度隔振系统

Figure 1. Vibration isolation system with high static stiffness and low dynamic stiffness

下载: 全尺寸图片幻灯片

永磁体嵌套在电磁线圈中，减小了整体结构的尺寸. 如图1（b）所示，底部磁环厚度大于上部和中部磁环，会在平衡位置产生一个向上的磁力差. 如图2所示，当隔振器承载与差值相等的隔振质量m时，中间磁环刚好回归平衡位置. 图中：F为轴向力，Δz为在平衡位置处的微小位移，z为轴向位移.

图 2 隔振系统轴向力示意

Figure 2. Axial force of vibration isolation system

下载: 全尺寸图片幻灯片

1.2 工作原理

永磁环间磁力差与隔振质量的重力保持平衡，使中间磁环保持在平衡位置. 未加载电流时，依靠永磁环间非线性斥力实现被动隔振，此时系统刚度为永磁体间磁力刚度. 当线圈加载额定的正向激励电流时，线圈与磁环耦合将产生排斥的电磁力，与两端磁环产生的斥力叠加，施加在中间磁环上的合力和系统的等效刚度都将增加，如图3（a），图中绿色箭头为斥力，红色箭头为吸引力. 当线圈加载额定的反向激励电流时，线圈与磁环耦合将产生吸引的电磁力，与两端磁环产生的斥力抵消，施加在中间磁环上的合力和系统的等效刚度都将减小，如图3（b）.

图 3 系统工作原理

Figure 3. System working principle

下载: 全尺寸图片幻灯片

2. 结构参数设计

该隔振系统涉及的结构参数主要包括磁环的内（外）半径R_m（r_m）和轴向厚度L_m，线圈的内（外）半径R_c（r_c ，匝数决定）和轴向厚度L_c. 考虑到磁环需要连接轴，中间需要设计圆孔，永磁铁嵌套在线圈中，所以磁环内径r_m，线圈内径r_c值固定，不再多做讨论. 所以最终优化设计参数为R_m、L_m、R_c、L_c. 以L_m为基准定义优化变量，磁环径厚比A=R_m/L_m，变化范围为[1.0,4.0]时，系统刚度变化曲线如图4.

图 4 A对系统刚度的影响

Figure 4. Effect of A on system stiffness

下载: 全尺寸图片幻灯片

系统刚度随着动磁环位移增加呈非线性规律增加，位移较小时，系统刚度较小；随着A的增大，系统的等效刚度整体增大，刚度变化平稳区域减小，非线性变化规律逐渐增强.

由工作原理可知，当通入负电流时可降低系统的等效刚度，所以，在线圈与磁环的厚度比变化和线圈径厚比变化的仿真时给定电流−8 A，观察系统刚度降低的变化.

线圈与磁环的厚度比B=L_c/L_m，变化范围为[1.0,4.0]. 线圈匝数不变，电流一定时，系统刚度的变化如图5所示. 起初系统起始刚度随着B增大而减小，但当B增加到一定值时，系统的起始刚度会增加. 然而，系统的刚度变化平稳区域却始终减小，非线性变化趋势增强，且随着位移增加，B=3.5和B=4.0的最大刚度会超过B=1.0.

图 5 B 比对系统刚度的影响

Figure 5. Effect of B on system stiffness

下载: 全尺寸图片幻灯片

影响系统刚度变化的另一指标就是线圈的匝数，将匝数与线圈的半径关联，定义优化变量X=R_c/L_c，取值范围为[1.5,4.0]，以线圈厚度L_c为基准，则X增加，线圈的匝数增加. 如图6所示，随着X增大，系统刚度降低，刚度变化平稳区域减小，随着X继续增大，X对系统刚度的影响逐渐减小，当X增加到一定值时，其对系统刚度影响作用趋于饱和.

图 6 X对系统刚度的影响

Figure 6. Effect of X on system stiffness

下载: 全尺寸图片幻灯片

通过仿真计算并结合工程实际情况，选取表1中的参数作为非线性隔振系统的结构尺寸.

表 1 结构参数

Table 1. Structural parameters mm

参数	数值
磁环内径	5
磁环外径	12
磁环厚度	8（上、中），12（下）
线圈内径	15
线圈外径	30
线圈厚度	16
工作气隙 Z	15
线径 s	1

下载: 导出CSV

| 显示表格

3. 力学模型

3.1 磁环间轴向磁力

根据分子电流法，假设永磁环由2组电流大小相等、方向相反的等效表面电流构成^[15-16]. 图7（a）为一对永磁环的分子电流模型，图中：“1”和“2”分别为永磁环1的内外表面，“3”和“4”分别为永磁环2的内外表面；L_g为两磁环间的气隙；L_m1、L_m2分别为磁环1、2的厚度. R₁、R₂、R₃、R₄分别为面1、2、3、4的半径，h₁、h₂分别为面“2”“3”上分子电流环的轴向位置.

图 7 1组磁环的分子电流模型

Figure 7. Molecular current model of magnetic rings

下载: 全尺寸图片幻灯片

以磁环1的“2”面和磁环2的“3”面为例，计算其作用力F₂₃. 图7（b）为“2”面上任意环电流对“3”面上任意环电流之间的作用力示意. 图中：P、Q为在“2”面和“3”面电流环上任意取的电流微元，θ_P、θ_Q分别为P、Q点在O-xy坐标平面上相对于x的夹角，i₁、i₂分别为“2”面、“3”面的等效电流强度，r为点P到点Q的标量距离.

根据毕奥-萨瓦尔-拉普拉斯定律可以得到载流环上点P在点Q处产生的磁感应强度为

$\text{d}{\boldsymbol{B}}=\frac{{\mu }_{0}{i}_{1}}{4{\text{π}} {{\boldsymbol{r}}}^{3}}{\text{•}}\left[\text{d}{{\boldsymbol{l}}}_{1}\times {\boldsymbol{r}}\right]{\text{•}}\text{d}{h}_{1} \text{，}$

(1)

式中：dl₁、dh₁为“2”面上假想电流微元的长度和厚度，μ₀为真空磁导率，r为Q点相对于P点的矢量位置.

由几何关系可以推出

$\begin{split} & {\boldsymbol{r}} = ( - {R_2}\cos \;{\theta _{\text{P}}} + {R_3}\cos \;{\theta _{\text{Q}}}){\boldsymbol{i}} + \\ &\quad ( - {R_2}\sin \;{\theta _{\text{P}}} + {R_3}\sin \;{\theta _{\text{Q}}}){\boldsymbol{j}} + \\ &\quad ({L_{{\mathrm{m}}1}} - {h_1} + {L_{\mathrm{g}}} + {h_2}){\boldsymbol{k}} , \end{split}$

(2)

${\text{d}}{{\boldsymbol{l}}_1} = - {R_2}\sin \;{\theta _{\text{P}}}{\text{d}}{\theta _{\text{P}}}{\boldsymbol{i}} + {R_2}\cos \;{\theta _{\text{P}}}{\text{d}}{\theta _{\text{P}}}{\boldsymbol{j}} + 0{\boldsymbol{k}} \text{，}$

(3)

$\begin{split} & r = [{( - {R_2}\cos \;{\theta _{\text{P}}} + {R_3}\cos \;{\theta _{\text{Q}}})^2} + \\ &\quad {( - {R_2}\sin \;{\theta _{\text{P}}} + {R_3}\sin \;{\theta _{\text{Q}}})^2} + \\ &\quad {({L_{{\mathrm{m}}1}} - {h_1} + {L_{\mathrm{g}}} + {h_2})^2}]^{1/2} , \end{split}$

(4)

式中：i、j、k分别为x轴、y轴、z轴的正向单位矢量.

对dB积分，可以得到磁环1环面“2”上的分子电流在点Q处的磁感应强度为

$\begin{split} & {\boldsymbol{B}} = \frac{{{\mu _0}{i_1}{R_2}}}{{4{\text{π}}}}\int_0^{2{\text{π}}} {\int_0^{{L_{{\mathrm{m}}1}}} {\frac{{({L_{{\mathrm{m1}}}} + {L_{\mathrm{g}}} - {h_1} + {h_2})\cos \;{\theta _{\text{P}}}}}{{{r^3}}}} } {\boldsymbol{i}} + \\ &\quad \frac{{({L_{{\mathrm{m1}}}} + {L_{\mathrm{g}}} - {h_1} + {h_2})\sin \;{\theta _{\text{P}}}}}{{{r^3}}}{\boldsymbol{j}} + \\ &\quad \frac{{{R_2} - {R_3}(\sin \;{\theta _{\text{P}}}\sin \;{\theta _{\text{Q}}} + \cos \;{\theta _{\text{P}}}\cos \;{\theta _{\text{Q}}})}}{{{r^3}}}{\boldsymbol{k}}{\text{d}}{\theta _{\text{P}}}{\text{d}}{h_1} . \end{split}$

(5)

由电磁场理论可知，点Q处电流微元所受到的作用力为

${\text{d}}{{\boldsymbol{F}}_{23}} = \left[ {{\text{d}}{{\boldsymbol{l}}_{\text{2}}} \times {\boldsymbol{B}}} \right]{i_2}{\text{d}}{h_2} \text{，}$

(6)

式中： ${\text{d}}{{\boldsymbol{l}}_{{2}}} = - {R_3}\sin \;{\theta _{\text{Q}}}{\text{d}}{\theta _{\text{Q}}}{\boldsymbol{i}} + {R_3}\cos \;{\theta _{\text{Q}}}{\text{d}}{\theta _{\text{Q}}}{\boldsymbol{j}} + 0{\boldsymbol{k}}$ ，dl₂、dh₂分别为“3”面上假想电流微元的长度和厚度.

对dF₂₃进行积分，得到“2”面和“3”面之间的磁力为

$\begin{split} & {{\boldsymbol{F}}_{23}} = \frac{{{\mu _0}{i_1}{i_2}{R_2}{R_3}}}{{4{\text{π}} }}\int_0^{2{\text{π}} } {\int_0^{2{\text{π}} } {\int_0^{{L_{{\mathrm{m1}}}}} {\int_0^{{L_{{\mathrm{m2}}}}} {} } } } \\ &\quad \frac{{\cos\; \beta ({R_2} - {R_3}\cos ({\theta _{\text{P}}} - {\theta _{\text{Q}}}))}}{{{r^3}}}{\boldsymbol{i}} + \\ &\quad \frac{{\sin\; \beta ({R_2} - {R_3}\cos ({\theta _{\text{P}}} - {\theta _{\text{Q}}}))}}{{{r^3}}}{\boldsymbol{j}} + \\ &\quad \frac{{({L_{{\mathrm{m1}}}} + {L_g} - {h_1} + {h_2})\cos ({\theta _{\text{P}}} - {\theta _{\text{Q}}})}}{{{r^3}}}{\boldsymbol{k}}{\text{d}}{\theta _{\text{P}}}{\text{d}}{\theta _{\text{Q}}}{\text{d}}{h_1}{\text{d}}{h_2}. \end{split}$

(7)

该隔振系统只有轴向磁力起作用，所以式（7）中只需要保留k项，“2”面和“3”面之间的轴向磁力为

$\begin{split} & {F_{{\mathrm{m}}23}} = \frac{{{\mu _0}{i_1}{i_2}{R_2}{R_3}}}{{4{\text{π}}}}\int_0^{2{\text{π}}} {\int_0^{2{\text{π}} } {\int_0^{{L_{{\mathrm{m}}1}}} {\int_0^{{L_{{\mathrm{m}}2}}} {} } } } \\ &\quad \frac{{({L_{{\mathrm{m}}1}} + {L_g} - {h_1} + {h_2})\cos ({\theta _{\text{P}}} - {\theta _{\text{Q}}})}}{{{r^3}}}{\text{d}}{\theta _{\text{P}}}{\text{d}}{\theta _{\text{Q}}}{\text{d}}{h_1}{\text{d}}{h_2} . \end{split}$

(8)

同理可分别推导出“1”“3”面、“1”“4”面和“2”“4”面之间的轴向磁力F_m13、F_m14、F_m24，可得两磁环间轴向磁力为

${F_{{\mathrm{m}}1}} = {F_{{\mathrm{m}}13}} + {F_{{\mathrm{m}}14}} + {F_{{\mathrm{m}}23}} + {F_{{\mathrm{m}}24}} .$

(9)

同理可求另一侧永磁环间的轴向磁力F_m2，系统中永磁环产生的总轴向磁力为

${F_{\mathrm{m}}} = {F_{{\mathrm{m}}1}} + {F_{{\mathrm{m}}2}}.$

(10)

3.2 线圈与磁环间的轴向磁力

将线圈等效为多层同向表面电流叠加，永磁环依然等效为2组电流大小相等但方向相反的等效表面电流.

如图8所示，为线圈-磁环分子电流模型示意，图中：“1”和“2”分别为磁环的内、外表面，h为磁环面“1”上分子电流环轴向位置，b、v分别为线圈上分子电流环的轴向与径向位置，L_cg为磁环与线圈间的气隙.

图 8 磁环与多层载流线圈的分子电流模型

Figure 8. Molecular current model of a magnetic ring and current-carrying coil with multiple layers

下载: 全尺寸图片幻灯片

通过模型可知，线圈等效电流i₃为

${i_3} = \frac{{{N_{\mathrm{c}}}I}}{{{L_{\mathrm{c}}}({R_{\mathrm{c}}} - {r_{\mathrm{c}}})}} \text{，}$

(11)

式中：N_c为线圈匝数，I为实际通入电流.

同永磁环分子电流模型推导步骤相同，线圈与永磁环“1”面之间的轴向磁力为

$\begin{split} & {F_{{\mathrm{c}}1}} = \frac{{{\mu _0}{i_3}{i_2}{r_{\mathrm{m}}}}}{{4{\text{π}}}}\int_0^{2{\text{π}}} {\int_0^{2{\text{π}}} {\int_0^{{L_{\mathrm{c}}}} {\int_0^{{L_{\mathrm{m}}}} {\int_{{r_{\mathrm{c}}}}^{{R_{\mathrm{c}}}} {} } } } } \\ &\quad \frac{{v({L_{\mathrm{c}}} + {L_{\mathrm{cg}}} - b + h)\cos ({\theta _{\text{P}}} - {\theta _{\text{Q}}})}}{{{r_{{\mathrm{cm}}}}^3}}{\text{d}}{\theta _{\text{P}}}{\text{d}}{\theta _{\text{Q}}}{\text{d}}b{\text{d}}h{\text{d}}v , \end{split}$

(12)

$\begin{split} & {r_{{\text{cm}}}} = [{( - v\cos \;{\theta _{\text{P}}} + {R_2}\cos \;{\theta _{\text{Q}}})^2} + \\ &\quad {( - v\sin \;{\theta _{\text{P}}} + {R_2}\sin \;{\theta _{\text{Q}}})^2} + {({L_{\mathrm{c}}} - b + {L_{\mathrm{cg}}} + h)^2}]^{1 /2}, \end{split}$

(13)

式中：r_cm为线圈上分子电流环与磁环上分子电流环的标量距离；db、dv分别为线圈上假想电流微元的厚度和径向厚度，dh为永磁体“1”面上假想电流微元的长度.

同理可求得线圈与面“2”之间的轴向磁力F_c2，因此，系统中单个线圈与永磁环间的轴向磁力为

${F_{{\mathrm{c}}12}} = {F_{{\mathrm{c}}1}} + {F_{{\mathrm{c}}2}} .$

(14)

同理可得另一线圈对永磁环的轴向磁力F_c21，线圈对磁环的总轴向磁力为

${F_{\mathrm{c}}} = {F_{{\mathrm{c}}12}} + {F_{{\mathrm{c}}21}} .$

(15)

通过3.1、3.2节的推导，系统的总轴向磁力为

${F_{\textit{z}}} = {F_{\mathrm{m}}} + {F_{\mathrm{c}}} .$

(16)

3.3 简化表达式

为得到显式表达式用于刚度建模和动力学方程求解，需要对式（16）进行近似处理. 式（16）在零刚度点z₀附近是连续的，可以对式（16）在z₀进行关于轴向位移z的三阶泰勒展开，如式（17）.

${F_{\textit{z}}} = {k_0} + {k_1}({\textit{z}} - {{\textit{z}}_0}) + {k_2}{({\textit{z}} - {{\textit{z}}_0})^2} + {k_3}{({\textit{z}} - {{\textit{z}}_0})^3} + o({\textit{z}}) \text{，}$

(17)

式中： ${k_0} \;=\; {\left. {{F_{\textit{z}}}} \right|_{{\textit{z}} = {{\textit{z}}_0}}},{k_1} \;=\; {\left. {{k_{\textit{z}}}} \right|_{{\textit{z}} \;=\; {{\textit{z}}_0}}},{k_2} \;=\; {\left. {\dfrac{1}{2}\dfrac{{\partial {k_{\textit{z}}}}}{{\partial {\textit{z}}}}} \right|_{{\textit{z}} \;=\; {{\textit{z}}_0}}}, {k_3} \;= {\left. {\dfrac{1}{6}\dfrac{{{\partial ^2}{k_{\textit{z}}}}}{{\partial {{\textit{z}}^2}}}} \right|_{{\textit{z}} = {{\textit{z}}_0}}}$ ，k_z为轴向位移， $o({\textit{z}})$ 是z大于3阶的高阶小项，在这里进行省略.

将理论解与−8 A电流时的三阶泰勒展开曲线求得的解进行比较，如图9所示，计算结果与三阶泰勒展开结果吻合较好，误差小于10%，可以基本满足隔振系统工程实际需求.

图 9 对比结果

Figure 9. Comparative results

下载: 全尺寸图片幻灯片

4. 动力学分析

4.1 动力学方程

假设隔振系统承受额定隔振质量m时，支撑平台的静平衡位置为零刚度点z₀. 此时，隔振器可以看作为单自由度系统，在平衡位置的静力学方程为

${k_0} + {k_1}({\textit{z}} - {{\textit{z}}_0}) + {k_2}{({\textit{z}} - {{\textit{z}}_0})^2} + {k_3}{({\textit{z}} - {{\textit{z}}_0})^3} - mg = 0 .$

(18)

为便于表达和计算，设置零刚度点z₀为坐标原点0，式（18）可以简化为

${k_0} + {k_1}{\textit{z}} + {k_2}{{\textit{z}}^2} + {k_3}{{\textit{z}}^3} - mg = 0 \text{.}$

(19)

当底座受到激励时，隔振系统的动力学简化模型，如图10所示，q为激励幅值，c为系统阻尼. 对其进行动力学分析，在激励频率ω的外部位移激励z_q=qcos ωt作用下，隔振系统在基座处的动力学方程为

图 10 动力学模型

Figure 10. Dynamic model

下载: 全尺寸图片幻灯片

$m\ddot {\textit{z}} + c\dot {\textit{z}} + {k_0} + {k_1}{\textit{z}} + {k_2}{{\textit{z}}^2} + {k_3}{{\textit{z}}^3} - mg = - m{\ddot {\textit{z}}_q} \text{.}$

(20)

由式（21）可知，系统为高阶非线性时变系统，通过四阶龙格库塔方法求得系统相平面图，如图11所示，随着时间t的不断增加，所有轨线最后都收敛于中间的一个极限环处，则该系统是稳定的.

图 11 相平面图

Figure 11. Phase plane

下载: 全尺寸图片幻灯片

在动力学求解中，以往的做法通常是省略k₂项，只保留一次项和三次项系数，但本系统轴向力两侧不对称，二次项系数不能省略. 为更好地进行对比，引入一个线性弹簧刚度k_s，所以线性弹簧的固有频率为 ${\omega _{\mathrm{n}}} = \sqrt {{{{k_{\mathrm{s}}}} / m}}$ ，将其引入到无量纲方程中. 因此，运动方程可以无量纲化为

$\hat {\textit{z}}'' + 2\zeta \hat {\textit{z}}' + \eta \hat {\textit{z}} + \alpha {\hat {\textit{z}}^2} + \gamma {\hat {\textit{z}}^3} = {\beta ^2}\cos\; \beta \tau \text{，}$

(21)

式中：求导是对无量纲时间τ进行的， $\hat {\textit{z}}$ 为无量纲相对位移幅值， $\xi$ 为无量纲阻尼比， $\;\beta$ 为无量纲频率比， $\eta 、\alpha 、\gamma$ 为无量纲刚度系数.

$\begin{split} \hat {\textit{z}} = \frac{{\textit{z}}}{q},\zeta = \frac{c}{{2m{\omega _{\mathrm{n}}}}},\beta = \frac{\omega }{{{\omega _{\mathrm{n}}}}}, \eta = \frac{{{k_1}}}{{{k_{\mathrm{s}}}}},\alpha = \frac{{{k_2}q}}{{{k_{\mathrm{s}}}}}, \gamma = \frac{{{k_3}{q^2}}}{{{k_{\mathrm{s}}}}}, \end{split}$

$\tau = {\omega _{\mathrm{n}}}t$ .

通常谐波平衡法（HB）只考虑基波项进行求解，无法考虑二次非线性项，如果想要正确反映各刚度系数的影响，就需要取更多的谐波项，求解比较困难. 所以本文采用增量谐波平衡法（IHB）求解，考虑了二次和三次非线性项的影响^[17-18].

设 ${\hat{{\textit{z}}}}_{0}、{\beta }_{0}$ 是式（21）的解，则其临近点可表示为 $\hat{{\textit{z}}}={\hat{{\textit{z}}}}_{0} + \Delta \hat{{\textit{z}}}，\beta ={\beta }_{0} + \Delta \beta$ ，其中， $\Delta \hat{{\textit{z}}}、\Delta \beta$ 为增量. 令 $\beta \tau = \phi$ ，将 $\hat{{\textit{z}}}、\beta$ 代入式（21），并略去高阶小量后得到增量方程为

$\begin{split} & \beta _0^2\Delta \hat {\textit{z}}'' + 2\xi {\beta _0}\Delta \hat {\textit{z}}' + \left( {\eta + 2\alpha {{\hat {\textit{z}}}_0} + 3\gamma \hat {\textit{z}}_0^2} \right)\Delta \hat {\textit{z}} = \\ &\quad R - \left( {2{\beta _0}{{\hat {\textit{z}}''}_0} + 2\xi {{\hat {\textit{z}}'}_0} - 2{\beta _0}\cos \;\phi } \right)\Delta \beta , \end{split}$

(22)

式中： $R = \beta _0^2\cos \;\phi - \left( {\beta _0^2{{\hat {\textit{z}}''}_0} + 2\xi {\beta _0}{{\hat {\textit{z}}'}_0} + \eta {{\hat {\textit{z}}}_0} + \alpha \hat {\textit{z}}_0^2 + \gamma \hat {\textit{z}}_0^3} \right)$ .

设式（22）的解为

${\hat {\textit{z}}} = \sum\limits_{n = 0}^N {\left( {{a_n}\cos \;n\phi + {b_n}\sin \; n\phi } \right)} \text{，}$

(23)

$\Delta \hat {\textit{z}} = \sum\limits_{n = 0}^N {\left( {\Delta {a_n}\cos \;n\phi + \Delta {b_n}\sin \;n\phi } \right)} .$

(24)

将式（23）、（24）代入式（22）中，求得以 $\Delta {a}_{n}、 \Delta {b}_{n}$ 的代数方程组，取项数N=5，求得a₁，…，a_n、b₁，…，b_n的值.

4.2 位移传递率

一般计算位移传递率 T的方法是忽略高次谐波项的影响，只取基波项，由质量位移幅值与基座位移幅值之比计算，存在一定误差. 在实际测量质量加速度时，波形含有高次谐波项，系统在13.5 Hz下采集的加速度图像如图12所示，所以，使用隔振质量位移响应对外部激励的均方根比来定义位移传递率.

图 12 质量与基座在13.5 Hz下的加速度图像

Figure 12. Acceleration image of mass and base at 13.5 Hz

下载: 全尺寸图片幻灯片

通过4.1节求解，幅值 $\hat{{\textit{z}}}$ 可表示为

$\hat{{\textit{z}}}={\displaystyle \sum _{n=0}^{5}({a}_{n}\mathrm{cos}\;n\phi + {b}_{n}\mathrm{sin}\;n\phi })={A}_{n}\mathrm{cos}(\beta \tau + {\theta }_{n}) \text{，}$

(25)

式中： ${A}_{{n}}=\sqrt{{a}_{{n}}{}^{2} + {b}_{{n}}{}^{2}}$ ， ${\theta }_{n}=\mathrm{arctan} \left(-\dfrac{{b}_{n}}{{a}_{n}}\right)$ .

基座无量纲位移为 $\hat{{\textit{z}}}_q$ =cos βτ，所以质量的无量纲位移 $\hat{{\textit{z}}}_m$ 为 $\hat{{\textit{z}}}$ 和 $\hat{{\textit{z}}}_q$ 的和，如式（26）.

${\hat{{\textit{z}}}}_{m}={\hat{{\textit{z}}}}_{q} + \hat{{\textit{z}}}=\mathrm{cos}\;\beta \tau + {A}_{n}\mathrm{cos}(\beta \tau + {\theta }_{n}).$

(26)

故位移传递率T表示为

$T = \frac{{{\mathrm{RMS}}({{\hat {\textit{z}}}_m}(t))}}{{{\mathrm{RMS}}({{\hat {\textit{z}}}_q}(t))}}.$

(27)

4.3 系统隔振效率分析

由于存在二次项与三次项等非线性项，该隔振系统表现出较为明显的非线性. 为了进一步分析该系统的隔振特性，通过数值仿真对比了不同电流下激励幅值q=1 mm与q=2 mm时系统的传递率曲线，如图13.

图 13 不同电流下的位移传递率

Figure 13. Simulated displacement transmissibility with different currents

下载: 全尺寸图片幻灯片

在数值仿真过程中隔振质量m=0.4 kg，磁环在闭合线圈中运动将会产生涡流阻尼. 但在本文中，动磁环的运动在线圈的外部，产生涡流阻力与主动控制的电磁力不在一个数量级，建模计算时将其忽略不计. 系统阻尼主要为直线轴承摩擦阻尼与空气阻尼构成，采用自由振动法，将系统阻尼估算为常数c.

由图13中可以看出：当q=1 mm时系统固有频率与最大传递率均随着负电流的增大而减小，整个系统表现出较强的线性特征；当q=2 mm时，系统固有频率与最大传递率同样随着电流的增大而减小，传递率曲线向右弯曲，表现出典型的渐硬弹簧特性.

5. 实验验证

搭建隔振系统实验测量系统，如图14所示. 隔振系统刚性垂直连接在VT6300-40型振动平台上，振动台的外部激励信号由VCS振动控制器和SA3K型功率放大器控制. 在振动平台和承载质量块上分别布置型号为1A213E和1A113E的加速度传感器.

图 14 试验系统示意图.

Figure 14. Test setup

下载: 全尺寸图片幻灯片

加载不同电流信号，测量所设计的非线性隔振系统的位移传递率. 隔振质量为0.4 kg，加载正弦激励信号，基座位移幅值为1 mm，进行线性扫频实验. 忽略其他部分的影响，系统阻尼近似为常数.

隔振器位移传递率实验测量结果如图15所示. 当电磁线圈通入负电流时，可以降低起始隔振频率，改善隔振性能，且随着负电流的增大，隔振性能改善越明显；当通入−8 A电流时，系统固有频率从10.62 Hz降低到8.64 Hz，起始隔振频率从14.86 Hz降低到12.00 Hz，传递率峰值从18.57 dB降到16.81 dB. 在此区间内可以通过改变电流值实现对不同振源的适应，实验测量结果与计算分析结果基本吻合.

图 15 不同电流下的位移传递率实验结果.

Figure 15. Experimental results of displacement transmissibility with different currents

下载: 全尺寸图片幻灯片

6. 结　论

本文基于电磁线圈嵌套永磁体的结构，提出一种高静-低动刚度的新型可变刚度隔振系统. 通过数值仿真和实验分析了激励、电流等因素对系统位移传递率的影响规律.

由仿真得到了系统参数对其刚度特性的影响规律，确定了结构的尺寸参数. 运用分子电流法构建了系统的力学模型及动力学模型. 采用增量谐波平衡法求解位移传递率，考虑了传统方法分析二次非线性项影响较困难的问题.

搭建了实验测试系统，由结果可知，当通入−8 A电流时，系统固有频率相对于未通入电流降低18.64%，起始隔振频率降低19.25%，且位移传递率的峰值降低了1.76 dB，使隔振系统的隔振效果进一步提升.

图 1 驾驶模拟实验平台及附加设备

Figure 1. Driving simulator and relevant equipment

下载: 全尺寸图片幻灯片

图 2 试验道路的线形和路段编号

Figure 2. Experimental roads and segment number

下载: 全尺寸图片幻灯片

图 3 闭眼眼睑开度阈值的计算示意

Figure 3. Schematic diagram of eyelid closure threshold

下载: 全尺寸图片幻灯片

图 4 疲劳预警提示逻辑流程

Figure 4. Flowchart of driving drowsiness alarming

下载: 全尺寸图片幻灯片

图 5 分级疲劳预警提示显示

Figure 5. Drowsiness alarming displays for different levels

下载: 全尺寸图片幻灯片

表 1 一般疲劳与严重疲劳预警前后的变量描述性统计指标汇总

Table 1. Descriptive statistics summary of normal and serious levels of drowsiness warning in before and after the warning period

报警类型	变量	均值	标准差	最小值	最大值
一般疲劳 N = 69	P_{clos_pre}	0.174 5	0.115 6	0.015 6	0.574 4
	P_{clos_post}	0.152 5	0.106 5	0	0.608 9
	M_{avg_pre}	0.222 4	0.577 6	–2.331 9	1.480 8
	M_{avg_post}	0.254 0	0.336 8	–0.618 5	0.946 4
	L_{sd_pre}	0.457 0	0.320 1	0.106 1	2.020 4
	L_{sd_post}	0.327 9	0.241 8	0.036 0	1.490 0
	D_{area_pre}	11.714 4	34.631 3	0	224.340 9
	D_{area_post}	3.140 0	10.091 8	0	60.062 8
严重疲劳 N = 96	P_{clos_pre}	0.359 9	0.131 7	0.063 3	0.713 3
	P_{clos_post}	0.315 0	0.142 7	0.062 2	0.757 8
	M_{avg_pre}	0.172 5	0.537 7	–3.381 4	0.861 0
	M_{avg_post}	0.191 4	0.399 3	–0.916 3	1.328 4
	L_{sd_pre}	0.433 0	0.242 8	0.125 1	1.668 7
	L_{sd_post}	0.488 7	0.826 4	0.091 7	8.138 6
	D_{area_pre}	7.994 6	37.169 4	0	351.571 1
	D_{area_post}	5.623 3	17.266 1	0	117.700 1

下载: 导出CSV

表 2 疲劳预警提示前后响应指标的Wilcoxon signed-rank配对检验

Table 2. Drowsiness alarming before and after the paired Wilcoxon signed-rank test of the responsive variables

报警类型	变量	均值	标准差	V值	P值
一般疲劳 N = 69	P_{clos_paired}	0.022 0	0.096 5	1 448.5	0.150 4
	M_{avg_paired}	–0.031 6	0.544 7	1 187.0	0.904 8
	L^*_{sd_paired}	0.129 1	0.313 6	1 879.0	0.000 1
	D^*_{area_paired}	8.574 4	33.605 3	1 064.0	0.002 1
严重疲劳 N = 96	P^*_{clos_paired}	0.044 9	0.138 6	3 153.5	0.000 1
	M_{avg_paired}	–0.018 9	0.560 8	2 383.0	0.842 1
	L_{sd_paired}	–0.055 7	0.811 5	2 644.0	0.248 9
	D_{area_paired}	2.371 4	34.455 4	1 520.0	0.477 0
注：*表示该指标在提示前后存在显著差异.

下载: 导出CSV

参考文献(21)

国务院滨保高速天津" 10•7”特别重大交通事故调查组. 滨保高速天津" 10•7”特别重大交通事故调查报告[R]. 北京:中华人民共和国国务院, 2012

国务院包茂高速陕西延安" 8•26”特别重大道路交通事故调查组. 包茂高速延安" 8•26”特大道路交通事故调查报告[R]. 北京: 中华人民共和国国务院, 2013

国务院办公厅. 国务院关于加强道路交通安全工作的意见 [EB/OL].（2012-07-27） [2017-02-05]. http://www.gov.cn/zwgk/2012-07/27/content_2193042.htm

中华人民共和国工业和信息化部. 关于进一步提高大中型客货车安全技术性能加强车辆《公告》管理和注册登记管理工作的通知[EB/OL].（2011-12-31）[2017-02-05]. http://www.miit.gov.cn/n1146295/n1652858/n1652930/n3757018/c3757341/content.html

中华人民共和国交通运输部. 关于加强道路运输车辆动态监管工作的通知[EB/OL].（2011-03-19）[2017-02-05]. http://www.chinasafety.gov.cn/Contents/Channel_5330/2011/0413/128291/asset000010003301211_0_1302663969013.html

G7智慧物联网公司. 中国物流大数据报告[EB/OL].（2017-01-13）[2017-02-05]. http://www.zqcn.com.cn/hongguan/201701/13/c491618.html

WIERWILLE W, WREGGIT S, KIRN L, et al. Research on vehicle-based driver status/performance monitoring; development, validation, and refinement of algorithms for detection of driver drowsiness[R]. Washington D. C.: US National Hightway Traffic Safety Adiministration, 1994

HORREY J, NOY I, FOLKARD S, et al. Research needs and opportunities for reducing the adverse safety consequences of fatigue[J]. Accident Analysis & Prevention, 2011, 43(2): 591-594

王福旺,王宏. 长途客车驾驶员疲劳状态脑电特征分析[J]. 仪器仪表学报,2013,34(5): 1146-1152 doi: 10.3969/j.issn.0254-3087.2013.05.027

WANG Fuwang, WANG Hong. EEG characteristic analysis of coach bus drivers in fatigue state[J]. Chinese Journal of Science Instrument, 2013, 34(5): 1146-1152 doi: 10.3969/j.issn.0254-3087.2013.05.027

程如中,赵勇,戴勇,等. 基于Adaboost方法的车载嵌入式疲劳驾驶预警系统[J]. 北京大学学报(自然科学版),2012,48(5): 719-726

CHENG Ruzhong, ZHAO Yong, DAI Yong, et al. An on-board embedded driver fatigue warning system based on Adaboost method[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2012, 48(5): 719-726

张希波,成波,冯睿嘉. 基于方向盘操作的驾驶人疲劳状态实时检测方法[J]. 清华大学学报(自然科学版),2010,50(7): 1072-1076

ZHANG Xibo, CHENG Bo, FENG Ruijia. Real-time detection of driver drowsiness based on steering performance[J]. Journal of Tsinghua University (Science and Technology), 2010, 50(7): 1072-1076

谭小强. 基于DSP技术的疲劳驾驶预警系统中车道偏离识别方法的研究[D]. 西安: 长安大学, 2010

胥川,王雪松,陈小鸿. 无侵入测量指标的驾驶疲劳检测性能评估[J]. 西南交通大学学报,2014,49(4): 720-726 doi: 10.3969/j.issn.0258-2724.2014.04.025

XU Chuan, WANG Xuesong, CHEN Xiaohong. Evaluating performance of non-intrusive indicators on drowsy driving detection[J]. Journal of Southwest Jiaotong University, 2014, 49(4): 720-726 doi: 10.3969/j.issn.0258-2724.2014.04.025

叶柠,孙宇舸. 基于EEG小波包子带能量比的疲劳驾驶检测方法[J]. 东北大学学报(自然科学版),2012,33(8): 1088-1092

YE Ning, SUN Yuge. A fatigue driving detection method based on wavelet packet sub-band energy ratio of EEG[J]. Journal of Northeastern University (Natural Science), 2012, 33(8): 1088-1092

赵晓华,许士丽,茉建,等. 基于ROC曲线的驾驶疲劳脑电样本熵判定阈值研究[J]. 西南交通大学学报,2013,48(1): 178-183 doi: 10.3969/j.issn.0258-2724.2013.01.028

ZHAO Xiaohua, XU Shili, RONG Jian, et al. Discriminating threshold of driving fatigue based on the electroencephalography sample entropy by receiver operating characteristic curve analysis[J]. Journal of Southwest Jiaotong University, 2013, 48(1): 178-183 doi: 10.3969/j.issn.0258-2724.2013.01.028

王连震,裴玉龙. 基于贝叶斯网络的驾驶疲劳程度识别模型[J]. 城市交通,2014,12(3): 66-74

WANG Lianzhen, PEI Yulong. Driving fatigue recognition model based on bayesian network[J]. Urban Transport of China, 2014, 12(3): 66-74

李家文,成波. 驾驶人状态适应式疲劳预警方法的研究[J]. 汽车工程,2011,33(8): 694-700

LI Jiawen, CHENG Bo. A study on the adaptive warning method for driving persons fatigue[J]. Automotive Engineering, 2011, 33(8): 694-700

胥川,王雪松,陈小鸿,等. 基于决策树的驾驶疲劳等级分析与判定[J]. 同济大学学报(自然科学版),2015,43(1): 75-81

XU Chuan, WANG Xuesong, CHEN Xiaohong, et al. Driver drowsiness level analysis and prediction based on decision tree[J]. Journal of Tongji University (Natural Science), 2015, 43(1): 75-81

DINGES D F, GRACE R. PERCLOS: A valid psychophysiological measure of alertness as assessed by psychomotor vigilance[R]. Washington D. C.: US Federal Highway Administration, 1998

胥川. 疲劳驾驶行为特征及提示有效性研究[D]. 上海: 同济大学, 2014

肖献强,殷延杰,王家恩. 基于个体特性的驾驶行为操纵模式建模方法[J]. 中国机械工程,2016,27(19): 2681-2686,2692 doi: 10.3969/j.issn.1004-132X.2016.19.021

XIAO Xianqiang, YIN Yanjie, WANG Jiaen. Driving behavior operation pattern modeling method based on individual characteristics[J]. China Mechanical Engineering, 2016, 27(19): 2681-2686,2692 doi: 10.3969/j.issn.1004-132X.2016.19.021

施引文献

期刊类型引用(7)

1.	Fengxiang Guo，Yiwen Zhou，Xiaoan Wang，Wenxuan Li，Jing Cai. Literature review of driving fatigue research based on bibliometric analysis. Journal of Traffic and Transportation Engineering(English Edition). 2024(06): 1401-1419 . 必应学术
2.	许永生，李丽丽，吴尤荻，支锦亦，向泽锐，张润芝. 基于注意力分配的高铁操纵台界面布局优化设计. 西南交通大学学报. 2022(02): 401-409 . 本站查看
3.	尉志礼，胥川，王雪松，王岳川. 基于隐马尔可夫模型的无侵入疲劳等级判别. 公路交通科技. 2022(05): 112-121 . 百度学术
4.	王朝侠，廖小菊. 基于用户需求的车载疲劳预警产品交互设计. 包装工程. 2022(14): 90-99 . 百度学术
5.	徐陶祎，张翼. 基于小波能量熵的司机疲劳驾驶行为响应算法. 计算机仿真. 2021(11): 133-137 . 百度学术
6.	李海舰，赵国强，杨艳芳，常鑫，赵晓华. 急刹预警作用下的车辆时空图特性. 华南理工大学学报(自然科学版). 2020(07): 76-84 . 百度学术
7.	王艳丽，杨烁，吴兵. 基于驾驶模拟的本科实验教学内容设计. 实验室研究与探索. 2020(11): 171-174 . 百度学术

其他类型引用(10)

附加材料(0)

访问统计

点击查看大图

图(5) / 表(2)

计量

文章访问数: 589
HTML全文浏览量: 288
PDF下载量: 33
被引次数: 17

1. 隔振系统结构与工作原理
1.1 结　构
1.2 工作原理
2. 结构参数设计
3. 力学模型
3.1 磁环间轴向磁力
3.2 线圈与磁环间的轴向磁力
3.3 简化表达式
4. 动力学分析
4.1 动力学方程
4.2 位移传递率
4.3 系统隔振效率分析
5. 实验验证
6. 结　论

疲劳预警分级提示下的驾驶行为响应特征

doi: 10.3969/j.issn.0258-2724.20180254

作者简介: 胥川（1987—），男，讲师，博士，研究方向为交通安全、数据科学等，E-mail： xuchuan@swjtu.edu.cn

通讯作者: 王雪松（1977—），男，教授，研究方向为交通安全、交通统计分析、交通规划、驾驶模拟器应用等，E-mail： wangxs@tongji.edu.cn

计量

出版历程

Driver Behavior Response to Drowsiness Alarming at Different Levels

1. 隔振系统结构与工作原理

1.1 结 构

1.2 工作原理

2. 结构参数设计

3. 力学模型

3.1 磁环间轴向磁力

3.2 线圈与磁环间的轴向磁力

3.3 简化表达式

4. 动力学分析

4.1 动力学方程

4.2 位移传递率

4.3 系统隔振效率分析

5. 实验验证

6. 结 论

期刊类型引用(7)

其他类型引用(10)

计量

出版历程

目录

1. 隔振系统结构与工作原理

1.1 结 构

1.2 工作原理

2. 结构参数设计

3. 力学模型

3.1 磁环间轴向磁力

3.2 线圈与磁环间的轴向磁力

3.3 简化表达式

4. 动力学分析

4.1 动力学方程

4.2 位移传递率

4.3 系统隔振效率分析

5. 实验验证

6. 结 论

作者简介:
胥川（1987—），男，讲师，博士，研究方向为交通安全、数据科学等，E-mail： xuchuan@swjtu.edu.cn

通讯作者:
王雪松（1977—），男，教授，研究方向为交通安全、交通统计分析、交通规划、驾驶模拟器应用等，E-mail： wangxs@tongji.edu.cn

1.1 结　构

6. 结　论

1.1 结　构

6. 结　论