简介围绕电动汽车充电站建模与电网谐波影响展开的PDF文献面向新能源汽车、电力电子与电网规划方向的师生及工程技术人员适合作为参考文献与专业指导材料。内容从大功率充电机整流装置的非线性特征切入分别建立单台充电机与充电站仿真模型重点分析各次谐波电流含有率、电流总谐波畸变率和功率因数随充电功率及充电机台数变化的规律并讨论适应大量电动汽车充电需求的合理充电技术与充电规模问题。压缩包内仅1个PDF文件约374KB篇幅精炼便于在电脑或移动端快速查阅与引用。已有188人学习下载可作为电动汽车充电站设计、电能质量评估及电力系统规划时的推导依据与数据参照也能帮助读者建立谐波分析的仿真思路。1. 从一台 60 kW 直流快充桩说起充电站谐波为什么不能只看单机某园区装了 10 台 60 kW 直流快充桩单桩出厂测试并网电流总谐波畸变率不到 5%但十台同时满功率输出时并网点电流畸变明显上升无功补偿柜温度异常。单机合格不等于站级合格多台充电机的前级整流挂在同一段母线上谐波电流的相位关系、电网侧阻抗和补偿电容支路的谐振点一起决定了最终畸变水平。下面几章围绕电动汽车充电站仿真模型的搭建与电网谐波影响的量化展开从单桩三级拓扑拆到站级并联再落到 FFT 提取谐波、THD 计算和阻抗谐振排查。内容面向做充电站配电设计、电能质量评估和充电桩 EPC 的工程人员也适合需要把仿真结论写进并网评估报告的人。2. 充电站仿真模型的三级拓扑与建模边界2.1 从充电模块到站级三级拓扑怎么拆电动汽车充电站仿真模型的第一步不是打开 Simulink而是决定哪一级进模型、哪一级用等效替代。直流快充桩的功率链常见为三级前级三相整流把 380 V 交流变成约 700 V 直流母线中间级是隔离 DC/DC多见于全桥 LLC 谐振变换器或移相全桥把母线电压变换为 2001000 V 可调输出末级是输出滤波、采样与充电控制逻辑。从电网谐波的角度看只有前级整流级是与电网直接交换电流的非线性环节250 次谐波电流几乎全部由它产生。中间级开关频率通常在几十千赫到一百千赫以上其谐波落在高频段超出工频谐波评估的关注范围但它决定直流母线电压纹波和整流级实际工作点所以建模时应该用等效功率负载替代而不是整段删掉。站级模型则是把 N 个单桩模型挂到同一台变压器低压侧母线上母线上再并联无功补偿电容器和站内其他负荷。站级模型的重点从器件开关过程转向谐波电流叠加关系与电网侧阻抗特性建模粒度可以比单桩粗但电网侧参数必须准。2.2 不控整流、PWM 整流与 Vienna 整流的选型对比不同整流拓扑对电网谐波的影响差别很大仿真前先确认实际桩用的是哪一种否则后面的 THD 结果没有参考价值。工程上常见的四类方案对比如下。整流拓扑并网电流 THD 量级功率因数能量双向仿真复杂度三相不控整流 电容滤波高常达数十个百分点低否低三相不控整流 LC 滤波/PFC中等较高否中三相 PWM 整流VSR低可调接近 1是高Vienna 整流低接近 1否高选型判断的逻辑是站里如果是 PWM 整流桩谐波治理压力小仿真重点放在高频开关纹波和共模路径如果是老旧不控整流桩谐波电流对交流侧电感和电网阻抗非常敏感必须把这两个参数建准。Vienna 整流在不控器件数量上有优势但三相桥臂的开关函数建模需要处理中点电位平衡仿真步长要求更高。2.3 全桥 LLC 谐振级要不要在谐波仿真里保留结论是做 250 次工频谐波评估时不保留开关级细节但要保留它的功率传输特性。常见做法是用基波近似FHA单独算一遍 LLC 增益曲线确认在所设母线电压和输出功率下它能正常工作然后在站级模型里把 LLC 换成受控功率负载或用效率曲线折算的等效电阻。% 全桥LLC谐振腔基波近似(FHA)增益用于判断输出调节能力 Vin 700; % 直流母线电压 V n 1.2; % 变压器原副边变比 Lr 30e-6; Cr 33e-9; Lm 150e-6; % 谐振电感、谐振电容、励磁电感 Rac 8*2.5^2/pi^2; % 副边负载折算到原边的等效交流电阻欧姆 fr 1/(2*pi*sqrt(Lr*Cr)); % 谐振频率 fsw linspace(0.5*fr, 2*fr, 200); M zeros(size(fsw)); for k 1:numel(fsw) w 2*pi*fsw(k); w0 2*pi*fr; Q sqrt(Lr/Cr)/Rac; % 品质因数负载越重Q越大 H 1/sqrt((1 Lr/Lm*(1 - w0^2/w^2))^2 Q^2*(w/w0 - w0/w)^2); M(k) n*H; % 归一化电压增益 end plot(fsw/fr, M); grid on; xlabel(归一化频率 f_{sw}/f_r); ylabel(增益 M);这段代码里Lr、Cr决定谐振频率frLm与Lr的比值影响增益曲线的陡度Q由负载决定。当fsw接近fr时增益趋近 1低于fr进入升压区。把这条曲线和实际充电电压范围对照就能确认 DC/DC 级不会成为限制条件之后在谐波模型里就可以放心地把它简化掉。2.4 仿真步长、采样率与仿真时长怎么定步长和时长直接决定 FFT 结果的可用性。步长太大会把高次谐波的谱线抹平时长太短则频率分辨率不够相邻谐波会糊在一起。关注最高谐波次数最高频率建议最大步长建议仿真时长25 次1250 Hz2e-5 s≥1 s约 50 个工频周期50 次2500 Hz1e-5 s≥2 s含开关级细节100 kHz 以上1e-7 s0.20.5 s提示先按上表定步长再用Fs/N检查频率分辨率确保能分开 3 次和 5 次谐波否则后面算出的 THD 会偏大。3. 在 Simulink 里把站级充电站仿真模型跑起来3.1 单桩最小可运行模型与关键参数最小可用模型包含五个部分理想三相电源与等效电网阻抗、三相整流桥、直流母线电容、等效恒功率负载、并网电流测量。Simscape Electrical 里有现成的三相桥和 PLL 模块直接用即可不必从开关器件搭起。参数取值示例作用交流线电压380 V决定整流桥直流侧电压上限电网等效电感0.5 mH影响换相角和电流上升率电网等值电阻0.02 Ω影响谐振阻尼直流母线电容2200 μF抑制母线电压纹波单桩额定功率60 kW决定等效负载电流开关频率PWM 整流时10 kHz与载波相位分散相关搭建顺序建议是先建理想电源加不控整流确认直流侧达到约 1.35 倍线电压再把电网阻抗串进去观察换相缺口是否出现。如果先把阻抗加满很容易因为换相过程导致仿真报错而不容易定位。3.2 站级并联多桩叠加与载波相位分散站级模型有三种搭法。第一种是复制 N 个单桩模型并到同一母线最直观但 N 大于 5 之后仿真速度下降明显第二种是把 N 个桩用同一套拓扑表示负载功率乘 N适合所有桩型号一致的场景第三种是把每台桩折算成诺顿等效电流源只做谐波域叠加速度最快。多桩并联时谐波电流的叠加结果强烈依赖各桩载波相位。PWM 整流桩如果载波完全同步特征次谐波会直接算术叠加把各桩载波相位按2πk/N分散开低次谐波会出现明显抵消。不控整流桩没有载波概念但各桩触发角、变压器接线组别和馈线长度差异同样会引入相位差。3.3 用 FFT 提取并网电流谐波并计算 THD仿真是手段出指标才是有用结果。下面这段脚本读取并网电流波形做加窗 FFT再算各次谐波和 THD。% 并网电流谐波分析与THD计算 f0 50; % 基波频率 Hz Fs 1e5; % 采样率 Hz建议不低于最高关注谐波的10倍 N 2^16; % 采样点数取2的幂便于FFT t (0:N-1)/Fs; ia current_sim; % 从Simulink输出或CSV读入的A相电流需整周期截断 w hann(N); % 汉宁窗抑制非整周期导致的频谱泄漏 X fft(ia(:).*w); P abs(X)/sum(w); % 幅值归一化恢复真实幅值 f (0:N-1)*Fs/N; fund interp1(f, P, f0, linear); % 基波幅值 hmax 50; H zeros(hmax,1); for n 2:hmax fh n*f0; idx f fh-2 f fh2; % 谐波谱线附近取窗内峰值 H(n) max(P(idx)); end THD sqrt(sum(H(2:end).^2))/fund*100; fprintf(基波 %.3f ATHD %.2f%%\n, fund, THD);Fs/N决定频率分辨率本例为 1.5 Hz足够分开相邻谐波。hann窗把频谱泄漏压下去代价是主瓣变宽所以取峰值时用 ±2 Hz 的窗口。interp1是为了让基波幅值不受谱线离散化影响直接取最近谱线在非同步采样时会有可观误差。3.4 批量扫描桩数与负载率生成谐波结果表实际评估要看“不同桩数 不同负载率”的组合手工改参数不现实用脚本批量跑。nPiles [1 2 4 6 8 10]; % 充电桩台数 loadRate [0.3 0.5 0.8 1.0]; % 负载率 result zeros(numel(nPiles), numel(loadRate)); for i 1:numel(nPiles) for j 1:numel(loadRate) set_param(cs_model/PileCount,Value,num2str(nPiles(i))); set_param(cs_model/LoadRate,Value,num2str(loadRate(j))); out sim(cs_model,SaveOutput,on); result(i,j) thd_calc(out.iabc); % 调用上一节的THD函数 end end disp(array2table(result,RowNames,string(nPiles), ... VariableNames,Lstring(loadRate)));set_param直接改模型工作区里对外的可见变量比重新打开模型快得多。内层每跑完一次就调用一次 THD 计算避免把所有波形都存内存。结果表按桩数为行、负载率为列能直观看出畸变随规模的上升趋势。4. 充电站对电网谐波影响的量化与谐振排查4.1 谐波指标的定义与仿真口径指标定义仿真里怎么算谐波含有率第 h 次谐波有效值与基波有效值之比FFT 后各次谱线幅值除以基波幅值总谐波畸变率 THD全部谐波含量有效值与基波有效值之比各次谐波平方和开方后除以基波总需求畸变 TDD谐波电流有效值与最大需求电流之比站级评估按容量考核时更合适谐波电流允许值按公共连接点电压等级和容量查表与 GB/T 14549 的表值逐次对照注意仿真给出的 THD 是站点注入值不等于公共连接点的最终畸变。背景谐波和系统阻抗会共同作用报告里要把这两项分开说明。4.2 多桩同时充电的谐波叠加规律谐波电流叠加不是算术相加。同一型号的不控整流桩5 次谐波电流在相位接近时接近算术和载波相位分散或馈线阻抗不一致时会出现部分抵消。下表是示例量级实际值必须由仿真或实测确定。桩数5 次谐波算术和实际叠加量级原因22.01.61.9相位接近接近算术和55.02.53.5载波相位分散产生抵消1010.035投切随机性进一步增强抵消4.3 电网阻抗扫频与并联电容谐振点排查站级谐波超标最常见的机理是谐振放大补偿电容器与系统等效电感在某个频率形成并联谐振若该频率正好落在充电桩特征谐波附近谐波电流被放大数倍。% 电网侧阻抗-频率扫描定位并联谐振点 f logspace(1, log10(5000), 500); % 10 Hz ~ 5 kHz w 2*pi*f; Ls 0.5e-3; % 系统等效电感 H R 0.02; % 等值电阻 C 20e-6; % 并联电容器组折算电容 Zg R 1j*w*Ls; Zc 1./(1j*w*C); Z (Zg.*Zc)./(ZgZc); % 并联阻抗 [~, idx] max(abs(Z)); fprintf(并联谐振频率约 %.1f Hz对应 %d 次谐波\n, f(idx), round(f(idx)/50));Ls和C决定谐振频率高低R决定谐振峰的尖锐程度。实际操作时把电网的短路容量换算成等效电感代入电容器组的实际容量就能得到一条阻抗频率特性曲线再看曲线峰值落在哪一次谐波上。若峰值与 5 次或 7 次重合优先考虑调整电容器分组投切策略或加装滤波支路。4.4 仿真结果与实测偏差的来源偏差主要来自四个方面一是背景谐波配电网本身已有畸变仿真是从零开始叠加二是电网阻抗时变运行方式一变谐振点就漂移三是多桩投切随机实际没有整齐的满载时刻四是变压器接线组别Dyn 接法能阻断 3 次及倍数次谐波向高压侧传播仿真里若用理想变压器就会丢掉这个效果。把接线组别建对、把背景谐波作为戴维南等效电源加入仿真和实测的偏差通常能压到可接受范围。5. 进阶把站级谐波仿真做薄的三条路子5.1 用诺顿等效把站级模型压成谐波电流源当充电站接入的是中压配电网、需要和几十个节点一起算潮流时逐个建单桩模型不现实。做法是先用单桩仿真得到各次谐波电流的幅值和相位把它写成一个频率相关的诺顿电流源再并联一个等效导纳代表充电站自身阻抗。% 由单桩仿真结果构造谐波域诺顿等效 h (1:50); % 谐波次数 Ih amp_h .* exp(1j*angle_h); % 各次谐波电流相量来自FFT Yh 1 ./ (Req 1j*2*pi*50*h*Leq); % 充电站等效导纳 I_inj Ih - Yh .* Vbus_ph; % 注入配网的等效谐波电流源Ih和angle_h来自第 3 章的 FFT 结果Req、Leq是充电站从并网点看进去的等效阻抗。这样配网级仿真里只需要一个谐波源模型节点数不再受桩数限制。5.2 用随机投切序列做概率性谐波评估单次仿真是确定性结果工程上更关心 95% 概率水平。把充电桩的启动时刻、功率等级按实际使用规律随机抽样跑几百次后统计各次谐波的累积分布。rng(7); % 固定随机种子保证可复现 nRun 300; thdSet zeros(nRun,1); for k 1:nRun tStart sort(rand(nPiles,1)*3600); % 1小时内随机投切 pSet randi([20 60], nPiles, 1); % 每桩功率 kW thdSet(k) run_station(tStart, pSet); % 单次站级仿真返回THD end fprintf(THD 95%% 分位值 %.2f%%最大值 %.2f%%\n, ... prctile(thdSet,95), max(thdSet));关键是把rng固定住否则每次报告的数值都不同评审时无法复现。nRun取 300 次左右通常够用再多边际收益很小。5.3 与充电站功率指令时间尺度耦合站级谐波不是恒定的。充电站和车辆之间的通信与控制标准如 ISO 15118 系列决定了功率指令下发的时间尺度功率指令一台阶式变化谐波电流的相位和幅值随之跳变。做长时间尺度评估时可以把功率指令序列作为站级模型的输入用 1 秒级步长跑时序潮流再把每个时段的稳态谐波电流按诺顿等效注入形成“时序 谐波”联合结果。落地的判断标准很直接如果联合仿真给出的 95% 分位 THD 与实测相差在一个百分点以内这套薄模型就足够支撑并网评估不需要回到全开关级模型重跑。本文还有配套的精品资源点击获取