扭曲尾翼弹箭的马格努斯数值研究

为提升掠飞攻顶弹箭较高转速下的飞行稳定性,运用数值计算方法研究了弹体-扭曲尾翼组合体在飞行过程中的马格努斯效应气动机理,并应用标准尾翼弹(BFM)模型的实验数据对数值方法进行了验证。分别研究了带有平板尾翼和扭曲尾翼弹体模型的马格努斯力和力矩随攻角的变化规律,并针对弹体弹翼组合体产生马格努斯效应的机理深入分析。结果表明,扭曲尾翼可有效改善翼面的压力分布,并降低弹体对翼面马格努斯效应的干扰,在大攻角时其表现更胜一筹;弹体所受马格努斯力较大,主要集中在受到涡对称畸变的尾锥部;尾翼主要由于弹体干扰以及几何外形的影响马格努斯力集中在尾部,两者产生的马格努斯力矩数值相差不大,但方向相反。     

腹支撑对飞行器高超声速试验结果影响的数值计算

用数值求解N-S方程模拟翼身组合体腹支撑干扰流场,研究了高超声速条件下腹支撑对飞行器流场、物面压力和整弹气动特性的影响。研究表明,腹支撑支杆对翼身组合体外形气动特性的影响,比我们通常概念中的影响要小得多,原因在于腹支杆不但在支杆附近产生激波,而且支杆的尾迹对后体及后体弹翼附近流场和表面压力带来很大的影响。腹支杆在支杆附近产生的激波对弹体的轴向力和法向力也产生正的影响,而支杆的尾迹对轴向力和法向力产生负的影响,它们之间的相对大小与来流攻角和Mach数有关,因此在一定的攻角条件下,反而引起轴向力和法向力下降。在本文计算外形和来流条件下,腹支杆的存在使α>1.5°以后的轴向力和α>4°以后的法向力反而下降,α<1.5°的轴向力和α<4°的法向力增大。α>1.5°以后由于支杆引起后体弹翼和弹身上载荷的下降,因此压力中心前移。Mach数越高,支杆的影响稍大。     

基于面元法的栅格翼翼身组合体高超声速气动特性研究

栅格翼是由外部框架和内部栅格布置形成的空间多升力面系统,其在高超声速状态具有突出优点。文中结合高超声速流动特点,考虑头部激波及弹身上洗干扰,基于面元法建立栅格翼翼身组合体高超声速气动特性计算模型,根据研究对象特点自动划分面元,进而完成典型栅格翼翼身组合体高超声速气动特性计算,并与数值计算结果进行对比分析。结果表明,在研究范围内,栅格翼翼身组合体高超声速气动特性面元法计算结果与数值计算结果吻合较好,前者较后者相对误差绝对值随攻角增加而增大。     

翼身组合弹箭马格努斯效应数值模拟研究

为研究翼身组合弹箭马格努斯特性及产生机理,采用完全时间相关的非定常N-S方程,对带翼弹箭在高速旋转状态下的绕流场进行数值模拟,得到马格努斯力和力矩系数随攻角变化的规律,所得结果与阿诺德工程发展中心（AEDC）试验及陆军研究实验室（ARL）计算结果吻合很好。分析表明压力差是产生马格努斯力的主因,切应力产生的马格努斯力只占压力产生马格努斯力的1%;弹身马格努斯力除α40°外皆为负,舵马格努斯力始终为正;α5°~30°每个舵的马格努斯力不是正弦变化规律,但合力呈现正弦变化规律。     

一种翼身组合体的气动设计及优化

翼身组合体具有较高的升阻比,可进行较大范围的机动,而且还可以提高落点精度、扩大再入走廊、降低热流峰值并降低过载。采用模线设计方法设计横截面控制点,借鉴航天飞机气动力工程计算方法发展了一套可以预估翼身组合体飞行器纵横向气动力的工程计算方法。提出并建立了翼身组合体飞行器的优化设计模型并进行了计算,获得了带后掠下反翼的翼身组合体优化方案。本方案在5°攻角时升阻比可达6.5,并给出了飞行器稳定配平的质心布置条件。在纵向稳定配平时,组合体飞行器在偏航及滚转方向均为静/动稳定的。研究表明,本方案可在较小攻角时获得较大升阻比,并具有纵横向稳定性,是高超声速机动的潜在可行方案。     

双翼末敏弹尾翼弹性对气动特性影响分析

受安装条件限制,无伞末敏弹尾翼通常很薄,在末敏弹下落过程中,厚度不超过1mm的末敏弹尾翼将在气动力作用下发生一定程度的变形,变形后的尾翼又将对末敏弹气动特性产生影响。文章采用双向流固耦合方法对S-S型双翼末敏弹弹性尾翼进行研究,分析末敏弹气动参数随攻角和速度的变化趋势;采用计算流体力学方法对末敏弹刚性尾翼进行研究。将弹性尾翼与刚性尾翼的末敏弹气动参数进行对比。结果显示,刚性翼、弹性翼末敏弹的阻力系数均随攻角的增大呈准线性递减趋势;刚性翼、弹性翼末敏弹的升力系数均随攻角的增大呈递增趋势。自由飞行试验结果显示,末敏弹的气动参数与弹性翼末敏弹的仿真结果更为贴切。     

超声速飞行器气动伺服弹性稳定性分析

基于分枝模态法,用修正活塞与细长体理论计算超声速飞行器升力面和机身的非定常气动力,考虑了翼身干扰和翼面对舵面的下洗与阻滞作用,通过控制系统输入参数与结构模态参数的关系,输出舵偏产生的非定常气动力的激励关系,建立由模态坐标运动微分方程和控制系统传递函数联合表示的飞行器气动伺服弹性模型,在时域内求解.某飞行器颤振和气动伺服弹性的计算证明该法可行、有效.另外,讨论了结构滤波器参数、敏感元件位置和及弹身固有频率对气动伺服弹性稳定性的影响.     

大长径比卷弧尾翼火箭弹气动特性数值研究

采用CFD方法对大长径比卷弧尾翼火箭弹流场进行了模拟,对比风洞实验数据,验证了数值方法可靠性.对比单个弹体和翼身组合体的气动特性,得到了卷弧尾翼对全弹气动特性的影响.利用正交设计方法,建立以卷弧尾翼6个几何参数为因素的正交数值实验表,得到卷弧尾翼几何参数对各个气动系数的影响.分析零安装角卷弧尾翼在超音速下零度攻角时流场,指出翼面压力分布差异,得到自诱导滚转力矩产生机理,同时分别指出超音速时其随攻角、曲率半径、展弦比和马赫数的变化趋势.     

变结构近空间飞行器大飞行包络控制特性研究

针对近空间飞行器大包络、多任务模式飞行运动,结合通用高超声速飞行器确立了近空间飞行器不同飞行阶段的机体结构,详细建立了气动参数由攻角、马赫数与控制舵面偏转最表示的12状态动力学方程,绘制了不同飞行阶段气动参数随攻角、马赫数的变化曲线,建立了近空间飞行器的大包络运动控制模型,而后研究了不同机体结构与飞行状态下的控制特性和各通道耦合性质,表明该系统可以有效用于全程强鲁棒稳定飞行控制系统的设计与仿真测试,充分体现了近空间飞行器非线性、时变、耦合等飞行运动特点,仿真验证表明了研究结果符合现有分析结果.     

弹性飞行器气动伺服弹性耦合动力学仿真

针对常见的轴对称串置翼布局外形,结构方面以分枝模态法为基础,气动力方面分别采用修正活塞理论和细长体理论计算升力面和旋成体机身的非定常气动力,干扰因子法考虑翼身干扰和串置翼下洗,通过控制系统敏感元件输入参数与结构振型的关系,输出控制力对结构振动的激励关系,考虑弹性振动产生的非定常气动力,建立由模态坐标运动微分方程和控制系统传递函数联合表示的受控弹性飞行器耦合系统动力学的数学模型。在此基础上,用状态空间法将结构运动微分方程和控制系统传递函数转换成状态方程,在时域内用龙格-库塔法实现系统动力学响应的数值仿真,并借以判断系统的动力学稳定性,辨识稳定性的临界参数。应用本文方法完成了算例飞行器无控状态的经典颤振及受控状态的气动伺服弹性分析,证明了方法的可行性和有效性。通过改变控制系统参数和结构参数的仿真计算,总结出这些参数的影响规律。
     

基于压力变送器及智能仪表的箱压自动控制技术

目前发动机地面试验过程中的泵入口压力控制主要是通过控制介质贮箱内压力来实现的.介绍了大型液体火箭发动机研制试验中介质贮箱压力自动控制的一种新方法,这种方法集试验过程信号采集、动态工艺参数显示、上下限设定值显示、报警显示输出及自动控制为一体,减小了手动调节箱压继电器的操作误差,使系统的可靠性得到了很大程度的提高.     

水击问题的Fourier谱方法计算

给出了推进剂供应系统管路内流体瞬变流动的数学模型,提出了采用Fourier谱方法求解瞬变流非线性偏微分方程的新方法。以一段两端分别连接贮箱和阀的直管道为例,利用该方法对阀门关闭后管道中形成的水击和压力振荡特性进行了求解,给出了相应的仿真结果,与已发表的采用特征线法和有限元法求解结果进行了比较。对数值计算中的非物理振荡问题进行了讨论。     

采用简化滤波器思想的卫星姿态确定系统故障诊断

针对星上计算机运算资源有限的问题,为了降低卫星姿态确定系统故障诊断的运算量,提出一种基于卡尔曼滤波器的故障检测与分离方法。该方法首先基于卫星姿态运动方程设计了一种加性卡尔曼滤波器,然后将卡尔曼滤波器与简化观测器思想相结合,进一步提出一种采用简化滤波器思想的姿态敏感器故障诊断律。所提出的故障诊断方法既可以实现对陀螺故障的检测与分离,又能够诊断星敏感器的故障。此外,该方法只利用一个滤波器即可实现故障检测与分离,其计算量小,有利于在轨实施。最终采用一个由三正交一斜装陀螺组件和星敏感器构成的姿态确定系统对所提出的方法进行了仿真校验,仿真结果表明了所提方法的有效性。     

一种相控阵通道的快速测量方法

随着相控阵系统在通信领域的大规模应用,需要供应商提供能够长期稳定工作的设备。而相控阵通道幅相特性在使用的过程中会随着时间产生变化,最终会导致波束性能下降。为了简化相控阵后期维护,降低维护时间和成本,需要在其维护阶段对相控阵通道进行幅相特性测量。这种测量必须由相控阵系统独立完成,且不应该依赖外部环境。传统的测量方法是依次对每个通道进行独立的测量,这种测量方法效率低下,大规模相控阵的测量时间一般都在数十分钟以上,会使通信业务长时间中断,不利于系统快速维护的需求。目前,对相控阵通道幅相的快速测量方法主要是在相控阵天线位于特定的测试环境下进行,目的是加快相控阵的生产周期,不适用于后期维护。将多载波和系统同步结合,提出了一种相控阵通道的快速测量方法。该方法在相同测量精度下,测量所消耗的时间大约比传统串行测量方法少两个数量级,相比于已有的快速测量方法测量时间大幅缩短。最后通过仿真验证了方法的有效性,得出了测量精度和测量时间的关系。并在相同测量精度条件下与传统串行测量方法和已有快速测量方法的测量时间进行对比。     

基于轮廓特征的X射线脉冲星信号多普勒估计

探测器速度引起的多普勒效应会影响脉冲星信号累积轮廓特征,反之,这种轮廓特征改变可用于多普勒估计。提出了一种基于轮廓特征的X射线脉冲星信号多普勒估计方法。通过建立累积轮廓非齐次泊松模型,获得多普勒效应对累积轮廓的影响。构造轮廓特征函数对该影响进行度量,并利用函数值与周期误差的关系搜索信号周期。在搜索过程中,将离散观测数据转换为连续能量密度函数,消除相位间隔对累积周期的限制。根据周期搜索结果,结合脉冲星先验知识完成多普勒估计。仿真实验验证了该方法的有效性,其多普勒估计精度优于频域方法。     

复合材料层合板的可靠性分析方法

基于复合材料层合板的一阶横向剪切变形理论,提出了同时考虑层合板面内和分层破坏的可靠性分析方法。该方法考虑了层间应力对层合板分层的影响,结合Tsai-Hill理论和层合板分层判据,给出了安全余量的表达形式,并考虑了各失效模式之间的相关性。在失效分析过程中,采用蔡氏所提出的刚度退化规律进行刚阵的减缩;利用随机有限元方法对安全余量进行敏度分析,结合改进的一次二阶矩法求解可靠性指标;用改进的分枝限界法寻找主要失效路径;用PNET法计算系统失效概率。计算表明,当考虑分层失效时结构系统失效概率有所增加,这是符合工程实际情况的。因此,设计过程中考虑分层失效是必要的。     

A numerical method to optimize the design of a space inflatable membrane reflector

A numerical method is proposed to optimize the design of a space inflatable membrane reflector. The initial geometry is expressed by polynomial series weighted by a set of shape parameters. The problem is formulated as a minimization of a cost function representing the difference between the effective shape of the reflector and a perfect parabolic surface. The minimization is performed using the Nelder–Mead method or downhill simplex method. The cost function is computed at each vertex of a simplex defined in the space of optimization parameters by solving direct problem thanks to a finite element method. The finite element model handles geometrical non-linearities and takes into account phenomena like membrane wrinkling and torus buckling which may affect the reflector shape when inflated.     

多约束下考虑控制饱和的自适应末制导方法

对地攻击飞行器在末制导过程中需满足导引头视场、攻角、落角等多种约束。针对这一问题提出一种多约束下的制导方法。将多约束下的制导问题转化为具有控制饱和的线性系统的镇定问题。针对该问题,利用基于参量Lyapunov方程的低增益反馈设计方法,推导带低增益参数的纵向制导律。借助Lambert W函数,研究闭环系统状态收敛特性与低增益参数的关系,给出低增益参数的自适应整定方法。通过该方法在线调整参数,使各约束条件得到满足。通过数值仿真对提出的制导方法进行校验。仿真结果表明,提出的制导方法满足约束条件且获得较高的制导精度。     

防空导弹翼面弹性支承边界参数识别

利用系统可测自由度上获得的不完全模态来识别尾翼的边界支承参数。直接由特征方程出发求解边界的未知位移和内力，借此计算边界动刚度，同时对带有一定误差的已知试验模态进行叠代修正，然后通过最小二乘法求解边界的支承刚度和附加质量。本文方法简捷。精度和效率都较高。在防空导弹及其他相关领域都有一定的应用价值。