计算机仿真实验报告.pdfVIP

计算机仿真技术作业一

题目：转速反馈单闭环直流调速系统仿真

直流电机模型框图如下图所示，仿真参数为

R=0.6，T=0.00833，T=0.045，Ce=0.1925。本次仿真采用算法为ode45，仿真时间5s。

lm

I

d

uR1n

d01/R

Ts1TsCe

m

l

图1直流电机模型

一、开环仿真

模型建立：

其中Ud0=220VR=0.6，T=0.00833，T=0.045，Ce=0.1925

lm

其中0~2.5s，电机空载，即I=0；

d

2.5s~5s，电机满载，即I=55A

d

采用算法：

ODE45:

5秒

前2.5秒：

后2.5秒：

空载转速：1143rpm负载转速：971rpm

空载静差率s=0负载静差率s=0.1505

仿真时有波动

ODE23:

整体：

前2.5秒

后2.5秒

ODE23s:

整体：

前2.5秒

后2.5秒

ODE23t:

整体：

前2.5秒

后2.5秒

ODE23tb:

整体：

前2.5秒

后2.5秒

算法分析比较：从上可以看出ODE45与ODE23算法较差，仿真结果与理论不符合，电机

转速有纹波。ODE23sODE23tODE23tb效果较好，基本满足仿真需要，波形基本符合理论。

原因在于：ODE45、与ODE23都是一步解法，即只要知道前一时间点的解y(tn-1)，就

可以立即计算当前时间点的方程解y(tn)。后三种算法适用于刚性系统的解法，而前两种不

可。其中ODE23tb最适合电力电子系统仿真，它采用TR-BDF2算法，即在龙格.库塔法的

第一阶段用梯形法，第二阶段用二阶的BackwardDifferentiationFormulas算法。

2、闭环仿真

模型建立：

1、比例校正（K=5）

转速指令为1130rpm

为了便于比较不同k值时转速波形，简便框图

p

先进行模块封装：

1、确定输入输出与变量K

2.选中整个模块右键选择createsubsystem

即得：

3、右键选择createmask

得：

定义变量名K，并使变量K与封装模块中的变量K相互对应。

4、此时点击该模块可得：

即可更改变量K

然后将多个模块集中在一起显示：

仿真K=5，K=10，K=20三种，采用算法ODE23tb

得：

整体：（黄、红、绿分别代表K=5，K=10，K=20）

我们可以看出，在一定范围内，K越大静态误差就越小，但是比例环节无法消除静态误差。

前2.5秒

空载时：K=5，n=1088rpms=0.0372

K=10,n=1108.66rpms=0.0189

K=20,n=1119.23rpms=0.00953

可以看出比例环节无法消除稳态误差，但是适当增加比例系数大小却可以减小误差。

后2.5秒：

加入负荷时：

K=5n=1081.75rpms=0.0427

K=10n=1105.42rpms=0.0218

K=20n=1117.60rpms=0.