用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

刚度导数法在含裂纹铆接加筋板应力强度因子计算中的应用

本文探讨了刚度导数法用于含裂纹铆接加筋板应力强度因子计算的问题。提出了选取对刚度导数矩阵有贡献元素的一种新方法,成功地解决了裂尖位于筋条轴线附近时,采用刚度导数法计算应力强度因子的困难。算例结果表明,本文提出的方法是正确有效的。文中还阐述了刚度方程降阶的方法,这种方法对于应用小容量计算机求解较大的结构分析问题,有一定的参考价值。     

有限平板绕流Monte—Carlo方法仿真

本文采用直接仿真Monte-Carlo方法求解有限长平板绕流等问题。本方法是通过计算机跟踪仿真分子的运动来实现数值模拟的。仿真分子间的碰撞计算由统计抽样确定。碰撞模型分别选用了硬球分子模型和负幂律分子模型。仿真分子与固壁作用采用由完全扩散反射和镜面反射按比例混合组成的模型。为了检验方法的可靠性,还计算了激波结构和Rayleigh问题等一维流动,二维计算采用同步并行程序。数值结果表明直接仿真Monte-Carlo方法能够较好地模拟稀薄气体力学中的一些问题。对于二维计算,所花费的机时和所需的内存均在国内机器所允许的范围内。     

应用等效系统法评价飞行品质

本文介绍等效系统方法在评价飞行品质中的应用,包括等效系统方法的原理、所采用的数学模型和进行等效系统拟配的数学方法。最后,介绍了具体拟配过程中应当考虑的若干问题。     

计算旋转薄壳自由振动的传递矩阵法

本文介绍了计算旋转薄壳自由振动的一种传递矩阵法。以计算截锥壳振动的传递矩阵公式为基础,通过适当的变换,进一步对一般旋转薄壳结构的固有频率及振型进行计算。     

叶片盘的旋转失速振动响应分析

稳定气动力状态下的叶片盘的模型由连续系统组成。利用弹性和波动理论以获得叶片盘响应。响应是带有波动模态阻尼的模态解。由模态解,可以得到共振相对转速的表示方法和解释试验中的某些现象,此外,讨论了在气流稳定和失速状态下的广义共振情况。最后得出的结论是:在均匀稳定气流中不出现波动响应,但在失速流动下,波动响应呈随机的波动振动,而此时不产生驻波共振。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

具有陶瓷涂层的火焰筒壁温和热流计算

对具有陶瓷隔热涂层的气膜冷却式火焰筒壁面温度和热流提出~种计算方法。在分别建立陶瓷涂层和金属壁面热平衡方程的基础上,给出两者之间的导热耦合关系,使计算模型更加符合实际传热过程。计算过程中始终考虑到金属壁面的轴向导热。通过算例研究了陶瓷涂层对气膜冷却式火焰筒壁面温度和热流的影响。有涂层时壁面温度低于无涂层时壁面温度,但它们的变化规律相似。燃气温度越高陶瓷涂层的隔热效果越好。     

多跨加筋板的动态裂纹扩展

提出了一种加筋板数值分析中筋条与板相互作用的计算模型,用有限差分法对多跨加筋板结构进行了动态裂纹扩展分析计算。当筋条位于板的边界时,筋条的作用可以作为边界力处理;当筋条位于板的内部时,其作用可作为局部作用的类似体积力进行分析。按照这种模型对多跨加筋裂纹板进行静态与动态断裂分析。计算结果表明:筋条具有降低裂纹尖端应力强度因子作用,在裂纹动态扩展情形还能降低裂纹扩展速度甚至止裂。     

EBE技术在结构分析中的应用——EBE—CG方法及其并行实现

共轭梯度法(CG法)是求解大型稀疏有限元方程组的有效迭代解法。针对大规模结构分析问题,提出了不形成总刚度矩阵时的EBE-CG法及其并行处理方法。在国内规模最大的大型MIMD分布式系统上的初步数值试验结果表明,并行的EBE-CG法是很有效的。