关于几何规划迭代法的讨论

本文利用正定几何规划原问题与对偶问题的关系以及迭代结果,提出了具体估计理论解存在范围的方法;针对缩并函数性态不良时,变量在迭代过程中易来回振荡,不易收敛,甚至不收敛的问题,提出了在迭代中对对偶变量增加一修正步,以减小变量振荡的方法;并由此使从可行域外,用一系列正定几何规划逐次逼近原问题最优解的方法获得成功,从而扩大了Avriel-Williams方法的应用范围,改善了缩并函数性态不良时的收敛性.     

摄动法结合Pad逼近在结构拓扑重分析中的应用

讨论了矩阵摄动理论结合Pad啨逼近在结构拓扑修改重分析中的应用,利用分步迭代的方法来取得高精度的近似解。定义过渡方程并利用原始结构信息得到其精确解;利用正交化的摄动基作Pad啨逼近,并采用迭代的方法得出对过渡方程解的增量,从而得到修改结构的近似解。     

摄动法结合Padé逼近在结构拓扑重分析中的应用

讨论了矩阵摄动理论结合Padé逼近在结构拓扑修改重分析中的应用,利用分步迭代的方法来取得高精度的近似解.定义过渡方程并利用原始结构信息得到其精确解;利用正交化的摄动基作Padé逼近,并采用迭代的方法得出对过渡方程解的增量,从而得到修改结构的近似解.     

弹性接触问题线性互补法的摄动原理

为使弹性摩擦问题的线性互补法得到收敛性保证,提出一个摄动原理,证明了通过“摆脱摄动误差”可得到原离散模型的精确解。     

弹性地基上圆板在轴对称载荷作用下的迭代解

提出了弹性地基上圆板承受轴对称载荷作用弯曲问题的一种新的解法－迭代法。文中导出了迭代过程的一般公式，并给出了关于收敛性的一般说明。未方法简单可靠，能收敛于精确解，可直接计算弹性地基圆板在轴对称载荷作用下的弯曲问题。文后给出算例，只二次迭代即得到了满意的结果。     

非线性结构动力学瞬态分析的摄动法

在用直接积分法求解非线性结构的动力响应时,常常需要做迭代运算。本文引入摄动方法后,加快了收敛速度,提高了计算效益。     

多自由度系统复模态理论的摄动方法 (三)非对称系统复特征值问题的摄动迭代法

基于实模态理论和相应特征值计算方法,本文在位形空间中用摄动方法求解非对称系统的复特征值和复模态。为提高精度,增加了一个迭代格式。最后给出算例并讨论了计算精度和收敛性。     

矩形板耦合热冲击问题的摄动解

本文通过对薄板耦合热弯曲问题的完备方程的无量纲化,引出了关于薄板的无量纲热弹性耦合系数,并以此系数为摄动参数,运用奇异摄动方法,导出了其摄动方程,得到了关于矩形薄板耦合热冲击问题的一致有效的渐近解。然后,在对该解讨论及计算的基础上,获得了一些关于该类问题的规律性的结论。     

椭圆抛物面扁壳某些应力集中问题

本文采用摄动法讨论了椭圆抛物面扁壳在集中载荷附近和圆孔附近的应力集中问题,所取的参数为ε=(k_x-k_y)/(k_x+k_y)(这里k_x和k_y是壳面的主曲率).按摄动法求解时,其初始解就是曲率为(k_x+k_y)/2的球壳的解,而各次修正解也归结为球壳方程问题的解.文中采用Hankel变换方法求得了各次修正解,并给出了集中力问题最初几次修正解的数字结果.数字计算表明,摄动法的收敛性相当良好.     

密集模态摄动的新方法

本文提出了一种密集模态结构系统(M_0,K_0)振动分析的矩阵摄动新方法.它将密集模态结构系统特征解的摄动问题转化为重特征值的摄动问题.文中给出了一个数值例子.     

双相介质波动方程孔隙率反演的同伦方法

从材料响应的理论合成应与实际测量数据相拟合这一出发点，将双相介质波劝方程参数的反演问题转化为非线性算子方程的零点求解问题，从而应用一种大范围收敛的同伦方尘土注来解非线性算子方程，并把这种方法用于Simon(1984)给出的具有解析的一维双相介质模型的数值模拟，最后的数值结果表明，给出的算法是十分有效的。     

On time independent Schr?dinger equations in quantum mechanics by the homotopy analysis method

A general analytic approach, namely the homotopy analysis method(HAM), is applied to solve the time independent Schr?dinger equations. Unlike perturbation method, the HAM-based approach does not depend on any small physical parameters, corresponding to small disturbances.Especially, it provides a convenient way to gain the convergent series solution of quantum mechanics. This study illustrates the advantages of this HAM-based approach over the traditional perturbative approach, and its general validity for the Schr?dinger equations. Note that perturbation methods are widely used in quantum mechanics, but perturbation results are hardly convergent. This study suggests that the HAM might provide us a new, powerful alternative to gain convergent series solution for some complicated problems in quantum mechanics, including many-body problems, which can be directly compared with the experiment data to improve the accuracy of the experimental findings and/or physical theories.     

处理一种非线性问题迭代法的计算机方法

本文依非线性问题的积分方程组,按迭代法得出一般迭代表达式.由此,任意次迭代解的所有待定系数可用计算机计算,使研究复杂的非路性问题成为可能.     

Comparison between the HAM and HPM solutions of thin film flows of non-Newtonian fluids on a moving belt

In this paper, we prove in general that the homotopy perturbation method (HPM) proposed in 1998 is only a special case of the homotopy analysis method (HAM) profound in 1992 when ħ = −1. Besides, by using the thin film flows of Sisko and Oldroyd 6-constant fluids on a moving belt as examples, we show that the solutions given by HPM (Siddiqui, A.M., Ahmed, M., Ghori, Q.K.: Chaos Solitons and Fractals (2006) in press) are divergent, and thus useless. However, by choosing a proper value of the auxiliary parameter ħ, we give convergent series solution by means of the HAM. These two examples also show that, different from the HPM and other traditional analytic techniques, the HAM indeed provides us with a simple way to ensure the convergence of the solution.     

Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media

Abstract. The Picard and modified Picard iteration schemes are often used to numerically solve the nonlinear Richards equation governing water flow in variably saturated porous media. While these methods are easy to implement, they are only linearly convergent. Another approach to solve the Richards equation is to use Newton's iterative method. This method, also known as Newton–Raphson iteration, is quadratically convergent and requires the computation of first derivatives. We implemented Newton's scheme into the mixed form of the Richards equation. As compared to the modified Picard scheme, Newton's scheme requires two additional matrices when the mixed form of the Richards equation is used and requires three additional matrices, when the pressure head-based form is used. The modified Picard scheme may actually be viewed as a simplified Newton scheme.Two examples are used to investigate the numerical performance of different forms of the 1D vertical Richards equation and the different iterative solution schemes. In the first example, we simulate infiltration in a homogeneous dry porous medium by solving both, the h based and mixed forms of Richards equation using the modified Picard and Newton schemes. Results shows that, very small time steps are required to obtain an accurate mass balance. These small times steps make the Newton method less attractive.In a second test problem, we simulate variable inflows and outflows in a heterogeneous dry porous medium by solving the mixed form of the Richards equation, using the modified Picard and Newton schemes. Analytical computation of the Jacobian required less CPU time than its computation by perturbation. A combination of the modified Picard and Newton scheme was found to be more efficient than the modified Picard or Newton scheme.     

Numerical solution of a singularly perturbed elliptichyperbolic partial differential equation on a nonuniform discretization mesh

In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order convergent, uniformly in the small parameter ε, to the solution of problem (1.1). Numerical results are finally provided.     

Using reproducing kernel for solving a class of partial differential equation with variable-coefficients

How to solve the partial differential equation has been attached importance to by all kinds of fields. The exact solution to a class of partial differential equation with variable-coefficient is obtained in reproducing kernel space. For getting the approximate solution, give an iterative method, convergence of the iterative method is proved. The numerical example shows that our method is effective and good practicability.     

概率约束评估的功能度量法的混沌控制

功能度量法是基于可靠度的结构优化设计中评估概率约束的一种方法,其改进均值(AMV)迭代格式具有简洁、高效的优点,但对一些非线性功能函数搜索最小功能目标点时可能陷入周期振荡或混沌解,本文利用混沌反馈控制的稳定转换法对功能度量法的AMV迭代格式实施收敛控制.首先展示一些功能函数应用功能度量法AMV格式迭代计算产生了周期解和混沌解现象,并对迭代算法进行了混沌动力学分析.然后利用稳定转换法对功能度量法迭代失败的参数区间进行混沌控制,使嵌入周期和混沌轨道的不稳定不动点稳定化,获得了稳定收敛解,实现了迭代解的周期振荡、分岔和混沌控制.     

A simple and efficient BEM algorithm for planar cavity flows

This paper presents a boundary element formulation for solution of planar Riabouchinsky cavity flow problems. An iterative procedure for adjusting the free surface position is developed and shown to be stable and convergent. Numerical results are compared with finite difference and finite element solutions, showing the superior accuracy of the BEM models.     

Using reproducing kernel for solving a class of partial differential equation with variable-coefficients

How to solve the partial differential equation has been attached importance to by all kinds of fields. The exact solution to a class of partial differential equation with variable-coefficient is obtained in reproducing kernel space. For getting the approximate solution, give an iterative method, convergence of the iterative method is proved. The numerical example shows that our method is effective and good practicability.