用有限元法计算边界称动的喷管温度场

以固体火箭发动机复合结构喷管为对象，用有限元法进行了二维轴对称瞬态热分析。计算中考虑了因烧蚀引起的边界移动，且实现了网格的自适应，给出了考虑与不考虑烧蚀情况下的温度场计算结果，并将计算结果与试验结果做了比较，二者吻合较好，说明该预估方法既具有理论意义又有实用价值。     

用有限元法计算边界移动的喷管温度场

以固体火箭发动机复合结构喷管为对象，用有限元法进行了二维轴对称瞬态热分析。计算中考虑了因烧蚀引起的边界移动．且实现了网格的自适应，给出了考虑与不考虑烧蚀情况下的温度场计算结果，并将计算结果与试验结果做了比较，二者吻合较好，说明该预估方法既具有理论意义又有实用价值。     

固液混合火箭发动机燃烧的边界层计算

针对固液混合火箭发动机典型的边界层燃烧理论，用边界层的方法对固液混合火箭发动机燃烧进行了计算，计算采用二维轴对称的层流边界层方程，化学模型采用有限速率模型。计算考虑了不同氧化剂流率，不同吹入参数的多种情况，得到了固体混合发动机边界层的温度场分布，并利用计算结果，拟合了固体燃料的燃速公式，与相关文献的比较说明计算正确，结果符合燃烧机理和流动规律，为进一步研究固液混合火箭发动机的燃烧问题打下了基础。     

二维粒子分离器的流场及分离效率的数值模拟

数值模拟了某型惯性粒子分离器(简称IPS)的二维有粘轴对称流场分布、及各种不同尺寸砂粒的轨迹及分离效率,采用贴体坐标网格,并将多连通域转换成单连通域,在任意非正交曲线坐标系下求解IPS的分叉流动、粒子运动轨迹和分离效率。计算与实验结果吻合良好,表明:本文所采用的方法和编制的程序能反映粒子分离器的主要特征。      

长尾喷管两相流流场计算

为了给长尾喷管热防护结构的烧蚀和温度场计算提供准确的热边界条件，从二维轴对称N－S方程出发，采用颗粒轨道模型和有限体积的Jameson格式计算了长层喷管纯气相流场和两相流流场，分析和比较表明：在长尾喷管喉部附近以及扩散段的相同轴向位置，两相流情况下轴线上的燃气马赫数小于纯气相情况下的燃气温度；在进行长尾喷管热防护结构设计时应以两相流流场作为热边界条件。     

二维有限元模型在燃气轮机转子-支承系统临界转速计算中的应用

应用专业的转子动力特性分析软件SAMCEF／ROTOR，建立了燃气轮机转子的二维有限元分析模型，计算了该转子系统的临界转速，并与其它经典方法的计算结果以及试验结果进行了比较，验证了用二维轴对称有限元模型计算转子临界转速的工程实用性和可靠性。     

气氢液氧同轴式单喷嘴燃烧室热态流场计算

建立了欧拉－拉格朗日坐标系下液体火箭发动机内部工作过程二维轴对称两相湍流化学反应流的完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ模型。应用该模型对气氢液氧同轴式单喷嘴燃烧室热态流场进行了数值模拟，得到了全流场的速度矢量分布、燃气组分等值线、流场等温线、流场等压线和流场等马赫线。计算是在任意斜交曲线坐标系下进行的，采用等参数正变换（代数变换）把物理平面上不规则计算域变换成了计算平面上的矩形域，控制方程也进行了相应的变换。计算结果表明，用计算流体力学（ＣＦＤ）的方法研究和分析液体火箭发动机内部工作过程是可行的。     

一种计算再生冷却推力室温度场的方法

为了能够快速而准确地得到再生冷却推力室的温度分布,建立了一种计算再生冷却推力室温度场的方法。首先建立了轴对称推力室的一维冷却模型,并使用换热经验公式,得到了推力室壁面在轴线方向上的温度分布;其次建立了推力室的冷却套二维导热模型,使用数值模拟的方法和一维计算的结果,得到了冷却套的温度场。然后使用这种方法研究了气壁材料、气壁厚度和冷却液流量对推力室再生冷却的影响,获得了比较满意的结果。从计算时间和准确性来说,这种方法能够为推力室的优化设计和性能估算提供参考。      

串装双燃速药柱发动机的内流场计算

对一维的加质燃烧室流场、轴对称的二维喷管内流场和高燃速药柱燃烧完成后的燃烧室二维内流场以及三维药柱几何做了一体化计算。应用轴对称的N－S方程、TVD有限体积法的显式MacCormack格式，对燃烧室后段和长尾喷管跨声速流场进行求解，得到了双燃速内弹道性能和轴对称二维流场中燃气参数的分布。结果具有较高的精度，并为进一步的燃烧室和喷管热结构分析提供数值依据。     

发动机短舱内外流场与结构温度场耦合计算

采用CFX软件,利用SST k-ω湍流模型,同时引入了辐射模型考虑热辐射的影响.采用CFX的共轭传热方法,在流体域求解N-S方程,固体域求解热传导方程,然后在流固耦合交界处交换数据,实现了流场与结构温度场的耦合数值分析.为了验证此方法的正确性,对文献中轴对称喷管进行了数值模拟,所得喷管壁温与试验数据吻合良好;另外进行了高超声速钝体流场与结构温度场的耦合计算,计算所得激波位置与试验结果一致,固壁非定常温度场分布与文献结果吻合.在此基础上,进行了航空发动机短舱内外流场与结构温度场耦合的数值研究,并得出了几个短舱外表面蒙皮冷却的可行性措施,为飞行器红外隐身设计提供参考.     

孔冷挤压有限元仿真中的铰削分界面位置确定方法

开展孔冷挤压过程有限元仿真计算是残余应力分布获取和疲劳寿命预测的前提。在有限元建模阶段，设置铰削层单元与基体材料单元之间的分界面，是模拟铰制终孔工艺过程的关键。通过弹塑性力学分析，建立了挤压强化过程芯棒、衬套和被挤压强化连接孔的应力分析方法；基于分析中得到的不同位置处微单元的径向位移量，建立了铰削分界面相对位置计算模型。并开展了关键参数的敏感性分析，定量研究了关键参数变化对残余应力分布和径向位移量的影响程度。本工作为孔冷挤压强化有限元模型建立中，铰削层单元与基体材料单元分界面相对位置确定，提供了便捷可靠的方法。     

斜喷管收敛段壳体三维应力场的有限元法计算

对固体火箭发动机斜喷管收敛段壳体应力场采用有限元法进行了分析和计算。提出了该方法的力学模型，基本公式，基本方法和算例。具有理论意义和工程应用价值。     

高、低周复合疲劳及其试验技术

文中介绍了一种简单易行的复合疲劳试验机,在此试验机上进行了高、低周复合疲劳试验。结果表明,随着次循环的频率和载荷幅增大,铝合金(LY12—CZ)试件的裂纹萌生寿命显著下降。然而,当次循环频率很低时,寿命也下降,这似乎属于“变形损伤”,因缺口附近的应变出现“棘轮”现象。复合疲劳的损伤模型尚有待研究。     

圆柱壳机匣热应力试验研究

试验设备与温度场的形成 试验设备与测量系统(图1):包括加温系统、温度测量系统和热应力测量系统。试件是1Cr18Ni9Ti材料制成的圆柱壳机匣,固定在底座上。加温器由周向均布的300V碘钨灯组成。热电偶是温度测量的传感器。在热应力测量中,采用中温应变片作为传感器,把电讯号经电阻平衡箱传到YJ-5应变仪,录取应变数据。 形成所需要的温度场是开展热应力试验研究的前提。      

2.5D C/SiC复合材料带孔板拉伸破坏的多尺度模拟

为了研究2.5D编织陶瓷基复合材料带孔板的拉伸破坏行为,提出一种可以模拟带孔板细观破坏过程的多尺度计算方法。该方法根据2.5D编织陶瓷基复合材料的细观结构建立细观模型,通过子模型法将平板的宏观有限元模型和孔周围区域的细观有限元模型耦合在一起,然后采用渐进损伤计算方法完成带孔板破坏的多尺度模拟。计算结果表明,带孔板在拉伸载荷较低时出现初始损伤,随着载荷的加大经纱发生轴向拉伸破坏,纬纱发生横向的破坏。裂纹从孔边沿横向扩展至板的两端,最终整个板完全断裂失效。失效时的应变为0.375%,最大加载应力为221.7MPa。     

固体火箭发动机的声学分析及燃烧稳定性预估

提出一种预估固体火箭发动机燃烧稳定性的方法。文中着重研究了燃烧室内腔声场的有限元数值分析问题,并应用声学有限元素法对一个二维轴对称装药发动机燃烧室内腔进行了固有频率及固有振型计算。随后,在声学分析的基础上对该发动机做了燃烧稳定性预估。其结果与多次试车的结果一致。最后还就声学分析的验证,减少机时和内存,微粒尺寸的确定,预估方法及计算程序的通用性等问题进行了讨论。     

基于非协调单元的液体火箭发动机推力室热结构分析

用巴兹公式和发动机的热力计算数据,得到了推力室内的燃气壁面对流传热系数和燃气温度的分布,并用有限单元法计算了推力室的温度场.在温度场计算基础上,采用Wilson非协调单元,对推力室的应力场进行了计算,其中位移函数采用二阶形式,将计算精度提高到二阶.结果表明:相比常规单元,Wilson非协调单元在计算结构比较复杂,温度梯度比较大的结构时,精度比较高,结果合理.      

复合材料叠层板固化成型后的变形分析计算

分析复合材料铺层的微观结构，研究复合材料成型过程变形特点，以此建立两类针对性的有限元模型，并采用典型复合材料叠层板固化成型温度载荷，进行叠层板翘曲变形和板内应力的计算分析。在此基础上，将计算结果与实际试验制件实测结果进行对比。结果表明，本文所提供的有限元分析计算模型可对实际情况进行简便和较为精确的模拟。     

弧线齿面齿轮应力过程分析

为了给弧线齿面齿轮的齿面接触强度和齿根弯曲强度设计提供理论依据,研究了弧线齿面齿轮的齿面接触应力和齿根弯曲应力随载荷和安装误差的变化规律.在齿面接触分析和承载接触分析的基础上应用弹性理论计算了弧线齿面齿轮副的齿面接触应力和应用有限元应力影响矩阵法计算了该齿轮副的齿根弯曲应力.给出了数字计算实例,计算结果表明:齿面接触强度和齿根弯曲强度在重载时的接触强度和弯曲强度由单齿啮合区的强度决定,轴向安装误差和轴夹角安装误差分别会增加齿面接触应力和齿根弯曲应力,轴夹角安装误差和轴间距安装误差对齿面接触应力影响甚小,而轴向安装误差和轴间距安装误差可以降低齿根弯曲应力,与直齿面齿轮相比,弧线齿面齿轮的接触和弯曲应力明显减小.      

焊接顺序对TA15钛合金壁板焊接变形的影响

建立了TA15钛合金壁板焊接有限元模型,利用热弹塑性有限元法分析了焊接顺序对壁板变形的影响。通过对比焊缝截面宏观形貌和残余应力,建立了TIG穿透焊焊接工艺的双椭球热源模型,进而模拟了5种焊接顺序下壁板变形情况。结果表明:按从一侧开始的顺序进行焊接时,下壁板中部边缘凸起,最大变形量为3.9mm;采用对称施焊,可有效地控制壁板变形量。