燃烧室内三维层流的数字模拟

本文以SIMPLE方法为基础,发展了一种模拟燃烧室内三维层流的数值计算方法和计算机程序.对存在冷却槽、主燃孔等不同流动条件下的流态进行了数字模拟.所得结果与已有的数值计算结果的一致性良好.     

两相紊流燃烧过程的数值计算

 本文描写了作者发展的两相紊流燃烧理论及数值方法;简介了作者自编的计算机程序GEMCHIP;提供了将本方法应用于几种燃烧室流动的计算结果。     

一种根据燃气成分计算燃气温度的方法

叙述了一种根据有限的燃气成分推算燃烧室出口高温燃气温度的计算方法，并分析了燃气成分测量误差、燃料低热值、燃料成分及燃料温度等对燃气温度的影响，该方法考虑了高温燃气解离的影响，利用数值解技术开展计算。根据此方法编制的计算机程序已应用于某高热容燃烧室出口燃气温度测量中，获得了燃烧室出口温度场.     

估计飞行器羽流影响的点源流动模型和计算机程序

目前和今后使用火箭发动机进行轨道和姿态控制的空间飞行器都要遇到发动机排气羽流影响问题.本文讨论了一种迅速估算轴对称的羽流远流场对飞行器影响的近似分析方法和相应的计算机程序.这种点源分析方法考虑了排气气体和喷管的特性.通过计算发动机排气羽流场,计算排气羽流干扰力矩和推力损失,这些影响在飞行器设计中必须给予考虑和控制.利用该方法和计算机程序对一个卫星模型进行了羽流影响的计算.     

三维瞬态导热温度场边界元计算及一阶数学奇异的处理

推导了三维瞬态导热温度场边界元解法的积分方程和离散代数方程，并编制了相应的计算机程序，采用坐标变换的解析消去法，处理边界元的一阶数学奇异点，从而提高了算法的数值精度。用本程序计算了正方体和圆柱环瞬态导热温度场，数值解与解析解的吻合程度令人满意。     

多级轴流压气机非设计性能的数值预估

采用LU-SGS隐式解法以及叶排间掺混面模型数值模拟压气机三维粘性流动,研究了多级轴流压气机的设计与非设计状态性能的预估问题。对一台五级高压压气机进行了详细的流场数值计算,包括了设计转速的100%,90%及80%等多个转速及各种工作状态。结果表明,本文所建立的数值方法与计算机程序对多级轴流压气机各种非设计状态的计算都具有收敛容易和快速的优点,计算得到的压气机总性能和各分级性能与设计参数进行了对比,证明该方法和程序的预估结果是比较准确的。根据设计转速下的计算结果,预估了压气机设计点的喘振裕度。此外,还对不同状态下叶片进口气流角的变化进行了分析。      

雷击建筑物时防雷系统中雷电流的分布研究

主要讨论雷直击建筑物时雷电流在防雷系统中的分布特性，介绍利用计算机进行数值计算的方法。通过典型算例、传统的电磁暂态计算程序——EMTP的计算结果和利用清华大学高压实验室的冲击大电流发生装置，对此结构模型进行模拟试验，将3种结果进行比较，从而验证了本研究所编制的计算机程序的正确性。为进一步研究室内电磁场的暂态分析做好基础。     

涡扇发动机加力燃烧室热阻损失的计算

加力燃烧室内热阻损失的计算虽然在许多著作里都有论述,但都作了一定的简化,作为工程计算,误差较大。本文在斯贝MK202发动机计算机程序中有关加力燃烧室热阻损失计算方法的基础上,探讨一种能同时计算出涡扇发动机加力燃烧室内涵燃烧区及外涵燃烧区各自热阻损失的一种工程上实用的计算方法。     

跨音涡轮中三维流动的数值试验

 本文从积分型的基本方程组出发,以我们的关于平面叶栅和任意旋成面叶栅的工作为基础,将文献（3、4）中的计算方案推广到三维流动情况。编制了相应的计算机程序,进行了一些数值试验,得到了初步的计算结果。 数值试验表明,拟三维计算结果可以导致相当大的误差,从而说明了完全三维流动计算的必要性。另一方面,如果取2000个左右的计算点,则在国产的013机上算一个例子,约需100分钟。从实际应用的观点来看也是可行的。     

不可压缩粘弹性有限元法及其对推进剂药柱应力分析的应用

由于固体推进剂具有近似不可压缩的特点 ,利用位移法有限元求解这类材料的力学问题不可能得到满意的结果。本文根据 Herrmann变分原理导出一种适用于不可压缩和近似不可压缩粘弹性物质的新的本构关系式 ,介绍了相应有限元公式 ,并开发了用于平面问题分析的计算机程序 VFAPINP。计算结果表明 ,本方法和程序可以用于所有泊松比的粘弹性问题的计算 ,尤其适用于固体推进剂药柱的应力分析。     

靶板抗轮盘碎片撞击能力的数值模拟

为研究机匣抵抗轮盘碎片撞击的能力和破坏方式,找到机匣结构优化设计的方法.用显式非线性动力学软件Dytran计算1/3轮盘碎片与单层和双层靶板的撞击过程.碎片及靶板均选用Johnson-Cook材料本构模型,结合Gruneisen状态方程.撞击过程基于Lagrange算法采用显式有限元计算.结果发现,靶板主要破坏方式为整体塑性变形、剪切撕裂和拉伸断裂;间距较大时双层靶板的弹道极限速度随间距的增加而增大;前靶板厚度比例较小或较大时弹道极限速度较大;单层靶板的抗击穿能力优于厚度均布的双层靶板.因此,通过对双层靶板的厚度和间距进行合理的搭配,能达到比同等厚度的单层靶板更好的抗击穿性能.      

涡轮盘疲劳寿命评估流程及应用

建立了较完整的涡轮盘疲劳寿命评估流程，并以一个涡轮转子模型的分析实例给出了寿命评估流程的实施过程，同时探索了涡轮盘考核部位模拟试验件设计的方法。结果表明:所建立的涡轮盘寿命评估流程具有较强的工程实用性，所发展的寿命方程参数确定方法可有效利用现有材料手册中各种材料基础试验数据并能获得一组物理意义明确的参数，此寿命方程预测的寿命精度在2倍分散带之内；通过涡轮盘计算分析发现应力梯度是影响涡轮盘中心孔、螺栓孔边以及过渡圆角等部位寿命的关键因素，需在设计过程中重点关注；通过带缺口的平板模型算例模拟了涡轮盘危险部位的应力梯度特征，实现试验件缺口部位的应力梯度与涡轮盘考核部位的相对梯度特征接近，并提供了缺口模拟试验件设计方法，为使用简单缺口试验件评估涡轮盘考核点的寿命提供了技术途径。      

采用碟形砂轮的面齿轮磨齿方法理论分析

为了制造出高精度硬齿面面齿轮和获得抛物线传动误差并提高传动稳定性,提出一种采用碟形砂轮加工面齿轮的磨齿方法.分析了碟形砂轮磨削面齿轮的展成原理和碟形砂轮的运动,根据展成原理推导了碟形砂轮的齿面方程,使用渐开线失配的碟形砂轮和改变砂轮的运动,推导出双向修形面齿轮的齿面方程.建立了双向修形面齿轮和常规渐开线小齿轮啮合的齿面接触分析模型,齿面计算和齿面接触分析实例表明,采用碟形砂轮加工双向修形面齿轮的磨齿方法是可行的,获得了面齿轮抛物线传动误差,避免了边缘接触并提高了传动的稳定性.      

基于表面缺陷特征的疲劳寿命预测方法

在含表面缺陷试样的疲劳数据的基础上，提出了表面缺陷对疲劳寿命影响的尺寸参数，将其引入Walker寿命方程，建立了可以考虑表面缺陷尺寸特征的疲劳寿命预测方程。将该方程的寿命预测结果同考虑应力梯度的寿命预测方法的计算结果进行对比，两者在±3倍以内，验证了方法是准确可靠的。进而，将该方程应用于粉末高温合金涡轮盘的疲劳寿命预测中，获得了不同尺寸的表面缺陷对涡轮盘寿命的影响规律，其工程意义在于:依据涡轮盘危险位置的应力特征，能够给出存在缺陷时的疲劳寿命，可作为使用过程中的重要参考数据，一旦出现漏检的表面缺陷，也能够保证涡轮盘的安全工作。      

压电网络对叶盘结构气动阻尼的影响

首先通过引入气动影响系数,建立了振动相关气动力与系统振动位移成比例的线性模型。基于此导出了叶盘结构气弹耦合的线性振动方程。接着在叶盘结构中引入压电网络,建立了流 机 电耦合动力学方程；并采用集总参数等效模型研究了压电网络对叶盘结构气动阻尼的影响,对流 机 电耦合系统的特征值问题进行求解,获得了可作为系统气弹稳定性判据的模态气动阻尼比。研究表明,压电网络可以有效改善叶盘结构主要振动形式——节径型振动的气弹稳定性；通过对压电网络中电学参数进行优化设计,可以在较宽范围内提高叶盘结构的颤振边界；压电网络对弯曲振动模态气动阻尼的改善效果优于弯扭耦合振动模态。      

飞力轮试验盘槽底裂纹损伤容限分析

用有限元和J积分相结合的方法, 确定飞力轮试验盘裂纹尖端参量随裂纹尺寸的变化规律, 并根据材料的断裂韧性确定临界裂纹长度。最后算出裂纹疲劳扩展寿命。计算结果与试验结果基本一致。      

处于非均匀温度场中含裂纹变厚度旋转盘的应力强度因子广义J积分解法

本文首先从直法线假设出发导出了处于非均匀温度场中变厚度各向同性钣元素本构方程并证明该方程与等厚度正交各向异性钣元素的本构方程相同。然后本文从上述本构方程出发,引用自由能与熵函数以及固体力学所有基本方程导出了与积分路径无关的广义J积分J_8,该广义J积分可用于处在非均匀温度场中的含裂纹变厚度旋转盘。最后本文应用这个广义J积分与有限元素法计算了涡轮盘应力强度.因子K_1并得到了令人满意的结果。     

盘用材料缺口特征尺寸对缺口强度的影响

现代航空发动机轮盘工况复杂,常常在榫槽、通气孔、螺栓孔等几何不连续处遭受破坏。这些几何不连续性特征可以视为广义的缺口,为了分析缺口附近几何特征尺寸对航空发动机轮盘复杂结构件抗拉强度的影响规律,基于典型盘用合金GH4169设计了双边缺口平板试样、双边不等缺口平板试样、不等厚双边缺口平板试样以及不等厚双边不等缺口平板试样进行试验研究。试验结果表明:缺口平板试样的应力集中系数与缺口抗拉强度存在一定的相关性,缺口抗拉强度无量纲化常数与应力集中系数的试验数据点存在一定的分布规律,二者的拟合曲线能很好地描述试验结果分布规律。包括GH4169、GH738和TC11合金等常用航空塑性金属材料在内的缺口平板试验数据点落在拟合曲线±10%以内的分散带内。     

小失谐叶轮的固有振动分析

失谐的带有 N个叶片的叶轮, 由于失去了 N阶回转对称性, 用于谐调叶轮的群论方法不再适用, 给其固有振动分析造成了困难。本文根据摄动理论, 以谐调叶轮的主模态作为 Ritz基, 提出了一种适用于小失谐叶轮固有振动分析的模态综合技术。算例表明, 文中方法是切实可行的。      

计算机翼跨声速压力分布的格林函数方法

本文给出一种用格林定理将跨声速小扰动速势方程化成积分方程式,然后进行数值求解的方法,文中对场元强度进行松弛迭代,一般只需20次左右即达收敛。算例表明,该方法结果与差分法结果一致。