高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

预测叶片失速颤振的一种半激盘法

本文利用一种半经验方法——变形半激盘法预测了轴流压气机（风扇）的亚/跨音速失速颤振。该模型的基本方法是:在变形激盘模型的基础上,发展简化的槽道一维非定常流动解。该模型改进了纯激盘模型中的一些不合理性。     

ARMAX模型参数估计的新两步法

 本文提出了一种ARMAX模型参数估计的新两步法。这种方法与其它递推估计方法（如增广矩阵法,极大似然法,Durbin两步法、三步法等）比较,具有较好的收敛特性,对参数的估计,特别是对C参数的估计比较精确,这是利用其它方法,一直未能很好解决的一个难题。     

基于配点法的轨道在线优化方法研究

针对传统轨道优化方法无法应用于轨道在线优化的问题,提出了三阶辛普森配点法结合乘子法求解轨道在线优化问题的方法.首先建立变轨的最优化模型;然后,用三阶辛普森配点法将该最优控制问题转化为受约束的非线性规划问题;最后,采用增广拉格朗日乘子法对转化得来的受约束的非线性规划问题进行求解.仿真结果表明,所提出的方法可以得到高精度的轨道优化结果,并且对状态量和控制量的初值选取不敏感,且仿真具有实时性,计算速度快,可以达到在线进行轨道优化的要求.     

空地导弹滑模制导方法研究

针对空地导弹的制导规律设计问题,给出了滑模制导方法,建立了纵向平面内的弹-目相对运动方程。基于滑模控制理论,设计了空地导弹的制导规律,以同时保证脱靶量和弹道倾角要求,并针对固定目标,给出了相对距离和相对速度的解算方法,介绍了针对移动目标的相对距离和相对速度估算方法。仿真结果表明,该制导方法具有较高的制导精度,并且满足弹道倾角约束。     

关于雨流计数法的讨论及修正

在疲劳载荷的统计处理中,雨流计数法一直被认为是比较合理的循环计数方法,在国内外得到了广泛应用。本文结合对应力-应变响应的研究指出了它的不足,提出了修正雨流计数法。此外,对于现有的两种利用有效系数模拟应力-应变响应的方法也进行了讨论,并提出了修正有效系数法。     

面内弯曲载荷作用下两边简支两边固支复合材料层合板的屈曲

为扩充航空复合材料结构稳定性设计手段,基于Levy形式级数和有限差分法,提出了面内纯弯曲载荷作用下两边简支两边固支（SSCC）复合材料层合板临界屈曲载荷的数值解法。首先,应用有限差分法将挠度函数进行离散化处理;然后利用边界条件扩充方程组,最后将临界屈曲载荷的求解转化为广义特征值问题。通过数值算例分别在各向同性板和复合材料层合板上验证了本文方法,并与有限元方法（FEM）进行了比较。结果表明,本文方法具有很高的精度,为复合材料层合板的屈曲分析开辟了新的途径。     

结构动力学与颤振约束下无人机机翼优化设计

提出了同时考虑结构动力学约束和颤振约束的优化设计方法，建立了利用规划优化算法求解问题的流程，并将其融入大型CAE软件系统HAJIF2013中，研究了优化问题中目标函数和约束条件的选取方式。以某型无人机机翼颤振模型为算例，利用拉盖尔迭代方法自动计算颤振速度，并对优化算法的有效性进行了验证。研究表明，所提出的优化设计方法有效，实现了整个优化求解流程的自动化，对工程设计具有指导意义。     

一种新颖的直升机舵回路故障诊断方法

研究了基于案例推理的直升机舵回路在线故障诊断方法。针对传统的CBR方法一般采用离散数据和符号向量来建立案例的不足,提出了动态系统多源异类信息的案例表达方法,利用系统特性集、征兆集和特征集混合模型作为案例的条件属性。在此基础上,研究了案例检索的综合相似性度量方法。最后,将上述方法在某型直升机的俯仰舵回路进行了试验验证,结果表明,该方法可以诊断出舵回路的控制参数漂移故障,为复杂动态系统的故障诊断提供了一种有效的新途径。     

混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

测试性研制阶段数据评估验证方法

分析了现有测试性验证方法应用情况,阐述了测试性研制阶段数据评估验证方法的必要性。通过与现有测试性验证方法的对比分析,提出了一种新的测试性验证方法,并说明了该验证方法的含义和特点。该验证方法可用于故障检测率、故障隔离率和虚警率的验证。建立了该验证方法的工作框架,包括验证要求、验证程序、验证技术和验证组织。重点阐述了该验证方法的两种工作方式实施流程、适用的最小样本量范围、测试性数据收集表格和测试性数据判据。其中,测试性数据收集表格可以记录故障检测、故障隔离和虚警信息。进行了案例应用,给出了数据分析结果和经验总结。案例应用表明该方法能够达到测试性验证的效果。     

基于有限元全耦合方法的耳片疲劳损伤力学分析

精确分析和预测连接耳片结构的疲劳寿命对保证飞机安全具有重要意义。采用损伤力学—有限元全耦合方法对飞机上典型连接耳片结构进行分析,得到该结构的疲劳损伤演化过程和疲劳寿命;将分析结果与疲劳寿命工程估算方法——名义应力法进行比较,证明损伤力学—有限元全耦合方法应用于结构寿命预测的可行性,并将损伤力学—有限元全耦合方法与线性Miner准则、损伤全解耦方法进行对比分析。结果表明：损伤力学有限元全耦合方法,其计算结果是唯一的、确定的,比运用线性Miner准则和损伤全解耦方法所得的结果更准确,弥补了名义应力法的不足。     

基于蒙特卡罗模拟的飞机结构可靠度分析

飞机可靠度的分析要综合考虑各种因素的影响,由于各因素的影响较为复杂,很难通过解析的方法来计算飞机的可靠度.随着计算机速度的不断提升,通过蒙特卡罗模拟计算飞机结构的可靠性逐渐变得可行.本文基于蒙特卡罗法,将损伤、工作环境、检修周期等作为核心要素,提出了面向工程应用的飞机结构可靠度分析方法,最终得到了飞机全机结构失效概率的计算方法.本方法在计算机上得到了初步实现.     

单子样寿命试验评估的综合概率分布法

提出了大型复杂结构体系或机构系统只做单子样1次寿命试验的评估方法。即据单子样的1次试验值创造性地提出确定当前概率分布及据此作出评估的方法,继而给出了结合半经验算法与当前概率分布算法的加权综合算法。附有算例阐明所提出的方法。     

Implicit discontinuous Galerkin method on agglomerated high-order grids for 3D simulations

High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is inves-tigated in combination with discontinuous Galerkin (DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry represen-tation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomer-ated high-order mesh works well with DG method in dealing with flow problems of curved geome-tries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant effi-ciency improvements compared to its explicit counterparts.     

改进参数化方法及在低空风切变中的应用

研究对参数化方法进行必要的改进,从而使该方法更好地用于求解最优控制问题。改进的基本思路是设法降低利用该方法转换而得的非线性规划子问题的维数。通常在利用参数化方法求解最优控制问题时,由于对控制变量和状态变量同时参数化,从而导致所转化而得的非线性规划子问题维数高,求解难。通过分析认为,只需对控制变量参数化,而对状态变量在时间节点处的值,则可通过对最优控制问题的微分约束进行分段积分求得,从而降低了转换而得的非线性规划子问题的维数,有利于该算法进一步推广应用。利用改进算法求解了飞机在低空风切变中的最优着陆轨迹问题,求解结果表明,改进算法可以顺利求解一类最优控制问题。     

适用于混合网格的约束最小二乘重构方法

为强化边界条件对计算流体力学(CFD)计算中流场重构过程的影响,基于三维任意多面体混合网格,对一种约束最小二乘法进行了研究。通过对不同边界类型的透明化处理,统一了约束方程组的形式,从而简化了该方法的应用。为了充分吸收消去法和加权法各自求解约束方程组的优点,提出了混合采用两种方法的思路。针对层流平板研究了加权法中加权系数如何取值问题,数值实验表明加权系数取到5基本已经足够大。最后通过亚声速层流平板和跨声速湍流ONERA-M6算例对比了混合法、加权法与原始最小二乘法,计算结果表明:约束最小二乘法相比原始最小二乘法的流场重构更加准确,特别是对边界质量较差网格的计算优势更为突出;混合消去加权法相比单独采用加权法的计算结果也有所改善。     

空间柔性网—膜结构分析的数学规划算法

 提出了以参数变分原理为基础的数学规划算法, 以计算分析太空用品中网- 膜结构起皱现象。针对网状柔性结构分析, 算法的原理是将结构化成多个具有一定非线性本构关系的杆单元, 并在此基础上进行杆系结构非线性分析的规划算法构造, 可以较好地模拟网状柔性结构的起皱性质。而对于一般的膜结构, 则构造了用于直接分析的膜单元相对应的本构模型, 并在此基础上构造问题的求解算法。不同于常规的数值计算方法, 提出的算法不依赖于应力的迭代过程。参数变分原理和数学规划算法的应用, 使本算法可以很有效地分析预测网- 膜起皱的区域、模式和应力分布, 从而为解决网- 膜结构起皱问题提供依据。给出的数值算例说明了方法正确与有效性。     

用太阳源法测量遥测系统的品质因数G/T

阐述利用太阳源法测量遥测系统品质因数G/T的原理,提出可用于工程实际的实现方法,重点对太阳流量密度、K1因子、K2因子等参数进行分析,给出工程应用计算方法。最后,以某遥测系统为例,给出了采用太阳源法和信标塔法的测量值,并对测试结果进行比对分析。实测结果分析表明,太阳源法测得的G/T重复测试精度小于0.4 dB,与信标塔法测试结果相差小于0.6 dB。因此,太阳源法合理可行,可广泛应用于遥测系统指标测试。     

自助法中若干问题研究及其在命中精度评估中的应用

在小样本的试验评估中,自助法是一种较优的小子样方法。本文在简要阐述自助法的基础上,深入分析了自助法应用中存在的若干问题,如再抽样样本的获取方式、再抽样样本量、各样本数据容量和适用范围,并给出了自助法的实施步骤,最后比较了自助法和落入概率点小样本方法在命中精度评估中的应用效果,指出了自助法的优点。