基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法,建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题,旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法,通过添加冗余约束,避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块,对基于运动嵌套网格 的空间流场进行求解,湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块,基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解,验证了网格变形模块和CFD模块的有效性,然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上,计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷,并与飞行测试数据进行了对比。计算结果表明,文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析,提高了旋翼气弹载荷的预测精度。     

一种基于余代数单子的Web服务形式化模型

为拓展服务计算的形式化研究视野、手段和方法,建立并实现了一种针对Web服务的服务计算形式化模型。鉴于开放环境下的服务实体主要来源于不同的第三方提供者,将软件实体抽象成余代数单子,从而以一种黑盒方式给出软件服务的语义模型。给出了余代数单子的一般性定义,在此基础上对软件服务进行单子描述,进而提出一种基于余代数方法和单子技术的Web服务参考模型。最后,实现了一个基于单子的Web服务计算平台原型系统,可支持从遗留系统中进行服务抽取、发布、发现和度量等。     

0.8m 风洞风扇系统的设计和运转

为提高低速风洞的最大风速,根本的措施是设计高转速大功率的风扇系统。风扇系统的设计首先在于有一个良好的设计方法,选择高升阻比的风扇翼型,根据需要恰当地确定设计工况点,合理选择桨叶和反扭导流片的各个参数,如桨毂比、实度、升力系数和圆弧角等。设计中针对已经设计的外形,进行性能计算,可以迅速鉴别设计的优劣,检查是否能满足设计要求,节约研制成本和缩短周期。     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

用腐蚀损伤计算金属日历寿命的原理和模型

通过对金属腐蚀T-H曲线和温湿谱的腐蚀损伤研究,找到一种用腐蚀损伤计算金属日历寿命的原理和模型。这种计算金属日历寿命的模型,只要给出真实使用环境下的腐蚀T-H曲线和温湿谱,就可计算出金属真实使用日历寿命。通过对30CrMnSiA钢的日历寿命实例计算结果与真实环境下的腐蚀试验结果的比较得到,相对误差是17.5%,这说明该日历寿命计算模型和方法是有效的。     

金属任意腐蚀损伤量的日历寿命计算模型和曲线

国际机械日历寿命研究,都是在假设已知腐蚀损伤容限Dc的前提下,求在这个Dc下的金属日历寿命。为了便于应用,本文不再假设腐蚀损伤容限Dc为已知数,而是把Dc当做变量,求金属任意腐蚀损伤量Di的日历寿命。本文通过多方面深入研究,发现腐蚀试验溶液的浓度dt与试验时间Ht的乘积除以试验获得的腐蚀损伤Dt等于一个常数。由此给出一种简单易行的金属任意腐蚀损伤量日历寿命的计算模型和相应的求解曲线。由此模型或曲线,可求得金属任意腐蚀损伤量的日历寿命。     

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

弧齿锥齿轮齿面拓扑修形及加工参数计算

为了能够实现对齿面啮合性能的灵活控制,针对弧齿锥齿轮小轮提出一种齿面拓扑修形方法,即借助二阶曲面对齿面偏差拓扑的近似表达,将齿面拓扑修形分解为5个方向:螺旋角修正、压力角修正、齿长曲率修正、齿廓曲率修正及齿面挠率修正,通过改变5个方向的修形系数对小轮齿面拓扑结构进行自由控制。在此基础上,建立齿面偏差与机床加工参数之间的修正数学模型,通过构建敏感性矩阵并采用最小二乘法求解,反求出获得修形齿面的小轮加工参数,以便指导加工。以一对弧齿锥齿轮副为例进行修形啮合分析,仿真结果表明:选取齿长曲率修形系数为0.0001,齿廓曲率修形系数为0.0005,齿面挠率修形系数为0.0003,对齿面进行拓扑修形后传动误差幅值为-25.60″,接触迹线倾斜角度变为54.7°,相比原始结果啮合性能得到改善。滚检接触区与理论仿真结果一致,验证了修形方法的有效性。      

基于图论的燃油管路稳态流动计算方法

基于节点分析法,提出将局部损失并入相应管段的处理方法,推导阻力计算公式,从而得到兼顾沿程和局部损失的稳态流动计算方法。将该计算方法应用于实际的航空发动机燃油管路,采用VB.NET语言编写计算软件;同时,对燃油管路开展三维数值模拟计算。结果表明:局部损失对管路的流量分配影响较大;该计算方法得到的流量分配与数值计算结果趋势吻合很好,差别在±10%。