二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

分区和有限容积法计算长尾喷管中的湍流流动

为了给长尾喷管热防护结构设计提供热边界条件，从轴对称二维雷诺平均N－S方程出发，采用分区算法和有限容积法计算了长尾喷管中的湍流流动，分析了流场的结构特性，提出了长尾喷管中燃气温度，压强和马赫数分布的某些规律以及燃气温度沿壁面的分布，计算和分析结果提出了长尾喷管热防护结构设计应注意的事项。     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

带壁面射流燃烧室燃烧数值模拟

利用基于非结构化网格有限体积法对三维有壁面射流的燃烧室内两相流动和燃烧进行了数值研究.对气相流动在Euler坐标系下求解,而对液滴相则利用Lagrange方法进行追踪求解.计算区域采用四面体网格进行划分,气相流场用SIMPLEC计算方法,对液滴相采用了欧拉隐式方法.考虑了液滴相与气相的完全双向耦合作用,分别采用了Spalding液滴蒸发模型和涡破碎(EBU)燃烧模型,数值计算结果与文献中实验数据吻合较好.      

直升机附件舱温度场的数值模拟

为考查某型直升机附件舱的热控制能力,采用有限元软件分析了其在强制冷却和自然冷却情况下的温度场分布.应用非结构化网格和有限体积法进行计算区域和控制方程的离散,采用标准k-ε模型封闭湍流控制方程.针对3种冷却系统设计方案,分析了不同冷却气流进口尺寸、相对位置以及不同附件散热量对附件舱温度分布的影响,计算了自然冷却情况下附件舱密闭时的温度场.计算和试验结果均表明,所采用的冷却方案可以满足附件舱冷却的需要.      

S型进气道隐身与气动综合特性数值分析

针对S型进气道弯度对其隐身气动特性影响程度的问题,利用CAD软件对进气道进行参数化建模,通过改变进、出口的中心偏移量实现对进气道弯度的改变。采用数值仿真方法分别研究了进气道前向电磁散射特性和进气道流场及气动特性。仿真结果表明,进气道弯度增加对隐身特性和气动特性的影响是完全相反的,并进一步分析出这一现象的产生机理,对进气道的隐身气动一体化设计具有指导性意义。     

Determination of thermal expansion coefficients for unidirectional fiber-reinforced composites

In the present work, the coefficients of thermal expansion（CTEs） of unidirectional（UD）fiber-reinforced composites are studied. First, an attempt is made to propose a model to predict both longitudinal and transverse CTEs of UD composites by means of thermo-elastic mechanics analysis. The proposed model is supposed to be a concentric cylinder with a transversely isotropic fiber embedded in an isotropic matrix, and it is subjected to a uniform temperature change. Then a concise and explicit formula is offered for each CTE. Finally, some finite element（FE） models are created by a finite element program MSC. Patran according to different material systems and fiber volume fractions. In addition, the available experimental data and results of other analytical solutions of CTEs are presented. Comparisons are made among the results of the cylinder model,the finite element method（FEM）, experiments, and other solutions, which show that the predicted CTEs by the new model are in good agreement with the experimental data. In particular, transverse CTEs generally offer better agreements than those predicted by most of other solutions.     

跨声速机翼流场的N—S方程计算

本文给出了三维机翼粘性流场的数值模拟过程,控制方程为时均化的三维可压缩薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型,空间离散采用中心有限体积格式,时间采用Ｒｕｎｇｅ－Ｋｕｔｔａ多齿格式进行式积分。     

陶瓷基复合材料力学性能计算及涡轮导叶宏观响应分析方法

针对涡轮导叶宏观热力学响应,将化学气相沉积工艺2D编织SiCf /SiC陶瓷基复合材料完成拉伸试验并将该试件进行分 割,对分割后的试样进行断面电镜扫描,测量统计细观编织结构参数,建立考虑缺陷位置和形状的代表性体积单元模型。采用细 观有限元方法,应用该模型计算得到SiCf /SiC陶瓷基复合材料的宏观弹性常数。结果表明：面内拉伸模量预测值比试验值低 4.4%,剪切模量预测值比试验值高7.8%。结合陶瓷基复合材料涡轮导叶的加工工艺,建立了材料分布映射模型,对导叶典型特征 结构件在给定热载荷下的宏观响应进行预测。结果表明：涡轮导叶典型结构加热段尾缘第1主应变预测值比试验值低5.5%。     

纤维增强涡轮轴结构失效模式分析方法及试验验证

针对连续纤维增强复合材料涡轮轴结构失效模式分析问题,基于宏-细观力学跨尺度分析方法,建立细观力学代表性体积元(RVE)模型,通过编程模拟实现模型的周期性边界条件,计算纤维增强复合材料应力响应,将其均值应力转化为真实应力,确定失效包线。建立连续纤维增强轴结构力学模型,计算轴结构在扭转载荷下的应力响应。通过复合材料层合板主偏轴关系应力转化,将危险单元各方向宏观应力响应计算结果转化到细观力学RVE模型上,即为细观力学RVE模型受载情况。结合细观力学失效边界确定复合材料轴结构危险位置失效模式,当扭转载荷达到5 000～5 500 N·m之间,复合材料最外层即层6(+45°)首先达到基体拉伸失效载荷。开展复合材料轴结构失效模式试验,在扭转载荷达到6 000 N·m时,声发射信号相互叠加,大部分均为中频信号,中频信号多为基体、界面开裂信号。与模拟仿真计算结果对比分析,验证连续纤维增强复合材料涡轮轴结构失效模式分析方法的有效性。利用所建立模型预测了某型发动机低压涡轮轴的失效载荷及失效模式。