超燃冲压发动机燃烧室主动冷却设计研究

超燃冲压发动机燃烧室主动冷却模型制造和热考核试验的费用高昂、周期长,为了降低试验风险,采用三维计算和经济的验证试验相结合的方法开展了超燃冲压发动机燃烧室主动冷却设计。冷却结构基本参数设计和侧壁冷却流动设计是确定设计方案的两大基础,前者采用三维传热计算结合冷却面板传热验证试验完成,后者采用超临界燃料流动三维并行计算结合水流动验证试验完成。在此基础上,经过多轮的结构设计与三维传热及强度计算评估迭代,确定了最终的燃烧室主动冷却结构。设计的主动冷却燃烧室在来流马赫数2.5,总温1700K条件下成功通过200s热考核试验,表明所采用的设计方法、验证试验和计算工具是有效和可信的。     

基于MATLAB/Simulink的飞机燃油箱内燃油温度仿真计算

燃油温度是飞机燃油箱可燃性评估的关键参数之一。为了研究飞机燃油箱内温度的分布规律及传热模型的正确性,采用MATLAB/Simulink软件平台搭建出了客机的燃油箱热仿真模型。在输入飞行实验所对应的的边界条件后,通过数值模拟获得燃油箱内部各计算节点处的燃油温度。结果表明,在三种不同的航行条件下,燃油箱热模型仿真结果均能较好地与飞行实验结果吻合,能够将计算误差控制在一定范围内。采用该燃油箱热模型进行热仿真数值模拟可以获得较为准确的飞机燃油箱热特性,能够在飞机设计阶段,用于对燃油箱结构以及内热源部件布置的优化。     

Gray bootstrap method for estimating frequency-varying random vibration signals with small samples

During environment testing, the estimation of random vibration signals （RVS） is an important technique for the airborne platform safety and reliability. However, the available meth- ods including extreme value envelope method （EVEM）, statistical tolerances method （STM） and improved statistical tolerance method （ISTM） require large samples and typical probability distri- bution. Moreover, the frequency-varying characteristic of RVS is usually not taken into account. Gray bootstrap method （GBM） is proposed to solve the problem of estimating frequency-varying RVS with small samples. Firstly, the estimated indexes are obtained including the estimated inter- val, the estimated uncertainty, the estimated value, the estimated error and estimated reliability. In addition, GBM is applied to estimating the single flight testing of certain aircraft. At last, in order to evaluate the estimated performance, GBM is compared with bootstrap method （BM） and gray method （GM） in testing analysis. The result shows that GBM has superiority for estimating dynamic signals with small samples and estimated reliability is proved to be 100% at the given confidence level.     

Accuracy enhancement of the spherical actuator with a two-level geometric calibration method

This paper presents a two-level geometric calibration method for the permanent magnet （PM） spherical actuator to improve its motion control accuracy. The proposed actuator is com- posed of a stator with circumferential coils and a rotor with multiple PM poles. Due to the assembly and fabrication errors, the real geometric parameters of the actuator will deviate from their design values. Hence, the identification of such errors is critical for the motion control tasks. A two-level geometric calibration approach is proposed to identify such errors. In the first level, the calibration model is formulated based on the differential form of the kinematic equation, which is to identify the geometric errors in the spherical joint. In the second level, the calibration model is formulated based on the differential form of torque formula, which is to calibrate the geometric parameters of the magnetization axes of PM poles and coils axes. To demonstrate the robustness and availability of the calibration algorithm, simulations are conducted. The results have shown that the proposed two-level calibration method can effectively compensate the geometric parameter errors and improve the positioning accuracy of the spherical actuator.     

Flutter analysis of an airfoil with nonlinear damping using equivalent linearization

The equivalent linearization method （ELM） is modified to investigate the nonlinear flut- ter system of an airfoil with a cubic damping. After obtaining the linearization quantity of the cubic nonlinearity by the ELM, an equivalent system can be deduced and then investigated by linear flut- ter analysis methods. Different from the routine procedures of the ELM, the frequency rather than the amplitude of limit cycle oscillation （LCO） is chosen as an active increment to produce bifurca- tion charts. Numerical examples show that this modification makes the ELM much more efficient. Meanwhile, the LCOs obtained by the ELM are in good agreement with numerical solutions. The nonlinear damping can delay the occurrence of secondary bifurcation. On the other hand, it has marginal influence on bifurcation characteristics or LCOs.     

Aeroelastic two-level optimization for preliminary design of wing structures considering robust constraints

An aeroelastic two-level optimization methodology for preliminary design of wing struc- tures is presented, in which the parameters for structural layout and sizes are taken as design vari- ables in the first-level optimization, and robust constraints in conjunction with conventional aeroelastic constraints are considered in the second-level optimization. A low-order panel method is used for aerodynamic analysis in the first-level optimization, and a high-order panel method is employed in the second-level optimization. It is concluded that the design of the abovementioned structural parameters of a wing can be improved using the present method with high efficiency. An improvement is seen in aeroelastic performance of the wing obtained with the present method when compared to the initial wing. Since these optimized structures are obtained after consideration of aerodynamic and structural uncertainties, they are well suited to encounter these uncertainties when they occur in reality.     

大型回流边界层风洞的气动与结构设计

大型边界层风洞是开展风工程研究的必备装备。以浙江大学ZD-1边界层风洞的研制为背景,详细介绍了大型回流边界层风洞气动设计和立式结构设计中的关键问题,在风洞气动设计时采用了收缩比为4∶1的单回路单试验段气动轮廓,在试验段中设置了0．22°的当量扩散角,对拐角导流片外形作了特殊处理,并采用钢结构与混凝土结构相结合的立式结构。流场校验结果表明,大型回流边界层风洞的气动与结构设计能满足设计要求,某些指标甚至达到航空风洞的标准,在试验段中设置扩散角有利于降低轴向静压梯度,立式结构设计对提高试验段气流的水平均匀性有一定的作用,可为今后类似风洞的研制提供参考。     

Transient simulation of regression rate on thrust regulation process in hybrid rocket motor

The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB(hydroxyl-terminated polybutadiene) hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.     

喉部长度和扩张角对离心压气机管式扩压器影响研究

针对带楔形扩压器的高负荷离心压气机,利用数值模拟手段对其进行了研究,计算获得的离心压气机特性与试验结果符合较好。在此基础上,设计了具有不同喉部长度和扩张角的管式扩压器,利用经过校核的数值模拟手段对管式扩压器的喉部长度和扩张角进行了参数化研究。研究发现较长的喉部会使压气机设计点的性能下降,但是却可以抑制近失速点扩压器内的流动分离,喉部长度与喉部直径之比为0.5时综合性能较优;较大的扩张角会造成扩压器内流动分离,增加流动损失,扩张角由4°增加到12°时压气机效率降低了约3%。     

LHT-100霍尔推力器热特性模拟分析

为了对LHT-100霍尔推力器提出热设计优化措施,采用有限元仿真软件进行LHT-100霍尔推力器的稳态、瞬态及空间在轨环境模拟热分析研究,并通过热平衡试验进行了结果比对。分析及试验结果表明,处于工作状态时霍尔推力器的高温部件主要是放电腔、阳极和导磁底座,而受高温影响薄弱部件内线圈、气路组件的温度则分别达到了约401~421℃和141~381℃。热设计优化建议为,在放电腔与内线圈之间增加独立热屏结构后可以有效降低内线圈温度约80~90℃,在阳极气路组件上存在的热应力会是影响霍尔推力器可靠性的重要因素,需要在热设计中得到充分考虑。