高速列车绕流的数值模拟

给出了用低速N－S方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Train3D，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

基于ANSYS的薄壁件反重力铸造充填形态的数值模拟

介绍了CFD的基本原理以及有限体积法的概念,采用ANSYS 6.0进行网格划分和控制方程的计算,模拟了不同壁厚、不同充型速度下薄板铸件的反重力充型形态,得出了壁厚,充型速度对充型形态的影响规律。     

尖头细长体中小迎角流动结构的数值模拟和实验研究

利用可压缩层流Navier-Stokes方程模拟了尖头细长体中小迎角的流动结构。给出了4个迎角状态的物面流谱,分析了极限流线随迎角的发展过程;给出了25°迎角的横截面流谱,分析了它们的拓扑特征。展示了由主涡涡对、二次涡对和Tertiary涡对等三重涡组成的完整的涡系结构,强调了Tertiary涡在涡系演化中的意义,及其沿轴向发展过程中迹线的合并与分叉现象。计算与实验结果定性一致。     

激光扫描路径对直接金属激光烧结温度场的影响

利用有限元分析软件AN SY S,对激光烧结温度场进行了数值模拟,建立了激光烧结薄板的温度模型;运用APDL语言控制热源的热流密度、移动速度以及扫描轨迹;研究了激光行程对烧结温度场的影响,扫描端点温度场的不对称及较大的温度梯度造成了端点球化现象;激光扫描线间的耦合作用使端点球化现象随扫描线的增加而逐渐显著。根据该程序得出的模拟结果和实验结果基本一致。     

粘性流体大幅晃动的ALE有限元模拟

采用有限元法数值求解了具有自由液面大幅移动边界的Navier－Stokes方程。对流体区域采用了任意的拉洛朗一欧拉（ALE）运动学描述，网格结点可以任意移动而不依赖于流体的运动，结合了拉格朗目描述易于处理移动边界与区域变形的优点和欧拉描述可克服单元缠结的优点，提出了简单合理的网格更新方法，能精确地跟踪运动的自由液面。为了更精确地处理强对流项的作用，采用了迎风流线Petrov－Galerkin（SUPG）加权余量法建立有限元方程。算例表明本方法对贮箱内流体的非稳态大幅晃动过程的数值模拟非常成功。     

尾支杆干扰非线性修正方法的初步研究

利用CFD方法 ,采用两套网格重叠的形式 ,研究了一种对尾支杆干扰进行修正的非线性方法。运用本方法 ,研究了尾支杆等直段长度、直径和锥角对模型气动特性的影响 ,研究了四个模型有攻角时的尾支杆干扰 ,并对GBM 0 4A模型的零阻进行了尾支杆干扰的修正。     

惩罚函数法在发动机性能计算中的应用

按照涡喷发动机工作时遵守的气动热力学定律，建立了发动机的稳定数学模型，并针对双轴涡轮发动机性能计算的非线性数学模型，利用优化和最优控制理论中的惩罚函数法，对发动机性能进行了基于约束的最优求解程序设计，利用程序对某型发动机的稳定飞行特性进行了计算，计算结果表明，使用惩罚函数法求出的特性曲线变化规律正确，发动机平衡方程的解的精度更高，在此基础上得出的发动机飞行特性也更准确。另外惩罚函数法对初值的要求很低，迭代运算不易发散。     

钛合金高压气瓶表面裂纹断裂试验及断裂应力预计

本文扼要地介绍气瓶基体金属疲劳表面裂纹断裂试验的结果,并与几种常用工程计算方法的计算值进行比较。对韧带未完全屈服的断裂情况,我们导出了一个适合计算表面裂纹断裂应力的公式,在此基础上还推荐了一种预计表面裂纹断裂应力的简便方法。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

应力极值分布、强度威布尔分布的可靠性

本文运用应力—强度干涉理论,推导了应力为Ⅰ型极小值分布,强度为威布尔分布的可靠度计算公式,并对冗长的计算公式进行简化,在简化公式的基础上,运用一定的数学技巧,改变积分公式中的积分变量和上下限。将被积函数化成在某一区域内的可积函数。采用de Boor编制的一种严谨的自适应Romberg外推格式的FORTRAN程序进行数值积分。对应予不同的组合参数,给出应力服从Ⅰ型极小值分布,强度服从威布尔分布的可靠度数值。本文最后讨论了服从这两种分布的组合参数的变化对可算度数值变化的影响。