简介面向压气机气动设计与叶轮机械方向的学习者和工程人员这份PDF文献围绕多级轴流式压气机展开借助Matlab编程替代繁琐的手工迭代梳理交换叶型设计与平均直径校核的先后次序帮助读者理解如何减少设计循环、提高前期效率。内容按通流尺寸计算与平均直径校核、气动参数确定、截面叶形叠加三步展开并讨论等环量级与等反动度级的径向平衡方程、约束条件选择及相对马赫数沿半径的变化差异截面叠加部分对比前缘叠加与重心叠加的受力特点且给出可导入Pro/E生成三维叶片的空间点坐标思路兼具matlab学习资料与专业参考文献价值。资源包为1个PDF文件约246KB轻量便于查阅。目前已有332人学习适合希望把程序计算与气动设计流程结合起来的中高级读者参考。1. 从一台 4 级轴流式压气机的设计任务说起手里只拿到一组进口条件和性能指标进口总温 288.15 K、总压 101325 Pa、流量 20 kg/s、总压比 4.0要求 12000 r/min 驱动。接下来的活是回答三个问题需要几级、各级速度三角形长什么样、叶片沿叶高怎么扭。多级轴流式压气机的气动设计几乎全部工作都压在这样一个一维骨架上后面再谈三维 CFD 和叶栅试验都是在骨架上去补细节。用 Matlab 做这件事的好处是参数全在脚本里改一个反力度或者转速整条压比链、速度三角形、叶栅稠度能一起重算还能顺手用 matlab 画图把沿叶高分布画出来对比。这篇按设计点核算、逐级速度三角形、稠度与损失、叶片型线、敏感性扫描的顺序走一遍目标是让第一次接触轴流压气机气动设计的人能照着把一维方案算出来也让做过的人交手时能看到参数边界在哪。2. 多级轴流式压气机气动设计的一维参数骨架2.1 设计点输入与总焓升核算一维设计从总焓升下手而不是从几何下手。轴流压气机每一级做的是对气流加功加功量直接由总温升反映几何只是实现手段。所以第一步是把压比和效率折算成总温升再把总温升按级数摊开。这个顺序不能反先画叶型再算压比是没法收敛的。设计点输入我一般固定成一张表脚本里也按这张表定义常量避免后面改参数时漏掉某个量。参数符号取值说明进口总温T01288.15 K标准状态进口总压p01101325 Pa标准状态质量流量mdot20 kg/s设计点总压比pr4.0设计点绝热效率eta0.88初值后续迭代转速n12000 r/min转速固定轴向速度ca150 m/s各级近似常值级数z4先估后校轴向速度取 150 m/s 是一个折中值。取得太高叶片相对马赫数上去了激波损失压不住取得太低通道面积变大前几级叶片过长叶尖切线速度又不够。工程上中小型机组 130170 m/s 是常见区间。% stage_enthalpy.m 设计点总焓升核算与级间分配 clear; clc cp 1005; % 定压比热 J/(kg*K) gam 1.4; % 比热比 Rg 287; % 气体常数 J/(kg*K) T01 288.15; p01 101325; mdot 20; pr 4.0; eta 0.88; n 12000; ca 150; z 4; % 等熵温升除以效率得到真实总温升 dT0_tot T01*(pr^((gam-1)/gam) - 1)/eta; dh0_tot cp*dT0_tot; dh0 dh0_tot/z; % 等焓升分配 T0 T01 (0:z)*dh0/cp; % 各级进口总温长度 z1 p0 zeros(1,z1); p0(1) p01; for k 1:z % 逐级按级效率缩放先按 0.88 试算 p0(k1) p0(k)*(1 eta*dh0/(cp*T0(k)))^(gam/(gam-1)); end fprintf(总温升 %.2f K每级焓升 %.2f kJ/kg\n, dT0_tot, dh0/1000); fprintf(出口总压 %.0f Pa实际压比 %.3f\n, p0(end), p0(end)/p01);等焓升分配是最简单的分法实际常用的是前几级多分一点、后几级少分一点原因是后面级气流密度大、叶高短二次流损失占比高。等焓升先跑通再切成递增或递减分布只需要把dh0从标量改成向量dh0 w.*dh0_tot权重w归一化即可。按上面参数跑出来总温升约 159.2 K每级约 39.8 kJ/kg出口总压回到 405000 Pa 附近说明分配自洽。到这里还没有任何几何信息但整机的能量账已经平了。2.2 载荷系数、流量系数与反力度的解析关系级几何的核心是速度三角形。对常轴向速度、进口无预旋的级三个无量纲量就能定形载荷系数 ψ Δh0/U²流量系数 φ ca/U反力度 R 1 − (cθ1 cθ2)/(2U)。其中 U 是级平均半径处的切线速度cθ 是绝对速度切向分量。取级平均直径 0.42 m12000 r/min 对应 U ≈ 263.9 m/s于是 ψ ≈ 0.574φ ≈ 0.568R ≈ 0.713。ψ 落在 0.50.6 属于单级压比 1.4 左右的常规水平R 到 0.71 偏高意味着静叶承担的扩压任务轻、动叶承担的扩压任务重动叶吸力面逆压梯度大容易提前分离。想把反力度压到 0.550.65最直接的办法是给进口加预旋。由 R 和 Δcθ 反解% velocity_triangle.m 由 R 与载荷系数反解速度三角形 Dm 0.42; % 级平均直径 m U pi*Dm*n/60; % 切线速度 m/s psi dh0/U^2; % 载荷系数 phi ca/U; % 流量系数 dCth psi*U; % Δcθ Δh0/U R 0.62; % 目标反力度 cth1 U*(1-R) - dCth/2; % 进口绝对切向分速 cth2 U*(1-R) dCth/2; % 出口绝对切向分速 beta1 atand((U - cth1)/ca); % 动叶进口相对气流角 beta2 atand((U - cth2)/ca); % 动叶出口相对气流角 alpha2 atand(cth2/ca); % 静叶进口绝对气流角 fprintf(U%.1f m/s, psi%.3f, phi%.3f\n, U, psi, phi); fprintf(beta1%.1f deg, beta2%.1f deg, alpha2%.1f deg\n, beta1, beta2, alpha2);几个输出要盯住beta1 一般在 55°65°beta2 在 30°45°两者之差就是动叶弯角超过 30° 就要警惕。alpha2 超过 50° 时静叶进口过于偏转扩压任务集中在静叶前段。这三条边界比任何优化目标都先用得上先卡边界再谈效率。2.3 级数分配与自由涡沿叶高规律级数不是拍脑袋定的。每级焓升除以 U² 得到的 ψ 超过 0.6 就要加级因为再往上加载荷损失会陡增。反过来 ψ 掉到 0.35 以下说明级数多了重量和成本白放。前面算 ψ 0.574 落在区间内4 级是合理的但注意这是级平均半径处的值叶根处 U 小、ψ 大叶尖处 U 大、ψ 小真正卡住的往往是叶根。沿叶高规律常用自由涡设计即 cθ·r 常数。这时反力度沿叶高按R(r) 1 − (cθ1m cθ2m)·rm / (2·ω·r²)分布半径越大反力度越大。按叶根 0.175 m、平均 0.21 m、叶尖 0.245 m 代入叶根反力度约 0.59叶尖约 0.79。叶根反力度偏低是自由涡设计的典型特征如果叶根算出来接近 0 甚至为负说明载荷系数过载要么加级要么降低 Δh0 的径向不均匀度。% free_vortex.m 自由涡沿叶高反力度与气流角分布 rm 0.21; Rm 0.62; omega 2*pi*n/60; r linspace(0.175, 0.245, 9); Rr 1 - (1-Rm)*(rm./r).^2; % 自由涡反力度分布 Ur omega*r; % 各半径切线速度 psi_r dh0./Ur.^2; alpha1_r atand((1-Rr).*Ur/ca*2 - atand(cth2/ca)); % 进口绝对角 disp(table(r, Ur, Rr, psi_r, VariableNames, {r_m,U_ms,R,psi}));这段的关键是(rm./r).^2这一项它决定了反力度径向变化有多剧烈。轮毂比越小叶高越大这项影响越强叶根越危险。中小型机组轮毂比常在 0.50.7叶高比大的低压级更容易在这儿翻车。3. Matlab 逐级速度三角形求解与级叠加计算3.1 用 struct 数组组织多级压气机的级参数单级速度三角形算一遍不难难的是让每一级都能独立调整参数、又能自动串成整机。用 struct 数组按级存参数最省事每级记录平均直径、焓升、反力度、总温总压进出口值循环里只读不写全局变量改哪一级改哪一级。% build_stages.m 级参数结构体初始化 function stg build_stages(z, dh0, T01, p01, Dm0, dDm) % z 级数 % dh0 每级焓升 J/kg标量或长度 z 向量 % Dm0 第一级平均直径 m % dDm 平均直径逐级缩小的量 m if isscalar(dh0), dh0 repmat(dh0, 1, z); end for k 1:z stg(k).id k; stg(k).Dm Dm0 - (k-1)*dDm; stg(k).dh0 dh0(k); stg(k).R 0.62; % 反力度初值 stg(k).T0in T01; % 占位循环里更新 stg(k).p0in p01; end end平均直径逐级缩小是通道收缩的一种实现方式。另一条路是保持平均直径不变、只缩叶高结构上更简单但叶根切线速度不变后几级的载荷系数会偏高。两种做法的取舍点在加工成本和气动收益之间方案阶段先用缩小直径结构评审时再定。3.2 主循环从进口到出口的逐级叠加级叠加stage stacking是这个脚本的核心每一级吃前一级的出口总温总压算自己的速度三角形、损失和出口状态把结果交出去。循环体里把所有随半径变化的量都留在平均半径上算径向分布留到后面单独处理。% stage_stacking.m 逐级叠加求解 cp 1005; gam 1.4; Rg 287; n 12000; ca 150; omega 2*pi*n/60; z 4; dh0 39986; % 每级升 J/kg stg build_stages(z, dh0, 288.15, 101325, 0.42, 0.008); for k 1:z U omega*stg(k).Dm/2; psi stg(k).dh0/U^2; dCth stg(k).dh0/U; cth1 U*(1-stg(k).R) - dCth/2; cth2 U*(1-stg(k).R) dCth/2; w1 sqrt(ca^2 (U-cth1)^2); % 动叶进口相对速度 w2 sqrt(ca^2 (U-cth2)^2); % 动叶出口相对速度 stg(k).sigma dCth/(2*w1*(0.6 - 1 w2/w1)); % 由 D0.6 反推稠度 % 损失系数叶型项随扩散因子增长二次流项随叶高弦长比增长 D 1 - w2/w1 dCth/(2*stg(k).sigma*w1); wprof 0.004 0.008*(D/0.6)^2; wsec 0.018*(0.04/(0.035*stg(k).sigma)); % 弦长与叶高比 wloss wprof wsec; stg(k).T0out stg(k).T0in stg(k).dh0/cp; p0id stg(k).p0in*(stg(k).T0out/stg(k).T0in)^(gam/(gam-1)); stg(k).p0out p0id*(1 - wloss); stg(k).eta ((stg(k).p0out/stg(k).p0in)^((gam-1)/gam) - 1) / ... (stg(k).T0out/stg(k).T0in - 1); if k z stg(k1).T0in stg(k).T0out; stg(k1).p0in stg(k).p0out; end end三个参数决定这段代码的输出形态。R控制速度三角形的形状调它就是在动叶和静叶之间搬扩压任务0.6是扩散因子上限调小会让稠度增大、叶栅变密wprof与wsec的系数是方案阶段的工程拟合值用来快速估效率定型阶段必须换成叶栅试验数据或者经过校核的损失模型否则算出来的效率只能当趋势看。build_stages里第一级入口总温总压写死后面级在前一级出口赋值这个链条不能断。如果发现某级出口总压低于进口先查stg(k).R是不是被改到了 0.9 以上反力度过高会把动叶压得喘不过气。3.3 结果校验流量、压比与稠度是否自洽算完不等于对了。第一件事是核对总压比是不是回到 4.0 附近第二件事是按质量流量反算各级通道面积看叶高是不是递减得合理。级平均直径 mU m/sψ稠度出口总压 Pa级效率10.420263.90.5741.151418000.89120.412258.90.5971.191980000.88330.404253.80.6211.242766000.87540.396248.80.6461.303865000.867压比 3.81比目标 4.0 差了一点说明损失估得偏保守典型处理是把反力度往低压一档或者把稠度调密一点让二次流损失下降再迭代一轮。载荷系数逐级上升是因为平均直径在缩这时第二级的 ψ 已经接近 0.6把dDm从 0.008 降到 0.005 会明显缓解。面积校核按% area_check.m 按连续性与状态方程反算通道面积 Rg 287; for k 1:z T stg(k).T0in - ca^2/(2*cp); p stg(k).p0in*(T/stg(k).T0in)^(gam/(gam-1)); rho p/(Rg*T); A 20/(rho*ca); fprintf(第%d级 面积 %.4f m2 静温 %.1f K\n, k, A, T); end面积逐级减小是正常的缩小比例大致等于密度的上升比例。如果算出来某一级面积反而变大说明轴向速度在那级被写错了或者总温分配出现了倒退。4. 从速度三角形到叶栅稠度、损失与沿叶高扭曲4.1 Lieblein 扩散因子反推叶栅稠度速度三角形只给出气流角叶片要做多密、弦长多长由稠度 σ c/s 决定。工程上用 Lieblein 扩散因子把稠度和分离风险连起来D 1 − w2/w1 Δcθ/(2σ·w1)D 经验上限取 0.6超过就容易在动叶吸力面后段分离。把 D 0.6 代进去解 σ得到的稠度就是满足扩压要求的最小值实际取值再往上留 10% 余量。第一级算下来 σ ≈ 1.15落到轴流压气机常见的 1.01.5 区间。稠度定下来弦长由栅距反推。栅距 s 2πr/Z叶片数 Z 先按经验给一个数比如 40得到 s ≈ 33 mm弦长 c σ·s ≈ 38 mm。叶片数不是随便给的它同时影响栅距和叶栅损失一般让 s/c 落在 0.71.1太密则摩擦损失上升、太重太稀则单叶负荷过大。提示稠度、弦长、叶片数这三个量在方案阶段互相耦合改一个另外两个都得跟着动。建议把叶片数做成外层循环变量稠度由扩散因子反推弦长作为输出观察别反过来先定弦长。4.2 损失系数表与逐级效率修正方案阶段用一组简化损失系数能快速给出效率趋势定型阶段必须换成经过校核的模型。常用的分解方式是叶型损失、二次流损失、叶尖间隙损失三项相加。损失项经验表达式量级主要影响因素叶型损失0.004 0.008(D/0.6)²0.010.03扩散因子、冲角二次流损失0.018·c/h0.010.04弦长叶高比、端壁叶尖间隙损失0.02·τ/h0.010.03相对间隙、叶高叶尖间隙损失随叶高减小而放大所以后几级的效率天然比前几级低这也是为什么级焓升不宜等分——后级分少一点效率损失反而摊得开。把三项损失加起来得到总压损失系数再按级出口理想总压乘上 (1 − ω) 得到实际总压级效率就出来了。% stage_efficiency.m 三项损失合成与级效率 function [eta, wloss] stage_efficiency(D, sigma, c, h, tau) wprof 0.004 0.008*(D/0.6)^2; wsec 0.018*(c/h); wtip 0.02*(tau/h); wloss wprof wsec wtip; eta []; % 由主循环用总温总压计算 endc/h是弦长与叶高比后几级叶高掉到 20 mm 以下时这项会明显抬头此时可以通过减小弦长加大叶片数来压一压但叶片数增加又会推高摩擦损失需要权衡。相对间隙tau/h取决于加工和装配水平方案阶段取 0.010.02 之间的值比较稳妥。4.3 中弧线与厚度分布生成叶片型线速度三角形和气流角最终要落成型线。工程上做方案设计不需要全三维造型用中弧线叠加厚度分布就能得到可用于初步评估的叶型。中弧线用单圆弧近似给定弦长 c 和弯角 θ由 beta1 与 beta2 之差扣掉冲角和落后角得到圆心在弦中垂线上% blade_profile.m 圆弧中弧线 NACA 厚度分布 function [xu, yu, xl, yl] blade_profile(c, theta, t, N) % c 弦长 m % theta 弯角 rad % t 最大相对厚度例如 0.08 % N 单侧离散点数 Rc c/(2*sin(theta/2)); % 中弧线圆弧半径 d Rc*cos(theta/2); % 圆心到弦的距离 phi linspace(-theta/2, theta/2, N); xc c/2 Rc*sin(phi); yc Rc*(cos(phi) - cos(theta/2)); % 最大弯度在 phi0 % 中弧线切线角用于沿法线叠加厚度 dx gradient(xc); dy gradient(yc); alpha atan2(dy, dx); xn linspace(0, 1, N); yt (t*c/0.2)*(0.2969*sqrt(xn) - 0.1260*xn - 0.3516*xn.^2 ... 0.2843*xn.^3 - 0.1015*xn.^4); xu xc - yt.*sin(alpha); yu yc yt.*cos(alpha); xl xc yt.*sin(alpha); yl yc - yt.*cos(alpha); endtheta的取值来自速度三角形动叶的弯角大致是 beta1 减 beta2但要扣掉进口冲角 i通常 −2°2°和出口落后角 δ。落后角常用 Carter 关系估算δ m·θ·σ^(−1/2)m 在压气机叶栅中取 0.20.25。薄壁圆弧中弧线适合大弯角、低马赫数高亚声速级更适合双圆弧或者多圆弧中弧线把弯度往通道后段搬降低前段加速区的峰值马赫数。厚度分布用 NACA 四位翼型的厚度公式最大相对厚度从叶根的 0.10 到叶尖的 0.06 递减兼顾结构强度和气动性能。生成的型线可以直接用patch或fill画出来对照检查前后缘是否闭合、厚度是否出现负值——sqrt(xn)在 xn 0 处导数为无穷离散点不够密时前缘会翘起来把 N 从 50 提到 200 就能压下去。5. 设计点校核与参数敏感性扫描的几个实用技巧方案跑通一次不算完真正决定成败的是偏离设计点还能不能站住。多级轴流压气机最怕的是前面级先堵、后面级先喘一维阶段能做的就是扫参数看趋势。把转速和反力度做成外层循环每换一组算一次级叠加看效率、扩散因子和反力度分布怎么变。% sensitivity_sweep.m 转速与反力度双参数扫描 nList 10800:300:13200; % 转速扫描%N 从 90 到 110 RList 0.55:0.02:0.70; effGrid zeros(numel(nList), numel(RList)); DmaxGrid zeros(size(effGrid)); for i 1:numel(nList) for j 1:numel(RList) [etaVec, Dvec] run_stage_stacking(nList(i), RList(j)); effGrid(i,j) prod(etaVec)^(1/numel(etaVec)); % 整机效率 DmaxGrid(i,j) max(Dvec); % 最大扩散因子 end end % 落在 D0.6 的区域才算可行再从中挑效率最高点 feasible DmaxGrid 0.60; effFeasible effGrid; effFeasible(~feasible) NaN; [bestEff, idx] max(effFeasible(:)); [ni, rj] ind2sub(size(effFeasible), idx); fprintf(可行域内最优转速 %.0f r/min反力度 %.2f整机效率 %.4f\n, ... nList(ni), RList(rj), bestEff);这一步的价值在于把可行域画出来而不是找那个最优点。单点最优的效率值受损失模型精度限制可信度有限可行域的形状却相当稳定能告诉你反力度调到 0.68 以上扩散因子立刻越界转速掉到 10500 r/min 以下时第一级载荷系数会冲上 0.7。这些边界才是交给结构和总体去谈的依据。如果想把叶片数、反力度、弦长一起调靠嵌套循环效率太低可以交给 matlab 优化工具箱里的fmincon把整机效率作为目标函数扩散因子、反力度下限、叶片数取整作为约束。写目标函数时注意把所有物理量归一到同一量级否则求解器在叶片数这种量级为 10 的变量和反力度这种量级为 1 的变量之间会走不动步这也是很多人用 matlab 优化工具箱做参数寻优时最先踩的坑。画图这一步别省。用imagesc把effFeasible画出来加上colorbar用plot把每一级的反力度沿叶高叠在一张图上看叶根那条线有没有贴到 0.4 以下。图形化之后的判断速度比看数字表快一个数量级特别是要同时比较三四种级数方案的时候。用 matlab 画图把设计点和两个典型非设计点叠在一起喘振边界和堵塞边界的相对位置一眼就能看出来再回头调第一级的稠度和叶尖间隙方向就明确了。本文还有配套的精品资源点击获取
Matlab多级轴流压气机气动设计：一维参数、速度三角形与级叠加
简介面向压气机气动设计与叶轮机械方向的学习者和工程人员这份PDF文献围绕多级轴流式压气机展开借助Matlab编程替代繁琐的手工迭代梳理交换叶型设计与平均直径校核的先后次序帮助读者理解如何减少设计循环、提高前期效率。内容按通流尺寸计算与平均直径校核、气动参数确定、截面叶形叠加三步展开并讨论等环量级与等反动度级的径向平衡方程、约束条件选择及相对马赫数沿半径的变化差异截面叠加部分对比前缘叠加与重心叠加的受力特点且给出可导入Pro/E生成三维叶片的空间点坐标思路兼具matlab学习资料与专业参考文献价值。资源包为1个PDF文件约246KB轻量便于查阅。目前已有332人学习适合希望把程序计算与气动设计流程结合起来的中高级读者参考。1. 从一台 4 级轴流式压气机的设计任务说起手里只拿到一组进口条件和性能指标进口总温 288.15 K、总压 101325 Pa、流量 20 kg/s、总压比 4.0要求 12000 r/min 驱动。接下来的活是回答三个问题需要几级、各级速度三角形长什么样、叶片沿叶高怎么扭。多级轴流式压气机的气动设计几乎全部工作都压在这样一个一维骨架上后面再谈三维 CFD 和叶栅试验都是在骨架上去补细节。用 Matlab 做这件事的好处是参数全在脚本里改一个反力度或者转速整条压比链、速度三角形、叶栅稠度能一起重算还能顺手用 matlab 画图把沿叶高分布画出来对比。这篇按设计点核算、逐级速度三角形、稠度与损失、叶片型线、敏感性扫描的顺序走一遍目标是让第一次接触轴流压气机气动设计的人能照着把一维方案算出来也让做过的人交手时能看到参数边界在哪。2. 多级轴流式压气机气动设计的一维参数骨架2.1 设计点输入与总焓升核算一维设计从总焓升下手而不是从几何下手。轴流压气机每一级做的是对气流加功加功量直接由总温升反映几何只是实现手段。所以第一步是把压比和效率折算成总温升再把总温升按级数摊开。这个顺序不能反先画叶型再算压比是没法收敛的。设计点输入我一般固定成一张表脚本里也按这张表定义常量避免后面改参数时漏掉某个量。参数符号取值说明进口总温T01288.15 K标准状态进口总压p01101325 Pa标准状态质量流量mdot20 kg/s设计点总压比pr4.0设计点绝热效率eta0.88初值后续迭代转速n12000 r/min转速固定轴向速度ca150 m/s各级近似常值级数z4先估后校轴向速度取 150 m/s 是一个折中值。取得太高叶片相对马赫数上去了激波损失压不住取得太低通道面积变大前几级叶片过长叶尖切线速度又不够。工程上中小型机组 130170 m/s 是常见区间。% stage_enthalpy.m 设计点总焓升核算与级间分配 clear; clc cp 1005; % 定压比热 J/(kg*K) gam 1.4; % 比热比 Rg 287; % 气体常数 J/(kg*K) T01 288.15; p01 101325; mdot 20; pr 4.0; eta 0.88; n 12000; ca 150; z 4; % 等熵温升除以效率得到真实总温升 dT0_tot T01*(pr^((gam-1)/gam) - 1)/eta; dh0_tot cp*dT0_tot; dh0 dh0_tot/z; % 等焓升分配 T0 T01 (0:z)*dh0/cp; % 各级进口总温长度 z1 p0 zeros(1,z1); p0(1) p01; for k 1:z % 逐级按级效率缩放先按 0.88 试算 p0(k1) p0(k)*(1 eta*dh0/(cp*T0(k)))^(gam/(gam-1)); end fprintf(总温升 %.2f K每级焓升 %.2f kJ/kg\n, dT0_tot, dh0/1000); fprintf(出口总压 %.0f Pa实际压比 %.3f\n, p0(end), p0(end)/p01);等焓升分配是最简单的分法实际常用的是前几级多分一点、后几级少分一点原因是后面级气流密度大、叶高短二次流损失占比高。等焓升先跑通再切成递增或递减分布只需要把dh0从标量改成向量dh0 w.*dh0_tot权重w归一化即可。按上面参数跑出来总温升约 159.2 K每级约 39.8 kJ/kg出口总压回到 405000 Pa 附近说明分配自洽。到这里还没有任何几何信息但整机的能量账已经平了。2.2 载荷系数、流量系数与反力度的解析关系级几何的核心是速度三角形。对常轴向速度、进口无预旋的级三个无量纲量就能定形载荷系数 ψ Δh0/U²流量系数 φ ca/U反力度 R 1 − (cθ1 cθ2)/(2U)。其中 U 是级平均半径处的切线速度cθ 是绝对速度切向分量。取级平均直径 0.42 m12000 r/min 对应 U ≈ 263.9 m/s于是 ψ ≈ 0.574φ ≈ 0.568R ≈ 0.713。ψ 落在 0.50.6 属于单级压比 1.4 左右的常规水平R 到 0.71 偏高意味着静叶承担的扩压任务轻、动叶承担的扩压任务重动叶吸力面逆压梯度大容易提前分离。想把反力度压到 0.550.65最直接的办法是给进口加预旋。由 R 和 Δcθ 反解% velocity_triangle.m 由 R 与载荷系数反解速度三角形 Dm 0.42; % 级平均直径 m U pi*Dm*n/60; % 切线速度 m/s psi dh0/U^2; % 载荷系数 phi ca/U; % 流量系数 dCth psi*U; % Δcθ Δh0/U R 0.62; % 目标反力度 cth1 U*(1-R) - dCth/2; % 进口绝对切向分速 cth2 U*(1-R) dCth/2; % 出口绝对切向分速 beta1 atand((U - cth1)/ca); % 动叶进口相对气流角 beta2 atand((U - cth2)/ca); % 动叶出口相对气流角 alpha2 atand(cth2/ca); % 静叶进口绝对气流角 fprintf(U%.1f m/s, psi%.3f, phi%.3f\n, U, psi, phi); fprintf(beta1%.1f deg, beta2%.1f deg, alpha2%.1f deg\n, beta1, beta2, alpha2);几个输出要盯住beta1 一般在 55°65°beta2 在 30°45°两者之差就是动叶弯角超过 30° 就要警惕。alpha2 超过 50° 时静叶进口过于偏转扩压任务集中在静叶前段。这三条边界比任何优化目标都先用得上先卡边界再谈效率。2.3 级数分配与自由涡沿叶高规律级数不是拍脑袋定的。每级焓升除以 U² 得到的 ψ 超过 0.6 就要加级因为再往上加载荷损失会陡增。反过来 ψ 掉到 0.35 以下说明级数多了重量和成本白放。前面算 ψ 0.574 落在区间内4 级是合理的但注意这是级平均半径处的值叶根处 U 小、ψ 大叶尖处 U 大、ψ 小真正卡住的往往是叶根。沿叶高规律常用自由涡设计即 cθ·r 常数。这时反力度沿叶高按R(r) 1 − (cθ1m cθ2m)·rm / (2·ω·r²)分布半径越大反力度越大。按叶根 0.175 m、平均 0.21 m、叶尖 0.245 m 代入叶根反力度约 0.59叶尖约 0.79。叶根反力度偏低是自由涡设计的典型特征如果叶根算出来接近 0 甚至为负说明载荷系数过载要么加级要么降低 Δh0 的径向不均匀度。% free_vortex.m 自由涡沿叶高反力度与气流角分布 rm 0.21; Rm 0.62; omega 2*pi*n/60; r linspace(0.175, 0.245, 9); Rr 1 - (1-Rm)*(rm./r).^2; % 自由涡反力度分布 Ur omega*r; % 各半径切线速度 psi_r dh0./Ur.^2; alpha1_r atand((1-Rr).*Ur/ca*2 - atand(cth2/ca)); % 进口绝对角 disp(table(r, Ur, Rr, psi_r, VariableNames, {r_m,U_ms,R,psi}));这段的关键是(rm./r).^2这一项它决定了反力度径向变化有多剧烈。轮毂比越小叶高越大这项影响越强叶根越危险。中小型机组轮毂比常在 0.50.7叶高比大的低压级更容易在这儿翻车。3. Matlab 逐级速度三角形求解与级叠加计算3.1 用 struct 数组组织多级压气机的级参数单级速度三角形算一遍不难难的是让每一级都能独立调整参数、又能自动串成整机。用 struct 数组按级存参数最省事每级记录平均直径、焓升、反力度、总温总压进出口值循环里只读不写全局变量改哪一级改哪一级。% build_stages.m 级参数结构体初始化 function stg build_stages(z, dh0, T01, p01, Dm0, dDm) % z 级数 % dh0 每级焓升 J/kg标量或长度 z 向量 % Dm0 第一级平均直径 m % dDm 平均直径逐级缩小的量 m if isscalar(dh0), dh0 repmat(dh0, 1, z); end for k 1:z stg(k).id k; stg(k).Dm Dm0 - (k-1)*dDm; stg(k).dh0 dh0(k); stg(k).R 0.62; % 反力度初值 stg(k).T0in T01; % 占位循环里更新 stg(k).p0in p01; end end平均直径逐级缩小是通道收缩的一种实现方式。另一条路是保持平均直径不变、只缩叶高结构上更简单但叶根切线速度不变后几级的载荷系数会偏高。两种做法的取舍点在加工成本和气动收益之间方案阶段先用缩小直径结构评审时再定。3.2 主循环从进口到出口的逐级叠加级叠加stage stacking是这个脚本的核心每一级吃前一级的出口总温总压算自己的速度三角形、损失和出口状态把结果交出去。循环体里把所有随半径变化的量都留在平均半径上算径向分布留到后面单独处理。% stage_stacking.m 逐级叠加求解 cp 1005; gam 1.4; Rg 287; n 12000; ca 150; omega 2*pi*n/60; z 4; dh0 39986; % 每级升 J/kg stg build_stages(z, dh0, 288.15, 101325, 0.42, 0.008); for k 1:z U omega*stg(k).Dm/2; psi stg(k).dh0/U^2; dCth stg(k).dh0/U; cth1 U*(1-stg(k).R) - dCth/2; cth2 U*(1-stg(k).R) dCth/2; w1 sqrt(ca^2 (U-cth1)^2); % 动叶进口相对速度 w2 sqrt(ca^2 (U-cth2)^2); % 动叶出口相对速度 stg(k).sigma dCth/(2*w1*(0.6 - 1 w2/w1)); % 由 D0.6 反推稠度 % 损失系数叶型项随扩散因子增长二次流项随叶高弦长比增长 D 1 - w2/w1 dCth/(2*stg(k).sigma*w1); wprof 0.004 0.008*(D/0.6)^2; wsec 0.018*(0.04/(0.035*stg(k).sigma)); % 弦长与叶高比 wloss wprof wsec; stg(k).T0out stg(k).T0in stg(k).dh0/cp; p0id stg(k).p0in*(stg(k).T0out/stg(k).T0in)^(gam/(gam-1)); stg(k).p0out p0id*(1 - wloss); stg(k).eta ((stg(k).p0out/stg(k).p0in)^((gam-1)/gam) - 1) / ... (stg(k).T0out/stg(k).T0in - 1); if k z stg(k1).T0in stg(k).T0out; stg(k1).p0in stg(k).p0out; end end三个参数决定这段代码的输出形态。R控制速度三角形的形状调它就是在动叶和静叶之间搬扩压任务0.6是扩散因子上限调小会让稠度增大、叶栅变密wprof与wsec的系数是方案阶段的工程拟合值用来快速估效率定型阶段必须换成叶栅试验数据或者经过校核的损失模型否则算出来的效率只能当趋势看。build_stages里第一级入口总温总压写死后面级在前一级出口赋值这个链条不能断。如果发现某级出口总压低于进口先查stg(k).R是不是被改到了 0.9 以上反力度过高会把动叶压得喘不过气。3.3 结果校验流量、压比与稠度是否自洽算完不等于对了。第一件事是核对总压比是不是回到 4.0 附近第二件事是按质量流量反算各级通道面积看叶高是不是递减得合理。级平均直径 mU m/sψ稠度出口总压 Pa级效率10.420263.90.5741.151418000.89120.412258.90.5971.191980000.88330.404253.80.6211.242766000.87540.396248.80.6461.303865000.867压比 3.81比目标 4.0 差了一点说明损失估得偏保守典型处理是把反力度往低压一档或者把稠度调密一点让二次流损失下降再迭代一轮。载荷系数逐级上升是因为平均直径在缩这时第二级的 ψ 已经接近 0.6把dDm从 0.008 降到 0.005 会明显缓解。面积校核按% area_check.m 按连续性与状态方程反算通道面积 Rg 287; for k 1:z T stg(k).T0in - ca^2/(2*cp); p stg(k).p0in*(T/stg(k).T0in)^(gam/(gam-1)); rho p/(Rg*T); A 20/(rho*ca); fprintf(第%d级 面积 %.4f m2 静温 %.1f K\n, k, A, T); end面积逐级减小是正常的缩小比例大致等于密度的上升比例。如果算出来某一级面积反而变大说明轴向速度在那级被写错了或者总温分配出现了倒退。4. 从速度三角形到叶栅稠度、损失与沿叶高扭曲4.1 Lieblein 扩散因子反推叶栅稠度速度三角形只给出气流角叶片要做多密、弦长多长由稠度 σ c/s 决定。工程上用 Lieblein 扩散因子把稠度和分离风险连起来D 1 − w2/w1 Δcθ/(2σ·w1)D 经验上限取 0.6超过就容易在动叶吸力面后段分离。把 D 0.6 代进去解 σ得到的稠度就是满足扩压要求的最小值实际取值再往上留 10% 余量。第一级算下来 σ ≈ 1.15落到轴流压气机常见的 1.01.5 区间。稠度定下来弦长由栅距反推。栅距 s 2πr/Z叶片数 Z 先按经验给一个数比如 40得到 s ≈ 33 mm弦长 c σ·s ≈ 38 mm。叶片数不是随便给的它同时影响栅距和叶栅损失一般让 s/c 落在 0.71.1太密则摩擦损失上升、太重太稀则单叶负荷过大。提示稠度、弦长、叶片数这三个量在方案阶段互相耦合改一个另外两个都得跟着动。建议把叶片数做成外层循环变量稠度由扩散因子反推弦长作为输出观察别反过来先定弦长。4.2 损失系数表与逐级效率修正方案阶段用一组简化损失系数能快速给出效率趋势定型阶段必须换成经过校核的模型。常用的分解方式是叶型损失、二次流损失、叶尖间隙损失三项相加。损失项经验表达式量级主要影响因素叶型损失0.004 0.008(D/0.6)²0.010.03扩散因子、冲角二次流损失0.018·c/h0.010.04弦长叶高比、端壁叶尖间隙损失0.02·τ/h0.010.03相对间隙、叶高叶尖间隙损失随叶高减小而放大所以后几级的效率天然比前几级低这也是为什么级焓升不宜等分——后级分少一点效率损失反而摊得开。把三项损失加起来得到总压损失系数再按级出口理想总压乘上 (1 − ω) 得到实际总压级效率就出来了。% stage_efficiency.m 三项损失合成与级效率 function [eta, wloss] stage_efficiency(D, sigma, c, h, tau) wprof 0.004 0.008*(D/0.6)^2; wsec 0.018*(c/h); wtip 0.02*(tau/h); wloss wprof wsec wtip; eta []; % 由主循环用总温总压计算 endc/h是弦长与叶高比后几级叶高掉到 20 mm 以下时这项会明显抬头此时可以通过减小弦长加大叶片数来压一压但叶片数增加又会推高摩擦损失需要权衡。相对间隙tau/h取决于加工和装配水平方案阶段取 0.010.02 之间的值比较稳妥。4.3 中弧线与厚度分布生成叶片型线速度三角形和气流角最终要落成型线。工程上做方案设计不需要全三维造型用中弧线叠加厚度分布就能得到可用于初步评估的叶型。中弧线用单圆弧近似给定弦长 c 和弯角 θ由 beta1 与 beta2 之差扣掉冲角和落后角得到圆心在弦中垂线上% blade_profile.m 圆弧中弧线 NACA 厚度分布 function [xu, yu, xl, yl] blade_profile(c, theta, t, N) % c 弦长 m % theta 弯角 rad % t 最大相对厚度例如 0.08 % N 单侧离散点数 Rc c/(2*sin(theta/2)); % 中弧线圆弧半径 d Rc*cos(theta/2); % 圆心到弦的距离 phi linspace(-theta/2, theta/2, N); xc c/2 Rc*sin(phi); yc Rc*(cos(phi) - cos(theta/2)); % 最大弯度在 phi0 % 中弧线切线角用于沿法线叠加厚度 dx gradient(xc); dy gradient(yc); alpha atan2(dy, dx); xn linspace(0, 1, N); yt (t*c/0.2)*(0.2969*sqrt(xn) - 0.1260*xn - 0.3516*xn.^2 ... 0.2843*xn.^3 - 0.1015*xn.^4); xu xc - yt.*sin(alpha); yu yc yt.*cos(alpha); xl xc yt.*sin(alpha); yl yc - yt.*cos(alpha); endtheta的取值来自速度三角形动叶的弯角大致是 beta1 减 beta2但要扣掉进口冲角 i通常 −2°2°和出口落后角 δ。落后角常用 Carter 关系估算δ m·θ·σ^(−1/2)m 在压气机叶栅中取 0.20.25。薄壁圆弧中弧线适合大弯角、低马赫数高亚声速级更适合双圆弧或者多圆弧中弧线把弯度往通道后段搬降低前段加速区的峰值马赫数。厚度分布用 NACA 四位翼型的厚度公式最大相对厚度从叶根的 0.10 到叶尖的 0.06 递减兼顾结构强度和气动性能。生成的型线可以直接用patch或fill画出来对照检查前后缘是否闭合、厚度是否出现负值——sqrt(xn)在 xn 0 处导数为无穷离散点不够密时前缘会翘起来把 N 从 50 提到 200 就能压下去。5. 设计点校核与参数敏感性扫描的几个实用技巧方案跑通一次不算完真正决定成败的是偏离设计点还能不能站住。多级轴流压气机最怕的是前面级先堵、后面级先喘一维阶段能做的就是扫参数看趋势。把转速和反力度做成外层循环每换一组算一次级叠加看效率、扩散因子和反力度分布怎么变。% sensitivity_sweep.m 转速与反力度双参数扫描 nList 10800:300:13200; % 转速扫描%N 从 90 到 110 RList 0.55:0.02:0.70; effGrid zeros(numel(nList), numel(RList)); DmaxGrid zeros(size(effGrid)); for i 1:numel(nList) for j 1:numel(RList) [etaVec, Dvec] run_stage_stacking(nList(i), RList(j)); effGrid(i,j) prod(etaVec)^(1/numel(etaVec)); % 整机效率 DmaxGrid(i,j) max(Dvec); % 最大扩散因子 end end % 落在 D0.6 的区域才算可行再从中挑效率最高点 feasible DmaxGrid 0.60; effFeasible effGrid; effFeasible(~feasible) NaN; [bestEff, idx] max(effFeasible(:)); [ni, rj] ind2sub(size(effFeasible), idx); fprintf(可行域内最优转速 %.0f r/min反力度 %.2f整机效率 %.4f\n, ... nList(ni), RList(rj), bestEff);这一步的价值在于把可行域画出来而不是找那个最优点。单点最优的效率值受损失模型精度限制可信度有限可行域的形状却相当稳定能告诉你反力度调到 0.68 以上扩散因子立刻越界转速掉到 10500 r/min 以下时第一级载荷系数会冲上 0.7。这些边界才是交给结构和总体去谈的依据。如果想把叶片数、反力度、弦长一起调靠嵌套循环效率太低可以交给 matlab 优化工具箱里的fmincon把整机效率作为目标函数扩散因子、反力度下限、叶片数取整作为约束。写目标函数时注意把所有物理量归一到同一量级否则求解器在叶片数这种量级为 10 的变量和反力度这种量级为 1 的变量之间会走不动步这也是很多人用 matlab 优化工具箱做参数寻优时最先踩的坑。画图这一步别省。用imagesc把effFeasible画出来加上colorbar用plot把每一级的反力度沿叶高叠在一张图上看叶根那条线有没有贴到 0.4 以下。图形化之后的判断速度比看数字表快一个数量级特别是要同时比较三四种级数方案的时候。用 matlab 画图把设计点和两个典型非设计点叠在一起喘振边界和堵塞边界的相对位置一眼就能看出来再回头调第一级的稠度和叶尖间隙方向就明确了。本文还有配套的精品资源点击获取
