流线图（`streamline`）

流线图（streamline）

05_场与矢量　5 张　MATLAB R2023b　100 mm 宽 @ 300 dpi

这个方法适合什么

从种子点出发积分出流线，看出场的方向结构。适合：变压器油流、变电站风场、电场线，需要指定种子点时用它。

电气场景

1　变压器油流种子点流线

2　棒板间隙电场线

3　变电站风场种子点流线

4　开关柜气流种子点流线

5　电缆隧道通风流线

代码

脚本 p5_field_vector.m 中对应 streamline 的一段，SEL 用于选择只跑某几个方法。

if any(SEL==7)
    % ---- 7.1 变压器油流种子点流线 ----
    try
        [X,Y,U,V] = gen_oil2(161,121);
        sx = linspace(0.12,3.88,9); sy = 0.10*ones(size(sx));
        verts = stream2(X,Y,U,V,sx,sy,[0.010 20000]);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamline(ax,verts);
        set(h,'Color',C(1,:),'LineWidth',0.65);
        xlabel(ml_tx('油道长度 <i>x</i> / m'));
        ylabel(ml_tx('油道高度 <i>y</i> / m'));
        xlim(ax,[0 4]); ylim(ax,[0 2.4]);
        save_png(f,FIGD,'变压器油流种子点流线');
    catch ME
        fprintf('[FAIL] 7-1 %s\n',ME.message); close all
    end

    % ---- 7.2 棒板间隙电场线 ----
    try
        [X,Y,Z,U,V,W] = gen_rod_E(41,31,25);
        th = linspace(0,2*pi,10); th = th(1:end-1);
        sx = 1.50 + 0.28*cos(th); sy = 1.50 + 0.28*sin(th); sz = 2.00*ones(size(th));
        verts = stream3(X,Y,Z,U,V,W,sx,sy,sz,[0.008 20000]);

        f  = ml_fig(PP,100,90);
        ax = ax3(f,false);
        h = streamline(ax,verts);
        set(h,'Color',C(2,:),'LineWidth',0.65);
        xlabel(ml_tx('横向距离 <i>x</i> / m'));
        ylabel(ml_tx('纵向距离 <i>y</i> / m'));
        zlabel(ml_tx('气隙高度 <i>z</i> / m'));
        xlim(ax,[0 3]); ylim(ax,[0 3]); zlim(ax,[0 2.4]);
        save_png(f,FIGD,'棒板间隙电场线');
    catch ME
        fprintf('[FAIL] 7-2 %s\n',ME.message); close all
    end

    % ---- 7.3 变电站风场种子点流线 ----
    try
        [X,Y,U,V] = gen_wind2(161,121);
        sy = linspace(2,28,9); sx = 0.8*ones(size(sy));
        verts = stream2(X,Y,U,V,sx,sy,[0.05 20000]);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamline(ax,verts);
        set(h,'Color',C(4,:),'LineWidth',0.65);
        xlabel(ml_tx('站区长度 <i>x</i> / m'));
        ylabel(ml_tx('站区宽度 <i>y</i> / m'));
        xlim(ax,[0 40]); ylim(ax,[0 30]);
        save_png(f,FIGD,'变电站风场种子点流线');
    catch ME
        fprintf('[FAIL] 7-3 %s\n',ME.message); close all
    end

    % ---- 7.4 开关柜气流种子点流线 ----
    try
        [X,Y,U,V] = gen_cab2(101,141);
        sx = linspace(0.10,1.10,7); sy = 0.06*ones(size(sx));
        verts = stream2(X,Y,U,V,sx,sy,[0.004 20000]);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamline(ax,verts);
        set(h,'Color',C(3,:),'LineWidth',0.65);
        xlabel(ml_tx('柜宽 <i>x</i> / m'));
        ylabel(ml_tx('柜高 <i>z</i> / m'));
        xlim(ax,[0 1.2]); ylim(ax,[0 1.8]);
        save_png(f,FIGD,'开关柜气流种子点流线');
    catch ME
        fprintf('[FAIL] 7-4 %s\n',ME.message); close all
    end

    % ---- 7.5 电缆隧道通风流线 ----
    try
        [X,Y,Z,U,V,W] = gen_tun(49,17,13);
        [sy,sz] = meshgrid([0.30 0.90],[0.25 0.75]);
        sx = 0.4*ones(size(sy));
        verts = stream3(X,Y,Z,U,V,W,sx,sy,sz,[0.05 20000]);

        f  = ml_fig(PP,100,90);
        ax = ax3(f,false);
        h = streamline(ax,verts);
        set(h,'Color',C(1,:),'LineWidth',0.65);
        xlabel(ml_tx('隧道长度 <i>x</i> / m'));
        ylabel(ml_tx('隧道宽度 <i>y</i> / m'));
        zlabel(ml_tx('高度 <i>z</i> / m'));
        xlim(ax,[0 12]); ylim(ax,[0 1.2]); zlim(ax,[0 1.0]);
        save_png(f,FIGD,'电缆隧道通风流线');
    catch ME
        fprintf('[FAIL] 7-5 %s\n',ME.message); close all
    end
end