用有限元法分析接触问题的一种新的力学模型及计算方法

阐明了用有限元法计算物体接触问题的一种新的力学模型:主动接触体网格中一个节点与被动接触体单元中任意点相接触的模型,放弃了过去常用的节点对模型,为网格划分带来了便利。在这种模型下,给出了粘式、摩擦滑移、开式等接触边界条件。按照修正势能原理导出了对称的接触控制方程,可求解接触体具有非线性的静态/动态的接触问题。本文给出了在一次求解迭代后判断各节点新的接触状态的方法,可进一步决定是否和如何修改接触控制方程。      

宇宙线福布什下降的数值模拟

本文从银河宇宙线的太阳调制方程出发,认为激波对宇宙线的影响是由一扰动区产生的。在此扰动区中太阳风速度增加,扩散系数下降。由此进行了数值模拟,模拟结果表明:扩散系数的下降在产生福布什下降时要比太阳风速度的增加更为有效;福布什下降在近日球层内向外传播时,其幅度随径向距离的增大而衰减;两个无相互作用的激波同时存在时所产生的福布什下降为每个激波单独存在时的福布什下降的简单线性迭加。      

偏置动量系统中的挠性天线指向控制研究

从工程实践角度,探讨了在目前广泛采用的偏置动量轮三轴姿态稳定卫星上,进行大尺寸（３～５ｍ）和高指向精度（０．０５°）的挠性天性指向控制问题。在建立具有挠性天线的星体动力学模型,以及在借鉴和改进有关抑制挠性天线振动研究的基础上,设计了三轴稳定卫星平台上的挠性天线指向控制系统的方案。最后,根据数值仿真结果,给出了影响挠性天线指向精度的主要因素。     

太阳电子事件粒子源区的位置

对ISEE-3的引个太阳电子事件资料,进行了系统地分析。通过Saha平衡模型分别计算了高Fe事件和贫Fe事件源区的平衡温度大小,指出两类事件粒子源区不相重合,高Fe事件的粒子源区位于高色球层,贫Fe事件粒子源区在色球层-日冕之间的过渡区。对高Fe事件和贫Fe事件的形成作了初步的设想。      

模拟高度复杂流动的自适应网格算法

提出一种基于对等弧分布律进行扭转修正的自适应网格算法。该方法将三维网格系统分裂为３个单方向自适应网格子系统,在计算中与主控方程的求解格式相耦合,交替地使用。应用所构造的算法结合ＭａｃＣｏｒｍａｃｋ显式格式对喷流／绕流干扰流场进行了雷诺平均Ｎ－Ｓ方程模拟,给出了清晰的剪切层与激波波系结构。计算结果表明,该自适应网格算法可较大程度地提高ＭａｃＣｏｒｍａｃｋ格式的分辨率。     

平均应变效应的循环粘塑性描述

提出了一个考虑具有平均应变循环变形行为的粘塑性本构模型,在该模型中,引入了可以计及平均应变效应的应变幅值记忆参数,构造了相应的背应力演化方程,定义了恰当的非比例度;还提出了对先前加载历史循环变形行为衰减记忆的各向同性变形阻力演化方程。模型用于１Ｃｒ１８Ｎｉ９Ｔｉ不锈钢的平均应变效应的本构描述,结果表明模型合理。     

N-S方程及其简化形式在圆柱绕流分离点邻域的数学性质

证明了在圆柱绕流分离点的邻域上 ,N-S方程和高智的简化 N-S方程是正则的 ,即存在级数解 ;而 Davis和 ГКП的简化 N-S与经典的边界层方程一样存在 Goldstein奇性     

翼-身组合体绕流的Euler方程数值模拟

用代数方法生成翼-身组合体H-O型网格,并用有限体积法研制出翼-身组合体绕流三维Euler方程计算程序。该方法的特点是改进了机身与机翼表面网格点分布,机翼后缘有后掠时也能保证后缘与网格线一致。程序除能按常规提供横流截面展向压强分布外,还能提供弦向压强分布。对NASA TND-712的翼-身组合体模型的计算结果与实验符合很好。     

飞机电网零序阻抗的数值计算

从麦克斯韦方程出发,推导了飞机电网零序阻抗参数的计算公式。对几种不同敷设方式载流导线的零序阻抗进行了计算,其结果与标准值较为接近。这种数值计算方法比传统的解析近似估算方法更为合理和准确。     

基于非交错有限体积格式和 k-ω 湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。