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

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

机构磨损可靠性

首先简述了磨损的常规失效判据,然后提出了磨损可靠性的安全裕量方程,给出了可靠性计算公式。用算例阐明基本概念及方法思路。最后对防止粘着磨损咬卡提出了合理的设计准则。对极限磨损量的代表性作了说明。     

用二维N—S方程计算汽车流场的地面效应

采用有限体积方法求解了Navier-Stokes方程,计算了绕二维类汽车物体的绕流并研究了地面效应。在非正交网格系统上以强守恒形式构造出有限体积方程,采用压力校正算法来求解所得到的方程。对于湍流情况,采用了双方程模型的κ_ε模型。 为了验证采用本程序计算移动壁面效应的正确性,本文给出了Couette流动的计算结果。 本文也给出了汽车绕流的计算结果,它包括层流情况,从汽车头部开始由层流转为湍流的混合流动情况及湍流情况。对于每一种情况都比较了移动壁面的计算结果和固定壁面的计算结果。     

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

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

用全位势方程计算机翼的亚声速,跨声速和超声速绕流

对大后掠小展弦比细长机翼,本文对机翼纵轴垂直的每一横流截面生成O型网格,形成对机翼流场的H-O型网格,用守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕机翼的可压缩位流。自由流可从亚声速直到低超声速的全部跨声速范围。本算法要求机翼前缘有大后掠角,后缘可稍许后掠或前掠。本文算例表明,所研制的计算程序已可提供工程实用。     

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

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

三维C-H网格下跨音压气机转子Euler流场数值模拟

介绍了三维C-H网格的生成和全三维转子流场的数值解。以法国宇航院设计的一个高负荷跨音转子为算例, 并与他们的计算结果作了比较。应用这套数值模拟技术得以计算带进口预旋的某跨音压气机转子流场, 并对进口预旋对其下游的转子流场的影响作了数值模拟, 得出了若干有意义的结论。      

二维可压湍流边界层改进的积分方法

采用含高阶项的边界层动量积分方程,同时对常规边界层卷吸方程及延迟方程进行相应的高阶影响修正,得到改进的边界层积分方程组。应用此方程组对多种高亚音速及跨音速翼型边界层流动作了计算,并与一阶Ｇｒｅｅｎ方法计算结果及实验结果作了比较。结果表明,在边界层积分方程中保留法向压力梯度项及雷诺法向应力项,明显改进了翼面边界层接近分离区域处参数的计算精度。     

跨声速开槽壁超临界翼型洞壁干扰的N－S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一步，其结果就符合得更好，这为下一步用Ｎ－Ｓ方程进行洞壁干扰修正打下了基础。     

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

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