一种多点定位的目标位置精确解算方法

多点定位是民航飞机导航系统的新技术,精确的目标位置解算方法是多点定位的关键.通过设定参考站和变换,将非线性到达时间(TOA)方程组转化为线性到在时间差(TDOA)方程组,提出了一种两步求解目标位置的闭式算法.首先忽略TDOA测量误差,获得了目标位置的粗解,在粗解的基础上做泰勒级数展开,克服了未知量不独立对精度的影响,再...     

三类非线性代数方程系统解的存在性

将大量的应用问题归纳成3类非线性代数方程组解的存在性问题.罗列了目前获得的主要结果,提出了一些待研究的问题.     

湍流的通用物理方程

利用变形张量及摄动变形位移的点积分解湍流脉动速度,并使用密度加权的侧偏系综平均方法取代 Reynolds平均方法,在不引入任何物理假设的情况下,对 Navier-Stokes方程进行平均。未知相关量在三阶精度上才出现,截断之后,获得具有二阶精度的封闭的湍流方程组。该方程组无相关量,不引入任何经验系数,无需构造任何湍流模型,即可描述 Ma<3.6流速的任意复杂湍流场并解释湍流的物理本质。本文仅对推导思路及方程组做概要介绍。      

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

多钉连接件钉传载荷计算的一个解析方法

对多钉连接件钉传载荷的计算问题提出了一个解析分析方法,推导了求解钉载的线性代数方程组并给出了若干算例。该方法使用方便,并具有工程分析所要求的精度。     

明渠水流三维数值模型研究

本文在具有无尺度铅垂坐标的直角坐标系中采用完全的三维雷诺方程组描述明渠中的水流运动,利用二维沿水深平均的连续方程确定自由表面形状。文中给出的隐式求解方法保持了经典SIMPLE算法的全部优点,并具有高的收敛速度。算例表明,计算结果与实测值十分一致。     

Navier—Stokes（NS）方程组差分计算中的物理和网格尺度效应及NS方程组的简化

对于高Re数流动计算，在通常二阶精度NS差分格式和网格数条件下，存在某些粘性项落入修正微分方程截断误差项的问题。这类NS方程组计算实际是计算某种简化NS方程组，而且重复计算误差物理粘性项既浪费机时和内存，误差积累又会对数值解产生不可预测的影响，避免外述缺陷的办法一个提高NS差分格式的精度，另一个是丢掉可能落入截断误差项的物理粘性项，把NS方程组简化为广义NS方程组，广义NS计算避免了误差物理粘性项误差积累对数值解的不可知影响，又可节省内存和机时，对高Re数流体工程计算很有好处。利用广义NS方程组计算超声速绕前向和后向台阶流动的结果表明：广义NS方程组与NS方程组的数值结果很好相符符。     

混合初边值问题解的渐近性质及其差分解

 本文讨论了当t→∞时拟线性方程组初边值解的一些渐近性质。在假设非定常问题的解唯一、定常问题的解存在时,对同一边界条件,证明了当t→∞时非定常方程组的解趋于定常问题的解,讨论了用差分方法求解这种问题时,边界点上一些量的精确算法;为了消除在物体尖角附近出现的局部非物理解,得到正确解,我们还给出了一个简单的耗散型物面边界格式;最后给出了一些尖角凹形体的算例。     

基于多小波的Volterra级数非定常气动力建模方法

非定常气动力建模除了要准确描述气动力的非定常特性,还要反映其非线性特性。Volterra级数因为对系统非线性具有很强的描述能力正日益受到重视。一阶Volterra核只能表达线性特性,要建立反映非线性特性的非定常气动力模型,需要引入二阶核甚至更高阶核的影响。高阶Volterra核辨识的主要困难在于待辨识参数的数量随着核阶次的增加而呈指数增加导致计算难度急剧加大,即出现所谓维数灾难问题。以一种分段二次多小波为基函数将Volterra核展开,求解一个高维病态方程组来计算展开系数,利用小波的多分辨分析在时间和频率两个维度的分解特性将方程降维,最终将问题转化为求解一个低维方程组得到稳定解。通过NACA0012翼型在马赫数0.8下作沉浮运动时,升力系数、阻力系数和俯仰力矩系数的二阶核和三阶核的辨识构建非定常气动力模型,然后由此计算不同减缩频率下的气动力并与CFD结果进行比较,验证了Volterra级数对非线性非定常气动力的描述能力和多小波处理方法的有效性。     

二维非结构网格上的高精度有限体积WENO格式

构造了非结构网格上二维双曲型守恒律的一类新的高精度有限体积WENO格式。其主要思想是:根据格式精度的要求,按照谱体积方法对三角形单元网格进行剖分,通过选取适当的子单元组成模板,利用WENO重构方法重构二阶和三阶多项式,利用有限体积公式和高阶Runge-Kutta TVD时间离散方法,构造了非结构网格上二维双曲型守恒律的一致二阶和三阶精度的有限体积WENO格式。然后,推广到二维Euler方程组。最后,给出几个数值算例,验证了格式的稳定性、高阶精度和高分辨捕捉激波等间断的能力。     

基于灵敏度分析的结构动力模型修改

本文基于灵敏度分析和Taylor级数展开式提出了一种参数型结构模型修改方法。为了改进模型修改精度,本文从能量的观点出发引入一权矩阵,并给出模型修改方程的加权最小二乘解。数字算例表明本文方法具有很好的精度,并且在利用结构模态试验数据进行修改方面具有较好的适应性。     

线性有延迟的微分方程的解

在解决一些有关生物数学的问题时，经常会遇到求解具有延迟的微分方程。到目前为止，还没有非常有效的方法求解这类方程，只是对于一小部分具有特殊形式可以求解。介绍一种此类方程的解法，其主要思想是把微分方程转化为一般的代数方程求解，然后研究解的稳定性，技巧性很强。     

解决欧拉方程奇异性的方法探讨

提出了一种解决欧拉方程奇异性的工程算法,从飞行器运动方程入手,推导出具体计算公式,并进行了数值计算验证。该方法无论是理论分析还是工程计算都证明是可行的,没有发现原理缺陷和方法误差,不会产生奇异性,在飞行器仿真、模拟和姿态控制中具有普遍的应用和参考价值。     

大展弦比飞机运动稳定性建模及分析

随着飞机设计技术的提高和新型材料的应用，现代飞机气动弹性效应越来越显著，刚体运动和弹性体振动互相解耦的关系不再明显，在分析弹性飞机相关问题时，应该综合地、统一地考虑飞行动力学与气动弹性力学问题。本文在一般体轴系下综合考虑了飞机刚体运动自由度和弹性振动自由度，建立弹性飞机刚弹耦合运动的方程，通过在配平状态下引入小扰动运动假设，以及利用有理函数拟合技术将偶极子格网法计算所得到的频域非定常气动力转化为时域形式，建立了大柔性飞机刚弹耦合小扰动状态空间方程，并对纵向运动稳定性进行分析计算。针对文中构型飞机，运用本文所采用的方法计算所得结果与传统的线性及刚体计算所得结果进行对比分析，非线性计算方法所得颤振速度更小，而且失稳模态亦发生变化，对于飞行力学模态而言，长周期模态稳定性变差，在计算速度范围内出现失稳的现象，由此说明机翼大变形对飞机稳定性所带来的影响。     

用杂交应力元及等参元计算应力强度因子

 本文介绍一种用奇异性杂交应力元及等参元计算平面裂纹的应力强度因子和J积分的方法。实际计算表明,此种方法的优点是精度高,计算时间短,便于实际应用。文中给出三个算例。各例分别采用216～256个自由度。计算结果同解析解及其他奇异元的解进行了比较。各例算得的应力强度因子值,均同解析解相差仅1～2％。每例中由4～6条路径求得的J积分值,其分散度均在±O.7％以内。 计算在TQ-16机上进行,程序编译约6～7分钟,每例的计算时间约6～8分钟。     

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

A Discrete Method for Parameter Identification in Linear Systems with Transport Lags

Linear dynamical systems with transport lags are characterized by linear differential-difference equations. The task of identifying unknown parameters in such systems from the input-output data is difficult due to mathematical complications associated with differentialdifference equations. This paper presents a method which solves the identification problem. The method is digitally oriented and shows how a continuous-time system can be identified by discrete techniques. The solution is based on Kalman's least square method. The identification procedure essentially involves two steps: 1) discretizing the continuous system via finite difference approximation, and 2) estimating the parameters through the identification of the resulting discrete model. Experimental results have verified the validity of the proposed method.     

几种生长模型参数的最小二乘估计方法及应用

提出一种Logistic模型、Gompertz模型、Mitscherlich模型和Bertallanffy模型的分布参数估计方法，这种方法以加权最小二乘法为基础，推导求解分布参数的方程组。利用这类曲线模型中参数形式的特点，使同时求解三个参数的问题简化为同时求解两个参数的问题，并给出求解参数耦合方程组的二重二分迭代法。预测地基沉降的计算实例表明，这种方法对于迭代初值的选取要求较低、求解速度快、精度高，而且通过选取权因子可以对已有数据实现“重近轻远”的处理。     

On Identifying Transfer Functions and State Equations for Linear Systems

Two methods are established for identifying constantcoefficient C2n-type noise-free linear systems if the time response data of the input-output or of all states are known. 2n response data are required to identify an nth-order transfer function or state equation for an unknown linear system. The order of the unknown system can be identified by checking a sequence of determinants. The Z transform and its inversion are mainly used.     

复合材料叠层板壳有限元的渐近解

 本文针对复合材料叠层板壳结构建立了一个分层计算的有限元模型,此模型的特点是以参考面的位移和各层横向剪应变作为独立的自变量,导出了叠层板壳的有限元模式,这是和原有分层模型的本质区别,从而具有以下优点:(1)刚度矩阵具有良好特性,可采用渐近法求解,大大提高了计算效率,克服了分层模型中的最大困难;(2)完全避免了较薄的板壳在有限元计算中常发生的“剪切自锁”现象;(3)较精确地考虑了各层的剪切效应,无须再引入剪切修正系数。还采用此模型较好地解决了任意铺层的叠层板弯曲、叠层板脱层屈曲和承受任意载荷的轴对称叠层壳的弯曲等问题,证明了此模型的优点和可靠性。此模型还极易推广应用于叠层板壳的动力、稳定和层间应力等问题。