第四周仿真作业

经过第四周关于交流电机的学习，我们通过对起调速特性的了解，进行仿真：

  • 控制电机带重物上升，从静止加速到800r/min
  • 保持800r/min匀速运动0.5s，
  • 减速到静止，保持静止状态0.5s，
  • 带重物下降，从静止达到600r/min
  • 保持600r/min匀速运动0.6s，
  • 减速到静止。
    (为了便于仿真，匀速和静止持续时间较短)

分析：这里，由于自耦降压启动具有K值方便直接调节的特点，所以我们采用此种启动方法。由于变频调速具有调速范围广、平滑性好、能耗小无级调速，故采用变频调速。制动采用相对方便的反接制动。

参数：我们需要设置K，a,b,c 4个参数，由于不方便计算，都是通过观察图像进行调试，从而得出近似值。我是将K先固定，通过调节a,使转速稳定时接近800，调节b,使其接近制动，调节c，使转速稳定时接近600.

从而我这里得出的参数K=0.8，a=0.54,b=0.07,c=0.405.

代码：

model SACIM "A Simple AC Induction Motor Model"
type Voltage=Real(unit="V");
type Current=Real(unit="A");
type Resistance=Real(unit="Ohm");
type Inductance=Real(unit="H");
type Speed=Real(unit="r/min");
type Torque=Real(unit="N.m");
type Inertia=Real(unit="kg.m^2");
type Frequency=Real(unit="Hz");
type Flux=Real(unit="Wb");
type Angle=Real(unit="rad");
type AngularVelocity=Real(unit="rad/s");

constant Real Pi = 3.1415926;

Current i_A"A Phase Current of Stator";
Current i_B"B Phase Current of Stator";
Current i_C"C Phase Current of Stator";
Voltage u_A"A Phase Voltage of Stator";
Voltage u_B"B Phase Voltage of Stator";
Voltage u_C"C Phase Voltage of Stator";
Current i_a"A Phase Current of Rotor";
Current i_b"B Phase Current of Rotor";
Current i_c"C Phase Current of Rotor";
Frequency f_s"Frequency of Stator";
Torque Tm"Torque of the Motor";
Speed n"Speed of the Motor";

Flux Psi_A"A Phase Flux-Linkage of Stator";
Flux Psi_B"B Phase Flux-Linkage of Stator";
Flux Psi_C"C Phase Flux-Linkage of Stator";
Flux Psi_a"a Phase Flux-Linkage of Rotor";
Flux Psi_b"b Phase Flux-Linkage of Rotor";
Flux Psi_c"c Phase Flux-Linkage of Rotor";

Angle phi"Electrical Angle of Rotor";
Angle phi_m"Mechnical Angle of Rotor";
AngularVelocity w"Angular Velocity of Rotor";

Torque Tl"Load Torque"; 

Resistance Rs"Stator Resistance";
parameter Resistance Rr=0.408"Rotor Resistance";
parameter Inductance Ls = 0.00252"Stator Leakage Inductance";
parameter Inductance Lr = 0.00252"Rotor Leakage Inductance";
parameter Inductance Lm = 0.00847"Mutual Inductance"; 
parameter Frequency f_N = 50"Rated Frequency of Stator";
parameter Voltage u_N = 220"Rated Phase Voltage of Stator";
parameter Real p =2"number of pole pairs";
parameter Inertia Jm = 0.1"Motor Inertia";
parameter Inertia Jl = 0.1"Load Inertia";
parameter Real K=0.8"starting rate";
parameter Real a=0.54"frequency rate";
parameter Real b=0.07"stable frequency rate"; 
parameter Real c=0.405"another frequency rate";
parameter Real P=0.7"stoping rate"; 

initial equation

Psi_A = 0; 
Psi_B = 0;
Psi_C = 0;
Psi_a = 0; 
Psi_b = 0;
Psi_c = 0;
phi = 0;
w = 0;

equation

u_A = Rs * i_A + 1000 * der(Psi_A);
u_B = Rs * i_B + 1000 * der(Psi_B);
u_C = Rs * i_C + 1000 * der(Psi_C);
= Rr * i_a + 1000 * der(Psi_a);
= Rr * i_b + 1000 * der(Psi_b);
= Rr * i_c + 1000 * der(Psi_c);

Psi_A = (Lm+Ls)*i_A + (-0.5*Lm)*i_B + (-0.5*Lm)*i_C + (Lm*cos(phi))*i_a + (Lm*cos(phi+2*Pi/3))*i_b + (Lm*cos(phi-2*Pi/3))*i_c;
Psi_B = (-0.5*Lm)*i_A + (Lm+Ls)*i_B + (-0.5*Lm)*i_C + (Lm*cos(phi-2*Pi/3))*i_a + (Lm*cos(phi))*i_b + (Lm*cos(phi+2*Pi/3))*i_c;
Psi_C = (-0.5*Lm)*i_A + (-0.5*Lm)*i_B + (Lm+Ls)*i_C + (Lm*cos(phi+2*Pi/3))*i_a + (Lm*cos(phi-2*Pi/3))*i_b + (Lm*cos(phi))*i_c;

Psi_a = (Lm*cos(phi))*i_A + (Lm*cos(phi-2*Pi/3))*i_B + (Lm*cos(phi+2*Pi/3))*i_C + (Lm+Lr)*i_a + (-0.5*Lm)*i_b + (-0.5*Lm)*i_c;
Psi_b = (Lm*cos(phi+2*Pi/3))*i_A + (Lm*cos(phi))*i_B + (Lm*cos(phi-2*Pi/3))*i_C + (-0.5*Lm)*i_a + (Lm+Lr)*i_b + (-0.5*Lm)*i_c;
Psi_c = (Lm*cos(phi-2*Pi/3))*i_A + (Lm*cos(phi+2*Pi/3))*i_B + (Lm*cos(phi))*i_C + (-0.5*Lm)*i_a + (-0.5*Lm)*i_b + (Lm+Lr)*i_c;

Tm =-p*Lm*((i_A*i_a+i_B*i_b+i_C*i_c)*sin(phi)+(i_A*i_b+i_B*i_c+i_C*i_a)*sin(phi+2*Pi/3)+(i_A*i_c+i_B*i_a+i_C*i_b)*sin(phi-2*Pi/3));

w = 1000 * der(phi_m);

phi_m = phi/p;
n= w*60/(2*Pi);

Tm-Tl = (Jm+Jl) * 1000 * der(w);
Tl = 10;

if time <= 100 then
u_A = 0;
u_B = 0;
u_C = 0;
f_s = 0;Rs = 0.531;
elseif time<=106 then
f_s = f_N*a;

Rs = 0.531;
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*K*a; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*K*a;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*K*a; 
elseif time<=1706 then
f_s = f_N*a;

Rs = 0.531;
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*a; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*a;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*a; 
elseif time<=1756 then
f_s = f_N*a;Rs = 3;
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*a; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*a;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*a;
elseif time<=2679 then
f_s = f_N*a;Rs = 0.531;
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*a; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*a;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*a;
elseif time<=3196 then
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*b; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*b;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*b; 
f_s = f_N*b;Rs = 0.531;
elseif time<=3310 then
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*K*c; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*K*c;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*K*c; 
f_s = f_N*K*c;Rs = 0.531;

elseif time<=4560 then
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*c; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*c;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*c; 
f_s = f_N*c;Rs = 0.531;
elseif time<=4595.5 then
f_s = f_N*P*a;Rs = 2;
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*a*P; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*a*P;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*a*P; 
elseif time<=5560 then
f_s = f_N*a;

Rs = 0.531;
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*a; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*a;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*a; 
else
u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000)*b; 
u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3)*b;
u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3)*b; 
f_s = f_N*b;Rs = 0.531;
end if;

end SACIM;

得出的仿真图像：

 

转载于:https://www.cnblogs.com/cyykhyql/p/5305839.html

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

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

相关文章

GVim中以十六进制方式打开文件

GVim中以十六进制方式打开文件

在windows下&#xff0c;有时我们想用GVim以十六进制形式来打开一个文件&#xff0c;下面就是不错的方法&#xff0c;将它作为一个右键菜单方式加入&#xff1a; 1. 将"GVim Hex"加入到右键菜单, 将以下内容添加到一个xxx.reg文件&#xff0c;再双击运行一次即可&am…
阅读更多...
高性能的智能日志

高性能的智能日志

【编者按】本文作者是 Archanaa Panda &#xff0c;从 2000 以来一直在软件开发&#xff08;构架、设计和编程&#xff09;团队担任 Java / JavaEE 构架师&#xff0c;目前立志于做一个与时俱进的独立的顾问架构师。在本篇文章中&#xff0c;作者通过多个方面为生产环节的日志提…
阅读更多...
Android, WindowsPhone7, IOS ,vc2010平台40多套图片滤镜开源

Android, WindowsPhone7, IOS ,vc2010平台40多套图片滤镜开源

历经三个多月开发&#xff0c;终于完成了这三个库的编码工作。用到的开发语言主要有使用了C&#xff08;含部分Object-c&#xff09;,JAVA,C#。经过检测&#xff0c;目前已正式投入到生产环境中。 为了实现三种语言一种命名方式&#xff0c;我采用了相同的类名&#xff0c;接口…
阅读更多...
找到一个比较熟悉Go语言的国内博客

找到一个比较熟悉Go语言的国内博客

2019独角兽企业重金招聘Python工程师标准>>> 他写了一些不错的Go安装和使用的文章&#xff0c;还参与了一本Go语言书籍的翻译工作&#xff0c;对Go语言相当熟悉。博客地址&#xff1a; http://www.mikespook.com/category/golang-2/ 转载于:https://my.oschina.net/…
阅读更多...
React开发（281）:一定是路由错误

React开发（281）:一定是路由错误

阅读更多...
React开发（282）:公共组件可以提升一下文件层级

React开发（282）:公共组件可以提升一下文件层级

阅读更多...
解决EditText不能撑满全屏的问题及EditText你应该知道的属性

解决EditText不能撑满全屏的问题及EditText你应该知道的属性

一般我们要实现去下图一的效果很简单&#xff1a;两个EditText就搞定效果图一&#xff1a;但是我们想让第二个EditText撑满剩余空间怎么做&#xff1f;如效果图二效果图二&#xff1a;解决&#xff1a;使用了ScrollView嵌套LinearLayout&#xff0c;将ScrollView中android:fill…
阅读更多...
Active Directory授权还原

Active Directory授权还原

Technorati 标签: Active Directory,授权还原额外域控制器有很多好处&#xff0c;例如可以平衡用户对AD的访问压力&#xff0c;有利于避免唯一的域控制器损坏所导致域的崩溃。域内所有的域控制器都有一个内容相同的Active Directory&#xff0c;而且Active Directory的内容是动…
阅读更多...
PHP的SESSION使用，文件上传下载等

PHP的SESSION使用，文件上传下载等

2019独角兽企业重金招聘Python工程师标准>>> PHP中Session的使用&#xff1a; 1. session_start(): 初始化session,如果页面要使用session&#xff0c;必须先使用这句话&#xff0c;否则会报错。 告诉服务器使用session。一般来说&#xff0c;php是不会主动使用sess…
阅读更多...
React开发（284）:注意点击事件大小写

React开发（284）:注意点击事件大小写

阅读更多...
[数学最安逸][UVa1638改编][第一类斯特林数+组合数]杆子的排列

[数学最安逸][UVa1638改编][第一类斯特林数+组合数]杆子的排列

有高为1,2,3,...,n的杆子各一根排成一行。从左边能看到l根&#xff0c;从右边能看到r根&#xff0c;求有多少种可能。 (l,r < 200,n < 200000) 给出T 组数据 (T < 500000) 对于每一组数据输出可能的个数&#xff0c;为避免写高精&#xff0c;将答案模 1e9 7 (它为质…
阅读更多...
Servlet的生命周期 与CGI的区别

Servlet的生命周期 与CGI的区别

2019独角兽企业重金招聘Python工程师标准>>> 对于Servlet&#xff0c;服务器只创建每个servlet的单一实例&#xff0c;每个用户请求都会引发新的线程-----将用户请求将会给相应的doGet/doPost进行处理。那么&#xff0c;servlet是如何进行创建和销毁的&#xff1f; …
阅读更多...
自动化打包资源混淆集成python实践----资源混淆

自动化打包资源混淆集成python实践----资源混淆

前面自动化打包资源混淆集成python实践----打包一文讲述了四种打包方案&#xff0c;以及美团打包方案、apk注释添加渠道号方案的实现。这里讲集成资源混淆。 1、资源混淆带来的好处&#xff1a; 1&#xff09;对资源文件起一定的保护作用&#xff0c;使其不能知其名不知意&…
阅读更多...
[win7] 去除将窗口拖到屏幕边缘时“自动最大化”

[win7] 去除将窗口拖到屏幕边缘时“自动最大化”

在win7中将一个窗口拖到屏幕的边缘&#xff0c;这个窗口将会自动以最大化方式显示&#xff0c;这对于那些需要同时查看两个窗口中内容的用户来说就是个恶梦&#xff01;因为屏幕大小是有限的&#xff0c;为了同时查看两个窗口中的内容&#xff0c;就要将这两个窗口以合适的方式…
阅读更多...
iOS 应用程序的国际化

iOS 应用程序的国际化

什么事应用程序的国际化? 就是不同国家之间使用的语言不一样,中国人的应用程序显示的汉字,美国人的应用程序是英文. 接下来我们就来实现应用程序的国际化: 1.首先我们新建一个工程,命名为:Internationalization-Demo; 2.添加语言(配置语言):(英文和中文) 3.目前&#xff0c;应…
阅读更多...
获得H.264视频分辨率的方法

获得H.264视频分辨率的方法

From: http://www.cnblogs.com/likwo/p/3531241.html 在使用ffmpeg解码播放TS流的时候&#xff08;例如之前写过的UDP组播流&#xff09;&#xff0c;在连接时往往需要耗费大量时间。经过debug发现是av_find_stream_info&#xff08;已抛弃&#xff0c;现在使用的是avformat_fi…
阅读更多...
Web服务器性能估算

Web服务器性能估算

1 【引题】 但凡写过技术方案的都知道&#xff0c;在技术方案最终落实到工程实施部署时&#xff0c;必须编制出当前解决方案需要部署的IT设备及环境&#xff0c;包括&#xff1a;需要的网络环境、端口、带宽、组网方式、网络安全保障措施&#xff1b;需配置的服务器设备性能…
阅读更多...
用实例分析H264 RTP payload

用实例分析H264 RTP payload

From: http://blog.csdn.net/zblue78/article/details/5948538 H264的RTP中有三种不同的基本负载&#xff08;Single NAL,Non-interleaved,Interleaved) 应用程序可以使用第一个字节来识别。 在SDP中也说明了本次会话的属性 SDP 参数 下面描述了如何在 SDP 中表示一个 H.264 流…
阅读更多...
实验4 颜色、字符串资源的使用

实验4 颜色、字符串资源的使用

课程名称 基于Android平台移动互联网开发 实验日期 3月25 实验项目名称 颜色、字符串资源的使用 实验地点 S3010 实验类型 □验证型 √设计型 □综合型 学 时 1学时 一、实验目的及要求&#xff08;本实验所涉及并要求掌握的知识点&#xff09; 掌握Androi…
阅读更多...
LAMP平台下构建Postfix邮件服务器

LAMP平台下构建Postfix邮件服务器

前言&#xff1a; 本人之前已经做过这个实验&#xff0c;只是版本不同&#xff0c;今天本来想用全新版本的源码包做实验&#xff0c;但是以httpd-2.4.2.tar.gz的源码包在配置虚拟主机这一块我不知道如何配置&#xff0c;按照以前的配置总是报错&#xff1a;AH00548: NameVirtua…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023