1、空间涡轮发电系统设计及性能特性仿真分析韩 灿，李育隆，田 林，于新刚，王鹏程（北京航空航天大学能源与动力工程学院，北京；北京空间飞行器总体设计部，北京）摘要：为解决空间飞行器应急发电问题，利用空间飞行器携带的推进剂为能量源，设计了一套空间涡轮发电系统，建立了稳态及启动过程的数学模型，对系统的稳态和启动特性进行了仿真分析。结果表明：空间涡轮系统可以利用推进剂实现高功率发电，系统可以启动和稳定运行；系统在设计点稳定运行时，推进剂流量为.，燃气发生器室压为.，燃气发生器出口温度为，输出电功率达到；当系统运行参数发生扰动偏离设计点时，系统性能会相应变化，混合比不变且氧化剂贮箱增压压力大于设计点.时，涡
2、轮输出功率变为设计值的.倍；系统启动过程仿真得到的性能参数和稳态设计值相对误差小于，验证了系统设计的合理性。关键词：空间飞行；涡轮发电；系统设计；性能仿真；能量利用中图分类号：文献标识码：文章编号：（）收稿日期：；修回日期：第一作者：韩灿，女，硕士研究生，研究方向为航天器能源利用技术。：，（，；，）：，.，.，：；第 卷 第 期 年 月 载 人 航 天 DOI:10.16329/ki.zrht.2023.01.015 引言 空间飞行器的生命保障、轨道与姿态控制等分系统均需持续可靠的电能供给才能稳定工作。但当空间飞行器发生故障导致空间电源系统不能正常工作时，如何紧急提供电能维持各分系统的运行将是
3、空间飞行器空间电源系统设计中的一项重要任务。如阿波罗 载人飞船地月出发后 天，服务舱的氧气罐爆炸导致燃料电池无法工作，造成电力系统严重降级。因此，有必要寻求空间飞行器特别是载人飞行器应急供电措施。空间飞行器携带的液体推进剂在燃气发生器中将会燃烧产生高温高压燃气，因此当电源系统发生故障时，利用这种高温高压燃气作为应急电源系统的能量来源将是一种可行的应急方案。例如阿波罗 在发生服务舱爆炸时，飞行器内还有大量未使用的推进剂，这些推进剂最终都被浪费。因此，可以使用剩余液体推进剂作为电能来源进行有效利用。涡轮是一种能量转化装置，在航空航天、发电等领域均有广泛应用。应急供电系统可以以涡轮作为转化高温高压气
4、体能量的主要部件，将液体推进剂通入燃气发生器中燃烧产生高温高压燃气，燃气再驱动微型涡轮机做功发电。热力发电系统的稳态及启动特性直接决定着系统性能的优劣。白杰对水下航行器的热电联合闭式循环动力系统进行了稳态参数设计，给出了设计工况稳态性能参数，建立了系统动态模型，通过系统级的计算，得到了系统动态特性，是系统性能预示的手段；利用模块化建模的思想建立了微型燃气轮机的动、静态模型，是目前公认的比较经典的微燃机模型，这种模块化建模的方法直观得到了微燃机标准工况下的特性。因此，针对于热力发电系统有必要进行稳态及启动特性的研究。综上，为满足空间飞行器能源系统的应急供电需求，同时解决液体推进剂有效利用问题，本
5、文提出空间涡轮发电系统并进行详细设计，建立稳态及启动过程数学及仿真模型，分析系统工作过程中的性能特性及参数变化情况。系统设计和数学模型的建立.系统设计 设计得到的空间涡轮发电系统构型如图 所示。系统包括的主要部件为：高压气瓶、减压阀、燃料贮箱、开关阀、燃气发生器、涡轮、对称喷管、行星齿轮箱、发电机、整流器等。图 空间涡轮发电系统构型图 .推进剂供应部分 系统采用的推进剂供应部分结构来自于空间飞行器已有的挤压式推进系统，主要功能是提高推进剂的压力以及以设计的质量流率向燃气发生器供给推进剂。由高压气体气瓶、减压阀、燃料贮箱、阀门和管路组成。燃气发生器的室压设计值为.，为了按设计值供给推进剂，氧化剂
6、贮箱增压压力为.，燃料贮箱增压压力为.。系统采用再生冷却的设计使燃料在进入燃气发生器燃烧之前先在冷却通道内循环流动，从而冷却高温部件并且增大燃料进入燃烧室前的初始能量，提高系统能量转换效率，系统在设计状态下的系统热效率为。.燃气发生器涡轮部分 燃气发生器为涡轮装置提供做功工质，燃气温度越高，涡轮做功能力越强，同时为了适应涡轮叶片材料的许用温度，选择设计状态下燃气绝热燃烧温度为 ，此时推进剂氧燃质量混合比为.，燃气发生器室压为.。图 为当推进剂总质量为 ，其中四氧化二氮和一甲基肼的质量比为设计混合比时，发电功率和发电时间的关系图，选择设计发电功率为 ，持续供电时间为。当整流器最终输出发电功率为
7、时，考虑到涡轮输出功率转化为输出电功率过程中耗散的功率，设计点时涡轮输出功率为 。发电涡轮只需输出功率带动外部负载即可，不需要高通流的设计去提高推力。设计状态燃气流量为.，流量越小需要涡轮的膨胀比越大，也即涡轮的负载越高，涡轮设计状态下膨胀比第 期 韩 灿，等 空间涡轮发电系统设计及性能特性仿真分析图 系统稳定工作时间与发电功率的关系 为。高转速是高负荷微型涡轮设计的关键，系统涡轮设计转速选择为 ，低于用于水下无人航行器燃气涡轮和微机电系统的微型涡轮转速，设计涡轮转速合理可行。系统采用的涡轮具有工质流量小、高负荷的特点。废气可以通过涡轮后的对称管路排出，不对空间飞行器产生额外的推力和力矩影响。
8、涡轮进口和燃气发生器出口可以通过法兰、石墨密封圈进行连接和密封，连接处需要解决密封、振动和耐热问题。气浮轴承具有高速、高寿命、耐热、耐冷等优点，本系统的涡轮转速较高，选择气浮轴承用于支撑转子，具有润滑、减振、冷却等作用。高压气瓶的高压气体可作为气浮轴承中的气体工质。涡轮发电机转子系统的转速较高，为了防冷焊，需要在金属部件表面镀高性能固体润滑薄膜。.发电部分 发电机使用微型永磁同步发电机，选择行星齿轮箱连接涡轮和电机。整流器可以将发电机产生的交流电转化为用电设备所需的直流电输出，设计状态下输出直流电压为 。发电部分与空间飞行器中的电能储存、变换、调节和分配装置相连，在选用连接接口器件时，需要考虑
9、工作电压等参数与抗振动、抗冲击、防电磁干扰性能及良好的高真空性能。综上，该系统设计工况下主要参数如表 所示。常规空间飞行器的推进系统中有推进剂供应系统和燃气发生器等装置，因此本文提出的空间涡轮发电系统可以在其基础上进行集成。以阿波罗 号为例，只需在飞行器上添加涡轮、行星齿 表 涡轮发电系统主要设计参数 参数数值氧化剂贮箱增压压力.燃料贮箱增压压力.燃气发生器室压.燃气发生器出口温度 燃气流量（）.推进剂混合比.涡轮转速（）涡轮输出功率 涡轮膨胀比涡轮出口温度.输出直流电压 输出电功率 系统热效率 系统稳定工作时间 推进剂使用总量 轮箱、发电机、整流器及相应的连接管路即可组成应急供电系统，额外补
10、充部件的总重量大约为，质量比功率为 。对于空间飞行器来说，目前常用的空间电源有燃料电池、太阳能电池以及空间核反应堆电源。燃料电池的质量比能量为 ，但燃料电池成本高，寿命短；比功率最大的太阳能电池可以达到 ，但发电受光照强度的影响，同时太阳能电池阵尺寸非常大；核反应堆电源为防核辐射污染需加厚重屏蔽防护装置，在重量比功率上并不具有优势。本文设计的空间涡轮发电系统相对于其他空间电源来说，重量和尺寸较小，比功率较大，在阿波罗 发生故障时，大量推进剂被抛弃，因此使用本系统进行应急供电还可以提高推进剂的利用效率。.稳态工作建模 建立稳态模型的基本思路为：首先建立系统相关部件的数学模型，再通过系统的平衡关系
11、建立非线性方程组求解，主要部件的稳态数学模型如下。.推进剂供应部分 推进剂供应部分简要介绍液路部分的静态数学模型，其余部件的数学模型可通过文献获得。载 人 航 天第 卷使用集中参数方法描述液路部分的静态数学方程如式（）所示。（）其中，为管道损失压力，为损失系数，为管路液体流量，为液体密度。.燃气发生器涡轮部分 燃气发生器中的热力计算如式（）所示。（，）（，）（）该式由文献数据拟合公式得到。为燃气热值，为燃气绝热指数，为燃气定压比热容，为推进剂混合比，为燃气发生器室压，为液体推进剂进入燃气发生器的初温。对于超声速涡轮级，需要建立其流量、效率和功率等数学模型，可由文献得到。.发电部分 永磁发电机出
12、口线电压 为式（）：（）（）（）其中，为发电机出口线电压的幅值，为电压系数，为发电机转动角速度。对于三相全波整流桥，考虑换相重叠角时，输出的直流电压 可以表示为式（）：（）其中，为发电机定子绕组的漏感，为整流器直流侧的直流电流。.系统参数平衡模型 对于主流路，流量平衡方程如式（）、（）所示。（）（）其中，为液体推进剂总质量流量，为氧化剂质量流量，为燃烧剂质量流量，为燃气质量流量。对于冷却回路，流量平衡方程如式（）所示。（）其中，为冷却液的质量流量。贮箱增压压力由减压阀的出口压力决定，认为减压阀为定值减压阀，出口压力保持稳定，在给定贮箱增压压力的情况下计算系统及液路各组件的静态特性。氧化剂流路和
13、燃料流路的压力平衡方程为式（）和式（）：（）（）（）上式中，为贮箱增压压力，为贮箱增压压力与贮箱出口压力之差，为推进剂通过开关阀的压降，为推进剂通过液体管道的压力损失，为推进剂通过喷注器的压降，为燃气发生器室压，表示冷却管路压降拟合系数，下标后加 表示燃料，下标后加 表示氧化剂。在系统工作时，涡轮发电部分的功率和转速需要达到平衡，平衡方程如式（）（）所示。（）（）其中，为涡轮输出的轴功率，为整流器输出的电磁功率，为涡轮功率转化为输出电功率过程中耗散的功率，为涡轮的转速，为行星齿轮箱传动比，为永磁发电机转速。.启动过程建模 建立系统启动过程模型的基本思路是：首先建立部件的动力学数学模型，再利用各
14、部件之间的压力、流量、温度等信号传递关系将各部件的动力学数学模型组合构成整个系统的数学模型。.推进剂供应部分 推进剂供应部分只介绍液路的动态数学模型，其 他 部 件 的 动 态 数 学 模 型 可 在 文 献 中查阅。利用集总参数模型建立管路流动模型，同时考虑液路流体的惯性、流阻和压缩性后，液路方程如式（）所示。|()()|（）其中，为管路的截面积，为管路的长度，为流体声速，为体积弹性模量，为管路的体积，为损失系数。.燃气发生器部分 为了对燃气发生器建立一个简单的动力学模型，做出如下假设：忽略燃烧时滞；任何瞬间，燃烧室内的压力、混合比、温度均匀分布；燃烧产物为理想气体；燃气发生器绝热。燃气发生
15、器压力的方程如式（）所示。第 期 韩 灿，等 空间涡轮发电系统设计及性能特性仿真分析()（）()（）其中，为燃气发生器体积，为给燃气发生器供应的氧化剂流量，为给燃气发生器供应的燃料流量，为燃气发生器出口燃气流量。燃烧室内混合比的方程如式（）所示。（）()（）.涡轮发电机转子部分 从系统角度仿真计算时，并不需要准确反映涡轮机、发电机和整流器内部工作过程，只需要求解反映主要性能的参数。因此，用准稳态模型对其进行计算。涡轮提供动力，带动发电机转动，涡轮发电机转子系统动力学模型如式（）所示。（）其中，为涡轮转矩，为发电机负载转矩，为摩擦转矩，为转动惯量，为涡轮转动角速度。以燃气发生器模块为例，使用 建
16、立该模块仿真模型，如图 所示。该模块可以计算输出燃气发生器的室压、燃气混合比、燃气热值、燃气比热比随时间的变化。建立完各部件模块的仿真模块后，连接输入输出端口，组成系统的仿真模型。图 燃气发生器的 仿真模块 性能特性仿真与分析.稳态工作性能分析 空间涡轮发电系统复杂，参数耦合性强，运行过程中系统性能参数可能偏离额定设计状态，并且重新平衡在一个新的稳定状态。需研究系统在各个稳定状态下，各部件的工作参数及系统性能随影响因素变化的关系。.固定混合比，改变贮箱压力 图 是保持氧化剂和燃料混合比 不变时，燃料贮箱增压压力、燃气发生器室压、涡轮输出功率 和整流器输出直流电压 标准化后随氧化剂贮箱增压压力 变化的关系图。图中纵坐标各参数的值为标准化后的相对值，用“”标出了设计点位置。仿真结果显示，随着氧化剂贮箱增压压力的增大，燃料贮箱增压压力、燃气发生器室压、涡轮输出功率和整流器输出直流电压均会升高，涡轮输出功率的相对值变化最大；当氧化剂贮箱增压压力大于设计点.时，涡轮输出功率大约为设计值的.倍；在氧化剂贮箱增压压力的变化范围内，种参数大致呈线性变化趋势。图、随氧化剂贮箱压力变化、图 是燃气流量、涡
