基于工程应变的轴对称塑性大变形几何方程及其应用

从工程线应变和工程切应变的物理概念出发，根据“轴对称”和“大变形”的特点，推导出相应的几何方程，另导出大变形条件下的体积不变定律。与现有“有限变形”条件下的几何方程相比，不但更为精确，而且其表达形式在具体应用中更为简单。在一些典型的工艺问题中，将小变形条件下导出的几何方程应用于大变形情况也可能得出非常好的结果，不但运算简单方便，而且其精度还优于“有限变形”条件下的几何方程。     

基于弹载捷联惯性导航系统精确导航的双欧拉全姿态方法

垂直发射的地空导弹在大机动飞行时,其四元数法姿态角输出存在奇异性。为解决该问题,提出了正、反欧拉角微分方程切换的双欧拉全姿态解算方法。对正、反欧拉角微分方程进行了详细推导,论述了正、反欧拉角微分方程精华区与奇异区的倒置关系。在此基础上,建立了正、反欧拉角的转换公式。根据正欧拉角的俯仰角对正、反欧拉角无缝切换,有效解决了俯仰角在±90°附近时姿态奇异性问题。仿真和半物理试验的结果表明:双欧拉全姿态算法解算的实际姿态与理论姿态保持一致,在奇异点附近的姿态角仍能正常输出,未出现突变现象,该方法得出的导航精度比四元数法得出的导航精度提高了9.81%。     

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

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

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

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

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

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

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

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

Nonlinear electron-acoustic rogue waves in electron-beam plasma system with non-thermal hot electrons

The properties of nonlinear electron-acoustic rogue waves have been investigated in an unmagnetized collisionless four-component plasma system consisting of a cold electron fluid, non-thermal hot electrons obeying a non-thermal distribution, an electron beam and stationary ions. It is found that the basic set of fluid equations is reduced to a nonlinear Schrodinger equation. The dependence of rogue wave profiles on the electron beam and energetic population parameter are discussed. The results of the present investigation may be applicable in auroral zone plasma.     

一种提高空间实验室定轨精度的方法

以近地航天器轨道动力学为基础,建立变阻力系数大气摄动模型,设计了求解变阻力系数的算法。然后利用天宫一号飞行器的测轨数据进行计算,分析了空间实验室飞行高度的轨道特性,其中包括：大气模式密度误差、变阻力系数与空间环境关系、定轨残差和星历误差。在空间环境平静和磁暴的条件下,制定了多种求解变阻力系数的策略,解决了空间实验室长弧段定轨精度受限的问题,并在空间环境平静条件下实现了优于10米的定轨精度,在磁暴条件下实现了优于20米的定轨精度。     

叠片式箔片气体动压推力轴承的静动态特性

针对一种叠片式箔片气体动压推力轴承,提出相应的结构刚度模型。运用有限差分法(FDM)求解可压缩气体雷诺方程,并与箔片弹性变形方程耦合求解得到该类型轴承的静动态特性。结果表明:由于支承结构在径向和周向都具有变刚度特点,导致单瓣箔片上出现二次收敛楔形效应,因此形成两个压力峰值,与传统类型轴承相比能在一定程度上提高承载能力。随着最小气膜厚度的增加,轴承承载能力和摩擦力矩都不断下降,且动态刚度也逐渐下降,但对动态阻尼的影响较小。存在最优楔形高度值以实现较大承载能力的同时,也能达到较小摩擦力矩,减小功率损耗,并且还能获得较高的动态刚度和动态阻尼系数。