褶皱结构成形时的形状偏差计算

解决了褶皱结构成形时相邻的平行四边形平面元素之间连接区棱的修圆半径值预测问题。把毛坯模拟成薄板，定义了平行四边形平面元素的边界条件，采用能量法导出挠度和力函数的非线性微分方程组，利用积分－差分法实现其数值解，并提供计算结果和分析应力－应变场。采用的计算方法可用于研究毛坯的刚度参数和工艺参数对褶皱结构的形状影响。     

轴对称锻件成形过程的热力耦合有限元分析技术的研究

锻件成形过程中,非稳态、不均匀的温度场对金属的塑性流动有很大的影响，尤其是高温成形过程。本文对轴对称段件成形过程的热力耦合有限元分析技术进行了研究。论述了刚塑性有限元分析方法，建立了热力耦合分析模型，开发了轴对称锻件成形过程的热力耦合有限元分析软件。通过将圆柱体镦粗过程的模拟结果与有关文献中的实验结果和模拟结果的比较，说明了该方法和软件的正确性。     

有进动的高速旋转体的应力分析

本文讨论了采用有限元法对具有进动的高速旋转的轴对称体进行应力分析的方法，首先，推导了轴对称体中哥氏惯性力的分布公式，在柱坐标系下，哥氏力中含有沿径向的零阶对称不载荷和沿轴向的一阶对称及一载荷，首次，采用了具有计算精度，对边界适应性强的八节点四边等参数环形单元，并将受哥氏惯性力的轴对称体看作非轴对称问题，对有限元分析中位移模式的建立和哥氏惯性力的处理作了讨论。最后，简要地说明了所研制的计算机程序，并     

板料应力成形极限判据的有限元分析程序

为了解决在板料成形有限元数值模拟中以依赖于加载路径的应变成形极限为判据所存在的问题,以板料塑性变形时所遵循的应力应变转换关系、H ill二次厚向异性准则为计算模型,基于有限元商用分析软件DYNAFORM,开发了用于板料塑性成形工艺的以应力成形极限为判据的有限元分析程序。利用该程序分析了两种拉深件的有限元数值模拟结果,实现了用应力成形极限图作为判据对板料成形极限的预测分析。所开发的以应力成形极限作为判据的有限元分析程序为分析多道次板成形的成形极限提供了新的手段。     

基于细以方法的颗粒增强复合材料弹塑性分析

针对高体积含量的颗粒增强复合材料，提出了一种细观结构模型：将颗粒简化为同质、同尺寸的弹性圆球，两颗粒之间的连接基本简化为一弹塑性短圆柱体，并假设细观应力、应变和塑性区均为轴对称分布。其体和颗粒的变形行为分别简化为类似于弹性地基和弹性半空是，以此为基础，建立了反映一对颗粒法向和切向变形协调关系的两个积分方程。数值求解这两个方程，可得到一对颗粒间基体中细观应力的分布形式，从而建立颗粒对的细观弹塑性本构，采用平均化方法，进一步推导出材料宏观的应力－应变关系。本文利用该模型对一种金属基复合材料的拉伸实验进行了模拟，理论预测与实验结果吻合较。     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

随机材料板结构的传递函数分析方法

本文将二维随机材料场作Karhunen-Loeve正交展开，应用传统有限条理论和摄动随机分析思想建立起静态随机板结构－系列微分方程，然后用传递函数方法求解这些微分方程得到结构的随机响应，该方法是一种半解析数值计算方法。     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

基于数值模拟的超塑成形压力—时间曲线计算机辅助设计

确定压力随时间变化的关系曲线对于超塑成形工艺设计非常重要,本文在超塑成形压力计算时引入能量消耗控制策略,提出了一种压力－时间（Ｐ－Ｔ）曲线设计方法即恒能耗率法,实现了基于有际元数值模拟的超塑成形Ｐ－Ｔ曲线计算机辅助设计,文中对有关的压力控制算法及算式进行了较详细的介绍和推导,分别给出了半球壳自由胀形和圆筒形凹模约束成形过程的Ｐ－Ｔ曲线设计计算实例,并对计算结果进行了讨论,本文的研究拓宽了刚粘塑性有     

结构固有振动特征值问题矩阵摄动法的一个拓广

本文探讨振动结构的物理、几何参数与固有模态之间的关系。处理这类问题的常见方法是特征值问题摄动法。但该方法得到的解只有当参数作小扰动时才有效。为克服这一缺陷,本文将前人已得到的模态关于参数的导数表达式看成是一组常微分方程式,从而运用常微分方程的定性研究方法和数值计算方法来研究这类问题。文中用这一拓广的方法探讨了模态在参数空间上的连续性及重频点对模态演变的影响,并对一个轮盘叶片模型进行了数值计算和实验分析。     

Roe格式在多组元燃烧流场数值模拟中的应用

对Roe通量分裂格式进行了详细讨论,并将其推广应用到基于非结构网格的高超声速多组元燃烧流场的数值模拟中。控制方程采用三维多组元Euler方程,重点针对Roe格式开展研究;考虑到燃烧流场求解过程中出现的刚性问题,对化学反应源项采用点隐式方法处理。对超声速条件下钝头体激波诱导燃烧流场进行了数值模拟,比较分析了3种经典化学动力学模型对流场结构的影响,获得与数值模拟、试验数据相符的结果。结果表明发展的方法可以用于多组元燃烧流场的数值分析,且化学动力学模型的选取直接影响激波诱导燃烧流场中诱导区厚度。     

含化学反应流动的TVD格式及其在燃气射流中的应用

以Harten标准TVD格式为基础，结合固体火箭燃气射流的特点，以数学方法系统地推导了适合于高温、高压和高速流体流动的数值格式，给出了压力偏导数的合理计算公式；利用化学动力学知识，对火箭燃气射流流场中存在的多组分、含有限速率的化学反应问题进行了论述，阐明了化学反应质量源项的求解方法。以12组分9反应方程模型为例，利用编制的计算程序，对某火箭燃气自由射流流场进行了模拟。通过对结果的分析，肯定了数值格式的正确性。     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。     

二维突扩管流的有限分析解

运用有限分析法求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，获得了二维突扩管流的数值解。对于雷诺数为２５，５０，１００，１５０，２００的来流，预报的流动图谱具有相似性，回流区域及再附距离的计算也获得了比较满意的结果     

机载蒸发循环制冷系统板翅式蒸发器换热特性

从叉流板翅式蒸发器翅片传热和隔板热平衡方程出发建立了蒸发器的分布参数稳态数学模型,用有限差分法计算了蒸发器内的温度分布和换热特性。模型中制冷剂侧采用均相流动假设,空气侧考虑了湿工况的影响。为了求解两流程板翅式相变换热器的数学模型,本文提出了一种以第一、第二流程分界隔板上温度分布为收敛依据的计算方法,得出了蒸发器中空气温度场及制冷剂两相流干度和温度场的分布规律。计算结果表明,蒸发器中温度场和干度的分布规律是由叉流式换热器的结构特点及制冷剂和空气的相对流动方向决定的,制冷剂在一个流动通道中沿垂直流动方向干度分布很不均匀,即一边是过热蒸气,而另一边则含有未蒸发的制冷剂液体。由于过热蒸气与空气的换热系数较低,蒸发器中存在大面积过热区势必影响其换热效率,同时,蒸发器出口含有制冷剂液滴不仅减少了系统制冷量,降低了系统的性能系数,严重时可能对压缩机的工作造成不良影响。因此在设计蒸发器时必须充分考虑干度不均匀性的影响,采取必要措施尽量使出口干度均匀。计算结果与相关文献中实验结果比较证明了本文所建立的模型是合理的。     

鸭式布局战斗机非常规机动的流场机理数值分析

以先进战斗机的非常规机动为对象,发展了一种适用于大幅度运动变化的非结构嵌套网格生成方法,建立了一整套非定常流场N-S方程数值求解方法。在对三角翼动态气动特性计算验证的基础上,模拟了飞机过失速机动条件下飞行姿态和来流速度的变化特征,对鸭式布局战斗机"眼镜蛇机动"的非定常涡结构、非定常气动力效应和气动特性进行了数值研究,揭示了鸭式布局战斗机"眼镜蛇机动"的非定常流场机理。     

平面COUETTE流动的湍流斑研究

湍流斑是近壁剪切流动中的重要现象,它们的生成和发展与流动转捩和湍流的形成密切相关.采用Navier-Stokes方程直接数值模拟和壁面脉冲模型研究平面Cou-ette流动中湍流斑的生成和演化.高精度高分辨率的三维耦合紧致差分格式被用于数值计算.探讨了平面Couette流动中湍流斑的重要特征,包括高频脉动的形成、雷诺应力的产生、扰动的急剧增长和湍流斑形状的变化,特别是流向涡的复杂演化过程.结果表明层流流动中的湍流斑具有湍流流动中一些基本特征.     

Development of a Synchronization Method for Fluid-Thermal Study of Hypersonic Flow

A synchronization method is developed for the fluid-thermal study of hypersonic flow.Different from conventional loosely/tightly coupled methods which separately deal with the flow field and the structure temperature field,the presented method expresses the governing equations in a unified framework so that the two fields can be calculated simultaneously.For efficiently solving the unified equations,the finite volume method together with the dual-time stepping approach is employed.Like in the flow field,the local time step is also used in the temperature field,which is determined from thermal conductivity spectral radii.In order to treat the fluid-structure interface more conveniently,an expanded virtual boundary is introduced.For validation,several fluid-thermal hypersonic flow problems are simulated.The computed results are compared with those obtained from the coupled methods and the experiment.In the continuous heating problems,the stagnation temperatures predicted by both the coupled and synchronization methods are in good agreements with the experimental data.In the unsteady flowthermal hypersonic flows,the stagnation heat fluxes predicted by the presented method and tightly coupled method are basically the same,which agree better with the experimental data than those predicted by the loosely coupled method.In terms of prediction of the stagnation temperature,the synchronization method shows better accuracy than the tightly coupled method.     

一种解Dirichlet问题的虚拟区域法及其在广义Stokes问题中的应用(英文)

讨论了求解线性 Dirichlet问题的虚拟区域法及其在广义 Stokes问题中的应用 .这种方法通过 Lagrange乘子技术来处理 Dirichlet边界条件 ,因而非常适用于粘性流动问题中的无滑移边界条件。为提高解的精度 ,我们根据事后误差估算对网格进行自适应加密。Mini-element离散被应用于广义 Stokes问题的求解。最后 ,我们给出一些数值结果以验证这种针对 Dirichlet边界条件的偏微分程解法。     

并联式TBCC排气系统模态转换过程的流场特性研究

为了研究并联式TBCC（Turbine Based Combined Cycle）排气系统在模态转换过程中出口流场的变化特性及其气动性能的变化规律,采用动网格等技术完成了某并联式TBCC组合排气系统在整个模态转换过程中的非定常数值模拟。同时,对该排气系统模态转换过程中若干工况时的冷流进行了风洞试验,并将试验与数值仿真结果进行比较。研究结果表明:模态转换过程中,并联式TBCC排气系统出口流场的波系结构十分复杂,分流板出口激波对排气系统的气动性能产生了一定影响;TBCC排气系统的推力系数始终保持在0.9以上,但其产生的升力变化较大;风洞试验获得的壁面压力分布及流场纹影与数值模拟结果吻合较好,从而证明了本文数值模拟方法的可行性。