解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

一种通用的冲击分析模型及低速冲击响应计算

 提出了一种通用的冲击分析模型,该模型采用连续介质理论描述,用拖带坐标定义摩擦冲击约束,用虚功原理定义冲击控制方程,可以作为不可刺穿型冲击问题的有限元求解基础。在此模型基础上,针对低速冲击问题,给出了采用自适应罚函数增广 Lagrange迭代处理接触力的有限元求解方法,既保证计算精度,又有效抑制算法震荡,并用算例进行了验证。     

基于短时谱分析的SISAR运动参数估计及侧影像定标

针对解目标运动参数、多普勒频率和角度之间的非线性方程运算量大且极易陷入局部极值的问题,提出一种估计阴影逆合成孔径雷达(SISAR)目标运动参数的新方法.该方法通过短时谱分析与抽取叠代相结合,提高了估计精度,减少了运算量.在对目标侧影像进行定标时,利用恢复的侧影轮廓中线相位的差分均值消除目标相对于雷达的非线性转动引起的附加相位,获得几乎不受目标运动参数影响的侧影轮廓中线的高度和长度.仿真验证了所提方法的可行性.     

有限元线性方程组改进的全稀疏存贮方法

有限元线性方程组的系数矩阵一般具有稀疏性和对称性的特点,全稀疏存贮方法就是利用这些特点,只存贮对称部分的非零元素,采用链表式管理,即节省存贮空间,又便于动态更改;但在实际求解计算时,查找元素花费的时间较多;这里给出了一种改进的方法,是在链表的基础上生成相应的线性表,能够方便快速查询元素,使存贮结构同时具有链表和线性表的优点.数值算例证明该存贮方案在时间和存贮上都较为占优,可靠高效,能够应用于有限元线性方程组的求解.     

用于失谐叶盘动力学特性分析的减缩计算方法

提出了一种用于失谐叶盘动力学特性分析的减缩计算方法.该方法采用基于周期对称的谐波平衡法,建立非线性的迭代方程,能够基于单扇区的矩阵求解整个失谐叶盘的强迫响应,且不存在任何截断误差,因而效率高于部件模态综合法但精度并不降低.从周期对称结构理论出发给出了方法原理的数学推导,在此基础上通过一个失谐叶盘有限元模型算例对该方法的效率与精度进行了讨论.结果表明:相较部件模态综合法,周期对称谐波平衡法矩阵维度减小59.5%,计算复杂度减少93.4%,计算时间节省57.4%.     

A Direct Solution to GPS-Type Navigation Equations

One solution to the navigation equations involves iteration on the 4 by 4 augmented range-direction-cosine matrix beginning with an assumed position and so assumed direction cosines, of which there are 12 for 4 satellites. An algebraic, direct solution to this same basic equation set has recently been published. Both of these methods are reviewed. We offer a direct solution using modified functions of the range magnitude data from four satellites to yield user's clock bias correction, user's position, and true range vectors if desired. The highest order of matrix inversion used is 2 by 2. The highest order, nonlinear equation is a numeric square root. The principle of the formulation is use of differences among the range magnitudes and range magnitudes squared. An additional auxiliary difference equation is formed. A computation basis uses the ephimeride differences and an orthogonal vector. The method offers convenience, speed, simplicity, low dimensionality, and precision, with no operational constraints.     

航空发动机非线性模型实时计算的迭代方法研究

针对航空发动机非线性模型将作为未来机载应用而需要严格实时性的问题,研究了求解模型的迭代方法,综述其历史发展情况,并讨论收敛性问题,同时提出一种改进应用方案,即模型在过渡态仿真时可将上一状态点雅可比矩阵直接应用于下一点计算,能够在典型应用方案基础上直接减少气动热力模型计算次数,从而提高整体模型执行效率.通过实际算例对这些应用方案进行仿真比较与分析,结果表明:①Broy-den秩1法执行效率最高,更适宜用于模型实时计算;②改进应用方案简单、有效,为进一步满足非线性模型的严格实时性要求奠定基础.      

基于双耦合Riccati方程的混合鲁棒飞行控制器

针对混合控制器设计中耦合Riccati方程难于得到解析解的问题,采用一种有效、快速的李雅普诺夫迭代法,并将此算法应用到飞行控制器的设计中,抑制了飞行控制系统中大量的非线性和不确定性因素对飞机的干扰.仿真结果表明,所设计的控制器不仅使闭环系统满足鲁棒稳定性,还满足一定的动态性能,从而证明了它的有效性和实用性.     

基于实验复模态参数的有限元模型修正

提出了一种利用非完备实验复模态参数,修正带有非比例阻尼矩阵结构有限元模型的迭代修正方法。首先,对处理困难的阻尼矩阵表达为质量与刚度矩阵的多项式组合,并给出了物理解释,被修正参数为该多项式的系数和结构的某些局部几何、材料参数,其修正结构物理意义明确,与结构的动态设计相适应,实例表明了本文修正方法的有效性。     

求解涡扇发动机数学模型的有限域搜索方法

提出一种基于N-R迭代方法的有限域优化搜索法,用于求解涡扇发动机数学模型。该方法不仅具有N-R优化迭代求解特色,还具有将初猜解附近划定域进行泛搜能力,使迭代运算具有强收敛性和平稳性。该方法所定义的迭代因子可通过编程自动调整,进而可以确定最小迭代次数及此意义下的最优迭代因子。      

模型变形视频测量的相机位置坐标与姿态角确定

高速风洞试验中的模型变形视频测量（VMD）要求双（多）相机大角度大重叠的测量方式,而通过现有共线方程的线性化模型难以取得高精度的相机位置与姿态角,故推导包含共线方程泰勒展开二次项的非线性误差模型,建立3控制点的VMD相机位置与姿态确定技术。多个工程实例表明该技术能取得高精度的相机位置坐标与姿态角,有实用价值。     

欧拉方程隐式迭代矢通量分裂跨音流场计算

 本文用矢通量分裂格式的有限体积法求解时间相关欧拉方程组。文中证明,矢通量分裂的有限体积法不能再采用物面反射边界条件。二维问题用隐式迭代法解代数方程组,第三维方向上则用强隐格式往返扫描。此法CFL数可达100的数量级,收敛快,占内存少。跨音流场计算结果与实验符合良好,激波只占一两个网格。     

高速发动机油膜惯性对活塞裙润滑影响

利用一种新的迭代算法 ,分析了高速发动机油膜惯性对活塞裙润滑特性的影响 .该迭代方法是利用有限元法和差分法交替求解Navier Stocks方程和雷诺方程 ,据此导出了适用于高速发动机活塞裙的润滑计算的混和润滑雷诺模型 .新的模型借助惯性系数 ,引入了油膜惯性项 ;同时给出求解含有油膜惯性项的迭代步骤和有限元表达式。计算结果表明 :随着惯性系数和活塞裙的长径比的不同 ,油膜惯性会对油膜的摩擦力、压力和承载力产生不同的影响 ,这种影响对承载力尤为明显 .该润滑模型也可用于中、低速发动机的活塞裙润滑计算以及不计入惯性项 (惯性系数置为零 )的某些润滑问题求解 .     

Kepler方程的一种迭代加速算法

采用Steffensen迭代法求解椭圆轨道下的Kepler方程,克服了不动点迭代法收敛速度慢以及Newton迭代法需要重复计算导数值和迭代结构相对复杂等方面的不足。试验表明利用该方法求解Kepler方程既满足了高精度要求,而且收敛速度快。同时它还能满足不同的偏心率和不同的精度的要求,因而是求解Kepler方程的一种较理想算法。     

直升机结构承载能力研究

介绍用迭代方法求解直升机结构承载能力的过程,求解非线性方程。计算时,首先假定结构是给定的,求出结构最大承载能力;以某一直升机模拟试验件为例,算出各试验件每一元素的应力、位移、屈服情况、失稳部位随着加载过程而变化情况以及最大承载能力等。计算结果与试验结果吻合很好。     

考虑库仑阻尼时叶盘系统振动响应的求解方法

本文研究了具有干摩擦阻尼的带冠叶盘系统振动响应的求解问题。解决了时频法中由于叠混和泄漏而使迭代发散的问题。最后用时频法获得了令人满意的响应计算结果。     

桥联重积分方程一种新的数值方法

纤维增强复合材料断裂力学中的桥联问题常常可归结为一个二重积分方程。常见的求解方法有矩阵迭代法和积分迭代法，但这两种方法都产复杂。本文给出一种新的分段插值方法，对此积分方程进行求解。数值计算结果表明，该方法效果令人满意。     

一种简化算法的航空发动机全状态数学模型

通过一种简化的燃气空间模型对某型双转子混排加力涡扇发动机各部件进行了建模和仿真计算,并给出了详细的计算流程以及数学公式.在仿真过程中结合实际工程应用经验,对涡轮处采用一个局部的小循环迭代以减小涡轮部件内燃气空间计算的时间步长,提高计算精度.与其他数值方法相比,该简化算法避免了求解非线性方程组时大范围的循环迭代,提高了全状态数学模型的计算速度.通过仿真结果与发动机实际试车数据的对比,证明该方法有足够的计算精度.      

Multiple target angle tracking using sensor array outputs

The use of the output of an array of sensors to track multiple independently moving targets is reported. The output of each sensor in the array is the sum of signals received from each of the targets. The results of direction-of-arrival estimation by eigenvalue analysis are extended to derive a recursive procedure based on a matrix quadratic equation. The solution of this matrix quadratic equation is used to provide updated target positions. A linear approximation method for estimating the solution of the matrix equation is presented. The algorithm is demonstrated by the simulated tracking of two targets. The main advantage of the algorithm is that a closed-form solution for updating the target angle estimates has been obtained. Also, its application is straightforward, and the data association problem due to uncertainty in the origin of the measurements is avoided. However, it requires the inversion of an N×N as well as other linear operations, so that the computational burden becomes substantial as N becomes very large     

全变差正则化方法对一类热传导反问题的应用

热传导终端反问题是典型的不适定问题,应用基于全变差的正则化方法来获得它的稳定数值解,将近似解的求解范围由连续函数空间扩大到了有界变差函数空间,并将固定点迭代法用于求解非线性的Euler方程,并进行了理论分析和数值实验。结果表明:该方法与算法具有高效稳定全局收敛等优点。