汽车工程

参数	数值
初始容积 $V 0$ /m³	0.000 79
有效承载面积 $A e$ /m²	0.003 5
初始压力 $p 0$ /MPa	1.48
实验加载 $F$ /N	4 900

参数	数值
初始容积 $V 0$ /m³	0.000 79
有效承载面积 $A e$ /m²	0.003 5
初始压力 $p 0$ /MPa	1.48
实验加载 $F$ /N	4 900

参数	仿真	实验	相对误差
滞回曲线面积/（N·mm）	0.8 Hz	10 mm	12 930	13 913	7.6%
12 mm	21 022	21 534	2.4%
1.0 Hz	10 mm	17 924	19 152	6.8%
12 mm	25 872	28 275	9.2%
滞回曲线刚度/（N·mm^-1）	0.8 Hz	10 mm	78.05	84.25	8.0%
12 mm	73.94	80.18	8.4%
1.0 Hz	10 mm	86.87	96.74	7.6%
12 mm	72.72	73.48	1.1%

参数	仿真	实验	相对误差
滞回曲线面积/（N·mm）	0.8 Hz	10 mm	12 930	13 913	7.6%
12 mm	21 022	21 534	2.4%
1.0 Hz	10 mm	17 924	19 152	6.8%
12 mm	25 872	28 275	9.2%
滞回曲线刚度/（N·mm^-1）	0.8 Hz	10 mm	78.05	84.25	8.0%
12 mm	73.94	80.18	8.4%
1.0 Hz	10 mm	86.87	96.74	7.6%
12 mm	72.72	73.48	1.1%

参数	数值
簧上质量 $m s$ /kg	1 730
质心至前轴距离 $a$ /mm	1 189
质心至后轴距离 $b$ /mm	1 548
绕 $x$ 轴转动惯量 $I x$ /（kg·m²）	700
绕 $y$ 轴转动惯量 $I y$ /（kg·m²）	3 480
绕 $z$ 轴转动惯量 $I z$ /（kg·m²）	3 710
前轮距 $t f$ /mm	1 644
后轮距 $t r$ /mm	1 635
质心高度 $h$ /mm	704
车轮质量 $m w$ /kg	64.5
轮毂电机质量 $m d$ /kg	31.5
轮胎刚度 $k t$ /（N·mm^-1）	218 900
减振器阻尼系数 $C s$ /（N·s·m^-1）	1 200

参数	数值
簧上质量 $m s$ /kg	1 730
质心至前轴距离 $a$ /mm	1 189
质心至后轴距离 $b$ /mm	1 548
绕 $x$ 轴转动惯量 $I x$ /（kg·m²）	700
绕 $y$ 轴转动惯量 $I y$ /（kg·m²）	3 480
绕 $z$ 轴转动惯量 $I z$ /（kg·m²）	3 710
前轮距 $t f$ /mm	1 644
后轮距 $t r$ /mm	1 635
质心高度 $h$ /mm	704
车轮质量 $m w$ /kg	64.5
轮毂电机质量 $m d$ /kg	31.5
轮胎刚度 $k t$ /（N·mm^-1）	218 900
减振器阻尼系数 $C s$ /（N·s·m^-1）	1 200

车辆响应	侧倾角加速度	空气弹簧力
60 km/h	RMS	0.094 8	0.093 5	534.6	570.1
max	0.287	0.291	1 114	1 117
min	-0.253	-0.286	-1 197	-1 205
75 km/h	RMS	0.137 2	0.129 6	604.7	641.9
max	0.357	0.379	1 289	1 294
min	-0.380	-0.410	-1 376	-1 386
90 km/h	RMS	0.189	0.177	685.0	752.1
max	0.416	0.432	1 416	1 460
min	-0.561	-0.573	-1 701	-1 716

车辆响应	侧倾角加速度	空气弹簧力
60 km/h	RMS	0.094 8	0.093 5	534.6	570.1
max	0.287	0.291	1 114	1 117
min	-0.253	-0.286	-1 197	-1 205
75 km/h	RMS	0.137 2	0.129 6	604.7	641.9
max	0.357	0.379	1 289	1 294
min	-0.380	-0.410	-1 376	-1 386
90 km/h	RMS	0.189	0.177	685.0	752.1
max	0.416	0.432	1 416	1 460
min	-0.561	-0.573	-1 701	-1 716

车辆响应	俯仰角加速度	空气弹簧力
3 m/s²	RMS	0.162	0.149	410.2	452.6
max	1.033	1.035	645.2	670.1
min	-0.563	-0.560	-1 318	-1 319
2.5 m/s²	RMS	0.146	0.136	366.8	399.9
max	0.899	0.901	543.7	591.2
min	-0.521	-0.521	-1 238	-1 235
2 m/s²	RMS	0.116	0.108	281.7	309.1
max	0.689	0.690	385.6	391.3
min	-0.440	-0.443	-1 025	-1 042

车辆响应	俯仰角加速度	空气弹簧力
3 m/s²	RMS	0.162	0.149	410.2	452.6
max	1.033	1.035	645.2	670.1
min	-0.563	-0.560	-1 318	-1 319
2.5 m/s²	RMS	0.146	0.136	366.8	399.9
max	0.899	0.901	543.7	591.2
min	-0.521	-0.521	-1 238	-1 235
2 m/s²	RMS	0.116	0.108	281.7	309.1
max	0.689	0.690	385.6	391.3
min	-0.440	-0.443	-1 025	-1 042

参数类别	系统参数
系统控制输入	空气弹簧质量流量
系统可测输出	俯仰角加速度、侧倾角加速度、悬架动挠度、簧上质量加速度、簧下位移
不可测输出	轮胎动载荷
可测干扰	制动力矩、前轮转角
不可测干扰	路面激励

参数类别	系统参数
系统控制输入	空气弹簧质量流量
系统可测输出	俯仰角加速度、侧倾角加速度、悬架动挠度、簧上质量加速度、簧下位移
不可测输出	轮胎动载荷
可测干扰	制动力矩、前轮转角
不可测干扰	路面激励

车辆响应	被动悬架	PID控制	MPC控制
前悬架动挠度/mm	未考虑分数阶	RMS	7.8	6.8	7.3
max	9.6	8.6	0
min	-28.3	-25.3	-15.7
考虑分数阶	RMS	7.9	6.9	7.4
max	10.0	8.8	0
min	-28.9	-25.8	-23.6
后悬架动挠度/mm	未考虑分数阶	RMS	9.3	6.9	7.9
max	41.2	30.8	16.7
min	-10.4	-7.05	0
考虑分数阶	RMS	9.7	7.2	9.2
max	44.1	33.0	19.2
min	-11.5	-8.2	0
俯仰角加速度/（rad·s^-2）	未考虑分数阶	RMS	0.135	0.114 5	0.082 97
max	0.903 7	0.903 6	0.903 9
min	-0.521 9	-0.464 5	-0.453 4
考虑分数阶	RMS	0.120	0.101 7	0.062 28
max	0.911 2	0.910 8	0.911 4
min	-0.512 1	-0.457 3	-0.455 3

车辆响应	被动悬架	PID控制	MPC控制
前悬架动挠度/mm	未考虑分数阶	RMS	7.8	6.8	7.3
max	9.6	8.6	0
min	-28.3	-25.3	-15.7
考虑分数阶	RMS	7.9	6.9	7.4
max	10.0	8.8	0
min	-28.9	-25.8	-23.6
后悬架动挠度/mm	未考虑分数阶	RMS	9.3	6.9	7.9
max	41.2	30.8	16.7
min	-10.4	-7.05	0
考虑分数阶	RMS	9.7	7.2	9.2
max	44.1	33.0	19.2
min	-11.5	-8.2	0
俯仰角加速度/（rad·s^-2）	未考虑分数阶	RMS	0.135	0.114 5	0.082 97
max	0.903 7	0.903 6	0.903 9
min	-0.521 9	-0.464 5	-0.453 4
考虑分数阶	RMS	0.120	0.101 7	0.062 28
max	0.911 2	0.910 8	0.911 4
min	-0.512 1	-0.457 3	-0.455 3

车辆响应	被动悬架	PID控制	MPC 控制
前左悬架动挠度/mm	未考虑分数阶	RMS	14.0	10.9	7.8
max	42.7	31.7	17.3
min	-23.8	-19.8	-20.7
考虑分数阶	RMS	16.5	12.8	10.7
max	43.4	32.1	26.8
min	-24	-19.9	-30.5
前右悬架动挠度/mm	未考虑分数阶	RMS	14.5	11.2	10.5
max	40.9	30.5	26.8
min	-24.2	-20.2	-21.0
考虑分数阶	RMS	17.3	13.6	11.2
max	41.9	31.1	30.8
min	-24.6	-20.5	-24.8
侧倾角加速度/（rad·s^-2）	未考虑分数阶	RMS	0.145 0	0.105 5	0.054 7
max	0.365 6	0.278 1	0.143 5
min	-0.377 5	-0.239 6	-0.144 1
考虑分数阶	RMS	0.125 8	0.089 3	0.049 39
max	0.380 5	0.287 9	0.145 8
min	-0.414 8	-0.254 6	-0.141 1

车辆响应	被动悬架	PID控制	MPC 控制
前左悬架动挠度/mm	未考虑分数阶	RMS	14.0	10.9	7.8
max	42.7	31.7	17.3
min	-23.8	-19.8	-20.7
考虑分数阶	RMS	16.5	12.8	10.7
max	43.4	32.1	26.8
min	-24	-19.9	-30.5
前右悬架动挠度/mm	未考虑分数阶	RMS	14.5	11.2	10.5
max	40.9	30.5	26.8
min	-24.2	-20.2	-21.0
考虑分数阶	RMS	17.3	13.6	11.2
max	41.9	31.1	30.8
min	-24.6	-20.5	-24.8
侧倾角加速度/（rad·s^-2）	未考虑分数阶	RMS	0.145 0	0.105 5	0.054 7
max	0.365 6	0.278 1	0.143 5
min	-0.377 5	-0.239 6	-0.144 1
考虑分数阶	RMS	0.125 8	0.089 3	0.049 39
max	0.380 5	0.287 9	0.145 8
min	-0.414 8	-0.254 6	-0.141 1

悬架动挠度	急减速工况	双移线工况
max/mm	3	2	9	9
min/mm	-10	-10	-6	-5
RMS/mm	3.73	3.67	3.72	3.59

悬架动挠度	急减速工况	双移线工况
max/mm	3	2	9	9
min/mm	-10	-10	-6	-5
RMS/mm	3.73	3.67	3.72	3.59

基于分数阶的空气弹簧建模及电动汽车主动悬架控制研究

PDF下载

冯桂珍 ¹^,² , 赵东鹏 ² , 李韶华 ¹

汽车工程 | 2024,46(7): 1282-1293

收起

汽车工程 | 2024, 46(7): 1282-1293

基于分数阶的空气弹簧建模及电动汽车主动悬架控制研究

全屏

冯桂珍¹^,², 赵东鹏², 李韶华¹

作者信息

^1. 石家庄铁道大学，省部共建交通工程结构力学行为与系统安全国家重点实验室，石家庄 050043

^2. 石家庄铁道大学机械工程学院，石家庄 050043

通讯作者:

赵东鹏，硕士研究生，E-mail：zdpbaoding@163.com。

Research on Air Spring Modeling Based on Fractional Order and Electric Vehicle Active Suspension Control

Guizhen Feng¹^,², Dongpeng Zhao², Shaohua Li¹

Affiliations

^1. Shijiazhuang Tiedao University，State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures，Shijiazhuang 050043

^2. School of Mechanical Engineering，Shijiazhuang Tiedao University，Shijiazhuang 050043

出版时间: 2024-07-25 doi: 10.19562/j.chinasae.qcgc.2024.07.015

文章导航

摘要

收起

电控空气悬架系统（electrically controlled air suspension， ECAS）具有调节悬架刚度和车身高度的功能，可有效改善车辆乘坐舒适性和操纵稳定性。以某乘用车ECAS为对象，利用分数阶理论描述橡胶气囊的黏弹性阻尼特性，考虑等效阻尼及滞回特性对其热力学模型进行了优化，结果与实验数据吻合良好，验证了优化后的空气弹簧模型的精确性。在此基础上，考虑车辆纵横向动力学特性与Dugoff轮胎模型，建立14自由度整车ECAS动力学模型，提出模型预测（model predictive control， MPC）主动悬架控制方法，以可测变量为控制器输入，实现直线及转向行驶工况下的主动控制。仿真与整车台架实验研究表明，分数阶修正模型可以很好地反映ECAS变刚度特性，基于MPC的主动悬架控制策略能实时调整空气弹簧刚度，控制车身姿态，有效改善电动汽车行驶时的平顺性与稳定性。论文的研究方法为车辆悬架系统建模及主动控制提供了一种新思路。

关键词

空气悬架 / 热力学 / 黏弹性阻尼 / 分数阶 / 模型预测控制

Abstract

收起

Electrically controlled air suspension （ECAS） has the function of adjusting suspension stiffness and body height， which can effectively improve vehicle ride comfort and handling stability. Taking a passenger car ECAS as an example， the viscoelastic damping characteristics of rubber airbag are described by fractional theory， and the thermodynamic model is optimized considering the equivalent damping and hysteretic characteristics， which is in good agreement with the experimental data， and the precision of the optimized air spring model is verified. On this basis， considering the longitudinal and lateral dynamic characteristics of the vehicle and the Dugoff tire model， a 14-degree-of-freedom vehicle ECAS dynamic model is established， and a Model Predictive Control （MPC） active suspension control method is proposed， with measurable variables as the input of the controller， to realize the active control under straight and turning driving conditions. Simulation and vehicle bench test show that the fractional correction model can well reflect the variable stiffness characteristics of ECAS， and the active suspension control strategy based on MPC can adjust the air spring stiffness in real time， control the body posture， and effectively improve the ride comfort and stability of the electric vehicle. The research method in this paper provides a new idea for vehicle suspension system modeling and active control.

Key words

air suspension / thermodynamics / viscoelastic damping / fractional order / model predictive control

引用本文

冯桂珍, 赵东鹏, 李韶华. 基于分数阶的空气弹簧建模及电动汽车主动悬架控制研究. 汽车工程, 2024 , 46 (7) : 1282 -1293 . DOI: 10.19562/j.chinasae.qcgc.2024.07.015

Guizhen Feng, Dongpeng Zhao, Shaohua Li. Research on Air Spring Modeling Based on Fractional Order and Electric Vehicle Active Suspension Control[J]. Automotive Engineering, 2024 , 46 (7) : 1282 -1293 . DOI: 10.19562/j.chinasae.qcgc.2024.07.015

正文

收起

前言

收起

随着电子技术的发展和制造成本的降低，电控技术在汽车行业得到了广泛应用，电控空气悬架系统（electrically controlled air suspension，ECAS）是其典型代表。作为当前最先进的悬架系统之一，ECAS能实时调节悬架刚度^［1］，具有更低的固有频率、更理想的非线性特性等^［2-3］，使得车辆能适应不同的路况、车辆负载和行驶速度，极大地改善了行驶性能，而如何准确、快速地调节空气弹簧刚度，成为当前研究的重点。

空气弹簧的刚度特性受气压、激励频率和幅值等诸多因素的影响，常见的空气弹簧的非线性模型有多项式模型及热力学简化模型。多项式模型一般采用三次多项式对空气弹簧刚度的幅频特性进行描述，冯霏等^［4］、盛云等^［5］都曾采用三次多项式，分别利用多尺度法、增量谐波平衡法等分析了汽车悬架系统的2自由度非线性振动系统的振动特性。尹万建等^［6］指出汽车系统中，空气弹簧刚度的频变现象并不明显，故可用三次多项式对空气弹簧刚度特性进行描述，但这类模型没有考虑橡胶气囊黏弹性特性的影响。热力学简化模型是应用最广的建模方法，将空气弹簧的力学模型分解为弹性力及黏性阻尼力两部分，其中黏性阻尼力存在非线性阻尼力以及惯性力，Zhu等^［7］基于Berg模型增加了对空气弹簧摩擦特性的描述，提出摩擦力由减振杆和气缸之间的相对运动以及橡胶材料颗粒产生。邬明宇等^［8］针对多腔室空气弹簧提出了一种考虑空气气囊热交换与阻尼滞回特性的力学模型，考虑滞回特性的精确解模型可以很好反映空气弹簧的动力学特性。Oda等^［9］研究建立了一种附加气室空气弹簧力学模型，具有一定的频率适应性。Docquier等^［10］研究了不同热交换系数对附加气室空气弹簧内部特性的影响，分析了传热对动态刚度的影响。陈俊杰等^［11］将分数导数Maxwell中建立空气弹簧模型，用来描述橡胶气囊的黏弹性特性。目前，针对空气弹簧力学特性及其模型的研究较多，而将空气弹簧应用于整车动力学模型控制悬架系统的研究尚不多见。

主动悬架控制一直是研究热点，常见的控制方法有天棚控制^［12-13］、最优控制^［14-16］、滑模控制^［17］及PID控制^［18］等。车辆在行驶过程中存有诸多不确定性，模型预测控制无须复杂的辨识过程，并且在优化过程利用实测信息不断矫正，对不确定性环境具有很强适应能力，因此，模型预测控制方法适合用于ECAS。同时，由于空气弹簧具有响应慢、时滞长等特点，目前多针对汽车行驶时的俯仰运动、侧倾运动以及车身高度控制等问题进行研究。Huo等^［19］提出将抑制俯仰运动的控制器作为补充控制模块，在特定工况下开启，并基于天棚控制设计PID控制器控制俯仰角。Eltantawie 等^［20］将半车模型分成两个独立的1/4模型，对每个模型设计分散神经模糊控制，优化车辆的垂向响应与俯仰角。Chen等^［21］设计了一款有高流量软管的空气悬架系统，可改进货车在车道变化和稳态转弯时的操控性和稳定性。上述文献均取得了一定的成果，但多针对阻尼减振器^［22］、油气悬架^［23］等，将空气弹簧作为控制对象的研究尚不多见。此外，现有关于ECAS研究多侧重于空气弹簧力学模型优化和车身高度调节等，涉及车辆行驶时侧向和纵向的车身姿态控制的研究较为少见。

鉴于此，本文基于热力学模型推导空气弹簧的弹性力，利用分数阶描述橡胶气囊的黏弹性阻尼力，以此修正热力学模型与工程实际中空气弹簧的差异，通过实验验证了模型的准确性。在此基础上，建立电动汽车整车ECAS垂向动力学模型，并考虑纵横向动力学特性与Dugoff轮胎模型，针对直线与转向行驶工况，以降低俯仰角加速度与侧倾角加速度为控制目标，搭建整车模型预测（MPC）控制器，并通过仿真分析和整车ECAS实验验证了控制策略的有效性。

1 空气弹簧力学模型与实验验证

收起

1.1 空气弹簧力学模型

空气弹簧在工作过程中，气囊内气体为多变过程，气体质量固定时，理想气体多变过程方程有：

p e V e k = p 0 V 0 k

（1）

式中：

p 0

为静止状态时空气弹簧内部的绝对压力，Pa；

V 0

为静止时空气弹簧气囊内容积，m³；

p e

为任意时刻空气弹簧内部的绝对压力，Pa；

V e

为任意时刻空气弹簧气囊内容积，m³；

k

为气体多变指数，绝热过程时

k

取1.4。

空气弹簧在充放气调节刚度的过程中，气体在气囊内流动时会有做功和焓变，一般充放气的时间较短，可视为绝热过程，根据热力学第一定律有^［24］：

W · + h i n q i n - h o u t q o u t = U a s

（2）

式中：

W

为空气弹簧做功；

h i n

、

h o u t

分别为流入、流出空气弹簧的比焓；

q i n

、

q o u t

分别为流入、流出空气弹簧的质量流量；

U a s

表示空气弹簧的内能。

结合理想气体状态方程

p e V e = m R T e

，式（2）可表达为

C v R (p ˙ e V e + p e V e ·) = k R T k - 1 q m - p e V e ·

（3）

式中：

R

为理想气体常数，

R

=287 J/（kg·K）；

T e

为气囊内气体温度；

q m

为质量流量，

q m = q i n - q o u t

；

C v

为定容比热容，

C v = R / (k - 1)

。

式（3）进一步整理为

p ˙ e V e = k R T q m - k p e V e ·

（4）

空气弹簧运动过程中，有效容积及气囊内压作用力为

V e · = A e ⋅ (z 2 · - z 1 ·) F k = (p e - p 0) ⋅ A e

（5）

式中：

A e

为空气弹簧有效承载面积，m²；

z 1

、

z 2

分别为簧上质量、簧下质量位移，m；

F k

为气囊内压作用力，N。

三次多项式模型、热力学模型没有考虑空气弹簧橡胶气囊材料的黏弹性阻尼特性，然而空气弹簧在工作频率范围内，由橡胶气囊黏弹性产生的阻尼效应会影响空气弹簧的力学响应，因此采用分数阶微分对橡胶气囊黏弹力进行表达^［25］，以此修正热力学模型与实际空气弹簧之间的差异，空气弹簧力学模型对比如图1所示。

分数阶微积分的定义形式很多，但都较为复杂，一般采用传递函数G（S）逼近分数阶模型^［26］。

G S = K ∏ k = 1 N s + ω ˜ k s + ω k

（6）

ω k = ω b ω u (2 k - 1 + α) / N ω ˜ k = ω b ω u (2 k - 1 - α) / N

（7）

ω u = ω h / ω b K = ω h α

（8）

式中：

K

为增益；

ω ˜ k

、

ω k

为极点或零点。

F c = d α' x (t) d t α' C F a s = F k + F c

（9）

式中：

F c

为黏弹性力；

α'

为空气弹簧垂向变形量对时间阶数；

C

为分数阶等效阻尼系数；

F a s

为空气弹簧力。

1.2 空气弹簧力学特性实验

为了验证基于分数阶的空气弹簧力学模型的准确性，需要对空气弹簧的动特性进行实验。某空气悬架减振总成由空气弹簧和CDC减振器串联组成，相关参数如表1所示，将减振总成安装在振动台上，在减振总成上盖板与实验台架之间安装力传感器，如图2所示。调节空气弹簧高度至226 mm，向空气弹簧充气达到初始压强后关闭气源。对空气弹簧施加正弦激励，振幅分别为10、12、15 mm等，频率为0.2~1.5 Hz，记录力与位移数据。本实验目的是测试减振器总成中空气弹簧的力学性能，为了消除与其串联的CDC减振器的影响，实验中还测量了相同工况下CDC减振器不通电时的阻尼力。

空气弹簧刚度和分数阶阻尼系数的计算方法如图3所示。空气弹簧的刚度

K

可由滞回曲线的最高点与最低点得出，滞回曲线的面积

S = ∫ F d x

为空气弹簧垂向振动过程中消耗的能量。

在消除CDC减振器的影响后，得到空气弹簧的滞回曲线，调整仿真模型参数与实验工况一致，限于篇幅，通过仿真得出在载荷为4 900 N时，仅呈现振幅为10、12 mm，频率

f

为0.8、1.0 Hz的正弦激励时的空气弹簧滞回曲线，如图4所示，仿真与实验结果及其相对误差如表2所示。

由图4和表2可知，仿真得到的空气弹簧分数阶力学模型滞回曲线与实验结果重合度较高，其中，滞回曲线面积

S

与实验相对误差在10%以内，刚度

K

与实验结果相对误差在10%以内，同时，与三次多项式模型及传统热力学模型相比，经分数阶修正后的力学模型具有滞回特性，说明此模型能有效反映空气弹簧的力学特性，具有较高的精度，可满足工程实际应用需求。

2 基于ECAS的电动汽车建模与仿真

收起

2.1 纵横垂向耦合动力学模型

本文主要聚焦直线和转向行驶时的车身姿态，须考虑车辆纵横向运动与整车垂向力学模型的耦合，反映直线行驶急减速及转向工况下，车辆的纵向、俯仰与侧倾运动以及悬架动挠度、车轮转动等运动响应。考虑车辆纵向速度及前轮转角，将轮毂电机作为车轮的附加质量，建立电动汽车纵-横-垂向耦合的14自由度整车ECAS动力学模型，如图5所示。

（1）车辆纵-横向动力学方程

车辆纵向、侧倾和横摆运动微分方程为

(F x f l + F x f r) c o s δ - (F y f l + F y f r) s i n δ + F x r l + F x r r = m a x

（10）

(F x f l + F x f r) s i n δ + (F y f l + F y f r) c o s δ + F y r l + F y r r =

m (v ˙ y + v x γ)

（11）

[(F x f r - F x f l) c o s δ + (F y f l - F y f r) s i n δ] a 2 + [(F x f r + F x f l) s i n δ + (F y f l + F y f r) s i n δ] a -

(F x r l - F x r r) b 2 - (F y r l + F x r r) b = I z γ ˙

（12）

式中：

δ

为前轮转角；

a x

为纵加速度；

v x

、

v y

分别代表车辆纵向速度、侧向速度；

m

为整车质量；

F x

、

F y

分别为轮胎纵向力与侧向力；

a

、

b

为质心至前后轴的距离；

γ

为横摆角速度；

I z

为整车绕

z

轴方向的转动惯量；下标

f l

、

f r

、

r l

、

r r

分别代表车辆左前轮、右前轮、左后轮和右后轮。

轮胎模型采用Dugoff模型，轮胎受力如图6所示，纵、侧向轮胎力为

F x = C x σ 1 + σ f (λ) F y = C α σ 1 + σ f (λ)

（13）

式中：

σ

为滑移率；

C x

、

C α

分别为轮胎纵向刚度、侧向刚度；

λ

与路面和轮胎间摩擦因数

μ

有关；

α

为侧偏角；函数

f (λ)

为

f (λ) = (2 - λ), λ < 1 1, λ ≥ 1

（14）

λ = μ F z (1 + σ) 2 (C x σ) 2 + (C α t a n α) 2

（15）

各车轮垂直载荷：

F z i = m g b 2 l ∓ m a y h b t f l - m a x h 2 l (i = f l, f r) F z i = m g a 2 l ∓ m a y h a t r l + m a x h 2 l (i = r l, r r)

（16）

同时，各车轮的纵向力与侧向力均不能超出轮胎的附着椭圆方程：

(F x i u F z i) 2 + (F y i u F z i) 2 ≤ 1 (i = f l, f r, r l, r r)

（17）

各车轮侧偏角：

α i = δ - a r c t a n v y + a γ v x ∓ t f γ / 2 (i = f l, f r) α i = - a r c t a n v y - b γ v x ∓ t r γ / 2 (i = r l, r r)

（18）

各车轮纵向速度：

v f l = (v x ∓ t f γ / 2) c o s δ + (v y + α γ) s i n δ (i = f l, f r) v r l = v x ∓ t r γ / 2 (i = r l, r r)

（19）

车轮旋转方程：

J w ω ˙ i = T d i - T b i - F x i R (i = f l, f r, r l, r r)

（20）

制动时滑移率：

σ i = w i R v i - 1 (i = f l, f r, r f, r r)

（21）

式中：

t f

、

t r

为前后轴轮距；

l

为前后轴距离；

F z

为轮胎垂向力；

h

为质心高度；

v

为车轮纵向速度；

J w

为车轮转动惯量；

T b

、

T d

分别为制动力矩与驱动力矩；

ω

为车轮角速度；

R

为车轮半径。

（2）车辆垂向动力学方程

由牛顿第二定理和车辆系统动力学原理，簧下质量运动可表示为

(m w i + m d i) ⋅ x ¨ w i = - k t i (x w i - q t i) + C s i (x ˙ i - x ˙ w i) - F a s i

(i = f l, f r, r l, r r)

（22）

式中：

m w

为车轮质量；

m d

为轮毂电机质量；

x

为车身4个悬架连接处的位移；

x w

为簧下质量的位移；

k t

为轮胎刚度；

q t

为路面激励；

C s

为减振器阻尼系数；

F a s

为空气弹簧力。

车身垂向、俯仰和侧倾运动表达式为

m s x ¨ s = F a s f l - C s f l (x ˙ f l - x ˙ w f l) + F a s f r - C s f r (x ˙ f r - x ˙ w f r) + F a s r l - C s r l (x ˙ r l - x ˙ w r l) + F a s r r - C s r r (x ˙ r r - x ˙ w r r) - m s g I y θ · · = b [F a s r l - C s r l (x ˙ r l - x ˙ w r l) + F a s r r - C s r r (x ˙ r r - x ˙ w r r)] - a [F a s f l + F a s f r - C s f l (x ˙ f l - x ˙ w f l) + F a s f r - C s f r (x ˙ f r - x ˙ w f r)] + M θ f l + M θ r l + M θ f r + M θ r r I x φ · · = c [- F a s f l + C s f l (x ˙ f l - x ˙ w f l) - F s r l + C s r l (x ˙ r l - x ˙ w r l)] + d [- F a s f r + C s f r (x ˙ f r - x ˙ w f r) - F s f r + C s r r (x ˙ r r - x ˙ w r r)] + M φ f l + M φ r l + M φ f r + M φ r r M θ i = F x i (h + x i - x w i) M φ i = F y i (h + x i - x w i) (i = f l, f r, r l, r r)

（23）

式中：

m s

为簧上质量；

x s

为簧上质量的位移；

θ

为俯仰角；

φ

为侧倾角；

c

、

d

为

1 / 2

轮距；

F x

为制动力矩不为零时地面对车轮的摩擦力；

M θ

、

M φ

分别为俯仰力矩与侧倾力矩。

2.2 仿真分析

为了分析所建立的整车ECAS模型能否有效反映车辆在行驶过程中的车身姿态变化，选取直线急减速和双移线转向工况分别进行仿真，车辆参数如表3所示。

（1）直线急减速工况

车辆匀速直线行驶，第7 s开始减速，车速由120减速到60 km/h，车辆纵向减速度a分别为3、2.5、2 m/s²，仿真结果如图7所示，相同工况下分数阶与未考虑分数阶的ECAS对车辆响应的影响，如表4所示。

由图7可知，随着纵向减速度的增大，车轮纵向力、俯仰角加速度和空气弹簧垂向力都增加，车辆的平顺性与稳定性随之降低。在开始减速一瞬间，俯仰角加速度达到第1个峰值，后逐渐稳定，这是因为急减速过程中制动力过大，导致车身俯仰角快速变化，在结束制动后俯仰角加速度出现了第2个峰值，失去制动力后，由于空气悬架的作用，车身趋于平稳。由表4可知，在纵向减速度分别为3、2.5、2 m/s²时，相比于未考虑分数阶，分数阶模型的空气弹簧力均方根值分别增加10.2%、9.0%、9.7%，俯仰角加速度均方根值降低8.1%、6.9%、6.9%，这是因为考虑了分数阶后，空气弹簧具有一定的阻尼特性，对减速时产生的俯仰具有一定的抑制。

（2）双移线转向工况

双移线转向行驶仿真，车速分别取60、75和90 km/h，第3 s开始转向，转向过程约8 s，仿真结果如图8所示。相同工况下考虑分数阶与未考虑分数阶的ECAS对车辆响应的影响对比如表5所示。

由图8可知，在双移线的过程中随着车速的增加，车侧向力、侧倾角加速度和空气弹簧力均增大，车辆的平顺性与稳定性降低。在开始变道时，侧倾角加速度达到第1个峰值，这是由于车身侧向速度增加引起的，变更回原车道后，车身逐渐平稳。由表5可知，在车速分别取60、75、90 km/h时，考虑分数阶后，空气弹簧力均方根值增加6.7%、6.2%、9.7%，侧倾角加速度降低1.4%、5.5%、6.4%，与直线减速行驶工况的仿真结果一致。

综上可见，随着减速度与车速的增加，车辆的平顺性与稳定性降低。基于第1节建立分数阶空气弹簧模型，通过纵-横-垂向整车ECAS动力学模型分析，可以很好地反映车辆在直线行驶时不同减速度与不同速度转向行驶过程中的车身姿态及空气弹簧垂向力的变化，进一步验证了所建立的空气弹簧模型的有效性。在此基础上，可对其进行控制策略的设计与验证。

3 模型预测（MPC）控制器设计及实验验证

收起

模型预测控制（MPC）的控制动作是在每一采样瞬间通过求解有限时域开环的最优控制问题而获得，其控制流程如图9所示。

3.1 模型预测（MPC）控制器设计

为便于MPC控制器的求解，对整车ECAS模型进行线性离散化处理，处理后的状态空间方程可写为

X · = A X + B u Y = C X + D u

（24）

其中：

X = (x w f l, x w f r, x w r l, x w r r, x ˙ w f l, x ˙ w f r, x ˙ w r l, x ˙ w r r, x s, x ˙ s, φ, φ ˙, θ, θ ˙) Y = (x f l - x w f l, x f r - x w f r, x r l - x w r l, x r r - x w r r, x s, x ˙ s, x ¨ s, φ, φ ¨, θ, θ ¨) u = (F a s f l, F a s f r, F a s r l, F a s r r, q f l, q f r, q r l, q r r, F x r (h + x s - q r), F x f (h + x s - q f), M f l, M r l, M f r, M r r)

（25）

式中 A 、 B 、 C 、 D 为状态矩阵。

预测模型的功能是根据被控对象的历史信息

[u k - j, y (k - j)]

，预测其未来输入

u k + j

与未来输出

y (k + j)

（

k

为当前时刻，

j ≥ 1

），具体输入输出参数如表6所示。

滚动优化是在每个时间段内的寻求最优，优化的性能指标在内容上也是实时变化的，预测控制的优化过程是反复在线进行的，即滚动优化。优化性能指标为

J = ∑ i = 1 p q i (y (k + i) - w (k + i)) 2 + ∑ j = 1 m r j u 2 (k + j - 1)

（26）

式中

q i

和

r j

分别为输出加权矩阵和受控输入加权矩阵，可根据期望设定相应权重系数。

在预测模型的基础上考虑被控对象的约束条件并对其进行处理是预测控制的一个主要特征，综合考虑车辆悬架系统的性能要求以及实际行驶情况，需要对滚动优化设定约束条件，即同时限定悬架动挠度不可超出限制（乘用车一般为7~9 cm），轮胎动载荷不能超出车身重力，且车辆滑移率与侧偏角均在合理区间，纵、侧向轮胎力满足式（17）的附着椭圆约束。

被控对象

t

时刻接收到预测控制器控制信号，被控系统

t + 1

时刻的实际输出

y (t + 1)

与预测模型预测输出

y ˜ (t + 1)

不相等，构成预测误差：

e (t + 1) = y (t + 1) - y^(t + 1) y ˜ p = y^p + h ⋅ e (t + 1)

（27）

式中：

y ˜

为预测系统输出；

h

为误差校正矢量。

3.2 模型预测（MPC）控制器仿真分析

在Simulink中搭建ECAS整车模型及MPC控制器，假设制动力矩恒定，将俯仰力矩作为状态空间外部输入。仿真工况为：车辆以120 km/h速度匀速直线行驶，在第7 s开始减速，减速过程持续6.4 s，车速减到60 km/h，得到两种工况下不同控制器对ECAS整车模型时域响应的影响，如图10所示。对比MPC控制与PID控制下的ECAS与被动空气悬架的前悬架动挠度、后悬架动挠度与俯仰角加速度的均方根值和极值，如表7所示。

由图10和表7可知，在对ECAS施加MPC控制后，前悬架动挠度、后悬架动挠度与俯仰角加速度的极值相比于被动悬架分别下降了18%、54%和-0.02%，相比于PID分别降低8%、42%和-0.06%，均方根值相比于被动悬架分别下降了9%、6%和48%，相比于PID分别降低7%、27%和39%，同时，相较于未考虑分数阶修正的ECAS，考虑分数阶修正的被动悬架、PID控制、MPC控制的车辆前悬架动挠度均方根值增大1.3%、1.5%、1.4%，俯仰角加速度减小11.2%、11.3%、25.1%。说明MPC控制器可有效减小悬架动挠度，改善车身姿态，提高舒适性。由于制动力矩为阶跃信号，MPC控制器针对俯仰角加速度极值的优化并不明显。

为进一步验证MPC控制器的有效性，进行双移线工况的仿真分析：车速为60 km/h，仿真时长15 s，在第3 s进行双移线转向，8 s后换回原车道，得到两种工况下不同控制器对ECAS时域响应的影响，如图11所示。MPC控制、PID控制与被动空气悬架的前悬架动挠度、后悬架动挠度与俯仰角加速度的均方根值和极值，如表8所示。

由图11和表8可知，施加MPC控制的前左悬架动挠度、前右悬架动挠度与侧倾角加速度的极值相比于被动悬架分别下降了38%、27%和66%，相比于PID分别降低17%、1%和45%，均方根值相比于被动悬架分别下降了24%、23%和61%，相比于PID分别降低2%、0%和45%，相较于未考虑分数阶修正的ECAS，考虑分数阶修正的被动悬架、PID控制、MPC控制的车辆前左悬架动挠度增大17.9%、17.4%、37.2%，侧倾角加速度减小13.7%、15.4%、9.8%。PID控制在急减速与双移线中车辆的响应指标除急减速时前、后悬架动挠度均方根略优于MPC控制外，其余指标均劣于MPC控制，并且由图10（d）和图11（d）可看出MPC控制下的电磁阀开关次数明显低于PID控制。由此可见，对ECAS施加MPC控制后，有效降低了车辆在进行急减速与双移线时的悬架动挠度、俯仰角加速度和侧倾角加速度，并且优化效果优于PID控制，可以使车辆更平稳，保证乘坐舒适性。

3.3 实验验证

3.3.1 实验系统及控制流程

ECAS实验台由空气弹簧与CDC减振器组成，主要包括上位机及测控系统、配电控制柜、ECAS悬架实验台，如图12所示。

本实验台架省去簧下质量部分，从PC机导出簧下质量位移

x w

作为外部激励，由伺服电动缸实现作动，控制流程图如图13所示，通过加速度传感器采集ECAS实验台簧上质量加速度信号和车身位移信号，发送至MPC5744P控制器，进行信号转换以及积分运算，经过滤波后作为控制器输入。ECAS控制器计算出理想质量流量，通过PWM处理后作为的控制信号。在Labview搭建上位机，可采集实验中被动与施加控制效果时的传感器信号，并对采集的信号进行示波显示及保存。

3.3.2 实验结果分析

由Simulink导出被动悬架的动挠度，采样间隔0.07 s，急减速工况及双移线工况的采样时间分别为7和12 s，每次实验时保持空气弹簧高度在同一固定值，且CDC减振器不通电，以便实现基于空气弹簧的控制效果的验证。采集悬架动挠度，并计算控制前与控制后的极值进行对比分析，以验证控制算法对ECAS的控制效果。

由图14和表9可以看出在急减速及双移线工况下，施加控制后的悬架动挠度有了一定改善，其中，急减速工况中控制后的悬架动挠度极大值下降了1 mm，均方根值下降了2%，双移线工况中控制后的悬架动挠度极小值下降了1 mm，均方根值下降了3.5%。证明本文所提出的MPC控制在整车ECAS控制中的有效性，可改善车辆在行驶时的平顺性和舒适性。

4 结论

收起

基于热力学理论和分数阶建立了空气弹簧力学特性模型，通过实验验证了模型的准确性。在此基础上，考虑车辆纵横向动力学特性与Dugoff轮胎模型，建立了整车ECAS动力学模型，针对直线与转向行驶工况，设计了整车模型预测（MPC）控制器，并进行了仿真分析和整车ECAS实验验证，研究结果表明：

（1）建立的分数阶空气弹簧热力学模型，可以很好地表达空气弹簧的力学特性，修正了空气弹簧的黏弹性阻尼特性，实验证明该模型可以较精确描述空气弹簧的力学特性。基于所提出的分数阶空气弹簧模型，通过整车ECAS纵-横-垂向动力学模型分析，可以很好地反映车辆在行驶过程中的车身姿态变化，进一步表明所建立的空气弹簧模型的可行性。

（2）以降低车身俯仰角加速度与侧倾角加速度为控制目标，设计提出了一种基于整车ECAS的MPC控制器，与PID控制和被动空气悬架控制相比，可有效改善车辆悬架动挠度、俯仰角加速度和侧倾角加速度。其中，急减速工况下，前悬架动挠度、后悬架动挠度与俯仰角加速度的极值相比于被动悬架分别下降了18%、54%和-0.02%，均方根值相比于被动悬架分别下降了9%、6%和48%，双移线工况下，前左悬架动挠度、前右悬架动挠度与侧倾角加速度的极值相比于被动悬架分别下降了38%、27%和66%，均方根值相比于被动悬架分别下降了24%、23%和61%。

（3）整车台架实验表明，采用MPC控制器控制的ECAS在车辆响应方面相较于未施加控制的空气悬架有一定优化作用，证明了所设计的控制策略应用于ECAS的有效性。

基金

收起

国家自然科学基金(U22A20246)
河北省高等学校科学技术研究项目(JZX2023014)
河北省省级科技计划资助项目(225676162GH)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

陈黎卿，郑爽，曹恺. 考虑EPS与ESP耦合作用的ECAS系统控制策略［J］. 汽车工程， 2017， 39（11）： 1261-1267.

CHEN Liqing， ZHENG Shuang， CAO Kai. Control strategy of ECAS system considering the coupling effect of EPS and ESP ［J］. Automotive Engineering， 2017，39 （11）： 1261-1267.

LV Y Y J.Kinematics & compliance analysis of double wishbone air suspension with frictions and joint clearances［J］.Mechanism and Machine Theory： Dynamics of Machine Systems Gears and Power Trandmissions Robots and Manipulator Systems Computer-Aided Design Methods， 2021， 156（1）.

王相平，李星，王剑锋，等. 中低速磁浮车辆空气弹簧动力学建模及其应用研究［J/OL］. 振动工程学报： 1-12. http：//kns.cnki.net/kcms/detail/32. 1349. TB. 20220812. 1607. 002.

WANG Xiangping， LI Xing， WANG Jianfeng， et al. Aerodynamic modeling and application research of air springs for medium and low speed maglev vehicles ［J/OL］. Journal of Vibration Engineering： 1-12. http：//kns.cnki.net/kcms/detail/32.1349.tb.2220812.1607.002.

冯霏，李凌轩，陈亚哲，等. 汽车非线性悬架系统共振特性的研究［J］. 东北大学学报（自然科学版）， 2010， 31（3）： 415-417，426.

FENG Fei， LI Lingxuan， CHEN Yazhe， et al. Study on resonance characteristics of automotive nonlinear suspension system ［J］. Journal of Northeastern University （Natural Science Edition）， 2010， 31 （3）： 415-417，426.

盛云，吴光强. 基于 IHBM 的汽车非线性悬架系统定量研究［J］.同济大学学报（自然科学版）， 2011， 39（3）：405-410.

SHENG Yun， WU Guangqiang. Quantitative study on automobile nonlinear suspension system based on IHBM ［J］. Journal of Tongji University （Natural Science Edition）， 2011，39 （3）： 405-410.

尹万建，韩鹰，杨绍普. 空气弹簧悬架系统在强迫振动下的动力学分析［J］. 中国公路学报， 2006（3）： 117-121.

YIN Wanjian， HAN Ying， YANG Shaopu. Dynamic analysis of air spring suspension system under forced vibration ［J］. China Journal of Highway and Transport， 2006（3）： 117-121.

ZHU H， YANG J， ZHANG Y， et al. Nonlinear dynamic model of air spring with a damper for vehicle ride comfort［J］. Nonlinear Dynamics， 2017， 89（4）： 1545-1568.

邬明宇，李雪冰，尹航，等. 双腔室空气弹簧动刚度理论模型与实验研究［J］. 振动工程学报， 2022， 35（4）： 834-840.

WU Mingyu， LI Xuebing， YIN Hang， et al. Theoretical model and experimental study on dynamic stiffness of double-chamber air spring ［J］. Journal of Vibration Engineering， 2022，35（4）：834-840.

ODA N， NISHIMURA S. Vibration of air suspension bogies and their design［J］. JSME International Journal， 2008， 13（55）：43-50.

DOCQUIER N， FISETTE P， JEANMART Hervé.Influence of heat transfer on railway pneumatic suspensions dynamics［J］. Iavsd Symposium on Dynamics of Vehicles on Roads & Tracks， 2009.

陈俊杰，郭孔辉，殷智宏，等. 囊式空气弹簧垂向刚度统一模型研究［J］. 机械工程学报， 2022，58（12）：180-187.

CHEN Junjie， GUO Konghui， YIN Zhihong， et al. Study on unified vertical stiffness model of bladder air spring ［J］. Journal of Mechanical Engineering， 2022，58（12）：180-187.

YI K， SONG B S. A new adaptive sky-hook control of vehicle semi-active suspensions［J］. Proceedings of the Institution of Mechanical Engineers Part D： Journal of Automobile Engineering， 1999， 213（3）：293-303.

杨艺，陈龙，汪若尘，等. 车辆半主动悬架广义天棚理论控制研究［J］. 振动与冲击， 2021（22）：040.

YANG Yi， CHEN Long， WANG Ruochen， et al. Study on generalized canopy theoretical control of vehicle semi-active suspension ［J］. Vibration and Impact， 2021（22）：040.

KOCK J， ROHDE A. Large potential in tuning and optimization of electronic air suspension systems with variable damping/encompassing a wide range of tuning parameters and limits［J］. VDI BERICHTE， 2003， 1791： 243-264.

DU H， SZE K Y， LAM J. Semi-active H∞ control of vehicle suspension with magneto-rheological dampers［J］. Journal of Sound & Vibration， 2005， 283（283）： 981-996.

BHANGAL N S， MAURYA A. Optimal control of vehicle active suspension system［J］. Journal of Automation and Control Engineering， 2018， 6（1）：22-26.

秦武，朱钢，上官文斌，等. 具有扰动观测器的汽车主动悬架滑模控制［J］. 振动工程学报， 2020， 33（1）：158-167.

QIN Wu， ZHU Gang， SHANGGUAN Wenbin， et al. Sliding mode control of automotive active suspension with disturbance observer ［J］. Journal of Vibration Engineering， 2020，33 （1）： 158-167.

程振扬，杨明，吴心杰，等. 基于BP神经网络PID控制的主动悬架仿真分析［J］. 农业装备与车辆工程，2022，60（2）：130-134.

CHENG Zhenyang， YANG Ming， WU Xinjie， et al. Simulation analysis of active suspension based on BP neural network PID control ［J］. Agricultural Equipment and Vehicle Engineering， 2022， 60（2）： 130-134.

HUO S， YU L， MA L， et al. Ride comfort improvement in post-braking phase using active suspension ［C］. International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. American Society of Mechanical Engineers， 2015， 57168： V006T10A072.

ELTANTAWIE M A. Decentralized neuro-fuzzy control for half car with semi-active suspension system［J］. International Journal of Automotive Technology， 2012， 13（3）：423-431.

CHEN Y， PETERSON A W， AHMADIAN M. Achieving anti-roll bar effect through air management in commercial vehicle pneumatic suspensions［J］. Vehicle System Dynamics， 2018.

李子先，潘世举，朱愿，等.基于状态反馈和预瞄前馈的智能车半主动悬架控制［J］.汽车工程， 2023， 45（5）：735-745.

LI Zixian， PAN Shiju， ZHU Yuan， et al. Semi-active suspension control of smart car based on state feedback and preview feedforward ［J］. Automotive Engineering， 2023，45 （5）： 735-745.

孙涛，张振东. 油气主动悬架不确定性对重载车辆平顺性的影响分析［J］.振动工程学报， 2007（6）：627-632.

SUN Tao， ZHANG Zhendong. Analysis of the influence of the uncertainty of hydro-pneumatic active suspension on the ride comfort of heavy-duty vehicles ［J］. Journal of Vibration Engineering， 2007（6）：627-632.

LIU H， LEE J C. Model development and experimental research on an air spring with auxiliary reservoir［J］. International Journal of Automotive Technology， 2011， 12（6）：839-847.

ALONSO A， GIL-NEGRETE N， NIETO J， et al. Development of a rubber component model suitable for being implemented in railway dynamic simulation programs［J］. Journal of Sound & Vibration， 2013， 332（12）：3032-3048.

TEPLJAKOV A， PETLENKOV E， BELIKOV J. FOMCON： a MATLAB toolbox for fractional-order system identification and control［J］. International Journal of Microelectronics & Computer Science， 2011.

2024年第46卷第7期

PDF下载

335

145

引用本文

BibTeX

文章信息

doi: 10.19562/j.chinasae.qcgc.2024.07.015

接收时间：2024-01-12
首发时间：2025-07-29
出版时间：2024-07-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2024-01-12
修回日期：2024-02-24

基金

国家自然科学基金(U22A20246)

河北省高等学校科学技术研究项目(JZX2023014)

河北省省级科技计划资助项目(225676162GH)

作者信息

^1. 石家庄铁道大学，省部共建交通工程结构力学行为与系统安全国家重点实验室，石家庄 050043

^2. 石家庄铁道大学机械工程学院，石家庄 050043

通讯作者:

赵东鹏，硕士研究生，E-mail：zdpbaoding@163.com。

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/qcygc/CN/10.19562/j.chinasae.qcgc.2024.07.015

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

车辆响应		侧倾角加速度		空气弹簧力
车辆响应		未考虑分数阶	考虑分数阶	未考虑分数阶	考虑分数阶
60 km/h	RMS	0.094 8	0.093 5	534.6	570.1
	max	0.287	0.291	1 114	1 117
	min	-0.253	-0.286	-1 197	-1 205
75 km/h	RMS	0.137 2	0.129 6	604.7	641.9
	max	0.357	0.379	1 289	1 294
	min	-0.380	-0.410	-1 376	-1 386
90 km/h	RMS	0.189	0.177	685.0	752.1
	max	0.416	0.432	1 416	1 460
	min	-0.561	-0.573	-1 701	-1 716

车辆响应		俯仰角加速度		空气弹簧力
车辆响应		未考虑分数阶	分数阶	未考虑分数阶	分数阶
3 m/s²	RMS	0.162	0.149	410.2	452.6
	max	1.033	1.035	645.2	670.1
	min	-0.563	-0.560	-1 318	-1 319
2.5 m/s²	RMS	0.146	0.136	366.8	399.9
	max	0.899	0.901	543.7	591.2
	min	-0.521	-0.521	-1 238	-1 235
2 m/s²	RMS	0.116	0.108	281.7	309.1
	max	0.689	0.690	385.6	391.3
	min	-0.440	-0.443	-1 025	-1 042

悬架动挠度	急减速工况		双移线工况
悬架动挠度	控制前	控制后	控制前	控制后
max/mm	3	2	9	9
min/mm	-10	-10	-6	-5
RMS/mm	3.73	3.67	3.72	3.59