孔冷挤压有限元仿真中的铰削分界面位置确定方法

开展孔冷挤压过程有限元仿真计算是残余应力分布获取和疲劳寿命预测的前提。在有限元建模阶段，设置铰削层单元与基体材料单元之间的分界面，是模拟铰制终孔工艺过程的关键。通过弹塑性力学分析，建立了挤压强化过程芯棒、衬套和被挤压强化连接孔的应力分析方法；基于分析中得到的不同位置处微单元的径向位移量，建立了铰削分界面相对位置计算模型。并开展了关键参数的敏感性分析，定量研究了关键参数变化对残余应力分布和径向位移量的影响程度。本工作为孔冷挤压强化有限元模型建立中，铰削层单元与基体材料单元分界面相对位置确定，提供了便捷可靠的方法。     

再入飞行器鼻锥逆向喷流对流场及气动热的影响

 使用计算流体力学（CFD）方法研究逆向喷流热防护系统对降低再入飞行器鼻锥物面热流的效果,获得了流场参数,回流再附点位置,物面压力分布以及热流分布。分析了逆向喷流对降低物面热流的物理机理,喷流通过与来流相互作用形成马赫盘,将来流导流到四周,不与物面直接作用形成气动加热,同时喷流回流形成低温区,降低物面与接触气体的温差,进而降低了物面热流。随着总压比率增大,这种效果越明显,气动加热越轻。为更合理分析喷流强度对流场及传热量的影响,将总压比率和流量相结合,提出了新的参数R PA。分析该参数的应用效果,结果发现不同的流量与总压比率组合成相同的参数R PA,可以实现相同的激波位置、再附点位置、表面热流峰值位置和总传热量。这说明该参数可用于表征喷流强度,用以分析喷流对流场及传热量的影响。     

基于表面压力信息的空间流向涡识别方法

旋涡与表面的相互作用广泛存在于各类飞行器的绕流中，旋涡影响飞行器表面的压力分布，引起其气动特性的改变。而表面压力分布同样反映飞行器绕流中复杂涡系的空间流动特征，结合迎角、侧滑角等来流参数，可以判断飞行器受力状态和运动趋势。以平面点涡和“镜像涡”理论为基础，建立基于表面展向压力分布曲线的空间涡识别方法。根据截面展向压力分布曲线是否存在因主涡诱导下洗气流产生的正压区，定义流向旋涡的近物面流动和远物面流动；根据二次涡相对于主涡的位置，定义压力分布曲线的“近二次涡侧”和“远二次涡侧”。利用展向压力分布曲线“远二次涡侧”的1/4峰值、峰值及其展向位置，识别流向旋涡空间位置特征和强度特征。搭建涡-面相互作用试验平台，比较空间流场测量结果和表面压力信息识别结果，验证该方法有效性。研究结果表明：通过表面压力分布曲线可以辨识流向旋涡的空间位置和强度特征，空间流场测量与基于表面压力信息的旋涡识别结果的关联性分析验证了该方法的有效性。为重构飞行器周围旋涡流动结构，以及实现飞行器气动力的预测奠定了重要的技术基础。     

涡方法分析并列圆柱的旋涡脱落现象

利用快速离散涡方法及流场可视化技术分析二维粘性不可压缩均匀流中并列圆柱体的旋涡脱落现象，识别作用在圆柱上的气动力及旋涡脱落频率。使用边界元方法及引入一个新颖的涡量边界条件计算表面涡量，利用系数影响矩阵技术获得无网格方法，采用自适应快速多极子算法快速计算涡元速度。数值算例再现了不同间距并列圆柱体的同步同相旋涡脱落、同步反相旋涡脱落、尾流的不对称性及间隙流的半稳定性等现象，数值结果与早期实验结果吻合得很好。     

卫星姿控发动机混合物羽流场分区耦合计算研究

研究求解喷管内流场N-S方程数值计算方法,发展基于N-S方程物面边界滑移流理论计算技术。提出求解羽流核心区轴对称DSMC模拟方法与远场三维DSMC仿真方案,发展多组元混合物羽流DSMC仿真方法。研究求解卫星姿控发动机内外近场、远场、倒流区和物面相互作用影响区多流域流场分区耦合计算技术,建立了一套用于求解混合物燃气羽流及对太阳电池帆板与卫星体表面撞击污染影响数值模拟方法。通过对分别安装于某在轨卫星不同位置两个典型姿控发动机燃气五组元混合物羽流计算研究及相关结果对比分析,证实本文数值方法可靠性。     

一种可任意给定环量分布的螺旋桨设计方法

提出了一种螺旋桨快速设计方法，该方法可以根据任意给定的环量分布及工况快速设计出高效率螺旋桨的几何外形（桨叶弦长、扭转角分布）。为对比有叶素阻力与无叶素阻力的最佳环量分布形式，首先对计入和不计叶素阻力下的最佳环量分布表达式进行了推导，分析了各参数对最佳环量分布形状的影响。然后基于片条理论进行了逆向推导，建立了以环量分布为输入的螺旋桨快速设计方法，并分析了螺旋桨理想效率的影响因素，发现理想效率的水平与来流动压、桨盘载荷有关。最后进行了螺旋桨设计实践，对比分析了不同环量分布下螺旋桨性能的差异。结果表明，所提出的设计方法能够根据给定环量分布进行精确设计，相对误差不超过7%，设计性能与参考文献中设计结果的差别不超过2%。      

水下单信标导航航路规划和最优航路接近方法

相比于传统的长基线和超短基线等导航方式，水下单信标导航具有布放简单的优点，但其导航精度有待进一步提高。为此，提出了单信标导航的航路规划方案，通过泰勒级数展开推导了水平位置精度因子的表达式，并分析了导航点和声信标的相对几何位置关系对导航精度的影响，最终提出了航路规划方案。在此基础上，还提出了自主水下航行器(autonomous underwater vehicle，AUV)从其他位置接近最优航路的方法，包括两部分：一是建立以时延为观测量的滤波模型，利用滤波算法实时获取AUV的位置估计值；二是基于可观测度分析结果对最优航路接近过程的轨迹进行了设计。二者的结合使得AUV高效率且高精度地逼近最优航路。仿真证明了采用所提出的航路规划方案和最优航路接近方法可以提高导航精度，航路规划方案使得AUV位置估计的均方根误差近似为2 m。     

不同来流条件对涡轮叶片表面颗粒沉积影响的实验研究

为了研究微颗粒在涡轮叶片表面的沉积特性，采用熔融态石蜡代替熔融态微颗粒，通过石蜡在平板表面的沉积来模拟真实涡轮叶片表面的微颗粒沉积，从而开展较低温度条件下微颗粒的沉积实验。通过实验观察石蜡在试验件表面附着固化分布，同时测量石蜡在平板表面不同位置的沉积量，研究了不同来流攻角、颗粒浓度对平板表面颗粒物沉积的影响规律，分析了等效加速实验的颗粒物沉积效果。结果表明：颗粒物沉积主要发生在平板前缘及压力面的中后部，前缘的颗粒沉积物质量远大于压力面和吸力面的。随着攻角的增大，压力面和吸力面沉积量显著增大，前缘变化不明显；随着来流颗粒浓度的增加，前缘和压力面的沉积量增加；对于已有等效加速实验准则，保持总粒子喷射质量相同，在颗粒浓度升高，实验时间变短的情况下，平板前缘位置的沉积量增大。     

一种基于输出误差观测的冗余MEMS-IMU标定技术研究

针对传统六位置标定存在标定步骤复杂、标定时间长等问题，同时随着MEMS器件冗余数目增加和冗余配置结构复杂化，利用分立标定技术实现器件误差参数辨识的难度进一步增大，且不同配置结构所采用的标定方法存在通用性较差的问题。因此提出了一种基于Kalman滤波的冗余MEMS-IMU分立标定方案。该方案首先采用小角度建模法实现安装误差的精确建模；然后针对直接以转台三轴角速率为观测值，导致部分状态量不可观的问题，提出了以器件的输出误差值作为观测量、以器件误差参数作为状态量设计Kalman滤波器；最后设计了高精度三轴转台转位编排方式，并利用四陀螺冗余结构进行标定仿真试验。仿真结果表明：该标定方法与六位置标定方案相比，标定精度平均提高了11.37%，可实现MEMS器件误差参数的快速辨识，对实际工程实践具有一定参考价值。     

基于缺陷概率特点的粉末冶金材料寿命预测概率模型

粉末冶金材料由于其制造工艺的特点,强度和寿命对微缺陷（夹杂、气孔、表面划伤）十分敏感,导致了粉末冶金材料的破坏具有较大的分散性,使得寿命预测更为困难。本文根据国内粉末材料中缺陷的分布特点,基于文献[79]中概率断裂分析的思路,对原有方法进行了修正和推广,重新给出了缺陷在表面、亚表面以及内部时的定义及缺陷在这些不同位置时出现概率断裂的确定方法,建立了一个可考虑缺陷形状、大小、位置等分布特征的粉末冶金材料寿命预测概率模型,并基于此计算了不同尺寸缺陷位于不同位置时材料的失效概率及总失效概率。分析表明:所给出的方法可以很好的表征国内工艺条件下粉末冶金材料缺陷的概率特征对强度寿命的影响,方法是有效的。      

侧向喷流与飞行器三维绕流干扰流场特性的DSMC/EPSM数值模拟

采用DSMC/EPSM混合算法,研究了带有侧向喷流的复杂外形高超声速飞行器稀薄区和连续区三维混合流场的干扰特性。采用一种飞行器物面网格与DSMC计算域网格分别标识的方法,通过判断模拟分子与表面碰撞来完成飞行器物面网格与DSMC计算域网格间的信息传递和信息存贮,对于复杂外形飞行器的精确描述的物面网格不需做进一步处理,直接应用于不依赖于飞行器外形的DSMC计算的通用子程序中。引入碰撞数做为混合参数,对流动区域进行划分,并通过网格中碰撞数的统计和处理,有效地将DSMC和EPSM方法结合在一起。仿真计算了三维复杂外形飞行器带喷流流场压力、热流分布量,飞行器表面气动力、气动力矩和气动热参数,对喷流与高超声速空气绕流相互作用以及它们与物面之间的相互扰动进行了分析,证明了采用的方案和技术的有效性。     

Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor

A robust unsteady rotor flowfield solver CLORNS code is established to predict the complex unsteady aerodynamic characteristics of rotor flowfield. In order to handle the difficult problem about grid generation around rotor with complex aerodynamic shape in this CFD code,a parameterized grid generated method is established, and the moving-embedded grids are constructed by several proposed universal methods. In this work, the unsteady Reynolds-Averaged Navier-Stokes(RANS) equations with Spalart-Allmaras are selected as the governing equations to predict the unsteady flowfield of helicopter rotor. The discretization of convective fluxes is accomplished by employing the second-order central difference scheme, third-order MUSCL-Roe scheme, and fifth-order WENO-Roe scheme. Aimed at simulating the unsteady aerodynamic characteristics of helicopter rotor, the dual-time scheme with implicit LU-SGS scheme is employed to accomplish the temporal discretization. In order to improve the computational efficiency of holecells and donor elements searching of the moving-embedded grid technology, the ‘‘disturbance diffraction method" and ‘‘minimum distance scheme of donor elements method" are established in this work. To improve the computational efficiency, Message Passing Interface(MPI) parallel method based on subdivision of grid, local preconditioning method and Full Approximation Storage(FAS) multi-grid method are combined in this code. By comparison of the numerical results simulated by CLORNS code with test data, it is illustrated that the present code could simulate the aerodynamic loads and aerodynamic noise characteristics of helicopter rotor accurately.     

直接统计模拟位置元方法中的分子表面反射确定论判据

在直接模拟Ｍｏｎｔｅ－Ｃａｒｌｏ（ＤＳＭＣ）方法的位置元方法中准确判断分子是否与某一表面元相撞,直接关系着气动力热的计算精度,是方法需要解决的关键问题。本文发展的分子表面反射判据,使这个问题得到了有效的解决。此外,还讨论了分子空间位置坐标的记录方法、分子碰撞对的极限距离、网格自适应等实际模拟技术。与圆球自由分子流准确解和过渡领域贴体网格ＤＳＭＣ模拟结果的比较,证明了上述判据和技术的有效性。     

非结构网格生成方法的改进及气膜冷却三维静子流场的求解

讨论五面体网格下,涡轮静子叶片含多排冷却孔时三维粘性流场雷诺平均Navier-Stokes方程组的数值求解问题。考虑到叶片含多排冷却孔的工程背景,本文发展了生成五面体非结构网格的快速有效方法即仅在一个拟S₁面上生成非结构网格,而沿叶高方向仍采用结构网格的生成办法。文中还采用了新的网格优化措施以改善生成网格的质量。流场的计算采用了非结构网格下的高分辨率迎风格式。几个典型算例表明:本文发展的网格生成方法快速有效,而且网格生成质量较高;发展的计算格式能够获得较为满意的三维流场数值解,并具有较高的激波分辨率。      

基于闭路制导的弹道导弹发射诸元的仿真算法

提出了一种基于闭路制导的弹道导弹的发射诸元算法.给定发射点与瞄准点的地理位置,以压低弹道为目标,通过一种迭代方法求解关机点弹道倾角,并基于地球引力级数展开模型提出了一种虚拟目标的快速计算方法,以主动段弹道平面和被动段弹道平面二面角为指标完成了发射方位角的迭代计算,最后通过仿真算例论证了该算法的可行性.     

基于DSP+FPGA的电机伺服控制

舵机控制系统的数字伺服控制具有高可靠性、调试灵活、易于扩展等优点,通过对某舵机控制系统中数字伺服控制器功能分析,提出采用DSP+FPGA的设计方案,由DSP完成双余度无刷直流电动舵机的控制率解算、余度切换;FPGA逻辑设计完成HALL信号解码、PWM输出、A/D采集。设计有效地减少了DSP的计算和管理任务,具有极高的通用性和可扩展性,并极大地简化了系统开发过程,对舵机系统的数字伺服控制设计具有一定的指导意义。     

分析计算航天飞机气动系数的边界元局部方法

空间运载工具设计阶段为分析比较不同外形高超音速各流动领域的性能,广泛应用局部方法。在其原始形式下,形状函数计算繁复,并限于简单几何外形组成的轴对称体。用边界元局部方法进行形状函数计算,可免去原局部方法面积分的繁复和藉助表格的局限。对STS-1轨道飞行器外形进行计算得到丁形状函数,并根据实验数据确定领域系数。由此计算出飞行器的法向力、轴向力和俯仰力矩系数,计算结果与实验数据相符。     

气冷叶片表面局部对流换热系数的数值计算

本文提出了联立求解二维可压缩边界层方程组和稳定导热方程的迭代法,从而修正了文献[1]假定物体表面温度等于来流总温的局限;同时提供了用这种方法计算气冷叶片表面局部对流换热系数的完整的FORTRAN程序。文中用近似积分法求解边界层方程组,用有限元素法计算叶片的稳定温度场。算例的结果与有关的理论及实验结果相接近。     

Numerical Simulation of Rotor Flow Field Based on Overset Grids and Several Spatial and Temporal Discretization Schemes

A numerical method based on solutions of Euler/Navier-Stokes (N-S) equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method (AUSM) scheme, upwind AUSM/van Leer scheme, AUSM+ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points (IGBPs). Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.     

飞机重心位置的计算分析方法

飞机重心位置对飞机的安定性和操纵性影响很大.不在飞机重心允许范围内装货或载客时,会使飞机操纵性和安定性变差,影响飞行安全.飞机重心位置计算分析是新型飞机研制、现役飞机加改(换)装设计等工程中必不可少的一个重要环节.以某型现役飞机进口短波单边带电台换装国产短波自适应/跳频电台工程作为研究对象,研究了基于力矩平衡法的飞机重心位置计算方法,计算了该型飞机换装短波电台后空机重心位置.计算结果表明,该型飞机进口短波电台换装国产短波电台后,空机状态下的重心位置满足重心允许范围要求.