滑动式Thiele型连分式插值方法在GPS精密星历中的应用

在动态GPS精密定位中,必须使用高精度的GPS卫星轨道数据,但是IGS组织只提供15min间隔的精密星历,无法满足间隔时间较短的动态定位要求,用多项式逼近效果很差,另外对于不同的星历插值没有高效的方法能确定最佳多项式阶数。因此,利用Thiele型连分式建立有理函数,并在此基础上提出滑动式Thiele型连分式插值的方法,简化了方法又提高了内插精度,并通过算例与Lagrange多项式和Chebyshev多项式进行了分析和比较,结果表明该插值方法可以更加有效地改进插值精度。     

卫星位置和速度的Lagrange插值算法分析

导航接收机中通常采用开普勒轨道参数法计算卫星的三维位置和速度.这种算法计算量较大,在资源有限的廉价用户机或要求高频度输出卫星位置和速度的高动态接收机中会使处理器资源非常紧张.采用Lagrange插值算法可以获得相同精度的卫星位置和速度,且大大降低计算量.利用GPS数据对不同插值区间长度和插值阶数的Lagrange插值法的计算精度和处理时延进行了定量分析.结果表明,Lagrange插值算法计算卫星位置和速度可以使处理时延降低十几倍,这在每个历元的导航解算需要同时计算多颗卫星位置和速度的接收机中具有重要意义.     

不带稳定翼飞船返回舱俯仰动稳定性研究

联盟号飞船返回舱,如果不存在稳定翼,从高超声速到低Mach数时,力矩曲线有由一个平衡迎角(α1),经二个(a1、α'2)到三个(α1、α2、α3)平衡迎角的情况出现.当有一个平衡迎角时,可能有一临界Mach数Mcr,当M∞＞Mcr时,平衡迎角处是动稳定的;M∞＜Mcr,变为动不稳定的,同时出现极限环.也可能不存在这个Mcr,这时随Mach数降低到M∞=M'cr,即2个平衡迎角的情况,在力矩曲线与α轴相切的平衡点α'2,会出现鞍、结点分叉.M∞＜M'cr后,对于存在三个平衡点的情况,在α和α的相图上,它是结点-鞍点-结点的结构.当α1处是稳定结点形态时,只有在很小的扰动情况下,它才能返回α1的状态.如果α1为不稳定结点状态,在α1附近,微小的扰动,使迎角由α1迅速变成α3或者更大.这是混沌的先兆.     

水下超声速气体射流线性稳定性研究

水下超音速气体射流稳定性是水下推进、水下焊接技术研究的焦点领域。通过建立线化小扰动气液混合流体控制方程,开展水下超音速气体射流线性稳定性研究。控制方程的形式表明射流基本流分布(速度与浓度分布)、雷诺数、无量纲密度、无量纲粘度能够影响水下超音速气体射流稳定性。使用Chebyshev配置点法对控制方程进行求解。计算结果表明射流基本流分布以及无量纲密度是影响射流稳定性的主要因素。射流基本流分布越平缓,射流最大扰动增长率越小;无量纲密度越小,即气液密度差越小,射流最大扰动增长率越小。随雷诺数增加,射流最大扰动增长率先减小后增加,当雷诺数大于105时,射流最大扰动增长率趋于常数。对不同无量纲粘度,射流扰动增长率-特征波数曲线高度重合。同时分析了射流不同截面的稳定性特征,计算结果表明射流最大扰动增长率随距喷口距离增加而减小。这些结论有助于理解水下超音速射流物理过程,为实际应用提供理论指导。     

KdV方程的数值解及数值模拟

文章在已构造的一种分片光滑的Lagrange型插值多项式空间中利用Galerkin法对非线性KdV方程进行求解。分析了数值格式的线性稳定性并进行了数值模拟,结果表明格式精确有效。     

台阶流模型流场的局部DQ-Lagrange求解方法

以Lagrange多项式作为Differential Quadrature(DQ)方法的基函数,应用局部DQ-Lagrange方法对Navier-Stokes方程进行了数值求解,建立了基于局部迎风DQ-Lagrange方法的复杂台阶流流场求解模型.并在此研究基础上,提出了本质和自然边界条件的处理方法.分析了支持域的选择、边界条件处理方法、迎风机制等对求解精度的影响.结果表明:支持域较大时,Lagrange高阶插值易出现振荡现象,可以采用局部替代全局的插值方法;迎风支持域的加入可以提高求解精度;建模时速度支持域的选择应该大于压力支持域.      

构造含结点旋转单元的坐标变换法

Ａｌｌｍａｎ提出了带角结点旋转的三角形单元及相应的四边形单元，法向位移的二次插值保证了该单元的精度，但带来了构成单元的复杂。本文把Ａｌｌｍａｎ单元看成是对一个有边中结点单元进行坐标变换的结果，大大提高了计算效率。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

一种改进的Dolph—Chebyshev方法在阵列天线方向图综合中的应用

通过对传统Chebyshev多项式增加修正因子得到一种修正的Chebyshev多项式，并以此为基础提出了一种新的阵列天线方向图综合方法。与传统Chebyshev综合法相比，该方法可以在保持副瓣电平不变的情况下对主瓣进行任意展宽，待求参数少，综合过程简单，有很强的实用性。     

集中涡的细微结构和破裂机理初探

本文利用微分方程定性理论和分叉的基本概念,证明了在轴向逆压梯度时集中涡截面流态是一个附着螺旋点,且存在至少一个极限环。破裂截面上的流态对应于一个退化结点。由此推出破裂时有V_θ/r|_(r=o)=0。通过分析认为,涡破裂是由于涡轴上轴向逆压梯度的扰动自动放大造成的。     

阶梯函数逼近阶与若干函数类的特征

利用阶梯函数逼近给出了连续函数类的特征,并在此基础上分别导出阶梯函数逼近李普希兹函数类和狄尼-李普希兹函数类的逼近阶。     

Stress-strain state analysis of shells of revolution based on the triangular discretization element with the Lagrange multipliers

The present paper sets forth an algorithm of analyzing thin shells with significant gradients of meridian curvature based on the Kirchhoff-Love hypothesis and finite element method. As a discretization element use is made of a triangular fragment of the middle surface with the Lagrange multipliers in additional nodes arranged in the triangle side middles. We apply the vector method of displacement interpolation that has a number of advantages over the conventional interpolation procedure. The high efficiency of the algorithm developed is confirmed by the numerical examples.     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

一种考虑截面任意变形模式梁的有限元模型

提出了一种考虑截面完整变形的截面插值梁模型,并用于机翼结构的建模与分析。首先引入截面插值函数——拉格朗日函数描述截面形状,以位移向量为未知变量描述截面位移,在此基础上依据插值理论构造梁单元位移场。不同于传统梁单元通过假定的中性轴的挠度和转角来确定梁截面各点位移,该梁模型摒弃了中性轴假设与平截面假设,通过截面插值函数得到梁截面面内、面外变形;然后通过有限元理论推导了梁单元刚度矩阵与质量矩阵;最后利用截面插值梁单元对机翼各部件进行有限元建模,并展开典型工况下的静力学分析以及动力学分析,与Nastran实体单元计算结果进行了对比分析,验证了该梁模型的有效性,为机翼的结构设计和强度分析提供了一种新的简化建模方法。     

斜齿轮齿面柔度矩阵与修形的有限元计算

提出了一套斜齿轮修形用的有限元计算程序。该程序自动生成轮齿有限元网格和齿面接触线, 计算得到接触线上结点的柔度矩阵;然后考虑支承结构与齿轮的综合变形, 由数学规划获得最佳齿面修形。计算与实验结果对比证明该程序操作简便、计算精度高, 具有工程应用价值。      

非线性热传导方程混合问题插值逼近

以Legendre-Gauss-type积分节点为插值节点,构造插值基函数展开数值解,逼近有界杆上的非线性热传导方程Dirichlet边界条件的正确解。给出算法格式和相应的数值算例,表明所提算法格式的有效性和高精度。所给算法适合于非线性问题求解。     

多轴联动条件下插补速度实时可调的前瞻控制算法

针对目前多轴联动条件下S形加减速算法复杂,一条插补线段内最大进给速度不能实时可调的问题,提出了一种插补速度实时可调的前瞻控制算法。算法首先根据机床特性和线段夹角建立了小线段衔接处进给速度的约束条件,然后采用加速度跟随原理提出了最大进给速度连续可调的S形加减速控制方法,在此基础上设计了一种加速度自适应调整的前瞻处理算法。该算法在不降低轨迹插补精度的前提下,能以最大的速度通过线段的转接点,从而使整个线段插补过程具有高度的柔性和快速性,能满足现代数控系统对前瞻处理的实时性要求。结果表明该算法有效降低了机床运动时的振动,与传统速度规划算法相比,同等加工条件下,加工效率明显提高,工件表面质量也得到改善。     

Direct trajectory optimization based on a mapped Chebyshev pseudospectral method

 In view of generating optimal trajectories of Bolza problems, standard Chebyshev pseudospectral (PS) method makes the points' accumulation near the extremities and rarefaction of nodes close to the center of interval, which causes an ill-condition of differentiation matrix and an oscillation of the optimal solution. For improvement upon the difficulties, a mapped Chebyshev pseudospectral method is proposed. A conformal map is applied to Chebyshev points to move the points closer to equidistant nodes. Condition number and spectral radius of differentiation matrices from both methods are presented to show the improvement. Furthermore, the modification keeps the Chebyshev pseudospectral method's advantage, the spectral convergence rate. Based on three numerical examples, a comparison of the execution time, convergence and accuracy is presented among the standard Chebyshev pseudospectral method, other collocation methods and the proposed one. In one example, the error of results from mapped Chebyshev pseudospectral method is reduced to 5% of that from standard Chebyshev pseudospectral method.     

学科间信息传递的参数化空间散乱数据插值法

耦合信息传递是实现耦合系统求解的关键技术之一,本课题研究了参数空间散乱数据插值方法.为了避免三维插值空间平整性导致的插值误差,将三维耦合界面进行参数化处理,在参数空间中构造插值函数进行插值传递.针对径流式叶片的特点,进行了温度插值传递.从插值效果看,这种基于参数空间的数据插值方法保证了较高的插值精度.比较了反距离加权平...