简介：「电动汽车充电站仿真模型及其对电网谐波影响」是一篇刊于《电工技术学报》的期刊论文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 四个阶段，每阶段持续约 5～6min；
- 各阶段 6k±1 次电流谐波含有率均超过 GB/Z 17625.6-2003 规定的接入条件，输出功率越大，特征谐波电流 6k+1 次的幅值也越大；
- 四个阶段中 5 次、7 次是最主要的谐波分量，各阶段偶次及 3 次谐波含有率均未越限。
注意，这里有一个容易反直觉的点：幅值随功率增大而增大，但含有率（相对基波的比例）反而随功率增大而降低。后文扫参时这个现象会反复出现。
2.3 恒功率工作点带来的负阻抗特性
如果把高频 DC-DC 之后的电池充电过程放在几个工频周期的时间尺度上看，电池侧的充电电压 Uo 和充电电流 Io 基本不变，也就是 DC-DC 环节工作在恒功率状态。恒功率 = 输入功率固定；当输入电压升高时，输入电流必然降低，这意味着它的输入阻抗呈现负阻抗特性。
在工频及以下频段内，这个负阻抗环节可以用一个非线性电阻近似模拟——不是随便挑一个电阻，而是让它满足功率与电源电压的约束关系。论文给出的等效关系是：
R = η·Uo² / Po
其中 η 是充电机效率，Uo 是电池侧等效直流电压，Po 是电池侧充电功率。整理一下物理含义：小电阻对应大功率充电，大电阻对应小功率充电。
提示：用电阻等效不是“偷懒”，而是把恒功率环节从仿真中解耦出去。这样可以绕开高频开关器件的开关损耗和步长问题，把整机仿真的时间尺度从微秒放宽到工频周期量级，扫参一次几百个工况才跑得动。
2.4 充电功率与等效电阻的映射
论文使用的 1～20Ω 电阻与对应充电功率如下表（节选几个关键点，完整表可自行扩展）：
| 等效电阻 / Ω | 充电功率 / kW | 工况含义 |
|---|---|---|
| 1 | 262.8 | 大功率快充 |
| 2 | 131.4 | 快充 |
| 5 | 52.6 | 中等功率 |
| 10 | 26.3 | 常规充电 |
| 20 | 13.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_eq | 1～20Ω 扫参 |
| powergui | Simulation mode | Continuous |
| Solver | 类型、最大步长 | 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 个工频周期，稳态取后面几个周期做 FFT；MaxStep 设 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-N+1:end); Y = fft(ia)/N; mag = abs(Y(1:N/2+1))*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(f>=50,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 大功率与小功率两种台数增长曲线
论文给出了两个典型工况的对比：
- 大功率充电 R=1Ω（P=262.8kW）：随充电机台数增加，各次谐波电流含有率呈减小趋势；相同台数下，5、7 至 31 次谐波含有率依次递减。
- 小功率充电 R=20Ω（P=13.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 个百分点，先回查直流滤波电容的取值，再检查等效电阻是不是没把充电机效率 η 算进去——这两处是差得最多的地方。
本文还有配套的精品资源，点击获取
