液体火箭发动机试验推力测量系统主要由试车架、力源加载装置和测量系统组成。试车架主要由定架、动架、预紧力装置及弹性约束件等组成。典型的大型液体火箭发动机垂直试车架结构形式如图1所示。其中,定架安装在试车台土建基础上,动架和发动机连接,测量传感器安装在动架和定架之间。动、定架力分界面有弹性约束件和水平波纹管。弹性约束件垂直于推力轴线,并能承受一定的横向力。动架、弹性约束件、波纹管组成的动架系统一般要求99%以上的被测推力传递到测量传感器上^[1-2]。预紧力装置的作用是拉起动架及发动机,保证测量传感器预先受到一定预紧力。试验前通过力源加载装置和基准力传感器,对推力测量传感器进行原位校准^[7]。

图1 大型垂直式发动机试车架
Fig.1 Large vertical thrust frame of engine test

图2 推力参数测量原理框图
Fig.2 Principle diagram of the thrust measurement system

压力测量系统由压力传感器、传输电缆、信号调理器、数据采集装置及数据处理分析装置等组成。压力传感器通过测压导管或直接安装在被测位置,一般力源校准装置将基准压力传感器和压力加载装置集成在一起^[8],在试验现场用力源校准装置对测量压力传感器进行现场校准后安装于被测位置。发动机自带的压力传感器一般用传感器的计量校准数据,在现场传感器连接处或数据采集装置输入端用标准电压源加载,进行等效电校准。

转速测量系统主要由转速传感器、传输电缆、信号调理器及数据采集处理装置等组成,如图3所示。转速传感器输出的是频率信号,信号波形为不规则的近似正弦波。转速参数测量的关键是信号调理及抗干扰技术^[9-14],信号调理的主要作用是:滤波、放大、整形、隔离、信号匹配、硬件分频及F-V变换等。

图3 转速信号调理器原理框图
Fig.3 Principle diagram of rotating speed signal converter

目前大型液体火箭发动机试验起动加速性用转速或推力测量值进行计算,轨姿控发动机用燃烧室压力或推力计算。控制和测量系统设置统一的起动计时基准,所有测量参数基准0 s和关机时刻都是统一的。从发动机工作机理分析,只要试验系统满足发动机起动的时序、入口压力、温度等条件,起动加速性由发动机自身性能所决定,与试验系统无关。但由于试验系统的结构和计算参数不同,虽然发动机起动过程推力、燃烧室压力、转速密切相关^[12],但计算的实际结果不尽相同。同一台发动机用不同参数计算结果有一定差异,同一台发动机在不同的试验系统进行试验,加速性计算值也有一定的差异。大型发动机的差异比轨姿控发动机的差异更大一些。差异较大时对研究发动机的起动特性产生一定影响。表1统计了某型号发动机在同一试验台试验时,推力、喷前压力、涡轮转速3个参数在发动机起动时测量值开始变化时刻,用推力计算的起动加速性和压力、转速分别计算的加速性数值之差。

表1 不同参数变化时刻与加速性差异
Tab.1 Starting acceleration of different parameter changing time 单位:s

由表1数据可以看出,同型号发动机在同一试验台试验计算的发动机起动加速性数值,转速最快,推力最慢。转速平均比推力计算的加速性快0.25 s,比喷前压力快0.12 s。喷前压力平均比推力快0.13 s。表1用不同参数计算发动机起动加速性差异的原因主要有两个方面:一是大型发动机试验系统庞大的推力测量机架结构及推进剂供应管路连接形式等综合因素影响推力参数响应慢; 压力参数安装的测压导管对动态响应有一定影响(传感器直接安装时,无此因素)^[8]。第二是测量系统对不同测量参数配置的信号调理器滤波特性的影响。转速采用绝对法直接测量,起动瞬间只需克服涡轮转子的惯性,信号仅受到调理器的滤波特性和门槛电平的影响,起动过程比推力和压力参数响应快。发动机起动的0 s时刻统一,在采样速率相同的条件下,图4所示发动机起动时推力比转速数值发生变化约晚0.32 s。喷前压力参数和推力参数在采样速率和信号调理器的滤波特性一致的条件下,图5所示发动机起动时推力比喷前压力晚0.2 s。

图4 推力和转速参数起动过程曲线
Fig.4 Starting curve of thrust and rotating speed

图5 推力和喷前压力参数起动过程曲线
Fig.5 Starting curve of thrust and pressure in front of the nozzle

f=1/(2π)(k/m)^1/2(1)

式中:m为动架总等效质量,kg; k为动架总等效刚度, N/m; f为动架固有频率,Hz。

图1所示的某发动机试验系统动架总质量约m=50 800 kg,f=3.7 Hz,可以推算出动架的总等效刚度k=2.742*107 N/m。如果动架总等效质量m减小到25 400 kg和12 700 kg,则动架的系统振动频率f提高到5.2 Hz和7.4 Hz,说明减小动架结构质量,提高了动态响应性能。图6用试验获得的实际推力进行激励,仿真不同动架等效质量的系统响应延迟情况。仿真结果证明推力架越重,固有频率越低,响应越慢,振幅越小^[10]。

图6 等效质量与推力系统响应的关系
Fig.6 Response relation between the equivalent mass and thrust

上述讨论的发动机加速性从快到慢依次为转速、喷前压力、推力。但有些小型发动机试验起动时,推力开始反应比较迟缓,一定时间后发生突变,如图7所示。加速性t₁₀时相差0.51 s,到t₉₀时相差仅有0.11 s。从t₁₀到t₉₀逐步变化情况如表2所示。有些小型发动机试验系统采用水平式试车架质量较小,燃烧室压力安装了不充填隔离液的较长测压导管,起动时推力及t₉₀均比室压反应快,如图8所示。说明研究发动机起动加速性时,应针对特定的试验系统和同型号发动机进行分析研究。

图7 推力和转速参数起动过程曲线
Fig.7 Starting curve of thrust and rotating speed

表2 起动过程不同参数加速性对比
Tab.2 Starting acceleration of different parameters in the engine start

图8 推力和燃烧室压力参数起动过程曲线
Fig.8 Starting curve of thrust and combustion chamber pressure

影响测量参数动态响应的另一个因素是测量系统电性能。电性能主要包括传感器频响、信号调理器的滤波器(含采集装置)频响两方面。采集装置和信号调理器分离时,应单独分析各自的频响特性。上述研究的大型试验系统用标准电压源对太平洋6000采集装置加载阶跃电压信号(将前置的调理器滤波频率设置为直通),响应时间t₉₀=0.3～0.4 ms。推力传感器的型号是GKCY15-1A,量程500 kN,频响500 Hz,用砝码瞬间加载,t₉₀=2.0～2.3 ms。压力传感器型号是YB1D,量程15 MPa,频响1.5 kHz,用落锤法加载,t₉₀=0.7～0.9 ms。由上述分析看,采集装置和传感器的动态响应在分析大型液体火箭发动机起动加速性时,影响因素可以忽略,主要分析信号调理器滤波特性的影响。

转速(频率型参数)和推力、压力(电压型参数)配置的信号调理器不同,滤波特性的差异对被测参数的动态响应产生了一定影响。以配置信号调理器的美国太平洋公司6000采集系统为例,说明信号调理器的滤波特性对参数动态响应影响程度。6000采集系统为电压型参数配置了6035采集板卡^[13]。传感器输出电信号经过两个VCVS二阶有源滤低通波器后,再进入采集装置。有源滤低通波器前并联了直通开关,可以将滤波器旁路,满足高频响应的特殊要求。6035采集板前端配置了巴特沃斯(Butterworth)四阶低通滤波电路(两个二阶的Sallen-Key低通滤波电路串联)。6000采集系统为频率型参数配置了6048采集板^[13]。该板卡每个通道预留一个一阶滤波器用于降噪,但出厂时未连接该滤波器。如果直接用6048测量转速参数,则开机和关机延时在转速参数上可以忽略。由于6048采集卡要求的输入信号是100 mV～100 V的过零信号,液体火箭发动机试验中部分转速传感器和涡轮流量计输出信号的幅值小于100 mV,有的转速信号幅值仅4 mV。一般在信号进入6048前研制了多功能信号调理器,原理如图3所示^[9],调理器中对测量有延时作用的主要是前端低通RC滤波器。通过仿真^[10],太平洋公司6000采集系统配置的6035采集板设置10 Hz滤波器与直通相差约70 ms,如图9所示。试验现场用阶跃信号对6035板滤波器进行测试,将采样速率提高到2 000点/秒,6035采集板仍然10 Hz滤波,现场加载结果如图 10 所示,带滤波器的信号比直通的信号晚72 ms,验证了仿真结果的正确性。由于6048板前设计了多功能信号调理器,经现场试验验证6035采集通道实际比6048采集通道延时40～50 ms,统计某型号发动机喷前压力起动峰值比转速起动峰值晚40～58 ms,平均值约48 ms。试验结果证明,测量参数时间基准统一的前提下,不同测量参数信号调理器(滤波器)引起的测量延时应在数据处理环节进行修正,确保不同参数在同一时刻的测量值反映发动机的固有特性。如太平洋6000采集系统采集的推力和压力参数数据处理时,时间坐标应相对于转速参数提前40～50 ms左右。实际工作中,标准低通滤波的延迟时间^[15]

T_τ=1/8·n/(f_c)(2)

式中:T_τ为延时时间,s; f_c为截止频率,Hz; n为滤波器阶数。

图9 阶跃信号输入仿真
Fig.9 Simulation of step signal input

图 10 6035滤波器延时测试
Fig.10 The delayed measurement of the 6035 filter