基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

基于应变路径非比例度的多轴疲劳寿命预测

在分析多轴疲劳几种常用非比例度定义的基础上,提出了一种非比例度定义方法,进而以American Society of Mechanical Engineers(ASME) 规范案例中非比例加载多轴疲劳设计准则采用的应变参量作为基本损伤参量,发展了一种新的多轴疲劳寿命预测模型.结果表明:①所提出的非比例度定义可以描述任意已知轮廓的、非周期的、变幅的非比例加载路径;②与两种常用的多轴非比例加载疲劳寿命模型的预测结果对比可知,新的寿命预测模型对14种比例和非比例加载路径下304不锈钢材料的寿命预测与试验吻合更好,预测结果基本位于2倍分散带以内.      

幂律型流体撞击射流破碎特性直接数值模拟

凝胶推进剂的双股射流撞击雾化广泛应用于液体火箭发动机的燃烧室中,其破碎特征及雾化效果直接影响燃烧效率。为探究雾化特性的发展规律,采用直接数值模拟DNS方法,对射流速度为100m/s的剪切稀化非牛顿液体正交撞击产生的雾化特征、液体表面积、表面波、涡特性以及非牛顿特性开展研究。结果表明,射流下形成的雾化流场迅速扩张形成液膜,液膜两侧边缘破碎成大量的液丝与液滴,核心部分产生撞击波后在气体力的作用下逐步发展为带有凸起和褶皱的不稳定表面波,其撞击波波长最大可达2.46倍射流直径。液体表面积不断增长,但无量纲表面积总体呈现先下降再上升的趋势。气体中的涡量分布则分为有序附着区和无序爆炸区两类,并且涡量主要集中分布于气相区域。此外,射流撞击时产生强剪切使该液体内部的粘性系数下降,最低仅为初始粘性系数的0.3倍。     

基于PaSR模型的低旋流燃烧大涡模拟

为了深入理解低旋流流场特征和燃烧稳定性,基于OpenFOAM平台,采用动态k方程模型和有限速率PaSR模型对甲烷/空气预混低旋流燃烧进行了大涡模拟,研究了气流入口速度、当量比和压力等流场参数对流场结构和燃烧非稳态特性的影响,分析了流场大尺度结构与火焰相互作用。结果表明,流场结构和火焰抬升高度受入口速度影响较小,流场和火焰形态能够保持自相似性;随着当量比和压力提高,流场扩张性增强并在燃烧区下游产生回流区,火焰稳定不依赖回流区,根部火焰锋面形状由U形转变为W形,火焰抬升高度降低。火焰锋面稳定在剪切层,剪切层产生的周期性有序涡结构引起当地流场速度脉动和火焰表面褶皱,反映了流场非稳态特性;通过剪切层监测点瞬时轴向速度分析,涡结构特征频率随速度增大而提高,由250Hz提高至300Hz,随当量比和压力提高而降低,由250Hz降低至125Hz。     

民用航空发动机引气系统部件寿命建模

针对飞机维修工作中发动机引气系统故障率较高问题,以寿命数据分析为基础,运用威布尔分布模型定量评估航空发动机部件可靠性,建立了民用航空发动机引气系统部件可靠性寿命模型。以V2500发动机引气系统为例,收集了该系统的可靠性数据,分析了部件的非计划拆换情况。在数据计算过程中采用3参数威布尔分布模型对引气系统部件的寿命数据进行拟合。假设检验结果表明,在对部件进行寿命可靠性分析时,运用威布尔分布建立的数学模型符合客观规律,且能够取得良好的拟合效果。依据计算出的部件可靠性特征量对部件的可靠性状况进行定量评估,结果表明:在评估复杂机械电子系统的可靠性时,威布尔分布模型具有较大应用价值。     

基于非接触式测量的旋转叶片动应变重构方法

基于叶端定时非接触式测量和振动响应传递比的概念,开展高速旋转叶片动应变重构方法的研究。在频域内推导了叶片任意测点位移与任意测点动应变的传递比,给出了单模态共振下响应传递比关于位移和应变模态振型的解析表达式;建立旋转叶片的三维(3D)有限元模型,开展考虑旋转预应力效应的叶片模态分析,提取位移和应变模态振型,获得任意转速下叶端位移与叶根关键点动应变的传递比。开展高速旋转叶片叶端定时非接触式测量实验,采用周向傅里叶算法对叶端定时信号进行处理,获得叶片在不同转速单模态共振下的叶端位移,结合响应传递比,重构5个旋转叶片的关键点动应变。结果表明:旋转叶片在9000r/min和13000r/min转速下发生1阶共振时,与应变片实测结果相比,叶根处应力最大点、次大点和边缘点3个关键点的动应变平均重构误差均小于15%,验证了旋转叶片动应变重构方法的有效性。     

非正交弧线齿面齿轮齿面设计及根切研究

非正交弧线齿面齿轮是以端面渐开线弧线齿圆柱齿轮为假想刀具包络展成的面齿轮。根据面齿轮啮合原理推导了非正交弧线齿面齿轮齿面方程,结合Matlab数学模型和根切理论研究了该种面齿轮的根切现象,计算出面齿轮不产生根切的最大内径,并通过Catia仿真滚齿来验证计算的准确性。该研究旨在提出一种新型非正交面齿轮,通过理论推导及仿真模拟来研究该种面齿轮的根切现象,在此基础上探究非正交弧线齿面齿轮的内径影响因素。研究结果表明:在同等参数下,当位置参数增大时,刀具渐开线截面高度和内径逐渐减小,当刀具圆弧半径增大时,内径逐渐减小,刀具渐开线截面高度逐渐增大。     

TGO非均匀增长对热障涂层应力演化和破坏机理的影响

采用内聚力模型和热生长氧化层(TGO)非均匀增长子程序,数值模拟了在热循环载荷作用下热障涂层(TBC)内部应力演化规律和开裂行为。涂层失效过程首先是源自陶瓷层(TC)内近波峰位置的拉伸和切应力共同主导的陶瓷层Ⅰ、陶瓷层Ⅱ混合型裂纹;随着循环数增加,则转向由TC内近波峰位置的切应力主导的Ⅱ型裂纹和波峰波谷中间的涂层厚度方向拉伸应力主导的Ⅰ型裂纹。整体非均匀增长和波谷均匀增长模式下的最大拉伸应力经过一定循环数后几乎不再随循环数而增加;而在波峰均匀增长和整体均匀增长模式下,最大拉伸应力则会随着循环数增加持续增长。整体非均匀增长、波谷非均匀增长模式下,20个循环后最大切应力出现在近波峰位置,分别为-16241 MPa和-15428 MPa;而整体波峰均匀增长和整体均匀增长模式下,最大切应力为-11382 MPa和-11198 MPa。对于波谷均匀增长和整体非均匀增长模式,在9个循环后出现界面裂纹。而对于波峰均匀增长和整体非均匀增长模式,在第17个循环出现界面裂纹。     

基于截尾概率-非概率混合模型的可靠性优化算法

针对工程中截尾概率变量与非概率变量同时存在的情况,给出一种新的截尾概率与非概率混合可靠性模型。在该混合可靠性模型基础上,按照可靠性指标(RIA)法给出双层嵌套可靠性优化模型,并采用改进搜索策略后的ST-Powell优化算法在外层搜索设计变量的最优值,内层采用能保证收敛的改进的有限步长迭代法求解混合可靠性指标。数值算例表明,改进搜索策略后的ST-Powell优化算法的全局寻优性得到显著提升;改进搜索策略后的ST-Powell优化算法与改进的有限步长迭代法相结合求解双层嵌套混合可靠性优化模型的正确性得到验证,且对于非线性程度较高的极限状态函数同样能够得到满足截尾概率与非概率混合可靠性模型指标要求的最优解,并对工程结构算例具有很好的适应性。     

低雷诺数下50°后掠三角翼的旋涡流动

采用数值模拟和流动显示的方法研究了50°后掠角三角翼在低雷诺数下的旋涡流动,结果表明:低雷诺数下,非细长三角翼在5°攻角时就形成了稳定的前缘涡,较小攻角时前缘主涡就开始破裂,并观察到泡型和螺旋型两种旋涡破裂方式。另外,在一定的攻角范围内,前缘主涡的外侧又生成一对新的集中涡,构成双涡结构;随着攻角的增大,前缘涡涡核不断升高,主再附线向中心移动,二次分离区扩大。