本发明属于汽车动力学性能仿真优化设计领域，特别涉及一种汽车动力学仿真分析中操纵稳定性和平顺性协同优化方法。
背景技术：
在汽车研发设计领域，汽车动力学性能是汽车开发过程中极其重要的组成部分，其内容主要包括操纵稳定性、平顺性等汽车动力学性能的设计及试验。在以往的汽车设计过程中，往往通过样车的反复试制、试验及改进来实现汽车动力学性能的开发，这一过程往往周期长、成本高、效率低，无法满足快速响应市场需求的要求。如何在最短的时间内开发出质优价廉的汽车产品，已成为各大汽车公司一直追求的目标。基于计算机仿真分析技术的虚拟样机技术的出现，使得这一目标正在成为现实。通过虚拟样机技术的运用，可以在样车制造之前，对汽车产品进行数字化设计和仿真分析，提前对汽车的动力学性能进行测量、评价和优化设计，以达到减少开发周期、降低开发成本、提高产品质量的目的。
随着汽车技术的不断发展，人们对汽车各项性能的要求也越来越高，为了满足这些要求，就需要对汽车的动力学性能进行优化设计。操纵稳定性和平顺性是汽车的两个基本动力学性能，这两项性能的研究对于提高汽车行驶安全性和舒适性具有重要意义。在以往汽车动力学性能优化设计过程中，大多只是对操纵稳定性和平顺性中某一单项指标性能进行单独优化，然而它们是相互影响、相互矛盾、密不可分的，因此，有必要提供一种操纵稳定性和平顺性的协同优化方法，以同时对这两项性能指标进行优化。
技术实现要素：
本发明的目的在于克服现有技术的不足，提供一种将正交试验设计和矩阵型权重分析法运用到汽车动力学仿真分析中，通过对少量正交仿真试验方案结果的分析，快速确定试验因素的最优水平组合方案及其主次顺序，原理简单、可操作性强、可靠性高的汽车动力学仿真分析中操纵稳定性和平顺性协同优化方法。
本发明的目的是通过以下技术方案来实现的：汽车动力学仿真分析中操纵稳定性和平顺性协同优化方法，包括以下步骤：
s1、确定仿真试验考察指标，即汽车的操纵稳定性和平顺性；
s2、确定仿真试验因素及其水平值，根据试验因素及水平值的设置选择合适的正交表；
s3、根据正交表确定的试验方案在adams仿真平台下进行整车操纵稳定性和平顺性试验，并记录各项试验结果，利用极差分析法获取正交试验各因素极差分析结果；
s4、利用矩阵型权重分析法对正交试验极差分析结果进行分析，选取试验因素优化水平组合方案。
进一步地，所述步骤s1具体实现方法为：确定的仿真试验考察指标为各种行驶工况下涉及汽车性能较多的操纵稳定性和平顺性，其中操纵稳定性考察其各项试验的综合评分值，平顺性考察其随机不平路面输入试验中车身在纵向、侧向和垂向三个方向振动的总加权加速度均方根值。
进一步地，所述步骤s2具体实现方法为：确定仿真试验因素时，应考虑那些对整车操纵稳定性和平顺性影响较大的汽车设计参数，且这类参数不会因改变整车结构形式而导致一系列工艺成本的增加，同时基于参数在设计中易于实现修改的原则，选取汽车悬架弹簧刚度、减振器阻尼和整车质心高度作为试验因素，每个试验因素在原始模型参数的基础上均设置三个水平值，并根据试验因素及其水平数确定三因素三水平的正交表。
进一步地，所述步骤s3具体实现方法为：选取的正交表为三因素三水平正交试验中l9(34)的前3列，根据正交表确定的试验方案进行整车操纵稳定性和平顺性随机不平路面输入仿真试验，并记录各项试验结果，利用极差分析法获取正交试验各因素极差分析结果。
进一步地，所述步骤s4具体实现方法为：利用矩阵型权重分析法对步骤s3中的正交试验极差分析结果进行分析，根据极差分析结果中各试验因素的权重大小直接获取其优化水平组合方案，并验证参数优化方案的有效性；其中，矩阵型权重分析法的具体步骤如下：
s41、建立仿真试验考察指标层矩阵p：正交试验中确定的试验因素个数及其水平数分别记为m和n；定义正交试验m个因素中的第i个因素在第j个水平上的试验考察指标的算术平均值为kij；定义试验考察指标值为kij，试验考察指标值的取值根据步骤s1中定义方式进行取值：若试验考察指标值kij越大表示性能越好，则令kij＝kij；若试验考察指标值kij越小表示性能越好，则令kij＝1/kij；构造的试验指标矩阵p如下所示：
s42、建立试验因素层矩阵l：定义
s43、建立试验因素水平层矩阵r：用ri表示正交试验中每个试验因素ai的极差，定义
s44、定义影响试验考察指标值的各因素水平的权重矩阵为：
ω＝plr，ωt＝[ω11ω12…ω1nω21ω22…ω2n…ωm1ωm2…ωmn]t；
s45、由权重矩阵得到影响试验考察指标的各试验因素水平的权重值，其值越大表示该因素该水平对指标的影响程度和优化作用越大，通过比较其大小得到影响试验考察指标的试验因素的主次顺序以及各个因素对应的最大权重值，并将最大权重值作为该因素水平正交试验的最优水平组合方案。
本发明的有益效果是：本发明将正交试验设计和矩阵型权重分析法运用到汽车动力学仿真分析中，通过对少量正交仿真试验方案结果的分析，快速确定试验因素的最优水平组合方案及其主次顺序，以实现汽车操纵稳定性和平顺性的协同优化。该优化方法具有原理简单、可操作性强、设计灵活、可靠性高及适用面广等特点，对汽车动力学性能优化和动力学控制等研究具有重要指导和借鉴意义。
具体实施方式
下面结合具体实施例进一步说明本发明的技术方案。
本发明的一种汽车动力学仿真分析中操纵稳定性和平顺性协同优化方法，包括以下步骤：
s1、确定仿真试验考察指标，即汽车的操纵稳定性和平顺性；具体实现方法为：确定的仿真试验考察指标为各种行驶工况下涉及汽车性能较多的操纵稳定性和平顺性，其中操纵稳定性考察其各项试验的综合评分值，平顺性考察其随机不平路面输入试验中车身在纵向、侧向和垂向三个方向振动的总加权加速度均方根值。
s2、确定仿真试验因素及其水平值，根据试验因素及水平值的设置选择合适的正交表；具体实现方法为：确定仿真试验因素时，应考虑那些对整车操纵稳定性和平顺性影响较大的汽车设计参数，且这类参数不会因改变整车结构形式而导致一系列工艺成本的增加，同时基于参数在设计中易于实现修改的原则，选取汽车悬架弹簧刚度、减振器阻尼和整车质心高度作为试验因素，每个试验因素在原始模型参数的基础上均设置三个水平值，并根据试验因素及其水平数确定三因素三水平的正交表。
本实施例采用的实验因素及其水平值如表一所示。
表一
减振器的阻尼特性为非线性，其特性参数用d0表示。
s3、根据正交表确定的试验方案在adams仿真平台下进行整车操纵稳定性和平顺性试验，并记录各项试验结果，利用极差分析法获取正交试验各因素极差分析结果；具体实现方法为：选取的正交表为三因素三水平正交试验中l9(34)的前3列，根据正交表确定的试验方案进行整车操纵稳定性和平顺性随机不平路面输入仿真试验，并记录各项试验结果。
根据表二设定的正交试验方案进行整车操纵稳定性和平顺性随机不平路面输入仿真试验，记录得到的试验结果如表三所示。
表二
表三
利用极差分析法得到各试验因素及其水平对这两项指标性能的影响情况如表四所示：
表四
由表四分析可得：影响操纵稳定性的试验因素主次顺序为：刚度>质心高度>阻尼，若单独考虑操纵稳定性的优化问题，根据考察指标值越大越好的原则，得到的参数水平优化组合方案为：刚度1、阻尼3、质心高度1；同理，可得到影响平顺性随机输入的试验因素主次顺序为：质心高度>刚度>阻尼，若单独对平顺性随机输入进行优化，根据考察指标值越小越好的原则，得到的参数水平优化组合方案为：刚度1、阻尼3、质心高度3。纵上所述，两个考察指标的参数水平优化组合方案并不一致，为了确定最终的参数水平优化组合方案，需要进一步采用本发明中使用的矩阵型权重分析法。
s4、利用矩阵型权重分析法对正交试验极差分析结果进行分析，选取试验因素优化水平组合方案；具体实现方法为：利用矩阵型权重分析法对步骤s3中的正交试验极差分析结果进行分析，根据极差分析结果各试验因素的权重大小直接获取其优化水平组合方案，并验证参数优化方案的有效性；其中，矩阵型权重分析法的具体步骤如下：
s41、建立仿真试验考察指标层矩阵p：正交试验中确定的试验因素个数及其水平数分别记为m和n；定义正交试验m个因素中的第i个因素在第j个水平上的试验考察指标的算术平均值为kij(第“i个因素在第j个水平上”的意思是第i个因素取第j个水平时所对应的参数值，比如本实施例中的“第一个因素在第二个水平上”即表示刚度取值为125n/mm)；定义试验考察指标值为kij，试验考察指标值的取值根据步骤s1中定义方式进行取值：若试验考察指标值kij越大表示性能越好，则令kij＝kij；若试验考察指标值kij越小表示性能越好，则令kij＝1/kij；构造的试验指标矩阵p如下所示：
s42、建立试验因素层矩阵l：定义
s43、建立试验因素水平层矩阵r：用ri表示正交试验中每个试验因素ai的极差，定义
s44、定义影响试验考察指标值的各因素水平的权重矩阵为：
ω＝plr，ωt＝[ω11ω12…ω1nω21ω22…ω2n…ωm1ωm2…ωmn]t；
s45、由权重矩阵得到影响试验考察指标的各试验因素水平的权重值，其值越大表示该因素该水平对指标的影响程度和优化作用越大，通过比较其大小得到影响试验考察指标的试验因素的主次顺序以及各个因素对应的最大权重值，并将最大权重值作为该因素水平正交试验的最优水平组合方案。
为方便起见，本实施例中分别用a和b表示试验中两个考察指标：操纵稳定性和平顺性随机输入，d、e、f分别表示试验因素刚度、阻尼、质心高度；其中a指标值越大越好，构造试验因素层矩阵时取kij＝kij；而b指标值追求越小越好，构造其试验因素层矩阵时取kij＝1/kij。则构造的两个指标的六个矩阵因子如下：
依据权重矩阵计算公式，得到两个试验考察指标的因素水平权重矩阵ωa和ωb，该正交试验的总权重矩阵为两个试验考察指标因素水平独立权重矩阵的算术平均值：
(ωa+ωb)/2＝ω总
得到的三个权重矩阵分别为：
由总权重矩阵可以得到影响试验考察指标的各试验因素水平的权重值，其值越大表示该因素该水平对指标的影响程度和优化作用越大。通过比较其大小可得到影响试验考察指标的三个试验因素的主次顺序为：质心高度>刚度>阻尼，且对应的最大权重值的水平为：f3、d1、e3。因此，最后确定的该三因素三水平正交试验的最优水平组合方案为：d1、e3、f3，即刚度为100n/mm，阻尼为1.2d0，质心高度为306mm。该参数优化水平组合方案即为正交试验设计中的方案3，该优化方案的两个考察指标分别为：操纵稳定性得分84.642，平顺性随机输入总加权加速度均方根值为0.129m/s2，相对原模型的83.505和0.140m/s2分别提高了1.36％和降低了7.8％。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。