【永磁同步电机(PMSM)】 4. 坐标变换的 Matlab 仿真

【永磁同步电机(PMSM)】 4. 坐标变换的 Matlab 仿真

    • 1. Clarke 变换的模型与仿真
      • 1.1 Clarke 变换
      • 1.2 Clarke 变换的仿真模型
    • 2. Park 变换的模型与仿真
      • 2.1 Park 变换
      • 2.2 Park 变换的仿真模型
    • 3. Simscape标准库变换模块
      • 3.1 abc to Alpha-Beta-Zero 模块
      • 3.2 abc to dq0 模块
    • 3. 基于 S函数的仿真模型


由于电机处于高速旋转运动状态,在 abc 三相静止坐标系(自然坐标系)所建立的 PMSM 数学模型,存在复杂的变量耦合,难以分析和求解。

通过数学变换将 abc 三相静止坐标系转换到 αβ两相静止坐标系 和 dq 同步旋转坐标系,可以对模型进行解耦和简化。本节讨论坐标系变换的模型与仿真。


1. Clarke 变换的模型与仿真

1.1 Clarke 变换

通过Clarke变换(有些文献也写作 Clark 变换)将 abc 三相静止坐标系转换到 αβ两相静止坐标系,记为 T 3 s / 2 s T_{3s/2s} T3s/2s

两相静止坐标系中的信号相互正交,更容易分离基波分量和谐波分量。对于 PMSM,假设虚拟电机有两个虚拟绕组,呈90°轴距分布,在定子槽中呈正弦分布,这两个绕组分别为直轴绕组和交轴绕组,即α,β轴。为简化分析过程,可以在建立αβ静止坐标系时将α轴与abc静止坐标系的a相轴重合。

Clarke 变换公式如下。通过 Clarke 变换,将abc三相静止坐标系下的 fa,fb,fc 分别映射到α-β轴上,得到αβ两相静止坐标系下的虚拟物理量 fα, fβ。

[ f α , f β , f 0 ] T = T 3 s / 2 s [ f a , f b , f c ] T [f_{\alpha}, f_{\beta}, f_0]^T = T_{3s/2s} [f_a, f_b,f_c]^T [fα,fβ,f0]T=T3s/2s[fa,fb,fc]T
T 3 s / 2 s = 2 3 [ 1 − 1 2 − 1 2 0 3 2 − 3 2 2 2 2 2 2 2 ] T_{3s/2s} =\frac{2}{3} \begin{bmatrix} 1 & -\frac{1}{2} & -\frac{1}{2}\\ 0 & \frac{\sqrt{3}}{2} & -\frac{\sqrt{3}}{2}\\ \frac{\sqrt{2}}{2} & \frac{\sqrt{2}}{2} & \frac{\sqrt{2}}{2} \end{bmatrix} T3s/2s=32 1022 2123 22 2123 22

f f f 表示电压 u u u,电流 i i i,磁链 ψ \psi ψ等时变变量,例如:

[ i α i β i o ] = T 3 s / 2 s [ i a i b i c ] = 2 3 [ 1 − 1 / 2 − 1 / 2 0 3 / 2 − 3 / 2 2 / 2 2 / 2 2 / 2 ] [ i a i b i c ] \begin{bmatrix}i_{\alpha} \\i_{\beta}\\i_o\end{bmatrix} = T_{3s/2s} \begin{bmatrix}i_a \\i_b \\ i_c\end{bmatrix}= \frac{2}{3} \begin{bmatrix} 1 & -1/2 & -1/2\\ 0 & \sqrt{3}/2 & -\sqrt{3}/2\\ \sqrt{2}/2 & \sqrt{2}/2 & \sqrt{2}/2 \end{bmatrix} \begin{bmatrix}i_a \\i_b \\ i_c\end{bmatrix} iαiβio =T3s/2s iaibic =32 102 /21/23 /22 /21/23 /22 /2 iaibic

式中, i 0 i_0 i0 为零序分量。

Clarke变换的逆变换,称为反Clarke变换(用T2s/3s表示)。

[ u a u b u c ] = T 2 s / 3 s [ u α u β u o ] = 2 3 [ 1 0 2 / 2 − 1 / 2 3 / 2 2 / 2 − 1 / 2 − 3 / 2 2 / 2 ] [ u α u β u o ] \begin{bmatrix}u_a \\u_b\\u_c\end{bmatrix} = T_{2s/3s} \begin{bmatrix}u_{\alpha}\\u_{\beta}\\u_o\end{bmatrix} =\frac{2}{3} \begin{bmatrix} 1 & 0 & \sqrt{2}/2\\ -1/2 & \sqrt{3} /2 & \sqrt{2}/2\\ -1/2 & -\sqrt{3} /2 & \sqrt{2}/2 \end{bmatrix} \begin{bmatrix}u_{\alpha} \\u_{\beta} \\u_o\end{bmatrix} uaubuc =T2s/3s uαuβuo =32 11/21/203 /23 /22 /22 /22 /2 uαuβuo


1.2 Clarke 变换的仿真模型

根据 Clarke 变换的数学公式,使用Matlab/Simulink建立 Clarke 变换的仿真模型。
(1)新建模型:打开Matlab软件,在Simulink模型编辑界面中新建“空白模型”。
(2)添加模块:打开库浏览器,从Simulink、Simscape\Electrical\Specialized Power Systems等标准库中依次选取模块。按照设计计算结果设置模块参数。用户也可以直接在模型编辑界面中,双击鼠标左键调出“搜索模块”弹窗,输入模块名称,直接选择和添加模块。
(3)搭建模型:按照电路原理图连接各模块,搭建Buck变换电路的仿真模型。
(4)信号监测:使用信号标记模块goto、信号分解模块Demux、总线选择模块Bus Select,提取和选择需要观测的信号,作为示波器的输入信号。
(5)模型设置:选择“模型配置参数”,在求解器中选择仿真算法ode23tb(stiff/TR-BDF2),仿真时间为0.02s。
按照以上步骤,建立 Clarke 变换和反变换的仿真模型(Clarke01.slx),如下图所示。

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述


仿真结果如下图所示。


2. Park 变换的模型与仿真

2.1 Park 变换

通过 Park 变换将 αβ两相静止坐标系 转换到 dq 同步旋转坐标系,记为 T 2 s / 2 r T_{2s/2r} T2s/2r

Park 变换公式如下。通过 Park 变换,将αβ两相静止坐标系fα, fβ 分别映射到d-q轴上,得到 dq 同步旋转坐标系下的虚拟物理量 fd, fq。

[ f d , f q ] T = T 2 s / 2 r [ f α , f β ] T [f_d, f_q]^T = T_{2s/2r} [f_{\alpha}, f_{\beta}]^T [fd,fq]T=T2s/2r[fα,fβ]T
T 2 s / 2 r = [ c o s ( θ ) s i n ( θ ) − s i n ( θ ) c o s ( θ ) ] T_{2s/2r} = \begin{bmatrix} cos(\theta)&sin(\theta)\\ -sin(\theta)&cos(\theta)\\ \end{bmatrix} T2s/2r=[cos(θ)sin(θ)sin(θ)cos(θ)]

f f f 表示电压 u u u,电流 i i i,磁链 ψ \psi ψ等时变变量,例如:

[ u d u q ] = T 2 s / 2 r [ u α u β ] = [ c o s ( θ ) s i n ( θ ) − s i n ( θ ) c o s ( θ ) ] [ u α u β ] \begin{bmatrix}u_d\\u_q\end{bmatrix} = T_{2s/2r} \begin{bmatrix}u_{\alpha} \\u_{\beta}\end{bmatrix} = \begin{bmatrix} cos(\theta)&sin(\theta)\\ -sin(\theta)&cos(\theta)\\ \end{bmatrix} \begin{bmatrix}u_{\alpha}\\u_{\beta}\end{bmatrix} [uduq]=T2s/2r[uαuβ]=[cos(θ)sin(θ)sin(θ)cos(θ)][uαuβ]

式中, i 0 i_0 i0 为零序分量。

Park 变换的逆变换,称为反Park 变换(用T2r/2s表示)。

[ u α u β ] = T 2 r / 2 s [ u d u q ] = [ c o s ( θ ) − s i n ( θ ) s i n ( θ ) c o s ( θ ) ] [ u d u q ] \begin{bmatrix}u_{\alpha}\\u_{\beta}\end{bmatrix} = T_{2r/2s} \begin{bmatrix}u_d\\u_q\end{bmatrix} = \begin{bmatrix} cos(\theta) & -sin(\theta)\\ sin(\theta) & cos(\theta) \end{bmatrix} \begin{bmatrix}u_d\\u_q\end{bmatrix} [uαuβ]=T2r/2s[uduq]=[cos(θ)sin(θ)sin(θ)cos(θ)][uduq]


2.2 Park 变换的仿真模型

根据 Clarke 变换的数学公式,使用Matlab/Simulink建立 Clarke 变换的仿真模型。

在这里插入图片描述

在这里插入图片描述


仿真结果如下图所示。

在这里插入图片描述

由仿真结果可见,经过坐标变换,在dq 同步旋转坐标系中的 Vd, Vq 成为解耦的直流量,便于进行控制。


3. Simscape标准库变换模块

Simscape标准库提供了 abc to Alpha-Beta-Zero模块和 abc to dq0 模块。

3.1 abc to Alpha-Beta-Zero 模块

abc to Alpha-Beta-Zero 模块对三相abc信号进行 Clarke 变换。Alpha-Beta-Zero to abc 模块对 αβ0 分量执行 Clarke 逆变换。

在这里插入图片描述

假设ua、ub、uc量表示三个正弦平衡电压:

[ u a u b u c ] = U [ c o s ( ω t ) c o s ( ω t − 2 π / 3 ) c o s ( ω t + 2 π / 3 ) ] \begin{bmatrix}u_a\\u_b\\u_c\end{bmatrix} = U \begin{bmatrix}cos{(\omega t)}\\cos{(\omega t-2\pi/3)}\\cos{(\omega t+2\pi/3)}\end{bmatrix} uaubuc =U cos(ωt)cos(ωt2π/3)cos(ωt+2π/3)

则 iα和iβ分量表示旋转空间矢量Is在α轴与相位a轴对齐的固定参考系中的坐标。振幅与三个电流产生的旋转磁动势成正比。


3.2 abc to dq0 模块

abc to dq0 模块 使用Park变换将三相(abc)信号变换为dq0旋转参考系。旋转框架的角位置由输入 wt 给出,单位为rad。
dq0 to abc 模块 使用逆Park变换将dq0旋转参考系变换为三相(abc)信号。旋转框架的角位置由输入 wt 给出,单位为rad。
当wt=0处的旋转帧对齐在相位A轴后90度时,Mag=1且phase=0度的正序信号产生以下dq值:d=1,q=0。

在这里插入图片描述

该模块支持用于 Park 变换:

  • 当旋转框架在t=0时与相位A轴对齐时,即t=0时,d轴与A轴对齐。这种类型的Park变换也称为余弦型Park变换。
  • 当旋转框架在相位A轴后对齐90度时,即t=0时,q轴与A轴对齐。这种类型的公园改造也被称为正弦型Park 变换。在具有三相同步和异步电机的Simscape™Electrical™专用电力系统模型中使用此转换。

通过在三相静止参考系中执行abc到αβ0-Clarke变换,从abc信号中推导出dq0分量。然后在旋转参考系中执行αβ0到dq0的变换,即通过对空间向量Us=uα+j·uβ执行−ωt旋转。

abc-to-dq0变换取决于t=0时的dq帧对齐。旋转坐标系的位置由ω.t给出,其中ω表示dq坐标系的旋转速度。
在这里插入图片描述


3. 基于 S函数的仿真模型

Matlab/Simulink在Simscape标准库中提供了多种坐标转换模块,使用abc_to_dq0模块可以完成从abc三相静止坐标系到dq旋转参考系的转换。

基于abc_to_dq0模块直接建立的仿真模型更加简洁。对于复杂算法也可以使用Simulink中的S函数建立仿真模型,如下图所示。

在这里插入图片描述


参考文献:袁雷等,现代永磁同步电机控制原理及MATLAB仿真,北京航空航天大学出版社,2016

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/429160.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

更换硬盘后,电脑装完系统进不去?或PE能识别硬盘但开机/启动/BIOS识别不了硬盘解决办法

由于现在的电脑主板,默认都是UEFI启动,硬盘只有使用GUID分区表,主板BIOS才找得到系统引导! 而当我们拿到一块新硬盘,使用分区工具默认类型分区,默认是MBR类型,所以这种分区的硬盘,B…

14.面试算法-字符串常见算法题(三)

1. 字符串回文问题 1.1 LeetCode.125. 验证回文串 回文问题在链表中是重点,在字符串中同样是个重点。当初我去美团面试第一轮技术面的第一个算法题就是让写判断字符串回文的问题。 这个本身还是比较简单的,只要先转换成字符数组,然后使用双…

肺结节检测系统源码分享

肺结节检测检测系统源码分享 [一条龙教学YOLOV8标注好的数据集一键训练_70全套改进创新点发刊_Web前端展示] 1.研究背景与意义 项目参考AAAI Association for the Advancement of Artificial Intelligence 项目来源AACV Association for the Advancement of Computer Visio…

uniapp uview扩展u-picker支持日历期间 年期间 月期间 时分期间组件

uniapp uview扩展u-picker支持日历期间 年期间 月期间 时分期间组件 日历期间、年期间、月期间及时分期间组件在不同的应用场景中发挥着重要的作用。这些组件通常用于表单、应用程序或网站中,以方便用户输入和选择特定的日期和时间范围。以下是这些组件的主要作用&a…

C++:日期类的实现

目录 一、前言 二、头文件 三、各个函数的实现 打印、检查日期及获取日期 、、-、-、 、<、<、>、>、 &#xff01; 日期-日期 >>、<< 一、前言 前面几篇讲了关于类和对象的一些知识&#xff0c;本篇就来实现一下前面用到的日期类。 二、头文…

TryHackMe 第3天 | Pre Security (中)

该学习路径讲解了网络安全入门的必备技术知识&#xff0c;比如计算机网络、网络协议、Linux命令、Windows设置等内容。上一篇中简短介绍了计算机网络相关的知识&#xff0c;本篇博客将记录 网络协议 部分。 How the web works? DNS in detail DNS (Domain name system&…

Java面试篇-AOP专题(什么是AOP、AOP的几个核心概念、AOP的常用场景、使用AOP记录操作日志、Spring中的事务是如何实现的)

文章目录 1. 什么是AOP2. AOP的几个核心概念3. AOP的常用场景4. 使用AOP记录操作日志4.1 准备工作4.1.1 引入Maven依赖4.1.2 UserController.java4.1.3 User.java4.1.4 UserService.java 4.2 具体实现&#xff08;以根据id查询用户信息为例&#xff09;4.2.1 定义切面类&#x…

SkyWalking 环境搭建部署

架构简介 skywalking agent : 和业务系统绑定在一起,负责收集各种监控数据skywalking oapservice : 是负责处理监控数据的,比如接受skywalking agent的监控数据,并存储在数据库中;接受skywalking webapp的前端请求,从数据库查询数据,并返回数据给前端。Skywalking oapserv…

华为HarmonyOS地图服务 7- 在地图上绘制标记

场景介绍 本章节将向您介绍如何在地图的指定位置添加标记以标识位置、商家、建筑等。 点标记用来在地图上标记任何位置,例如用户位置、车辆位置、店铺位置等一切带有位置属性的事物。Map Kit提供的点标记功能(又称 Marker)封装了大量的触发事件,例如点击事件、长按事件、…

【RabbitMQ】消息分发、事务

消息分发 概念 RabbitMQ队列拥有多个消费者时&#xff0c;队列会把收到的消息分派给不同的消费者。每条消息只会发送给订阅该队列订阅列表里的一个消费者。这种方式非常适合扩展&#xff0c;如果现在负载加重&#xff0c;那么只需要创建更多的消费者来消费处理消息即可。 默…

linux网络编程5

24.9.21学习目录 一.TCP1.TCP流程2.TCP相关函数3.三次握手 一.TCP 1.TCP流程 服务器流程&#xff1a; 创建套接字socket&#xff08;&#xff09;将套接字与服务器网络信息结构体绑定bind&#xff08;&#xff09;将套接字设置为监听状态listen&#xff08;&#xff09;阻塞等…

进程间的通信4 共享内存

共享内存 1.共享内存简介 共享内存是将分配的物理空间直接映射到进程的用户虚拟地址空间中&#xff0c;减少数据在内核空间缓存共享内存是一种效率较高的进程间通讯的方式在 Linux 系统中通过 ipcs -m 查看所有的共享内存 共享内存模型图 2.共享内存的创建 1.函数头文件 #…

Java算法专栏

专栏导读 在当今这个技术日新月异的时代&#xff0c;Java算法作为软件开发的核心&#xff0c;对于提升程序性能和解决复杂问题至关重要。本“Java算法”专栏旨在帮助读者深入理解Java编程语言中的算法原理和应用&#xff0c;通过实战案例和深入分析&#xff0c;使读者能够掌握…

Java汽车销售管理

技术架构&#xff1a; springboot mybatis Mysql5.7 vue2 npm node 功能描述&#xff1a; 针对汽车销售提供客户信息、车辆信息、订单信息、销售人员管理、财务报表等功能&#xff0c;提供经理和销售两种角色进行管理 效果图&#xff1a;

Python基础学习(3)

目录 一&#xff0c;函数 1&#xff0c;函数的定义 2&#xff0c;函数的参数 1&#xff0c;默认值 2&#xff0c;传参 3&#xff0c;返回值 4&#xff0c;变量的作用域 5&#xff0c;函数的调用 二&#xff0c;常用数据结构 1&#xff0c;列表 列表的定义 列表的特性…

【Geoserver使用】REST API调用(工作空间部分)

文章目录 前言一、Geoserver REST API(GeoServer Workspace)二、GeoServer Workspace接口使用1.GET请求 /workspaces2.POST请求 /workspaces3.GET请求 /workspaces/{workspaceName}4.PUT /workspaces/{workspaceName}5.DELETE /workspaces/{workspaceName} 总结 前言 根据Geos…

C++ | Leetcode C++题解之第423题从英文中重建数字

题目&#xff1a; 题解&#xff1a; class Solution { public:string originalDigits(string s) {unordered_map<char, int> c;for (char ch: s) {c[ch];}vector<int> cnt(10);cnt[0] c[z];cnt[2] c[w];cnt[4] c[u];cnt[6] c[x];cnt[8] c[g];cnt[3] c[h] - …

YOLOv10 简介

YOLOv10&#xff0c;由清华大学的研究人员基于 Ultralytics Python 包构建&#xff0c;引入了一种全新的实时目标检测方法&#xff0c;该方法解决了以往 YOLO 版本中后处理和模型架构方面的不足。通过消除非极大值抑制&#xff08;NMS&#xff09;并优化各种模型组件&#xff0…

【解决】chrome 谷歌浏览器,鼠标点击任何区域都是 Input 输入框的状态,能看到输入的光标

chrome 谷歌浏览器&#xff0c;鼠标点击任何区域都是 Input 输入框的状态&#xff0c;能看到输入的光标 今天打开电脑的时候&#xff0c;网页中任何文本的地方&#xff0c;只要鼠标点击&#xff0c;就会出现一个输入的光标&#xff0c;无论在哪个站点哪个页面都是如此。 我知道…

十四、运算放大电路

运算放大电路 1、理想运算放大器的概念。运放的输入端虚拟短路、虚拟断路之间的区别; 2、反相输入方式的运放电路的主要用途&#xff0c;以及输入电压与输出电压信号的相位 3、同相输入方式下的增益表达式(输入阻抗、输出阻抗)