基于Simulink建模逻辑的双馈异步风力发电机(DFIG)数学仿真系统
项目介绍
本系统是一个基于MATLAB脚本编写的高集成度双馈异步风力发电机(DFIG)数学仿真环境。系统通过纯代码形式模拟了风力发电的核心物理过程,包括风能捕获子系统、传动链力学模型、双馈感应电机(WRIG)的dq轴坐标系数学模型,以及双PWM变流器的矢量控制逻辑。该仿真旨在研究风机在多变风速条件下的动态响应,重点展示了定子磁场定向控制(SFO)和最大功率点跟踪控制(MPPT)的实现机理,能够精确反映电磁转矩、功率流向及直流母线电压的稳态与动态特性。
功能特性
- 全流程数学建模:涵盖了从自然风速到叶片扫掠、齿轮箱升速、感应电机发电直至并网功率控制的全过程。
- 矢量控制实现:完整复现了机侧变流器(RSC)的转速-电流双闭环解耦控制和网侧变流器(GSC)的电压控制逻辑。
- 动态风速模拟:支持包括恒定风速、阶跃风速和渐变风速在内的复合环境输入,模拟风场真实工况。
- MPPT追踪能力:通过在线实时计算最优叶尖速比,动态调整电机转速以维持最高风能利用系数Cp。
- 电力电子闭环仿真:包含基于能量守恒原理的直流母线电压动态更新模型,体现了背靠背变流器的功率协调过程。
- 可视化分析:系统可自动生成风速、转速、有功/无功功率、转矩平衡、Cp曲线以及定子三相电流局部放大图。
系统要求
- 软件环境:MATLAB R2018b 或更高版本。
- 基础模块:MATLAB内置函数库。
- 硬件性能:由于采用数值迭代(欧拉法),普通办公笔记本电脑即可执行毫秒级步长的长周期仿真。
实现逻辑说明
仿真核心程序遵循严密的物理反馈回路,主要逻辑如下:
- 参数与状态初始化:程序开始时预设1.5MW等级发电机的定转子电阻、感抗、极对数及惯量,并初始化直流母线电压为1200V、发电机转速为0.8倍同步速。
- 仿真主循环:采用时间步长为0.2ms的循环进行迭代计算,总时长为5秒。
- 风力机物理模拟:通过非线性Cp近似公式计算风能利用系数。根据即时风速、叶片半径和转速得出机械转矩,作为发电机驱动源。
- 定子磁场定向(SFO):在dq旋转坐标系下,将定子磁链矢量定向于d轴,实现定子有功与无功的解耦。
- 控制策略执行:
- 机侧控制:转速环输出q轴转子电流指令以控制转矩(MPPT),d轴指令控制无功,通过PI调节器输出转子电压补偿项。
- 网侧控制:电压环维持中间直流母线恒定,d轴电流控制保证单位功率因数并网。
- 电机状态演化:基于欧拉法求解转子电流的微分方程,根据转子磁链与定子磁链的电磁耦合关系更新定子电流。
- 转轴动力学模型:通过电磁转矩与机械转矩的差值,更新发电机的角加速度、角速度及转子物理位置角度。
- 功率与母线更新:计算定子端输出的瞬时有功和无功,并根据机侧与网侧变流器的功率差额更新直流电容的电荷状态。
关键算法与实现细节分析
- PI调节器算法:
实现了一个带有输出限幅功能的通用比例积分控制器。该函数通过累加误差的历史值实现积分项,并结合当前误差的比例项输出控制量,有效地抑制了系统稳态误差,同时通过限幅处理防止控制量饱和带来的系统失稳。
- 坐标变换与简化:
在坐标系处理上,模型简化了定子电压方程,假设定子侧电网电压稳定且定向在q轴。通过3/2坐标变换逻辑,将定子侧的d、q轴电流映射回静止坐标系下的三相平衡电流,从而生成直观的ias、ibs、ics波形图。
- MPPT查找算法:
利用一个预设的最优叶尖速比函数(lambda_opt),在任意风速下计算出发电机应达到的目标转速。这一逻辑是机侧变流器速度环的核心参考源,确保了风机在非额定风速下始终运行在最大效率区间。
- 感应电机状态空间求解:
代码通过推导出的转子电压方程和磁链关系,建立了双馈电机的电流微分方程组。在每一采样瞬间,通过当前的电压输入和磁链状态计算电流导数,实现了对电机内部复杂电磁瞬变的数值仿真。
- 功率平衡逻辑:
直流母线电压的演化取决于网侧吸收功率与机侧转子交换功率的平衡。代码通过能量转换方程 $dvdc = Delta P / (Vdc cdot C)$ 模拟了电容充放电过程,体现了背靠背变流器之间电压支撑功能的硬连接特性。