梁壳组合结构的一种非线性有限元计算方法

提出了梁壳垂直组合结构在有限元计算中存在的问题,并以几何大变形问题为依托提供了解决该问题的计算原理和有限元方法。数值算例证明了该理论和方法的正确性。该法不仅适用于梁壳组合结构的大变形非线性计算,也适用于静动力小变形问题的线弹性分析计算。     

某直升机“空中共振”分析

直升机“空中共振”问题,是一种危及直升机飞行安全的自激振动问题,是直升机设计中必须解决的问题。某直升机“空中共振”计算,是指该直升机在悬停状态的“空中共振”计算。在计算过程中,应用“八五”重点预研课题《星形柔性桨毂直升机“空中共振”计算方法研究》的研究成果,建立了“空中共振”计算模型,研制了一套实用的计算软件,分别采用临界稳定性分析法和特征根分析法进行“空中共振”计算,并进行了计算结果的对比,得出了明确的计算结论。     

粘性极低速流动的预处理方法研究

讨论了采用预处理方法计算低速流动时遇到的参数选择、计算格式、舍入误差等问题,提出了一种适用于粘性计算的预处理参数计算方法.计算结果表明,使用高速流的计算格式计算低速问题时,即使残值收敛,也有可能得不到正确的计算结果;AUSM+-up格式适用于低速计算;马赫数小于等于10-6时,舍入误差会造成数值解无法收敛.     

一维喷管化学反应流的吉尔方法研究

建立了一维喷管化学反应流计算模型,并进行了数值求解,给出了喷管内一维化学反应流计算结果。采用吉尔方法解决控制方程组的刚性问题获得满意效果。计算中解决了化学反应流计算中初值难于给定和喉部条件难于确定等问题。并通过对氢氧发动机喷管的分析计算,证明了这种方法的正确性。     

求解航空发动机非线性方程组的变步长牛顿法

详细分析了航空涡轮发动机部件共同工作的非线性方程组的求解收敛性问题,并对计算中常见的迭代不收敛问题的机理进行了研究总结。针对常见的不收敛问题,提出了以传统牛顿迭代法为基础的变步长牛顿法加以解决。采用变步长牛顿法后,对某型定几何单轴涡喷发动机进行了一系列的验证计算,计算结果表明计算收敛性得到较大幅度的改善。变步长牛顿法对改善发动机总体性能计算收敛性具有很好的参考价值。     

公差原则中相关要求的计算问题

关于公差原则中的相关要求的计算问题 ,本文提出了两种方法 :计算法和图解法。运用这两种方法可以简化此类问题的计算过程 ,准确求得计算结果。     

多面体网格在静气动弹性计算中的应用

介绍了多面体网格技术在静气动弹性问题中的应用.基于M6算例和径向基函数插值法,探究了多面体网格的CFD计算效率和变形能力,使用AGARD445.6机翼进行了静气动弹性数值计算.结果表明,多面体网格具有计算效率高、变形能力强、计算结果准确等优点,可适用于静气动弹性问题的数值计算.     

起落架坠毁性能计算研究

对目前国内外起落架坠毁性能计算情况进行了分析,针对分析中发现的一些问题和疑问,对起落架坠毁性能计算中轮胎压死后计算和起落架结构刚度影响两个问题进行了研究,通过研究确定起落架结构刚度是起落架坠毁性能计算不可忽略的重要因素。     

直升机桨叶扬起下坠碰撞动响应计算

针对铰接式旋翼桨叶扬起下坠碰撞问题,采用中等变形梁模型,同时引入有限转角模拟桨叶绕铰的刚性运动。运用Hamilton原理建立桨叶扬起下坠碰撞动力学方程,采用隐式Newmark法计算桨叶在有初始挥舞角下的下坠碰撞响应。计算结果与试验数据的对比表明:如果采用定步长的方法计算响应,计算时间较长,计算结果会产生明显的延迟。为了解决上述问题,在碰撞前后设置了一碰撞缓冲区,采用了变步长的方法强制计算点在跳过此缓冲区时必须进入该区。此法在计算桨叶扬起下坠碰撞动响应时,能较好的解决计算延迟问题,计算时间大大减少,计算结果与试验数据吻合较好。     

一种飞机交流电网计算新方法

提出了一种飞机交流电网计算新方法。该方法基于汇流条流入功率等于流出功率得到的功率平衡方程,将交流电网计算问题转化为求解非线性方程组的问题,采用牛顿-拉夫逊法进行求解计算,并采用MATLAB编制了相应的交流电网计算程序。以民用飞机交流电网为例,将本文算法和常规方法的计算结果进行对比,结果表明本文算法可以提高计算结果的准确性。     

罚函数法在求解航空发动机非线性方程组中的应用

针对求解分排涡扇发动机非线性方程组过程中出现的喷管进入无意义工作范围导致计算发散的问题.在传统的牛顿法中结合罚函数法,提出了一种新的非线性方程组求解方法.对某大涵道比分排涡扇发动机进行了仿真计算验证,并且将计算结果、收敛情况和传统的牛顿法进行了对比.结果表明,能够较好地处理喷管进入无意义工作范围导致计算发散的问题,改善了计算的收敛性.     

多级轴流压气机二维性能预测方法

为了满足多级轴流压气机性能预估需求,在一个流线曲率程序基础上开展了经验、半经验关系式研究。通过对文献中的损失模型进行校验及融合,建立了相匹配的损失计算模块,并改进了端壁损失计算方法;研究了利用S₁流面计算修正S₂正问题的方法,解决了基于传统平面叶栅试验数据的攻角、落后角模型与先进技术叶型之间不匹配的问题,继而发展出了一个高精度的S₂正问题计算方法。为了验证计算方法,利用3个不同负荷水平的、经试验验证的多级压气机进行了校验计算。对比表明,发展的程序对多级压气机具有很高的计算精度和稳定性,可用于多级轴流压气机性能分析。     

外挂物连续投放和非对称投放动载荷计算

 以求解广义动力学方程为基础,分析了飞机在外挂物弹射投放时所引起的动载荷问题,实现了用直接矩阵提取方法（DMAP）计算外挂物连续投放动载荷;用半机模型计算全机外挂物非对称投放动载荷;以及未投放外挂物的多自由度动载荷计算问题。计算结果表明,该方法能有效地解决飞机的各种动载荷问题。     

轮盘振动特性计算

本文提出了一种计算轮盘振动特性的数值解法。采用薄板和厚板两种理论,将问题化成一阶微分方程组进行数值计算。编制了计算机程序,得出了令人满意的结果。由于没有涉及到复杂的数学问题,计算机程序大为简化,计算时间少。      

重型燃气轮机透平气动性能研究

为提高某重型燃气轮机单轴多级透平气动性能,利吲航空发动机设计及校核程序,对某重型燃气轮机透平的设计进行了反演算,利用1维平均半径流场的计算方法进行了计算分析,得出了气动参数的展向分布,并进行了透平冷、热态流道的计算对比,分析了大尺寸透平尺寸换算的变化规律,通过采用变比热、多损失模拟的准3维正问题计算方法,对结果进行了校核计算。结果表明：采用现有的航宅发动机设计平台进行重型燃机单转子多级透平反问题设计和正问题校核计算误差较小,其技术可靠,为进一步利用该体系完成透平性能改进打下了基础。     

基于两步优化策略的结构网格负载平衡

国家数值风洞（NNW）工程旨在发展完全自主知识产权的计算流体力学（CFD）软件,结构网格负载平衡问题研究是该工程中的一个重要组成部分。本文发展了两步优化策略以求解结构化网格的负载平衡问题。第1步优化采用传统的贪婪算法,完成对大块网格的剖分和以进程计算时间为指标的网格块分配;第2步采用遗传算法（GA）,目标函数兼顾进程计算时间和通信时间,在第1步优化结果的基础上,对网格块在进程上的分配开展二次优化。为准确计算GA的目标函数,构建了一套计算时间和通信时间的建模方法,包括样本生成、模型建立和模型验证,整体方法具有一定的通用性。根据负载平衡问题以及两步优化策略的特点,对GA的编码、交叉、变异和种群初始化进行了研究,详细分析了交叉操作的递归问题及解决方法。算例验证说明建立的进程计算时间和通信时间模型具有较高的计算精度,能够用于GA的目标函数计算;两步优化策略能够在第1步优化的基础上进一步改善优化结果,从而减少CFD问题的整体计算时间,对于计算量巨大的工程问题具有较大的实用价值。     

有热解的轴对称烧蚀体传热有限元计算

本文采用有限元法,分析和计算了有热解的再入烧蚀体的传热问题。计算方法和程序先用两个算例进行了考证,最后给出了沿真实弹道的计算结果。     

复约束方法计算斜齿圆柱齿轮结构特征值问题

根据齿轮结构的循环对称性，以仅含有一个完整轮齿的扇形子结构作为齿轮有限元计算模型，用计算循环对称结构特征值问题的复约束方法，通过求解一个扇形子结构的特征值问题得到整个齿轮结构的特征值和特征向量。计算了４个不同参数齿轮的低阶固有频率和振型，并与实验结构数据比较，获得了满意的结果。这证明复约束法既极大地减少了计算所需内存和机时，又可建立包含轮齿和轮体在内的计算模型，保证了齿轮结构特征值问题计算的完整性和准确性。     

低雷诺数粘性流动的数值解

本文采用格林函数方法,使问题降低一个维度,计算了低雷诺数粘性流动,给出边界积分中几何系数的解析表示式,并讨论了计算区域选取对计算的影响。     

亚音速机翼有限基本解方法中的最佳展向控制点位置及诱导阻力计算

本文提出一种最佳展向控制点分布,用以提高有限基本解方法（method offinite elementary solutions,MFES）的精确度,同时降低计算结果对网格数目的依赖性,从而使网格数目大为减少。由于在诱导阻力计算中要求对升力分布计算得非常准确,所以诱阻计算是个非常困难的问题。在采用了最佳展向控制点分布后,我们顺利地解决了诱阻计算问题。当对侧滑机翼计算滚转力矩或对细长机翼计算侧缘吸力时,最佳展向控制点分布也呈现其高度的优越性。