三栅极离子推力器电子反流失效影响参数的敏感性研究

电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC (Particle in Cell-Monte Carlo Collision)方法,构建了三栅极系统数值仿真模型。采用模型研究了地面真空舱本底压力、屏栅电压、加速栅电压、屏栅与加速栅间距、屏栅上游等离子体密度和放电室工质利用率等参数的影响敏感度对比。研究结果显示,真空舱本底压力可以作为加速寿命试验的首选加速应力,在推力器结构和工作本征参数中工质利用率为最敏感应力,其次是屏栅电压、屏栅上游等离子体密度、加速栅电压、屏栅和加速栅间距。     

30cm离子推力器三栅极组件设计参数对寿命的影响研究

为了研究30cm离子推力器三栅极组件设计参数对预估寿命的影响,在完成失效模式分析的基础上,通过PIC-MCC方法对离子推力器三栅极组件的离子溅射速率进行了计算,建立起栅孔二维寿命预估模型,并针对栅极设计参数对预估寿命的影响进行研究。结果显示:导致三栅极组件的主要失效模式为5kW高功率模式下的离子直接轰击所造成的栅极早期结构失效,且减速栅的过快离子溅射腐蚀成为影响三栅极组件寿命的关键,而不同工作模式不会产生新的失效方式,仅影响栅极的离子溅射速率以及寿命;在现有三栅极设计参数条件下,当推力器工作时,栅极引出的离子束流处于明显欠聚焦状态,且加速栅寿命预估值约为9062h,而减速栅约为2642h;通过PIC-MCC方法得到的栅极三个关键设计参数对寿命的影响模拟结果显示,降低加速栅电压对提升减速栅寿命的作用较小;缩小加速栅与减速栅冷态间距后,离子溅射速率会随着冷态间距的缩小逐渐降低,冷态间距由1mm缩小至0.6mm后,减速栅在5kW工况下的工作寿命可提升至10726h,且经试验验证该间距可满足推力器力学环境试验要求;缩小屏栅孔径对改变离子束流引出形状具有显著作用,单孔束流发散角度随着屏栅孔径的缩小出现了明显降低,且束流离子几乎不会再直接轰击至减速栅上游区域,当屏栅孔径由1.9mm缩小至1.6mm后,减速栅工作寿命可提升至9259h;分析结果对后续开展栅极组件的寿命优化设计提供了参考。     

离子推力器栅极系统电子反流阈值的数值分析

阻止束流等离子体中电子反流到加速栅上游区域是离子推力器加速栅负电压的主要作用之一,能够阻止电子反流的加速栅电压最小值称为电子反流阈值。加速栅电压的选择直接影响到离子推力器的工作性能和运行寿命,电子反流阈值电压是确定加速栅电压的重要参考参数。基于PIC方法计算了20cm氙离子推力器加速栅电子反流阈值,并分析了加速栅孔径、栅间距、单孔引出束流电流大小对加速栅电子反流阈值电压的影响,计算结果与试验测量值符合较好。该数值模型为加速栅参数的选择和降低电子反流失效风险方法提供了参考,为下一步电子反流现象对加速栅寿命的预测分析奠定了基础。     

不同环境温度对30cm离子推力器三栅极组件的离子刻蚀速率影响分析

为了对不同环境温度造成的30cm离子推力器三栅极组件离子刻蚀速率的影响进行分析,采用有限元仿真与试验验证相结合的方法,计算并试验验证了不同环境温度下的三栅极组件热平衡温度以及栅极间的相对位移变化,采用流体方法模拟了不同环境温度(20℃,-70℃,-120℃和-170℃)对三栅极组件的刻蚀影响,并结合短期寿命试验结果进行验证。结果显示:随着环境温度的降低,屏栅达到温度平衡的时间无变化,而加速栅温度平衡所需的时间则明显延长,20℃下的屏栅和加速栅热仿真结果与室温下推力器热平衡试验结果比对误差分别为7%和5%;其次环境温度的降低,会导致屏栅与加速栅的中心间距和边缘间距均缩小,而加速栅和减速栅的边缘间距却逐渐拉大,仿真结果与栅极热间距摄像测量结果符合性较好;根据三栅极组件的栅孔径扩大率随环境温度变化的计算结果来看,加速栅中心和减速栅边缘是离子刻蚀的主要位置,轰击至加速栅中心区域的离子数速率约是边缘的3倍,而轰击至减速栅边缘区域的离子数速率是中心的2.5倍,且环境温度的降低对加速栅中心区域离子刻蚀的影响更为强烈;经2100h的寿命试验验证,仿真结果与试验结果基本符合,误差经分析认为主要来自于流体方法的参数设置过程以及栅孔壁面均匀刻蚀的计算假设。     

30cm离子推力器在轨环境下的工作点及热控措施研究

栅极间距变化是影响离子推力器在轨环境下从冷态条件正常点火启动的重要因素,同时也决定了离子推力器的在轨工作时机和热控实施策略。本文采用有限元仿真与地面热平衡试验验证相结合的方法,建立起30cm离子推力器有限元分析模型并进行了模型校验,之后对离子推力器在轨受太阳光照影响的栅极温度场分布和间距变化,以及推力器在5kW工况下的三个典型温度点所对应的栅极间距变化进行了仿真分析,最后考虑了主动热控干预对推力器最恶劣工作点的栅极间距变化影响。结果显示：纯太阳光照影响下的栅极组件存在周期性温度变化,栅极最大温差可达到100℃,栅间距缩小量在0.06mm~0.16mm范围内波动;在太阳光照基础上实施60W的主动热控后,栅极最大温差降低至60℃,栅间距缩小量波动范围则变为0~0.03mm;栅极最高温度点和最低温度点分别是推力器冷态启动最容易和最困难的两个工作时机点,两点所对应的启动后屏栅和加速栅最小间距分别为0.22mm和0.04mm;在10W、70W和120W的热控加热功率下,从最低温度点启动后的屏栅和加速栅最小间距分别为0.06mm、0.20mm和0.29mm;采取主动热控措施能够有效降低推力器工作过程中的栅极热形变位移峰值,且加热功率为120W即温控点温度为50℃的主动热控可以满足30cm离子推力器在轨冷态启动时的0.25mm安全栅极间距要求。     

射频自偏压效应离子推力器束流特性实验研究

为了获得不同推力器工况对射频自偏压效应离子推力器的自偏压幅值和束流特性的影响，本文通过地面实验研究了栅极射频功率、线圈放电功率、工质种类对自偏压幅值和羽流区等离子体参数的影响，同时对直流栅极工况和射频栅极工况下的束流特性进行了对比。研究结果表明：栅极射频功率的增大会提高自偏压幅值并提升束流强度，但在较高栅极射频功率下，栅极下游区域将发生自持放电并形成等离子体；放电腔内放电模式转换会通过改变等离子体阻抗的方式大幅影响栅极直流自偏压幅值和栅极电压的射频分量，进而影响羽流区等离子体参数；与直流栅极工况相比，射频栅极能同时引出并加速离子和电子，并在栅极下游实现自中和，且在Ar,Kr,Xe三种工质下均具有自中和能力。     

配电方式对Kaufman离子推力器性能的影响研究

目前Kaufman离子推力器主要有两种最具代表性的配电方式:屏栅极电源正端分别连接阳极电源正、负端的配电方式。为了研究配电方式对Kaufman离子推力器工作性能的影响,基于等离子体理论和推力器工作原理,分析两种主要配电方式下放电室电极电势及电流平衡关系,推导了放电室等离子体特性表达式,理论分析了配电方式对离子推力器多种性能参数的影响。结合兰州空间技术物理研究所自研的LIPS300离子推力器在两种配电方式下工作在3kW和5kW的性能试验,通过解析方法对离子推力器多种工作参数和性能参数进行分析,试验结果与理论分析结果具有良好的一致性。研究表明:采用屏栅极电源正端连接阳极电源负端的配电方式能够获得更大的推力和比冲,并能提高离子对栅极透明度,减少离子对屏栅极的溅射,从而提高栅极寿命,但束离子产生成本稍高。研究结果可为离子推力器配电方式的设计与优化提供依据。     

离子推力器欠聚焦冲击电流的数值模拟

离子光学系统的离子束引出过程是离子推力器重要的物理过程,该过程直接关系到推力器的推力、比冲、效率等参数。为研究离子在离子推力器光学系统中的运动特性,使用了基于IFE-PIC（Immersed Finite Element Particle-In-Cell）的离子推力器光学系统离子束引出过程的三维数值计算模型,计算了栅极间电场分布、电荷密度,栅极冲击电流及欠聚焦极限。计算结果表明,当屏栅极电压不同时,发生欠聚焦的等离子束电流也不同。在欠聚焦工况下,一部分离子与栅极碰撞,产生冲击电流。冲击电流随电离室等离子体数密度增加而增大。     

射频离子推力器多元工质束流调节试验研究

为了获得射频离子推力器离子束流随放电参数的变化规律,采用试验研究的方法,就推力器引出束流与射频功率强度、工质种类、工质流量之间的调节规律开展了研究,搭建了射频离子推力器束流调节试验系统。研究结果表明:屏栅电压1200V,加速电压-250V,射频功率200W～700W,工质流量0.2mg/s～4.76mg/s,Xe,Ar,O_2,N_2四种工质下能够可靠放电并稳定引出,实现束流从54mA～467mA的调节,电离效率XeArO_2N_2,离子束流随射频功率和工质流量线性增加,在1.01mg/s的氙工质下,推力、比冲随射频功率从100W～400W线性增加实现推力7.35mN～27.5mN,比冲1191s～3696s大范围连续可调,工质利用率为21.1%～78.8%,并在射频功率为276W时工质利用率和功耗之间存在明显拐点,在应用中要根据任务选择最佳工作区间,合理控制工作参数可以提高推力器工作性能和效率。     

离子推力器加速栅寿命概率性分析

交换电荷离子对加速栅极的溅射腐蚀是离子推力器的关键失效模式之一,基于交换电荷离子对加速栅溅射腐蚀的物理机理,对离子推力器加速栅工作寿命进行了概率性建模。利用该模型对20cm Xe离子推力器加速栅寿命和其达到预期寿命的可靠度进行了评估。结果显示加速栅的寿命近似服从高斯分布,当推力器工作环境压力近似6.7×10-3Pa时,加速栅工作寿命达到3kh的可靠度为0.9352。     

ν-gap度量及其在飞行控制律评估中的应用

 传统控制律评估方法主要用于单输入单输出（SISO）系统，且对模型参数摄动考虑不够全面,针对这些不足，研究了ν-gap度量方法。在介绍系统广义稳定裕度相关概念的基础上，给出了ν-gap度量的定义、特点和性质以及近似摄动模型的计算，提出ν-gap度量评估控制律的步骤。实例结果表明，该方法不仅克服了上述传统评估方法的缺陷，而且还有根据所求的各摄动影响情况忽略影响小的元素，以减少计算量及可以找到最坏情况下的参数摄动组合等优点。     

Integrated Entry Guidance for Reusable Launch Vehicle

A method for the implementation of integrated three-degree-of-freedom constrained entry guidance for reusable launch vehicle is presented. Given any feasible entry conditions, terminal area energy management interface conditions, and the reference trajectory generated onboard then, the method can generate a longitudinal guidance profile rapidly, featuring linear quadratic regular method and a proportional-integral-derivative tracking law with time-varying gains, which satisfies all the entry corridor constraints and meets the requirements with high precision. Afterwards, by utilizing special features of crossrange parameter, establishing bank-reversal corridor, and determining bank-reversals according to dynamically adjusted method, the algorithm enables the lateral entry guidance system to fly a wide range of missions and provides reliable and good performance in the presence of significant aerodynamic modeling uncertainty. Fast trajectory guidance profiles and simulations with a reusable launch vehicle model for various missions and aerodynamic uncertain-ties are presented to demonstrate the capacity and reliability of this method.     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

航空发动机管路测量数据分割方法

 为了提高航空发动机管路测量数据的反求建模效率,提出了一种区域增长分割算法。该方法主要是利用管路表面在曲面索引系数映射、主曲率映射、高斯映射上的特性,并基于均值漂移算法和遗传算法提出了区域分割算法中种子区域的选择策略。然后利用区域增长分割方法实现了对管路测量数据的分割。经过仿真和实测数据验证,所提管路分割算法具有较好的分割质量和效率。     

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

Interfacial Properties Modification of Carbon Fiber/ Polyarylacetylene Composites

This work was dedicated to performing surface oxidation and coating treatments on carbon fibers (CF) and investigating the changes of fiber surface properties after these treatments, including surface composition, relative volume of functional groups, and sur- face topography with X-ray photoelectron spectroscopy (XPS) and atom force microscopy (AFM) technology. The results show that, after oxidation treatments, interfacial properties between CF and non-polar polyarylacetylene (PAA) resin are remarkably modified by removing weak surface layers and increasing fiber surface roughness. Coating treatment by high char phenolic resin solution after oxida- tion makes interface of CF/PAA composites to be upgraded and the interfacial properties further bettered.     

Thermal Stresses and Cracks During the Growth of Large-sized Sapphire with SAPMAC Method

本文利用有限元法对冷心放肩微量提拉法(SAPMAC)生长大尺寸蓝宝石单晶过程中的热应力进行了模拟计算,应用临界缺陷理论解释了裂纹的萌生与扩展机理。研究表明,晶体内的热应力主要与晶体生长速率,环境温度和结晶取向有关;较大热应力多出现在结晶界面,放肩、收尾以及直径急剧变化等位置;最大热应力总是出现在籽晶与新生晶体的界面附近。裂纹将在临界缺陷位置产生,并在应力作用下沿a或m面扩展。计算结果与实验结果基本吻合,通过晶体生长系统和生长工艺的改进,对晶体的开裂问题得到了有效的抑制。     

Optimization Design System for Composite Structures Based on Grid Technology

To solve the topology optimization of complicated multi-objective continuous/discrete design variables in aircmit structure design, a Parallel Pareto Genetic Algorithm （PPGA） is presented based on grid platform in this paper. In the algorithm, the commercial finite element analysis （FEA） software is integrated as the calculating tool for analyzing the objective functions and the filter of Pareto solution set based on weight information is introduced to deal with the relationships among all objectives. Grid technology is utilized in PPGA to realize the distributed computations and the user interface is developed to realize the job submission and job management locally/remotely. Taking the aero-elastic tailoring of a composite wing for optimization as an example, a set of Pareto solutions are obtained for the decision-maker. The numerical results show that the aileron reversal problem can be solved by adding the limited skin weight in this system. The algorithm can be used to solve complicated topology optimization for composite structures in engineering and the computation efficiency can be improved greatly by using the grid platform that aggregates numerous idle resources.