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

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

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

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

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

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

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

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

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

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

管道热边界层方程的迎风有限元分析

对于管道热边界层方程,除了采用动量积分方法求得理论解析解外,也可以用数值方法求解,如有限差分、有限体积、有限元等方法.理论解析解是采用一定的简化并忽略若干项之后得到的,因此,也只是一种近似解,数值解可以考虑完整的方程和各种边界条件,因而其解较为全面.采用伽辽金有限元方法求解,管道热边界层方程为标准的对流扩散方程,当对流项较强时,需要采用迎风方法,因而也给出了迎风有限元方法的模型.     

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

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

梁、板非线性动力响应分析的样条有限点法

本文用样条有限点法分析梁和基于Berger假设的板的几何非线性动力响应,推导了以样条结点参数表示的基本方程,用wilson-θ法求解动力增量方程,求得挠度响应曲线。本文的方法与有限元法比较,共计算时间省,程序简单,可在微机上计算。     

二维非线性正交各向异性材料的瞬态热传导反问题数值方法

针对正交各向异性材料的二维非线性热传导反问题,采用顺序函数法进行表面热流辨识问题的研究。在求解反问题时,采用有限体积法、牛顿-拉夫逊法并引入未来时间步的概念。在每个时间步内,将待辨识热流视为非线性方程组的未知量,通过一个迭代过程进行求解。算例的研究表明,热流辨识结果与真实热流相近,从而证明了本方法在辨识二维非线性热传导反问题时是准确、稳定、有效的。     

截顶圆锥薄壳振动特性的计算

用有限单元法计算,单元为轴对称截锥壳。法向位移函数w 和周向位移函数v 采用三次多项式,纵向位移函数u 为线性多项式。其节圆自由度选取节圆上三个位移,一个转角和一个剪应变。计算结果与文献[1]、[2]和[3]作了比较。这方法和程序可以应用于计算截顶圆锥壳的振动特性。     

双曲线法抗跟踪干扰研究

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

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

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

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

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

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

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

OGY法控制混沌研究

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

螺旋桨梢涡卷起数值计算

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

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

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

一种齿轮箱振动控制方法

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