充电站谐波仿真：单机等效电阻到多机并联抵消规律

简介「电动汽车充电站仿真模型及其对电网谐波影响」是一篇刊于《电工技术学报》的期刊论文PDF面向电气工程、电力电子及新能源汽车充电设施方向的高校师生与科研人员聚焦大功率充电机作为非线性用电设备接入电网后引发的谐波问题。文件包仅含1个PDF文档约374KB篇幅紧凑适合快速通读与引用参考。论文从合理充电技术与合理充电规模出发分别搭建单台充电机与整站仿真模型重点分析各次谐波电流含有率、电流总谐波畸变率及功率因数随充电功率和充电机台数变化的规律并给出大功率、小功率充电下谐波变化的差异化结论可为充电站规划建设与电力系统设计提供数据参考。目前已有188人学习下载。读者可从中获得整流装置谐波来源、仿真建模思路与功率因数耦合机制等完整论述便于直接用于课题研究、论文写作或工程方案论证。1. 充电站谐波不是“单台乘以台数”做充电站电能质量评估时最省事的估算方式是拿一台充电机的电流频谱乘上台数直接当作整站注入电网的谐波。示范站实测数据经常把这种算法的结果打脸——实测谐波反而更小。原因在于每台充电机启动时刻、电池荷电状态、整流触发时刻都不一致谐波电流相位角各不相同多台并联时向量和会出现互相抵消。换句话说充电站谐波与台数不是线性叠加关系而是带相位多样性的耦合关系。这篇重庆大学的《电动汽车充电站仿真模型及其对电网谐波影响》就是围绕这个问题来的先建单台高频充电机的等效模型再扩到多机并联的充电站模型用 Matlab/Simulink 沿充电功率和充电机台数两个维度扫参看各次谐波电流含有率、电流总谐波畸变率THD和功率因数怎么变。它适合做配电网规划、充电站电气设计、电能质量评估的工程师也适合刚上手 Simulink 电力电子仿真、想找一个可复现案例练手的人。2. 充电机为什么能用一个电阻等效2.1 高频充电机的拓扑与非线性源头当前主流的大功率充电机结构是三相桥式不可控整流电路把 380V 交流整成直流经过直流滤波后再送入高频 DC-DC 功率变换电路最后经输出滤波给动力电池充电。整条链路里整流环节是谐波的根源——三相桥式二极管整流属于强非线性环节交流侧电流不是正弦而是带凹陷的准方波谐波含量可观。后面的高频 DC-DC 因为开关频率远高于工频它在工频尺度上看不出“纹波式”的非线性只是决定充电机在哪个工作点上取多少功率。理解这一点很关键谐波的“形状”由整流桥决定谐波的“幅度”由 DC-DC 决定的功率决定。仿真模型只要把整流桥画对再把后级等效成一个合适的负载就能把谐波幅频特性复现出来。2.2 6 脉波整流留下的谐波指纹三相桥式不可控整流输出为 6 脉波理想情况下只含 6k±1 次特征谐波即 5、7、11、13、17、19、23、25、29、31……次。偶次谐波和 3 次谐波因为三相对称性被消掉。对某示范站的单台充电机实测报告印证了这一点一次完整充电过程输出电流呈阶梯下降曲线直流侧充电电流依次经历 320A、240A、160A、80A 四个阶段每阶段持续约 56min各阶段 6k±1 次电流谐波含有率均超过 GB/Z 17625.6-2003 规定的接入条件输出功率越大特征谐波电流 6k1 次的幅值也越大四个阶段中 5 次、7 次是最主要的谐波分量各阶段偶次及 3 次谐波含有率均未越限。注意这里有一个容易反直觉的点幅值随功率增大而增大但含有率相对基波的比例反而随功率增大而降低。后文扫参时这个现象会反复出现。2.3 恒功率工作点带来的负阻抗特性如果把高频 DC-DC 之后的电池充电过程放在几个工频周期的时间尺度上看电池侧的充电电压 Uo 和充电电流 Io 基本不变也就是 DC-DC 环节工作在恒功率状态。恒功率 输入功率固定当输入电压升高时输入电流必然降低这意味着它的输入阻抗呈现负阻抗特性。在工频及以下频段内这个负阻抗环节可以用一个非线性电阻近似模拟——不是随便挑一个电阻而是让它满足功率与电源电压的约束关系。论文给出的等效关系是R η·Uo² / Po其中 η 是充电机效率Uo 是电池侧等效直流电压Po 是电池侧充电功率。整理一下物理含义小电阻对应大功率充电大电阻对应小功率充电。提示用电阻等效不是“偷懒”而是把恒功率环节从仿真中解耦出去。这样可以绕开高频开关器件的开关损耗和步长问题把整机仿真的时间尺度从微秒放宽到工频周期量级扫参一次几百个工况才跑得动。2.4 充电功率与等效电阻的映射论文使用的 120Ω 电阻与对应充电功率如下表节选几个关键点完整表可自行扩展等效电阻 / Ω充电功率 / kW工况含义1262.8大功率快充2131.4快充552.6中等功率1026.3常规充电2013.1小功率慢充如果换成脚本化扫参用 MATLAB 一行就能生成这张表R_eq 1:20; % 等效电阻 Ω Uo 513; % 电池侧等效直流电压 V eta 0.95; % 充电机效率 P_kW eta * Uo^2 ./ R_eq / 1e3; % 充电功率 kW disp([R_eq(:), P_kW(:)]);逻辑说明R_eq 是扫参自变量Uo、eta 是充电机本身的常量P_kW 是因变量。用这个映射把“功率扫描”转换为“电阻扫描”Simulink 里只需要改一个负载电阻值就能换工作点。注意Uo、η 的取值会直接改变功率尺度。做你自己的项目时不要照抄 513V 和 0.95用实测电池包成组电压和充电机铭牌效率替换。3. 在 Simulink 里把单台充电机跑通3.1 主电路模块与关键参数论文的仿真在 Matlab/Simulink 加 SimPowerSystems现在叫 Simscape Electrical环境下完成。单台充电机模型的关键模块与典型参数模块关键参数推荐取值Three-Phase Source线电压、频率、短路容量380V、50Hz、按母线实际值Three-Phase Series RLC线路电阻、电感模拟电源内阻抗Universal Bridge桥臂数、器件类型3 臂、Diodes直流滤波电容 C、电感 L按滤波截止频率选等效负载电阻 R_eq120Ω 扫参powerguiSimulation modeContinuousSolver类型、最大步长ode23tb、1e-5 s两个参数最容易踩坑电源短路容量决定总线电压刚性短路容量太小会让谐波电压畸变反过来影响充电机电流扫出来的谐波含有率会偏高直流侧滤波电容决定整流桥的导通角电容过大导通角过窄电流波形更尖THD 会虚高。做对比实验时这两个量必须固定否则台数和功率的灵敏度分析没有意义。3.2 脚本化搭模型在 GUI 里点几千次扫参是不现实的。常见做法是用脚本新建模型、加模块、连信号、批量仿真。核心片段model EVCharger; if bdIsLoaded(model), close_system(model, 0); end new_system(model); open_system(model); add_block(powerlib/powergui, [model /powergui], ... SimulationMode, Continuous); add_block(powerlib/Electrical Sources/Three-Phase Source, ... [model /Source], ... PhaseToPhaseVoltage, 380, Frequency, 50); add_block(powerlib/Power Electronics/Universal Bridge, ... [model /Bridge], Arms, 3, PowerElectronicDevice, Diodes); add_block(powerlib/Elements/Series RLC Branch, ... [model /Rload], BranchType, R, Resistance, 1); add_block(simulink/Sinks/To Workspace, [model /ia], ... VariableName, ia, SaveFormat, Structure With Time); % 连线并配置求解器 add_line(model, Source/1, Bridge/1); add_line(model, Bridge/2, Rload/1); set_param(model, Solver, ode23tb, MaxStep, 1e-5, ... StopTime, 0.2);逻辑说明add_block 按 Simulink 库路径添加模块set_param 在添加时直接写入参数add_line 完成端口连线最后 set_param 配置求解器和仿真时长。这个脚本只画出“电源 整流桥 等效电阻”的最小可跑版本还需要补全三相电源的 RLC 内阻抗支路和直流滤波支路。参数说明StopTime 取 0.2s 是给整流启动暂态留出 10 个工频周期稳态取后面几个周期做 FFTMaxStep 设 1e-5 是为了让整流换相瞬间被采到太大会导致高频毛刺被平掉谐波分布失真。3.3 FFT 提取各次谐波含有率与 THD仿真跑完从 ia 里截稳态的一整段做 FFT这是把时域波形转成谐波指标的必经一步Fs 1e4; % 采样频率需与 MaxStep 匹配 N 2000; % 窗口点数2000/1e4 0.2s t (0:N-1)/Fs; ia out.ia.Data(end-N1:end); Y fft(ia)/N; mag abs(Y(1:N/21))*2; % 单边幅值谱 f (0:N/2)*Fs/N; I1 mag(find(f 50, 1)); % 基波幅值 h [5 7 11 13 17 19 23 25 29 31]; hr zeros(size(h)); for k 1:numel(h) hr(k) mag(find(f h(k)*50, 1)) / I1 * 100; % 含有率 % end non_fund mag(f 100); % 2 次及以上 THD_i sqrt(sum(non_fund.^2)) / I1 * 100;逻辑说明FFT 归一化后 abs 求幅值乘 2 还原为单边谱find(f50,1) 取最接近工频的谱线作为基波hr 依次算出 6k±1 特征次含有率THD_i 用二次方和的平方根除以基波。参数说明N 必须覆盖整数个工频周期否则会引入频谱泄漏如果无法保证整周期就要加窗函数或者用 Simulink 自带的 powergui FFT Analysis 工具代替。3.4 单机扫参结果怎么读把 R_eq 从 1Ω 扫到 20Ω会看到三组趋势5 次、7 次谐波含有率随充电功率下降电阻增大而明显上升11 次、13 次及以上高次谐波含有率随功率变化较平缓且谐波次数越高含有率越低电流总谐波畸变率随充电功率增大而减小功率因数随充电功率增大而增大。第二条容易理解第三条里功率因数的走向其实和整流桥导通角有关大功率时直流母线电压升高、二极管导通角变宽电流波形向正弦收紧位移因数和畸变因数都改善总功率因数自然上去了。4. 多机并联充电站模型的搭建与抵消规律4.1 把单机封装成可复用的子系统多机模型不是把单机模型复制 15 遍就完事需要先把单机里的电源内阻抗拿出来放到母线上单机本体只保留“整流桥 滤波 等效电阻”三部分封装成 Subsystem。封装时给一个参数外露口 R_eq每个 Subsystem 实例可以设为不同的电阻值模拟不同 SOC 阶段的充电机。常见的母线模型是一个共享的三相电源 母线短路阻抗母线下挂 N 个并联充电机子系统。每台充电机的初始投入时刻错开几个毫秒让各机的相位在统计上独立这是让“相位多样性”出现的关键做法。4.2 相位多样性与谐波抵消的仿真验证把多机模型的 A 相总电流和 N 个“单机电流乘台数”的等效电流做对比可以看到明显的差别低频段 5、7 次谐波差异最大高次谐波因为幅值小差异相对不明显。论文的仿真结论与这个现象一致——不同充电机产生的谐波存在补偿甚至互相抵消所以多台并联时各次谐波含有率和 THD 不是单机值的线性叠加。4.3 大功率与小功率两种台数增长曲线论文给出了两个典型工况的对比大功率充电 R1ΩP262.8kW随充电机台数增加各次谐波电流含有率呈减小趋势相同台数下5、7 至 31 次谐波含有率依次递减。小功率充电 R20ΩP13.1kW随台数增加各次谐波含有率变化较平缓几乎不随台数显著下降。两个工况的差异本质上是相位多样性能不能发挥作用。大功率时单机谐波幅值大各机之间只要相位差上几十度向量和的合成就能明显低于代数和小功率时单机谐波幅值本来就小向量和的波动被“淹没”在基波附近的噪声里曲线上看不出明显趋势。4.4 THD 随台数衰减与最优台数区间把充电机台数 N 从 1 扫到 15、等效电阻从 1Ω 扫到 20Ω绘制 THD 的三维曲面能看到两个规律THD 随台数增加整体呈减小趋势且充电功率越大减小幅度越显著台数从 10 台开始各功率下的 THD 减小趋势都逐渐变缓。也就是说多装充电机带来的谐波改善存在明显的边界收益递减。充电站里不同充电机设定的充电功率不同最优台数也会随之变化。对规划人员来说仿真给出的不是唯一的“最优台数”而是一张 台数-功率-THD 曲面最优点需要结合变压器容量和负荷预测来切。5. 功率因数与 THD 耦合曲线的工程读法论文最值得反复看的一句结论是功率因数的变化由充电功率与充电机数目的耦合机制决定。这意味着不能只看单机功率因数数据就乘个台数去估算整站。实操上怎么做第一步用 4.4 的 THD 曲面做 Pareto 筛选。把 (N, R) 组合按 THD 和台数两目标排序从 Pareto 前沿里再选。下面这段脚本读入仿真生成的三维表找出每个功率下 THD 最小的台数N 1:15; R 1:20; THD load(thd_surface.mat); % 行为 R列为 N THD THD.THD; [minTHD, k] min(THD, [], 2); % 每个功率下的最优台数 figure; plot(R, N(k), o-, LineWidth, 1.2); xlabel(等效电阻 / \Omega); ylabel(THD 最优台数); % 功率因数切片 PF load(pf_surface.mat); PF PF.PF; plot(R, PF(:, 1), k--, R, PF(:, 15), r-); legend(1 台, 15 台);逻辑说明min(THD,[],2)沿列方向取每行最小值对应每个功率水平下 THD 最小的台数第二个 plot 找两个极端台数下的功率因数曲线做对比观察耦合走向。注意如果仿真数据只到 15 台曲线在 15 台处可能还没收敛规划时不要直接把最优台数取在扫描范围边界上至少要再多扫 5 台确认趋势。第二步做短路容量敏感性。把母线的短路容量从基准值改到 0.8 倍和 1.2 倍重跑仿真观察同一 THD 水平下允许的最大台数。短路容量越小同样谐波电流下引起的母线电压畸变越大充电机之间相位耦合更强THD 曲面的形状会明显变形。这一步常被跳过但它决定了前面算出来的“最优台数”能不能直接搬到实际配电站使用。第三步把结果落到滤波器选型。5、7 次是主导谐波分量无源滤波优先按这两个频次设计11 次以上含有率低但次数高用电抗率 6% 或 7% 的串联电抗器搭配并联电容器抑制即可。功率因数补偿容量不要按单机额定功率因数反算要用 4.4 的曲面在目标台数下取点否则补偿过头会把容性无功倒送到上级电网。第四步验证方法。仿真结果最好和示范站实测对比一次对比口径统一到同一功率水平、同一采样窗口、同一 FFT 参数。如果 5、7 次含有率的偏差超过 2 个百分点先回查直流滤波电容的取值再检查等效电阻是不是没把充电机效率 η 算进去——这两处是差得最多的地方。本文还有配套的精品资源点击获取