高精度有限体积法与间断有限元法的比较

通过数值算例,比较了高精度有限体积法和间断有限元法在求解不同问题时的表现。研究发现:在精度相同的条件下,间断有限元法的计算误差要明显小于有限体积法;间断有限元法的重构过程与高精度有限体积法相比较为简单,但高阶情形下解多项式的自由度较多并且需要计算体积分,因此整个求解时间较长。降低时间积分时解多项式的自由度数目是实现高精度算法在实际问题中应用的重要手段。     

基于迎风型有限体积法的非结构化网格流场格心数据高质量体绘制方法

在计算流体力学领域,对复杂流场进行数值模拟时常用非结构化网格划分空间,并广泛采用迎风型有限体积法求解流场.该方法将数值解存储在非结构化网格单元(控制体)中心,产生格心数据,并假设物理量在相邻单元邻接面上形成"间断".然而,在可视化此类数据时,目前的重构方法破坏了间断约束,导致"流间断"信息丢失,严重制约了绘制质量.为此...     

应用Groebner基方法求解代数方程组的解

在字典序下计算方程组的多项式生成的理想的Groebner基G,根据Groebner基G中单变元多项式的解,依次递推求出多项式方程组的解.通过求解多项式函数条件极值问题的一般步骤和具体实例说明该方法的计算过程.     

高精度有限体积格式在三维曲线坐标系下的应用

为了构造在光滑区具有较高分辨率并且可以无振荡捕捉激波的高精度有限体积格式,同时降低格式在模板选择上遇到的困难,基于逐维重构方法,发展了结构网格下的高精度有限体积格式,并将这一格式推广到三维曲线坐标系下,从而可以适应相对复杂外形下的计算。为充分验证格式的有效性,选取一系列典型算例进行计算:在等熵涡输运和二维喷管流动中验证了格式的精度可以达到设计精度;在双马赫反射问题中格式也表现出良好的捕捉激波的能力。数值计算表明,上述格式在曲线网格上具有较高的数值精度和鲁棒的激波捕捉能力,适用于流体力学方程的计算。     

运动载体涡流磁场数值计算与特性分析

为提高运动载体涡流磁干扰建模与补偿的精度,对载体运动引起的涡流磁场的产生机理和分布特性进行分析研究。首先,对涡流磁场的产生机理进行了溯源分析;然后,采用有限元法对涡流磁场进行数值计算研究。考虑到计算空间包含结构的磁导率单一问题,采用了有限元节点法矢量磁位计算方式,并针对舱体壁计算单元,考虑了耦合积分电势自由度,同时考虑到磁场周期变化,采用了瞬态有限元计算方式;最后,通过数值计算分析,揭示了运动载体涡流磁场与外磁场变化频率的关系,并分析了涡流磁场沿径向和轴向的空间分布特性和传播规律,为运动载体涡流磁场干扰的高精度建模与补偿奠定了基础。     

基于表面电荷法的像管静电场计算

与传统的有限差分法和有限元法相比,表面电荷法可降低单元剖分难度并能处理开放边界,更适于计算静电透镜的静电场。针对轴对称及非轴对称两类结构,采用样条函数及奇异形状函数改进电极表面的电荷密度估计,讨论了二维及三维表面电荷法的实现。与具有解析解的单位圆盘进行比较,分析了表面电荷法的电位计算精度,并用于轴对称及非轴对称像管计算。结果表明：对轴对称系统,在文中给定步长下,二维及三维表面电荷法与有限差分法计算相对误差不超过10^-3,且二维方法精度整体高于三维方法,二维、三维表面电荷法对3种像面位置的计算结果与有限差分法结果相比,最大误差不超过0.6mm;对非轴对称系统,3种方案下轴外电子运行轨迹的差异明显,非对称因素所产生的干涉场影响是不可忽略的。     

非局部薄板动力学特性分析的广义有限积分变换方法

基于应力梯度非局部薄板理论模型,推导了非局部薄板动力学特性求解的广义有限积分变换方法.通过选取适应边界条件的积分核函数并构建广义积分变换对,应用积分变换将非局部薄板的高阶偏微分方程变换成线性方程组,直接求解得到固有频率.将广义有限积分变换方法的计算结果和有限元法及已有文献的结果进行对比,验证了本文方法的正确性.在此基础...     

线性弹道模型建立与仿真

对6自由度刚体弹道模型进行了线性化,得到了线性弹道方程组。求解线性弹道方程组得到其解析解,并给出详细求解过程。以6自由度刚体弹道模型为基准,在不同射角条件下用6自由度刚体弹道模型和线性弹道模型进行了弹道仿真。仿真结果表明:该线性弹道模型能满足工程精度要求,它模型简单,弹道计算方便快速,在弹道滤波处理中可作为卡尔曼滤波模型,也可作为弹道快速模型用于简易制导弹箭预测控制系统设计。     

一种求解飞行器廓线间断处超音速复杂流场的方法

用有限体积法求解三维完全气体的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,并将通量守恒的分区计算方法和Ｂ－Ｌ湍流模型修正结合在一起,捕捉到了间断处的分离涡、激波、膨胀波等流场信息     

一种GNSS卫星轨道高精度实时插值方法

针对GNSS导航信号模拟源中多星座轨道计算实时性问题,从计算效率、精度等角度分析和对比了牛顿多项式和埃尔米特多项式插值算法的性能,通过对牛顿插值的加“窗”改进,使得卫星位置、速度的插值精度显著提高.算例结果表明,该方法能够在等间距时间点上多个卫星位置已知、速度未知的条件下,高精度实时地内插出卫星的位置和速度,其中卫星位置插值精度为mm级,速度插值精度小于1e-5 m/s,计算量为广播星历直接计算的1/7,位置和速度的拟合精度均优于三阶埃尔米特插值结果.     

基于Gauss伪谱法和直接打靶法结合的月球定点着陆轨道优化

将一种求解最优控制问题的新方法—高斯伪谱法( Gauss Pseudospectral Method-GPM)和传统的直接打靶法有效结合,对月球着陆器定点软着陆轨道快速优化问题做出了研究.推导了高精度模型下着陆动力学方程.针对优化方法各自的特点和多约束条件下最优月球软着陆轨道设计的难点,提出了问题求解的串行优化策略:将控制变量和终端时间一同作为优化变量,同时离散控制变量与状态变量,取较少的Gauss节点,利用GPM求解初值,初值的求解采用从可行解到最优解的串行优化策略;在Gauss节点上离散控制变量,利用直接打靶法求解精确最优解.仿真结果表明,本文提出的轨道优化方法具有较强的鲁棒性和快速收敛性.     

双曲线法抗跟踪干扰研究

针对中、低速跳频电台跟踪干扰问题,从通信抗干扰的角度出发,通过对干扰椭圆方程的变形分析,提出了抗跟踪干扰的双曲线法,推导了抗干扰双曲线的方程,并计算出了跟踪干扰无效区域。最后,分析了该双曲线的特点及影响因素,介绍了该方法的运用。与现有的抗跟踪干扰理论依据相比,根据双曲线法选择抗干扰措施更直观、更简单。     

失效概率函数求解的高效算法

结合基于分数矩约束的极大熵方法和替代模型法,发展了一种失效概率函数求解的高效算法。所提算法的基本思路是利用自主学习的迭代Kriging方法来构造失效概率函数,即采用较少的训练样本来构造粗糙的失效概率函数,在此基础上通过添加新的违反学习函数约束的样本来更新失效概率函数,直到达到精度要求。对于每一个分布参数的训练样本点,所提方法采用分数矩约束的极大熵法来求解相应的失效概率样本。由于分数矩的计算采用了高效的降维积分,并且由于分数矩约束下极大熵法中优化策略高效地逼近了响应的概率密度函数,从而使得失效概率样本能够被高效高精度地估计出来。为了检验所提方法的精度及效率,给出了两个算例,对比了所提方法与已有的失效概率函数求解的Bayes公式法及Monte Carlo法等,结果表明,所提方法适用于求解复杂的功能函数问题,且在满足精度要求的基础上大大降低了计算量。     

基于备件需求预测的两步法理论?

在对备件需求时间序列研究的基础上,结合指数平滑法和Croston法的特点,分析了指数平滑法与两步法的原理,通过对指数平滑法和两步法方差的研究,得到两步法是指数平滑法的一般形式的结论,为两步法的进一步研究提供一定的理论支撑。     

基于广义标准轨迹的平衡滑翔状态反馈制导方法

对多约束条件下远程助推滑翔飞行器再人滑翔飞行问题,提出了一种基于广义标准轨迹的平衡滑翔状态反馈制导方法.建立了远程助推滑翔飞行器的动力学模型,确定了飞行轨迹约束条件,详细阐述了基于广义标准轨迹的平衡滑翔状态反馈制导方法的制导原理,设计了远程助推滑翔飞行器的侧向和纵向制导律,并采用LQR( linear quadratic regular)方法设计了纵向制导参数,仿真验证了该方法的可行性.与以往再人滑翔制导方法不同,该制导方法主要利用飞行攻角的变化来调节飞行轨迹,飞行过程中飞行器的速度倾侧角较小.仿真结果表明,该制导方法能满足远程助推滑翔飞行器的再入滑翔制导问题,并且具有较好的鲁棒性和自适应性.     

OGY法控制混沌研究

OGY法是一种有效的混沌控制方法,其利用混沌运动对微小扰动极端敏感及遍历特性,通过系统参数的微小摄动将混沌控制到期望的轨道.文中对OGY法及其两种改进方法作了介绍,进行数值仿真,比较3种方法的差异,为今后的混沌控制研究方法指出了新的方向.     

螺旋桨梢涡卷起数值计算

通过解析公式计算得到面源和面偶的诱导速度,在面偶附近区域,用涡环代替面偶,计算得到螺旋桨的速度场,根据尾涡面必须和当地流体流速相切的原则,修正原来尾涡面的形状,逐步迭代直至尾涡形状收敛,用这种方法预报了三维水翼和螺旋桨的尾涡面的形状和梢涡卷起,与试验结果有较好的一致性。     

用面元法预报吊舱推进器水动力性能

面元法的各个环节有不同的数值处理方法,不同的方法计算结果有所不同,通过对各种方法的计算结果进行比较,得到适合螺旋桨计算的数值处理方法。用面元法对吊舱螺旋桨的水动力性能进行了计算,得到了吊舱推进器的敞水曲线和轴承力。     

一种齿轮箱振动控制方法

从机械结构的振动控制理论出发 ,提出了一种齿轮箱振动控制方法 ,研究并讨论几种齿轮箱箱体以及齿轮本体的结构参数对动态特性的影响 ,其结论可以为齿轮箱的振动设计提供参考依据