一种用于飞机结冰的海绵冰模型及应用

基于两相流理论,采用欧拉法对飞机结冰进行数值模拟,在空气流场速度分布的基础上,求解水滴流场,得到水滴收集率和撞击极限.结果表明:欧拉法比拉格朗日法在计算水滴收集效率和撞击极限方面有优势,尤其是在较大的水滴平均等效直径情况下.对翼型的冰形和圆柱表面的温度预测数据和实验数据进行对比,圆柱冰形厚度误差在15.0%之内,计算结果与实验结果吻合较好,同时表明结冰模型是有效可行的.      

悬停状态直升机桨叶扭转分布的优化数值计算

建立了一套基于高精度计算流体力学(CFD)技术和代理模型优化算法的旋翼气动外形设计方法。在该方法中,旋翼流场气动性能的计算采用了基于Navier-Stokes/Euler方程的CFD方法,并根据流场特点、精度和效率的要求采用Baldwin-Lomax(B-L)湍流模型,通量计算采用Roe-MUSCL格式进行。为提高网格生成质量和便于流场控制方程的求解,将流场分成两个区域,即围绕旋翼的黏性区和无黏的背景网格区。其中,桨叶网格使用了基于二维翼型网格的参数化方法生成,数值计算结果表明该方法有效地提高了网格生成质量及效率。在参考旋翼流场及桨叶细节流动分析的基础上给出设计变量及范围,有效减小了优化问题的规模;为满足优化和机理分析的需要,将基于置换遗传算法优化的拉丁超立方(PermGA LHS)方法和径向基函数(RBF)的代理模型优化方法引入到桨叶外形的优化设计中。首先以Helishape 7A旋翼为算例,检验了数值模拟方法的准确性。然后,应用所建立的优化方法针对旋翼负扭转分布进行了优化计算,结果表明优化后的旋翼悬停气动性能比优化前有了明显提高。     

三阶幻方问题的代数解法

介绍了中国古代数学家甄鸾和杨辉关于三阶幻方的一个构造方法。用线性代数的方法构造三阶幻方约束方程组时,当自由未知量分别取1,2,3,…,9时,对应的方程组共有56组整数解,为了找出幻方的全部解,需要确定自由未知量取值的限制条件,由此对方程组中的自由未知量的取值进行了分类探讨,进而求出三阶幻方的全部解。利用线性代数的方法去探求三阶幻方的解法种数显然不是最简捷的,但却有完备的理论做保证。此外,提供了一个如何利用线性代数知识来解决实际问题的实例,并简介了三阶幻方所具有的奇妙的特性。     

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

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

考虑背景孔隙的单开孔两空间结构的风致内压响应研究

风致内压对建筑结构的安全性有着显著的影响,在结构抗风设计时不容忽视。本文基于合理的假定,利用非定常伯努利方程、质量守恒定律和绝热气体状态方程推导了有背景孔隙的单开孔两空间结构的内压响应非线性控制方程组,并通过数值算例分析了背景孔隙对内压响应的影响。结果表明:本文导出的内压控制方程组可很好地分析有背景空隙的单开孔两空间结构的风致内压响应;无背景孔隙时,内侧房间的内压响应要高于外侧房间;背景孔隙所引入的附加阻尼使得两房间的内压响应均受到抑制。这些研究成果对结构风致内压有理论与实际应用的价值。     

过渡状态倾转旋翼气动力模拟的高效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的计算时间。     

不同温度和应变速率下7022铝合金流变应力行为

通过在293-773 K的温度范围内和应变速率为0.001-0.1 s-1下对7022铝合金薄板进行温拉伸试验,研究了7022铝合金温拉伸性能,以及该合金在升温条件下流变应力与变形温度和应变速率之间的关系.并利用改进了的Hooke law和Grosman方程建立了7022铝合金在温拉伸时应力-应变本构模型.研究结果表明:7022铝合金的流变应力随温度的升高而降低,随应变速率的升高而升高;温拉伸试样的延伸率随变形温度的升高而升高,随应变速率的增大而减小.     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

湍流模型对喷流计算精度的影响研究(英文)

采用三维RANS代码MBNS3D对包含了内流、外流及后体干扰的三种简化尾喷流,进行了数值模拟研究。通过壁面压力分布及喷流截面速度型与实验数据的对比,研究了SA和SST两种湍流模型及其修正(SACC,SADES,SSTCC)对喷流计算精度的影响。本文计算表明,SSTCC模拟效果最为精确;利用SA模型进行脱体涡模拟,在同样的网格下优势不明显。湍流模型可压缩修正改善了模拟喷流的膨胀趋势,但对附面层分离位置的预测,在本文的工作中,还不如原始模型。两类模型相比,SST模型的结果好于SA模型。     

梯形翼风洞试验模型数值模拟技术

基于雷诺平均Navier-Stokes (RANS)方程和结构网格技术,采用二阶空间离散精度的MUSCL格式,并结合k-ω剪切应力输运(SST)两方程湍流模型和γ-Re_θ转捩模型,研究了梯形翼风洞试验模型中前缘缝翼、后缘襟翼连接装置对气动特性的影响。简要介绍了本文采用的计算方法;介绍了梯形翼的风洞试验模型及风洞试验结果;在网格收敛性研究的基础上,采用"全湍流"方式和转捩模型研究了梯形翼试验模型连接装置对气动特性的影响。通过与不带连接装置的计算结果的对比,采用"全湍流"模拟方式,计算模型中考虑试验模型的连接装置引起升力系数下降、阻力系数下降、低头力矩减小以及失速迎角提前;通过与试验数据的对比,进一步考虑转捩影响可以提高梯形翼风洞试验模型气动特性的计算结果与试验结果的吻合程度,梯形翼风洞试验模型失速迎角附近的气动特性数值模拟技术还需要进一步的研究。