简介这份文档面向电气工程及其自动化专业学生与铁路供电方向技术人员聚焦AT牵引供电系统的计算与仿真帮助读者掌握从参数选取到牵引网特性分析的完整流程。资源为单个doc文件压缩包约589KB内容以理论推导、计算步骤与MATLAB程序实现为主适合课程设计、毕业设计及工程入门参考。文档目录涵盖AT牵引供电系统简介、计算参数、MATLAB仿真软件说明、计算步骤图与关键函数程序并进一步展开无机车运行、机车正常运行及T—R、F—R、T—F各类短路故障下的电压电流分布计算与特性分析还整理了计算过程中常见问题与解决方法。已有84人学习读者可借此理清牵引供电计算的建模思路、函数编写方法与故障分析框架为后续仿真优化和工程应用打下基础。1. 从一份“计算模板.doc”说起AT牵引供电系统到底要算什么很多电气工程及其自动化专业的同学第一次接触 ATAuto-Transformer自耦变压器牵引供电系统是在毕业设计或课程设计里拿到一份名为“AT牵引供电系统计算模板.doc”的文件。打开一看里面是一堆公式、表格和待填参数但很少有人告诉你这些参数从哪来、算出来对不对、为什么牵引网要用 AT 而不是直接 BT 供电。AT 供电的核心价值在于它把牵引网电压从 27.5 kV 提升到 2×27.5 kV 的馈电结构通过自耦变压器把供电电压抬高、把电流压低从而大幅延长供电臂长度、降低线路损耗。对电气化铁路、城市轨道交通的供电专业来说这套计算直接决定了变电所布点、导线截面选型和短路电流水平。这份模板要解决的就是把“牵引网阻抗—电压损失—短路电流—回流分配”这条链路算清楚让设计有据可依。适合正在做供电系统设计、需要把手工计算升级成可复现脚本的从业者和学生。2. AT牵引供电系统的计算模型与MATLAB建模准备2.1 AT牵引网的结构与等效电路怎么理解AT 牵引供电系统由牵引变电所、AT 所、接触网T、正馈线F、钢轨R和地G构成。接触网与正馈线对地电压各为 27.5 kVAT 所把两者串起来使接触网对正馈线的电压达到 55 kV。列车取流时电流从接触网流入机车再经钢轨和地返回AT 绕组按 1:1 变比把回流“吸”回正馈线这就是所谓“吸流效应”。做计算前必须先把多导体传输线简化成等值电路。常见做法是把接触网—钢轨、正馈线—钢轨、钢轨—地之间的互阻抗和自阻抗写成阻抗矩阵再用 AT 的边界条件绕组电流关系消元。很多模板直接给出“单位长度阻抗”经验值但真正要复现得从导线几何参数算起。参数符号典型取值说明接触网单位阻抗z_T0.15j0.55 Ω/km含接触线承力索正馈线单位阻抗z_F0.12j0.42 Ω/kmLGJ 型导线钢轨单位阻抗z_R0.10j0.35 Ω/km含回流线AT 漏抗z_AT0.05j0.30 Ω单台折算AT 间距D10~15 km决定供电臂长度注意不同线路的导线型号、轨型、大地电阻率差异很大上表只是量级参考实际必须按设计资料代入。2.2 用MATLAB搭一个最小可跑的AT计算脚本下面这段脚本把上面的等值电路落成矩阵运算算一段供电臂的电压损失。它不依赖任何工具箱纯基础语法方便对照模板里的公式。% AT牵引供电系统电压损失最小计算示例 clear; clc; % 单位长度阻抗 (Ω/km) zT 0.15 1j*0.55; % 接触网 zF 0.12 1j*0.42; % 正馈线 zR 0.10 1j*0.35; % 钢轨 % AT间距与供电臂长度 D 12; % km L 36; % 供电臂总长 km % 列车电流 (A)功率因数0.95 I 600 * exp(-1j*acos(0.95)); % 简化按AT分段叠加每段回流系数取0.6 k 0.6; % 逐段累加电压损失 dU 0; for x 0:D:L seg min(D, L - x); % 接触网压降减去正馈线回流补偿 dU dU (zT - k*zF) * I * seg; end fprintf(供电臂末端电压损失 %.2f V\n, abs(dU)); fprintf(折算到27.5kV侧百分比 %.2f%%\n, abs(dU)/27500*100);逻辑说明循环按 AT 间距把供电臂切成若干段每段用(zT - k*zF)表示接触网压降被正馈线回流部分抵消后的等效阻抗k是回流系数反映吸流效果。参数说明D改小会让分段更细、结果更接近分布参数k越接近 1说明 AT 吸流越充分电压损失越小I用复数是为了同时保留幅值和功率因数信息。跑完看百分比如果超过 10%通常就要缩短供电臂或加 AT 所。2.3 模板里的公式和脚本结果怎么对齐拿到“计算模板.doc”后别急着填数。先把模板里每个公式的假设抄出来是集中参数还是分布参数回流系数是常数还是随位置变短路计算用的是标幺值还是有名值然后把你脚本里对应的量改成同样假设跑一遍对比。常见差异来自三点一是模板把 AT 漏抗忽略脚本里却算了二是模板用平均电流脚本用末端最大电流三是模板的钢轨阻抗含大地回流修正脚本用了简化值。对齐这三点误差一般能压到 5% 以内。3. 牵引网短路电流与回流分配的MATLAB实现3.1 短路电流计算的边界条件怎么设AT 牵引网短路分接触网对钢轨短路、接触网对正馈线短路等类型。短路点不同AT 绕组提供的回路也不同。做短路计算时核心是把所有 AT 所当作理想变压器处理写出节点导纳矩阵再在短路点注入单位电流求电压分布。常见做法是把每个 AT 所看作一个理想变压器节点接触网节点和正馈线节点通过 AT 绕组耦合钢轨和地合并成一个回流节点。短路时在接触网节点和回流节点之间加一个很小的过渡电阻解网络方程得到短路电流。% AT牵引网短路电流计算节点法简化版 nAT 3; % AT所数量 Zat 0.05 1j*0.30; % AT漏抗 Zs 0.10 1j*0.35; % 钢轨回流阻抗 Zt 0.15 1j*0.55; % 接触网阻抗 D 12; % AT间距 km % 构造节点导纳矩阵接触网节点 回流节点 N nAT 1; Y zeros(2*N, 2*N); for i 1:N Y(i,i) Y(i,i) 1/Zt; Y(iN,iN) Y(iN,iN) 1/Zs; if i N Y(i,i1) Y(i,i1) - 1/Zt; Y(i1,i) Y(i1,i) - 1/Zt; end % AT耦合 Y(i,iN) Y(i,iN) - 1/Zat; Y(iN,i) Y(iN,i) - 1/Zat; end % 短路点设在第2个AT处接触网对回流 k 2; Y(k,k) Y(k,k) 1e6; % 近似理想短路 U Y \ zeros(2*N,1); fprintf(短路点电压 %.2f V\n, abs(U(k)));逻辑说明导纳矩阵按接触网节点和回流节点分块AT 耦合项放在对角块之间。1e6是人为加大短路点导纳模拟金属性短路。参数说明nAT决定矩阵规模实际线路可能有 5~8 个 AT 所Zat越小AT 提供的短路回路越强短路电流越大。算完电压后短路电流等于短路点两侧支路电流之和需要再回代一次。3.2 回流分配系数对结果的影响回流分配是 AT 系统最容易被忽略的部分。列车电流经钢轨返回时一部分被 AT 吸到正馈线一部分经大地流回变电所。回流系数k不是常数它随列车位置、AT 间距、大地电阻率变化。工程上常用 0.5~0.7 的经验值但做精确计算时应该把它作为位置函数。列车位置回流系数 k电压损失影响靠近变电所0.75损失小供电臂中部0.60损失中等靠近末端0.45损失最大提示如果模板里k取固定值而你用变值算末端电压损失可能差 15% 以上这直接关系到要不要加 AT 所。3.3 用循环扫描找出最不利短路点短路电流最大的点不一定在供电臂末端也不一定在变电所出口。把短路点从变电所逐段移到末端每次重算导纳矩阵记录短路电流就能画出短路电流沿线路的分布。这个扫描用 MATLAB 写起来很短Isc zeros(1, N); for k 1:N Yk Y; Yk(k,k) Yk(k,k) 1e6; U Yk \ zeros(2*N,1); Isc(k) abs(U(k)) / abs(Zt); % 近似支路电流 end plot(1:N, Isc, -o); xlabel(AT节点编号); ylabel(短路电流 (A)); title(AT牵引网短路电流分布);逻辑说明每次把短路点移到下一个节点重解一次网络。参数说明Isc存每个节点的短路电流plot用来直观找峰值。实际工程里峰值点就是选断路器开断容量的依据。4. 把计算模板变成可复用的MATLAB函数与参数化脚本4.1 函数化改造输入参数、输出结果模板最大的问题是“一次性”改一个参数要重填整张表。把它改造成函数输入导线参数、AT 间距、列车电流输出电压损失、短路电流和回流系数就能批量跑不同方案。function [dU, Isc_max, k_avg] at_traction(zT, zF, zR, D, L, I, nAT) % at_traction 计算AT牵引供电系统电压损失与最大短路电流 % 输入单位阻抗(Ω/km)、AT间距(km)、供电臂长(km)、列车电流(A)、AT所数 % 输出末端电压损失(V)、最大短路电流(A)、平均回流系数 k 0.6; % 回流系数初值 dU 0; for x 0:D:L seg min(D, L - x); dU dU (zT - k*zF) * I * seg; end % 短路扫描 N nAT 1; Y zeros(2*N); for i 1:N Y(i,i) Y(i,i) 1/zT; Y(iN,iN) Y(iN,iN) 1/zR; if i N Y(i,i1) Y(i,i1) - 1/zT; Y(i1,i) Y(i1,i) - 1/zT; end Y(i,iN) Y(i,iN) - 1/(0.051j*0.30); Y(iN,i) Y(iN,i) - 1/(0.051j*0.30); end Isc zeros(1,N); for kk 1:N Yk Y; Yk(kk,kk) Yk(kk,kk) 1e6; U Yk \ zeros(2*N,1); Isc(kk) abs(U(kk)) / abs(zT); end Isc_max max(Isc); k_avg k; end逻辑说明函数把第 2、3 章的代码合并输入输出清晰。参数说明zT/zF/zR用复数传入保留阻抗角nAT控制矩阵规模k目前是常数可以改成随位置变化的向量。调用时一行就能跑一个方案[dU, Isc, k] at_traction(0.151j*0.55, 0.121j*0.42, ... 0.101j*0.35, 12, 36, 600*exp(-1j*acos(0.95)), 3); fprintf(dU%.1f V, Isc%.1f A\n, abs(dU), Isc);4.2 参数敏感性分析哪些参数最值得调不是所有参数都同等重要。用循环扫一遍看电压损失对哪个参数最敏感能帮你把设计精力放在刀刃上。参数变化范围电压损失变化敏感度AT间距 D8~15 km-18%~12%高回流系数 k0.45~0.7522%~-15%高接触网阻抗 zT±10%±8%中列车电流 I±20%±20%高注意回流系数和 AT 间距是强耦合的单独调一个可能得出错误结论最好做二维扫描。4.3 和模板.doc的对照验证方法把函数跑出来的结果按模板的假设重新算一遍如果模板忽略 AT 漏抗就把函数里Zat设成 0如果模板用平均电流就把I换成平均值。然后逐项对比电压损失、短路电流、回流系数。差异超过 10% 时优先检查钢轨阻抗是否含大地回流修正、AT 变比是否按 1:1 处理、短路点过渡电阻是否一致。这套对照做完模板就从“黑箱”变成了可解释的计算工具。5. 进阶用MATLAB做AT牵引供电系统的批量仿真与结果可视化5.1 批量跑方案不同AT间距下的电压损失曲线实际设计往往要比较多个 AT 间距方案。把at_traction放进循环一次算出 8、10、12、15 km 四种间距下的电压损失和短路电流直接出图。D_list [8 10 12 15]; dU_list zeros(size(D_list)); for i 1:length(D_list) [dU, ~, ~] at_traction(0.151j*0.55, 0.121j*0.42, ... 0.101j*0.35, D_list(i), 36, 600*exp(-1j*acos(0.95)), 3); dU_list(i) abs(dU); end bar(D_list, dU_list); xlabel(AT间距 (km)); ylabel(电压损失 (V)); title(不同AT间距下的供电臂电压损失);逻辑说明循环调用函数把结果存进数组再画柱状图。参数说明D_list是要比较的间距dU_list存对应损失。从图上能直观看出间距从 8 km 放宽到 15 km 时损失涨了多少从而判断是否值得多设 AT 所。5.2 短路电流的热稳定校验算出最大短路电流后还要校验导线和设备的动热稳定。常见做法是用短路电流有效值Isc和短路持续时间t算热效应Isc^2*t再和导线允许值比。Isc 8000; % A t 0.5; % s Q Isc^2 * t; % A^2·s fprintf(短路热效应 %.2e A^2·s\n, Q);逻辑说明Q是焦耳积分用来查导线热稳定曲线。参数说明t取保护动作时间加断路器开断时间通常 0.3~0.8 s。如果Q超过导线允许值就要加大截面或缩短保护时间。5.3 结果导出与报告生成最后一步是把计算结果写成表格方便贴进设计说明书。MATLAB 的writetable可以直接导出 CSV再用模板里的表格格式整理。T table(D_list, dU_list, VariableNames, {AT间距_km,电压损失_V}); writetable(T, at_result.csv);逻辑说明把两个数组拼成 table 再写文件。参数说明VariableNames用中文列名导出后可直接进报告。跑完这一步一份“计算模板.doc”就真正变成了可复现、可批量、可追溯的计算流程改参数不用再手工重算短路电流和电压损失也能一起出结果。本文还有配套的精品资源点击获取