一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

混合边界子结构的航空发动机机匣模型修正

提出了一种基于混合边界模态综合的复杂结构有限元模型修正方法。其主要步骤包括:①子结构划分,根据结构形式划分待修正区域，得到子结构和残余结构；②缩聚和装配,利用混合边界模态综合法，将子结构内部自由度集缩聚至混合边界自由度集，得到子结构缩聚矩阵，并与残余结构的系统矩阵进行装配；③修正,基于灵敏度分析方法，对装配后的残余结构进行参数修正。将该方法应用于航空发动机外机匣的精细化有限元建模及模型修正研究，针对局部连接结构参数修正，子结构模型修正方法的参数收敛后最大误差为064%，其计算效率提高了515倍。算例结果表明，该方法在保证修正精确度的同时，能提高大规模复杂结构有限元模型修正的计算效率。      

高速风洞洞壁干扰减小与修正研究

亚跨声速风洞试验的洞壁干扰问题是影响风洞试验结果准确度的一个重要因素。南京航空学院ＮＨ－１高速风洞首先使用了两种可变开闭比风洞壁来减小洞壁干扰，然后发展了一系列用壁压信息法对剩余洞壁干扰效应进行修正的方法，对国内外大量高速风洞实验数据进行了洞壁干扰修正，将修正结果和ＮＡＳＡ的非线性修正法结果。自修正风洞实验结果及无洞壁干扰的Ｎ－Ｓ方程计算结果进行了比较。结果表明，本文的修正方法结构正确，而计算量远     

静止叶栅三维粘性数值模拟

应用Ｎ－ Ｓ方程压力修正算法和ｋ－ ε两方程湍流模型,对端弯叶片和基准（无端弯）叶片的平面叶栅进行了三维粘性数值模拟。在采用简单交错网格的情况下,通过对计算控制体通量的速度进行修正的方法解决了Ｎ－ Ｓ方程压力修法算中遇到的压力场和速度场呈锯齿状分布的问题。用实验测量值对本数值模拟方法的精度和可靠性进行了检验     

基于统计分析和响应面方法的有限元模型修正研究

通过统计学中的方差分析对待修正参数进行筛选,利用对修正参数的方差分析结果,选择对响应特征有显著影响的参数作为修正对象。再利用试验设计方法,通过有限元计算,获得一系列试验点的响应特征,拟合出响应特征关于修正参数的响应面模型。模型修正时,用响应面模型替代有限元模型进行迭代计算。由于只需要在试验点的响应特征计算时调用有限元计算程序,对于大型工程结构,可极大的提高修正效率。     

液体火箭发动机喷管简化建模及模型修正

为了准确掌握喷管的动力学特性,提出了一种等效简化建模结合试验数据的模型修正方法。首先,将喷管的原始几何模型经过几何处理成简化模型,建立了喷管的有限元模型,对有限元模型在自由条件下进行模态计算,将得到的计算数据与测试数据对比分析,再利用测试数据对有限元模型的弹性模量参数进行修正,修正后的喷管有限元模型前9阶模态计算结果与测试结果频差在5%以内,MAC值(模态置信准则)在0.8以上。表明此方法是一种高效可行的喷管简化建模方法,既保证了精度又提高了计算效率,对其动力学特性分析、振动响应预测等方面具有重要应用价值,对于液体火箭发动机其他部件的动力学建模及分析也具有普适性。     

高度修正实现刚度等效的梁结构有限元建模方法

民用飞机中对于梁结构的有限元建模,一般采用刚度等效方法将梁建成杆单元与剪切板单元的组合形式,这种有限元简化方法仅能保证截面刚度与真实截面刚度近似等效,但有限元中的截面面积与实际面积却相差较大,致使有限元分析的结果不能准确反映结构严重部位的应力情况。针对传统刚度等效有限元建模方法的不足,提出一种采用高度修正的刚度等效有限元建模方法;将该方法的计算结果与传统建模方法的计算结果、理论计算结果、试验测量结果进行对比。结果表明：本文提出的高度修正的有限元建模方法比传统杆板系建模方法更准确,与试验结果更接近。将该方法用于民用飞机结构有限元建模,能够对真实结构进行更准确地模拟,分析结果也更精确,对民用飞机结构的强度分析和安全保障具有重要的意义。     

亚声速大迎角模型试验洞壁干扰修正方法研究

通过测量洞壁附近的压力分布来模拟透气壁试验段的洞壁边界条件，采用数值求解Ｅｕｌｅｒ方程的方法模拟模型在风洞中的绕流场，然后将洞壁边界条件用远场边界条件替换，用同样方法算出模型在自由流中的绕流场，从而计算出洞壁干扰对模型的绕流场和气动力的影响。针对大迎角模型的洞壁干扰问题的特殊性，采用重叠网格技术，对飞机大迎角标模和ＳＢ－０３模型进行了亚声速大迎角的洞壁干扰计算与修正。     

某航空发动机机匣的动力学模型修正

利用动力学模型修正技术对某航空发动机机匣的有限元模型进行了修正.通过振动模态测试得到了实际机匣的模态数据用以作为有限元模型 修正的基准.利用频率对单元刚度的灵敏度分析选定了修正区域.在此基础上,应用1阶优化方法对机匣的有限元模型进行修正.研究结果表明:修正后机匣有限元模型的前10阶模态的计算值与实际测试的误差都在29%以内,可以应用在后续的发动机整机动力学分析等方面.      

计算流体力学中有限元的强间断处理

本文首次提出了无波动，无自由参数的有限元离散方法，它是基于矢通量分裂后其物理传播特性，利用差分格式构造过程中对激波上下游不同性质的分析，对一维Ｅｕｌｅｒ方程进行了仔细分析，并将其推广到二维，轴对称的无粘，有粘方中，从计算实践看，这种方法简便、实用，对激波，分离，激波与边界层干扰的复杂流场有较好的分辨能力和计算稳定性，比以往的迎风有限元，加入Ｌａｐｄｉｕｓ人工粘性修正的有限元有更好的流场分辨能力，使     

Navier-Stokes方程积分表达式的有限元解法

本文采用定常不可压Navier-Stokes（NS）方程的积分表达式方法计算粘性流动问题,用线性三角形单元划分计算区域,建立了适当的迭代求解过程。本文改进了压力计算方法,用全变量松驰加快了收敛速度,用矩阵分块计算减小了内存限制,并可以通过对矩阵特性的分析更好地控制迭代过程的稳定性。     

高速发动机油膜惯性对活塞裙润滑影响

利用一种新的迭代算法 ,分析了高速发动机油膜惯性对活塞裙润滑特性的影响 .该迭代方法是利用有限元法和差分法交替求解Navier Stocks方程和雷诺方程 ,据此导出了适用于高速发动机活塞裙的润滑计算的混和润滑雷诺模型 .新的模型借助惯性系数 ,引入了油膜惯性项 ;同时给出求解含有油膜惯性项的迭代步骤和有限元表达式。计算结果表明 :随着惯性系数和活塞裙的长径比的不同 ,油膜惯性会对油膜的摩擦力、压力和承载力产生不同的影响 ,这种影响对承载力尤为明显 .该润滑模型也可用于中、低速发动机的活塞裙润滑计算以及不计入惯性项 (惯性系数置为零 )的某些润滑问题求解 .     

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

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

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

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

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

钝头飞行器大攻角绕流流场计算

本文用有限差分法求解Ｎ－Ｓ方程，对具有大面积底部飞行器的高超声速大攻角绕流流场进行了计算，通过积分物面压力给出了飞行器的气动力特性，并研究了攻角的影响，物面压力系数和气动力系数与实验数据进行了比较，符合较好，结果可应用到返回舱和弹头的再入气动特性研究。     

线性GMRES算法在求解全速位方程中的应用

从全速位方程出发,利用压强极小积分得到有限元方程组。采用适用于有限元法的块结构网格生成技术,给出求解有限元总体方程组的线性GMRES算法。该算法比线松弛迭代法的收敛速度快得多。对M6机翼跨音速有升力情况进行了计算,其结果与实验值符合较好。     

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果     

粘弹性阻尼结构的动特性分析

本文研究了粘弹性阻尼结构动力特性的分析方法。根据分数导数模型表示的粘弹性材料本构方程,对一种阻尼材料用拟合法构造了它的力学性能解析表达式。建立了一种八自由度夹层梁有限元模型,从而给出了粘弹性阻尼结构的有限元动力学方程。为分析粘弹性阻尼结构的动力特性,提出了一种非线性特征值问题的迭代解法。应用上述方法对悬臂约束阻尼层梁的复模态参数进行了分析,并进行了实验验证。结果表明,分析与实验的结果一致,说明本文所提出的方法是可靠的且有较好的精度和实用性。     

绕升力机身跨音速流动的有限元素法数值分析

 <正> 本文应用全速位方程最小压强积分的有限元素法解绕升力机身的跨音速流动。机身头部可以是尖的或带有进气道的。用头部伸出无穷长圆柱来模拟头部进气。尖头机身在元素形状处理上具有一定的复杂性,从数值实验角度考虑我们分别按尖头和头部伸出无穷长细圆柱来近似尖头头部两种方法对尖头机身进行计算。采用人工密度法捕获激波,但由于绕机身流线形状比较复杂,直接采用机翼问题中的人工粘性公式,解往往不收敛。我们推导出一种较精确的,适用于复杂流动情况的人工粘性计算公式,为了加速收敛,应用网格逐次加密技术,可使收敛速度提高2～3倍。本文还给出绕升力机身跨音速流动的远场速位解析式。