简介一份基于Matlab/Simulink的垂直轴风力发电系统设计毕业论文PDF面向电气工程、新能源及自动化方向的本科或研究生可作为风力发电系统建模、仿真与控制策略研究的参考模板。压缩包内仅含1个PDF文件整体大小8.3MB内容为完整论文正文包含系统需求分析、垂直轴风机设计原理、控制策略、Simulink建模与仿真流程等章节阅读方便适合直接对照学习。已有127人学习下载属于较为专精的毕业设计类资料。读者既能学习垂直轴风力发电机组的气动与机械建模方法也能了解最大功率点追踪、桨距控制和偏航控制在Simulink中的实现思路并可借鉴其从需求分析、系统设计到仿真优化的论文组织方式对开展相关课题和撰写毕业论文均有直接帮助。1. 垂直轴风力发电系统仿真的切入点比想象中更靠近控制器做垂直轴风力发电系统的仿真很多人第一反应是先把风轮画出来或者急着找一台现成的风力机模型塞进 Simulink。实际上真正决定论文能不能往下写、仿真能不能收敛的是风力机的气动功率系数 Cp 和叶尖速比 TSR 之间的那条曲线以及后续变流器控制器的带宽设置。垂直轴风机和水平轴不一样它的叶片攻角在旋转周期内一直在变转矩脉动大启动特性差一点但也没有对风装置结构简单在城市分布式发电和岛屿微网里反而合适。做设计时要用 Simulink 搭一个完整的“风轮-永磁同步发电机-PWM 变流器-电网”闭环重点不是把每个器件画得像而是让气动输入、转速控制、电流环和直流母线电压控制互相匹配。这篇从建模说起最后一直走到 C 代码生成前的那一步中间涉及的参数都能直接抄进模型里改。2. 垂直轴风机的气动建模Cp-λ 曲线在 Simulink 里的三种落法2.1 垂直轴风轮的 Cp 特性为什么是控制系统设计的第一个难点水平轴风机在额定风速以下的控制目标很直观让叶尖速比 λ 保持在最优值附近逼近最大功率系数 Cpmax。垂直轴风机尤其是 H 型达里厄风机虽然同为升力型风机但叶轮旋转一圈时叶片在不同方位角下对来流攻角变化非常剧烈气动载荷不平稳所以 Cp-λ 曲线不像水平轴那么平滑最优叶尖速比附近有一个尖锐的峰偏离 0.5~1 个 λ功率系数就明显掉下来。这就意味着控制策略必须把转速精确锁在最优 TSR 附近。Simulink 里做这个环节惯性环节和采样时间都不能设得太粗糙否则两个控制周期内转速已经偏离最优工作点。做气动建模时我一般不会直接把某个现成的 Cp 函数当作黑盒塞进模型而是把 Cp 和 TSR 的映射关系拆成两种可复现的形式一种是用 MATLAB Function 块写解析公式一种是用 Lookup Table 直接查表。解析公式适合做参数化研究查表适合从 CFD 或风洞数据直接导入。最忌讳的是把 Cp 写成固定常数那样 MPPT 闭环就完全失去了意义。2.2 用 MATLAB Function 块写 Cp-TSR 关系的最小代码设 H 型达里厄风轮的半径 R扫掠面积 A 2RH风速 Vw风轮转速 ωr叶尖速比 λ ωr R / Vw。气动转矩 Ta 0.5 ρ A Vw^2 R Cp(λ) / λ。这里的 Cp(λ) 用升力线理论拟合形式表示常见的参数化写法如下function Cp cp_h_darrieus(lambda, sigma, alpha0) % 适用于中等实度 H 型达里厄风机的简化 Cp 模型 % lambda: 叶尖速比sigma: 实度alpha0: 参考攻角 if lambda 0 Cp 0; return; end % 用高斯形式近似达里厄风机的 Cp 峰 Cp 1.6 * (1 - 0.33 * sigma) * exp(-0.12 * (lambda - 2.8)^2) ... * (sin(alpha0 * pi / 180) 0.5 * cos(alpha0 * pi / 180)) / 1.5; % 限制物理上限 if Cp 0.45 Cp 0.45; end if Cp 0.02 Cp 0.02; end end这段代码放在 Simulink 的 MATLAB Function 块里输入是 λ输出是 Cp。实度 σ 取 0.25~0.35对应中等实度 H 型风机。需要说明的是这个解析式不是为了替代 CFD而是为了让模型在宽风速范围内有个合理的单调趋势。真正的风洞或 CFD 数据出来后这个函数可以直接替换成查表块。仿真时要注意MATLAB Function 块在定步长离散模型里输出是零阶保持的所以风速突变时气动转矩的更新会有一步延迟这对 MPPT 控制器的响应分析是有影响的。如果做快速暂态分析要改用 S-Function 并在内部维护连续状态。2.3 算例参数半径 2.5m、高度 5m 的风轮怎么配转速范围给定设计风速 Vw_design 10 m/s最优叶尖速比 λ_opt 2.8那么风轮最优角速度 ω_opt λ_opt Vw_design / R 2.8×10 / 2.5 11.2 rad/s换算成转速约 107 rpm。扫掠面积 A 2RH 2×2.5×5 25 平方米在空气密度 1.225 kg/m³ 下理论最大功率 Pmax 0.5ρAVw_design^3Cp_opt取 Cp_opt 0.38得到 Pmax 约 5.8 kW。这个功率级对应的小型直驱永磁同步发电机额定转速通常在 100~200 rpm属于低速大扭矩工况反电动势不高所以变流器的直流母线电压可以设计在 400V 左右而不用升压变压器。仿真模型里风轮转速输入用积分环节把气动转矩和电磁转矩差积分得到转动惯量 J 需要根据叶片结构和发电机转子估算一般取 2~6 kg·m²。转动惯量设置偏小转速环的响应会显得过快MPPT 会产生高频抖动偏大在阶跃风下转速跟踪变慢动态能量损失会很大。开始时可以用 J 3.5 kg·m² 跑通闭环再微调。2.4 三种气动建模方案选型对比方案实现方式参数可调性精度适用阶段解析公式法MATLAB Function 写 Cp(λ) 表达式高可直接调 Cp 峰值中趋势正确但峰型偏理想毕业论文前期、控制策略对比查表法Lookup Table 导入 CFD/实验数据中需要插值方式匹配高接近真实风轮设计定型后的精确校核S-Function 法在 C 或 MATLAB 中写完整气动模型高可加动态尾流高但开发量大研究动态失速、尾流影响2.4.1 查表法的一个关键细节用 Lookup Table 时断点数据必须是单调递增的。如果你的 CFD 数据在 λ 接近 0 时出现数值振荡导致 Cp 曲线不是严格单调的Simulink 会报错或者插值结果异常。处理办法是先对原始数据做滑动平均再保证第一个断点从 λ0.5 开始不要放 0。因为在 λ0 附近数值上 Cp 的定义本身已经病态。3. 最大功率追踪 MPPT最佳叶尖速比法在 Simulink 里的闭环实现3.1 转速外环需要的不是最大功率本身而是最优转速信号额定风速以下垂直轴风机的控制目标就是追踪最优叶尖速比。最佳叶尖速比法TSR 法的原理是实时测量风速计算当前风速下的最优转速 ω_ref λ_opt Vw / R然后通过转速闭环让风轮转速逼近这个参考值。这个方法的关键在于风速测量的准确性。小型垂直轴风机安装在城市环境湍流强度大机舱风速计受叶片尾流干扰后测量偏差可能超过 20%。所以工程上通常用风速估计即通过当前风轮转矩和转速反推等效风速而不是直接使用风速计读数。3.2 在 Simulink 中用 S-Function 实现 TSR 法 MPPT下面给出一个简单的 Level-2 MATLAB S-Function输入为转速和风速输出为转矩参考值。这个实现比 Simulink 内置的 MPPT 模块更容易解释给审稿人看因为表达式是显式的。function tsrMppt(block) setup(block); end function setup(block) block.NumInputPorts 2; block.NumOutputPorts 1; block.InputPort(1).Dimensions 1; % 转速 block.InputPort(2).Dimensions 1; % 风速 block.OutputPort(1).Dimensions 1; % 转矩参考 block.SampleTimes [0 1e-3]; % 固定步长 1ms block.RegBlockMethod(Outputs, Output); end function Output(block) w block.InputPort(1).Data; % 转速 rad/s V block.InputPort(2).Data; % 风速 m/s R 2.5; lambda_opt 2.8; Cp_opt 0.38; rho 1.225; A 25; w_ref lambda_opt * V / R; if w 1e-3 Te_ref 0.5 * rho * A * V^3 * Cp_opt / (w * lambda_opt); else Te_ref 0; end block.OutputPort(1).Data Te_ref; end这个 S-Function 相当于一个前馈转矩计算器。它没有用 PI 来跟踪转速而是直接用风速算出期望电磁转矩近似为力矩控制。实际模型中转速环 PI 仍然需要保留因为转矩前馈并不能保证转速对参考值的跟踪尤其在风机和发电机转动惯量大而阻尼小的时候。常见做法是把这个转矩前馈作为前馈通道加到转速 PI 的输出上这样既有快速响应又有稳态精度。参数方面的说明采样时间 1ms 对 MPPT 环足够了风电系统的机械时间常数通常在百毫秒量级。S-Function 里直接用大常数表达式没有做标幺化这在小功率模型里没问题但如果把模型扩展到兆瓦级建议改用标幺制否则数值范围跨度过大会影响定步长求解器误差。3.3 用阶跃风速验证 MPPT 动态性能的一组波形标准仿真里把风速从 8 m/s 阶跃到 10 m/s观察下面几个量是否达标风轮转速能否在 2~3 秒内过渡到新的最优转速超调量不超过 10%电磁转矩曲线没有明显抖振尤其在阶跃后的第一个周期内发电机输出功率在过渡后保持在理论最大功率的 95% 以上3.4 和功率信号反馈法以及爬山法的区别MPPT 方法需要传感器动态响应稳态振荡适用场景最佳叶尖速比法风速、转速快无风速平稳的并网型系统功率信号反馈法转速、功率较快无大型风机、无风速传感器爬山法转速、功率慢需扰动有小型独立系统、湍流强烈爬山法在垂直轴风机上有个额外问题扰动步长必须和风轮转动惯量匹配。步长太小收敛慢太大则在 Cp 曲线峰值附近来回振荡输出功率纹波变大。有一个经验值扰动步长设为额定转矩的 2%~5%扰动周期设为机械时间常数的 3~5 倍这样基本不会发散。4. PMSG 与背靠背变流器dq 解耦控制从模块接线到 PI 整定4.1 直驱垂直轴系统为什么选 PMSG以及它带来的控制频率问题垂直轴风机的转速通常在 100~300 rpm 范围属于低速直驱场景所以选择永磁同步发电机 PMSG 比异步电机更合适因为不需要齿轮箱没有励磁电流损耗也没有滑环维护问题。但 PMSG 的多极对数大所以电频率不是风轮转速本身而是极对数倍的机械频率。比如极对数 p20风轮转速 100 rpm对应机械频率约 1.67 Hz电频率就是 33.3 Hz。这个频率低于工频 50 Hz所以必须通过背靠背变流器先整流再逆变。如果电机极对数更多频率会高一些但电流环的采样频率总是要远高于电频率。Simulink 里用 Simscape Electrical 的 Permanent Magnet Synchronous Machine 模块参数表中要同时配置定子电阻 Rs、电感 Ld/Lq、磁链 psi、以及极对数 p。其中的定子电感参数很多人直接拿电机设计书的线电感填入但 Simulink 模型里的 Ld 和 Lq 是 dq 坐标系下的电感。如果电机是表贴式LdLq如果是内置式LdLq这时 MPPT 和弱磁控制都要分开设计。对于垂直轴直驱应用表贴式转子结构更常见因为不需要额外的磁阻转矩控制简单。4.2 机侧变流器转速环与电流环配合功率能稳住的接线顺序机侧变流器的主流控制策略是零 d 轴电流控制也就是 id_ref0用 q 轴电流控制电磁转矩。整个控制环的接线顺序从外到内是转速环输出 q 轴电流参考值电流环输出 d/q 轴电压参考值再经过反 Park 变换得到静止坐标系电压进入 PWM 发生器。Simulink 具体搭建时用 Demux 把测量到的三相电流拆开经过 Clarke 变换得到 α/β 轴分量再经过 Park 变换得到 d/q 轴分量。这一步角度用 PMSG 内部输出的转子角度不要自己用积分器积分风速或者转速来凑角度否则带负载后角度相位差会直接击穿电流环。机侧变流器通常用 SVPWM因为直流母线电压利用率比 SPWM 高约 15%。Simulink 里可以用 Simscape Electrical 自带的 Universal Bridge 模块选择 SVPWM也可以直接用 Mean Value Inverter 模块加速仿真。注意Mastering 论文时不要为了省事把逆变器替换成受控电压源这样电流谐波分析、直流母线电容纹波分析都做不真实答辩被问一句就漏了。4.3 网侧变流器直流母线电压恒定和单位功率因数必须同时做到网侧变流器的外环是直流母线电压环内环是网侧电流环。直流母线电压给定通常设在 400V电容容值根据允许的电压纹波和母线功率级估算。另一个重要控制目标是单位功率因数也就是网侧 q 轴电流参考值为 0。但如果想为电网提供无功支持可以让 q 轴电流参考不为 0。在论文仿真里先把 q 轴电流设为 0验证系统有功功率平衡后再尝试无功调节。网侧电流环的结构和机侧类似但需要锁相环 PLL 看电网电压相位。Simulink 中三相 PLL 模块输出的角度用于 Park 变换。用 Three-Phase Programmable Voltage Source 做电网时仿真步长要控制在 50μs 以内因为 PWM 载波频率通常在 10kHz 左右。步长太大会出现谐波混叠直流母线电容电压的计算也不准确。4.4 带宽法整定 PI 参数从电机参数到控制器参数的直接换算电流环 PI 参数用电流环带宽和电机电感、电阻换算工程公式如下Ki_rate R_s / (L_s * tau_i) Kp L_s * bw_i Ki Kp / tau_i这里 bw_i 是电流环带宽单位为 rad/s通常取电频率的 10~20 倍tau_i 是电流环期望时间常数取 1/bw_i。转速环的带宽要比电流环低 5~10 倍否则两个控制环互相激励。参数典型值说明电流环带宽 bw_i62.8~125.6 rad/s对应 10~20Hz至少 10 倍于电频率电流环 KpL×bw_i保证电流环不出现相位裕度不足电流环 KiKp×bw_i消除稳态误差的积分作用转速环带宽 bw_w6~12 rad/s低于电流环 5~10 倍转速环 KpJ×bw_w / K_tK_t 为转矩常数需查 PMSG 数据表整定完要在 Simulink 里用阶跃输入验证给 d 轴电流一个 0 参考q 轴电流一个阶跃给定观察电流实际值上升时间是否在 2~3 倍时间常数以内是否有超调。用线性分析器或 Simulink Control Design 提取电流环的闭环 bode 图相位裕度应大于 45 度。折腾完以上环节整个 Simulink 模型才算真正完整风力机气动部分、MPPT 部分、PMSG、背靠背变流器、电流环、转速环、直流母线电压环。下一步要考虑的就是把这个模型搬到嵌入式处理器之外去验证也就是现代控制开发流程中的模型到 C 代码生成。5. Simulink 模型走到 C 代码生成前必做的配置检查清单5.1 模型的求解器、步长和数据转换三处配置建议能跑通的仿真模型和能生成可执行 C 代码的模型差别很大。做 C 代码生成前Simulink 模型必须满足下面这组约束求解器必须是定步长离散常见用 ode4 或 ode3步长建议取 1e-3 秒。模型中不能有连续状态块和变步长求解器依赖所有连续模块要用离散近似替换。例如 PMSM 模块在 Simscape Electrical 中需要换成 Simscape Electrical 的离散化版本否则生成代码时直接报错。数据转换必须显式做所有 double 类型的信号要转换到 single 或 int16 等目标数据类型不然在嵌入式设备上运行效率非常低。Simulink 的 Data Type Conversion 块放在每个接口处用硬件实现的方式是 Embedded Coder 里配置硬件信息后再做自动类型转换。5.2 用一段 MATLAB 脚本批量修改模型配置逐个点菜单改模型配置容易漏掉某个关键选项。下面这段脚本可以直接把模型切换到适合代码生成的状态function configModelForCGen(modelName) % 将模型配置从仿真模式切换到适合代码生成的状态 load_system(modelName); cfg getActiveConfigSet(modelName); % 求解器配置定步长离散ode4 set_param(cfg, SolverType, Fixed-step); set_param(cfg, Solver, ode4); set_param(cfg, FixedStep, 1e-3); % 应用类型设置为基础代码生成 set_param(cfg, SystemTargetFile, ert.tlc); set_param(cfg, GenerateMakefile, off); set_param(cfg, GenCodeOnly, on); % 硬件信息假设目标为 ARM Cortex-M 系列 set_param(cfg, ProdHWDeviceType, ARM Compatible-ARM Cortex-M); % 代码生成时保留模型注释方便反查 set_param(cfg, GenerateComments, on); set_param(cfg, GenerateTraceInfo, on); save_system(modelName); disp(模型配置已更新可尝试生成代码); end脚本中使用 ert.tlc 表示生成嵌入式实时目标代码。如果不做硬件在环只做处理器在环测试可以先将 GenCodeOnly 保持为 on防止 Simulink 生成多余的文件先跑一个“干净”的 C 代码出来再手动拉进编译环境。注意 ert.tlc 这个目标不能在普通 Simulink 里直接选需要安装 Embedded Coder。没有这个工具箱时只能选择 slrt.tlc 或 grt.tlc但这样代码的可移植性差很多。5.3 外部模式验证从数据观测器到波形记录的一步到位生成的代码烧进实机之后怎么验证代码和模型行为一致Simulink 外部模式 External Mode 可以把目标机上的实时信号传回宿主机显示。具体操作是在模型配置中勾选 External mode设置传输接口为串口或 TCP/IP然后在代码生成选项中勾选 Generate code only 之外的那项 Genesis 接口代码。启动外部模式后在 Simulink 的 Scope 里观测到的数据和模型仿真数据完全是一套汉化界面但数据结构已经被简化成适合实时运行的版本。外部模式有一个限制它只能通过通信接口传输信号值不能实时在线调参。如果需要在实机上改控制参数建议加上校准接口用生成的 ASAP2 文件配合 INCA 或 CANape 工具做在线标定。论文阶段通常不需要走那么深只需证明生成的模型在外部模式下能跑、波形符合预期即可。一个很实用的验证技巧把风速阶跃和 MPPT 的测试序列预先存在目标机的数据采集卡里利用 Simulink 的 Data Store Memory 模块在外部模式下实现测试场景的自动切换。这样反复跑动态工况不需要每次重新编译下载可以看数据波形是否有不规则尖峰。这一步做完并且波形数据和模型仿真基本吻合基本可以说这个基于 Matlab/Simulink 的垂直轴风力发电系统设计已经从纸面理论走到可运行的设备控制阶段了。本文还有配套的精品资源点击获取