机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

扩散抛物化(DP)NS方程组的意义及其在计算流体力学中的应用

本文首先讨论扩散抛物化（DP）NS方程组的早期研究工作：它的提出、数学性质、意义和在CFD中的应用，然后讨论扩散抛物化理论的一些新发展。这些新发展是对NS方程组数值计算进行物理分析的基础上得到的，其中包括NS方程组差分计算时，粘性剪切流对网格间距和格式精度的要求；粘性项只保留剪切项的广义扩散抛物化（GDP）NS方程组，它的性质和应用。由于高Re数流动在NS方程组的差分计算中，网格Re数彼此相差悬殊的特点，产生了计算离散单元守恒方程组的新的算法思路，即离散流体力学（DFD）算法。在DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD格式的构造、它的优点和应用。并以超声速绕前后台阶流动为算例，来说明GDPNS方程组的用处和DFD算法的优点。DPNS方程组、GDPNS方程组、DFD算法是高智提出的，对这些问题，他的合作者从理论、算法、数值验证和某些应用又作了系统的研究。     

航空发动机状态监控亚定方程组有效数字算法

缺少状态监控的信息量是基于数学模型的航空发动机状态监控的主要困难。航空发动机状态监控亚定方程组缺少成熟可靠的求解算法。提出一种求解航空发动机状态监控亚定方程组的有效数字算法。此算法按实际工程应用所需的有效数字位数将解空间离散化,再加上亚定方程组本身的物理背景约束,然后用求解整数规划问题的方法将亚定方程组求解出来。仿真结果表明此算法是有效的,且2或3位有效数字的仿真结果曲线拟合较好,能较好地满足实际工程应用的需要。     

基于Broyden算法的航空发动机气路故障诊断

针对基于Kalman的故障诊断算法响应速度慢、多故障诊断及非设计点诊断精度低的问题,提出一种基于改进Broyden算法求解方程组的航空发动机气路故障诊断方法。针对涡轴发动机,以模型输出跟踪发动机输出为准则确定3个方程,结合发动机模型中的2个平衡方程,构建气路故障诊断方程组,通过改进Broyden算法求解方程组以获得部件性能退化因子及模型猜值。数字仿真结果表明,所提出的基于Broyden算法求解方程组的航空发动机气路故障诊断方法,在包线内的单故障和多故障诊断稳态误差均小于0.35%,且诊断过程算法单步运行最大耗时小于2ms,具有良好的实时性,远优于Kalman滤波方法,验证了算法的先进性。     

求解非线性方程及方程组的模拟退火算法

探讨了用模拟退火算法求解非线性方程及方程组的方法,该方法把方程组求解问题转化为函数优化问题.计算中不需要使用目标函数的导数信息,在Matlab环境下实现了该算法,数值实验结果表明了该算法的有效性.     

一种改进模拟退火算法在非线性方程组求解中的应用

将非线性方程组的求解问题转化为求多元连续函数极小值的无约束优化问题,提出一种求解该优化问题的改进模拟退火算法,从而得到非线性方程组问题的高精度解,数值仿真结果表明了该算法是可行而有效的.     

二维可压缩Euler方程组的SUPG有限元数值解

本文在文献[1]的基础上,构造了二维非定常可压缩Euler方程组的SUPG变分方程组。文中对由文献[1]提出的预估——校正算法作了分析,并用该算法对变分方程组进行解算,所做算例分别为:无粘性激波在固壁上的反射、零攻角10％抛物翼型亚、跨声速绕流、带攻角NACA 0012翼型跨声速绕流,取得了较好的数值结果。     

含损伤粘弹性本构及其在有限元分析中的实现

为研究固体发动机药柱的损伤粘弹性行为,针对一种耦合损伤的三维非线性粘弹性本构模型,采用积分算法对其进行了数值离散,导出了本构方程的增量有限元形式.针对该增量型本构方程对应的非线性方程组,建立了采用Newton-Raphson迭代法求解的应力更新方法,并基于Abaqus软件和Fortran语言开发了相应材料本构子程序,最后应用算例对所设计数值方法及程序进行了验证与分析.结果表明,本文有限元解与对应解析解吻合良好,保持了很高的精度,所设计方法及程序有效,可应用于实际型号发动机药柱损伤结构分析.     

螺旋锥齿轮齿面接触分析改进算法研究

针对现有齿面接触分析(TCA)的不足,提出了螺旋锥齿轮齿面接触分析的改进算法.基于空间曲面啮合原理、Enler-Rodrigues公式以及接触区几何形状的等量关系对齿面接触分析算法进行了改进,通过方程变换,减少了求解接触点位置的非线性方程组维数,由啮合方程求得到初值选取计算公式,使初始值求解有规律可循.以一对弧齿锥齿轮为例,对算法进行了对比分析与验证,结果表明:改进的齿面接触分析算法可减少求解非线性方程组时初值选取的繁复性与增强非线性方程组迭代求解的稳定性;并利用几何等量关系求解接触区几何参数,能避免求解复杂的齿面相对曲率等二阶几何参数.      

基于粒子群算法的发动机部件模型求解

传统求解航空发动机部件模型中的非线性方程和共同工作方程组的方法因对初值的依赖使得模型并不总能收敛。首先针对模型中的非线性方程求解的特点,在带邻域的粒子群算法的基础上,将收敛因子、被动聚集压力因子和自适应惯性权重引入算法,提出一种混合粒子群算法HPSO1。而对模型中的共同工作方程组的求解,则是在基本粒子群算法的基础上,将收敛因子和被动聚集压力因子引入基本粒子群算法,提出另一种混合粒子群算法HPSO2。实验仿真表明HPSO1和HPSO2均克服了对初值的依赖性,因而对发动机部件模型求解是有效的。     

有限元的一种并行算法及实现

当前,由于有限元素法大都使用串行算法,处理时间较长。本文根据有限元素法的特点和采用的并行计算机结构,提出了一种并行有限元算法,能大大提高其处理速度,加快处理周期。且算法简单,实现方便。     

一种无条件稳定的时域推进法与有限元混合分析瞬态热传导问题

本文将一种新的无条件稳定的时域推进法与有限元结合,用于分析瞬态热传导问题。根据微-积分型热传导控制方程,对时间变量在小区闻内插值积分,变初-边值问题为一系列离散时刻的边界值问题,再应用有限元法求解之。由于推进中的每一时刻的解,都严格满足原给定的初始条件,就消除了累积误差的影响。此外,这种时域推进法对时间变量的数值积分也优于数值微分。因而,与现行的基于微分型控制方程,采用时-空有限元模型同时离散,再逐步求解的直接积分算法相比,计算精度可大大提高。不难预料,当求解较长时间后的瞬态值时,本算法的优越性会更加明显。     

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

无粘跨声速绕流的有限元法和界积分方程的某些进展

本文对无粘跨声速绕流的有限元法和边界积分方程法的目前进展,特别是对提高跨声速流算法效率起着主要作用的成果,例如普适元素和块结构有限元网格生成技术有限元解法的多重网格技术、SUPG有限元法、GRMES算法和可压流的完全边界积分方程法等方面作了回顾和评介。     

结构有限元模型的试验建模

 对于复杂结构系统,用试验测取的模态参数来识别或修改系统的物理参数,以构造其动态数学模型。这种方法是建立振动系统数学模型的逆方法。 近些年来,国内外不少学者在这一领域做了很多工作。他们方法的共同点是识别系统的总质量矩阵和总刚度矩阵,因此未知量较多。考虑到结构系统总矩阵元素之间有一定关系,文献[3]讨论了一种建立动态有限元模型的逆方法。     

各向异性Chaboche黏塑性本构方程隐式应力积分算法

针对各向异性Chaboche黏塑性本构方程,利用有效参数控制法推导了该方程的隐式积分算法,该算法把所有未知量都用同一个控制参数表示,并形成控制方程,从而把本构方程应力积分的问题转化成一个非线性控制方程的求解问题.该积分算法程序,与通用有限元软件Marc结合,通过单个单元对单向拉伸、循环、松弛的响应的实验数据进行了数值模拟,通过与显式算法对比,该隐式算法的稳定性和收敛性有所提高.构件的对比分析也表明隐式算法比显式算法具有较高的效率.      

Navier-Stokes方程两层稳定有限元算法分析

分析了定常Navier-Stokes方程的两种两层稳定有限元算法。它们将局部Gauss积分稳定化技术和两层算法的思想充分结合,采用低次等阶有限元P1-P1或Q1-Q1对N-S方程进行数值求解。误差分析和数值算例结果表明,当粗、细网格尺度H=O（h1/2）时,它们与在细网格上的单层有限元算法具有相同的收敛速度,而两层算法却节省了大量的计算时间。相比之下,Simple算法具有更高的计算效率。而且进一步发现Oseen算法能够对小粘性系数N-S方程进行有效求解。     

考虑预应力和液体影响的圆柱壳颤振分析

 为了研究预应力和液体对圆柱壳颤振的影响,采用杂交元方法,建立了圆柱壳的气动弹性方程。由Sanders薄壳理论和经典有限元理论,从壳的精确解推导节点位移函数入手,并由一阶活塞理论得到结构的气动阻尼矩阵和刚度矩阵,最终推导出考虑预应力和内部液体影响的圆柱壳的混合有限元公式和气动弹性方程。通过特征值法验证了此种有限元方法的正确性,并重点研究了预应力和内部注有液体对圆柱壳稳定性的影响。计算结果表明:预应力和液体对圆柱壳的颤振特性有显著的影响。     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。