本发明涉及领域，具体是涉及一种基于python的纯电动汽车扭矩分配系数分析方法及系统。
背景技术：
：目前新能源汽车市场中，纯电动汽车驱动方式多为两驱，四驱的车辆占比较少，和两驱相比四驱可以最大限度的利用车辆全部的附着力，可以显著提升车辆动力性。但由于四驱车辆和两驱车辆相比增加了一套动力系统，所以四驱车辆要比两驱车辆明显重量大，大重量增加车辆附着力的同时，也增加了车辆阻力，相应的会增加车辆能耗。车辆能耗主要由克服车辆阻力损耗和车辆动力系统损耗组成，降低车辆动力系统损耗可以降低车辆能耗。可以通过保证动力系统效率最优来降低车辆动力系统损耗。目前纯电动汽车四驱前后电机扭矩的分配按照固定的比例进行分配，并没有结合实时的工况，因此无法达到系统最优。技术实现要素：本发明的目的是为了克服上述
背景技术：
的不足，提供一种基于python的纯电动汽车扭矩分配系数分析方法及系统，对前后电机进行扭矩分配，可以保证在满足车辆动力性需求时车辆动力系统效率时刻保持最优，降低车辆能耗，提升车辆续驶里程第一方面，提供一种基于python的纯电动汽车扭矩分配系数分析方法，包括以下步骤：通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点；获取采样点前电机效率和采样点后电机效率，采样点为所述车辆工况点中任意多个；通过python第三方库numpy模块、math模块、scipy模块及panda模块，根据所述采样点前电机效率和所述采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率；根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率。根据第一方面，在第一方面的第一种可能的实现方式中，所述“通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点”步骤，具体包括以下步骤：读取前电机外特性数据和后电机外特性数据；根据所述前电机外特性数据和所述后电机外特性数据分别确定前电机最大转速、后电机最大转速、任一转速下的前电机最大扭矩及任一转速下的后电机最大扭矩，选取所述前电机最大转速和所述后电机最大转速中的较小值作为动力系统最大转速，选取任一转速下的所述前电机最大扭矩和所述后电机最大扭矩作为所述任一转速下的动力系统最大扭矩；按照预设采样规则，结合所述动力系统最大转速和所述动力系统最大扭矩生成各车辆工况点，所述车辆工况点包括电机转速和及对应的动力系统扭矩。根据第一方面，在第一方面的第二种可能的实现方式中，所述“根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率”步骤之前，具体还包括以下步骤：选取任一目标车辆工况点，获取所述目标车辆工况点的目标动力系统扭矩；若所述目标动力系统扭矩小于等于第一扭矩值t1，则所述目标车辆工况点的扭矩分配系数范围设为[0,1]，所述第一扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较小值；若所述目标动力系统扭矩大于第一扭矩值t1、且小于等于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为所述第二扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较大值；若所述目标动力系统扭矩大于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为根据第一方面的第二种可能的实现方式，在第一方面的第三种可能的实现方式中，所述“根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率”步骤，具体包括以下步骤：在所述扭矩分配系数范围内设置多组效率比对组，所述效率比对组中将所述目标动力系统扭矩按照不同的扭矩分配系数细分为相应的目标前电机扭矩和目标后电机扭矩；根据所述目标车辆工况点的目标转速和所述目标前电机扭矩确定对应的目标前电机效率，根据所述目标车辆工况点的目标转速和所述目标后电机扭矩确定对应的目标后电机效率；获取前减速器效率和后减速器效率；根据所述目标前电机扭矩、目标后电机扭矩、所述目标前电机效率、所述目标后电机效率、所述前减速器效率和所述后减速器效率分析各效率比对组的动力系统效率；根据所述动力系统效率确定所述最优扭矩分配系数。根据第一方面的第三种可能的实现方式，在第一方面的第四种可能的实现方式中，，所述“根据所述目标前电机扭矩、目标后电机扭矩、所述目标前电机效率、所述目标后电机效率、所述前减速器效率和所述后减速器效率分析各效率比对组的动力系统效率”步骤，具体包括以下步骤：根据所述目标前电机扭矩得到前电机输出功率p01，根据所述目标后电机扭矩得到后电机输出功率p02；根据所述目标前电机效率η1和所述前电机输出功率p01计算前电机输入功率p11，根据所述目标后电机效率η2和所述后电机输出功率p02计算后电机输入功率p12，根据所述前电机输出功率p01、所述后电机输出功率p02、所述前电机输入功率p11、所述后电机输入功率p12、所述前减速器效率η3和所述后减速器效率η4分析各效率比对组的动力系统效率η，第二方面，提供一种基于python的纯电动汽车扭矩分配系数分析系统，包括：数据读取处理模块，用于通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点，获取采样点前电机效率和采样点后电机效率，采样点为所述车辆工况点中任意多个；工况点效率分析模块，与所述数据读取处理模块通信连接，用于通过python第三方库numpy模块、math模块、scipy模块及panda模块，根据所述采样点前电机效率和所述采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率；以及，最优分配系数计算模块，与所述工况点效率分析模块通信连接，用于根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率。根据第二方面，在第二方面的第一种可能的实现方式中，所述数据读取处理模块包括：数据读取单元，用于读取前电机外特性数据和后电机外特性数据；数据处理单元，与所述数据读取单元通信连接，用于根据所述前电机外特性数据和所述后电机外特性数据分别确定前电机最大转速、后电机最大转速、任一转速下的前电机最大扭矩及任一转速下的后电机最大扭矩，选取所述前电机最大转速和所述后电机最大转速中的较小值作为动力系统最大转速，选取任一转速下的所述前电机最大扭矩和所述后电机最大扭矩作为所述任一转速下的动力系统最大扭矩；按照预设采样规则，结合所述动力系统最大转速和所述动力系统最大扭矩生成各车辆工况点，所述车辆工况点包括电机转速和及对应的动力系统扭矩；数据获取单元，与所述数据处理单元通信连接，用于用于获取采样点前电机效率和采样点后电机效率，采样点为所述车辆工况点中任意多个。根据第二方面，在第二方面的第二种可能的实现方式中，所述系统还包括：扭矩分配约束模块，与所述工况点效率分析模块通信连接，用于：选取任一目标车辆工况点，获取所述目标车辆工况点的目标动力系统扭矩；若所述目标动力系统扭矩小于等于第一扭矩值t1，则所述目标车辆工况点的扭矩分配系数范围设为[0,1]，所述第一扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较小值；若所述目标动力系统扭矩大于第一扭矩值t1、且小于等于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为所述第二扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较大值；若所述目标动力系统扭矩大于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为根据第二方面的第二种可能的实现方式，在第二方面的第三种可能的实现方式中，所述最优分配系数计算模块，与所述扭矩分配约束模块通信连接，包括：分组单元，用于在所述扭矩分配系数范围内设置多组效率比对组，所述效率比对组中将所述目标动力系统扭矩按照不同的扭矩分配系数细分为相应的目标前电机扭矩和目标后电机扭矩；电机效率分析单元，与所述分组单元通信连接，用于根据所述目标车辆工况点的目标转速和所述目标前电机扭矩确定对应的目标前电机效率，根据所述目标车辆工况点的目标转速和所述目标后电机扭矩确定对应的目标后电机效率；效率获取单元，用于获取前减速器效率和后减速器效率；系统效率分析单元，与所述电机效率分析单元和所述效率获取单元通信连接，用于根据所述目标前电机扭矩、目标后电机扭矩、所述目标前电机效率、所述目标后电机效率、所述前减速器效率和所述后减速器效率分析各效率比对组的动力系统效率；最优系数分析单元，与所述系统效率分析单元通信连接，用于根据所述动力系统效率确定所述最优扭矩分配系数。根据第二方面的第三种可能的实现方式，在第二方面的第四种可能的实现方式中，所述系统效率分析单元包括：输入功率分析子单元，用于根据所述目标前电机扭矩得到前电机输出功率p01，根据所述目标后电机扭矩得到后电机输出功率p02；输出功率分析子单元，与所述输入功率分析子单元通信连接，用于根据所述目标前电机效率η1和所述前电机输出功率p01计算前电机输入功率p11，根据所述目标后电机效率η2和所述后电机输出功率p02计算后电机输入功率p12，系统效率分析子单元，与所述输入功率分析子单元和所述输出功率分析子单元通信连接，用于根据所述前电机输出功率p01、所述后电机输出功率p02、所述前电机输入功率p11、所述后电机输入功率p12、所述前减速器效率η3和所述后减速器效率η4分析各效率比对组的动力系统效率η，与现有技术相比，本发明使用python语言用于前后电机扭矩分配系数计算，找出所有工况点下的车辆动力系统效率最优的分配系数，以便将所有最优分配系数组成分配系数表，四驱模式下车辆在进行扭矩控制过程中可以实时查询分配系数表，对前后电机进行扭矩分配，可以保证在满足车辆动力性需求时车辆动力系统效率时刻保持最优，降低车辆能耗，提升车辆续驶里程。附图说明图1是本发明一种基于python的纯电动汽车扭矩分配系数分析方法一实施例的流程示意图；图2是本发明一种基于python的纯电动汽车扭矩分配系数分析方法另一实施例的流程示意图；图3是本发明一种基于python的纯电动汽车扭矩分配系数分析方法另一实施例的流程示意图；图4是本发明一种基于python的纯电动汽车扭矩分配系数分析系统一实施例的结构示意图。附图标号：100、基于python的纯电动汽车扭矩分配系数分析系统；110、数据读取处理模块；111、数据读取单元；112、数据处理单元；113、数据获取单元；120、工况点效率分析模块；130、最优分配系数计算模块；131、分组单元；132、电机效率分析单元；133、效率获取单元；134、系统效率分析单元；135、最优系数分析单元；140、扭矩分配约束模块。具体实施方式现在将详细参照本发明的具体实施例，在附图中例示了本发明的例子。尽管将结合具体实施例描述本发明，但将理解，不是想要将本发明限于所述的实施例。相反，想要覆盖由所附权利要求限定的在本发明的精神和范围内包括的变更、修改和等价物。应注意，这里描述的方法步骤都可以由任何功能块或功能布置来实现，且任何功能块或功能布置可被实现为物理实体或逻辑实体、或者两者的组合。为了使本领域技术人员更好地理解本发明，下面结合附图和具体实施方式对本发明作进一步详细说明。注意：接下来要介绍的示例仅是一个具体的例子，而不作为限制本发明的实施例必须为如下具体的步骤、数值、条件、数据、顺序等等。本领域技术人员可以通过阅读本说明书来运用本发明的构思来构造本说明书中未提到的更多实施例。参见图1所示，本发明实施例提供一种基于python的纯电动汽车扭矩分配系数分析方法，包括以下步骤：通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点；获取采样点前电机效率和采样点后电机效率，采样点为车辆工况点中任意多个；通过python第三方库numpy模块、math模块、scipy模块及panda模块，根据采样点前电机效率和采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率；根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，前电机效率包括采样点前电机效率和非采样点前电机效率，后电机效率包括采样点后电机效率和非采样点后电机效率。具体的，本实施例中，通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点，各个车辆工况点中转速和扭矩的取值阶梯设置，其阶梯值基于电机的参数、运转特性进行设置，或根据需要自主设置。选取车辆工况点中若干工况点作为采样点，获取采样点前电机效率和采样点后电机效率，采样点前电机效率和采样点后电机效率为该采样点工况下实测的前电机和后电机的效率。其中，采样点为车辆工况点中任意多个，且采样点中分配给前电机和后电机的扭矩的值也是唯一确定的，即采样点中前电机转速和扭矩值是确定的值，后电机转速和扭矩值也是确定的值，该条件下前电机和后电机的效率也是唯一确定的值。通过python第三方库numpy模块、math模块、scipy模块及panda模块，根据采样点前电机效率和采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率，插值方法采用”最近取值法”。采样点的转速、扭矩以及效率均已知，结合非采样点的转速和扭矩得到相应的效率。根据不同的转速以及扭矩下的所有的前电机效率和后电机效率分析各车辆工况点不同分配系数下的动力系统的效率，从而得到每个车辆工况点的最优扭矩分配系数，前电机效率包括采样点前电机效率和非采样点前电机效率，后电机效率包括采样点后电机效率和非采样点后电机效率。将每个车辆工况点的最优扭矩分配系数生成最优扭矩分配系数表格，后续在车辆行驶过程中，根据车辆的电机转速以及需求扭矩找到相应的匹配的车辆工况点，选取对应的最优扭矩分配系数对前电机和后电机进行扭矩分配，以获取当前状况最优的扭矩分配结果。本发明使用python语言编写脚本文件进行车辆动力系统效率最优前后电机扭矩分配系数计算。python语言是一个高层次的结合了解释性、编译性、互动性和面向对象的脚本语言，具有很强的可读性并且易于学习、易于阅读，另外python具有丰富的第三方库，可以协助处理各种工作，其中在科学计算和大数据处理方面使用python尤为简单方便。使用python语言用于车辆动力系统效率最优前后电机扭矩分配系数计算，可以达到简单、方便、快捷的效果。可选地，在本申请另外的实施例中，“通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点”步骤，具体包括以下步骤：读取前电机外特性数据和后电机外特性数据；根据前电机外特性数据和后电机外特性数据分别确定前电机最大转速、后电机最大转速、任一转速下的前电机最大扭矩及任一转速下的后电机最大扭矩，选取前电机最大转速和后电机最大转速中的较小值作为动力系统最大转速，选取任一转速下的前电机最大扭矩和后电机最大扭矩作为任一转速下的动力系统最大扭矩；按照预设采样规则，结合动力系统最大转速和动力系统最大扭矩生成各车辆工况点，车辆工况点包括电机转速和及对应的动力系统扭矩。具体的，本实施例中，读取前电机外特性数据和后电机外特性数据，依据前后电机实测外特性数据，明确前电机转速范围以及特定转速对应的扭矩范围，明确后电机转速范围以及特定转速对应的扭矩范围。根据上述明确的前后电机转速范围和特定转速对应的扭矩范围，计算车辆动力系统转速范围和特定转速对应扭矩范围，车辆动力系统包括前电机加减速器系统和后电机加减速器系统。根据上述明确的车辆动力系统转速范围以及特定转速对应的扭矩范围，进行车辆工况点(转速-扭矩)的设计。例如，划分原则(仅为了便于理解进行举例说明，并不限定为下列划分原则)如下：a)转速以500rpm为起始转速，以最高转速为终止转速，以500rpm为梯度得到全转速范围内的转速工况点。b)选取第a步某一转速工况点，扭矩以10nm为起始扭矩，以当前转速最大扭矩为终止扭矩，以10nm为梯度得到当前转速下所有扭矩工况点。c)重复第b步，得到所有转速下的扭矩工况点，进而得出所有扭矩分配系数计算工况点。本申请导入前后电机外特性数据，确定动力系统转速最大值和车辆动力系统转速范围；根据车辆动力系统转速范围，以及前后电机特定转速下最大扭矩设计车辆工况点更加合理。可选地，如图2所示，在本申请另外的实施例中，“根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，前电机效率包括采样点前电机效率和非采样点前电机效率，后电机效率包括采样点后电机效率和非采样点后电机效率”步骤之前，具体还包括以下步骤：选取任一目标车辆工况点，获取目标车辆工况点的目标动力系统扭矩；若目标动力系统扭矩小于等于第一扭矩值t1，则目标车辆工况点的扭矩分配系数范围设为[0,1]，第一扭矩值为目标转速下前电机最大扭矩和后电机最大扭矩的较小值；若目标动力系统扭矩大于第一扭矩值t1、且小于等于第二扭矩值t2，则目标车辆工况点的扭矩分配系数范围设为第二扭矩值为目标转速下前电机最大扭矩和后电机最大扭矩的较大值；若目标动力系统扭矩大于第二扭矩值t2，则目标车辆工况点的扭矩分配系数范围设为具体的，本实施例中，选取任一目标车辆工况点，获取目标车辆工况点的目标动力系统扭矩，即选取目标车辆工况点下的目标转速下和目标动力系统扭矩，根据与该目标转速下前电机最大扭矩和后电机最大扭矩进行比对，然后得到该目标车辆工况点的扭矩分配系数范围。因此本申请在进行分配系数范围约束时，方法较为简单且所选范围能覆盖所有满足要求的取值。可选地，在本申请另外的实施例中，“根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，前电机效率包括采样点前电机效率和非采样点前电机效率，后电机效率包括采样点后电机效率和非采样点后电机效率”步骤，具体包括以下步骤：在扭矩分配系数范围内设置多组效率比对组，效率比对组中将目标动力系统扭矩按照不同的扭矩分配系数细分为相应的目标前电机扭矩和目标后电机扭矩；根据目标车辆工况点的目标转速和目标前电机扭矩确定对应的目标前电机效率，根据目标车辆工况点的目标转速和目标后电机扭矩确定对应的目标后电机效率；获取前减速器效率和后减速器效率；根据目标前电机扭矩、目标后电机扭矩、目标前电机效率、目标后电机效率、前减速器效率和后减速器效率分析各效率比对组的动力系统效率；根据动力系统效率确定最优扭矩分配系数。具体的，本实施例中，各车辆工况点中只设有转速和对应的动力系统扭矩，而动力系统扭矩并没有具体分配到前电机和后电机。因此对于每一个车辆工况点，在扭矩分配系数范围内设置多组效率比对组，每组效率比对组内按照一定的扭矩分配系数将动力系统扭矩分别分配给前电机和后电机，也就是说在每组效率比对组内前电机和后电机的转速以及扭矩都是确定的，因此由上述实施例中的说明中采样点前电机效率、采样点后电机效率、非采样点前电机效率及非采样点后电机效率，可以确定每组效率比对组内目标前电机效率和目标后电机效率。其中，效率比对组之间扭矩分配系数的变化阶梯可以根据不同的需求进行设置，例如，扭矩分配系数变化梯度设定为0.01。由于车辆动力系统包括前电机加减速器系统和后电机加减速器系统，也就是减速器效率会对车辆动力系统产生影响，因此获取前减速器效率和后减速器效率。根据目标前电机扭矩、目标后电机扭矩、目标前电机效率、目标后电机效率、前减速器效率和后减速器效率分析各效率比对组的动力系统效率。将各效率比对组的动力系统效率进行比对，选取其中效率最高的效率比对组对应的扭矩分配系数作为最优扭矩分配系数。本申请设置多组扭矩分配系数接阶梯变化的效率比对组分别计算其动力系统效率，以便更加全面地选取出最优扭矩分配系数。可选地，在本申请另外的实施例中，“根据目标前电机扭矩、目标后电机扭矩、目标前电机效率、目标后电机效率、前减速器效率和后减速器效率分析各效率比对组的动力系统效率”步骤，具体包括以下步骤：根据所述目标前电机扭矩得到前电机输出功率p01，根据所述目标后电机扭矩得到后电机输出功率p02；根据所述目标前电机效率η1和所述前电机输出功率p01计算前电机输入功率p11，根据所述目标后电机效率η2和所述后电机输出功率p02计算后电机输入功率p12，根据所述前电机输出功率p01、所述后电机输出功率p02、所述前电机输入功率p11、所述后电机输入功率p12、所述前减速器效率η3和所述后减速器效率η4分析各效率比对组的动力系统效率η，具体的，本实施例中，根据扭矩值可以直接计算得到相应的功率值，目标前电机扭矩和目标后电机扭矩都是需要电机输出的扭矩值，因此根据目标前电机扭矩得到前电机输出功率p01，根据目标后电机扭矩得到后电机输出功率p02，即计算得到该车辆工况在该效率对比组扭矩分配系数下前电机输出功率和后电机输出功率，由于转速和扭矩一定时，前电机效率和后电机效率一定且已经获取，因此根据目标前电机效率η1和前电机输出功率p01计算前电机输入功率p11，根据目标后电机效率η2和后电机输出功率p02计算后电机输入功率p12，最后根据前电机输出功率p01、后电机输出功率p02、前电机输入功率p11、后电机输入功率p12、前减速器效率η3和后减速器效率η4分析各效率比对组的动力系统效率η，本发明依托上述实例编写计算车辆动力系统效率最优前后电机扭矩分配系数脚本程序如下：参见图1所示，本发明实施例提供一种基于python的纯电动汽车扭矩分配系数分析方法，包括以下步骤：s1、导入前后电机外特性数据，其中前电机外特性数据如表一所示，后电机外特性数据如表二所示；表一前电机外特性数据转速(r/min)扭矩(nm)50023210002331500235……1250091表二后电机外特性数据转速(r/min)扭矩(nm)50030810003051500305……1600088根据前后电机外特性数据可知，nmax_f＝12500,nmax_r＝16000,设定前后减速器速比一致，所以取动力系统转速最大值nmax＝min(nmax_f,nmax_r)＝12500,车辆动力系统转速范围为0-12500。s2、根据车辆动力系统转速范围，以及前后电机特定转速下最大扭矩tmax_n＝tmax_f_n tmax_r_n，其中n＝500,1000,1500…12500,设计车辆工况点如表三所示：表三车辆工况点转速(r/min)扭矩(nm)500105002050030…500540100010……12500198s3、计算使每个工况点车辆动力系统效率最优的前后电机扭矩分配系数，根据本发明所阐述的扭矩分配系数计算方法，首先要比较该工况点下扭矩值与该工况点对应的转速下前后电机最大扭矩的大小关系。如果t_n<t_max_f_nandt_n<t_max_r_n，则限定前后电机扭矩分配系数范围为[0,1]，通过如下语句实现该条件判断：iffloat(ws.cell(row＝i,column＝2).value)/torque_max_90<1:proportions＝list(range(0,101,1))设定本实施例中同一转速下，后电机最大扭矩始终大于前电机最大扭矩，所以只判断了t_n与t_max_f_n的大小关系。如果t_n>t_max_f_nandt_n<t_max_r_n，则限定前后电机扭矩分配系数范围为[0,t_max_f_n/tmax_n]，通过如下语句实现该条件判断：eliffloat(ws.cell(row＝i,column＝2).value)/torque_max_90>1andfloat(ws.cell(row＝i,column＝2).value)/torque_max_160<1:a＝round(float(torque_max_160/ws.cell(row＝i,column＝2).value))*100b＝round(float(torque_max_90/ws.cell(row＝i,column＝2).value))*100proportions_2＝list(range(0,b 1,1))如果t_n<t_max_f_nandt_n<t_max_r_n，则限定前后电机扭矩分配系数范围为[t_max_f_n/tmax_n,t_max_r_n/tmax_n]，通过如下语句实现该条件判断：eliffloat(ws.cell(row＝i,column＝2).value)/torque_max_160>1:a＝math.floor((float(torque_max_160/ws.cell(row＝i,column＝2).value)*100))b＝math.ceil(float(torque_max_90/ws.cell(row＝i,column＝2).value)*100)proportions_3＝list(range(b,a,1))限定前后电机扭矩分配系数范围的目的是，保证前后电机实际所分配的扭矩不超过当前转速下前后电机的最大扭矩能力。s4、计算每个车辆工况点下车辆动力系统效率，根据下式计算每个车辆工况点前后电机输出功率：p输出_前电机_n＝t_n*n/9550,其中n＝500,1000,1500…12500，p输出_后电机_n＝t_n*n/9550,其中n＝500,1000,1500…12500；根据下式计算前后电机输入功率：,p输入_前电机_n＝p输出_前电机_n/η_f_n,其中n＝500,1000,1500…12500，p输入_后电机_n＝p输出_后电机_n/η_r_n,其中n＝500,1000,1500…12500；根据下式计算车辆动力系统效率：η_n＝(p输出_前电机_n p输出_后电机_n)/(p输入_前电机_n p输入_后电机_n)*0.97，其中0.97为前后减速器效率，仅为了便于理解进行举例说明，实际减速器效率并不进行限定。上述公式中每个工况点下的前电机和后电机的效率，需要根据前后电机实测效率数据，通过插值的方法得出，本发明中编写的脚本程序中，如下语句为前后电机效率插值程序：result_ef_90＝griddata(points_90,values_90,points_new_90,method＝'nearest')result_ef_160＝griddata(points_160,values_160,points_new_160,method＝'nearest')插值方法采用“最近取值法”，其中插值程序中的前后电机样本点通过如下语句得出：ws_90＝wb['前电机外特性']ws_160＝wb['后电机外特性']rows＝ws.max_rowrows_90＝ws_90.max_rowrows_160＝ws_160.max_rowlist_speed_90＝[]list_torque_90＝[]list_ef_90＝[]foriinrange(2,rows_90 1):speed_90＝ws_90.cell(row＝i,column＝1).valuetorque_90＝ws_90.cell(row＝i,column＝2).valueeff_90＝ws_90.cell(row＝i,column＝3).valuelist_speed_90.append(speed_90)list_torque_90.append(torque_90)list_ef_90.append(eff_90)list_new_90＝list(zip(list_speed_90,list_torque_90))points_90＝np.array(list_new_90)values_90＝np.array(list_ef_90)list_speed_160＝[]list_torque_160＝[]list_ef_160＝[]foriinrange(2,rows_160 1):speed_160＝ws_160.cell(row＝i,column＝1).valuetorque_160＝ws_160.cell(row＝i,column＝2).valueeff_160＝ws_160.cell(row＝i,column＝3).valuelist_speed_160.append(speed_160)list_torque_160.append(torque_160)list_ef_160.append(eff_160)list_new_160＝list(zip(list_speed_160,list_torque_160))points_160＝np.array(list_new_160)values_160＝np.array(list_ef_160)得到前电机样本点为(points_90,values_90),后电机样本点为(points_160,values_160)。前电机插值点和后电机插值点，通过如下语句得出：points_new_90＝np.array([speed_90,torque_90])points_new_160＝np.array([speed_160,torque_160])在得出样本点和插值点之后，可以得出前后电机插值点的效率result_ef_90和result_ef_160，通过如下语句可以计算每个工况点的车辆动力系统效率：front_power_input＝((float(torque_90)*float(speed_90))/9550/float(result_ef_90))*100front_power_output＝(float(torque_90)*float(speed_90))/9550*0.97rear_power_input＝((float(torque_160)*float(speed_160))/9550/float(result_ef_160))*100rear_power_output＝(float(torque_160)*float(speed_160))/9550*0.97total_power_input＝front_power_input rear_power_inputtotal_power_output＝front_power_output rear_power_outputeff_n＝total_power_output/total_power_inputs5、限定前后电机扭矩分配系数范围之后，则需要在限定的系数范围之内，找出使车辆动力系统效率最优的点，并将所有点组成车辆动力系统效率最优前后电机扭矩分配系数表，通过如下语句可得出该分配表：eff_max.append(proportions[eff.index(max(eff))]/100)eff_max_coefficient.append(max(eff))通过上述步骤，本发明所依托的实例最后得出的每个车辆工况点的扭矩分配系数如表四所示；表四车辆工况点的扭矩分配系数转速(r/min)扭矩(nm)分配系数500101500200500300………5005400.43…125001980.45参见图4所示，本发明实施例提供一种基于python的纯电动汽车扭矩分配系数分析系统100，包括：数据读取处理模块110，用于通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点，获取采样点前电机效率和采样点后电机效率，采样点为车辆工况点中任意多个；工况点效率分析模块120，与数据读取处理模块110通信连接，用于通过python第三方库numpy模块、math模块、scipy模块及panda模块，根据采样点前电机效率和采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率；以及，最优分配系数计算模块130，与工况点效率分析模块120通信连接，用于根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，前电机效率包括采样点前电机效率和非采样点前电机效率，后电机效率包括采样点后电机效率和非采样点后电机效率。数据读取处理模块110包括：数据读取单元111，用于读取前电机外特性数据和后电机外特性数据；数据处理单元112，与数据读取单元111通信连接，用于根据前电机外特性数据和后电机外特性数据分别确定前电机最大转速、后电机最大转速、任一转速下的前电机最大扭矩及任一转速下的后电机最大扭矩，选取前电机最大转速和后电机最大转速中的较小值作为动力系统最大转速，选取任一转速下的前电机最大扭矩和后电机最大扭矩作为任一转速下的动力系统最大扭矩；按照预设采样规则，结合动力系统最大转速和动力系统最大扭矩生成各车辆工况点，车辆工况点包括电机转速和及对应的动力系统扭矩；数据获取单元113，与数据处理单元112通信连接，用于用于获取采样点前电机效率和采样点后电机效率，采样点为车辆工况点中任意多个。系统还包括：扭矩分配约束模块140，与工况点效率分析模块120通信连接，用于：选取任一目标车辆工况点，获取目标车辆工况点的目标动力系统扭矩；若目标动力系统扭矩小于等于第一扭矩值t1，则目标车辆工况点的扭矩分配系数范围设为[0,1]，第一扭矩值t1为目标转速下前电机最大扭矩和后电机最大扭矩的较小值；若目标动力系统扭矩大于第一扭矩值t1、且小于等于第二扭矩值t2，则目标车辆工况点的扭矩分配系数范围设为第二扭矩值t2为目标转速下前电机最大扭矩和后电机最大扭矩的较大值；若目标动力系统扭矩大于第二扭矩值t2，则目标车辆工况点的扭矩分配系数范围设为最优分配系数计算模块130，与扭矩分配约束模块140通信连接，包括：分组单元131，用于在扭矩分配系数范围内设置多组效率比对组，效率比对组中将目标动力系统扭矩按照不同的扭矩分配系数细分为相应的目标前电机扭矩和目标后电机扭矩；电机效率分析单元132，与分组单元131通信连接，用于根据目标车辆工况点的目标转速和目标前电机扭矩确定对应的目标前电机效率，根据目标车辆工况点的目标转速和目标后电机扭矩确定对应的目标后电机效率；效率获取单元1333，用于获取前减速器效率和后减速器效率；系统效率分析单元134，与电机效率分析单元132和效率获取单元133通信连接，用于根据目标前电机扭矩、目标后电机扭矩、目标前电机效率、目标后电机效率、前减速器效率和后减速器效率分析各效率比对组的动力系统效率；最优系数分析单元135，与系统效率分析单元134通信连接，用于根据动力系统效率确定最优扭矩分配系数。系统效率分析单元134包括：输入功率分析子单元，用于根据目标前电机扭矩得到前电机输出功率p01，根据目标后电机扭矩得到后电机输出功率p02；输出功率分析子单元，与输入功率分析子单元通信连接，用于根据目标前电机效率η1和前电机输出功率p01计算前电机输入功率p11，根据目标后电机效率η2和后电机输出功率p02计算后电机输入功率p12，系统效率分析子单元，与输入功率分析子单元和输出功率分析子单元通信连接，用于根据前电机输出功率p01、后电机输出功率p02、前电机输入功率p11、后电机输入功率p12、前减速器效率η3和后减速器效率η4分析各效率比对组的动力系统效率η，具体的，本实施例中各个模块的功能在相应的方法实施例中已经进行详细说明，因此不再一一赘述。显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。当前第1页1 2 3
技术特征：1.一种基于python的纯电动汽车扭矩分配系数分析方法，其特征在于，包括以下步骤：
通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点；
获取采样点前电机效率和采样点后电机效率，采样点为所述车辆工况点中任意多个；
通过python第三方库numpy模块、math模块、scipy模块及panda模块，根据所述采样点前电机效率和所述采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率；
根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率。
2.如权利要求1所述的方法，其特征在于，所述“通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点”步骤，具体包括以下步骤：
读取前电机外特性数据和后电机外特性数据；
根据所述前电机外特性数据和所述后电机外特性数据分别确定前电机最大转速、后电机最大转速、任一转速下的前电机最大扭矩及任一转速下的后电机最大扭矩，选取所述前电机最大转速和所述后电机最大转速中的较小值作为动力系统最大转速，选取任一转速下的所述前电机最大扭矩和所述后电机最大扭矩作为所述任一转速下的动力系统最大扭矩；按照预设采样规则，结合所述动力系统最大转速和所述动力系统最大扭矩生成各车辆工况点，所述车辆工况点包括电机转速和及对应的动力系统扭矩。
3.如权利要求1所述的方法，其特征在于，所述“根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率”步骤之前，具体还包括以下步骤：
选取任一目标车辆工况点，获取所述目标车辆工况点的目标动力系统扭矩；
若所述目标动力系统扭矩小于等于第一扭矩值t1，则所述目标车辆工况点的扭矩分配系数范围设为[0,1]，所述第一扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较小值；
若所述目标动力系统扭矩大于第一扭矩值t1、且小于等于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为所述第二扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较大值；
若所述目标动力系统扭矩大于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为
4.如权利要求3所述的方法，其特征在于，所述“根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率”步骤，具体包括以下步骤：
在所述扭矩分配系数范围内设置多组效率比对组，所述效率比对组中将所述目标动力系统扭矩按照不同的扭矩分配系数细分为相应的目标前电机扭矩和目标后电机扭矩；
根据所述目标车辆工况点的目标转速和所述目标前电机扭矩确定对应的目标前电机效率，根据所述目标车辆工况点的目标转速和所述目标后电机扭矩确定对应的目标后电机效率；
获取前减速器效率和后减速器效率；
根据所述目标前电机扭矩、目标后电机扭矩、所述目标前电机效率、所述目标后电机效率、所述前减速器效率和所述后减速器效率分析各效率比对组的动力系统效率；
根据所述动力系统效率确定所述最优扭矩分配系数。
5.如权利要求4所述的方法，其特征在于，所述“根据所述目标前电机扭矩、目标后电机扭矩、所述目标前电机效率、所述目标后电机效率、所述前减速器效率和所述后减速器效率分析各效率比对组的动力系统效率”步骤，具体包括以下步骤：
根据所述目标前电机扭矩得到前电机输出功率p01，根据所述目标后电机扭矩得到后电机输出功率p02；
根据所述目标前电机效率η1和所述前电机输出功率p01计算前电机输入功率p11，根据所述目标后电机效率η2和所述后电机输出功率p02计算后电机输入功率p12，
根据所述前电机输出功率p01、所述后电机输出功率p02、所述前电机输入功率p11、所述后电机输入功率p12、所述前减速器效率η3和所述后减速器效率η4分析各效率比对组的动力系统效率η，
6.一种基于python的纯电动汽车扭矩分配系数分析系统，其特征在于，包括：
数据读取处理模块，用于通过python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点，获取采样点前电机效率和采样点后电机效率，采样点为所述车辆工况点中任意多个；
工况点效率分析模块，与所述数据读取处理模块通信连接，用于通过python第三方库numpy模块、math模块、scipy模块及panda模块，根据所述采样点前电机效率和所述采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率；以及，
最优分配系数计算模块，与所述工况点效率分析模块通信连接，用于根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数，所述前电机效率包括所述采样点前电机效率和所述非采样点前电机效率，所述后电机效率包括所述采样点后电机效率和所述非采样点后电机效率。
7.如权利要求6所述的系统，其特征在于，所述数据读取处理模块包括：
数据读取单元，用于读取前电机外特性数据和后电机外特性数据；
数据处理单元，与所述数据读取单元通信连接，用于根据所述前电机外特性数据和所述后电机外特性数据分别确定前电机最大转速、后电机最大转速、任一转速下的前电机最大扭矩及任一转速下的后电机最大扭矩，选取所述前电机最大转速和所述后电机最大转速中的较小值作为动力系统最大转速，选取任一转速下的所述前电机最大扭矩和所述后电机最大扭矩作为所述任一转速下的动力系统最大扭矩；按照预设采样规则，结合所述动力系统最大转速和所述动力系统最大扭矩生成各车辆工况点，所述车辆工况点包括电机转速和及对应的动力系统扭矩；
数据获取单元，与所述数据处理单元通信连接，用于用于获取采样点前电机效率和采样点后电机效率，采样点为所述车辆工况点中任意多个。
8.如权利要求6所述的系统，其特征在于，所述系统还包括：
扭矩分配约束模块，与所述工况点效率分析模块通信连接，用于：
选取任一目标车辆工况点，获取所述目标车辆工况点的目标动力系统扭矩；
若所述目标动力系统扭矩小于等于第一扭矩值t1，则所述目标车辆工况点的扭矩分配系数范围设为[0,1]，所述第一扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较小值；
若所述目标动力系统扭矩大于第一扭矩值t1、且小于等于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为所述第二扭矩值为所述目标转速下前电机最大扭矩和后电机最大扭矩的较大值；
若所述目标动力系统扭矩大于第二扭矩值t2，则所述目标车辆工况点的扭矩分配系数范围设为
9.如权利要求8所述的系统，其特征在于，所述最优分配系数计算模块，与所述扭矩分配约束通信连接，包括：
分组单元，用于在所述扭矩分配系数范围内设置多组效率比对组，所述效率比对组中将所述目标动力系统扭矩按照不同的扭矩分配系数细分为相应的目标前电机扭矩和目标后电机扭矩；
电机效率分析单元，与所述分组单元通信连接，用于根据所述目标车辆工况点的目标转速和所述目标前电机扭矩确定对应的目标前电机效率，根据所述目标车辆工况点的目标转速和所述目标后电机扭矩确定对应的目标后电机效率；
效率获取单元，用于获取前减速器效率和后减速器效率；
系统效率分析单元，与所述电机效率分析单元和所述效率获取单元通信连接，用于根据所述目标前电机扭矩、目标后电机扭矩、所述目标前电机效率、所述目标后电机效率、所述前减速器效率和所述后减速器效率分析各效率比对组的动力系统效率；
最优系数分析单元，与所述系统效率分析单元通信连接，用于根据所述动力系统效率确定所述最优扭矩分配系数。
10.如权利要求9所述的系统，其特征在于，所述系统效率分析单元包括：
输入功率分析子单元，用于根据所述目标前电机扭矩得到前电机输出功率p01，根据所述目标后电机扭矩得到后电机输出功率p02；
输出功率分析子单元，与所述输入功率分析子单元通信连接，用于根据所述目标前电机效率η1和所述前电机输出功率p01计算前电机输入功率p11，根据所述目标后电机效率η2和所述后电机输出功率p02计算后电机输入功率p12，
系统效率分析子单元，与所述输入功率分析子单元和所述输出功率分析子单元通信连接，用于根据所述前电机输出功率p01、所述后电机输出功率p02、所述前电机输入功率p11、所述后电机输入功率p12、所述前减速器效率η3和所述后减速器效率η4分析各效率比对组的动力系统效率η，
技术总结本发明公开了一种基于Python的纯电动汽车扭矩分配系数分析方法及系统，涉及电动汽车领域，其方法包括：通过Python第三方库openpyxl模块，读取前电机外特性数据和后电机外特性数据，按照预设采样规则生成各车辆工况点；获取采样点前电机效率和采样点后电机效率；通过Python第三方库numpy模块、math模块、scipy模块及panda模块，根据采样点前电机效率和采样点后电机效率，插值得到各车辆工况点中的非采样点前电机效率和非采样点后电机效率；根据前电机效率和后电机效率分析各车辆工况点的最优扭矩分配系数。本发明对前后电机进行扭矩分配，可以保证在满足车辆动力性需求时车辆动力系统效率时刻保持最优，降低车辆能耗，提升车辆续驶里程。
技术研发人员：白志浩;赵征澜;佘建强;张丽;徐飞
受保护的技术使用者：东风汽车集团有限公司
技术研发日：2020.12.01
技术公布日：2021.03.12
