三角翼大迎角绕流的数值模拟

采用FINETM/HEXA非结构网格流场计算软件包对三角翼飞行器的低速大迎角绕流进行了数值模拟.结果表明,计算值与实验数据符合较好,采用Agglomeration 方法能大幅度提高收敛速度,而网格自适应方法用总体不大的网格单元数可较精细地模拟流场.本文还分析了前缘分离涡破裂前后的流动现象和旋涡横截面流线图谱的变化规律.     

模拟X型号飞机复合材料机翼有限元分析

基于建立的模拟半梁式全金属x型号飞机机翼模型，用等刚度设计方法建立半梁式半复材、全复材的模拟x型号飞机机翼模型。用软件ABAQUS分别采用其中的三角形以及四边形壳单元计算全复合材料机翼、半复合材料机翼、全金属机翼在同种气动外形和载荷下的有限元分析，研究位移和应力的收敛情况。数值计算表明：位移收敛速度较快，应力收敛较慢，用不同类型单元和逐渐加密有限元网格对收敛判断是必要的。半复材机翼比全金属机翼减重35％左右，全复材机翼比全金属机翼减重接近50％。     

气动热数值模拟中的网格相关性及收敛

针对气动热数值模拟中的网格相关性和收敛问题,采用计算流体力学(CFD)计算方法,以典型的钝头为算例,进行了数值模拟研究。研究结果表明:网格是气动热数值模拟中的关键因素,壁面附近法向网格间距最为敏感,热流结果随网格间距不同会出现数倍乃至数量级的差异,在得出正确的壁面压力结果的网格上不能保证得到正确的热流结果。分析认为壁面附近客观存在的巨大温度梯度是网格要求敏感的原因。计算表明:气动热计算的收敛过程比压力和流场收敛慢得多;在压力和流场完全收敛时,热流与最终的迭代收敛值偏差达20%以上;以方程残值的下降量级和压力收敛作为热流收敛的判别标准是不合适的。在气动热计算中,应直接观察热流数据的收敛,确保得到真正的收敛解。     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

不同类型搭接网格对周向槽处理机匣模拟结果影响研究

摘 要: 为深入研究不同类型搭接网格对周向槽处理机匣数值模拟结果的影响,以某单级跨声速压气机及其转子叶顶的周向槽处理机匣作为研究对象,划分了多种计算网格并采用不同类型搭接网格关联周向槽与压气机主流道。详细的三维定常数值模拟和分析表明,传统完全非匹配搭接网格,在计算上会导致流场数据通过搭接网格传递时明显失真,使转子尖部复杂流动呈现显著间断特征,给压气机内流和性能模拟结果带来较大不确定性;增强搭接网格匹配度,有利于提高流场数据传递精度,弱化转子尖部复杂流动的不连续性,降低模拟结果的不确定性;采用完全匹配类型搭接网格,可在计算上完全避免搭界网格上的流场数据传递失真,从而使压气机内流及性能模拟结果具有更高的可信度。     

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

基于工程应用的有限元网格划分研究

首先介绍和分析了划分有限元网格的基本原则,然后基于MSC.Patran有限元软件,结合3个不同的工程实际问题进行网格划分分析,总结了划分有限元网格的一般过程,最后对划分有限元网格进行了小结。     

CFD计算网格误差分析的一个算例

为了给一般工程问题的CFD仿真的网格划分提供参考,研究了光滑直圆管中粘性稳态等温流动在不同尺寸、不同纵横比网格情况下的CFD计算数值解的误差.采用FLUENT软件进行数值计算,应用CFD仿真验证和确认的基本方法,与精确解比较,计算了各种网格的数值计算结果的误差,得到了网格独立解的网格尺寸,给出了误差随网格尺寸变化的规律.     

通道中有稳定器的湍流流场分区计算方法

采用任意曲线坐标系对含有Ｖ型稳定器的通道内湍流流场进行数值计算。对含有Ｖ型槽的复杂流场中网格分区方法作了详细的分析,提出一个新的网格划分方法。数值计算采用曲线坐标系下非交错网格的ＳＩＭＰＬＥ方法,计算时对整个流场进行分区迭代,直到得到收敛结果。文中对一含有Ｖ型稳定器的直通道湍流流场进行了数值计算,结果合理。     

非结构同位网格SIMPLE类算法收敛性能比较

以非结构网格的SIMPLE算法为基础,将算法扩展为SIMPLEC算法.利用30°角的斜方腔流动计算成果,分析了非结构同位网格的SIMPLE/SIMPLEC算法的收敛性能;比较了因网格的非正交而引入的非正交项的取舍对该算法收敛性能的影响;并采用显式校正步法对SIMPLEC算法进行了显式校正.比较表明,在非结构同位网格SIMPLEC算法中可忽略非正交项,但有必要对压力作亚松弛.显式校正步法可显著地加速在非结构网格上求解N-S方程的收敛性能,而且在不同的松弛因子组合下,均有较好的收敛速率.     

预应力作用下弧齿锥齿轮的动频率计算

弧齿锥齿轮是航空发动机中的基本元件,常发生共振破坏.运用自主开发的弧齿锥齿轮设计分析系统建立了包含齿轮完整结构的有限元网格模型,并导人ANSYS软件中进行了考虑工作转速和啮合扭矩引起的预应力影响的弧齿锥齿轮动频率计算,结果表明工作转速引起的离心力和啮合扭矩对弧齿锥齿轮的振动频率有一定的影响.     

一种考虑齿轮副连续啮合过程的接触有限元动力学分析方法

在系统分析齿轮副连续啮合过程不同典型啮合状态特点的基础上,提出了一种可以考虑齿轮副连续啮合过程中啮合齿对变化,受载弹性变形以及滑动摩擦等影响因素的接触有限元分析方法.利用该方法对三种考虑不同因素的齿轮副模型进行了对比分析,并得到了各种因素对齿轮副连续啮合过程动态传递误差及动态接触力特性的影响.研究表明:该方法可以真实模拟齿轮副连续啮合过程中单、双齿对啮合及其交替啮合状态的动力学特性,包括滑动摩擦引起的节点冲击,受载弹性变形引起的啮入、啮出冲击以及时变刚度等激励特性;并可以得到啮入、啮出冲击的大小及作用时间,以及滑动摩擦和齿廓修形对动态啮合特性的影响;动态啮合特性分析结果与齿轮副连续啮合状态的特征完全吻合.      

纺织复合材料细观力学分析的一般性周期性边界条件及其有限元实现

针对纺织复合材料细观有限元分析中单胞网格的快速生成与顺利施加周期性边界条件之间的矛盾,提出了非周期性网格划分条件下,一般性周期性边界条件的数学表达形式.基于ABAQUS有限元软件平台,通过在单胞模型的相对面、相对边及相对角点施加多点约束(MPC)方程,实现了一般性周期性边界条件的施加.结合三维四向编织复合材料单胞模型,对比分析了周期性网格划分和非周期性网格划分情况下,单胞模型受载下的变形状态、应力分布及弹性性能的预测结果,验证了一般性周期性边界条件的正确性和有效性.研究表明:一般性周期性边界条件可以实现复杂细观结构单胞模型的自由网格划分,降低网格划分的难度,提高网格生成的质量.     

基于接触有限元的齿轮副扭转共振分析

以典型斜齿轮副为对象，基于接触有限元法，先模拟其准静态啮合过程并求解啮合刚度，并与国家标准偏差值进行了对比，证明了分析方法的正确性.再模拟其动态啮合过程，通过对啮合力响应进行频谱分析识别出扭转振动的固有频率，并与绘制频响曲线的识别结果一致.在此基础上，数值模拟了该齿轮副扭转共振的发生、发展和稳定过程，定量分析了共振状态下的典型动力学参数响应.结果表明:脱离啮合现象的发生使齿轮副的啮合冲击更为显著，共振状态下的最大应力值相比啮合频率为3000Hz时增大了20%.      

Correlation analysis of PCB and comparison of test-analysis model reduction methods

The validity of correlation analysis between finite element model（FEM） and modal test data is strongly affected by three factors, i.e., quality of excitation and measurement points in modal test,FEM reduction methods, and correlation check techniques. A new criterion based on modified mode participation（MMP） for choosing the best excitation point is presented. Comparison between this new criterion and mode participation（MP） criterion is made by using Case 1 with a simple printed circuit board（PCB）. The result indicates that this new criterion produces better results. In Case 2, 35 measurement points are selected to perform modal test and correlation analysis while 9 selected in Case 3.System equivalent reduction expansion process（SEREP）, modal assurance criteria（MAC）, coordinate modal assurance criteria（CoMAC）, pseudo orthogonality check（POC） and coordinate orthogonality check（CORTHOG） are used to show the error introduced by modal test in Cases 2 and 3. Case 2 shows that additional errors which cannot be identified by using CoMAC can be found by using CORTHOG.In both Cases 2 and 3, Guyan reduction, improved reduced system（IRS） method, SEREP and Hybrid reduction are compared for accuracy and robustness. The results suggest that the quality of the reduction process is problem dependent. However, the IRS method is an improvement over the Guyan reduction, and the Hybrid reduction is an improvement over the SEREP reduction.     

某预警飞机机头ESM天线罩模型的建立与响应分析

建立了某型预警飞机机头ESM天线罩有限元数值模型,对机头ESM天线罩结构进行了静态应力分析,得到了机头ESM天线罩受力情况,同时进行了机载天线罩模态分析以及响应计算。算例表明,罩体与机身对接约束附近区域应力最大,变形较大的区域集中在两肩附近,受撞击时有明显的局部振动,减小两肩附近的蜂窝尺寸可以改善天线罩的稳定性。     

波子结构法在失谐带分流叶片离心叶轮动态特性分析中的应用

为解决高保真失谐叶轮模型计算量大的问题，在Craig-Bampton部件模态综合法的基础上，结合波子结构法对界面自由度进一步减缩，发展了适用于带分流叶片离心叶轮的模型减缩方法，并给出其复数形式下的数学表达。采用该方法建立了工程实际叶轮高保真有限元模型的减缩模型，自由度数减缩率达99.8%。分别采用减缩模型和完整模型，通过计算谐调和不同程度失谐时的模态频率和振型，验证了减缩模型的模态计算精度；通过计算不同程度失谐和不同激励阶次下叶轮的归一化最大失谐位移响应幅值，进一步验证了减缩模型的响应计算精度。数值计算结果表明：在叶片弹性模量失谐标准差15%范围内，采用该方法建立的减缩模型具有较高的模态计算精度并且在各激励阶次下具有较高的响应计算精度，波减缩基对不同程度失谐具有良好的鲁棒性。     

基于能量等效的直齿锥齿轮时变啮合刚度计算

针对直齿圆锥齿轮啮合刚度的计算问题,从微元思想出发,将变截面齿廓划分为若干等截面微段齿段,基于能量等效建立微段齿段啮合刚度计算模型,并利用积分方法获得单齿啮合刚度。此外,基于力平衡和变形协调条件进一步提出了齿轮时变啮合刚度计算模型,同时根据几何关系导出了相应传动误差计算式。利用有限元分析对解析计算模型进行了验证,并分析了误差来源。结果表明:利用该模型不仅能够将直齿锥齿轮啮合刚度计算精度保证在2%以内,还达到了快速求解的目的。      

差曲面拓扑的齿轮啮合刚度计算与承载接触分析

基于齿条-齿轮等切共轭产形原理,构建齿面数值模型、ease-off差齿面,对ease-off蕴含的齿面啮合信息进行解析,获得了齿面接触路径、传动误差、接触线瀑布图;综合ease-off拓扑仿真与轮齿刚度非线性单元耦合解析,给出了修形拓扑齿面的啮合刚度、承载传动误差的计算方法。沿接触路径遍历接触线序列,获得了轮齿时变啮合刚度、承载传动误差与载荷分布图;给出了2阶抛物面对称与对角拓扑两种修形形式算例,求出了系列载荷作用下的啮合刚度、承载传动误差、齿面载荷分布。结果显示:随着载荷的增加,轮齿啮合刚度时变效应明显减弱;承载传动误差波动与啮合刚度、修形梯度密切相关;对角修形在啮合刚度、传动误差、载荷分布特性方面好于对称修形。      

盘一盘螺栓连接结构模态频率分析

基于ANSYSWorkbench软件,采用层单元法和多点约束技术,模拟了盘一盘螺栓连接结构的螺栓预紧力,对其进行了模态分析,并与试验结果进行了比较。研究表明：采用层单元法模拟螺栓预紧力的模态结果与试验结果吻合较好,为在预紧力作用下螺栓连接结构模态频率的有限元计算提供了1种有效方法。