CN114151216B - 一种基于缸压重构的主动减振控制方法及装置 - Google Patents

一种基于缸压重构的主动减振控制方法及装置 Download PDF

Info

Publication number
CN114151216B
CN114151216B CN202111273303.7A CN202111273303A CN114151216B CN 114151216 B CN114151216 B CN 114151216B CN 202111273303 A CN202111273303 A CN 202111273303A CN 114151216 B CN114151216 B CN 114151216B
Authority
CN
China
Prior art keywords
combustion
axis
cylinder
current
gas
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111273303.7A
Other languages
English (en)
Other versions
CN114151216A (zh
Inventor
孙进伟
张金羽
杨福源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202111273303.7A priority Critical patent/CN114151216B/zh
Publication of CN114151216A publication Critical patent/CN114151216A/zh
Application granted granted Critical
Publication of CN114151216B publication Critical patent/CN114151216B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D41/00—Electrical control of supply of combustible mixture or its constituents
    • F02D41/02—Circuit arrangements for generating control signals
    • F02D41/14—Introducing closed-loop corrections
    • F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W20/00—Control systems specially adapted for hybrid vehicles
    • B60W20/10—Controlling the power contribution of each of the prime movers to meet required power demand
    • B60W20/15—Control strategies specially adapted for achieving a particular effect
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02B—INTERNAL-COMBUSTION PISTON ENGINES; COMBUSTION ENGINES IN GENERAL
    • F02B75/00—Other engines
    • F02B75/06—Engines with means for equalising torque
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D41/00—Electrical control of supply of combustible mixture or its constituents
    • F02D41/02—Circuit arrangements for generating control signals
    • F02D41/14—Introducing closed-loop corrections
    • F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
    • F02D41/1486—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor with correction for particular operating conditions
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D41/00—Electrical control of supply of combustible mixture or its constituents
    • F02D41/02—Circuit arrangements for generating control signals
    • F02D41/14—Introducing closed-loop corrections
    • F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
    • F02D2041/1433—Introducing closed-loop corrections characterised by the control or regulation method using a model or simulation of the system
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/60—Other road transportation technologies with climate change mitigation effect
    • Y02T10/62—Hybrid vehicles

Abstract

本发明实施例公开一种基于缸压重构的主动减振控制方法及装置，该方法包括重构发动机燃烧模型，得到第一零维燃烧模型；根据曲轴上的转速信息和发动机机体信息构建第一往复惯性力矩模型，根据第一往复惯性力矩模型和曲轴传递函数获得第一修正数据；利用第一修正数据对第一零维燃烧模型进行修正，得到第二零维燃烧模型；根据曲轴转角和第二零维燃烧模型，获得发动机的实时缸压曲线；根据实时缸压曲线获得发电机q轴前馈电流。应用本发明可实现抑制混合动力系统发动机工作过程中的转矩和转速波动，从而达到改善混合动力系统NVH的效果。

Description

一种基于缸压重构的主动减振控制方法及装置
技术领域
本发明涉及新能源应用领域，具体而言，涉及一种基于缸压重构的主动减振控制方 法及装置。
背景技术
混合动力作为新能源车用动力重要组成部分承担着重要的节能减排任务。对于车用 领域来说，小型化和一体化是轻量化实现的重要手段，也是节能减排的重要方法之一。在传统的混合动力构成中，一般是常规发动机和永磁同步电机组合的方式。车用发动机 通常是循环往复式活塞发动机，通过缸内混合气急速燃烧产生爆发压力推动活塞下行， 然后通过曲柄连杆机构将化学能转化为机械动能。曲柄连杆机构在工作的时候会分为进 气、压缩、燃烧和排气四个过程，在这个四个过程中，尤其是在燃烧时缸内压力急速上 升，活塞迅速下移，循环往复从而完成做功。在曲轴旋转一圈的过程中自然的会产生一 个转速和扭矩的脉动，脉动造成的车辆振动，降低了驾驶和乘车的舒适感。如何对发动 机进行减振，提高驾乘舒适度，是新能源应用领域内需要解决的问题。
发明内容
本发明提供了一种基于缸压重构的主动减振控制方法及装置，以实现对车用混合动 力系统主动减振的作用。具体的技术方案如下：
第一方面，本发明实施例提供了一种基于缸压重构的主动减振控制方法，所述方法 包括：
根据热力学第一定律推导得到混合动力汽车发动机的初始零维燃烧模型，所述发动 机气缸为闭口系统，所述缸内气体为理想气体，所述初始零维燃烧模型公式如下：
Figure BDA0003329358660000011
其中，
Figure BDA0003329358660000012
为燃烧过程中曲轴瞬时转角，
Figure BDA0003329358660000013
为燃烧放热率，
Figure BDA0003329358660000014
为排气过程中随排气工 质加入到燃烧室的热量变化率，
Figure BDA0003329358660000021
Figure BDA0003329358660000022
为进气过程中随排气工质加入到燃烧室的热 量变化率，
Figure BDA0003329358660000023
Figure BDA0003329358660000024
为通过燃烧室壁散热的散热率，
Figure BDA0003329358660000025
为工质所做机械功率的变化率，
Figure BDA0003329358660000026
P为缸内气体压力，
Figure BDA0003329358660000027
为气缸容积对时间的变化率，
Figure BDA0003329358660000028
为工质内能的变 化率，
Figure BDA0003329358660000029
mg为缸内气体质量，且mg为常数，Cv为缸内气体定容比热容，Tgas为缸内气体温度，
Figure BDA00033293586600000210
为缸内气体温度对曲轴转角的变化率；
对理想气体状态方程进行求导，得到所述发动机的第一理想气体状态公式：
PdV+Vdp＝mgRadTgas
式中，P为缸内气体压力，V为缸内容积，mg为缸内气体质量，且mg为常数，Ra为 空气气体常数，Tgas为缸内气体温度；
根据所述初始零维燃烧模型、所述第一理想气体状态公式、第一关系式，得出第一零维燃烧模型，所述第一关系式是Cv和Ra之间的关系式，所述Cv和Ra之间的关系式为Cv＝Ra/(k-1)，k为绝热指数，所述第一零维燃烧模型公式如下：
Figure BDA00033293586600000211
其中，由Woschni提出的经验公式得出Qw＝hc(Tgas-Tw)，Tw为壁面温度，hc为热 传导系数，
Figure BDA00033293586600000212
ηu为燃烧效率，
Figure BDA00033293586600000213
为燃烧速度，Qf为燃烧放热量，Hu为燃 烧的低热值，x为已燃燃油百分比，
Figure BDA00033293586600000214
x由Weibo模型可得，c 为燃烧效率系数，m为燃烧品质指数，
Figure BDA00033293586600000215
为燃烧始点对应的曲轴转角，
Figure BDA00033293586600000216
Figure BDA00033293586600000217
为燃烧终点对应的曲轴转角，
Figure BDA00033293586600000218
Figure BDA00033293586600000219
为喷油或点火始点所对应的曲轴转角， τid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差，
Figure BDA00033293586600000220
为燃烧持续期；
根据曲轴上的转速信息和发动机机体信息构建第一往复惯性力矩模型，对所述第一 往复惯性力矩模型进行快速FFT变换，获得第二往复惯性力矩模型；
提取所述第二往复惯性力矩模型中的特征值，得到第一特征数据，提取曲轴传递函 数中的特征值，得到第二特征数据，分别将所述第一特征数据中的特征值与第二特征数据中与第一特征数据中对应的特征值相减，得到第一修正数据；
利用所述第一修正数据对所述第一零维燃烧模型进行修正，得到第二零维燃烧模型， 所述第二零维燃烧模型为参数化的燃烧模型；
根据曲轴转角和所述第二零维燃烧模型，获得发动机的实时缸压曲线；
根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通过查表方式获得，所述q 轴前馈电流公式如下：
Figure BDA0003329358660000031
式中，
Figure BDA0003329358660000032
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得， PDem为目标需求功率，PDC为直流母线上总功率。
可选的，所述根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通过查表方式 获得，所述q轴前馈电流公式如下：
Figure BDA0003329358660000033
式中，
Figure BDA0003329358660000034
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得， PDem为目标需求功率，PDC为直流母线上总功率，包括：
根据车辆实车参数和目标车速，得到目标需求功率PDem；
将缸内气体压力P、目标需求功率PDem、直流母线上总功率PDC、q轴前馈电流组成 多项式函数，对所述多项式函数进行离散化，获得MAP表，所述缸内气体压力由所述实 时缸压曲线获得；
查找所述MAP表，得到q轴前馈电流。
可选的，还包括根据q轴前馈电流和q轴原始电流，确定q轴目标电流，所述q轴 目标电流是q轴原始电流与q轴前馈电流抵消后的q轴电流；
其中，通过永磁同步电机控制矢量控制算法得到q轴原始电流，公式如下：
Figure BDA0003329358660000035
式中，uq为q轴电压，iq为q轴原始电流，Lq为q轴电感，Rs为转子电阻，ωr为电 角度，ψm为磁链系数。
可选的，所述永磁同步电机的d轴前馈电流为0。
可选的，还包括根据发动机往复惯性力和发电机转子机械转动惯量，确定发电机转 子目标转矩公式：
Figure BDA0003329358660000036
式中，J是发电机转子机械转动惯量，Te为发电机转子目标转矩，Bm为含有发动机往复惯性力的机械阻尼系数，ωr为电角度，ωm为机械角度。
第二方面，本发明实施例提供了一种基于缸压重构的主动减振控制装置，所述装置 包括；
第一建模模块，用于根据热力学第一定律推导得到混合动力汽车发动机的初始零维 燃烧模型，所述发动机气缸为闭口系统，所述缸内气体为理想气体，所述初始零维燃烧模型公式如下：
Figure BDA0003329358660000041
其中，
Figure BDA0003329358660000042
为燃烧过程中曲轴瞬时转角，
Figure BDA0003329358660000043
为燃烧放热率，
Figure BDA0003329358660000044
为排气过程中随排气工 质加入到燃烧室的热量变化率，
Figure BDA0003329358660000045
Figure BDA0003329358660000046
为进气过程中随排气工质加入到燃烧室的热 量变化率，
Figure BDA0003329358660000047
Figure BDA0003329358660000048
为通过燃烧室壁散热的散热率，
Figure BDA0003329358660000049
为工质所做机械功率的变化率，
Figure BDA00033293586600000410
P为缸内气体压力，
Figure BDA00033293586600000411
为气缸容积对时间的变化率，
Figure BDA00033293586600000412
为工质内能的变 化率，
Figure BDA00033293586600000413
mg为缸内气体质量，且mg为常数，Cv为缸内气体定容比热容， Tgas为缸内气体温度，
Figure BDA00033293586600000414
为缸内气体温度对曲轴转角的变化率；
第一获得模块，用于对理想气体状态方程进行求导，得到所述发动机的第一理想气 体状态公式：
PdV+Vdp＝mgRadTgas
式中，P为缸内气体压力，V为缸内容积，mg为缸内气体质量，且mg为常数，Ra为 空气气体常数，Tgas为缸内气体温度；
第二建模模块，用于根据所述初始零维燃烧模型、所述第一理想气体状态公式、第一关系式，得出第一零维燃烧模型，所述第一关系式是Cv和Ra之间的关系式，所述Cv和 Ra之间的关系式为Cv＝Ra/(k-1)，k为绝热指数，所述第一零维燃烧模型公式如下：
Figure BDA00033293586600000415
其中，由Woschni提出的经验公式得出Qw＝hc(Tgas-Tw)，Tw为壁面温度，hc为热 传导系数，
Figure BDA00033293586600000416
ηu为燃烧效率，
Figure BDA00033293586600000417
为燃烧速度，Qf为燃烧放热量，Hu为燃 烧的低热值，x为已燃燃油百分比，
Figure BDA00033293586600000418
x由Weibo模型可得，c 为燃烧效率系数，m为燃烧品质指数，
Figure BDA00033293586600000419
为燃烧始点对应的曲轴转角，
Figure BDA00033293586600000420
Figure BDA00033293586600000421
为燃烧终点对应的曲轴转角，
Figure BDA00033293586600000422
Figure BDA00033293586600000423
为喷油或点火始点所对应的曲轴转角， τid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差，
Figure BDA0003329358660000051
为燃烧持续期；
第三建模模块，用于根据曲轴上的转速信息和发动机机体信息构建第一往复惯性力 矩模型，对所述第一往复惯性力矩模型进行快速FFT变换，获得第二往复惯性力矩模型；
特征提取模块，用于提取所述第二往复惯性力矩模型中的特征值，得到第一特征数 据，提取曲轴传递函数中的特征值，得到第二特征数据，分别将所述第一特征数据中的特征值与第二特征数据中与第一特征数据中对应的特征值相减，得到第一修正数据；
模型修正模块，用于利用所述第一修正数据对所述第一零维燃烧模型进行修正，得 到第二零维燃烧模型，所述第二零维燃烧模型为参数化的燃烧模型；
第二获得模块，用于根据根据曲轴转角和所述第二零维燃烧模型，获得发动机的实 时缸压曲线；
第三获得模块，用于根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通过查 表方式获得，所述q轴前馈电流公式如下：
Figure BDA0003329358660000052
式中，
Figure BDA0003329358660000053
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得， PDem为目标需求功率，PDC为直流母线上总功率。
可选的，所述第三获得模块包括：
确定模块，用于根据车辆实车参数和目标车速，确定目标需求功率PDem；
建表模块，用于将缸内气体压力P、目标需求功率PDem、直流母线上总功率PDC、q 轴前馈电流组成多项式函数，对所述多项式函数进行离散化，获得MAP表，所述缸内气 体压力由所述实时缸压曲线获得；
查表模块，用于查找所述MAP表，得到q轴前馈电流。
可选的，还包括q轴目标电流确定模块，用于根据q轴前馈电流和q轴原始电流， 确定q轴目标电流，所述q轴目标电流是q轴原始电流与q轴前馈电流抵消后的q轴电 流；
其中，通过永磁同步电机控制矢量控制算法得到q轴原始电流，公式如下：
Figure BDA0003329358660000054
式中，uq为q轴电压，iq为q轴原始电流，Lq为q轴电感，Rs为转子电阻，ωr为电 角度，ψm为磁链系数。
可选的，所述永磁同步电机的d轴前馈电流为0。
可选的，还包括目标转矩确定模块，用于根据发动机往复惯性力和发电机转子机械 转动惯量，确定发电机转子目标转矩公式：
Figure BDA0003329358660000061
式中，J是发电机转子机械转动惯量，Te为发电机转子目标转矩，Bm为含有发动机往复惯性力的机械阻尼系数，ωr为电角度，ωm为机械角度。
由上述内容可知，本发明实施例提供的一种基于缸压重构的主动减振控制方法及装 置，该方法包括重构发动机燃烧模型，得到第一零维燃烧模型；根据曲轴上的转速信息和发动机机体信息构建第一往复惯性力矩模型，对第一往复惯性力矩模型进行快速FFT 变换，获得第二往复惯性力矩模型；提取第二往复惯性力矩模型中的特征值，得到第一 特征数据，提取曲轴传递函数中的特征值，得到第二特征数据，分别将第一特征数据中 的特征值与第二特征数据中与第一特征数据中对应的特征值相减，得到第一修正数据； 利用第一修正数据对第一零维燃烧模型进行修正，得到第二零维燃烧模型，第二零维燃 烧模型为参数化的燃烧模型；根据曲轴转角和第二零维燃烧模型，获得发动机的实时缸 压曲线；根据实时缸压曲线获得发电机q轴前馈电流。
应用本发明实施例，可基于现有的传感器信号，采用模型重构的方法得到发动机起 始燃烧时刻及燃烧爆发压力，从而控制电机输出谐波转矩方法来改善系统整个轴系上转 矩和转速波动，另外通过查表方式获得q轴前馈电流，简化了计算过程。当然，实施本 发明的任一产品或方法并不一定需要同时达到以上所述的所有优点。
本发明实施例的创新点包括：
1、车用混合动力系统一般由常规发动机和永磁同步电机组合方式构成，发动机中设 置有飞轮，飞轮在整个混合动力系统中起减振、储能和滤波的作用。本发明实施例将发动机飞轮取消，在发动机曲轴旋转的过程中自然的会产生转速和扭矩的脉动，取消飞轮 该脉动会尤为明显。本发明实施例通过在发电机的转子上施加一个与上述脉动相反的脉 动转矩，通过反向脉动转矩来抵消没有飞轮的发动机曲轴上的转矩和转速波动，使得整 个系统在三维空间中工作更加平顺，从而达到改善混合动力系统NVH(NoiseVibrationHarshness噪声、振动与声振粗糙度)的效果，起到主动减振的作用。
2、本发明实施例系统由无飞轮的循环往复式活塞发动机、永磁同步电机、曲轴位置 传感器、凸轮轴位置传感器、旋转变压器以及ECU(Electronic Control Unit发动机控制 器)和MCU(Micro Controller Unit电机控制器)等构成，本发明实施例基于现有的传感器，利用高精度的旋变信号和发动机相位信号等来重构发动机缸内燃烧模型，得到发动 机经过燃烧施加在扭矩上的转矩波动，将这个转矩作为ISG电机谐波转矩补偿的前馈值， 从而减少整个系统上的转矩波动。
3、本发明实施例通过重构发动机缸内燃烧模型，可得到发动机缸内压力的变化情况， 即发动机实时缸压曲线，通过缸压曲线可获得ISG电机转子转矩前馈目标值，从而来达到削峰填谷的效果。本发明实施例中通过实时缸压曲线获得缸内气体压力，将含有缸内 气体压力、目标需求功率PDem、直流母线上总功率PDC、q轴前馈电流的多项式函数进行 离散化，最终共同组成一张MAP表，所以q轴前馈电流可直接通过查找MAP的方式获 得，极大的简化了计算过程。
4、现有技术中，由整车控制器(HCU)发送给电机控制器(MCU)的发电机转子 转子转矩不含有抵消发动机产生的转速和扭矩的脉动，本发明实施例中由整车控制器 (HCU)发送给电机控制器(MCU)的发电机转子目标转矩数据Te则由原始发电机转子 转矩和综合机械阻尼获得，其中Bm中为含有发动机往复惯性力的机械阻尼系数，本发明 实施例抵消了发动机产生的转速和扭矩的脉动，应用本发明实施例可改善混合动力系统 整个轴系上转矩和转速的波动。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有 技术描述中所需要使用的附图作简单介绍。显而易见地，下面描述中的附图仅仅是本发明的一些实施例。对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可 以根据这些附图获得其他的附图。
图1为本发明实施例提供的基于缸压重构的主动减振控制方法的一种流程示意图；
图2为本发明实施例提供的某一转速下的电机转矩补偿示意图；
图3为本发明实施例提供的基于缸压重构的主动减振控制装置的一种结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整 的描述。显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，本发明实施例及附图中的术语“包括”和“具有”以及它们的任何变形， 意图在于覆盖不排他的包含。例如包含的一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或 可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
本发明提供了一种基于缸压重构的主动减振控制方法及装置。下面对本发明实施例 进行详细说明。
图1为本发明实施例提供的基于缸压重构的主动减振控制方法的一种流程示意图。 该方法可以包括如下步骤：
S101：根据热力学第一定律推导得到混合动力汽车发动机的初始零维燃烧模型，所 述发动机气缸为闭口系统，所述缸内气体为理想气体，所述初始零维燃烧模型公式如下：
Figure BDA0003329358660000081
其中，
Figure BDA00033293586600000814
为燃烧过程中曲轴瞬时转角，
Figure BDA0003329358660000082
为燃烧放热率，
Figure BDA0003329358660000083
为排气过程中随排气工 质加入到燃烧室的热量变化率，
Figure BDA0003329358660000084
Figure BDA0003329358660000085
为进气过程中随排气工质加入到燃烧室的热 量变化率，
Figure BDA0003329358660000086
Figure BDA0003329358660000087
为通过燃烧室壁散热的散热率，
Figure BDA0003329358660000088
为工质所做机械功率的变化率，
Figure BDA0003329358660000089
P为缸内气体压力，
Figure BDA00033293586600000810
为气缸容积对时间的变化率，
Figure BDA00033293586600000811
为工质内能的变 化率，
Figure BDA00033293586600000812
mg为缸内气体质量，且mg为常数，Cv为缸内气体定容比热容， Tgas为缸内气体温度，
Figure BDA00033293586600000813
为缸内气体温度对曲轴转角的变化率。
在传统的混合动力构成中，一般是常规发动机和永磁同步电机组合的方式，常规车 用发动机通常是循环往复式活塞发动机，循环往复式活塞发动机通过缸内混合气急速燃 烧产生爆发压力推动活塞下行，然后通过曲柄连杆机构将化学能转化为机械动能。曲柄连杆机构在工作的时候会分为进气、压缩、燃烧和排气四个过程，循环往复这四个过程 从而完成做功。在曲轴旋转一圈的过程中自然的会产生一个转速和扭矩的脉动，尤其是 在燃烧时缸内压力急速上升，活塞迅速下移时，发动机产生的扭矩最大。车用发动机中 的飞轮会抵消一部分该脉动，飞轮在发动机工作过程中起着减振、储能和滤波的作用。 本发明实施例将发动机飞轮取消，利用电机转子来承接一部分飞轮的作用。
在发动机曲轴旋转一圈的过程中自然的会产生一个转速和扭矩的脉动，本发明实施 例将发动机飞轮取消，取消飞轮之后使得这种脉动变得尤为明显。本发明实施例主要是通过在电机转子上面施加一个与发动机脉动相反的脉动转矩，通过反向脉动转矩来抵消没有飞轮的发动机曲轴上的转矩和转速波动，使得整个系统在三维空间中工作更加平顺。在一个可实现的方式中，本发明实施例将取消飞轮的循环往复式活塞发动机和ISG(IntegrateStartGenerator集成启动/发电一体化电机)电机直连。在ISG电机中注入谐波转矩，以抵消发动机曲轴上的转矩和转速波动。如图2所示的本发明实施例提供的某一 转速下的电机转矩补偿示意图。需要说明的是图2中的电机谐波转矩即为电机转子上面 施加的与发动机脉动相反的脉动转矩。
本发明实施例通过建立发动机零维燃烧模型，得到发动机缸内工作状态，对于发火 顺序为1-3-4-2四缸机来说，曲轴每转两圈发动机完成一个完整的循环，即曲轴每转一圈 就有两个缸完成一个进气、压缩、做功、排气四个循环。
本发明实施例中，根据热力学第一定律得到初始零维燃烧模型，热力学第一定律是 涉及热现象领域内的能量守恒和转化定律，反映了不同形式的能量在传递与转换过程中 守恒。本申请中，通过建立初始零维燃烧模型重新构建了发动机缸压，在缸压重构过程中为了减小算法的复杂程度，假设发动机气缸内为一个闭口系统，因此，
Figure BDA0003329358660000091
Figure BDA0003329358660000092
Figure BDA0003329358660000093
将缸内气体假设为理想气体，其内能的变化率为：
Figure BDA0003329358660000094
其中，mg为缸内气体质量，Cv为缸内气体定容比热容，Tgas为气体温度，
Figure BDA0003329358660000095
为缸内气体温度对曲轴转角的变化率。
S102：对理想气体状态方程进行求导，得到所述发动机的第一理想气体状态公式：
PdV+Vdp＝mgRadTgas
式中，P为缸内气体压力，V为缸内容积，mg为缸内气体质量，且mg为常数，Ra为 空气气体常数，Tgas为缸内气体温度。
本发明实施例中的气缸为闭口系统，闭口系统中气体质量mg为常数，可通过参数标 定的方式实现，Ra可通过查表得到。
S103：根据所述初始零维燃烧模型、所述第一理想气体状态公式、第一关系式，得出第一零维燃烧模型，所述第一关系式是Cv和Ra之间的关系式，所述Cv和Ra之间的关系 式为Cv＝Ra/(k-1)，k为绝热指数，所述第一零维燃烧模型公式如下：
Figure BDA0003329358660000096
其中，由Woschni提出的经验公式得出Qw＝hc(Tgas-Tw)，Tw为壁面温度，hc为热 传导系数，
Figure BDA0003329358660000097
ηu为燃烧效率，
Figure BDA0003329358660000098
为燃烧速度，Qf为燃烧放热量，Hu为燃 烧的低热值，x为已燃燃油百分比，
Figure BDA0003329358660000101
x由Weibo模型可得，c 为燃烧效率系数，m为燃烧品质指数，
Figure BDA0003329358660000102
为燃烧始点对应的曲轴转角，
Figure BDA0003329358660000103
Figure BDA0003329358660000104
为燃烧终点对应的曲轴转角，
Figure BDA0003329358660000105
Figure BDA0003329358660000106
为喷油或点火始点所对应的曲轴转角， ωid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差，
Figure BDA0003329358660000107
为燃烧持续期。
绝热指数k可以通过标定得到，壁面温度Tw也可通过数据标定得到。利用Weibo(韦柏)模型可得已燃燃油百分比：
Figure BDA0003329358660000108
其中，
Figure BDA0003329358660000109
为燃烧始点对应的曲轴转角；
Figure BDA00033293586600001010
为燃烧终点对应的曲轴转角；c为燃烧效率 系数；m为燃烧品质指数。
在实际发动机运行过程中，燃烧始点、燃烧终点、燃烧品质指数很难确定，因此，考虑到实际应用，可采用喷油或点火始点所对应的曲轴转角
Figure BDA00033293586600001011
替代上述参数，
Figure BDA00033293586600001012
可通过 电机旋变信号获得，因此，
Figure BDA00033293586600001013
可以通过以下关系得到：
Figure BDA00033293586600001014
Figure BDA00033293586600001015
其中，τid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差；
Figure BDA00033293586600001016
为燃烧持续期， 这两个参数都可以通过发动机台架标定得到。
S104：根据曲轴上的转速信息和发动机机体信息构建第一往复惯性力矩模型，对所 述第一往复惯性力矩模型进行快速FFT变换，获得第二往复惯性力矩模型。
往复式活塞式发动机作用在曲轴上的作用力主要为气体作用力、往复惯性力、摩擦 转矩、电机轴上的负载转矩等。根据研究发现摩擦转矩虽然也与发火频率为周期，但其幅值非常小，所以其低次谐波分量相对于气体作用力矩、往复惯性力矩的同次分量来说 可以忽略不计。因此，发动机总转矩只是气体往复作用力矩和惯性力矩的合成，是频率 为4Cz的周期性信号，所以发动机总转矩即是往复惯性力矩。利用曲轴上的转速信息和 发动机机体信息构建往复惯性力矩模型，对往复惯性力矩模型进行快速FFT变换(fast Fouriertransform快速傅里叶变换)后，提取特征值，提取到的特征值用于修正第一零维 燃烧模型。
S105：提取所述第二往复惯性力矩模型中的特征值，得到第一特征数据，提取曲轴传递函数中的特征值，得到第二特征数据，分别将所述第一特征数据中的特征值与第二 特征数据中与第一特征数据中对应的特征值相减，得到第一修正数据。
从第二往复惯性力矩模型中提取特征值，得到第一特征数据；依据发动机构建曲轴 传递函数，从曲轴传递函数中提取特征值，得到第二特征数据，需要说明的是，第一特征数据中的特征与第二特征数据中的特征一一对应，按照对应特征将两组特征数据中特征值相减，得到第一修正数据，该第一修正值用于修正第一零维燃烧模型。
S106：利用所述第一修正数据对所述第一零维燃烧模型进行修正，得到第二零维燃 烧模型，所述第二零维燃烧模型为参数化的燃烧模型。
利用特征值对第一零维燃烧模型进行修正，得到最终的发动机燃烧模型，即第二零 维燃烧模型。
S107：根据曲轴转角和所述第二零维燃烧模型，获得发动机的实时缸压曲线。
根据喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差ωid、喷油或点火始点所对 应的曲轴转角
Figure BDA0003329358660000111
确定燃烧始点对应的曲轴转角
Figure BDA0003329358660000112
根据燃烧持续期
Figure BDA0003329358660000113
喷油或点火始点所对应的曲轴转角
Figure BDA0003329358660000114
确定燃烧终点对应的曲轴转角
Figure BDA0003329358660000115
根据
Figure BDA0003329358660000116
确定燃烧放 热率
Figure BDA0003329358660000117
根据燃烧放热率、缸内容积，以及第二零维燃烧模型，可获得发动机的实时 缸压曲线。
S108：根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通过查表方式获得， 所述q轴前馈电流公式如下：
Figure BDA0003329358660000118
式中，
Figure BDA0003329358660000119
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得， PDem为目标需求功率，PDC为直流母线上总功率。
通过参数化的燃烧模型即第二零维燃烧模型，得到发动机的实时缸压曲线，离散化 实时缸压曲线，将离散化的缸压值P、需求功率PDem、总线功率PDC、q轴电流共同组成 一个多项式函数：
Figure BDA00033293586600001110
其中，
Figure BDA00033293586600001111
为q轴前馈电流，P为实时缸内气体压力，PDem为车辆实际需要的目标需求功率，PDC为直流母线上总功率。目标需求功率PDem通过车辆实车参数及目标车速可 以得到，属于一个已知变量。PDC可通过公式PDC＝uDC*iDC计算得到，uDC为直流母线 端电压，iDC为直流母线端和电流，PDC所以也属于已知变量。为了降低模型运算的复杂 度，将上述多项式函数进行离散化，共同组成一张MAP表，通过这个MAP可以得到q 轴电流的前馈值。
在一个可选的实施例中，所述根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电 流通过查表方式获得，所述q轴前馈电流公式如下：
Figure BDA0003329358660000121
式中，
Figure BDA0003329358660000122
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得， PDem为目标需求功率，PDC为直流母线上总功率，包括：
根据车辆实车参数和目标车速，确定目标需求功率PDem。
将缸内气体压力P、目标需求功率PDem、直流母线上总功率PDC、q轴前馈电流组成 多项式函数，对所述多项式函数进行离散化，获得MAP表，所述缸内气体压力由所述实 时缸压曲线获得。
查找所述MAP表，得到q轴前馈电流。
直流母线上总功率PDC可通过直流母线上电压和电流获得，对多项式函数进行离散化， 共同组成一张MAP表，通过查找MAP表的方式得到q轴前馈电流，极大的降低了运算 的复杂度，简化了缸压主动减振方法，需要说明的是，该多项式函数即为q轴前馈电流 公式。
在上述实施例的基础上，还包括根据q轴前馈电流和q轴原始电流，确定q轴目标电流，所述q轴目标电流是q轴原始电流与q轴前馈电流抵消后的q轴电流；
其中，通过永磁同步电机控制矢量控制算法得到q轴原始电流，公式如下：
Figure BDA0003329358660000123
式中，uq为q轴电压，iq为q轴原始电流，Lq为q轴电感，Rs为转子电阻，ωr为电 角度，ψm为磁链系数。
q轴原始电流是不向发电机的q轴施加q轴前馈电流时的q轴电流，uq是不向发电机的q轴施加q轴前馈电流时的q轴电压，q轴目标电流是q轴原始电流抵消q轴前馈电流 后实际的q轴上的电流。ωr电角度通过电机旋变信号获得。
在一个可选的实施例中，所述永磁同步电机的d轴前馈电流为0。
电机的交轴也叫q轴，直轴也叫d轴，q轴和d轴是坐标轴而不是实际的轴。在永磁同步电机控制中，为了能够得到类似直流电机的控制特性，在电机转子上建立了一个坐 标系，此坐标系与转子同步转动，取转子磁场方向为d轴，垂直于转子磁场方向为q轴， 将电机的数学模型转换到此坐标系下，可实现d轴和q轴的解耦，从而得到良好控制特 性。在本发明实施例中，默认d轴前馈电流为0。
通过永磁同步电机控制矢量控制算法可得到d轴原始电流，公式如下：
Figure BDA0003329358660000131
其中，ud为d轴电压，id为d轴电流，Ld为d轴电感，Rs为转子电阻，ωr为电角度。
d轴原始电流是不向发电机的d轴施加d轴前馈电流时的d轴电流，ud是不向发电机的d轴施加d轴前馈电流时的d轴电压，d轴目标电流是d轴原始电流抵消d轴前馈电流 后，实际的d轴上的电流。因为d轴前馈电流为0，所以d轴目标电流数值等于d轴原始 电流数值。
在上述实施例的基础上，还包括根据发动机往复惯性力和发电机转子机械转动惯量， 确定发电机转子目标转矩公式：
Figure BDA0003329358660000132
式中，J是发电机转子机械转动惯量，Te为发电机转子目标转矩，Bm为含有发动机往复惯性力的机械阻尼系数，ωr为电角度，ωm为机械角度。
其中，Bm是一个综合系数，它里面含有发动机往复惯性力的机械阻尼系数，Bm可通过经验公式和数据标定相结合的方式获得具体数据，Te是由整车控制器(HCU)发送给 电机控制器(MCU)的发电机转子目标转矩数据，该发电机转子目标转矩抵消了发动机 经过燃烧施加在扭矩上的转矩波动，起到了主动减振的作用。
应用本发明实施例，通过重构发动机缸内燃烧模型，得到发动机缸内压力的变化情 况，即发动机缸压曲线，将重构的缸压曲线作为ISG电机转子转矩前馈目标值，从而来达到削峰填谷的效果。
相应于上述方法实施例，本发明实施例提供了一种基于缸压重构的主动减振控制装 置，其特征在于，所述装置包括：
第一建模模块S301，用于根据热力学第一定律推导得到混合动力汽车发动机的初始 零维燃烧模型，所述发动机气缸为闭口系统，所述缸内气体为理想气体，所述初始零维燃烧模型公式如下：
Figure BDA0003329358660000133
其中，
Figure BDA0003329358660000134
为燃烧过程中曲轴瞬时转角，
Figure BDA0003329358660000135
为燃烧放热率，
Figure BDA0003329358660000136
为排气过程中随排气工 质加入到燃烧室的热量变化率，
Figure BDA0003329358660000137
Figure BDA0003329358660000138
为进气过程中随排气工质加入到燃烧室的热 量变化率，
Figure BDA0003329358660000139
Figure BDA00033293586600001310
为通过燃烧室壁散热的散热率，
Figure BDA00033293586600001311
为工质所做机械功率的变化率，
Figure BDA0003329358660000141
P为缸内气体压力，
Figure BDA0003329358660000142
为气缸容积对时间的变化率，
Figure BDA0003329358660000143
为工质内能的变 化率，
Figure BDA0003329358660000144
mg为缸内气体质量，且mg为常数，Cv为缸内气体定容比热容， Tgas为缸内气体温度，
Figure BDA0003329358660000145
为缸内气体温度对曲轴转角的变化率。
本发明实施例通过建立发动机零维燃烧模型，得到发动机缸内工作状态，对于发火 顺序为1-3-4-2四缸机来说，曲轴每转两圈发动机完成一个完整的循环，即曲轴每转一圈 就有两个缸完成一个进气、压缩、做功、排气四个循环。
本发明实施例中，根据热力学第一定律得到初始零维燃烧模型，热力学第一定律是 涉及热现象领域内的能量守恒和转化定律，反映了不同形式的能量在传递与转换过程中 守恒。本申请中，通过建立初始零维燃烧模型重新构建了发动机缸压，在缸压重构过程中为了减小算法的复杂程度，假设发动机气缸内为一个闭口系统，因此，
Figure BDA0003329358660000146
Figure BDA0003329358660000147
Figure BDA0003329358660000148
将缸内气体假设为理想气体，其内能的变化率为：
Figure BDA0003329358660000149
其中，mg为缸内气体质量，Cv为缸内气体定容比热容，Tgas为气体温度，
Figure BDA00033293586600001410
为缸 内气体温度对曲轴转角的变化率。
第一获得模块S302，用于对理想气体状态方程进行求导，得到所述发动机的第一理 想气体状态公式：
PdV+Vdp＝mgRadTgas
式中，P为缸内气体压力，V为缸内容积，mg为缸内气体质量，且mg为常数，Ra为 空气气体常数，Tgas为缸内气体温度。
本发明实施例中的气缸为闭口系统，闭口系统中气体质量mg为常数，可通过参数标 定的方式实现，Ra可通过查表得到。
第二建模模块S303，用于根据所述初始零维燃烧模型、所述第一理想气体状态公式、 第一关系式，得出第一零维燃烧模型，所述第一关系式是Cv和Ra之间的关系式，所述Cv和Ra之间的关系式为Cv＝Ra/(k-1)，k为绝热指数，所述第一零维燃烧模型公式如下：
Figure BDA00033293586600001411
其中，由Woschni提出的经验公式得出Qw＝hc(Tgas-Tw)，Tw为壁面温度，hc为热 传导系数，
Figure BDA0003329358660000151
ηu为燃烧效率，
Figure BDA0003329358660000152
为燃烧速度，Qf为燃烧放热量，Hu为燃 烧的低热值，x为已燃燃油百分比，
Figure BDA0003329358660000153
x由Weibo模型可得，c 为燃烧效率系数，m为燃烧品质指数，
Figure BDA0003329358660000154
为燃烧始点对应的曲轴转角，
Figure BDA0003329358660000155
Figure BDA0003329358660000156
为燃烧终点对应的曲轴转角，
Figure BDA0003329358660000157
Figure BDA0003329358660000158
为喷油或点火始点所对应的曲轴转角， τid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差，
Figure BDA0003329358660000159
为燃烧持续期。
绝热指数k可以通过标定得到，壁面温度Tw也可通过数据标定得到。利用Weibo(韦柏)模型可得已燃燃油百分比：
Figure BDA00033293586600001510
其中，
Figure BDA00033293586600001511
为燃烧始点对应的曲轴转角；
Figure BDA00033293586600001512
为燃烧终点对应的曲轴转角；c为燃烧效率 系数；m为燃烧品质指数。
在实际发动机运行过程中，燃烧始点、燃烧终点、燃烧品质指数很难确定，因此，考虑到实际应用，可采用喷油或点火始点所对应的曲轴转角
Figure BDA00033293586600001513
替代上述参数，因此，
Figure BDA00033293586600001514
可以通过以下关系得到：
Figure BDA00033293586600001515
Figure BDA00033293586600001516
其中，τid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差；
Figure BDA00033293586600001517
为燃烧持续期， 这两个参数都可以通过发动机台架标定得到。
第三建模模块S304，用于根据曲轴上的转速信息和发动机机体信息构建第一往复惯 性力矩模型，对所述第一往复惯性力矩模型进行快速FFT变换，获得第二往复惯性力矩模型。
往复式活塞式发动机作用在曲轴上的作用力主要为气体作用力、往复惯性力、摩擦 转矩、电机轴上的负载转矩等。根据研究发现摩擦转矩虽然也与发火频率为周期，但其幅值非常小，所以其低次谐波分量相对于气体作用力矩、往复惯性力矩的同次分量来说 可以忽略不计。因此，发动机总转矩只是气体往复作用力矩和惯性力矩的合成，是频率 为4Cz的周期性信号，所以发动机总转矩即是往复惯性力矩。利用曲轴上的转速信息和 发动机机体信息构建往复惯性力矩模型，对往复惯性力矩模型进行快速FFT变换后，提 取特征值，提取到的特征值用于修正第一零维燃烧模型。
特征提取模块S305，用于提取所述第二往复惯性力矩模型中的特征值，得到第一特 征数据，提取曲轴传递函数中的特征值，得到第二特征数据，分别将所述第一特征数据中的特征值与第二特征数据中与第一特征数据中对应的特征值相减，得到第一修正数据。
从第二往复惯性力矩模型中提取特征值，得到第一特征数据；依据发动机构建曲轴 传递函数，从曲轴传递函数中提取特征值，得到第二特征数据，需要说明的是，第一特征数据中的特征与第二特征数据中的特征一一对应，按照对应特征将两组特征数据中特征值相减，得到第一修正数据，该第一修正值用于修正第一零维燃烧模型。
模型修正模块S306，用于利用所述第一修正数据对所述第一零维燃烧模型进行修正， 得到第二零维燃烧模型，所述第二零维燃烧模型为参数化的燃烧模型。
利用特征值对第一零维燃烧模型进行修正，得到最终的发动机燃烧模型，即第二零 维燃烧模型。
第二获得模块S307，用于根据根据曲轴转角和所述第二零维燃烧模型，获得发动机 的实时缸压曲线。
根据喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差τid、喷油或点火始点所对 应的曲轴转角
Figure BDA0003329358660000161
确定燃烧始点对应的曲轴转角
Figure BDA0003329358660000162
根据燃烧持续期
Figure BDA0003329358660000163
喷油或点火始点所对应的曲轴转角
Figure BDA0003329358660000164
确定燃烧终点对应的曲轴转角
Figure BDA0003329358660000165
根据
Figure BDA0003329358660000166
确定燃烧放 热率
Figure BDA0003329358660000167
根据燃烧放热率、缸内容积，以及第二零维燃烧模型，可获得发动机的实时 缸压曲线。
第三获得模块S308，用于根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通 过查表方式获得，所述q轴前馈电流公式如下：
Figure BDA0003329358660000168
式中，
Figure BDA0003329358660000169
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得， PDem为目标需求功率，PDC为直流母线上总功率。
通过参数化的燃烧模型即第二零维燃烧模型，得到发动机的实时缸压曲线，离散化 实时缸压曲线，将离散化的缸压值P、需求功率PDem、总线功率PDC、q轴电流共同组成 一个多项式函数：
Figure BDA00033293586600001610
其中，
Figure BDA00033293586600001611
为q轴前馈电流，P为实时缸内气体压力，PDem为车辆实际需要的目标需求功率，PDC为直流母线上总功率。目标需求功率PDem通过车辆实车参数及目标车速可 以得到，属于一个已知变量。PDC可通过公式PDC＝uDC*iDC计算得到，uDC为直流母线 端电压，iDC为直流母线端和电流，PDC所以也属于已知变量。为了降低模型运算的复杂 度，将上述多项式函数进行离散化，共同组成一张MAP表，通过这个MAP可以得到q 轴电流的前馈值。
在一个可选的实施例中，第三获得模块S308包括：
确定模块，用于根据车辆实车参数和目标车速，确定目标需求功率PDem；
建表模块，用于将缸内气体压力P、目标需求功率PDem、直流母线上总功率PDC、q 轴前馈电流组成多项式函数，对所述多项式函数进行离散化，获得MAP表，所述缸内气 体压力由所述实时缸压曲线获得；
查表模块，用于查找所述MAP表，得到q轴前馈电流。
直流母线上总功率PDC可通过直流母线上电压和电流获得，对多项式函数进行离散化， 共同组成一张MAP表，通过查找MAP表的方式得到q轴前馈电流，极大的降低了运算 的复杂度，简化了缸压主动减振方法，需要说明的是，该多项式函数即为q轴前馈电流 公式。
在上述实施例的基础上，还包括q轴目标电流确定模块，用于根据q轴前馈电流和q轴原始电流，确定q轴目标电流，所述q轴目标电流是q轴原始电流与q轴前馈电流抵 消后的q轴电流；
其中，通过永磁同步电机控制矢量控制算法得到q轴原始电流，公式如下：
Figure BDA0003329358660000171
式中，uq为q轴电压，iq为q轴原始电流，Lq为q轴电感，Rs为转子电阻，ωr为电 角度，ψm为磁链系数。
q轴原始电流是不向发电机的q轴施加q轴前馈电流时的q轴电流，uq是不向发电机的q轴施加q轴前馈电流时的q轴电压，q轴目标电流是q轴原始电流抵消q轴前馈电流 后实际的q轴上的电流。
在一个可选的实施例中，所述永磁同步电机的d轴前馈电流为0。
电机的交轴也叫q轴，直轴也叫d轴，q轴和d轴是坐标轴而不是实际的轴。在永磁同步电机控制中，为了能够得到类似直流电机的控制特性，在电机转子上建立了一个坐 标系，此坐标系与转子同步转动，取转子磁场方向为d轴，垂直于转子磁场方向为q轴， 将电机的数学模型转换到此坐标系下，可实现d轴和q轴的解耦，从而得到良好控制特 性。在本发明实施例中，默认d轴前馈电流为0。
在上述实施例的基础上，还包括目标转矩确定模块，用于根据发动机往复惯性力和 发电机转子机械转动惯量，确定发电机转子目标转矩公式：
Figure BDA0003329358660000181
式中，J是发电机转子机械转动惯量，Te为发电机转子目标转矩，Bm为含有发动机往复惯性力的机械阻尼系数，ωr为电角度，ωm为机械角度。
其中，Bm是一个综合系数，它里面含有发动机往复惯性力的机械阻尼系数，Bm可通过经验公式和数据标定相结合的方式获得具体数据，Te是由整车控制器(HCU)发送给 电机控制器(MCU)的发电机转子目标转矩数据，该发电机转子目标转矩抵消了发动机 经过燃烧施加在扭矩上的转矩波动，起到了主动减振的作用。
上述系统、装置实施例与系统实施例相对应，与该方法实施例具有同样的技术效果， 具体说明参见方法实施例。装置实施例是基于方法实施例得到的，具体的说明可以参见 方法实施例部分，此处不再赘述。本领域普通技术人员可以理解：附图只是一个实施例的示意图，附图中的模块或流程并不一定是实施本发明所必须的。
本领域普通技术人员可以理解：实施例中的装置中的模块可以按照实施例描述分布 于实施例的装置中，也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然 可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换； 而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和 范围。

Claims (10)

1.一种基于缸压重构的主动减振控制方法，其特征在于，所述方法包括：
根据热力学第一定律推导得到混合动力汽车发动机的初始零维燃烧模型，发动机气缸为闭口系统，缸内气体为理想气体，所述初始零维燃烧模型公式如下：
Figure FDA0003329358650000011
其中，
Figure FDA0003329358650000012
为燃烧过程中曲轴瞬时转角，
Figure FDA0003329358650000013
为燃烧放热率，
Figure FDA0003329358650000014
为排气过程中随排气工质加入到燃烧室的热量变化率，
Figure FDA0003329358650000015
Figure FDA0003329358650000016
为进气过程中随排气工质加入到燃烧室的热量变化率，
Figure FDA0003329358650000017
为通过燃烧室壁散热的散热率，
Figure FDA0003329358650000018
为工质所做机械功率的变化率，
Figure FDA0003329358650000019
P为缸内气体压力，
Figure FDA00033293586500000110
为气缸容积对时间的变化率，
Figure FDA00033293586500000111
为工质内能的变化率，
Figure FDA00033293586500000112
mg为缸内气体质量，且mg为常数，Cv为缸内气体定容比热容，Tgas为缸内气体温度，
Figure FDA00033293586500000113
为缸内气体温度对曲轴转角的变化率；
对理想气体状态方程进行求导，得到所述发动机的第一理想气体状态公式：
PdV+Vdp＝mgRadTgas
式中，P为缸内气体压力，V为缸内容积，mg为缸内气体质量，且mg为常数，Ra为空气气体常数，Tgas为缸内气体温度；
根据所述初始零维燃烧模型、所述第一理想气体状态公式、第一关系式，得出第一零维燃烧模型，所述第一关系式是Cv和Ra之间的关系式，所述Cv和Ra之间的关系式为Cv＝Ra/(k-1)，k为绝热指数，所述第一零维燃烧模型公式如下：
Figure FDA00033293586500000114
其中，由Woschni提出的经验公式得出Qw＝hc(Tgas-Tw)，Tw为壁面温度，hc为热传导系数，
Figure FDA00033293586500000115
ηu为燃烧效率，
Figure FDA00033293586500000116
为燃烧速度，Qf为燃烧放热量，Hu为燃烧的低热值，x为已燃燃油百分比，
Figure FDA00033293586500000117
x由Weibo模型可得，c为燃烧效率系数，m为燃烧品质指数，
Figure FDA00033293586500000118
为燃烧始点对应的曲轴转角，
Figure FDA00033293586500000119
Figure FDA00033293586500000120
为燃烧终点对应的曲轴转角，
Figure FDA00033293586500000121
Figure FDA00033293586500000122
为喷油或点火始点所对应的曲轴转角，τid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差，
Figure FDA00033293586500000123
为燃烧持续期；
根据曲轴上的转速信息和发动机机体信息构建第一往复惯性力矩模型，对所述第一往复惯性力矩模型进行快速FFT变换，获得第二往复惯性力矩模型；
提取所述第二往复惯性力矩模型中的特征值，得到第一特征数据，提取曲轴传递函数中的特征值，得到第二特征数据，分别将所述第一特征数据中的特征值与第二特征数据中与第一特征数据中对应的特征值相减，得到第一修正数据；
利用所述第一修正数据对所述第一零维燃烧模型进行修正，得到第二零维燃烧模型，所述第二零维燃烧模型为参数化的燃烧模型；
根据曲轴转角和所述第二零维燃烧模型，获得发动机的实时缸压曲线；
根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通过查表方式获得，所述q轴前馈电流公式如下：
Figure FDA0003329358650000021
式中，
Figure FDA0003329358650000022
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得，PDem为目标需求功率，PDC为直流母线上总功率。
2.如权利要求1所述的方法，其特征在于，所述根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通过查表方式获得，所述q轴前馈电流公式如下：
Figure FDA0003329358650000023
式中，
Figure FDA0003329358650000024
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得，PDem为目标需求功率，PDC为直流母线上总功率，包括：
根据车辆实车参数和目标车速，确定目标需求功率PDem；
将缸内气体压力P、目标需求功率PDem、直流母线上总功率PDC、q轴前馈电流组成多项式函数，对所述多项式函数进行离散化，获得MAP表，所述缸内气体压力由所述实时缸压曲线获得；
查找所述MAP表，得到q轴前馈电流。
3.如权利要求2所述的方法，其特征在于，还包括根据q轴前馈电流和q轴原始电流，确定q轴目标电流，所述q轴目标电流是q轴原始电流与q轴前馈电流抵消后的q轴电流；
其中，通过永磁同步电机控制矢量控制算法得到q轴原始电流，公式如下：
Figure FDA0003329358650000031
式中，uq为q轴电压，iq为q轴原始电流，Lq为q轴电感，Rs为转子电阻，ωr为电角度，ψm为磁链系数。
4.如权利要求1所述的方法，其特征在于，所述永磁同步电机的d轴前馈电流为0。
5.权利要求1～4中任意一项所述的方法，其特征在于，还包括根据发动机往复惯性力和发电机转子机械转动惯量，确定发电机转子目标转矩公式：
Figure FDA0003329358650000032
式中，J是发电机转子机械转动惯量，Te为发电机转子目标转矩，Bm为含有发动机往复惯性力的机械阻尼系数，ωr为电角度，ωm为机械角度。
6.一种基于缸压重构的主动减振控制装置，其特征在于，所述装置包括：
第一建模模块，用于根据热力学第一定律推导得到混合动力汽车发动机的初始零维燃烧模型，发动机气缸为闭口系统，缸内气体为理想气体，所述初始零维燃烧模型公式如下：
Figure FDA0003329358650000033
其中，
Figure FDA0003329358650000034
为燃烧过程中曲轴瞬时转角，
Figure FDA0003329358650000035
为燃烧放热率，
Figure FDA0003329358650000036
为排气过程中随排气工质加入到燃烧室的热量变化率，
Figure FDA0003329358650000037
Figure FDA0003329358650000038
为进气过程中随排气工质加入到燃烧室的热量变化率，
Figure FDA0003329358650000039
Figure FDA00033293586500000310
为通过燃烧室壁散热的散热率，
Figure FDA00033293586500000311
为工质所做机械功率的变化率，
Figure FDA00033293586500000312
P为缸内气体压力，
Figure FDA00033293586500000313
为气缸容积对时间的变化率，
Figure FDA00033293586500000314
为工质内能的变化率，
Figure FDA00033293586500000315
mg为缸内气体质量，且mg为常数，Cv为缸内气体定容比热容，Tgas为缸内气体温度，
Figure FDA00033293586500000316
为缸内气体温度对曲轴转角的变化率；
第一获得模块，用于对理想气体状态方程进行求导，得到所述发动机的第一理想气体状态公式：
PdV+Vdp＝mgRadTgas
式中，P为缸内气体压力，V为缸内容积，mg为缸内气体质量，且mg为常数，Ra为空气气体常数，Tgas为缸内气体温度；
第二建模模块，用于根据所述初始零维燃烧模型、所述第一理想气体状态公式、第一关系式，得出第一零维燃烧模型，所述第一关系式是Cv和Ra之间的关系式，所述Cv和Ra之间的关系式为Cv＝Ra/(k-1)，k为绝热指数，所述第一零维燃烧模型公式如下：
Figure FDA0003329358650000041
其中，由Woschni提出的经验公式得出Qw＝hc(Tgas-Tw)，Tw为壁面温度，hc为热传导系数，
Figure FDA0003329358650000042
ηu为燃烧效率，
Figure FDA0003329358650000043
为燃烧速度，Qf为燃烧放热量，Hu为燃烧的低热值，x为已燃燃油百分比，
Figure FDA0003329358650000044
x由Weibo模型可得，c为燃烧效率系数，m为燃烧品质指数，
Figure FDA0003329358650000045
为燃烧始点对应的曲轴转角，
Figure FDA0003329358650000046
Figure FDA0003329358650000047
为燃烧终点对应的曲轴转角，
Figure FDA0003329358650000048
Figure FDA0003329358650000049
为喷油或点火始点所对应的曲轴转角，τid是喷射起始点对应曲轴转角和燃烧始点之间曲轴转角差，
Figure FDA00033293586500000410
为燃烧持续期；
第三建模模块，用于根据曲轴上的转速信息和发动机机体信息构建第一往复惯性力矩模型，对所述第一往复惯性力矩模型进行快速FFT变换，获得第二往复惯性力矩模型；
特征提取模块，用于提取所述第二往复惯性力矩模型中的特征值，得到第一特征数据，提取曲轴传递函数中的特征值，得到第二特征数据，分别将所述第一特征数据中的特征值与第二特征数据中与第一特征数据中对应的特征值相减，得到第一修正数据；
模型修正模块，用于利用所述第一修正数据对所述第一零维燃烧模型进行修正，得到第二零维燃烧模型，所述第二零维燃烧模型为参数化的燃烧模型；
第二获得模块，用于根据根据曲轴转角和所述第二零维燃烧模型，获得发动机的实时缸压曲线；
第三获得模块，用于根据所述实时缸压曲线，获得永磁同步电机的q轴前馈电流，所述q轴前馈电流为用于抵消发动机做功时产生振动的电流，所述q轴前馈电流通过查表方式获得，所述q轴前馈电流公式如下：
Figure FDA00033293586500000411
式中，
Figure FDA00033293586500000412
为q轴前馈电流，P是缸内气体压力，P由所述实时缸压曲线获得，PDem为目标需求功率，PDC为直流母线上总功率。
7.如权利要求6所述的装置，其特征在于，所述第三获得模块包括：
确定模块，用于根据车辆实车参数和目标车速，确定目标需求功率PDem；
建表模块，用于将缸内气体压力P、目标需求功率PDem、直流母线上总功率PDC、q轴前馈电流组成多项式函数，对所述多项式函数进行离散化，获得MAP表，所述缸内气体压力由所述实时缸压曲线获得；
查表模块，用于查找所述MAP表，得到q轴前馈电流。
8.如权利要求7所述的装置，其特征在于，还包括q轴目标电流确定模块，用于根据q轴前馈电流和q轴原始电流，确定q轴目标电流，所述q轴目标电流是q轴原始电流与q轴前馈电流抵消后的q轴电流；
其中，通过永磁同步电机控制矢量控制算法得到q轴原始电流，公式如下：
Figure FDA0003329358650000051
式中，uq为q轴电压，iq为q轴原始电流，Lq为q轴电感，Rs为转子电阻，ωr为电角度，ψm为磁链系数。
9.如权利要求6所述的装置，其特征在于，所述永磁同步电机的d轴前馈电流为0。
10.权利要求6～9中任意一项所述的装置，其特征在于，还包括目标转矩确定模块，用于根据发动机往复惯性力和发电机转子机械转动惯量，确定发电机转子目标转矩公式：
Figure FDA0003329358650000052
式中，J是发电机转子机械转动惯量，Te为发电机转子目标转矩，Bm为含有发动机往复惯性力的机械阻尼系数，ωr为电角度，ωm为机械角度。
CN202111273303.7A 2021-10-29 2021-10-29 一种基于缸压重构的主动减振控制方法及装置 Active CN114151216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111273303.7A CN114151216B (zh) 2021-10-29 2021-10-29 一种基于缸压重构的主动减振控制方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111273303.7A CN114151216B (zh) 2021-10-29 2021-10-29 一种基于缸压重构的主动减振控制方法及装置

Publications (2)

Publication Number Publication Date
CN114151216A CN114151216A (zh) 2022-03-08
CN114151216B true CN114151216B (zh) 2022-08-12

Family

ID=80458596

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111273303.7A Active CN114151216B (zh) 2021-10-29 2021-10-29 一种基于缸压重构的主动减振控制方法及装置

Country Status (1)

Country Link
CN (1) CN114151216B (zh)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7463968B2 (en) * 2007-05-03 2008-12-09 Gl Global Technology Operations, Inc. Method and apparatus to control engine stop for a hybrid powertrain system
US20090271057A1 (en) * 2008-04-29 2009-10-29 Ise Corporation Hybrid Vehicle Vibration Reduction System and Method
CN107701320B (zh) * 2017-09-30 2023-07-11 东风商用车有限公司 一种多缸发动机缸内燃烧状态的估计方法及系统
CN109795470B (zh) * 2018-12-13 2020-07-07 清华大学 一种基于发动机瞬时转矩观测器的动力系统主动减振方法
CN109488471B (zh) * 2018-12-20 2021-09-03 南京理工大学 一种自由活塞发动机稳定运行过程控制方法

Also Published As

Publication number Publication date
CN114151216A (zh) 2022-03-08

Similar Documents

Publication Publication Date Title
JP6842285B2 (ja) ハイブリッド車の能動型振動低減制御装置及び方法
US20130233268A1 (en) Engine starting apparatus
JP7431512B2 (ja) 内燃機関制御装置
JP3922277B2 (ja) 内燃機関の空気量推定装置
JP2010023790A (ja) 電動機の制御装置
CA2709022C (en) Method and apparatus for starting an internal combustion engine
JP2013163436A (ja) モータによるエンジン変動抑制装置
CN109641587B (zh) 发动机扭矩平滑
KR100494919B1 (ko) 병렬형 하이브리드 전기자동차의 모터 제어방법 및 그제어장치
CN109795470A (zh) 一种基于发动机瞬时转矩观测器的动力系统主动减振方法
JP7031028B1 (ja) 内燃機関の制御装置及び制御方法
CN110366512A (zh) 混合动力车辆的动力控制方法及动力控制装置
CN110325420A (zh) 混合动力车辆的动力控制方法及动力控制装置
CN110300689A (zh) 混合动力车辆的动力控制方法及动力控制装置
CN114151216B (zh) 一种基于缸压重构的主动减振控制方法及装置
TWI545876B (zh) 電動車延距系統及其充電方法、發電設備與發電設備控制方法
Dönitz et al. Validation of the fuel saving potential of downsized and supercharged hybrid pneumatic engines using vehicle emulation experiments
Liu et al. A study on dynamic torque cancellation in a range extender unit
US10576959B2 (en) Control apparatus for hybrid vehicle
CN113728158B (zh) 内燃机控制装置
JP6489509B2 (ja) ハイブリッド車両の動力制御方法及び動力制御装置
Rusman Charge air configurations for propulsion diesel engines aboard fast naval combatants
CN110466501A (zh) 混合动力汽车发动机启停控制方法
Kataoka et al. Crankshaft positioning utilizing compression force and fast starting with combustion assist for indirect injection engine
CN110290992A (zh) 混合动力车辆的动力控制方法及动力控制装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant