基于放大因子与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%.      

矩的显式积分算法研究与应用

矩的求解通常被用于求解有限元、体积、惯性矩等问题中.基于矩的叠加性,首先给出了在三维空间中计算域的离散方式,并推导了矩的显式积分公式,随后将其推广到n维空间中,该表达式易于在计算机上实现;设计了矩的并行计算算法,并通过Fortran和Python混编的方式,实现了矩的并行计算;对多重精度下的样例数据给出了一个算例,实现了零阶矩和二阶矩的计算,并和串行算法、逐次降维算法作出比较,进行了效率分析和误差分析.结果显示,矩的显式积分并行计算算法易于程序实现,并且在效率上高于串行算法,能够很容易推广到高维空间,该算法具有高度可并行性,误差主要来自计算域离散.      

基于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%,在他源数据上的安全预测率也很高,适用于航天器防护结构设计。探讨了撞击极限方程中的系数、速度分界值的优化思路,以提高撞击极限方程的预测率。     

小步长低阶Collocation方法轨道积分性能研究

介绍了Collocation方法的原理及误差构成,以Collocation方法在逐历元滤波轨道解算中的应用为切入点对小步长低阶Collocation方法轨道积分的性能进行了研究。分别就二体问题及真实GPS星座卫星轨道为对象进行了一系列验算,得到了一些有益结论：就积分方法本身性能而言,Collocation积分方法较传统数值积分法精度要高且稳定性要好;在一定长度的积分步长范围内,将小步长低阶Collocation方法用于滤波轨道解算是可行且有效的。 但在积分步长很小的情形下,从计算效率方面考虑,不推荐使用Collocation方法;利用小步长低阶Collocation方法取代大步长高阶方法能提高效率且可抑制蚀现象引起的积分结果发散。经试验表明：选用积分步长为900s的4阶方法是合适的。     

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

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

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

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