全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

斜齿轮齿面柔度矩阵与修形的有限元计算

提出了一套斜齿轮修形用的有限元计算程序。该程序自动生成轮齿有限元网格和齿面接触线, 计算得到接触线上结点的柔度矩阵;然后考虑支承结构与齿轮的综合变形, 由数学规划获得最佳齿面修形。计算与实验结果对比证明该程序操作简便、计算精度高, 具有工程应用价值。      

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

航空用有机玻璃(PMMA)的非线性热粘弹-塑性本构关系

实验测定了航空用有机玻璃T=-40～40℃,ε_t=1.0×10_-4～1.0×10_-1s_-1范围内单轴压缩循环的应力-应变曲线。提出一种考虑温度和应变率效应的非线性粘弹性变形和塑性变形相耦合的本构关系,给出了完备的微分和积分形式的关系方程组,用于拟合各不同温度和应变率的实测曲线,获得了良好吻合的结果。     

一种新的有限元线性方程组的求解算法

本文采用链表数据结构存储有限元结构总体刚度矩阵,对一有限元结构分析系统进行了改进。算例表明,改进后的系统可使存储总体刚度矩阵所需内存及分解总刚阵时间分别下降到原系统的1/2～1/4,大大提高了计算机的解题能力和执行效率。该系统可在微机上进行大型复杂结构的受力分析。     

复杂前机身—进气道组合体的流场计算

本文提出了一种能够对复杂前机身-进气道组合体进行流场计算与气动分析的方法。通过求解Euler方程来模拟前机身-进气道组合体在亚、跨、超声速,不同攻角和发动机流量条件下的各类流场。     

一个非线性拟协调退化壳有限元

 根据非线性协调有限元列式方法,构造了几何非线性拟协调9结点四边形退化壳有限元。首先给出9结点四边形退化壳单元的几何描述和位移模式;接着讨论了非线性协调有限元列式方法;假设应力和应变;然后用所构造的壳单元,进行了方板、圆拱、球壳和圆柱壳的数值分析,并同精确解以及其他单元的解答进行了比较。结果表明本文壳单元是准确的和有效的。     

航空交叉流式湿热交换器的热性能分析

 本文分析研究两种流体互不混合的湿空气交叉流式热交换器的热性能。按热质平衡原理,首先建立了一套描述其热湿性能的微分方程组,而后应用有限差分技术把它们化为差分格式进行计算,最后与实验结果进行了比较,获得较为满意的结果。