磁浮列车动态仿真

获奖团队介绍

学校:北京交通大学
指导老师:杜欣
团队成员:王战阳、陶浩鹏、潘雨冉、张雨笑

简介

磁浮列车实际上是依靠电磁吸力或电动斥力将列车悬浮于空中并进行导向,实现列车与地面轨道间的无机械接触,然后利用线性电机驱动列车运行。本案例介绍磁悬浮列车静态悬浮控制仿真分析,分别通过电流控制系统和气隙控制系统实现列车静态悬浮。
电磁悬浮就是对车载的、置于导轨下方的悬浮电磁铁通电励磁而产生电场,磁铁与轨道上的铁磁构件相互吸引,将列车向上吸起悬浮于轨道上,磁铁和铁磁轨道之间的悬浮间隙一般约为8—12mm。列车通过控制悬浮磁铁的励磁电流来保证稳定的悬浮间隙,通过直线电机来牵引列车前进。

使用说明

1.在打开主要模型前,首先需要在 Syslab 加载对应的.jl文件,进行路径规划,文件为 maglev.jl
代码如下:

using TyControlSystems
using TyMath
using TyPlot
kc = -3005915.12     # 位移刚度系数
ki = 801.577   # 电流刚度系数
M = 1226.9   # 单电磁铁(模拟磁浮列车)质量1226.9Kg
# 磁悬浮系统状态空间矩阵定义
A = [0 1
     -kc/M 0]
B = [0 -ki/M]'
C = [1 0]
D = 0
# 磁悬浮系统模型
G = ss(A, B, C, D)
M1 = ctrb(A, B) # 计算系统的能控性判别矩阵
N = obsv(A, C) # 计算系统的能观性判别矩阵
n_M1 = rank(M1) # 计算系统能控性判别矩阵的秩
n_N = rank(N) # 计算系统能观性判别矩阵的秩
# 系统重构
Ahat = [0 1 0
        0 0 1
        0 -kc/M 0]
Bhat = [0
        B]
# 计算理想主导极点位置并分配系统所有极点位置
AllPoles = [-42.9+34.4im, -42.9-34.4im, -500]
Khat = place(Ahat, Bhat, AllPoles)
# 状态增益值分解
K = reshape(Khat[1][1:3], 1, 3)

2.在Sysplorer加载 maglevmath.mo 和 Default 文件中的 package.mo 文件,模型如下图

该仿真模型能够在 0~0.2s 内,通过控制系统让磁浮列车从 16 mm 的初始位置运动到 8 mm 的平衡位置,同时在 3s 受到干扰,控制系统驱动电流变大,保证列车位置不变。

仿真结果分析

1.仿真在设定输入曲线和干扰的情况下,悬浮对象的位置从 16 mm 在 0.2s 内快速悬浮到了 8 mm,稳态误差为 0,满足指标要求,3s 前电磁铁实际电流为 30 A,3s 时由于磁浮系统受到方向向下的扰动电流大小由 30 A 上升到 183 A 以增加电磁吸力维持磁浮列车气隙平衡。

  • 输入平衡位置和干扰信号 fd 曲线:

  • 励磁线圈电流和气隙、电流双闭环系统位置响应曲线:

  • 0-0.6 秒位置响应曲线:

2.在 Sysplorer 仿真界面,在仿真界面选择图表>动画,绘制磁悬浮列车和轨道的三维动画,单击回放控制>播放按钮,实现动画播放。

版本说明

V0.0.1,2025-09-03 14:55

  • 初始版本

使用许可

本模型库版权由北京交通大学蓝桥杯参赛团队版权所有,未经原创作者许可,不得用于商业用途。