技术领域
[0001] 本
发明涉及伺服系统控制领域，具体涉及一种基于内模的电动汽车速度跟踪控制方法。
背景技术
[0002] 随着世界石油资源短缺和环境污染日益严重，寻找新形式的
能源代替石油为汽车提供动
力已成为必然的趋势，而电动汽车则是新能源汽车中的主力军。永磁同步
电机因其功率因数高、低损耗等优良性能而被广泛应用在电动汽车中。然而，把永磁同步电机用作电动汽车的驱动装置时，电机或汽车
传动系统参数的变化以及系统外界干扰的引入会直接影响电机的控制性能，进而影响车辆的行驶特性。
[0003] 此外，非线性输出调节理论近年来经历了快速的发展，而作为其常用控制策略的内模控制方法能够解决复杂系统的轨迹跟踪与干扰抑制问题，并具有良好的鲁棒性。针对电动汽车的整体复杂非线性系统，内模控制方法可以实现高
精度的速度跟踪与干扰抑制性能，并允许整体系统的所有参数未知。
发明内容
[0004] 基于背景技术存在的技术问题，本发明提出了一种基于内模的电动汽车速度跟踪控制方法。针对电动汽车的整体复杂系统，设计了内模
控制器，从而实现了高精度的速度跟踪与干扰抑制性能，并具有良好的鲁棒性。
[0005] 本发明的技术方案是：
[0006] 一种基于内模的电动汽车速度跟踪控制方法，其特征在于，包括以下步骤：
[0007] 步骤1：选择永磁同步电机作为
驱动电机并将电机模型与电动汽车传动系统整合，建立整体系统的非线性数学模型；
[0008] 步骤2：将电动汽车的速度跟踪和干扰抑制问题描述为一个全局鲁棒伺服控制问题；
[0009] 步骤3：设计内模，将整体系统的全局鲁棒伺服控制问题转化为由整体系统和内模组成的增广系统的全局鲁棒镇定问题；
[0010] 步骤4：采用反步法设计状态反馈控制器，解决增广系统的全局鲁棒镇定问题。
[0011] 进一步地，所述的一种基于内模的电动汽车速度跟踪控制方法，其特征在于，步骤1中，选择永磁同步电机作为驱动电机并将电机模型与电动汽车传动系统整合，建立整体系统的非线性数学模型，其过程如下：
[0012] 1.1，永磁同步电机数学模型如下：
[0013]
[0014]
[0015]
[0016]
[0017] 其中θM为电机
转子角度，ωM为电机转子
角速度，id,iq,ud,ud为dq轴
定子电流与
电压，JM为电机的
转动惯量，TL为电机的负载转矩，Φv为转子磁链，p为电机极对数，Rs,L为定子
电阻与电感，B为粘性
摩擦系数。
[0018] 1.2，通过分析电动汽车在行驶过程中的受力，可以得到电动汽车传动系统的模型如下：
[0019]
[0020] 其中v为汽车行驶速度，m为汽车的整体
质量，g为重力
加速度，FL为汽车受到的
牵引力，χroad为路面坡度，Fwind为
空气阻力，FR为
滚动阻力。
[0021] 1.3，建立空气阻力Fwind的模型如下：
[0022]
[0023] 其中cair为空气阻力系数，ρa为空气
密度，AL为车辆迎
风正面面积。
[0024] 1.4，建立滚动阻力FR的模型如下：
[0025] FR＝m(cr1+cr2rωW) (4)
[0026] 其中cr1,cr2是取决于轮胎和轮胎压力的常数，r为
车轮半径，ωW为车轮转速。
[0027] 1.5，将车轮转速ωW与汽车行驶速度v关系表示如下：
[0028]
[0029] 1.6，根据
牛顿第二定律得出车轮转速与车轮所受力矩之间的关系如下：
[0030]
[0031] 其中JW为车轮转动惯量，TW为连接车轮
传动轴的输出转矩，Tf为车轮的摩擦转矩。
[0032] 1.7，根据传动部分转矩关系与速度关系得出如下关系式：
[0033] TW＝ntnfTL,ωW＝ntnfωM (7)
[0034] 其中nt,nf分别为变速箱和
主减速器的变速比。
[0035] 1.8，结合公式(1)～(7)，得到永磁同步电机驱动电动汽车的整体系统数学模型如下：
[0036]
[0037]
[0038]
[0039] 其中n＝ntnf；δ＝JW+mr2+n2JM；
[0040] 进一步地，所述的一种基于内模的电动汽车速度跟踪控制方法，其特征在于，步骤2中，将电动汽车的速度跟踪和干扰抑制问题描述为一个全局鲁棒伺服控制问题，其过程如下：
[0041] 2.1，假定车轮的参考转速ωd和包含摩擦转矩的等效干扰βcr1+βgsinχroad+Tf可由如下的外部系统产生：
[0042]
[0043] 其中A1,G1,G2为定常矩阵。
[0044] 2.2，令：x1,1＝ωW,x1,2＝iq,x2,1＝id,u1＝uq,u2＝ud,a1 6＝ n p , a 22 ＝n p ,
将系统(8)写为如下形式：
[0045]
[0046]
[0047]
[0048]
[0049] 其中e1为车轮转速跟踪误差，e2为d轴电流跟踪误差。
[0050] 2 .3，考虑不确定因素所产生的系统参数摄动，定义不确定参数其中 为整体系统的标称值，w∈R13。将系统(9)与(10)相结合，得到如下的紧凑形式：
[0051]
[0052] v＝A1v,
[0053] e＝H(x,u,v,w). (11)
[0054] 其中x＝(x1,1,x1,2,x2,1)T,u＝(u1,u2)T，
[0055]
[0056] 2.4，此时系统(8)的全局速度跟踪控制问题已被描述为系统(11)的全局鲁棒伺服控制问题，其控制目标为在保证从任意初始值出发的
闭环系统轨迹有界的情况下，跟踪误差渐近趋于0。
[0057] 进一步地，所述的一种基于内模的电动汽车速度跟踪控制方法，其特征在于，步骤3中，设计内模，将整体系统的全局鲁棒伺服控制问题转化为由整体系统和内模组成的增广系统的全局鲁棒镇定问题，其过程如下：
[0058] 3.1，求解如下的调节器方程：
[0059]
[0060] 0＝H(x(v,w),u(v,w),v,w). (12)
[0061] 其中x(v,w),u(v,w)分别是稳态状态和稳态输入。得到状态和输入的稳态解如下：
[0062] x1,1(v,w)＝G1v,
[0063]
[0064] x2,1(v,w)＝0,
[0065]
[0066]
[0067] 3.2，令g(x,u)＝col(x1,2,u1,u2)，用gi(x,u)表示g(x,u)中第i个元素，其中i＝1,2,3。构建如下的稳态发生器来产生稳态解：
[0068]
[0069]
[0070] 其中Ti为任意非奇异矩阵，(Φi,Ψi)为一对能
观测矩阵。
[0071] 3.3，选择一对能控矩阵(Mi,Ni)，其中Mi为Hurwitz矩阵，使得Ti满足如下的Sylvester方程：
[0072] TiΦi-MiTi＝NiΨi. (14)
[0073] 3.4，设计内模为如下形式：
[0074]
[0075]
[0076]
[0077] 3.5，进行如下的坐标变换与输入变换：
[0078]
[0079]
[0080]
[0081]
[0082]
[0083]
[0084]
[0085]
[0086] 得到如下的误差方程：
[0087]
[0088]
[0089]
[0090]
[0091]
[0092]
[0093] 其中
[0094] b1＝Ψ1T1-1(M1+N1Ψ1T1-1),
[0095] b2＝Ψ1T1-1N1,
[0096]
[0097] c2(v)＝-a11-2a12G1v+b2,
[0098]
[0099]
[0100]
[0101] c6(v)＝-a16G1v,
[0102]
[0103]
[0104]
[0105] c10(v)＝a22G1vΨ1T1-1,
[0106]
[0107] c12(v)＝a22G1v,
[0108]
[0109] d2＝b11Ψ1T1-1,
[0110] d3＝-a12,
[0111]
[0112]
[0113]
[0114]
[0115]
[0116] d9＝-a14Ψ1T1-1-b1,
[0117]
[0118]
[0119] d12＝-a16,
[0120]
[0121]
[0122]
[0123]
[0124] d17＝b21Ψ3T3-1,
[0125] d18＝Ψ3T3-1N3-a21,
[0126] d19＝a22,
[0127] d20＝a22Ψ1T1-1,
[0128]
[0129] 3.6，此时，系统(11)的全局鲁棒伺服控制问题已经被转化为系统(17)的全局鲁棒镇定问题。
[0130] 进一步地，所述的一种基于内模的电动汽车速度跟踪控制方法，其特征在于，步骤4中，采用反步法设计状态反馈控制器，解决增广系统的全局鲁棒镇定问题，其过程如下：
[0131] 4.1，为了使用反步法设计控制器，首先定义如下的记号：
[0132]
[0133]
[0134]
[0135]
[0136]
[0137] 其中 为特定的非负光滑函数。
[0138] 4.2，令 针对X1子系统，令 其中P1为一个正定对称矩阵满足
I为实对称矩阵，m1,m2为特定的正数。选择 于是存在
一个足够大的增益k1满足如下的不等式：
[0139]
[0140] 其中l1为特定的正数。
[0141] 4.3，针对X2子系统，令 其中P2为一个正定对称矩阵满足 m3,m4为特定的正数。选择 于是存在
一个足够大的增益k2满足如下的不等式：
[0142]
[0143] 其中l2为特定的正数。
[0144] 4.4，最后，令 其中P3为一个正定对称矩阵满足 m5,m6为特定的正数。选择 于是存在一个足
够大的增益k3满足如下的不等式：
[0145]
[0146] 4.5，得到如下控制律解决系统(17)的全局镇定问题：
[0147]
[0148]
[0149]
[0150] 4.6，得到最终的控制器为如下形式：
[0151]
[0152]
[0153]
[0154]
[0155]
[0156]
[0157] 本发明的优点是：
[0158] 本发明提出的一种基于内模的电动汽车速度跟踪控制方法，针对电动汽车行驶过程中运行环境造成系统参数变化和系统外界干扰引入的现象，设计了基于内模的状态反馈控制器，实现了电动汽车速度跟踪控制，具有高精度的速度跟踪与干扰抑制性能，并具有良好的鲁棒性。
附图说明
[0159] 图1为电动汽车基本受力分析图；
[0161] 图3为车轮转速跟踪曲线；
[0162] 图4为车轮转速误差曲线；
[0163] 图5为id电流曲线；
[0164] 图6为iq电流曲线；
[0165] 图7为定子电压ud曲线；
[0166] 图8为定子电压uq曲线。
具体实施方式
[0167] 下面将结合本发明
实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
[0168] 实施例。
[0169] 如图1、2所示，一种基于内模的电动汽车速度跟踪控制方法，包含以下步骤：
[0170] 步骤1，选择永磁同步电机作为驱动电机并将电机模型与电动汽车传动系统整合，建立整体系统的非线性数学模型，其过程如下：
[0171] 1.1，永磁同步电机数学模型如下：
[0172]
[0173]
[0174]
[0175]
[0176] 其中θM为电机转子角度，ωM为电机转子角速度，id,iq,ud,ud为dq轴定子电流与电压，JM为电机的转动惯量，TL为电机的负载转矩，Φv为转子磁链，p为电机极对数，Rs,L为定子电阻与电感，B为粘性摩擦系数。
[0177] 1.2，通过分析电动汽车在行驶过程中的受力，可以得到电动汽车传动系统的模型如下：
[0178]
[0179] 其中v为汽车行驶速度，m为汽车的整体质量，g为
重力加速度，FL为汽车受到的牵引力，χroad为路面坡度，Fwind为空气阻力，FR为滚动阻力。
[0180] 1.3，建立空气阻力Fwind的模型如下：
[0181]
[0182] 其中cair为空气阻力系数，ρa为空气密度，AL为车辆迎风正面面积。
[0183] 1.4，建立滚动阻力FR的模型如下：
[0184] FR＝m(cr1+cr2rωW) (4)
[0185] 其中cr1,cr2是取决于轮胎和轮胎压力的常数，r为车轮半径，ωW为车轮转速。
[0186] 1.5，将车轮转速ωW与汽车行驶速度v关系表示如下：
[0187]
[0188] 1.6，根据牛顿第二定律得出车轮转速与车轮所受力矩之间的关系如下：
[0189]
[0190] 其中JW为车轮转动惯量，TW为连接车轮传动轴的输出转矩，Tf为车轮的摩擦转矩。
[0191] 1.7，根据传动部分转矩关系与速度关系得出如下关系式：
[0192] TW＝ntnfTL,ωW＝ntnfωM (7)
[0193] 其中nt,nf分别为变速箱和主减速器的变速比。
[0194] 1.8，结合公式(1)～(7)，得到永磁同步电机驱动电动汽车的整体系统数学模型如下：
[0195]
[0196]
[0197]
[0198] 其中n＝ntnf；δ＝JW+mr2+n2JM；β＝rm；
[0199] 步骤2，将电动汽车的速度跟踪和干扰抑制问题描述为一个全局鲁棒伺服控制问题，其过程如下：
[0200] 2.1，假定车轮的参考转速ωd和包含摩擦转矩的等效干扰βcr1+βgsinχroad+Tf可由如下的外部系统产生：
[0201]
[0202] 其中A1,G1,G2为定常矩阵。
[0203] 2.2，令：x1,1＝ωW,x1,2＝iq,x2,1＝id,u1＝uq,u2＝ud,a1 6＝ n p , a 22 ＝n p ,
将系统(8)写为如下形式：
[0204]
[0205]
[0206]
[0207]
[0208] 其中e1为车轮转速跟踪误差，e2为d轴电流跟踪误差。
[0209] 2 .3，考虑不确定因素所产生的系统参数摄动，定义不确定参数其中 为整体系统的标称值，w∈R13。将系统(9)与(10)相结合，得到如下的紧凑形式：
[0210]
[0211] v＝A1v,
[0212] e＝H(x,u,v,w). (11)
[0213] 其中x＝(x1,1,x1,2,x2,1)T,u＝(u1,u2)T，
[0214]
[0215] 2.4，此时系统(8)的全局速度跟踪控制问题已被描述为系统(11)的全局鲁棒伺服控制问题，其控制目标为在保证从任意初始值出发的闭环系统轨迹有界的情况下，跟踪误差渐近趋于0。
[0216] 步骤3，设计内模，将整体系统的全局鲁棒伺服控制问题转化为由整体系统和内模组成的增广系统的全局鲁棒镇定问题，其过程如下：
[0217] 3.1，求解如下的调节器方程：
[0218]
[0219] 0＝H(x(v,w),u(v,w),v,w). (12)
[0220] 其中x(v,w),u(v,w)分别是稳态状态和稳态输入。得到状态和输入的稳态解如下：
[0221] x1,1(v,w)＝G1v,
[0222]
[0223] x2,1(v,w)＝0,
[0224]
[0225]
[0226] 3.2，令g(x,u)＝col(x1,2,u1,u2)，用gi(x,u)表示g(x,u)中第i个元素，其中i＝1,2,3。构建如下的稳态发生器来产生稳态解：
[0227]
[0228]
[0229] 其中Ti为任意非奇异矩阵，(Φi,Ψi)为一对能观测矩阵。
[0230] 3.3，选择一对能控矩阵(Mi,Ni)，其中Mi为Hurwitz矩阵，使得Ti满足如下的Sylvester方程：
[0231] TiΦi-MiTi＝NiΨi. (14)
[0232] 3.4，设计内模为如下形式：
[0233]
[0234]
[0235]
[0236] 3.5，进行如下的坐标变换与输入变换：
[0237]
[0238]
[0239]
[0240]
[0241]
[0242]
[0243]
[0244]
[0245] 得到如下的误差方程：
[0246]
[0247]
[0248]
[0249]
[0250]
[0251]
[0252] 其中
[0253] b1＝Ψ1T1-1(M1+N1Ψ1T1-1),
[0254] b2＝Ψ1T1-1N1,
[0255]
[0256]
[0257]
[0258]
[0259]
[0260] c6(v)＝-a16G1v,
[0261]
[0262]
[0263]
[0264] c10(v)＝a22G1vΨ1T1-1,
[0265]
[0266] c12(v)＝a22G1v,
[0267]
[0268] d2＝b11Ψ1T1-1,
[0269] d3＝-a12,
[0270]
[0271]
[0272]
[0273]
[0274]
[0275] d9＝-a14Ψ1T1-1-b1,
[0276]
[0277]
[0278] d12＝-a16,
[0279]
[0280]
[0281]
[0282]
[0283] d17＝b21Ψ3T3-1,
[0284] d18＝Ψ3T3-1N3-a21,
[0285] d19＝a22,
[0286] d20＝a22Ψ1T1-1,
[0287]
[0288] 3.6，此时，系统(11)的全局鲁棒伺服控制问题已经被转化为系统(17)的全局鲁棒镇定问题。
[0289] 步骤4，采用反步法设计状态反馈控制器，解决增广系统的全局鲁棒镇定问题，其过程如下：
[0290] 4.1，为了使用反步法设计控制器，首先定义如下的记号：
[0291]
[0292]
[0293]
[0294]
[0295]
[0296] 其中 为特定的非负光滑函数。
[0297] 4.2，令 针对X1子系统，令 其中P1为一个正定对称矩阵满足
I为实对称矩阵，m1,m2为特定的正数。选择 于是存在
一个足够大的增益k1满足如下的不等式：
[0298]
[0299] 其中l1为特定的正数。
[0300] 4.3，针对X2子系统，令 其中P2为一个正定对称矩阵满足 m3,m4为特定的正数。选择 于是存在
一个足够大的增益k2满足如下的不等式：
[0301]
[0302] 其中l2为特定的正数。
[0303] 4.4，最后，令 其中P3为一个正定对称矩阵满足 m5,m6为特定的正数。选择 于是存在一个足
够大的增益k3满足如下的不等式：
[0304]
[0305] 4.5，得到如下控制律解决系统(17)的全局镇定问题：
[0306]
[0307]
[0308]
[0309] 4.6，得到最终的控制器为如下形式：
[0310]
[0311]
[0312]
[0313]
[0314]
[0315]
[0316] 为了验证所提方法的有效性，本发明对所提出控制器的控制效果进行仿真验证，所选用的永磁同步电机的标称值为：极对数p＝4，定子电阻 磁链定子电感 转动惯量 粘性摩擦系数
所选用的电动汽车传动系统的标称值为：车体总质量 车
轮转动惯量 迎风面积 空气阻力系数 空气密度
总变速比n＝7，轮胎常数 轮胎压力常数 重力
2
加速度g＝9.8m/s，路面坡度
[0317] 摩擦转矩设置为一个常值
信号Tf＝10N.m，参考速度设置为一个斜坡信号与常值信号的组合，当0
[0318] G1＝[1 0 0],G2＝[0 0 1].
[0319] 控制器参数设计为：
[0320]
[0321]
[0322] Ψ1T1-1＝[1 2 10000], k1＝800,k2＝10,k3＝300。
[0323] 系 统 不 确 定 参 数 选 择 如 下 ：设置初值为ωW(0)＝0rad/s,id(0)＝0.1A,
iq(0)＝0A,η(0)＝0。
[0324] 基于上述的一系列参数，图3为车轮转速跟踪曲线，图4为车轮转速误差曲线，图5为定子电流id曲线，图中曲线反映了所设计的控制器在外部干扰和整体系统参数摄动的条件下有良好的速度跟踪性能和电流跟踪性能，图6为定子电流iq曲线，图7和图8分别为永磁同步电机ud和uq曲线。从图3-图8的仿真结果可以看出，所设计的内模控制器可以实现高精度的速度跟踪与干扰抑制性能，并具有良好的鲁棒性。
[0325] 以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉
本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
