深埋隧洞围岩应力的精确解与近似解的对比分析

对不同断面形状的深埋隧洞进行了分析,比较了隧洞围岩应力解析解与通过当量半径方法得到的近似解之间的差别.首先,应用复变函数的基本理论,给出圆形、椭圆、矩形、直墙拱形等几种常见深埋隧洞围岩应力的解析表达式.其次,应用当量半径的折算形式,将其任意形状的边界转化为标准圆形断面,利用Lamé解答得到了各围岩应力分量.最后,考虑隧洞断面形状参数的变化,通过数值算例对精确解和近似解进行了比较,分析了当量半径折算形式的精确度.在此基础上,应用有限元方法验证了复变函数解析解的精确性,以椭圆、矩形和直墙拱形的复变函数解验证当量半径精确度.结果表明,当量半径的折算形式解答与精确解答之间相似程度与隧洞的断面形状和几何参数之间有着密切的关系.     

任意变系数微分方程的精确解析法

工程中的许多问题归结为求解任意变系数微分方程的解.本文首次提出精确解析法,用以求解任意变系数微分方程在任意边界条件下的解.文中还给出精确解析法的一般计算格式,得到了一致收敛于精确解及其任意阶导数的解析表达式,并给出收敛性证明.文末给出四个算例,均得到较好的结果,证明了本文理论的正确性.     

随机Navier-Stokes方程数值解的收敛性

Navier-Stokes方程在流体力学中有广泛的应用.通常情况下,大多数Navier-Stokes方程没有精确解,数值方法显得尤为重要.本文根据BDM法,利用It公式,Burkholder-Davis-Gundy不等式,Doob不等式和Gronwall引理对随机Navier-Stokes方程数值解的收敛性进行了讨论,得出数值解均方意义下收敛到解析解.     

同伦分析法求解Burgers方程的初边值问题

采用同伦分析法求解了Burgers方程的一初边值问题,得到了它的近似解析解.在不同粘性系数情形下,对近似解与精确解进行了比较,发现在粘性系数不是非常小的情况下,用此方法得到的解析解与精确解符合地很好.     

Winkler地基上四边自由正交各向异性矩形中厚板弯曲的有限积分变换解（英文）

本文应用有限积分变换法研究Winkler地基上四边自由正交各向异性矩形中厚板的弯曲问题.具体由正交各向异性矩形中厚板弯曲的基本方程组和边界条件出发,结合有限积分变换法及其对应的逆变换法推导出正交各向异性矩形中厚板弯曲问题的解析解.该解析解统一适用于计算各向同性/正交各向异性矩形薄板、中厚板和厚板的弯曲问题,并且通过具体算例验证了所得解析解的正确性.     

简支夹层矩形板的非线性弯曲

本文应用变分法导出了具有软夹心的夹层矩形板的非线性弯曲理论的基本方程和边界条件.然后,使用摄动法研究了均布横向载荷作作用下简支夹层矩形板的非线性弯曲问题,得到了相当精确的解析解.     

精确的磁流体动力学汇流解析解

就一个特殊的磁流体动力学(MHD)流动,即速度幂指数为-1时的汇流,得到著名的Falkner-Skan方程精确的解析解.解析解是封闭的,并有多重解分支.分析了磁场参数和壁面伸长参数的影响.发现了有趣的速度分布现象:即使壁面固定,回流区域依然出现.在一个罕见的FalknerSkan MHD流动中,得到了一组解,以精确封闭的解析公式表示,极大地丰富了著名的Falkner-Skan方程的解析解,也加深了对这重要又有趣方程的理解.     

复合荷载下波纹圆板的非线性分析*

本文按照各向同性和正交各向异性圆板的大挠度理论,研究了具有光滑中心的波纹圆板在均布和中心集中荷载联合作用下的非线性弯曲问题.应用修正迭代法,我们得到了夹紧固定和滑动固定两种边界条件下十分精确的解析解.     

单位圆内解析系数的高阶齐次线性微分方程解与小函数的关系

本文对系数为单位圆内的解析函数的某类高阶齐次线性微分方程解的性质进行研究,得到解的增长级和超级的精确估计以及解与小函数之间的关系.     

非均匀双向加肋圆柱壳非线性轴对称变形的精确解析法

本文首次利用精确解析法分析了环向和纵向加肋非均匀圆柱壳在任意载荷和边界条件下非线性轴对称变形问题.导出了一致收敛于精确解的位移和内力解析表达式,文中给出收敛性问题.问题最后归结为求解二元一次代数方程组,计算既简便又迅速.文末给出四个数值算例表明,本文提出的方法,可以得到满意的结果.     

On a direct method for the solution of nearly uncoupled Markov chains

Summary This note is concerned with the accuracy of the solution of nearly uncoupled Markov chains by a direct method based on the LU decomposition. It is shown that plain Gaussian elimination may fail in the presence of rounding errors. A modification of Gaussian elimination with diagonal pivoting and correction of small pivots is proposed and analyzed. It is shown that the accuracy of the solution is affected by two condition numbers associated with aggregation and the coupling respectively.This work was supported in part by the Air Force Office of Sponsored Research under Contract AFOSR-87-0188     

A new application of the homotopy analysis method in solving the fractional Volterra’s population system

This paper considers a Volterra’s population system of fractional order and describes a bi-parametric homotopy analysis method for solving this system. The homotopy method offers a possibility to increase the convergence region of the series solution. Two examples are presented to illustrate the convergence and accuracy of the method to the solution. Further, we define the averaged residual error to show that the obtained results have reasonable accuracy.     

Higher-order compact finite difference method for systems of reaction–diffusion equations

This paper is concerned with a compact finite difference method for solving systems of two-dimensional reaction–diffusion equations. This method has the accuracy of fourth-order in both space and time. The existence and uniqueness of the finite difference solution are investigated by the method of upper and lower solutions, without any monotone requirement on the nonlinear term. Three monotone iterative algorithms are provided for solving the resulting discrete system efficiently, and the sequences of iterations converge monotonically to a unique solution of the system. A theoretical comparison result for the various monotone sequences is given. The convergence of the finite difference solution to the continuous solution is proved, and Richardson extrapolation is used to achieve fourth-order accuracy in time. An application is given to an enzyme–substrate reaction–diffusion problem, and some numerical results are presented to demonstrate the high efficiency and advantages of this new approach.     

Analytical solution to nonlinear oscillation system of the motion of a rigid rod rocking back using max–min approach

This paper used an ideal periodic solution which is called max–min approach (MMA) to evaluate oscillation systems with nonlinearity terms such as motion of a rigid rod rocking back. This method introduces an alternative to overcome the difficulty of computing the periodic behavior of the oscillation problems in engineering. To assess the accuracy of solutions, the results were compared with the exact ones. The most significant features of this method are the simplicity and the excellent agreement with the exact results for the various parameters. Furthermore, the results reveal that one iteration leads to high accuracy of the solution. This solution may be useful for the explanation of some practical physical problems.     

搭接焊缝内应力分布规律的准确解

本文给出了搭接焊缝内应力分布规律的准确解，从而使评价以往的有关研究成果、改进迄今采用的设计方法和设计规范成为可能。     

The modified regularization method for identifying the unknown source on Poisson equation

This paper discusses the problem of determining an unknown source which depends only on one variable in two-dimensional Poisson equation from one supplementary temperature measurement at an internal point. The problem is ill-posed in the sense that the solution (if it exists) does not depend continuously on the data. The regularization solution is obtained by the modified regularization method. For the regularization solution, the Hölder type stability estimate between the regularization solution and the exact solution is given. Numerical results are presented to illustrate the accuracy and efficiency of this method.     

A dynamic meshless method for the least squares problem with some noisy subdomains

The interpolation method by radial basis functions is used widely for solving scattered data approximation. However, sometimes it makes more sense to approximate the solution by least squares fit. This is especially true when the data are contaminated with noise. A meshfree method namely, meshless dynamic weighted least squares (MDWLS) method, is presented in this paper to solve least squares problems with noise. The MDWLS method by Gaussian radial basis function is proposed to fit scattered data with some noisy areas in the problem’s domain. Existence and uniqueness of a solution is proved. This method has one parameter which can adjusts the accuracy according to the size of noises. Another advantage of the developed method is that it can be applied to problems with nonregular geometrical domains. The new approach is applied for some problems in two dimensions and the obtained results confirm the accuracy and efficiency of the proposed method. The numerical experiments illustrate that our MDWLS method has better performance than the traditional least squares method in case of noisy data.     

On the effect of numerical integration in the Galerkin boundary element method

Summary. We discuss the effect of cubature errors when using the Galerkin method for approximating the solution of Fredholm integral equations in three dimensions. The accuracy of the cubature method has to be chosen such that the error resulting from this further discretization does not increase the asymptotic discretization error. We will show that the asymptotic accuracy is not influenced provided that polynomials of a certain degree are integrated exactly by the cubature method. This is done by applying the Bramble-Hilbert Lemma to the boundary element method. Received May 24, 1995     

Analysis of a fourth-order compact ADI method for a linear hyperbolic equation with three spatial variables

This paper is concerned with a three-level alternating direction implicit (ADI) method for the numerical solution of a 3D hyperbolic equation. Stability criterion of this ADI method is given by using von Neumann method. Meanwhile, it is shown by a discrete energy method that it can achieve fourth-order accuracy in both time and space with respect to H ¹- and L ²-norms only if stable condition is satisfied. It only needs solution of a tri-diagonal system at each time step, which can be solved by multiple applications of one-dimensional tri-diagonal algorithm. Numerical experiments confirming the high accuracy and efficiency of the new algorithm are provided.     

The truncation method for identifying an unknown source in the Poisson equation

This paper deals with the problem of determining an unknown source which depends only on one variable in two-dimensional Poisson equation, with the aid of an extra measurement at an internal point. The problem is ill-posed in the sense that the solution (if it exists) does not depend continuously on the data. We obtain the regularization solution by the truncation method. For the regularization solution, the Hölder type stability estimate between the regularization solution and the exact solution is given. Numerical results are presented to illustrate the accuracy and efficiency of this method.