CN109002599B - 基于田口实验的汽车平顺性优化分析方法 - Google Patents

Abstract

本发明公开了一种基于田口实验的汽车平顺性优化分析方法，包括以下步骤：S1、基于田口方法实验，确定需要进行优化的影响汽车平顺性的可控因子和噪声因子；S2、根据可控因子和噪声因子的个数和其水平数生成组合正交表的表头；S3、在Matlab/Simulink中对各可控因子与噪声因子的组合依次进行仿真实验，算出座椅垂向加速度均方根值；S4、对步骤S3中所得数据进行统计学分析与计算，得到信噪比的方差分析表，并计算出各可控因子对汽车平顺性及其稳健性的贡献率；可将仿真分析与田口方法相结合，分析各可控设计变量与随机干扰因素对平顺性的影响，并可优选出最佳参数组合，提高正常平顺性及其稳定性。

Description

基于田口实验的汽车平顺性优化分析方法

技术领域

本发明涉及汽车测试技术领域，特别涉及一种基于田口实验的汽车平顺性优化分析方法。

背景技术

汽车的平顺性主要是保证汽车在行驶过程中产生的振动和冲击环境对乘员舒适性的影响在一定界限之内。悬架是影响汽车平顺性的主要因素，众多专家学者从事悬架系统的研究、优化与控制工作，使汽车平顺性得到不断地提高。对于悬架的优化，现有技术中主要存在这么几种方式，一种是考虑了簧上质量变化、执行机构等延迟等因素，设计了一个可靠的模糊H_∞控制器，并采用四分之一车辆悬架模型验证了设计方法的有效性；另一种是应用改进的多目标遗传算法优化了某种油气悬架的参数；再一种是结合蒙特卡罗方法和6σ稳健性优化技术，对悬架参数进行多目标稳健性优化，改善了设计变量在最优解附近波动而引起平顺性变差的状况。但是，这些优化都将路面、载荷和车速等因素视为不变的，没有考虑随机干扰因素对汽车平顺性的影响，因而不能保证汽车平顺性在复杂工况下的稳健性。

发明内容

有鉴于此，本发明提供一种基于田口实验的汽车平顺性优化分析方法，可综合考虑悬架参数、座椅参数等可控设计变量和路面、载荷、车速等随机干扰因素等汽车平顺性的影响，将仿真分析与田口方法相结合，分析各可控设计变量与随机干扰因素对平顺性的影响，并可优选出最佳参数组合，提高正常平顺性及其稳定性。

本发明的基于田口实验的汽车平顺性优化分析方法，包括以下步骤：

S1、基于田口方法实验，确定需要进行优化的影响汽车平顺性的可控因子(即可控设计变量)和噪声因子(即随机干扰因素)，其中，可控因子包括前后悬架的刚度、阻尼和驾驶员座椅的刚度、阻尼，噪声因子包括路面等级、载荷和车速，同时确定每个可控因子和噪声因子的水平数；

S2、根据可控因子的个数和其水平数生成田口正交表作为内表；

根据噪声因子的个数和其水平数生成田口正交表为外表；

由内、外表的设计得到各可控因子与噪声因子的组合正交表的表头；

S3、在Matlab/Simulink中对各可控因子与噪声因子的组合依次进行仿真实验，得到座椅垂向振动加速度图，并算出座椅垂向加速度均方根值，并计算各可控因子组合下的响应均值和信噪比；

S4、对步骤S3中所得数据进行统计学分析与计算，得到信噪比的方差分析表，并计算出各可控因子对汽车平顺性及其稳健性的贡献率。

进一步，步骤S1中，每个可控因子的水平数为五，并以第三个水平数代表车型初始设计参数值，其余水平对应的值分别为初始值的倍数；

步骤S2中，选取L₂₅(5⁶)正交表为内表。

进一步，步骤S1中，每个噪声因子的水平数为三；步骤S2中，选取L₉(3³)正交表为外表。

进一步，步骤S1中还有步骤S0，进行整车平顺性模型及仿真实验，以确定车型是否需要平顺性优化。

进一步，步骤S0，具体为：

建立八自由度整车动力学模型；该模型中，车身被简化为刚体，八自由度分别为垂直位移(Z_c)、侧倾角(Φ)和俯仰角(θ)3个自由度，驾驶员座椅垂直位移(Z_s)1个自由度，4个非簧载质量分别具有1个垂直位移(Z₁、Z₂、Z₃、Z₄)自由度；

取系统的广义坐标为：Z^T＝[Z₁ Z₂ Z₃ Z₄ Z_c Z_s Φ θ]；

由拉格朗日方程和达朗贝尔原理，可得8自由度整车平顺性模型振动微分方程为

式中：M为8×8阶的系统质量矩阵；C为8×8阶的系统阻尼矩阵；K为8×8阶的系统刚度矩阵；K_t为8×4阶的轮胎刚度矩阵；Q为4×1阶的路面激励矩阵；Z为8×1阶的位移矩阵；

利用Matlab/Simulink进行仿真实验，模拟该车在某级路面，空载状况下，以特定速度匀速行驶时的整车振动状况，得到驾驶员座椅垂向加速度随时间变化的曲线图；并在Matlab中计算得到座椅垂向加速度均方根值，座椅垂向加速度均方根值越小，代表汽车平顺性越好，否则，代表车型的汽车平顺性需要优化。。

本发明的有益效果：本发明的基于田口实验的汽车平顺性优化分析方法，可综合考虑悬架参数、座椅参数等可控设计变量和路面、载荷、车速等随机干扰因素等汽车平顺性的影响，将仿真分析与田口方法相结合，分析各可控设计变量与随机干扰因素对平顺性的影响，并可优选出最佳参数组合，提高正常平顺性及其稳定性。

附图说明

下面结合附图和实施例对本发明作进一步描述。

图1为汽车四轮所受路面随机激励曲线图；

图2为整车八自由度振动动力学模型示意图；

图3为座椅垂向振动加速度图；

图4为本发明方法的步骤图。

具体实施方式

下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围内。

本发明提出一种基于田口实验的汽车平顺性优化分析方法。参照图1-图4，图1为汽车四轮所受路面随机激励曲线图；图2为整车八自由度振动动力学模型示意图；图3为座椅垂向振动加速度图，图4为本发明方法的步骤图，如图所示：

本发明的基于田口实验的汽车平顺性优化分析方法，在本方法的具体步骤进行前，要建立路面激励模型，以为具体步骤施行奠定基础；具体为：

车辆振动的输入路面不平度主要采用路面功率谱密度来描述其特性。这里采用滤波白噪声法建立路面激励时域模型，其单轮路面激励q(t)的时域模型为

式中：n_∞为路面空间截止频率，n_∞＝0.011m^-1；u为汽车行驶速度；n₀为参考空间频率，G_q(n₀)为路面不平度的几何平均值；n₀＝0.1m^-1；ω(t)为均值为0的单位白噪声随机信号；

同侧前后车轮所受激励存在迟滞性，单轮辙激励模型为

式中：α为与路面等级相关的路面常数；q₁、q₂分别为汽车同侧前、后车轮的路面激励；l为汽车前后轴的距离；Δ为前后车轮路面激励的输入延迟，Δ＝l/u；

左右两侧车轮轨迹具有异辙相干性，假设左、右两侧轮辙的白噪声输入分别ξ为ξ_x、ξ_y，则两侧轮辙相干性方程为

右侧轮辙白噪声输入模型为

式中：a₀、a₁、a₂、b₀、b₁、b₂为道路品质常数，a₀＝3.1815、a₁＝0.2063、a₂＝0.0108、b₀＝3.223、b₁＝0.59、b₂＝0.0327；x₁、x₂为中间状态变量；

本发明所研究的汽车主要行驶在A级、B级和C级路面。对于A级、B级和C级路面，其路面不平度的几何平均值G_q(n₀)和路面常数α的数值如表一所示。

表一路面参数G_q(n₀)与α的取值

以B级道路路面不平度作为车轮激励，取车速为u＝20m/s，在Matlab/Simulink中进行仿真实验，得到汽车四轮所受路面随机激励如图1所示。

以下进入步骤S0，整车平顺性模型及仿真，以确定车型是否需要平顺性优化，具体为：

经过合理简化，建立8自由度整车动力学模型如图2所示。该模型中，车身被简化为刚体，有质心的垂直位移(Z_c)、侧倾角(Φ)和俯仰角(θ)3个自由度，驾驶员座椅垂直位移(Z_s)1个自由度，4个非簧载质量分别具有1个垂直位移(Z₁、Z₂、Z₃、Z₄)自由度；

取系统的广义坐标为Z^T＝[Z₁ Z₂ Z₃ Z₄ Z_c Z_s Φ θ]，由拉格朗日方程和达朗贝尔原理，可得8自由度整车平顺性模型振动微分方程为

式中：M为8×8阶的系统质量矩阵；C为8×8阶的系统阻尼矩阵；K为8×8阶的系统刚度矩阵；K_t为8×4阶的轮胎刚度矩阵；Q为4×1阶的路面激励矩阵；Z为8×1阶的位移矩阵；

然后进行仿真分析，以表二中所示汽车的动力学参数取值进行说明；

表二汽车动力学参数说明及取值

利用Matlab/Simulink进行仿真实验，模拟该车在C级路面，空载状况下，以25m/s的速度匀速行驶时的整车振动状况，得到驾驶员座椅垂向加速度随时间变化的曲线如图3所示。在Matlab中计算得到座椅垂向加速度均方根值为0.6057m/s²，根据提供的垂向加速度均方根值与人主观感觉之间的关系，座椅垂向加速度均方根值偏大，对于所研究车型，驾驶员的舒适性不够好，需要进行平顺性优化。

步骤S1基于田口方法实验，确定需要进行优化的影响汽车平顺性的可控因子(即可控设计变量)和噪声因子(即随机干扰因素)，其中，可控因子包括前后悬架的刚度、阻尼和驾驶员座椅的刚度、阻尼，噪声因子包括路面等级、载荷和车速，同时确定每个可控因子和噪声因子的水平数；

田口方法实验是基于损失模型的稳健优化设计实验，其原理是利用设计参数和响应之间的非线性关系来实现产品的稳健性，其目的是寻找可控因素的最佳组合，使得噪声因素的影响尽可能地小。田口方法的基本工具是信噪比和正交表；正交实验设计方法只需要做少数试验就能全面反映出试验条件完全组合的内在规律，以确定参数的最佳组合，可以缩短设计时间，提高设计效率；

悬架及座椅的动力学参数对汽车乘坐舒适性的影响较为明显。因此，这里选取前后悬架的刚度、阻尼和驾驶员座椅的刚度、阻尼作为6个可控因子；取示例车型，前悬架偏频范围为1.00～1.45Hz，后悬架偏频范围为1.17～1.58Hz，结合已知参数，可以计算出悬架参数的取值范围。为了详细分析可控因子各水平的影响，每个因子均选择5种水平，如表三所示。其中，水平3指该车初始设计参数值，其余水平对应的值分别为初始值的倍数。根据可控因子个数及其水平数，选取L₂₅(5⁶)正交表为内表。

表三可控因子水平表

依据汽车平顺性的外界影响因素分析，选取路面等级、载荷和车速为3个噪声因子，其水平的选择根据该车型的常态行驶环境来确定，如表四所示。根据噪声因子的个数及其水平数，选取L₉(3³)正交表为外表。

表四噪声因子水平表

步骤S2，田口方法中，内正交表是可控因子各水平组合的安排，外正交表是噪声因子各水平组合的安排，由内外正交表的设计得到正交表的表头。

步骤S3、在Matlab/Simulink中对各种可控因子与噪声因子的组合依次进行仿真实验，得到座椅垂向振动加速度图，编程算出座椅垂向加速度均方根值(单位：mm/s²)，并计算各可控因子组合下的响应均值和信噪比，得到实验结果如表五所示。

表五田口方法实验正交表

步骤S4、对表五的数据进行统计学分析与计算，得到信噪比的方差分析表，如表六所示。通过贡献率的计算，可以得出各可控因子对平顺性及其稳健性的影响程度。由表六可知：各可控因子中，前悬架刚度对该车平顺性及其稳健性的贡献率最高，达到71.11％；后悬架刚度次之，为10.26％；前、后悬架的阻尼和座椅的刚度、阻尼对应的贡献率较小。因此，对于该车型，悬架的刚度对汽车平顺性影响显著，其设计与加工精度应着重考虑。

汽车平顺性的客观评价常以座椅的垂向加速度均方根值作为指标。这里以悬架动挠度和轮胎动载荷为约束条件，以驾驶员座椅的垂向加速度均方根值最小为优化目标，其值越小，汽车平顺性越好。田口稳健设计要求目标函数及其在噪声因素干扰下的波动值越小越好。考虑望小特性，基于损失函数法的信噪比为

式中：yi为对应的目标函数响应值。

表六信噪比方差分析表

噪声因子的影响趋势如表七所示。由表七可知：汽车平顺性随路面不平度的增加而变差，随载荷的增加而改善，随速度的加快而变差。各噪声因子对汽车平顺性的影响程度由高到低依次为：路面、车速、载荷。

表七噪声因子反应表单位:mm/s²

最终优化结果分析：

稳健优化前后各可控因子的水平值及汽车在以上9种工况下的座椅垂向加速度均方根值(单位：mm/s²)如表八所示。由表八可知：优化后，汽车驾驶员座椅的垂向加速度均方根值在以上9种噪声因子的组合工况下均有所减小，其变化范围和均值均变小，其中，均值由345.1下降至267.9，减小了22.37％；信噪比变大；说明驾驶员座椅的舒适性显著提升，且其随噪声因子变化的波动范围变小，即提高了汽车的平顺性及其稳健性。

仿真实验结果显示：汽车各轮胎相对动载荷与各悬架动行程于优化前后的反应具有一致性。在以上9种工况下，右后轮相对动载荷均方根值gi与右后悬架动行程均方根值fi(单位：mm)如表九所示。由表九可知：优化后，右后悬架动行程均方根值与右后轮相对动载荷均方根值在各工况下均减小，汽车操纵稳定性提升。信噪比增大，说明悬架动行程与轮胎相对动载荷在各工况下的波动范围也变小，汽车操纵稳定性的稳健性也有所提高。由此，验证了田口方法汽车平顺性分析与稳健优化的可行性。

表八优化前后可控因子水平及响应分析

表九优化前后汽车操纵稳定性分析

通过本发明，建立了8自由度整车平顺性模型，利用Matlab/Simulink进行仿真实验，分析驾驶员座椅的乘坐舒适性，指出了所研究汽车平顺性的不足。通过田口方法试验，分析了汽车动力学参数变化与外界随机不可控因素对汽车平顺性的影响；优选出了汽车动力学参数的最佳水平组合，提高了汽车的操纵稳定性、平顺性及其稳健性。且通过本发明分析了悬架和座椅的动力学参数对汽车舒适性的影响和贡献率，对汽车的设计与生产标准提供了理论依据。本文的研究方法对汽车的振动、噪声等其它方面的研究具有一定的借鉴意义。

最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (3)

1.一种基于田口实验的汽车平顺性优化分析方法，其特征在于：包括以下步骤：

S1、基于田口方法实验，确定需要进行优化的影响汽车平顺性的可控因子和噪声因子，其中，可控因子包括前后悬架的刚度、阻尼和驾驶员座椅的刚度、阻尼，噪声因子包括路面等级、载荷和车速，同时确定每个可控因子和噪声因子的水平数；

S2、根据可控因子的个数和其水平数生成田口正交表作为内表；

根据噪声因子的个数和其水平数生成田口正交表为外表；

由内、外表的设计得到各可控因子与噪声因子的组合正交表的表头；

S3、在Matlab/Simulink中对各可控因子与噪声因子的组合依次进行仿真实验，得到座椅垂向振动加速度图，并算出座椅垂向加速度均方根值，并计算各可控因子组合下的响应均值和信噪比；

S4、对步骤S3中所得数据进行统计学分析与计算，得到信噪比的方差分析表，并计算出各可控因子对汽车平顺性及其稳健性的贡献率；

步骤S1中还有步骤S0，进行整车平顺性模型及仿真实验，以确定车型是否需要平顺性优化；

在步骤S0中，利用Matlab/Simulink进行仿真实验，模拟该车在某级路面，空载状况下，以特定速度匀速行驶时的整车振动状况，得到驾驶员座椅垂向加速度随时间变化的曲线图；并在Matlab中计算得到座椅垂向加速度均方根值，座椅垂向加速度均方根值越小，代表汽车平顺性越好，否则，代表车型的汽车平顺性需要优化；

所述步骤S0，具体为：

建立八自由度整车动力学模型；该模型中，车身被简化为刚体，八自由度分别为垂直位移(Z_c)、侧倾角(φ)和俯仰角(θ)3个自由度，驾驶员座椅垂直位移(Z_s)1个自由度，4个非簧载质量分别具有1个垂直位移(Z₁、Z₂、Z₃、Z₄)自由度；

取系统的广义坐标为：Z^T＝[Z₁ Z₂ Z₃ Z₄ Z_c Z_s φ θ]；

由拉格朗日方程和达朗贝尔原理，可得8自由度整车平顺性模型振动微分方程为

式中：M为8×8阶的系统质量矩阵；C为8×8阶的系统阻尼矩阵；K为8×8阶的系统刚度矩阵；K_t为8×4阶的轮胎刚度矩阵；Q为4×1阶的路面激励矩阵；Z为8×1阶的位移矩阵。

2.根据权利要求1所述的基于田口实验的汽车平顺性优化分析方法，其特征在于：所述步骤S1中，每个可控因子的水平数为五，并以第三个水平数代表车型初始设计参数值，其余水平对应的值分别为初始值的倍数；

步骤S2中，选取L₂₅(5⁶)正交表为内表。

3.根据权利要求2所述的基于田口实验的汽车平顺性优化分析方法，其特征在于：所述步骤S1中，每个噪声因子的水平数为三；步骤S2中，选取L₉(3³)正交表为外表。

Images