在烧蚀条件下固体火箭发动机复合喷管的温度预示

本文探讨了由石墨和高硅氧-酚醛组成的固体火箭发动机复合喷管,在考虑热解炭化和烧蚀引起喷管内型面变化条件下的瞬态轴对称二维温度场.借用当量热容的概念和坐标变换的方法,应用单元体能量平衡原理,在内点用显式格式,边界结点用隐式格式进行数值求解,这样的计算比较简捷,所得温度场与实验结果作了比较,说明本文的预示方法是可以符合工程精度要求的.     

硅基内衬喷管扩张段烧蚀和温度场的耦合计算方法

本文叙述固体火箭喷管扩张段硅基内衬烧蚀和温度场的耦合计算方法.烧蚀采用非定常的液体层烧蚀模型;温度场计算采用圆柱坐标系下的一维径向瞬态导热模型,考虑材料在高温下形成的炭化层和热解层,而由烧蚀产生的边界移动问题用坐标变换的方法处理.利用烧蚀和导热计算方程组中二者的结合部,通过充分耦合,获得了更准确的预示计算结果.     

烧蚀端头的瞬变外形及内部温度分布

本文通过对传热、烧蚀和形状变化的耦合计算取得了各向同性和各向异性石墨再入端头帽的温度场。直接从热传导偏微分方程出发,通过座标变换,建立三层无条件稳定的显式差分格式,计算简捷。在碳表面,包括十七组元的平衡烧蚀计算表明,壁面的主要化学成分有N_2、CO、C_2N、C_3、CN等几种。对非线性形状变化方程,采用可显式计算的隐式差分格式,稳定性良好。二种石墨端头的温度场差别很大,但算例条件下的总烧蚀量相差无几。     

平面温度场的延拓有限元分析

 平面温度场的延拓有限元分析 大多数复杂温度场问题的解析解往往是难以获得的,必须采用一些近似的数值计算方法。延拓有限元技术合理利用了单元邻近外点信息,以构造单元内高精度广义插值近似函数,从而在不增加网格剖分单元及结点数目的条件下,有效地改善温度场的分析精度。     

高超声速飞行器碳基头锥烧蚀外形计算

高超声速飞行器头锥部位的烧蚀外形会影响飞行器的升阻比和飞行器的稳定性,烧蚀外形的准确评估对于长时间滑翔式飞行器的综合设计具有重要意义。本文介绍了头锥烧蚀外形计算方法,基于该方法对飞行器头锥烧蚀外形进行了仿真研究,某飞行器头锥烧蚀的仿真结果与试验结果吻合良好。对比了飞行器进行俯仰运动和滚转运动时的烧蚀外形,总结了烧蚀外形随飞行状态变化的演变规律,俯仰运动影响头锥迎风面和背风面的烧蚀量,滚转运动影响头锥烧蚀位置。     

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

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

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

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

烧蚀环境下的圆柱绕流计算模型

基于多学科理论建立了高温高压气流环境下圆柱烧蚀、剥蚀的数理模型。利用离散涡方法计算流场与圆柱表面压力分布;采用三方程烧蚀模型计算热化学控制机制下的烧蚀速率,并将其与扩散控制机制下的结果相比较以确定烧蚀量;引入颗粒轨道模型求解剥蚀颗粒的运动。基于该模型对高温高压燃烧室内圆柱形烧蚀试件的绕流场、烧蚀量以及剥蚀颗粒的运动轨迹进行了编程计算分析。研究表明:低雷诺数条件下烧蚀环境对绕流流场的影响不大,而较高雷诺数条件下烧蚀绕流流场与无烧蚀绕流流场的涡街分布存在显著差别;在低雷诺数条件下流场对烧蚀速率的影响甚微,较高雷诺数条件下烧蚀速率略大于无流场烧蚀情况;涡量的存在决定着剥蚀颗粒的运动轨迹与分布,进而影响尾流热场与烧蚀速率。该模型可为热防护材料的烧蚀实验提供参考,并为烧蚀、剥蚀过程的非线性分析提供数据支持。     

高强化活塞烧蚀失效研究

运用软件AVL FIRE对某型柴油机进行了有近壁燃烧工况时的燃烧特性的数值模拟,同时确定了该型柴油机在有近壁燃烧工况时活塞有限元计算的热边界条件,运用软件ANSYS进行了活塞的烧蚀过程有限元分析,得到了近壁燃烧对活塞烧蚀的影响规律,进行了活塞烧蚀的实机试验,模拟计算结果与试验结果基本吻合.结果表明:所建立的烧蚀计算方法可以很好地模拟活塞的烧蚀过程;柴油机喷油器喷孔磨损变大将导致活塞喉口附近燃气温度和壁面传热系数在较大的曲轴转角内都保持较大的数值,从而使活塞发生烧蚀失效;烧蚀量随着喷油器孔径的增大先增大后减小,在孔径增大15%时达到峰值;烧蚀量随着烧蚀时间的增大而增大,但烧蚀的速度在缓慢减小.      

高超声速三维碳—碳烧浊流场的数值研究

本文针对高超声速再体的烧蚀现象 ,利用简单隐式TVD差分格式和激波捕捉法 ,数值求解三维化学非平衡Navier Stokes方程 ,其中化学模型是碳 碳 (C C)空气化学模型 ,考虑 12个化学组分和 31个化学反应过程 ,研究了C C烧蚀对再入体头部区域的壁面温度和热流分布的影响。为了计算效率和稳定性提出壁面条件显式处理的方法。对再入高度为 6 5km和速度为 8km/s的再入体头部区域烧蚀流场进行了数值模拟 ,用飞行迎角α =0°的计算结果与国外文献进行了比较 ,符合得较好。同时给出了三维小迎角α =5°的计算结果     

飞机外挂物离机过渡过程的有限元算法.

研究飞机外挂物在按给定轨迹下沉旋转的离机过渡过程中，流场非结构网格的调整和再生成方法以及过渡过程欧拉方程有限元数值解法。中指出外挂物运动时流场网格调整或再生成仅需在围绕在外挂物的一小区内进行以减少计算量，而且固定网格的欧拉方程的有限元数值解法也可适用于运动网格情况，只要将固定网格的守恒通量项Fk用运动网格的Ek来代替，并建立了Fk和Ek之间的关系。末给出了带外挂物飞机的实例计算结果。     

任意坐标系中的弹塑性矩阵

 本文推导了三维问题及平面应力、平面应变、轴对称问题在任意曲线坐标系中的弹塑性矩阵。由此最一般的表达式可以方便地得到任何特定坐标系中的弹塑性矩阵,供各类不同特点问题的有限元计算时使用。 文中推导得到的任意曲线坐标系中的本构方程,也同时为曲线差分法求解弹塑性问题提供了刚度系数显式。     

径向开槽药柱发动机内流场数值模拟

采用时间相关的显示ＭａｃＣｏｒｍａｃｋ预估－校正两步格式和可压缩全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了轴对称径向开槽药柱发动机内流场。计算中采用了贴体坐标系统和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。边界参数计算采用参考平面上的特征线法和物理边界条件。通过大量算例的分析表明，方法有效，计算结果合理。当有集中质量加入时，燃烧室头部压强升高，喷射强度和开槽口的倾斜方向对头部压强也有较大的影响。     

动载下的圆柱齿轮接触应力分析

本文研究了一般三维弹性接触问题的有限元局部网格密化技术 ,并在此基础上结合有限元—线性规划法首次提出了一种计算直齿和斜齿圆柱齿轮啮合过程接触应力分布的方法。文中以一对航空斜齿圆柱齿轮为例 ,计算了在系统动载和静载作用下该对斜齿轮不同接触线上的接触应力 ,并与航标进行了对比分析。计算结果表明本文方法及其微机用程序系统是快速可靠的。     

贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.      

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

采用近似边界求解壁板颤振的CBS有限元算法

为了简化Euler方程求解壁板颤振问题,将近似边界条件成功地推广到CBS(characteristic-based split)有限元方法求解壁板颤振领域,把运动壁面的影响转化为Euler方程的边界条件,避免了对动网格的处理.通过对近似边界方法、动网格方法以及3阶活塞理论的数值模拟结果进行对比发现:首先,加载近似边界的CBS有限元方法操作非常简便,且模拟结果与动网格方法完全吻合;其次,Euler方程可以准确捕捉到壁板前后两端不光滑尖点所引起的压力突跳,而3阶活塞理论则无法做到;最后,采用近似边界的CBS有限元方法在求解无黏可压缩流动中细长体小变形的流-固耦合问题上具有潜力和优势.