《机电传动控制》第四周作业

2023-11-02 21:59

本文主要是介绍《机电传动控制》第四周作业,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、题目要求:

结合本周学习的交流电机原理及启动、调速、制动特性,用Modelica设计和仿真一个用三相交流异步电机带动起重机起升机构运行。具体要求如下:

1)实现如下机械运动周期:

控制电机带重物上升,从静止加速到800r/min

保持800r/min匀速运动0.5s,

减速到静止,保持静止状态0.5s,

带重物下降,从静止达到600r/min

保持600r/min匀速运动0.6s,

减速到静止。
(为了便于仿真,匀速和静止持续时间较短)

2) 升降机构和重物折算到到电机转子轴上的等效负载惯量为1Kg.m^2,折算到到电机转子轴上的等效负载转矩是15N.m。

3)使用统一的电机模型,如果控制策略中用到转子串电阻,允许将该电机的转子改为绕线式转子(参数不变)。

4)参照教材中给出的交流电机启动、调速和制动方法,设计控制策略,用Modelica实现控制策略并与电机模型实现联合仿真。

5)可以采用定子串电阻、转子串电阻、定子调压、定子调频等手段,但必须具备工程上的可实施性。

6)评价指标:快速启动、制动,冲击转矩和冲击电流小,能耗小,兼顾实施的经济性。

二、仿真分析:

由于转速n与fs存在正比例的关系,所以采用改变fs来控制n的方法。在老师所给程序案例中,fs为50Hz时对应电机转速n为1500r/min,所以认为n=30fs。以下程序编写,均以控制不同时间段的fs为基础,来完成相应的转速要求,具体内容参考郑君政同学的参数设置。

三、仿真代码及结果:

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";

 

parameter Resistance Rs = 0.531+0.5 "Stator Resistance";

parameter Resistance Rr = 0.408+0.5 "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 = 1"Load Inertia";

parameter Frequency f1 = 27;

parameter Frequency f2 = 15;

parameter Frequency f3 = 0.2;

parameter Frequency f4 = 19.8;

parameter Frequency f5 = 19.8;

parameter Real t1 = 100+500+671; //加速 800r/min恒速

parameter Real t2 = t1+50; //减速

parameter Real t3 = t2+600; //静止0.5s

parameter Real t4 = t3+200; //反向加速

parameter Real t5 = t4+300+600; //600r/min恒速

 

 

Real time1(start=0);

Real time2(start=0);

Real time3(start=0);

 

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);

 

0 = Rr * i_a + 1000 * der(Psi_a);

0 = Rr * i_b + 1000 * der(Psi_b);

0 = 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);

 

if time <= 100 then

u_A = 0;

u_B = 0;

u_C = 0;

Tl = 0;

elseif time <= t4 then

u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000);

u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000-2*Pi/3);

u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000-4*Pi/3);

Tl = 15;

else

u_A = u_N * 1.414 * sin(2*Pi*f_s*time/1000);

u_B = u_N * 1.414 * sin(2*Pi*f_s*time/1000+2*Pi/3);

u_C = u_N * 1.414 * sin(2*Pi*f_s*time/1000+4*Pi/3);

Tl = 15;

end if;

 

 

algorithm

 

if time <= 100 then

f_s := 0;

elseif time <= t1 then

f_s := f1;

elseif time <= t2 then

f_s := f2;

 

elseif time <= t3 then

f_s := f3;

elseif time <= t4 then

 

f_s := (time-t3)/(t4-t3)*(f4-f3)+f3;

 

elseif time <= t5 then

f_s := f5;

else

f_s := 0;

end if;

 

if n > 800*0.98 and n < 800*0.99 then

time1 := time;

end if;

if n < 0.1 then

time2 := time;

end if;

if n > -600*0.98 and n < -600*0.99 then

time3 := time;

end if;

end SACIM;

四、总结

基本按照要求完成任务,曲线整体较为平滑。

转载于:https://www.cnblogs.com/yujin13/p/5308803.html

这篇关于《机电传动控制》第四周作业的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/333601

相关文章

Python异步编程中asyncio.gather的并发控制详解

《Python异步编程中asyncio.gather的并发控制详解》在Python异步编程生态中,asyncio.gather是并发任务调度的核心工具,本文将通过实际场景和代码示例,展示如何结合信号量... 目录一、asyncio.gather的原始行为解析二、信号量控制法:给并发装上"节流阀"三、进阶控制

使用DrissionPage控制360浏览器的完美解决方案

《使用DrissionPage控制360浏览器的完美解决方案》在网页自动化领域,经常遇到需要保持登录状态、保留Cookie等场景,今天要分享的方案可以完美解决这个问题:使用DrissionPage直接... 目录完整代码引言为什么要使用已有用户数据?核心代码实现1. 导入必要模块2. 关键配置(重点!)3.

SpringSecurity 认证、注销、权限控制功能(注销、记住密码、自定义登入页)

《SpringSecurity认证、注销、权限控制功能(注销、记住密码、自定义登入页)》SpringSecurity是一个强大的Java框架,用于保护应用程序的安全性,它提供了一套全面的安全解决方案... 目录简介认识Spring Security“认证”(Authentication)“授权” (Auth

python之流程控制语句match-case详解

《python之流程控制语句match-case详解》:本文主要介绍python之流程控制语句match-case使用,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐... 目录match-case 语法详解与实战一、基础值匹配(类似 switch-case)二、数据结构解构匹

Spring Security注解方式权限控制过程

《SpringSecurity注解方式权限控制过程》:本文主要介绍SpringSecurity注解方式权限控制过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、摘要二、实现步骤2.1 在配置类中添加权限注解的支持2.2 创建Controller类2.3 Us

Python中如何控制小数点精度与对齐方式

《Python中如何控制小数点精度与对齐方式》在Python编程中,数据输出格式化是一个常见的需求,尤其是在涉及到小数点精度和对齐方式时,下面小编就来为大家介绍一下如何在Python中实现这些功能吧... 目录一、控制小数点精度1. 使用 round() 函数2. 使用字符串格式化二、控制对齐方式1. 使用

Springboot控制反转与Bean对象的方法

《Springboot控制反转与Bean对象的方法》文章介绍了SpringBoot中的控制反转(IoC)概念,描述了IoC容器如何管理Bean的生命周期和依赖关系,它详细讲解了Bean的注册过程,包括... 目录1 控制反转1.1 什么是控制反转1.2 SpringBoot中的控制反转2 Ioc容器对Bea

浅析如何使用Swagger生成带权限控制的API文档

《浅析如何使用Swagger生成带权限控制的API文档》当涉及到权限控制时,如何生成既安全又详细的API文档就成了一个关键问题,所以这篇文章小编就来和大家好好聊聊如何用Swagger来生成带有... 目录准备工作配置 Swagger权限控制给 API 加上权限注解查看文档注意事项在咱们的开发工作里,API

Spring IOC控制反转的实现解析

《SpringIOC控制反转的实现解析》:本文主要介绍SpringIOC控制反转的实现,IOC是Spring的核心思想之一,它通过将对象的创建、依赖注入和生命周期管理交给容器来实现解耦,使开发者... 目录1. IOC的基本概念1.1 什么是IOC1.2 IOC与DI的关系2. IOC的设计目标3. IOC

Python实现局域网远程控制电脑

《Python实现局域网远程控制电脑》这篇文章主要为大家详细介绍了如何利用Python编写一个工具,可以实现远程控制局域网电脑关机,重启,注销等功能,感兴趣的小伙伴可以参考一下... 目录1.简介2. 运行效果3. 1.0版本相关源码服务端server.py客户端client.py4. 2.0版本相关源码1