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

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

电弧喷射推力器流动区域电磁场数值解法

为了有效地进行电弧喷射推力器流动区域全场(流场、电场、磁场和化学反应)耦合数值模拟并揭示其电磁场特征，建立了适当的电磁场模型，对电磁场控制方程的数值解法进行了详细研究。模型基于推力器稳态工作特征，控制方程由麦克斯韦方程组简化得到，采用有限控制容积积分方法离散椭圆型控制方程，采用9种不同的迭代方法求解离散方程。给出了不同迭代方法的收敛速度、数值稳定性和最终精度。研究表明，Gauss—Seidel逐线超松弛迭代法是推力器流动区域电磁场离散方程的一种快速有效的数值解法。     

超燃冲压发动机流量匹配的临界流量法

为了使用基于积分方程的超燃冲压发动机性能计算模型模拟跨燃、亚燃工作模态,用临界流量法建立用于表征燃烧室热壅塞程度的残差方程,并以隔离段出口静压为独立变量,通过牛顿迭代法迭代求解残差方程,模拟了燃烧室内存在临界截面（即喉道）但不出现热壅塞的实际跨燃、亚燃工作状态。结果表明：该方法能够确定有喉道的超燃冲压发动机燃烧室一维参数分布,并具有计算精度高、计算速度快、收敛性好等优点,可以为基于积分方程的超燃冲压发动机性能计算模型求解提供一定的参考。     

弹性飞机空气动力及动态特性的数值计算和研究

本文给出了一种弹性飞机动态特性的数值计算和分析方法。文中考虑了飞机受扰后弹性振动与飞机总体运动的耦合作用,并着重研究了前几阶低频结构振型对飞机总体运动的影响。本方法采用了格林定理所导出的速势积分方程方法计算定常和完全非定常空气动力。采用特征根迭代法求解弹性运动方程,只需经较少次数的迭代即可收敛。文中附有算例,其计算结果是合理的。     

二维对流—扩散方程的控制体有限元法数值模拟

推导了二维对流—扩散方程的控制体有限元公式，其主要特点是所得离散化方程适合于逐点迭代法及共轭梯度法进行求解，使得该方法更适合于求解复杂的几何域及节省更多的计算机内存。通过几个算例表明，该方法对于特别大的Ｐｅｃｌｅｔ数其精度也是很高的     

基于离散伴随方程的三维雷达散射截面几何敏感度计算

针对有限差分法计算雷达散射截面（RCS）梯度效率低，采用高精度雷达散射截面评估时计算代价高的问题，提出了一种基于麦克斯韦积分方程离散伴随方程的RCS梯度高效计算方法。基于伴随方程的梯度计算可以通过一次雷达散射截面求解、一次伴随方程求解获得RCS关于所有设计变量的梯度。其中麦克斯韦积分方程离散伴随方程的形式与原方程基本一致，可以采用矩量法（MOM）及多层快速多极子算法（MLFMA）求解。伴随方程求解计算量与直接雷达散射截面评估基本一致，存储量在直接雷达散射截面评估的基础上增加不明显。通过双椎体模型、导弹模型对基于矩量法、多层快速多极子算法的伴随梯度进行验证，证明了基于伴随方法的RCS梯度计算可以实现复杂外形中RCS梯度的高效、高精度求解，为基于梯度的高精度气动/隐身一体化优化提供了基础。     

翼型跨音速绕流的数值计算

本文选用精确速势方程作为翼型跨音速无粘绕流的数学模型。在变换过的直角座标中,用有限差分(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用列松弛迭代法求解。 为了计入粘性效应,利用精确速势方程和附面层动量积分方程进行联合迭代求解。这一点对超临界翼型的计算显得特别重要。 算例与其它数值解及实验作了比较,对于一般翼型和超临界翼型,本文结果是良好的。     

一种求解任意形状导体目标内谐振时表面电流的有效方法

用矩量法求解闭合导体目标的电场积分方程或磁场积分方程将得到不正确的表面电流。提出了一种确定正确表面电流的有效方法，即正确的表面电流由非谐振模电流和归一化的谐振模电流与一未知复因子的乘积组成，通过令闭合导体内给定点的总电场为零可得该未知复因子。计算了一无限长理想导体圆柱谐振时的表面电流，所得结果与解析解一致，进而证明了本方法的有效性和准确性。     

数值微分的小波最小二乘方法

给出了利用小波求解数值微分的方法.将微分问题转化为等价的第一类Fredholm积分方程,然后选取小波基,使用最小二乘求解该积分方程.我们给出的分析框架容许使用具紧支集的小波,利用正则化保证了计算的稳定性.数值计算实例表明,本方法稳定,运算速度快,计算精度高.     

机翼跨声速绕流的压强极小积分有限元算法及应用

本文从全速位方程出发。利用压强极小积分,给出了机翼跨声速绕流的有限元算法,并编制出相应的FORTRAN语言程序。文中给出了适于有限元法的网格生成技术,采用了线松弛迭代法和逐步加密法进行数值求解,提供了M6机翼和一个三角翼的计算结果。     

纤维增强复合材料币状I型裂纹桥联问题

利用平面等效方法，将轴对称币状裂纹问题转化平面裂纹问题，由此给出复合材料内币状Ⅰ型裂纹的桥联积分控制方程。对权函数法应用于求解币状裂纹问题的可行性提出质疑，并给出修正的建议。然后对币状Ⅰ型裂纹桥联问题给予了断裂力学分析。利用一种新的分段插值方法求解了桥联积分控制方程。针对该积分程提出了一种近似分析解法，其结果与数值结果符合得相当好。     

计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

边界元法分析回转体的自由扭振问题

 边界元法用于分析非线性问题和动力问题时,在最终得到的积分方程中同时包含边界积分和区域积分。通常的做法是对区域积分项采用常值有限元离散,这需要用很密的网格才能达到较高的精度。     

一种基于自适应粘滞Landweber算法的雷达杂波模拟方法

在Weibull分布雷达杂波模拟的模型中,针对隐函数迭代法求解无记忆非线性变换法(Zero Memory Nonlinear,ZMNL)中相关参数效率低下的问题,对Landweber迭代算法进行了优化,使之具备粘滞迭代格式和自适应步长,进而将相关系数求解问题转化为第一类算子方程问题。仿真结果表明,改进的算法收敛速度明显提升,显著提高了Weibull雷达杂波的模拟精度。     

弹塑性接触问题边界元缩聚法分析

给出了利用隐含对称条件的二维弹塑性边界元分析公式,并应用于弹塑性接触问题分析。根据接触条件列出了弹塑性接触问题迭代求解的耦合方程。利用接触问题局部非线性的特点,将缩聚法用于迭代求解。基于拟线性、双重迭代法的思想,解决了两重非线性相耦合的问题。最后给出了一个算例。     

一个n阶微分方程的两点边值问题

有一些n阶非线性微分方程的两点边值问题已被人们讨论过,大多是在微分方程的右端函数有界、满足李普希兹条件或在相应边值问题存在上下解的情况下讨论的。[3]、[4]在这篇文章里,讨论了一个n阶非线性微分方程的两点线性边值问题,是首先通过n-2次累次积分将原方程化成了一个二阶微积分方程边值问题,并用拓扑度理论讨论了解的存在性,同时给出了解的唯一性条件。推广了文献[1]的结果。     

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

Solving a radiation problem with forward speed using a lifting surface method with a Green''s function

The progress of seakeeping computations requires development of computating codes for unsteady flows around a ship or its elements. In this paper, we present a method of calculation concerning waves radiated by an oscillating surface-piercing flat plate with forward speed, with a yaw angle. By use of Green's third identity, the problem is transformed into the resolution of a Fredholm integral equation of the first kind by a panel method using Green's function. The Green's diffraction-radiation function with forward speed is used. Its numerical values are calculated by an adaptative integration procedure to reduce the computation time. The present method permits determination of the pressure jump distribution across the plate, the total forces and moments. The results obtained are compared with other numerical methods in hydrodynamics and in aerodynamics, and with experimental data obtained in a water tank.     

Optimum Signal and Filter Design in Underwater Acoustic Echo Ranging Systems

Optimization of the filter, the signal, and the signal and filter jointly are studied in the sonar environment under noise and reverberation limited conditions. The maximization of the receiver output signal-to-interference ratio is used as a performance criterion with unit energy constraint on both signal and filter. In the filter design problem, the optimum filter function is the solution of a linear integral equation. The kernel of the integral equation is a function of the target and medium scattering functions and the reverberation distribution. In the signal design problem, a similar type of integral equation is obtained as in the filter optimization problem. In the joint signal and filter design problem, it is shown that the optimum signal and filter functions are the solutions to a pair of linear integral equations with the largest (SIR)O. Several examples are investigated for different mediums and reverberation distributions with the finite matrix approximation method. An interative technique is used to compute the joint optimization of signal and filter.