Singular finite elements for the sudden-expansion and the die-swell problems

The singular finite element method is used to solve the sudden-expansion and the die-swell problems in order to improve the accuracy of the solution in the vicinity of the singularity and to speed up the convergence. The method requires minor modifications to standard finite element schemes, and even coarse meshes give more accurate results than refined ordinary finite element meshes. Improved normal stress results for the sudden-expansion problem have been obtained for various Reynolds numbers up to 100 using the singular elements constructed for the creeping flow problem. In addition, the normal stresses at the walls appear to be insensitive to the singularity powers used in the construction of the singular basis functions. The die-swell problem is solved using the singular elements constructed for the stick–slip problem. The singular elements accelerate the convergence of the free surface dramatically.     

Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.     

A uniform high-order method for a singular perturbation problem in conservative form

A uniform high order method is presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems (1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O(h~m+1)accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O(h~m+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.     

Experimental estimation of the load distribution in bearings by the method of caustics

The objective of this paper is the experimental estimation of the load distribution in roller-bearing by the metrhod of caustics. Contact problems have many practical and important applications¹. For the solution of such problems, besides mechanical analysis, the experimental method of caustics² can also be applied. The optical method of caustics is suitable for the experimental study of singularities in stress fields created either by discontinuities or by loading. Previously, caustics has been applied to the study of singular stress fields developed near concentrated or uniformly distributed loads which are applied along straight boundaries.^3,4 In this work, it is applied to study the load distribution in rollerbearings.     

边界积分方程中近奇异积分计算的一种变量替换法

准确估计近奇异边界积分是边界元分析中一项很重要的课题,其重要性仅次于对奇异积分的处理. 近年来已发展了许多方法,都取得了一定程度的成功,但这个问题至今仍未得到彻 底的解决. 基于一种新的变量变换的思想和观点,提交了一种通用的积分变换法, 它非常有效地改善了被积函数的震荡特性,从而消除了积分的近奇异性,在不增加计算量的情况 下, 极大地改进了近奇异积分计算的精度. 数值算例表明,其算法稳定,效率高, 并可达到很高的计算精度,即使区域内点非常地靠近边界,仍可取得很理想的结果.     

Interaction of a die with a layered elastic foundation

Plane and axisymmetric contact problems for a three-layer elastic half-space are considered. The plane problem is reduced to a singular integral equation of the first kind whose approximate solution is obtained by a modified Multhopp-Kalandiya method of collocation. The axisymmetric problem is reduced to an integral Fredholm equation of the second kind whose approximate solution is obtained by a specially developed method of collocation over the nodes of the Legendre polynomial. An axisymmetric contact problem for an transversely isotropic layer completely adherent to an elastic isotropic half-space is also considered. Examples of calculating the characteristic integral quantities are given. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 165–175, May–June, 2006.     

双曲偏微分方程的局部伪弧长方法研究

重点研究了局部伪弧长方法在处理偏微分方程,尤其是双曲型偏微分方程出现激波间断的奇异性问题,对比分析了全局伪弧长方法空间转化的形式及其网格自适应的性质。为提高求解效率,提出了局部伪弧长方法,利用激波间断的性质,给出了判断奇异点位置以及模板选择的方法,涉及如何处理激波振荡,如何引入弧长参数,以及怎样求解间断等问题。通过数值算例验证了局部伪弧长在激波捕捉和追踪方面的可行性,通过比较局部伪弧长方法与Godunov方法处理不同初值条件的双曲问题,显示出局部伪弧长方法处理双曲偏微分方程的优越性,为伪弧长方法应用到物理问题奠定基础。     

三维位势问题的梯度边界积分方程的新解法

三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题.已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算.本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法.该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解.值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差.辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.     

A new approach for the solution of singular optimum in structural topology optimization

In order to overcome the difficulties caused by singular optima, in the present paper, a new method for the solutions of structural topology optimization problems is proposed. The distinctive feature of this method is that instead of solving the original optimization problem directly, we turn to seeking the solutions of a sequence of approximated problems which are formulated by relaxing the constraints of the original problem to some extent. The approximated problem can be solved efficiently by employing the algorithms developed for sizing optimization problems because its solution is not singular. It can also be proved that when the relaxation parameter is tending to zero, the solution of the approximated problem will converge to the solution of the original problem uniformly. Numerical examples illustrate the effectiveness and validity of the present approach. Results are also compared with those obtained by traditional methods. The project supported by the National Natural Science Foundation of China under project No. 19572023     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

A numerical method for the evaluation of hydrodynamic forces of translating bodies under a free surface

This paper presents a numerical method to evaluate the hydrodynamic forces of translating bodies under a free surface. Both steady and unsteady problems are considered. Analytical and numerical studies are carried out based on the Havelock wave‐source function and the integral equation method. Two main problems arising inherently in the proposed solution method are overcome in order to facilitate the numerical implementation. The first lies in evaluating the Havelock function, which involves integrals with highly oscillatory kernels. Particular integration contours leading to non‐oscillatory integrands are derived a priori so that the integrals can be evaluated efficiently. The second problem lies in evaluating singular kernels in the boundary integral equation. The corresponding non‐singular formulation is derived using some theorems of potential theory, including the Gauss flux theorem and the property related to the equipotential body. The subsequent formulation is amenable to the solution by directly using the standard quadrature formulas without taking another special treatment. This paper also attempts to enhance the computational efficiency by presenting an interpolation method used to evaluate matrix elements, which are ascribed to a discretization procedure. In addition to the steady case, numerical examples consist of cases involving a submerged prolate spheroid, which is originally idle and then suddenly moves with a constant speed and a constant acceleration. Also systematically studied is the variation of hydrodynamic forces acting on the spheroid for various Froude numbers and submergence depths. Copyright © 2000 John Wiley & Sons, Ltd.     

A singular perturbation problem for periodic boundary partial differential equation

In this paper,we consider a singular perturbation elliptic-parabolic partial differentialequation for periodic boundary value problem,and construct a difference scheme.Using themethod of decomposing the singular term from its solution and combining an asymptoticexpansion of the equation,we prove that the scheme constructed by this paper convergesuniformly to the solution of its original problem with O(τ h~2).     

A singular finite element for Stokes flow: The stick–slip problem

Abrupt changes in boundary conditions in viscous flow problems give rise to stress singularities. Ordinary finite element methods account effectively for the global solution but perform poorly near the singularity. In this paper we develop singular finite elements, similar in principle to the crack tip elements used in fracture mechanics, to improve the solution accuracy in the vicinity of the singular point and to speed up the rate of convergence. These special elements surround the singular point, and the corresponding field shape functions embody the form of the singularity. Because the pressure is singular, there is no pressure node at the singular point. The method performs well when applied to the stick–slip problem and gives more accurate results than those from refined ordinary finite element meshes.     

Asymptotic solution of nonlocal nonlinear reaction-diffusion Robin problems with two parameters

In this paper, the nonlocal nonlinear reaction-diffusion singularly perturbed problems with two parameters are studied. Using a singular perturbation method, the structure of the solutions to the problem is discussed in relation to two small parameters. The asymptotic solutions of the problem are given.     

On interface crack in laminated anisotropic medium

In this paper, the interface crack problems of a multilayered anisotropic medium under a state of generalized plane deformation are considered within the framework of anisotropic theory. A general solution procedure is introduced such that it can be uniformly applied to media with transversely isotropic, orthotropic, monoclinic, etc. layers. The problem is reduced to the solution of a system of singular integral equations by means of Fourier transform method and the stiffness matrix formulation. A Jacobi polynomial technique is then used to solve the integral equations numerically. The stress intensity factors are provided. The stress intensity factors have been calculated numerically and displayed graphically.     

一种求解弹性接触问题的缩减二次规划方法

基于势能原理以节点位移为设计变量、以接触条件为约束方程构建了无摩擦弹性接触问题的二次规划数学模型,在此基础之上利用力平衡线性约束方程的特解和由基础解向量构成的奇异模态矩阵,提出一种新的基于奇异坐标变换的自由度缩减方法,大大降低了二次规划的规模,并使得二次规划模型不再含显性等式约束;根据弹性接触力学体系的特点,通过人为假定接触自由度位移模式,提出了一种简单高效的奇异模态矩阵的计算方法。通过两圆柱接触、轴孔间隙配合接触两个数值算例的对比分析,验证了对于弹性接触问题的求解,缩减二次规划方法有效克服了传统方法计算量大、对求解参数设置敏感、收敛困难的问题。     

Precise integration method for solving singular perturbation problems

This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.     

A REDUCED QUADRATIC PROGRAMMING METHOD FOR ELASTIC CONTACT PROBLEMS

基于势能原理以节点位移为设计变量、以接触条件为约束方程构建了无摩擦弹性接触问题的二次规划数学模型,在此基础之上利用力平衡线性约束方程的特解和由基础解向量构成的奇异模态矩阵,提出一种新的基于奇异坐标变换的自由度缩减方法,大大降低了二次规划的规模,并使得二次规划模型不再含显性等式约束;根据弹性接触力学体系的特点,通过人为假定接触自由度位移模式,提出了一种简单高效的奇异模态矩阵的计算方法。通过两圆柱接触、轴孔间隙配合接触两个数值算例的对比分析,验证了对于弹性接触问题的求解,缩减二次规划方法有效克服了传统方法计算量大、对求解参数设置敏感、收敛困难的问题。     

Integral transform solution of multi-material structure with finite crack perpendicular to the interface

The plane elasticity problem for layered elastic systems containing a finite crack perpendicular to the interface is considered. To derive the singular integral equations. Fourier transform in conjunction with dislocation is used. The singular integral equation is solved with the Lobatto-Chebyshev method commonly applied to such problems. In order to have an idea about the usefulness of the method described, a two-layer structure which contains a cut parallel toh is considered.     

域外奇点法解杆的弹塑性扭转问题

本文提出一种借助于沙丘比拟的求解杆弹塑性扭转问题的域外奇点法,这种方法可降低所求问题的维数,有效地避免解的奇异性。它具有方法简单,不需要数值积分,计算时间短和精度高等优点。