不同支承条件下弯曲振动梁的频率方程及特征值

本文推导出两端带有扭转弹簧的铰支梁的弯曲振动频率方程,并由此导出一些特殊支承条件下的梁的振动频率方程,并用数值方法计算出不同支承条件下梁的特征值,本文还通过绘制变化曲线的方法,直观地说明了特征值随梁的刚度EI,长度L和弹簧弹性系数K三个参数的变化及相互影响的情况。     

悬停状态倾转旋翼噪声试验及数值计算

结合消声室试验和数值计算研究了悬停状态倾转旋翼气动噪声特性。消声室试验采用孤立倾转旋翼模型,测试了不同总距角和桨尖马赫数状态的气动噪声数据。数值计算基于CFD(computational fluid dynamic)结合FW-H(Ffowcs Williams-Hawkings)声学方程的方法,采用试验数据验证了计算模型的可靠性。分析了倾转旋翼的噪声传播特点以及拉力系数和桨尖马赫数对噪声总声压级的影响,并对比了孤立旋翼和双旋翼状态的气动噪声特性。结果表明:倾转旋翼噪声随着拉力系数和桨尖马赫数的增加均有所增加,维持旋翼拉力不变时降低桨尖马赫数虽然使得拉力系数增加,旋翼噪声水平仍然降低;倾转双旋翼噪声相对纵向平面对称分布,在多个方位角区域存在着局部最大值,这和双旋翼噪声传播时的相互叠加以及双旋翼间气动干扰相关。      

低重环境下圆柱贮箱内液体晃动特性研究的一种解析法

考虑低重环境下圆柱贮箱内由于表面张力影响而呈现弯曲自由液面的情况,利用Fourier Bessel级数对贮箱受横向激励时的自由液面处的边界条件进行展开,得到液体晃动系统的广义状态方程,给出了固有频率、晃动波高、晃动力和晃动力矩等晃动特性的计算公式。通过数值算例与文献结果对比,在验证文中方法正确性的同时,具体研究了各晃动特性随充液深度、外激励频率和Bond数等参数变化的规律。     

带副翼偏转的三角翼自由滚转运动数值模拟

 通过耦合求解非定常Euler/Navier-Stokes方程和单自由度滚转运动方程,对带副翼偏转的65°后掠角尖前缘三角翼WI1-SLE自由滚转运动进行了研究,Navier-Stokes方程的求解采用基于Spalart-Allmaras湍流模型的脱体涡模拟(DES)。在多块结构网格上,应用基于弧长的无限插值理论(TFI)生成变形网格,实现副翼偏转,而三角翼的滚转运动则通过网格的整体旋转实现。结果表明:Euler方程和DES方法均准确地模拟出了三角翼在滚转运动过程中存在的3个平衡位置。出现平衡位置的原因分别是:①流动对称性;②机翼左侧发生涡破裂的分离涡与右侧分离涡相互平衡使得滚转力矩为0,并且平衡位置仅与三角翼两侧涡强的差有关;③副翼偏转和左右机翼不对称分离涡涡强差产生的滚转力矩相互平衡。此外,滚转运动对副翼偏角幅值很敏感,幅值的微小改变会影响最终的平衡位置和向平衡位置运动的路径。     

基于气密载荷作用的飞机舱门可靠性分析

在阐明飞机舱门系统启动工作原理的基础上,重点考虑舱门所需克服气密载荷产生的密封带摩擦力的大小及其分散性的影响,构建了舱门意外打开及启动的可靠性分析模型,提出了相应的基于气密载荷作用的飞机舱门可靠性分析方法。以某型飞机舱门机构为例,作了具体的算例阐明,分析方法与结果可为舱门设计提供参考。     

一组新的四元数轨道要素建模方法

针对传统以欧拉角为参数的轨道要素的奇异性、不确定性以及计算效率低等问题,提出了一种新的四元数轨道要素。建立了新轨道要素与经典轨道要素,以及新轨道要素与惯性系下位置、速度的相互转换关系,推导了基于新轨道要素的高斯摄动方程。以J₂项摄动下的轨道推演为例进行仿真验证,结果表明新轨道要素不仅在圆轨道与赤道平面轨道处不存在奇异性和不确定性,而且由于新轨道要素不涉及三角函数运算,新高斯摄动方程积分效率和计算精度明显提高。     

有限差分法的坐标变换诱导误差

有限差分法应用于具有复杂外形的网格时需要进行坐标变换,在此过程中经常会引入坐标变换诱导误差。在柱坐标系下使用均匀网格进行均匀流场计算,计算结果表明,即使物理坐标对计算坐标的变换函数连续可导、计算过程中坐标变换系数直接采用准确的解析式、网格完全正交并且充分光滑,也无法避免坐标变换诱导误差。理论分析表明,产生坐标变换诱导误差的机理是笛卡尔坐标系下的守恒型欧拉方程变换至贴体坐标系下后增加了源项。针对该问题,目前国内外学者通常采用几何守恒律,构建与差分格式相匹配的坐标变换系数计算方法来消除源项。本文介绍了从包含源项的离散等价方程基础上直接进行离散的新算法,在此基础上针对非等距网格条件下MUSCL类格式重构过程进行误差分析,理论推导表明重构中需要考虑非等距插值公式的影响系数,将变量转换至计算空间内进行MUSCL重构才能保证该过程具有均匀网格下的插值精度。通过理论分析及数值实验证明新算法对于均匀流场完全不会引入坐标变换误差。     

接收模式下加罩天线系统特性的快速电磁仿真

针对接收模式下电大尺寸的天线-天线罩系统电磁特性分析问题,提出了一种全波法和高频法结合的混合算法:积分方程/修正型表面积分+多层快速多极子(IE/MSI+MLFMA)方法.该方法将由高频方法MSI得到的天线罩内场分布作为天线阵列的激励场,再利用基于体-面混合积分方程(VSIE)和MLFMA的快速全波分析方法精确计算天线阵;同时应用球谐函数展开、预处理技术、混合并行等技术进一步改进算法的计算效率.同传统的全波数值法相比,该方法在保证计算精度、满足工程需要的前提下,解决了求解计算时间长、计算效率低的问题,实现了对电大尺寸天线-天线罩系统的快速、高效仿真.      

CFD数学模型的线性化方法及其应用

计算流体力学(CFD)方法不仅仅起到数值模拟的作用,它本身是一个复杂的非线性系统。在流动稳定性分析、气动弹性分析、优化设计以及流动控制等领域,从系统的角度出发,对CFD数学模型线性化后,可以对模型的系统矩阵进行定量分析获得更多的系统特性。但是CFD数学模型往往非常复杂且阶数很高,因此其线性化系统矩阵的获得比较困难。鉴于此,采用人工编程和自动微分相结合,构造有限体积法并行CFD模型的线性化系统矩阵。其中自动微分只被用来得到每个界面通量的局部雅可比矩阵,而采用人工编程方法来实现并行环境下的稀疏雅可比矩阵的组装。线性化系统的并行求解采用了块雅可比预处理的广义最小残量法,每个并行进程内部则采用零填充不完全LU分解预处理。为了验证这种线性化方法,上述方法被用于:① NACA 0012翼型的非定常绕流线性系统构造与求解;② NACA 0012翼型稳态流动的伴随方程构造与求解;③ AGARD wing 445.6机翼颤振问题降阶建模。上述三个算例的结果与CFD模拟的吻合一致。