基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.      

基于Deform软件二次开发和BP神经网络的TA15多向锻造微观组织预报

通过热压缩实验构建了TA15合金的热变形真应力-应变曲线,以此为基础分别建立了合金双相区及单相区温度区间的热变形本构方程;基于热压缩试样动态再结晶的统计数据建立了TA15合金的动态再结晶模型。借助Deform提供的二次开发功能实现相关数学模型的程序化,制定正交实验方案,实现了TA15合金多向锻造变形的微观组织仿真。通过正交实验分析得出各项因子的影响对象及强弱差异,提出了双相区及单相区温度区间内的多向锻造最佳因子组合。建立了TA15合金多向锻造变形微观组织的BP （Back Propagation）神经网络预测模型,将预测结果与有限元仿真结果进行比较,结果表明两种方法的预测结果基本一致,但神经网络具备有限元仿真难以实现的良好细节预测能力,能更为细致地实现对微观组织分布状态的划分。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

蜂窝夹层板超高速撞击极限方程分析

为分析撞击极限方程对蜂窝夹层结构的穿透特性的预测能力,调研得到了3类撞击极限方程的表达形式和等效方法,以及131个采用碳纤维复合材料(CFRP)面板的蜂窝夹层板结构的试验数据,并对撞击极限方程的预测能力进行了比较计算。结果发现,MET方程对他源数据的四种(未失效、失效、总体以及安全)预测率均大于80%,进行在轨航天器结构的失效分析时可优先选用;SRL方程对本源数据的安全预测率达到了100%,在他源数据上的安全预测率也很高,适用于航天器防护结构设计。探讨了撞击极限方程中的系数、速度分界值的优化思路,以提高撞击极限方程的预测率。     

各向异性Chaboche黏塑性本构方程隐式应力积分算法

针对各向异性Chaboche黏塑性本构方程,利用有效参数控制法推导了该方程的隐式积分算法,该算法把所有未知量都用同一个控制参数表示,并形成控制方程,从而把本构方程应力积分的问题转化成一个非线性控制方程的求解问题.该积分算法程序,与通用有限元软件Marc结合,通过单个单元对单向拉伸、循环、松弛的响应的实验数据进行了数值模拟,通过与显式算法对比,该隐式算法的稳定性和收敛性有所提高.构件的对比分析也表明隐式算法比显式算法具有较高的效率.      

叶片全三维反问题优化设计方法

基于数值求解Navier-Stokes方程,对叶轮机械叶片全三维反问题优化设计方法进行研究.反问题给定叶片表面的静压分布反求叶型,假设叶表网格点有虚拟移动速度,迭代过程中根据叶表的实际静压和目标静压之差来修改叶型.对单排叶栅进行了反问题设计,并给出了初步的优化算例,结果能很好地满足目标静压分布,验证了该方法的有效性,计算的时间和精度也适合进行工程应用.      

等极孔球形压力容器平面缠绕规律

基于空间解析几何理论及运动学分析方法,建立复合材料球形压力容器平面缠绕落纱点轨迹运动方程,推导出等极孔球形容器表面缠绕纤维落纱点空间位置、速度以及加速度的求解函数;通过对落纱点运动规律的分析,解出空间落纱点各方向运动速度的大小及其各方向夹角余弦值。最后利用Matlab软件对不同的丝嘴运动平面倾斜角的缠绕轨迹进行模拟,通过该计算方法模拟的轨迹结果与实际球形容器缠绕线型相符合,验证了该算法的有效性。     

射频干扰对InSAR干涉处理的影响分析

窄带干扰是低频段InSAR系统面临的主要射频干扰之一,它的存在会对干涉相位产生严重影响,进而导致高程反演误差。本文推导了窄带干扰影响下的InSAR干涉相位的解析表达式,给出了干涉相位对窄带干扰各参数的敏感度方程,仿真分析了窄带干扰对P波段InSAR系统干涉处理的影响。结果表明,窄带干扰经过成像后其幅度调制函数近似为方位时间的sinc函数,在SAR图像中表现为沿距离向的干扰条带。干涉相位对窄带干扰信号的功率较敏感,对其频率较不敏感,干涉相位和高程反演误差随着干信比的增加而增加。     

过渡状态倾转旋翼气动力模拟的高效CFD方法

为显著减少倾转旋翼过渡飞行时气动力CFD模拟的计算代价,提出并建立了适合倾转旋翼过渡状态气动特性分析的高效混合CFD方法。首先,提出了适合于过渡状态模拟的嵌套网格系统,并发展了相应的挖洞和贡献单元搜寻方法。在此基础上,结合叶素理论和动量理论建立了旋翼气动力模拟的简化虚拟桨盘模拟方法(Virtual rotor model,VBM)。为了能够分析倾转旋翼气动力的细节特性,多层运动嵌套网格系统和单指令多数据流(Single program multiple data,SPMD)并行技术被引入来建立精确的旋翼模拟方法(Real blade model,RBM)。然后,通过将VBM和RBM方法结合,构建了适合倾转旋翼过渡状态气动特性分析的高效(Hybrid blade model,HBM)方法。最后,通过对有试验值对比的悬停状态典型旋翼和7A旋翼分别验证了VBM和RBM方法的有效性。分别采用3种方法预测了过渡状态不同倾转角下旋翼的气动特性,VBM表现出最优的计算效率,能用于倾转旋翼总体气动性能的分析。HBM方法在保证流场求解精度的基础上,相对于高精度的RBM方法节省了1/3的计算时间。     

电推力器中复杂电场求解的超松弛并行算法

目前,电推力器中复杂电场的三维泊松方程求解仍存在速度偏慢、并行性欠佳等缺陷。通过分解三维求解区域,针对复杂电场提出一种超松弛并行迭代算法。将求解区域划分为多个子块,利用超松弛迭代格式构造出若干分组显式格式,分别给出不同迭代步数下的求解方程以进行并行计算。通过数值模拟,计算时间至少缩短一半以上,该P SOR方法较传统迭代格式在快速性方面取得较大进展,对电推进领域的数值仿真研究起到促进作用。