大范围收敛的摄动Lambert问题新型解法:拟线性化-局部变分迭代法

研究实时、高效、稳定性强的高性能空间轨道计算方法对于中国未来航天工程具有重大应用价值。针对强非线性系统的多维两点边值问题,提出了一种拟线性化-局部变分迭代法（QL-LVIM）,通过拟线性化（QL）思想,将非线性两点边值问题转化为一系列具有一定迭代格式,并且成对出现的初值问题,进而通过局部变分迭代法（LVIM）对其进行求解。利用拟线性化的大范围收敛特性和局部变分迭代法的快收敛、高精度特性,该方法能够在较大的时间和空间尺度下快速精确获得摄动Lambert问题的初速度和转移轨道,其收敛域远大于传统的牛顿打靶法,为航天器轨道转移提供了一种简便高效、稳定性强的新型计算方法。在不同轨道情形下,与几类参考方法对比,结果表明本方法能够在计算效率方面实现大幅提升,并且能够在大范围内实现快速收敛。方法的有效性在地-月系三体问题中得到了进一步验证。     

双基推进剂热分解反应动力学参数的计算

本文根据在CDR-1型热分析仪上测得的一条非等温DSC曲线,确定了双基推进剂热分解反应的最可几机理函数,提出了计算三个动力学参数(E,A和n)的数值方法.该法用线性最小二乘法、迭代法以及二分法与最小二乘法相结合的方法,以积分方程、微分方程和放热速率方程拟合实验数据.在逻辑选择建立了微分形式或积分形式的机理函数的最可几一般表达式后,用放热速率方程计算相应的动力学参数E,A和n的值.     

含裂纹半无限大板面内受集中力的应力强度因子估算

以含裂纹(边缘裂纹、内部裂纹、半无限裂纹)的半无限大板板面内受集中力为例,阐述了用交替迭代法求应力强度因子的过程,并将3种情况的应力强度因子的计算结果绘成曲线,以供参考.     

解大型结构特征值问题的同时迭代的一个新算法

本文提出了同时迭代法的一种变型算法,用来求解大型结构动力分析中的广义特征值问题。使用了E_k子空间和特征方向的概念,证明了该算法的收敛率,并讨论了如何在计算机上实现的某些细节问题。通过统计每步迭代所需运算量和大量的实际计算,表明该算法与当前几种常用的同时迭代法或子空间迭代法相比,收敛率相同,但减少了每步迭代所需的计算量,提高了计算效率。     

含裂纹半无限大板面内采集中力的应力强度因子估算

以含裂纹（边缘裂纹，内部裂纹，半无限裂纹）的半无限大板板面内受集中力为例，阐述了用交替迭代法求应力强度因子的过程，并将3种情况的应力强度因子的计算结果绘成曲线，以供参考。     

二元翼型结冰数值模拟研究

对二元翼型前缘结冰数值模拟进行了研究。采用中心有限体积法求解N—S方程；四阶龙格一库塔法求解水滴运动轨迹方程；逐次二分法寻找水滴撞击极限，然后计算局部水滴收集率；分析翼型表面的传质、传热项，建立相应质量、能量方程并迭代求解，得到翼型表面结冰质量；根据角度二分法生成结冰冰形。预测的NACA0012结冰冰形与实验及Lewice模拟结果吻合很好，结冰过程中的特征参数与冰型及相应冰形特征吻合。     

Kepler方程的一种迭代加速算法

采用Steffensen迭代法求解椭圆轨道下的Kepler方程,克服了不动点迭代法收敛速度慢以及Newton迭代法需要重复计算导数值和迭代结构相对复杂等方面的不足。试验表明利用该方法求解Kepler方程既满足了高精度要求,而且收敛速度快。同时它还能满足不同的偏心率和不同的精度的要求,因而是求解Kepler方程的一种较理想算法。     

减速器机匣安装边的疲劳寿命分析

讨论了减速器机匣安装边疲劳寿命分析的方法。提出了一种安装边载荷计算方法。采用半解析有限元法计算安装边的应力 ,然后采用局部应力应变法计算疲劳寿命。给出了两个算例 ,采用本文方法对一机匣试验件安装边的应力及疲劳寿命进行了计算并与试验结果进行了比较和分析。     

数值解法在燃烧室出口燃气温度计算中的应用

为了求解燃气温度,利用双变量迭代法和赫夫法求解根据物质守恒、化学平衡、压力和能量守恒方程建立的非线性方程组,通过算例验证了2种数值解法求解化学平衡成分和温度的计算结果基本一致,在燃烧室出口温度场试验中,2种数值解法都能很好地满足测试需求,但双变量迭代法的计算收敛更快。试验结果表明:燃气分析的温度计算结果与热电偶的测温结果沿着燃烧室周向变化的趋势基本一致,燃气分析计算的温度基本大于热电偶测试结果。     

多轴应力状态下的局部应力应变近似计算方法及其应用

根据弹塑性力学基本原理 ,建立了单轴与多轴应力应变之间的关系 ,应用多轴Neuber法和Glinka法 ,计算了圆棒结构件应力集中处的应力应变 ,并与弹塑性有限元结果进行了比较 ,分析表明 :多轴应力状态下直接采用单轴近似算法会产生较大误差 ,采用多轴Neuber法及Glinka法计算得到的局部应力应变结果与弹塑性有限元结果比较接近。在此基础上 ,对航空发动机轮盘榫槽部位进行了局部应力应变计算分析     

高维局部非线性转子-轴承动力系统的稳定性和分岔

采用模态缩减方法对具有非解析轴承的高维局部非线性转子系统进行自由度降阶。提出基于变分不等方程的有限元互补解的的数值方法,在几乎不增加计算量的情况下,使得实际轴承油膜力的Jacobian矩阵可与油膜力自身计算同时完成,并取得协调一致的精度。结合打靶法和Floquet理论对实际转子-轴承系统的非线性不平衡响应及其分岔行为进行计算分析,数值结果表明,所提出的方法不仅极大地降低了计算量,而且具有足够高的精度。     

某粉末冶金高温合金涡轮盘破裂转速分析

介绍了平均应力法、平均径向应力修正法、残余变形法、局部塑性应变法等破裂转速预测方法,并采用这些方法对某粉末冶金高温合金涡轮盘破裂转速进行了分析.通过理论预测与试验结果的对比得知:采用平均径向应力修正法、残余变形法、局部塑性应变法预测的破裂转速与试验结果吻合较好,而采用基于径向应力分布均匀假设的平均应力法推算的破裂转速与试验结果相差较大.      

考虑尺寸效应的缺口疲劳寿命预测方法

利用试样缺口根部归一化应力的分布规律，在考虑应力梯度的缺口疲劳寿命预测方法的基础上，着重考虑尺寸效应对缺口疲劳寿命的影响，构建了综合考虑平均应力、应力梯度和尺寸效应的缺口疲劳寿命预测方法.利用已有的TC4合金缺口疲劳试验数据对所发展的方法进行了验证，并与常用的几种局部应力应变法的寿命预测结果进行了对比分析.结果表明:所对比的几种局部应力应变法的寿命预测结果过于保守，而所发展的缺口疲劳寿命预测方法的预测结果几乎全在3倍分散带以内，且大多数预测结果在2倍分散带以内.      

离心式压气机耦合松弛多学科设计优化方法

为提高多学科设计优化的寻优效率,以离心式压气机为对象,开展了耦合松弛多学科设计优化方法的研究.针对松散耦合方法建立的离心式压气机的流-热-固耦合分析模型,为了避免反复迭代导致的分析、优化时间过长,构建了分层次的离心式压气机耦合松弛优化设计系统.在气动优化中采用全局寻优算法探索系统最优气动设计方案;在耦合优化阶段,综合考虑学科间耦合及相互作用,在气动优化的基础上局部寻优校验气动优化结果并得到最终设计方案.这种不同复杂层次模型的优化设计有效缩减了系统优化设计周期.     

航空发动机非线性模型实时计算的迭代方法研究

针对航空发动机非线性模型将作为未来机载应用而需要严格实时性的问题,研究了求解模型的迭代方法,综述其历史发展情况,并讨论收敛性问题,同时提出一种改进应用方案,即模型在过渡态仿真时可将上一状态点雅可比矩阵直接应用于下一点计算,能够在典型应用方案基础上直接减少气动热力模型计算次数,从而提高整体模型执行效率.通过实际算例对这些应用方案进行仿真比较与分析,结果表明:①Broy-den秩1法执行效率最高,更适宜用于模型实时计算;②改进应用方案简单、有效,为进一步满足非线性模型的严格实时性要求奠定基础.      

多场耦合计算平台与高超声速热防护结构传热问题研究

从有限元法(FEM)和有限体积法(FVM)的单元特性出发,提出一种具有局部守恒特性的界面载荷插值方法.采用共享内存技术开发适用于通用有限元和计算流体力学(CFD)软件的多场耦合计算平台,并基于分区耦合方式实现流固耦合传热计算.作为验证,分别将ANSYS与Fluent和CFD-FASTRAN软件进行耦合,计算外壁冷却的喷...     

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

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

基于改进烟花算法的资源分配

资源分配问题作为一个NP-Hard问题,在云计算、无线电、卫星调度、多无人机协同作业等领域皆有研究需求,是一个共性的数学问题。烟花算法作为一种智能优化算法,具有求解大规模资源分配问题的能力,但也存在求解精度低等问题。为了提高传统烟花算法的计算效率和全局寻优能力,提出一种改进烟花算法,用遗传算法中的变异算子替代高斯变异操作,并增加模拟退火流程。最后在多无人机协同作业任务分配数学模型上进行仿真验证,实验结果表明在收敛速度以及计算精度方面,该算法均优于其余3种烟花算法。