脉冲子结构与有限元刚-弹混合连接的子结构方法

航天器结构的日益复杂庞大对系统级的动力学建模仿真以及进一步的结构优化提出了巨大挑战。为提高动力学求解效率,通常引用动态子结构方法。本文利用适于处理瞬态冲击问题的脉冲子结构(IBS)方法,并对其进行改进,将基于脉冲响应函数(IRF)的子结构与有限元建立的子结构综合,同时考虑刚性、弹性以及刚-弹混合连接情形下的子结构综合格式。通过3个数值算例,验证了方法的正确性。最后将刚-弹混合连接下的子结构方法应用到月球探测器软着陆的动态响应预测,结果表明该方法适于对月球探测器软着陆动态响应进行高精度快速预测,并且可以应用于月球探测器的局部动力学结构优化。     

基于径向基函数插值的气动弹性计算方法

动网格与界面数据传递方法是气动弹性计算的主要关键技术.将径向基函数（RBF,Radial Basis Function）插值引入气动弹性计算,介绍了径向基函数的概念并给出其应用于动网格及界面数据传递的基本方法.对AGARD 445.6机翼颤振问题进行了计算流体动力学（CFD,Computational Fluid Dynamics）/计算结构动力学（CSD,Computational Structure Dynamics）耦合计算,在不同的来流速度条件下得到了振动收敛、等幅谐振及振动发散的现象,临界颤振速度及颤振频率与实验数据符合良好,验证了径向基函数插值方法用于气动弹性计算的有效性.      

多饱和约束拦截弹参数时变系统控制器设计

针对拦截弹参数时变和执行机构饱和非线性影响控制系统性能,提出基于齐次多项式参数相关 Lyapunov 理论和非线性矩阵微分方程凸多面体转换算法相结合的鲁棒最优控制系统设计方案。 通过引入标准饱和约束算子,建立控制系统数学模型并给出控制系统结构;引入饱和归一化因子,把非线性参数时变微分方程转化成线性多胞型微分方程;基于齐次多项式参数相关Lyapunov理论和拓展Pólya’s定理,把参数时变多胞型微分方程转换成齐次多项式参数相关微分方程,并给出齐次多项式参数相关控制器结构;利用迭代算法把非线性矩阵不等式转换成线性矩阵不等式(LMIs),进而基于LMI凸优化理论获得控制器解,设计出的控制器是关于齐次多项式阶次g和Pólya’s松弛度d的函数;通过定点和制导控制系统仿真校验了该算法,仿真结果表明随着齐次多项式阶次g和Pólya’s松弛度d的增加,设计出的控制器在确保系统稳定的前提下改善控制系统性能。     

一种并行组合采样系统信号重构方法

提出了在并行组合采样系统这种多通道组合采样系统中,对于时基非同步信号,使用插值算法对组合信号进行重构的方法。实际的并行组合采样中,由于各采样通道之间的采样时钟相位控制存在时基偏差,引起了采样信号的时间非同步,导致采样波形非线性失真。为补偿系统的无杂散动态范围和实际信噪比,有效地还原信号原貌,需要对采样序列进行信号重构。提出了用插值算法来重构信号,实现对时基偏差的校正,并通过仿真验证,分析了此方法的可行性。     

基于引射器函数法的程序优化设计及验证

针对考虑多因素综合影响的引射器优化问题,基于引射器函数法,深入研究了静压协调函数的数学特性并分析了其数学曲线上的奇异点。在此基础上,通过程序设计和基准推进归纳法得到了多分支工作特性曲线,分析了不同分支下解的特性。另外结合工作特性曲线重点研究了混合室背压、主次流总压比和混合不均匀度对引射性能的影响,提出了静压特性曲面、临界曲线的设计概念。通过试验数据对比分析验证了引射器函数法的可靠性。结论表明:靠近临界曲线工况时引射器设计性能较好;混合不均匀度对引射性能影响重大,不均匀度为1.5时对比理想状态,引射系数最大误差达到32.73％;考虑壁面摩擦设计时需对摩擦因数公式模型进行修正。研究结果为引射器优化设计提供了重要指导。      

基于火焰图像诊断的模型燃烧室燃烧不稳定特性

提出了一种基于火焰图像诊断的旋流燃烧室燃烧不稳定特性研究方法。分析甲基滤镜下高速CCD （charge coupled device）拍摄的火焰面振荡图像,发现不稳定工况下的旋流火焰在以较高频率抖动的同时伴随着周期性的火焰脱落,其脱落信息可由火焰面轴向峰值位置体现。对CCD图像进行数字图像处理获得轮廓线,统计其边缘轴向最远处坐标作为火焰抖动峰值位置。基于图像分析确定火焰脱落阈值,获得不同工况下火焰脱落发生时的脉冲时序信号,统计得各工况火焰脱落频率。研究发现该频率分别与火焰抖动时序、CH^*时序及燃烧室压力脉动时序信号各自3阶本征模态函数的傅里叶谱主峰值接近。在多工况下研究了不稳定火焰脱离频率在时序信号中的存在方式和规律,发现对于本模型燃烧室,火焰脱落频率随着油气比的降低而升高,在油气比约0.0096处到达火焰脱落频率临界点,之后呈现缓慢下降趋势。     

基于频响主元的卫星结构损伤定位与评估方法

鉴于目前国内卫星预、复振频响结果的评估主要依靠设计师的经验分析,缺乏高效准确的量化手段,文章基于频响函数(FRF)的主元分析(PCA),提出了两种针对卫星结构的损伤定位指标,分别是多主元综合指标和FRF主元残差均方根;并采用平方预测误差(SPE)指数对卫星结构的整体损伤程度进行评估。以卫星主承力部件中心承力筒为分析对象,验证了基于卫星结构频响主元的损伤定位和评估方法的有效性。文章对基于预、复振对比分析的损伤识别、量化试验评估标准、提高卫星研制效率有促进作用。     

基于物理量梯度修正的RBF数据传递方法

径向基函数是复杂外形飞行器多场耦合计算中重要的数据传递方法。复杂飞行器物理量分布通常与局部流动特征密切相关且在不同区域变化剧烈,如何在具备各向同性特征的径向基函数中考虑实际物理量各向异性分布因素的影响成为制约其精度进一步提升的关键。针对径向基函数的问题,提出了一种基于局部物理量梯度对径向基三方向权重进行自适应修正的数据传递新方法。采用高超声速舵面、翼身组合体算例对该方法在单次数据传递中的可行性及效果进行了验证,同时基于该方法对压缩拐角外形开展了气动力/热/结构多场耦合数值模拟,并与风洞试验结果进行了对比分析,以验证该方法在实际耦合问题计算长时间高频次数据传递中的适用性及可靠性。研究结果表明该方法能够很好地提升如薄板样条插值（TPS）、多重二次函数插值（MQ）等全域基函数的单次插值效果,并且有效地改进了MQ方法对形态参数的鲁棒性,而对于梯度修正后的紧支基函数,其改进效果更佳,以较少的点达到了TPS、MQ等全域基函数的插值效果,实现了数据传递效率和精度的兼顾。采用该方法获得的耦合计算结果与试验对比一致,表明了该方法的可用性以及在复杂外形多场耦合问题中的良好工程应用前景。     

奇异积分方程方法在接触压力分析中的应用

总结回顾了奇异积分方程的数值解法,对于第二类奇异积分方程使用分段连续函数法进行了求解;对于两个弹性体接触问题,通过接触体之间的滑移函数和间隙函数,建立求解接触压力的奇异积分方程.分别针对圆柱体与弹性半空间体、抛物线型压头与弹性半空间体、榫头与榫槽3类接触问题,确立奇异积分方程的具体表达式,而后使用分段连续函数方法进行求解,获得接触面上的接触压力.最后将计算所得的接触压力分别与理论解和有限元解进行了对比.对于圆柱体与弹性半空间体接触问题,奇异积分方程法的最大接触压力与理论解和有限元解的相对误差分别为0.3%和0.5%;对于榫头与榫槽接触问题,奇异积分方程方法计算所得的最大接触压力与有限元解的相对误差为1.8%,验证了奇异积分方程方法的有效性.      

基于改进多目标遗传算法的舰尾紊流模拟方法

为提高舰尾紊流自由大气紊流分量仿真的可信度,提出了一种紊流数值模拟的新方法.首先,使用带有修正系数的Euler前向差分格式表示各个方向的紊流序列,同时结合智能算法的思想,把紊流相关性检验中的均方差误差和相关函数误差作为待优化目标函数,将修正系数的选择看成一个多目标优化问题,并采用改进的多目标遗传算法进行求解.最后,通过仿真算例验证了本文方法的正确性与合理性,计算结果表明该方法可以根据不同的采样步长灵活地生成所需紊流.尤其在小步长情况下,亦可得到很好符合理论值的紊流序列,可以满足虚拟飞行实验的要求.