利索能及
我要发布
收藏
专利号: 2023101618087
申请人: 西安理工大学
专利类型:发明专利
专利状态:已下证
更新日期:2026-05-14
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.超高速永磁同步电机转子位置和转速估计方法,其特征在于:具体包括如下步骤:步骤1,建立考虑铁损电阻的超高速永磁同步电机全速域转子位置代价函数;

所述步骤1的具体过程为:

考虑铁损的超高速永磁同步电机在两相静止坐标系的电压方程如下公式(1)所示:2

其中,Ri=(ωr(k)) LdLq/ri是铁损,ωr(k)是第k拍转子电角频率,Ld是d轴电感,Lq是q轴电感,ri是铁损电阻,uα(k)、uβ(k)分别是第k拍定子电压在α轴和β轴的分量,iα(k)、iβ(k)分别是第k拍定子电流在α轴和β轴的分量,Rs是定子电阻,θr(k)是第k拍实际转子位置,L1=0.5(Ld+Lq),L2=0.5(Ld‑Lq),Ts是开关周期,Δiα(k)=iα(k)‑iα(k‑1),Δiβ(k)=iβ(k)‑iβ(k‑1),ψf是永磁体磁链;

为了分离与转子位置无关项,定义中间变量如下公式(2)所示:其中,eα(k)、eβ(k)分别是第k拍不包含转子位置信息的中间变量在α轴和β轴的分量;

由公式(1)和公式(2)可得转子位置相关函数如下公式(3)所示:其中,Fα(θe(k))、Fβ(θe(k))分别是第k拍转子位置相关函数在α轴和β轴的分量,F(θeT(k))=[Fα(θe(k))Fβ(θe(k))];

由公式(3)构建转子位置代价函数如下公式(4)所示:其中,h(θe(k))是转子位置代价函数;

为了保证转子位置的可解性,转子位置代价函数式(4)至少是局部凸函数,在静止及低速区和中高速区转子位置代价函数,简化为如下公式(5)所示:为了在静止和低速增加转子位置代价函数的凸度,在10%额定转速以下注入高频方波电压如下公式(6)所示:

其中,udh、uqh分别是在d轴和q轴注入的高频方波电压,Vhk是注入高频方波电压的振幅;

步骤2,由步骤1得到的全速域转子位置代价函数通过加速自适应凸优化估计初步转子位置;

所述步骤2的具体过程为:

步骤2.1,计算加速自适应凸优化的加速迭代步长;

所述步骤2.1的具体过程为:

根据公式(3)计算自适应因子 如下公式(7)所示:其中, 是自适应因子,λi是自适应系数,i表示第i次迭代;

根据公式(3)和公式(7)计算自适应加速凸优化的迭代步长如下公式(8)所示:其中,li是自适应加速凸优化的迭代步长,R(θe(i))是当前步雅可比矩阵或在上一次迭代中使用的雅可比矩阵,当前步雅可比矩阵I是单位矩阵;

根据公式(3)、公式(7)和公式(8)计算自适应加速凸优化的近似迭代步长如下公式(9)所示:

其中,gi是自适应加速凸优化的近似迭代步长;

根据公式(8)和公式(9)计算自适应加速凸优化的加速迭代步长如下公式(10)所示:ci=li+σigi(10);

其中,ci是自适应加速凸优化的加速迭代步长,σi为可调参数;

步骤2.2,加速自适应凸优化估计初步转子位置;

所述步骤2.2的具体过程为:

为验证当前迭代步的有效性,采用自适应评价指标如下公式(11)所示:其中,ei是自适应评价指标;

根据公式(11)计算的结果是否接受本次迭代计算的加速迭代步长如下公式(12)所示:其中,d0是加速迭代步长取舍指标,θe(i+1)是第i+1迭代得到的转子位置;

根据公式(11)计算的结果计算可调参数σi+1如下公式(13)所示:根据公式(11)计算的结果是否选择更新雅克比矩阵R(θe(i+1))如下公式(14)所示:其中,d1是更新雅克比矩阵的自适应评价指标ei的下限阀值,x是当前雅克比矩阵使用的次数,xmax是同一雅可比矩阵使用的最大次数;

根据公式(11)计算的结果是否选择更新自适应系数λi+1如下公式(15)所示:其中,λmin是自适应系数的最小值,d2是自适应因子时自适应评价指标ei的下限阀值,d3是自适应因子时自适应评价指标ei的上限阀值;

根据公式(11)计算的结果是否选择自适应因子 如下公式(16)所示:步骤3,由步骤2中得到的估计初步转子位置通过锁相环估计超高速永磁同步电机的转子位置和转速。

2.根据权利要求1所述的超高速永磁同步电机转子位置和转速估计方法,其特征在于:所述步骤3的具体过程为:

由加速自适应凸优化估计的初步转子位置θe(k)与锁相环上一拍输出的估计转子位置θer(k‑1)做差计算转子位置误差如下公式(17)所示:Δθ(k)=θe(k)‑θer(k‑1)(17);

其中,θe(k)是加速自适应凸优化估计的第k拍初步转子位置,θer(k‑1)是锁相环第k‑1拍输出的估计转子位置,Δθ(k)是第k拍转子位置误差;

第k拍转子位置误差Δθ(k)通过PI调节器调节得到估计的转速如下公式(18)所示:其中,ωer(k)是估计的第k拍转子转速,Kp是PI的比例增益,Ki是PI的积分增益;

第k拍转子转速ωer(k)通过积分得到第k拍转子位置如下公式(19)所示:θer(k)=θer(k‑1)+ωer(k)Ts(19);

其中,θer(k)是估计的第k拍转子位置。