Dirac矩阵和Pauli矩阵在塑性力学中的应用

本文在求解塑性力学问题中,采用了量子电动力学中著名的Dirac矩阵和Pauli矩阵,使求解塑性应力增量的问题变得十分简单.     

一种全速域的计算方法及其应用

针对原可压缩流动求解器不能用于低速不可压缩流动预测的缺点,采用预处理技术对控制方程特征系统、隐式求解方法进行修正,并采用预处理修正的AUSM+-up格式离散对流项.采用修正后的求解器对无粘鼓包流动、顶盖驱动粘性方腔流动以及Laval(拉瓦尔)喷管流动等算例进行数值仿真,并将数值仿真结果与基准解进行对比.结果表明将预处理技术应用于全速域流动问题的求解是可行的,经预处理修正后的求解器能够用于低速、亚音速、跨音速以及超音速流动问题的求解.     

关于拟线性抛物型方程的加速收敛分析

在工程实际问题中,求解非线性问题的有限元近似解,通常采用低次元,得到的解精度较低。为了达到较高的精度,就必须采用高次元,但这给实际计算带来了很大的困难,为此,文[1]作者提出了一种求解非线性问题的加速计算方法。首先用低次元求出非线性问题的有限元近似解,再利用这个近似解把非线性问题的求解化成相应的线性问     

“运价”为负的运输模型计算机求解问题的处理

为解决一些计算机软件求解"运价"既有正值又有负值运输模型时"不可求解"的问题,本文采用"运价同额增减法"决策模型转换的方法,将原模型的"运价"全部转换为正值后再用计算机软件求解,并分别编写了EXCEL求解模板和求解程序对该方法的计算加以印证。结果表明,采用该方法求解得出的最优解(最优决策方案)与原模型求得的最优解完全一样,而最优值(最优决策效果)减去虚增(或加上虚减)的部分就是原模型的最优值。采用这种方法能成功地解决一些计算机软件"不可求解"的问题。     

带平衡约束的离散网络平衡设计问题的遗传算法

本文研究了带平衡约束的离散网络设计问题及其求解算法.模型中上层是一个离散网络设计的数学规划模型,采用遗传算法来求解.下层是采用变分不等式描述的用户平衡配流问题,利用对角化方法直接求解.通过实例对算法进行验证,结果表明该算法是有效的.     

基于人工蜂群算法的配送中心选址问题求解

采用人工蜂群算法对配送中心选址问题进行求解,给出食物源的编码方法,通过整数规范化,使算法能在整数空间内对问题进行求解.应用算法进行了仿真实验,并将结果与其它一些启发式算法进行了比较和分析.计算结果表明人工蜂群算法可以有效求解配送中心选址问题,同时也为算法求解其它一些组合优化问题提供了有益思路.     

基于正弦型光滑打磨函数对0-1规划问题的连续化求解方法

传统的求解0-1规划问题方法大多属于直接离散的解法.现提出一个包含严格转换和近似逼近三个步骤的连续化解法:(1)借助阶跃函数把0-1离散变量转化为[0,1]区间上的连续变量;(2)对目标函数采用逼近折中阶跃函数近光滑打磨函数,约束条件采用线性打磨函数逼近折中阶跃函数,把0-1规划问题由离散问题转化为连续优化模型;(3)利用高阶光滑的解法求解优化模型.该方法打破了特定求解方法仅适用于特定类型0-1规划问题惯例,使求解0-1规划问题的方法更加一般化.在具体求解时,采用正弦型光滑打磨函数来逼近折中阶跃函数,计算效果很好.     

一类高阶牛顿迭代法及其在线性互补问题中的应用

通过等价转换,把线性互补问题转化为一个不可微的非线性方程组,进而采用光滑函数处理,得到一个光滑非线性方程组,利用高阶牛顿迭代法进行求解.该方法不再区分线性互补问题是否单调,因此扩大了线性互补问题的求解范围.计算结果表明,方法计算速度快,对线性互补问题求解较为有效.     

凸二次函数优化问题在快速多极边界元法中的应用

利用快速多极边界元法(FMM-BEM)求解大规模工程问题最终结为稀疏线性方程组的求解,因此,采用更好的方法求解线性方程组可以提高边界元法的计算效率.本文利用最优化数值技术处理,将稀疏线性方程组的求解等价为求解一个凸二次函数极小化的问题,并利用最优化理论及相关数学理论证明了其解的存在唯一性,为该理论的形成和发展奠定了理论基础.     

边界节点法计算二维瞬态热传导问题

采用边界节点法(BKM)结合双重互易法(DRM)求解二维瞬态热传导问题.采用差分格式处理时间变量,可将原瞬态热传导方程转化为一系列非齐次修正的Helmholtz方程.随后,方程的解可分为特解和齐次解两部分计算,引入双重互易法在区域内部配点求解方程的特解,采用边界节点法仅需边界配点求解方程的齐次解.给出的数值算例显示该方法计算精度高,适用性好,具有很好的稳定性和收敛性,适合求解瞬态热传导问题.     

研究伪随机数序列独立性的均匀偏度法

§1 引言 由具有在[0,1]上均匀分布的总体中产生的简单子样称为随机数序列,其中的每一样本称为随机数。所谓伪随机数序列,一般是指用数学递推公式所产生的随机数序列。如最常见的乘同余方法,它所产生的伪随机数序列就是对给定的正整数M和小于M的任意正整数初值β_1,用如下递推公式确定的:     

一种基于矩阵初等变换的Schmidt正交化方法

对向量组的Schmidt正交化法和合同变换法的关系进行了分析,指出Schmidt正交化法就是合同变换法中利用规范化初等变换后的一种特殊情况,由此给出一种基于矩阵初等变换的Schmidt正交化方法——Schmidt初等变换正交化法,以及这一方法在软件Matlab上实现的程序.     

自适应多重网格法与超松弛法的比较

多重网格法(Multiple Grid Method,简称M-G方法)是近年来出现的快速方法之一,本文在M-G方法中采用自适应控制层间转换的技术,并将自适应M-G方法与G-S迭代方法及SOR迭代方,法进行了比较。其计算结果表明,自适应M-G方法的计算量比G-S迭代及SOR迭代少得多,当M-G方法所用层数为4-6层,这种优越性就更加明显,且自适应M-G方法中选取控制参数有很大的灵活性。     

Finding roots of a real polynomial simultaneously by means of Bairstow's method

Aberth's method for finding the roots of a polynomial was shown to be robust. However, complex arithmetic is needed in this method even if the polynomial is real, because it starts with complex initial approximations. A novel method is proposed for real polynomials that does not require any complex arithmetic within iterations. It is based on the observation that Aberth's method is a systematic use of Newton's method. The analogous technique is then applied to Bairstow's procedure in the proposed method. As a result, the method needs half the computations per iteration than Aberth's method. Numerical experiments showed that the new method exhibited a competitive overall performance for the test polynomials.     

一种求解非线性互补问题的方法及其收敛性

本文将Newton方法和外梯度方法相结合,提出了一种求解非线性互补问题的方法,证明了此方法的全局收敛性和超线性收敛性,在适当的条件下给出了一个有限终止结果。数值实验表明,此方法是有效的。     

对称正定矩阵的多级迭代法

本文主要研究解对称正定矩阵的多级迭代法,并对其收敛性进行证明。最后用数值实验验证此方法的有效性．多级迭代法特别适用于并行计算,并且可以被理解为古典迭代法的扩展,或共轭梯度法的预处理子。     

Accelerated Imaginary-time Evolution Methods for the Computation of Solitary Waves

Two accelerated imaginary-time evolution methods are proposed for the computation of solitary waves in arbitrary spatial dimensions. For the first method (with traditional power normalization), the convergence conditions as well as conditions for optimal accelerations are derived. In addition, it is shown that for nodeless solitary waves, this method converges if and only if the solitary wave is linearly stable. The second method is similar to the first method except that it uses a novel amplitude normalization. The performance of these methods is illustrated on various examples. It is found that while the first method is competitive with the Petviashvili method, the second method delivers much better performance than the first method and the Petviashvili method.     

A Random Numerical Method with Application in Combustion

This paper presents a random numerical method which combines the random vortex method and the random choice method. A random choice method is used for the modeling reaction-diffusion system. The splitting of source terms in two dimensions for the random choice method is tested. A hybrid random vortex method is used for solving the Navier-Stokes equation which governs the fluid motion. With the assumption of incompressibility, the fluid motion can be uncoupled from the chemistry. The method is applied to a flow passing a circular cylinder which is kept cold or heated. In both cases the method demonstrates an ability to resolve turbulent effects on flame front propagation.     

拟Raabe判别法与拟对数判别法的强弱关系

拟Raabe判别法是新近提出的关于正项级数收敛性的一种比较细致的判别法.对通项递减的正项级数来说,此判别法强于传统的Raabe判别法与Gauss判别法.通过对拟Raabe判别法与另一个细致的判别法——拟对数判别法强弱关系的探讨,得出了后一判别法强于前者的结论.     

A COLUMN-SECANT UPDATE TECHNIQUE FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

In this paper,we present a column-secant modification of the SCC method,which is called the CSSCC method.The CSSCC method uses function values more efficiently than the SCC method,and it is shown that the CSSCC method has better local q-convergence and r-convergence rates than the SCC method.The numerical results show that the CSSCC method is competitive with some well known methods for some standard test problems.