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

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

跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

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

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

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

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

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

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

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

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

海杂波的多重分形判定及广义维数谱自动提取

海杂波中的微弱目标检测历来是一个难题,传统的检测方法一直不能得到较好的检测效果。多重分形理论的引入为海杂波中微弱目标检测提供了一个很好的途径。文章提出了多重分形无标度区间自动识别的改进算法,既保留了该算法的客观性,又降低了它的保守程度,使估计的区间更贴近真实无标度区间,同时实现了广义维数的自动计算。实验证明,该方法优于原方法,结果更为精确。     

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

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

跨音叶栅流动的多级有限元分析

将各向异性张力有限单元和Taylor-Galerkin有限元法在多级有限元法的思想下结合起来, 构成了收敛速度和稳定性均较好的适用于跨音流动计算的各向异性多级有限元法。利用这一方法, 基于Euler方程, 研究了平面跨音速叶栅流动, 并与实验结果作了比较。结果表明, 本方法是跨音速计算的强有力手段。      

A compositional method to model dependent failure behavior based on PoF models

In this paper, a new method is developed to model dependent failure behavior among failure mechanisms. Unlike the existing methods, the developed method models the root cause of the dependency explicitly, so that a deterministic model, rather than a probabilistic one, can be established. Three steps comprise the developed method. First, physics-of-failure(PoF) models are utilized to model each failure mechanism. Then, interactions among failure mechanisms are modeled as a combination of three basic relations, competition, superposition and coupling. This is the reason why the method is referred to as ‘‘compositional method". Finally, the PoF models and the interaction model are combined to develop a deterministic model of the dependent failure behavior. As a demonstration, the method is applied on an actual spool and the developed failure behavior model is validated by a wear test. The result demonstrates that the compositional method is an effective way to model dependent failure behavior.     

用于压气机流动计算的3种模型比较

为了更好地对压气机流动进行模拟,在课题组自行开发的结构化有限体积解算器上实现了用于压气机流场计算的混合平面法、谐波平衡法及相滞后法.以NASA Stage 35为例,对3种方法的计算结果进行了比较分析.结果表明:相滞后法的计算精度最高,混合平面法的计算精度最低;相滞后法与半环的双时间推进法结果相近,计算速度提高了20倍;相比混合平面法,谐波平衡法能准确地模拟动静叶间的非定常干涉及进出口参数变化;在谐波阶数达到5阶后,谐波平衡法计算得到的结果不随阶数变化,且与相滞后法的结果基本吻合;混合平面法的计算效率远高于另外两种方法,相滞后法与谐波平衡法在谐波阶数为5阶时的计算效率相当.      

一种适合迭代求解的反馈力浸入边界法

对Goldstein提出的反馈力浸入边界法进行了新的思考，改进了其对力源项的计算，拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时，含有对速度误差的时间积分项，只能用于含时间项的Navier-Stokes （N-S）方程的求解，且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项，避免了时间相关的参数，使其不仅能适合非定常隐式时间推进，还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性，对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算，计算结果均与文献结果符合较好，表明本文改进的方法是有效的。得出的结论为：可以直接基于迭代次数进行反馈力源项的计算，改进的反馈力浸入边界法不仅可与非定常N-S方程结合，进行隐式求解，还可以与定常N-S方程结合用于定常流动的模拟，可将改进的方法运用到更多的流动问题当中。     

Fluid structure interaction simulation of supersonic parachute inflation by an interface tracking method

An Arbitrary Lagrangian–Eulerian (ALE) approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation. A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished. In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact, a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces. In addition, the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method. Based on this method, fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained. Numerical results such as maximum Root Mean Square (RMS) drag, general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results. This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.     

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

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

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

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

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

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

谐波平衡法在非定常流场中的应用

谐波平衡法（HBM）适用于模拟周期性非定常流场。求解这类问题时,只需要一个周期中几个等距时刻的流场解,即可重建整个周期流动的时间历程,计算效率较高;同时还采用谱算子精确求解流场变量对物理时间的导数项,保障了算法的计算精度。鉴于此,详细介绍了谐波平衡法的原理和实现过程,并给出了如何在现有计算程序下实现谐波平衡法的流程。取跨声速翼型和超声速钝锥两个周期性强迫俯仰振荡绕流为例,以双时间步方法为参考,对谐波平衡法的效率和精度作了详细的考核和分析。计算结果表明,谐波平衡法的内存消耗较大,谐波数N_H=2时约为双时间步方法的4倍,但谐波平衡法取得与双时间步方法相近结果的计算时间仅为后者的1/5。因此,谐波平衡法是一种有应用前景的工程实用方法。对于长周期问题,谐波平衡法的优势更加明显。     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

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

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

考虑应力松弛的涡轮盘蠕变-疲劳寿命可靠性分析方法

提出了考虑材料蠕变应变模型参数分散性及应力松弛效应的涡轮盘蠕变-疲劳寿命可靠性分析方法。首先建立蠕变应变概率模型;然后采用响应面法分别回归蠕变和疲劳损伤函数;通过蒙特卡洛法模拟总损伤分布及寿命可靠度。同时提出给定可靠度的寿命计算方法和两种考虑应力松弛的涡轮盘蠕变-疲劳寿命可靠度分析简化方法。通过算例验证了方法的有效性。