找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 33332|回复: 12
打印 上一主题 下一主题
收起左侧

基于MATLAB的异步电机变频调速系统的设计文档(Simulink仿真模型)

  [复制链接]
跳转到指定楼层
楼主
基于MATLAB的异步电动机变频调速仿真实现
SIMULATION FOR FREQUENCY CONTROL SYSTEM OF ASYNCHRONOUS MOTOR BASED ON MATLAB


本文主要对交流异步电动机SPWM变频调速矢量控制系统进行建模与仿真。变频调速系统在异步电动机的各种调速方式中效率最高、性能最好,因此有着极其重要的地位。电气传动控制系统计算机仿真是应用现代软件工具对其工作特性进行研究的一种十分重要的方法。通过仿真试验,可以比较各种策略与方案,优化并确定相关参数。因此进行系统仿真是不可或缺的,为科学决策提供了可靠的依据。
本文介绍了交流调速系统概况、矢量控制的基本概念以及异步电动机变频调速系统在 MATLAB/Simulink仿真工具中模型建立以及特性研究。一方面,本文通过对交流异步电动机矢量控制调速系统各部分仿真,得出该系统各部分的运行特性;另一方面,通过对转矩内环的转速、磁链闭环矢量控制系统和转差频率控制的异步电动机矢量控制系统的仿真,熟悉了矢量控制系统的参数设置和工作特性。
本文通过仿真实验不仅了解和掌握了异步电动机运行特性,更重要的是得出的仿真数据,为新的实验设备的引进和进一步开发打下了坚实的基础。

仿真图:

摘要
目录
1 绪论
1.1 概述
1.2 交流变频调速技术的现状
1.3 仿真工具MATLAB/Simulink简介
1.4 毕业设计的研究内容及章节安排
2 基于动态模型的异步电动机调速系统工作原理
2.1 异步电动机的数学模型
2.2 坐标变换
2.3 异步电动机在两相坐标系上的数学模型
3 交流异步电动机性能的仿真研究
3.1在交流情况下异步电动机工作仿真
3.2 PWM变频器-电动机系统仿真
PWM
4 交流异步电动机矢量控制调速系统仿真
4.1 转矩内环的转速、磁链闭环矢量控制系统仿真及分析
4.2 转差频率控制的异步电动机矢量控制系统仿真
4.3 两种矢量控制系统的特点与存在的问题
结论
参考文献
致谢


1 绪论

1.1 概述

电动机作为在工农业生产过程中主要的动力来源,发挥着日益重要的作用。电动机调速性能的好坏,对劳动生产率、节能环保产品质量有着非常重要的影响。直流电动机转速与交流电动机相比,更加容易被控制和调节。当为额定转速以下,可以通过维持励磁电流恒定,改变电枢电压方法来实现恒转矩调速;当为额定转速以上,可以通过维持电枢电压恒定,调节励磁的方法来实现恒功率调速。因此采用转速和电流双闭环直流调速系统可获得优良动态和静态调速特性。
然而,随着科技的不断发展和前进,直流调速弱点也逐步显现出来。由于直流调速系统中使用换向器,直流电动机维护工作量重,最高转速、单机容量以及使用环境等都受到限制,因此直流调速系统的应用范围狭窄。从20世纪30年代以来,交流调速技术逐步成为人们研究的热点,但交流电动机自身结构特点制约了调速技术发展。70 年代以来,伴随着电机控制技术及电力电子技术的飞速发展,交流调速系统性能也有了长足进步,几乎可以达到与直流调速系统相媲美的程度,为满足生产发展的需要,世界各国都非常重视交流调速技术的研发。
目前,交流调速系统已逐步取代直流调速系统。交流调速系统广泛采用交-直-交变频器作为交流电源,交-直-交变频器主要包括整流部分和逆变部分,整流部分将工频交流电整流成直流电,逆变部分再将直流电逆变成频率可调整的交流电,从而满足生产过程中对电动机调速系统的需求。
系统仿真是指使用计算机设备来模拟或描述一个系统或过程的运行行为,从而研究该系统或过程的一种技术。仿真技术是集系统理论、计算机技术、相似原理、控制论于一体的综合性技术。电气传动控制系统的计算机仿真是应用现代科学手段对其进行科学研究的十分重要的手段之一。近年来,几乎所有电动机调速控制系统的高品质控制均离不开系统仿真方面的研究。
1.2 交流变频调速技术的现状

目前,交流电机变频调速技术经已逐步成为当代电机调速的主流方法,日益得到人们的关注和重视。交流电机调速也由曾经的电磁离合器调速、变极调速、串级调速、转子串电阻调速、调压调速等有级调速方式逐步转变为无极调速方式。变频调速技术飞速发展为异步电动机已经可以实现无极平滑调速提供了可能。由于异步电动机变频调速系统还具有调速范围宽的优点,在工业控制领域中逐渐替代有极调速方式。交流变频调速克服了直流调速中的一些缺点,并拥有造价低、结构简单、坚固耐用等优点,已被公认为最有前途的调速方式。电工领域中各种容量和电压等级的变压器的出现,使得改变交流电压大小很容易实现,从而使得电力系统行业得到快速的发展。长期以来,交流电的频率保持固定不变,但是变频调速技术的出现使得交流电的频率可以可控地变化,这样就使得交流电得到充分利用。
随着我国经济高速发展,工业自动化领域中变频调速形成了一个巨大的市场。日本三菱、富士、德国西门子、ABB等国外公司先后进入中国市场并在中国建厂,国外公司产品最初约占中国变频器市场80%市场份额。经过近30 年的发展,国内像阿尔法、佳灵、利德华福、英威腾、汇川等近300多家国产变频器开拓市场,已经取得了很好的业绩。在国内变频器市场品牌中,内资品牌约占70%左右,但市场份额却只有30%左右,业内人士分析未来十年后市场才能达到饱和状态,因此该行业发展潜力十分可观。
变频调速技术的研究现状及发展趋势大体可归结为以下几个方面:
1、功率器件的发展。电力电子器件从上世纪五十年代开始,由第一代半控晶闸管 SCR及相关的派生器件,发展到全控器件诸如GTO,GTR,电力MOSFET,IGBT,IGCT 以及智能功率模块(IPM),IEGT(耐压达4KV以上的IGBT系列电力电子器件),主要向低通态电压、高压大容量、模块封装化方向发展。
2、控制策略的发展。变频器由最初的恒压频比控制方式(V/F),转差频率控制方式(SF)发展到矢量控制方式(VC)和直接转矩控制方式(DTC),调速精度和动态性能上均有明显提高。
3、微电子技术的发展。高性能矢量控制系统,若没有高效能微处理器支持,将很难实现。16位、32位高速微处理器及含浮点处理单元的DSP和专用集成电路ASIC(Application Specific IC)技术的快速发展,为实现变频调速器高效多功能化提供了硬件支持。
4、产品化生产。基础工业和各种工业制造业的高速发展,为变频器相关配套件生产提供了有力支持。
5、高压大功率变频调速装置以及永磁同步电动机及其控制系统的发展。近些年国际上高压大电流及串、并联技术和多电平技术的发展,使高压变频调速得以实现,有着明显的节能效果。
1.3 仿真工具MATLAB/Simulink简介

MATLAB是矩阵实验室(Matrix Laboratory)的简称,是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。MATLAB可以进行矩阵运算、绘制函数图像和数据、实现算法、创建用户界面、连接其他编程语言程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。现在MATLAB软件不但广泛应用于控制领域,也应用于其它的工程和非工程领域。
MATLAB具有很多优点,例如:此高级语言可用于技术计算;此开发环境可对代码、文件和数据进行管理;交互式工具可以按迭代的方式探查、设计及求解问题;数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以及数值积分等 ;二维和三维图形函数可用于可视化数据;各种工具可用于构建自定义的图形用户界面;各种函数可将基于MATLAB的算法与外部应用程序和语言(如C、C++、Fortran、Java、COM 以及Microsoft Excel)集成。
Simulink是一个用来对动态系统进行建模、仿真和分析的软件包。它为用户提供了用方框图进行建模的图形接口,使得建模就像用纸和笔来画画一样容易。Simulink允许用户定制和创建自己的模块。Simulink是建模仿真可视化软件,它用框图表示系统的各个环节,用带方向的连线,来表示各环节的输入输出关系。电力系统(SimPower System)模块库是专用于RLC电路、电力电子电路、电机传动控制系统和电力系统仿真用的模型库。模型库中包含了各种交直流电源,大量电气元器件和电工测量仪表等。利用这些模块可以模拟电力电子系统、电力拖动控制系统和电力系统运行和故障的各种状态运行。
1.4 毕业设计的研究内容及章节安排

本毕业设计主要用MATLAB软件对在矢量控制条件下的异步电动机变频调速系统的进行仿真,文章的章节安排如下:
第一章:绪论
对毕业设计研究背景进行说明,并对异步电动机调速系统所属领域及相关技术的发展进行概述,并且简单介绍了MATLAB/Simulink仿真工具。
第二章:基于动态模型的异步电动机调速系统的工作原理
介绍异步电动机数学模型,在异步电动机数学模型的基础上推演了电机的坐标变换,最后介绍了矢量控制的原理。
第三章:异步电动机性能的仿真研究
在MATLAB/Simulink中,通过仿真研究异步电动机的特性。首先是介绍鼠笼型交流异步电动机在正弦电压下的工作情况;然后是对交流异步电动机在变频器供电时的工作情况进行分析。
第四章:交流异步电动机矢量控制变频调速系统的仿真
在MATLAB/Simulink中,分别对转差频率控制的异步电动机矢量控制系统和带转矩内环的转速、磁链闭环矢量控制系统,在不同的给定转速条件下的仿真,并作出比较。

2 基于动态模型的异步电动机调速系统工作原理
异步电动机具有非线性、强耦合、多变量的性质,要获得高动态调速性能,必须从动态模型出发,分析异步电动机的转矩和磁链控制规律,研究高性能异步电动机的调速方案。矢量控制系统和直接转矩控制系统是已经获得成熟应用的两种基于动态模型的高性能交流电动机调速系统。矢量控制系统通过矢量变换和按转子磁链定向,得到等效的直流电动机模型,然后模仿直流电动机的控制策略设计控制系统。
2.1 异步电动机的数学模型
1.在研究异步电动机的多变量数学模型时,常做如下假设:
(1)忽略空间谐波,设三相绕组对称,在空间互差(电角度),所产生的磁动势沿气隙圆周按正弦规律分布;
(2)忽略磁路饱,各绕组的自感和互感都是恒定的;
(3)忽略铁心损耗;
(4)不考虑温度和频率的变化对电机参数的影响。
无论电动机转子是绕线型的还是鼠笼型的,都将它等效成绕线转子,到定子侧,折算后的每相绕组匝数都相等。这样,实际电动机就被等效为图示的三相异步电动机的物理模型。在三相异步电动机物理模型中,定子三相绕组轴线A、B、C在空间是固定的,故定义为三相静止坐标系。设A轴为参考坐标轴,转子以速度旋转,转子绕组轴线为a、b、c随转子旋转。转子a轴和定子A轴间的电角度差为空间角位移变量。如图2-1所示,为三相异步电动机物理模型。

图2-1 三相异步电动机物理模型

2.异步电动机三相动态模型的数学表达式
(1)磁链方程
异步电动机的数学模型由下述的电压方程、磁链方程、转矩方程和运动方程组成,其中磁链方程和转矩方程为代数方程,电压方程和运动方程为微分方程。
异步电动机每个绕组的磁链是它本身的自感磁链和其他绕组对它的互感磁链之和,因此,六个绕组磁链如下所示。
式中,L—6x6阶的电感矩阵,其中对角线元素是各相关绕组的自感,其余各项则是绕组间的互感。对于每一项绕组来说,它所交链的磁通是互感磁通与漏磁通之和,因此,定子和转子各相自感为:
两相绕组之间只有互感,互感的情况较为复杂,定子和转子的六个绕组之间的互感可考虑有两类:一类是A、B、C相绕组及a、b、c相绕组之间因位置固定,故互感为常数;另一类是定子任一相与转子任一相之间的位置是变化的,互感是角位移θ的函数。由于三相绕组的轴线在空间的相位差是电角度,在假设气磁通为正弦分布的情况下,有:
定、转子间的互感为:

矩阵方程写成分块矩阵的形式:
定子自感矩阵:
转子自感矩阵:
定子、转子之间的互感矩阵:
两个分块矩阵互为转置,且与转子位置有关,是时变元素,这是因为系统非线性的根源,可以用坐标变换把变参数转换成参数。
(2)电压方程
三相定子绕组电压平衡方程式为:
三相转子绕组折算到定子侧后的电压方程式为:
上述各量己折算到定子侧,将电压方程用矩阵形式表示可写成:
(3)转矩方程
按照机电能量转换原理,可求出电磁转矩的表达式:
式中,——电磁转矩;
——电机的磁极对数。
(4)运动方程
作用在电动机轴上的转矩与电动机速度变化之间的关系可以用运动方程来表达,一般情况下,电气传动系统的运动方程式为:
对于恒转矩负载D=K=0,则:
上述的异步电动机动态模型是在线性磁路、磁动势在空间按正弦分布的假定条件下得出来的,对定、转子电压和电流未作任何假设。因此,该动态模型完全可以用来分析含有电压、电流谐波的三相异步电动机调速系统的动态过程。
2.2 坐标变换
三相绕组可以用相互独立的两相正交对称绕组等效替代,等效原则是产生相等是磁动势。所谓独立是指两相绕组间无约束条件,即不存在约束条件;所谓正交是指两相绕组在空间上互差;所谓对称是两相绕组的匝数和阻值相等。图2-2中绘出两相绕组α、β,通过两相平衡交流电流,也能产生旋转磁动势。当三相绕组和两相绕组产生的两个旋转磁动势大小和转速都相等时,即认为两相绕组与三相绕组等效,这就是3/2变换。三相坐标系和两相坐标系物理模型如图2-2所示。

图2-2 三相坐标系和两相坐标系物理模型

异步电动机三相原始动态模型相当复杂,分析和求解这组非线性方程十分困难。异步电动机数学模型之所以复杂,关键是因为有一个复杂的电感矩阵和转矩方程,它们体现了异步电动机的电磁耦合和能量转换的复杂关系。因此,需从电磁耦合关系入手。
两相匝数相等相互正交的绕组d、q,分别通以直流电流,产生合成磁势F,其位置相对于绕组来说是固定的。如果人为地让包括两个绕组在内的整个铁心以同步转速旋转,则磁势F自然也随之旋转起来,成为旋转磁动势。如果这个旋转磁动势的大小和转速与固定的交流绕组产生的旋转磁动势相等,那么这套旋转的直流绕组也就和前面两套固定的交流绕组等效了。旋转正交坐标系的物理模型如图2-3所示。

图2-3 旋转正交坐标系的物理模型

由此可见,以产生相同的旋转磁动势为准则,三相交流绕组、两相交流绕组和旋转的直流绕组彼此等效。或者说,在三相坐标系下和在两相坐标系下的以及在旋转正交坐标系下的直流产生的旋转磁动势相等。在图2-3中的d、q两个绕组而言,当观察者站在地面上看上去,它们是与三相交流绕组等效的旋转直流绕组;如果跳到旋转着的铁心上看,它们就的的确确是一个直流电动机的物理模型。
这样,通过坐标系的变换,可以找到与交流三绕组等效的直流电动机模型。现在的问题是如何求出之间准确的等效关系,这就是坐标变换的任务。
1.三相-两相变换(3/2变换)
A、B、C为三相静止绕组,每相绕组的有效匝数为,通以三相平衡的正弦电流,产生合成磁动势F,以同步转速旋转,A、B、C轴称为三相静止坐标系;α、β为两相静止绕组,每相绕组的有效匝数为,它们在空间互差,且通入时间上互差的两相电流,也产生与上相同的磁动势F,并以同步转速旋转,α、β轴称为两相静止坐标系。当进行三相/两相坐标变换时,三相总磁动势应该与两相总磁动势相等,两套绕组瞬时总磁动势在α、β轴上的投影都应相等,写成矩阵形式,得:
在变换前后总功率不变的前提下,得到匝数比为:
三相静止坐标系等效变换到两相静止坐标系的变换矩阵:
如果从两相静止坐标系变换到三相静止坐标系,可以利用增广矩阵的方法,把扩成方阵,求其逆矩阵之后,再除以增加的一列,得:

为了便于利用功率不变条件下的坐标变换矩阵,需将变换矩阵变为方阵,因此,在坐标系中增设零轴,得:
上述公式同样适合于电压和磁链的变换矩阵。
2.静止两相-旋转正交变换(2s/2r变换)
从静止两相正交坐标系αβ到旋转正交坐标系dq的变换,称作静止两相-旋转正交变换,简称2s/2r变换。为两相静止坐标系统(2S)与磁场定向的M-T两相同步旋转坐标系统(2r)之间的变换称为两相/两相变换,简称变换,如图2-4所示。

图2-4 静止两相正交坐标系和旋转正交坐标系中磁动势矢量

图中两相交流电流,两相直流电流,产生同样的以同步速度旋转的合成磁动势F。由于各绕组匝数相等,可以消去磁动势中的匝数,可以用电流表示,但是此时的电流是空间矢量,而非时间相量。轴和轴静止不动,分量的大小却随时间而变化,相当于绕组在的交流磁动势的瞬时值。M轴和轴之间的夹角是随时间而变化的,而M轴和转子总磁链空间矢量方向重合,由图可知:
写成矩阵形式为:
式中,
为两相同步旋转坐标系变换到两相静止坐标系的变换矩阵。
反之,由两相静止坐标系变换到两相同步旋转坐标系的变换矩阵为:
电压和磁链的旋转变换矩阵也与电流(磁动势)的旋转变换矩阵相同。
2.3 异步电动机在两相坐标系上的数学模型
1.在任意两相旋转坐标系(dq坐标系)上的数学模型
两相坐标系可以是静止的,也可以是旋转的,其中任意旋转的坐标系是最常见的一种情况,由此,求某一具体两相坐标系上的数学模型就相对比较容易了。
设两相坐标d轴与三相坐标A轴的夹角为为dq坐标系相对于定子的角速度;为dq坐标系相对于转子的角速度。要把三相静止坐标系上的电压方程、磁链方程和转矩方程都变换到两相旋转坐标系上来,可以先利用3S/2S变换将方程式中的定子和转子的电流、电压、磁链和转矩都转换到两相静止坐标系上,然后再用旋转变换矩阵将这些变量都变换到两相旋转坐标系dq上。定子各量用下标1表示,转子各量用下标2表示。具体的变换过程比较复杂,变换后得到的数学模型如下。
(1)dq坐标系中的电压方程:
(2)dq坐标系中的磁链方程
数学模型简化的根本原因可从磁链方程和dq坐标系物理模型上看出。其磁链方程为:
由于变换到dq坐标系上以后,定子和转子等效绕组都落在两根轴上,而且两轴相互垂直,它们之间没有互感的耦合关系,互感磁链只在同轴绕组之间存在,所以式中每个磁链分量只剩下两项了。
(3)dq坐标系中的转矩方程和运动方程
将坐标变换矩阵代入ABC三相坐标系中的转矩方程式,简化后,得到dq0坐标系中的转矩方程为:
所以在dq0坐标系中的运动方程为:
2.按转子磁场定向的两相同步旋转坐标系上的数学模型
设d轴沿着转子磁链的方向,称之M轴;q轴逆时针旋转,即垂直于M轴,称为T轴;这样两相同步旋转的坐标系就具体规定为MT坐标系,即按转子磁场定向的坐标系。坐标系旋转速度等于同步角速度转子旋转机械角速度为;MT轴相对于转子的角速度为,即为转差的角速度。MT坐标系上的数学模型:


1



安徽理工大学毕业设计
3 交流异步电动机性能的仿真研究
本章主要通过仿真研究异步电动机的特性,为下一章交流异步电动机矢量控制调速系统的仿真打下基础。主要内容是在MATLAB/Simulink工作环境下,首先是对鼠笼型交流异步电动机在正弦交流电压下的工作情况进行仿真,参数设置以及出波形的特性;然后是对对在变频器供电时异步电动机的工作情况进行仿真,并对相关波形进行分析。
3.1在交流情况下异步电动机工作仿真
1.工作原理
当三相异步电机接三相交流电源时,三相定子绕组流过三相对称电流产生的三相磁动势(定子旋转磁动势)并产生旋转磁场。该旋转磁场与转子导体之间有相对切割运动,根据电磁感应原理,转子导体产生感应电动势并产生感应电流。载流的转子导体在磁场中受到电磁力作用,形成电磁转矩,驱动转子旋转,当电动机轴上带机械负载时,便向外输出机械能。三相异步电动机的转速永远低于旋转磁场的同步转速,使转子和旋转磁场间有相对运动,从而保证转子的闭合导体切割磁力线,感生电流,产生转矩。
2.在Simulink中交流异步电动机仿真模型的建立
在Simulink中异步电动机连接三相正弦电源时仿真模型如图3-1所示。
图3-1 交流条件下异步电动机仿真模型
异步电动机连接三相正弦电源,电动机负载由常数模块TL设定,电动机参数通过电动机测量模块(Machines Measurement Demux)测量,通过示波器观测电动机定子三相电流(isa,isb,isc)、转子三相电流(ira,irb,irc)、转速speed和转矩Te,并且由XY图示仪(XY Graph)观测电动机的机械特性(Te-n特性)。模型参数如表3-1所示。
表3-1 异步电动机特性研究模型参数
三相电源
电压
220V
频率
50Hz
电动机(容量Pn=4.73e3kW)
电压Un
380V
频率fn
50Hz
定子绕组电阻Rs
0.68
定子绕组漏感L1s
0.0042H
转子绕组电阻Rs
0.45
转子绕组漏感L1s
0.0042H
互感Lm
0.1486H
转动惯量J
0.05kg·
摩擦系数F
0.081
极对数P
2
负载设定TL
加载时间
0.5s
加载值
132
仿真参数
算法
Ode23t
相对误差
3.交流异步电动机的仿真波形
当异步电动机连接三相正弦交流电源时,异步电动机仿真波形如图3-2所示。
图3-2 异步电动机的仿真波形
4.对异步电动机仿真波形的分析
电动机在额定电压下空载起动时,图3-2分别描述了定子电流、转子电流、电动机转速、电动机转矩和机械特性变化情况。
在电动机启动到空载运行和过载运行过程中,定子电流和转子电流如图3-2(a)、(b)所示,在起动中随着转速的上升定子电流减小,在0.5s加载后定子电流迅速增大,定子电流为50Hz的正弦波。转子电流的变化与定子电流相同,但是从转子电流的波形可以看出,转子电流频率随电动机转差率的变化而变化,在启动过程中随转速上升转差率变小,转子电流频率下降,当电动机达到理想空载转速1500r/min,转子电流的大小和频率都是0,加上负载后随转速的下降和反向后转差变大,转子电流频率又增加。
由图3-2(c)、(d)可知,起动时电动机转速迅速上升,在0.2s时能够达到稳定转速1500r/min左右。用阶跃信号TL模拟,在t=0.5s时给电动机加上负载转矩132N·m,此时观察电动机转速知,电动机转速下降,转差率变大,在1.45s时转速下降为0,因为该负载远大于电动机额定负载40N·m,1.45s后电动机转速变为负值,这时相当于电动机待位能性负载,负载过大使得电动机处于倒拉反转的状态。图3-2(d)是电动机的转矩响应,起动中交流电动机的转矩有波动的,严重过载时引起电动机反转时,电动机产生很小的转矩。
图3-2(e)是电动机的动态机械特性,该机械特性与理论上异步电动机起动、空载运行以及带负载运行的机械特性相一致。在异步电动机起动时起动转矩变化比较大,此时电动机的转速较低;当空载运行时,转差率较小,电动机转速较高;在t=0.5s时,给电动机加上负载转矩132N·m,远大于电动机额定负载40N·m,此时异步电动机的工作点会下移,越过转折点,异步电动机进入不稳定区。
3.2 PWM变频器-电动机系统仿真
1.在Simulink中PWM变频器-电动机系统仿真模型的建立
模型由三相交流电源模块、两个通用桥模块(Universal Bridge)、调制器(PWM Generator)、电感L2、电容C1、异步交流电动机模块(Asynchronous Machine SI Units)以及电机测量(Machines Measurement Demux)组成。
仿真电路将三相交流电通过由通用桥模块(Universal Bridge)组成的整流单元整流,并经过电感、电容滤波,得到近似直流电。该直流电经过调制器(PWM Generator)控制下的逆变单元给异步电动机供电,其中逆变单元也是由通用桥模块组成,选择的半导体器件为IGBT。整流单元和逆变单元组成了变频器的主电路,调制器(PWM Generator)相当于控制单元。交流异步电动机由PWM变频器供电的仿真模型如图3-3所示。

图3-3PWM变频器-交流异步电动机系统仿真模型
2.PWM变频器-电动机系统仿真波形及其分析
PWM变频器-电动机系统仿真波形如图3-4所示。图3-4分别描绘了逆变器输出电压、电动机定子电流、电动机转子电流以及转速的变化波形。
图3-4 逆变器供电电动机工作波形
图3-4(a)是逆变器输出线电压Uab的波形,PWM发生器给逆变器提供驱动信号,控制开关管IGBT的导通,从而控制Uab输出波形。
图3-4(b)、(c)分别是电动机定子电流和转子电流,空载启动运行时,刚刚启动
时定子电流和转子电流很大,随着电动机渐渐达到额定转速时,定子和转子电流都逐渐减小;当t=0.15s时,异步电动机转速达到空载转速,此时定子电流和转子电流为最小值;当在t=0.25s,通过阶跃信号模块TL突加负载转矩45N·m,定子电流和转子电流幅值分别逐渐上升,最终达到稳定状态。定子电流和转子电流变化的根本原因是:为了维持异步电动机转速的恒定和稳定运行,必须使电机产生的电磁转矩和负载转矩相等,由异步电动机的等效T型模型可知,电磁转矩Te与定子电流或者转子电流成正比关系。所以负载转矩的不断变化使得定子和转子电流的不断变化。开始阶段,异步电动机的起动转矩较大,导致了起动电流较大;空载运行时,负载转矩很小,导致了此时定子和转子电流较小;突加负载后,负载转矩上升,导致了定子和转子的电流上升。这样我们就看到了图3-4(b)、(c)定子电流和转子电流的变化情况。
图3-4(d)是异步电动机从空载启动达到额定转速再负载运行时转速变化情况。由图可知,异步电动机从空载启动,转速上升平稳,大约在0.2s时达到额定转速约为1500r/min;在t=0.25s突加负载后,转速微微下降,此时电磁转矩上升,定子电流和转子电流也有所上升;最终异步电动机转速达到稳定状态,波形中有微微波动,不是一条直线,这主要是因为交流电源经整流模块和滤波单元后,由于滤波不完全,中间存在各次谐波,导致输出直流电波形存在扰动情况。如果用三相电压型SPWM逆变器代替PWM逆变器充当电源供给,输出波形可以更好点。

1



安徽理工大学毕业设计
4 交流异步电动机矢量控制调速系统仿真
按转子磁链定向矢量控制的基本思想是通过坐标变换,在按转子磁链定向同步旋转正交坐标系中,得到等效的直流电动机模型,仿照直流电动机的控制方法来控制电磁转矩与磁链,然后将转子磁链定向坐标系中的控制量通过反变换得到三相坐标系的对应量,以实施控制。由于变换的是矢量,所以这样的坐标变换也可称作矢量变换,相应的控制系统称为矢量控制(Vector Control,VC)系统或按转子磁链定向控制(Flux Orientation,FOC)系统。
本章主要对带转矩内环的转速和磁链闭环异步电动机矢量控制系统、转差频率控制的异步电动机矢量控制系统以及无速度传感器的矢量控制系统仿真,并通过改变模型参数比较和分析交流异步电动机矢量控制调速系统的特性,加强对交流异步电动机矢量控制系统的理解。
4.1 转矩内环的转速、磁链闭环矢量控制系统仿真及分析
带转矩内环的转速、磁链闭环矢量控制系统主电路采用电流滞环控制型逆变器。在控制电路中,在转速环后增加了转矩控制内环,转速调节器ASR的输出是转矩调节器ATR的给定,而转矩的反馈信号,则通过矢量控制方程计算得到。电路中磁链调节器ApsiR用于电动机定子磁链的控制,并设置了电流变换和磁链观测环节。ATR和ApsiR的输出分别是定子磁链的控制,并设置了电流变换和磁链观测环节。ATR和APsiR的输出分别是定子电流的转矩分量和励磁分量经过2r/3s变换后得到三相定子电流的给定值,并通过电流滞环控制PWM逆变器控制电动机定子的三相电流。
1.带转矩内环转速、磁链闭环控制系统仿真模型
在带转矩内环转速、磁链闭环控制系统仿真模型中,包含220V三相交流电源、两个通用桥模块(Universal Bridge)、调制器(PWM Generator)、电感L2、电容C1、异步交流电动机模块(Asynchronous Machine SI Units)以及电机测量(Machines Measurement Demux)组成了模型的主电路,逆变器的驱动信号由滞环脉冲发生器产生。三个调节器ASR、ATR和ApsiR是带输出限幅的PI调节器(见图4-2~图4-4)。转子磁链观测使用两相同步旋转坐标系上的磁链模型(Current model),函数模块Fcn用于计算转矩,两相旋转坐标系/三相静止坐标系变换模块(dq0-to-abc)实现了矢量控制中的2r/3s的坐标变换。带转矩内环转速、磁链闭环控制系统仿真模型如图4-1所示。
图4-1 带转矩内环的转速、磁链闭环的矢量控制系统仿真模型
ASR为转速调节器,ATR为转矩调节器,ApsiR为磁链调节器。三个调节器结构相同,只是内部参数设置不同。带输出限幅的PI调节器ASR、ATR和ApsiR都含有饱和模块(Saturation),其主要作用是限制信号的范围,使信号变化在符合实际要求的范围之内。ASR为转速调节器如图4-2所示,ATR为转矩调节器如图4-3所示,ApsiR为磁链调节器如图4-4所示。
图4-2 ASR调节器结构
图4-3 ATR调节器结构
图4-4 ApsiR调节器结构
2.仿真模型参数设定
电动机参数同表3-1异步电动机仿真模型参数的设定。模型仿真算法为ode23tb,ode23tb在龙格-库塔法的第一阶段用梯形法,第二阶段用二阶的backward differentiation formulas算法。由于仿真模型中采用了双闭环的控制方式,系统抗扰动能力较强,因此采用该算法可以在不影响系统精度的前提下,较为迅速地得出仿真结果。在仿真模型中,给定转速为1400r/min,空载起动,在0.6s时突加负载60N·m。ASR、ATR、ApsiR调节器参数如表4-1所示。
表4-1 调节器参数参考值
调节器

比例放大器放大倍数
积分放大器放大倍数
积分器限幅
调节器输出限幅
上限
下限
上限
下限
转速调节器ASR
3.8(G1)
0.8(G2)
80
-80
75
-75
转矩调节器ATR
4.5(G3)
12(G4)
60
-60
60
60
磁链调节器ApsiR
1.8(G5)
100(G6)
15
-15
13
-13
3.仿真输出波形及分析
(1)在给定转速=1400r/min时,系统仿真波形如图4-5,4-6所示。
图4-5(a)带转矩内环转速和磁链闭环控制系统在给定转速=1400r/min时,转速响应,从波形中可以看出,在矢量控制下,在起动阶段,转速上升平稳;在t=0.35s达到给定转速,此时为空载状态下的转速;当t=0.6s时,给电机加负载,TL=60N·m后,电动机转速略有下降。
图4-5 =1400r/min时带转矩内环的转速和磁链闭环矢量控制系统仿真波形
图4-5(c)为异步电动机转矩波形,空载起动时,起动转矩较大,这时电动机产生较大的电磁转矩,使得异步电动机的定子电流也比较大;当起动结束后,空载转速接近同步转速,此时负载转矩很小接近为0;当t=0.6s时,加上负载后,转速不能突变,电磁转矩增加,电磁转矩Te近似等于负载转矩60Nm。所以所得仿真波形与理论上分析的波形基本一致。
图4-5(b)为A相定子电流,图4-5(d)经过2r/3s变换的A相电流给定值,比较两者波形可知观察给定值和测量值波形可知,两者波形近似相同,即波形始终围绕着值在规定的范围内波动。在起动过程中,由于定子电流的转矩和励磁分量都保持不变,定子电流的给定值也不变,所以在起动过程中,定子电流基本保持不变,实现恒电流起动;在t=0.35s时,此时电动机达到空载状态下的额定转速,当忽略异步电动机T型等效电路中的励磁支路时,有,所以t=0.35s到t=0.6s时由于电动机空载运行,所以A相定子电流ia值比较小;t=0.6s加上负载TL=60N·m后,电磁转矩增大,定子电流也随之增大。
=1400r/min时带转矩内环的转速、磁链闭环矢量控制系统中异步电动机定子磁链轨迹和转矩-转速曲线如图4-6所示。
图4-6 =1400r/min时系统仿真波形
图4-6(a)、(b)分别为异步电动机定子磁链轨迹和转矩-转速曲线。比较两者波形可以看出,在起动阶段,磁场建立过程比较平滑,磁链呈螺旋形增加,最终形成圆形旋磁场。同时电动机转矩也不断上升,这有赖于磁链闭环矢量控制方式的作用。当异步电动机转子磁链发生波动时,电磁转矩不是稳定的,电动机转速也受到影响。此时,转子磁链调节器力图使转子磁链恒定,而转速调节器则调节电流转矩分量,以抵消转子磁链的变化对电磁转矩的影响,最后达到平衡状态。
(2)在给定转速=1700r/min时,仿真输出波形
=1700r/min时闭环矢量控制系统仿真波形如图4-7、4-8所示。
图4-7 =1700r/min时系统仿真波形
图4-7、4-8系统仿真波形总体上和图4-5、4-6相仿,给定值=30*/,通过改变转速给定,来改变转子角速度给定值。图4-7(a)带转矩内环转速和磁链闭环控制系统在给定转速=1700r/min时,转速响应,从波形中可以看出,在矢量控制下,在起动阶段,转速上升平稳;在t=0.4s时,空载状态下的转速接近给定转速=1700r/min;当t=0.6s时,给电机加负载,TL=60N·m后,电动机转速略有下降。
图4-7(b)为A相定子电流,图4-7(d)经过2r/3s变换的A相电流给定值,比较两者波形可知观察给定值和测量值波形可知,两者波形近似相同,即波形始终围绕着值在规定的范围内波动。在起动过程中,由于定子电流的转矩和励磁分量都保持不变,定子电流的给定值也不变,所以在起动过程中,定子电流基本保持不变,实现恒电流起动;在t=0.4s时,此时电动机达到空载状态下的额定转速,当忽略异步电动机T型等效电路中的励磁支路时,有,所以t=0.4s到t=0.6s时由于电动机空载运行,所以A相定子电流ia值比较小;t=0.6s加上负载TL=60N·m后,电磁转矩增大,定子电流也随之增大。
图4-7(c)为异步电动机转矩波形,空载起动时,起动转矩较大,这时电动机产生较大的电磁转矩,使得异步电动机的定子电流也比较大;当起动结束后,空载转速接近同步转速,此时负载转矩很小接近为0;当t=0.6s时,加上负载后,转速不能突变,电磁转矩增加,电磁转矩Te近似等于负载转矩60N·m。
=1700r/min时带转矩内环的转速、磁链闭环矢量控制系统中异步电动机定子磁链轨迹和转矩-转速曲线如图4-8所示。
图4-8 =1700r/min时系统仿真结果
图4-8(a)、(b)分别为=1700r/min异步电动机定子磁链轨迹和转矩-转速曲线。比较两者波形可以看出,在起动阶段,磁场建立过程比较平滑,磁链呈螺旋形增加,最终形成近似六边形的旋转磁场。同时电动机转矩也不断上升,这有赖于磁链闭环矢量控制方式的作用。比较图4-8(a)和图4-6(a)知,图4-6(a)中旋转磁场更加接近圆形,这是因为图4-8(a)中PWM开关频率设定值较低,而PWM开关频率越高旋转磁场越接近为圆。所以在实际情况中提高PWM开关频率,有助于建立圆形旋转磁场。
通过比较=1400r/min和=1700r/min条件下,带转矩内环的转速和磁链闭环矢量控制系统仿真结果,我们可以知道,在两种条件下,实际电流都能快速的跟随给定值,充分发挥了闭环控制的优越性。因为根据异步电动机矢量变换及等效直流电动机模型可知,转子磁链定向实现了定子电流的转矩分量和励磁分量的解耦,转子的角速度主要受定子电流的转矩分量控制。采用电流闭环控制,通过改变转子角速度的给定值,通过反馈环节可以实现对定子电流转矩分量的控制。当系统受到扰动时,定子电流转矩分量虽然发生变化,但是由于反馈环节的作用,能够很快的恢复到扰动之前的状态。
纵观图4-4带转矩内环的转速和磁链闭环的矢量控制系统仿真模型,转速给定通过转速调节器(ASR)得到转矩调节器(ATR)的给定转矩,而仿真模型中增加了转矩内环,和矢量控制方程计算出的反馈信号在转矩调节器(ATR)中作用,得出定子电流的转矩分量给定值。转矩内环的作用是:当转子磁链发生波动时,通过转矩调节器及时调整电流转矩分量的给定值,以抵消磁链变化的影响,尽可能不影响或者少影响电动机转速。另一方面,转子磁链给定值为与通过转子磁链定向转子磁链电流模型得出的转子磁链比较,经过转子磁链调节器ApisR,得出电子电流的励磁分量。磁链调节器作用是:力图使转子磁链恒定。ATR和APsiR的输出分别是定子电流的转矩分量和励磁分量经过2r/3s变换后得到三相定子电流的给定值,并通过电流滞环控制PWM逆变器控制电动机定子的三相电流。
4.2 转差频率控制的异步电动机矢量控制系统仿真            
1.转差频率控制异步电动机矢量控制系统的工作原理
以上介绍的磁链闭环控制的矢量控制系统中,磁链幅值和位置信号均由磁链模型计算得到,都受电机参数变化的影响,造成控制系统的不不准确性。采用磁链开环的控制方式,无需转子磁链幅值,但对于矢量变换而言,仍然需要转子磁链的位置信号,转子磁链的计算仍然不可避免。如果利用给定值间接计算转子磁链的位置,可简化系统结构,这种方法称为间接定向。间接定向的矢量控制系统借助于矢量控制方程中转差公式,构成转差型的矢量控制系统。
转速采取了转差频率控制,即异步电动机定子角频率由转子角频率组成(=+),这样在转速变化过程中,电动机定子电流频率始终能随转子的实际转速同步升降,使转速调节更为平滑。
从矢量控制方程中可以看出,在保持转子磁链不变的控制下,电动机转矩直接收定子电流的转矩分量控制,并且转差可以通过定子电流转矩分量计算,转子磁链也可以通过定子电流的励磁分量来计算。在系统中以转速调节器ASR的输出为定子电流的转矩分量,并通过计算得到。如果采用磁通不变控制,则p=0,而
由于矢量控制方案得到的是定子电流的励磁分量和转矩分量,而本系统采用了电压型逆变器,需要相应的将电流控制转换为电压控制,其变换关系为:
式中,为定子电压的励磁分量和转矩分量;为漏磁系数,
经过两相旋转坐标系/三相静止坐标系的变换(2r/3s),得到SPWM逆变器的三相电压控制信号,并控制逆变器的输出电压。
2.转差频率控制异步电动机矢量控制系统参数设置及仿真
转差频率控制异步电动机矢量控制系统仿真模型如图4-9所示。

图4-9 转差频率控制的矢量控制系统仿真模型
系统的控制部分由给定、PI调节器、函数运算、两相/三相坐标变换、PWM脉冲发生器等环节组成。其中给定环节有定子电流励磁分量im*和转子速度n*。放大器G1、G2和积分器组成带限幅的转速调节器ASR.。电流电压模型转换由函数Um*、Ut*模块实现。函数运算模块ws*根据定子电流的励磁分量和转矩分量计算转差,并与转子频率相加得到定子频率,再经积分器得到定子电压矢量转角。模块sin、cos、dq0/abc实现了两相旋转坐标系至三相静止坐标系变换。dq0/abc的输出时PWM发生器的三相调制信号,因为调制信号幅度不能大于1,在dq0/abc输出后插入了衰减环节。PWM发生器设置为内部模式,然后运行模型,根据dq0/abc输出和PWM发生器的三相调制输入信号幅值小于1的要求,计算G4的衰减系数。各放大器的参数取值见表4-2。
表4-2 转差频率矢量控制仿真模型放大器参数
放大器
放大倍数
放大器
放大倍数
G1
35
G4
2
G2
0.15
G5、G6
9.55
G3
0.0076


仿真了给定转速为1400r/min时的空载起动的过程,在起动后0.45s加载TL=65N·m。该系统是比较复杂的系统,收敛是仿真计算中经常出现问题,经试用各种计算方法,选择了固定步长算法ode5,步长取
3.转差频率控制的矢量控制系统仿真波形及其分析
转差频率控制的矢量控制系统仿真波形如图4-10所示。
图4-10 转差频率矢量控制系统仿真结果
图4-10(a)是电动机的转速响应,图4-10(b)是定子电流A相电流响应,图4-10(c)是异步电动机的电磁转矩和负载转矩给定。
图4-10(a)是电动机的转速响应,电机在起动过程中,转速上升平稳,此时异步电动机产生的电磁转矩近似等于起动转矩,起动电流较大;当t=0.4s时,电动机转速稳定在空载运行转速上,此时转速接近同步转速,约1400r/min,异步电动机产生的电磁转矩下降为空载转矩,近似为零,定子电流也比较小;当t=0.45s是,突加负载TL=65N·m,此时电磁转矩上升,转速略有下降,但不是很明显,并且最终稳定下来。
图4-10(b)是定子电流A相电流响应,起动过程中,由于电磁转矩近似等于起动转矩,起动转矩较大,从而导致起动电流较大,即异步电动机的定子电流较大,幅值接近50A;当t=0.4s时,电磁转矩近似等于空载转矩,非常小,近似为零,从而使得A相定子电流很小,接近为0;当t=0.45s时,突加负载转矩TL=65N·m,为了维持电动机稳定运行,电磁转矩上升,等于负载转矩,这就要求定子电流上升,以产生足够大的电磁转矩。
图4-10(c)是异步电动机的电磁转矩和负载转矩给定,当t=0时,起动转矩较大,转矩最大值超过了100N·m;当t<0.05s时,异步电动机处于刚刚起动状态,起动转矩波动比较大,异步电动机运行不是很平稳,会出现抖动等不稳定的状况;当0.05s< t<0.2s时,电磁转矩逐步稳定,稳定在70 N·m左右;当t=0.4s时,电动机处于空载运行状态,负载转矩接近为0,电磁转矩也接近为0;当t=0.45s时,突加负载转矩TL=65N·m,负载转矩突然上升,电磁转矩也跟随上升,此时,电磁转矩近似等于负载转矩给定。
4.3 两种矢量控制系统的特点与存在的问题
对于带转矩内环的转速、磁链闭环矢量控制系统,在本仿真中主电路采用电流滞环控制型逆变器。在计算转子磁链的电流模型都需要实测的电流和转速信号,无论转速高低时都能适用,但都受电动机参数影响。例如在实际情况下,电动机温升和频率变化都会影响转子电阻,磁饱和程度将影响电感。这些影响都将导致磁链幅值与位置信号失真,而反馈信号的失真比然使磁链闭环控制系统的性能降低,这是带转矩内环的转速、磁链闭环矢量控制系统的不足。
对于磁链开环转差型矢量控制系统的磁场方向由磁链和电流转矩分量给定信号确定,靠矢量控制方向保证,并没有磁链模型实际计算转子磁链及其相位,所以属于间接的磁场定向。但由于矢量控制方程中包含了电动机转子参数,定向精度仍受参数变化影响,磁链和电流转矩分量给定值与实际值存在差异,将影响系统性能。
综上可得,两种矢量控制系统存在的共同特点:
(1)按转子磁链定向,实现定子电流励磁分量和转矩分量的解耦,需要电流闭环控制;
(2)转子磁链控制在系统的控制对象是稳定的惯性环节,可以采用磁链闭环控制,也可以采用开环控制;
(3)采用连续的PI控制,转矩与磁链变化平稳,电流闭环控制可有效限制起、制动电流。
同时,两种矢量控制系统也存在着一些共同的问题:
(1)转子磁链计算精度易受变化的转子电阻的影响,转子磁链的角度精度影响定向的准确性,使得系统的误差变大。
(2)需要进行矢量变换,系统结构复杂,运算量大。


结论
本文通过在Simulink环境下建立交流异步电动机矢量控制仿真模型,作为研究交流电动机变频调速的基础。本文第一章是绪论部分,主要是介绍毕业设计研究背景和总体思路;第二章主要介绍了基于动态模型的异步电动机调速系统的理论基础;第三章是对交流异步电动机矢量控制系统各部分单元进行仿真,第四章打下基础;第四章对两种矢量控制下的变频调速系统进行仿真,通过仿真实验研究了系统的性能。
仿真实验结果表明了通过矢量变换建立的电机模型的正确性和有效性。系统仿真表明,适当地选取调节参数和合适的仿真参数和仿真算法,可使该系统得到良好的动态性能,对给定值具有良好的动态跟踪能力,控制系统的各个部分协调运行,实现了优化,获得了具有理想性能的调速系统。同时也可以看到,矢量控制是控制学科发展历程中的一个崭新的阶段。目前,矢量控制的研究与应用己经深入到众多的领域;同样,它的发展也给电气传动系统的控制策略带来了新思想、新方法。SIMULINK软件在矢量控制交流调速系统中确实为优秀而便利的软件,必将获得更广泛的应用。
通过大量的理论计算和仿真实验研究知,实现矢量控制的方法不是惟一的,比如,为了实现磁场定向控制,需要测出转子磁链的幅值和相位,但是是电动机内部的物理量,直接检测在技术上难以实现。因此,多采用间接观测的方法,即检测与有关的电动机运行参数,如定子电压、电流或转速等,然后根据电动机的动态数学模型,通过运算,实时计算转子磁链的幅值和相位。此外也可以根据系统运行的指令值和检测到的转子位置信号或转子速度信号,通过计算求得的幅值和相位。同时利用能够实际测量的物理量的不同组合,可以得到不同形式的各有利弊的转子磁链观测模型。


以上的Word格式文档51黑下载地址:
基于matlab的异步电机变频调速系统的设计.doc (1.43 MB, 下载次数: 217)



分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏5 分享淘帖 顶 踩
回复

使用道具 举报

沙发
ID:395889 发表于 2018-9-10 10:05 | 只看该作者
我也是这个方向的的嘞,感谢分享
回复

使用道具 举报

板凳
ID:428261 发表于 2018-11-21 08:56 | 只看该作者
楼主你好,能不能分享一下磁链观测器和dq0-to-abc模块怎么做出来的嘛
回复

使用道具 举报

地板
ID:489394 发表于 2019-3-12 17:50 | 只看该作者
想看下~能下载吗
回复

使用道具 举报

5#
ID:502499 发表于 2019-4-1 16:00 | 只看该作者
您好,请问可以分享一下转差频率控制的矢量控制系统仿真模型(即图4-9)的原图吗?或者MATLAB文件,谢谢!
回复

使用道具 举报

6#
ID:532428 发表于 2019-5-10 13:28 | 只看该作者
hehe22 发表于 2018-9-10 10:05
我也是这个方向的的嘞,感谢分享

请问您做完了吗
回复

使用道具 举报

7#
ID:565832 发表于 2019-6-18 15:54 | 只看该作者
收获颇多,感谢!
回复

使用道具 举报

8#
ID:764967 发表于 2020-5-30 13:28 | 只看该作者
您好 请问在么
您的那篇基于MATLAB的异步电机变频调速系统的设计文档分析的很全面 内容也很到位   我目前的方向和这篇有很大关联  恳请您能把仿真部分的文件或代码借我参考  
回复

使用道具 举报

9#
ID:893984 发表于 2021-3-19 16:37 | 只看该作者
大佬,我看看,现在正在做这方面
回复

使用道具 举报

10#
ID:905007 发表于 2021-4-13 13:48 | 只看该作者
能给个仿真参数不?
回复

使用道具 举报

11#
ID:841381 发表于 2021-6-15 19:22 | 只看该作者
能把源码和仿真模型发附件到论坛给大家参考下吗!谢谢
回复

使用道具 举报

12#
ID:395359 发表于 2021-10-21 08:58 | 只看该作者
好东西呀!想学习好久了!谢谢!
回复

使用道具 举报

13#
ID:395359 发表于 2021-10-25 15:50 | 只看该作者
太牛了!但看不明白!真的想学习呀!
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

快速回复 返回顶部 返回列表