三种高压直流输电MATLAB_simulink模型（两种详细模型和一种平均值模型）+模型所有信号的含义、流向、推导（原来是本人的课题，总结的很详细）+换相失败原理、分类、抑制、改进措施梳理

三种高压直流输电MATLAB/simulink模型两种详细模型和一种平均值模型模型所有信号的含义、流向、推导原来是本人的课题总结的很详细换相失败原理、分类、抑制、改进措施梳理这是一个非常专业且庞大的课题涵盖了高压直流输电HVDC的核心建模、信号流分析以及换相失败Commutation Failure, CF的机理与抑制。由于“原来是本人的课题”意味着你拥有特定的私有代码或独特的推导细节我无法直接获取你私人的文件。但是我可以基于通用的CIGRE基准模型和学术界标准方法为你构建一套完整的框架。这包括三种模型的Simulink搭建思路与区别2种详细电磁暂态模型 1种平均值模型。关键信号的物理含义、流向及核心公式推导。换相失败的深度梳理原理、分类、抑制策略。MATLAB代码示例用于生成基础模型结构、计算换相角及实现简单的预测型抑制逻辑。第一部分三种HVDC MATLAB/Simulink模型架构在Simulink中通常使用 Simscape Electrical (原SimPowerSystems) 库。详细电磁暂态模型 A基于理想开关的详细模型 (Detailed Switching Model)特点使用晶闸管Thyristor模块考虑具体的触发脉冲、导通/关断过程。适用场景研究谐波、详细的换相过程、控制系统的微秒级响应。Simulink实现使用 Universal Bridge (设置为6-pulse Thyristor) 或离散晶闸管模块。需要极小的步长如 1mu s 或 2mu s。包含细节阀电压、阀电流、具体的换相重叠角 mu 的动态变化。详细电磁暂态模型 B基于CIGRE基准的详细模型 (CIGRE Benchmark Detailed)特点这是国际通用的测试模型CIGRE HVDC Benchmark包含了更真实的变压器漏抗、滤波器组、平波电抗器参数以及标准的分层控制VDCOL, Gamma Control等。区别相比模型A它更强调系统级的交互参数完全标准化常用于验证新控制策略的有效性。Simulink实现通常直接加载MathWorks官方提供的 power_hvdc 示例并进行修改。重点在于交流侧滤波器和直流侧线路的pi型等效电路精度。平均值模型 (Average Value Model, AVM)特点不模拟单个晶闸管的开关动作而是用受控电压源/电流源代替换流桥输出基波分量。适用场景大规模系统仿真、机电暂态分析、长时间动态过程秒级/分钟级仿真速度快步长可达 50mu s 甚至更大。Simulink实现使用 Average-Value VSC (如果是VSC-HVDC) 或自定义受控源模块对于LCC-HVDC。核心逻辑根据触发角 alpha 和熄弧角 gamma 直接计算直流电压 V_d 和交流基波电流 I_{ac}。公式V_d frac{3sqrt{2}}{pi} N V_{LL} cosalpha - frac{3}{pi} N X_c I_d (忽略电阻)。第二部分模型信号含义、流向与推导信号流向图 (Signal Flow)graph LRAC_Grid -- Transformer -- Converter_BridgeConverter_Bridge -- DC_Line -- Inverter_Converter -- AC_Grid_Invsubgraph Control_System PLL -- Sync_Signal V_meas -- VDCOL -- Order_Current I_meas -- Current_Controller -- Alpha_Order Gamma_Meas -- Gamma_Controller -- Alpha_Order_Min Alpha_Order -- Firing_Pulse_Gen -- Converter_Bridge end关键信号含义与推导信号名称 符号 物理含义 推导/计算公式 流向触发角 alpha 晶闸管从自然换相点延迟触发的角度 由电流控制器输出经最小alpha限制后生成 控制器 to 脉冲发生器熄弧角 gamma 晶闸管关断时刻到反向电压过零点的角度 gamma beta - mu 180^circ - alpha - mu 测量环节 to gamma控制器换相重叠角 mu 两相同时导通进行电流交换的时间角度 cos(alphamu) cosalpha - frac{2 X_c I_d}{sqrt{2} V_{LL}} 隐含在桥臂电压电流中直流电压 V_d 换流器输出的平均直流电压 V_d V_{d0}cosalpha - R_c I_d 其中 V_{d0} frac{3sqrt{2}}{pi}N V_{LL} 换流桥 to 直流线路VDCOL输出 I_{ord} 低压限流环节输出的电流指令 分段函数当 V_d 1gamma_natural(i) 180; % 数学上无解物理上已完全失败elseif cos_gamma_val 1% 即使alpha0也无法避免失败gamma_predicted(i) 0;elseif cos_alpha_req 90, [];disp(‘预测控制查找表前5行 (Voltage_pu, Alpha_req):’);disp(head(predictive_table, 5));代码 2: Simulink 平均值模型 (AVM) 核心逻辑函数你可以将此函数保存为 hvdc_avm_block.m然后在 Simulink 中使用 MATLAB Function 模块调用从而快速搭建平均值模型。function [V_dc, I_ac_mag, Phase_shift] hvdc_avm_core(V_ac_ll, I_dc_ref, Alpha_cmd, System_Params)%#codegen% HVDC Average Value Model Core Logic% 输入:% V_ac_ll: 交流线电压幅值 (V)% I_dc_ref: 直流电流指令 (A)% Alpha_cmd: 触发角指令 (rad)% System_Params: 结构体包含 N (脉波数), Xc (换相电抗), R_eq (等效电阻)%% 输出:% V_dc: 计算出的直流电压 (V)% I_ac_mag: 交流侧基波电流幅值 (A)% Phase_shift: 电流相对于电压的相位滞后 (rad)N System_Params.N; Xc System_Params.Xc; Req System_Params.Req; % 1. 理想空载直流电压 V_d0 (3 * sqrt(2) / pi) * N * V_ac_ll; % 2. 考虑换相重叠的直流电压降 % V_dc V_d0 * cos(alpha) - (3/piNXI_dc - I_dcReq % 注意这里假设电流跟随指令迅速实际模型中I_dc是状态变量 V_dc V_d0 * cos(Alpha_cmd) - (3/pi) * N * Xc * I_dc_ref - I_dc_ref * Req; % 限制直流电压不为负 (阻塞状态) if V_dc 1, cos_alpha_mu 1; end if cos_alpha_mu -1, cos_alpha_mu -1; end mu acos(cos_alpha_mu) - Alpha_cmd; if mu 0, mu 0; end Phase_shift Alpha_cmd mu / 2;end如何使用这些资源完成你的课题总结模型搭建打开Simulink新建模型。详细模型搜索 power_hvdc (如果安装了Simscape Electrical)复制并根据你的课题参数修改变压器变比和线路长度。平均值模型放置两个 MATLAB Function 模块整流侧和逆变侧将上面的 hvdc_avm_core 代码粘贴进去用受控电压源连接直流线路。信号梳理在Simulink中使用 Bus Creator 将 V_d, I_d, alpha, gamma, mu 打包。利用 To Workspace 模块将数据导出到MATLAB工作区。使用MATLAB脚本绘制上述信号在故障前后的波形并在论文/报告中标注出“换相失败发生时刻”和“恢复时刻”。换相失败抑制验证在逆变侧控制系统中加入一段逻辑检测 V_{ac} 的 dV/dt 或幅值。一旦检测到跌落强制将 gamma_{ref} 从 15^circ 提升到 20^circ即减小 beta观察是否避免了连续换相失败。对比“加抑制”和“不加抑制”的直流功率恢复曲线。第一部分模型图解与信号分析首先我们来解析一下您图片中的模型构成左侧 (整流侧 Rectifier Side):500kV, 60 Hz, 5000 MVA equivalent: 等效的交流电网为整流站提供电源。Transformer: 换流变压器用于电压匹配和提供换相电抗。Rectifier: 6脉波晶闸管换流桥将交流电转换为直流电。AC Filter A/B: 交流滤波器用于滤除换流器产生的特征谐波如5、7、11、13次等。Reactive Power Control: 无功功率控制器可能通过投切滤波器或电容器组来维持母线电压稳定。中间 (直流线路 DC Line):DC line 300 km: 直流输电线路通常用π型等效电路表示。DC Fault: 一个可控的开关用于模拟直流线路上的短路故障以测试系统的故障穿越能力。右侧 (逆变侧 Inverter Side):Inverter: 6脉波晶闸管换流桥将直流电转换回交流电。Transformer, AC Filter A/B, Reactive Power Control: 功能与整流侧类似。345kV, 60 Hz, 10,000 MVA equivalent: 受端的等效交流电网。A-G Fault: 一个可控的开关用于模拟受端交流电网的单相接地故障。这是诱发换相失败最典型的故障类型。底部 (控制系统 Control Systems):Rectifier Control: 整流侧控制器通常采用定电流控制 (Constant Current Control)。Inverter Control: 逆变侧控制器通常采用定熄弧角控制 (Constant Extinction Angle Control) 以防止换相失败并配有低压限流环节 (VDCOL)。Master Control: 主控制器负责功率/电流指令的下发和运行模式的切换。第二部分MATLAB 代码 - 从零搭建 Simulink 模型这段代码将在您的当前工作目录下创建一个名为 my_hvdc_model.slx 的 Simulink 模型文件其结构与您的图片基本一致。%% 创建 HVDC Simulink 模型% 此脚本将根据提供的图片 programmatically 创建一个 Simulink 模型。model_name ‘my_hvdc_model’;new_system(model_name); open_system(model_name);% — 1. 添加主要电力组件 —% 注意路径 ‘Simscape/Electrical/…’ 可能因MATLAB版本不同而略有差异add_block(‘powergui/Powergui’, [model_name ‘/powergui’], ‘Position’, [20 20 80 60]);% — 交流系统 —add_block(‘powerlib/Electrical Sources/Three-Phase Source’, …[model_name ‘/AC Source (Rect)’], …‘VoltageAmplitude’, ‘500e3/sqrt(3)’, …‘Frequency’, ‘60’, …‘Position’, [50 200 100 240]);add_block(‘powerlib/Electrical Sources/Three-Phase Source’, …[model_name ‘/AC Source (Inv)’], …‘VoltageAmplitude’, ‘345e3/sqrt(3)’, …‘Frequency’, ‘60’, …‘Position’, [850 200 900 240]);% — 变压器 —add_block(‘powerlib/Elements/Three-Phase Transformer (Two Windings)’, …[model_name ‘/Transformer (Rect)’], …‘Nominal power and frequency [Pn(VA) fn(Hz)]’, ‘[5000e6 60]’, …‘Winding 1 parameters [V1(P-P Phasor) R1(pu) L1(pu)]’, ‘[500e3 0.002 0.08]’, …‘Winding 2 parameters [V2(P-P Phasor) R2(pu) L2(pu)]’, ‘[230e3 0.002 0.08]’, …‘Position’, [150 200 200 280]);add_block(‘powerlib/Elements/Three-Phase Transformer (Two Windings)’, …[model_name ‘/Transformer (Inv)’], …‘Nominal power and frequency [Pn(VA) fn(Hz)]’, ‘[5000e6 60]’, …‘Winding 1 parameters [V1(P-P Phasor) R1(pu) L1(pu)]’, ‘[230e3 0.002 0.08]’, …‘Winding 2 parameters [V2(P-P Phasor) R2(pu) L2(pu)]’, ‘[345e3 0.002 0.08]’, …‘Position’, [750 200 800 280]);% — 换流器 (使用通用桥臂简化表示) —add_block(‘powerlib/Power Electronics/Universal Bridge’, …[model_name ‘/Rectifier’], …‘Number of bridge arms’, ‘6’, …‘Snubber resistance Rs (Ohms)’, ‘100000’, …‘Forward voltage [Von(V) Rd(Ohm)]’, ‘[1.5 0.001]’, …‘Position’, [250 200 300 300]);add_block(‘powerlib/Power Electronics/Universal Bridge’, …[model_name ‘/Inverter’], …‘Number of bridge arms’, ‘6’, …‘Snubber resistance Rs (Ohms)’, ‘100000’, …‘Forward voltage [Von(V) Rd(Ohm)]’, ‘[1.5 0.001]’, …‘Position’, [650 200 700 300]);% — 直流线路 (简化为电阻电感) —add_block(‘powerlib/Elements/Series RLC Branch’, …[model_name ‘/DC Line R’], …‘Branch type’, ‘R’, …‘Resistance R (Ohms)’, ‘5’, … % 300km线路的近似电阻‘Position’, [450 240 500 260]);add_block(‘powerlib/Elements/Series RLC Branch’, …[model_name ‘/DC Line L’], …‘Branch type’, ‘L’, …‘Inductance L (H)’, ‘0.5’, … % 300km线路的近似电感‘Position’, [510 240 560 260]);% — 故障模块 —add_block(‘powerlib/Extras/Discontinuities/Breaker’, …[model_name ‘/DC Fault’], …‘Initial status’, ‘closed’, …‘Position’, [480 300 530 340]);add_block(‘powerlib/Extras/Discontinuities/Breaker’, …[model_name ‘/A-G Fault’], …‘Initial status’, ‘closed’, …‘Position’, [820 300 870 340]);% — 滤波器 (简化为一个并联RLC支路) —add_block(‘powerlib/Elements/Three-Phase Series RLC Load’, …[model_name ‘/AC Filter A (Rect)’], …‘Nominal voltage Vln (Vrms)’, ‘230e3/sqrt(3)’, …‘Nominal frequency fn (Hz)’, ‘60’, …‘Active power P (W)’, ‘0’, …‘Inductive reactive power QL (var)’, ‘-300e6’, … % 容性无功‘Position’, [200 350 250 400]);% … (其他滤波器和无功补偿模块可以类似方式添加)% — 2. 连接模块 (简化连接示例) —% 连接整流侧add_line(model_name, ‘AC Source (Rect)/1’, ‘Transformer (Rect)/1’);add_line(model_name, ‘Transformer (Rect)/2’, ‘Rectifier/1’);add_line(model_name, ‘Rectifier/2’, ‘DC Line R/1’);add_line(model_name, ‘DC Line R/1’, ‘DC Line L/1’);add_line(model_name, ‘DC Line L/1’, ‘Inverter/2’);add_line(model_name, ‘Inverter/1’, ‘Transformer (Inv)/2’);add_line(model_name, ‘Transformer (Inv)/1’, ‘AC Source (Inv)/1’);% 连接故障add_line(model_name, ‘DC Fault/1’, ‘DC Line R/1’); % 一端接在直流线路上add_line(model_name, ‘DC Fault/2’, ‘DC Line R/1’); % 另一端也接在同一点通过闭合实现短路% — 3. 保存模型 —save_system(model_name);disp([模型已创建并保存为: , model_name, ‘.slx’]);disp(‘请注意此脚本创建了基本拓扑详细的参数、测量模块和控制回路需要您手动或在后续脚本中完善。’);第三部分核心控制与换相失败抑制代码这部分是课题的灵魂。我们将创建一个 MATLAB Function 模块的核心代码它实现了带预测功能的熄弧角控制器。您可以将此代码放入 Simulink 的 MATLAB Function 模块中。代码功能实时计算熄弧角 gamma。检测交流电压跌落判断是否有换相失败风险。启动预测控制一旦检测到风险立即增大触发超前角 beta (即减小 alpha)为换相过程预留更多时间从而保证 gamma gamma_{min}。function alpha_out predictive_gamma_controller(V_ac_meas, I_dc_meas, gamma_ref, sys_params)%#codegen% 预测型熄弧角控制器% 输入:% V_ac_meas: 测量的交流线电压幅值 (V)% I_dc_meas: 测量的直流电流 (A)% gamma_ref: 参考熄弧角 (rad), 例如 18 度% sys_params: 包含系统参数的结构体 (N, Xc, etc.)% 输出:% alpha_out: 计算出的触发角指令 (rad)% --- 1. 参数提取 --- N sys_params.N; % 脉波数 Xc sys_params.Xc; % 换相电抗 (Ohm) gamma_min 10 * pi / 180; % 最小允许熄弧角 (rad) Kp 5.0; % 比例增益 Ki 100.0; % 积分增益 % --- 2. 状态变量 (需要使用持久变量 persistent 来保持) --- persistent gamma_integrator alpha_prev; if isempty(gamma_integrator) gamma_integrator 0; end if isempty(alpha_prev) alpha_prev 15 * pi / 180; % 初始触发角 end % --- 3. 计算当前实际熄弧角 gamma_actual --- % 公式: cos(gamma) cos(alpha) - (2 * Xc * Id) / (sqrt(2) * Vll) % 我们用上一拍的 alpha_prev 来估算当前的 gamma cos_alpha_prev cos(alpha_prev); term (2 * Xc * I_dc_meas) / (sqrt(2) * V_ac_meas); cos_gamma_est cos_alpha_prev - term; % 边界处理 if cos_gamma_est 1 gamma_actual 0; elseif cos_gamma_est 1 alpha_out 0; % alpha 最小为0 elseif cos_alpha_req max_delta_alpha alpha_out alpha_prev sign(delta_alpha) * max_delta_alpha; end else % --- 正常运行模式 (定 gamma 控制) --- error gamma_ref - gamma_actual; % 简单的PI控制律 d_gamma_integrator Ki * error; gamma_integrator gamma_integrator d_gamma_integrator; alpha_out alpha_prev Kp * error gamma_integrator; % 限制 alpha 的输出范围 alpha_min 5 * pi / 180; alpha_max 160 * pi / 180; if alpha_out alpha_max alpha_out alpha_max; gamma_integrator 0; % 抗积分饱和 end end % --- 5. 更新状态 --- alpha_prev alpha_out;end您需要手动添加测量模块如 Three-Phase V-I Measurement来获取 V_{ac} 和 I_{dc}。 添加脉冲发生器如 Six-Pulse Generator并将上述 predictive_gamma_controller 的计算结果 alpha_out 连接到它的触发角输入端。 配置好所有模块的参数使其与您的课题数据匹配。集成控制代码在 Simulink 库中找到 User-Defined Functions - MATLAB Function。将其拖入您的模型双击打开将 predictive_gamma_controller 的代码粘贴进去。