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.     

The integrated singular basis function method for the stick-slip and the die-swell problems

We further develop a new singular finite element method, the integrated singular basis function method (ISBFM), for the solution of Newtonian flow problems with stress singularities. The ISBFM is based on the direct subtraction of the leading local solution terms from the governing equations and boundary conditions of the original problem, followed by a double integration by parts applied to those integrals with singular contributions. The method is applied to the stick-slip and the die-swell problems and improves the accuracy of the numerical results in both cases. In the case of the die-swell problem it considerably accelerates the convergence of the free surface profile with mesh refinement. The advantages and disadvantages of the ISBFM when compared to other singular methods are also discussed.     

Converged solutions of the Newtonian extrudate‐swell problem

Both the axisymmetric and the planar Newtonian extrudate‐swell problems are solved using the standard and the singular finite element methods. In the latter method, special elements that incorporate the radial form of the stress singularity are used around the exit of the die. The convergence of each of the two methods with mesh refinement is studied for various values of the Reynolds and the capillary numbers. The numerical results show that the singular finite elements perform well if coarse or moderately refined meshes are used, and appear to be superior to the standard finite elements only when the Reynolds number is low and the surface tension is not large. The standard finite elements perform better as the surface tension or the Reynolds number are increased. This implies that the effect of the stress singularity on the accuracy of the numerical solution in the neighborhood of the die exit becomes less significant when the Reynolds number is high or the surface tension is large. Copyright © 1999 John Wiley & Sons, Ltd.     

用拟协调元直接构造平面任意四边形单元——进入有限元的禁区

利用拟协调元方法,在直角坐标系下直接构造了一族平面任意四边形单元,对其收敛性进行了分析,并与平面等参元进行了对比研究.结果证明平面任意四边形单元可采用多项式基函数直接列式,并可以保障单元的收敛性;拟协调元列式可以使平面问题的有限元方法得到统一.与平面等参元相比,单元列式简单,性能稳定,具有显式的刚度阵,计算量小,这说明对于有限元平面问题拟协调元是一个更正确、有效的做法.      

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

Diffraction of waves by semi-infinite breakwater using finite and infinite elements

A finite and infinite element model is derived to predict wave patterns around a semi-infinite breakwater in water of constant depth. Both circular and square meshes of elements are used. The wave theory used is that of Berkhoff. The appropriate boundary conditions for finite and infinite boundaries are described. The singularity in the velocity at the breakwater tip is modelled effectively using the technique of Henshell and Shaw originally developed in elasticity. The results agree well with the analytical solution. In addition the problem of waves incident upon a semi-infinite breakwater and parabolic shoal, where both diffraction and refraction are present, is solved. There is no analytical solution for this case. The combination of finite and infinite elements is found to be an effective and accurate technique for such problems.     

弱连续条件下的九参三角形板元

首先对薄板弯曲平衡方程的弱形式进行了推导,导出保证单元收敛的弱协调条件,即三角形顶点函数值连续和三边的法向导数积分连续这两个条件;对比拟协调元、广义协调元和双参数法中所使用的3个积分连续条件,本条件更弱;再对这3个积分协调条件的构成方法进行了总结和分析,现有采用积分连续条件构造的有限元大都采用了这些构成方法.采用弱协调条件构造有限元,比原来的构造范围更广,井以此构造出几种单元作为算例.采用这种构成法还可构造多种单元,它们都具有采用最小势能原理法构成有限元的简便的优点,并在任意网格下收敛到真解.     

弹塑性力学问题的虚单元法

具有有限差分法特征的虚单元法,可视为是有限元法向任意多边形单元的扩展。在材料细观力学性能表征、非均质材料力学分析等非线性问题方面,传统的弹塑性有限元法具有网格数目多、效率低下等不足之处,而虚单元法使网格划分更加灵活,为材料的弹塑性力学分析等非线性问题提供了新的思路。基于增量法弹塑性力学原理和双线性投影算子,建立了弹塑性力学问题的虚单元法求解技术,提出了弹塑性力学问题虚单元法的应力更新方案,研究了弹性力学问题虚单元法的精度和收敛性,讨论了虚单元法求解弹塑性力学问题的网格依赖性。同时,开展了任意多边形和凹多边形单元的数值试验研究,结果表明,虚单元法无须分割多边形,仅需节点自由度便可求得单元刚度矩阵和应力等效荷载,程序实现简单,计算精度高,改善了传统有限元的网格依赖性和塑性区的网格奇异性。     

Analytical removal of singularity and direct evaluation of singular integrals in 3-D elastoplastic finite deformation analysis by BEM

As a further development of the present authors' research work [1,2], in this paper a method of the so-called quadratic pentahedron polar co-ordinate transformation and analytical removal of singularity of Cauchy principal value singular integrals is proposed to evaluate the strongly singular integrals in the sense of Cauchy principal values and the weakly singular integrals over quadratic internal cells in 3-D elastoplastic finite deformation analysis by BEM. First, a quadratic pentahedron polar co-ordinate transformation technique is used to reduce the order of singularity of the singular integrals. Then, a form of Gauss' theorem is introduced to remove the singularity in the Cauchy principal value singular integrals analytically. Therefore, the evaluation of all those strongly and weakly singular integrals can be carried out by standard Gaussian quadrature accurately and efficiently. Numerical examples of the 3-D elastoplastic problem and 3-D finite deformation problem are given to demonstrate that the method possesses good accuracy and numerical stability, and is convenient to implement. The method in this paper can be applied extensively to evaluating the singular integrals over cubic and higher order elements.     

一种确定结合材料界面端应力强度因子的数值方法

基于弹性力平面问题的基本方程，给出了结合材料界面端的应力奇异性特征方程以及位移场和奇异应力场。提出了一种确定结合材料界面端应力强度因子的数值外插方法。对界面端区域进行了有限元网格单元划分。经过具体实例检验进一步确定了求解应力强度因子的最佳方向，该数值外插法的计算结果精度符合工程应用的要求，为工程材料强度的评价提供了有效的计算途径。     

无奇异性边界积分方程法

边界元法中如何有效地处理奇异积分,一直是人们极为关心的课题。本文提出了建立由“线源”产生的边界积分方程,其优点是可任意阶地降低奇异性阶次。计算实例表明,本文所提出的方法优于常规的边界元法。尤其是提高了位移的计算精度。对边界单元网格不等长划分,更显示出该法的优越性。     

Three-dimensional fracture analysis using tetrahedral enriched elements and fully unstructured mesh

In the absence of automated and customized methods and tools, some of today’s existing methods for solving three-dimensional fracture problems require comprehensive finite element meshing, labor-intensive analysis and post-processing efforts. In this study, a tetrahedral enriched element method and related applications are presented that demonstrate employment of fully unstructured tetrahedral meshes for general mixed-mode three-dimensional fracture problems. As in the case of hexahedral enriched elements, the tetrahedral enriched elements also alleviate the needs of pre- and post-processing the finite element model, allowing direct computation of stress intensity factors in the solution phase. In addition, when tetrahedral enriched elements are used, the crack front region can also be meshed using unstructured elements allowing direct use of automatic free-meshing programs. The applications presented are plane-strain central crack problem, mode-I surface crack in a plate, inclined penny-shaped crack, edge-cracked bar under constant heat flux and lens-shaped crack embedded in a large elastic body. The results obtained are in good comparative agreement with those available in the literature. Thus, it is concluded that the enriched tetrahedral elements can be applied efficiently and accurately on a general three-dimensional fracture problem allowing usage of fully unstructured finite element meshes.     

A finite element method for liquid with memory

Steady, two-dimensional flows of viscoelastic liquid with memory are analyzed by means of a nonlinear integral constitutive equation, novel streamlined finite elements, Galerkin's method of weighted residuals, and Newton iteration. By their relation to streamlines the new elements allow, for the first time, full Newton iteration of the algebraic equation set to which the governing integrodifferential system reduces. Streamlines are computed simultaneously with velocity components and pressure, the primary unknowns; to track the history of liquid particles the deformation equations are solved analytically in a Protean system of coordinates that conforms to the streamlines. In the cases studied the Newton iteration converges quadratically to a creeping flow state of Weissenberg number (product of upstream wall shear rate and relaxation time of the liquid) up to about 20 in channel flow, 10 in film flow, and 2 in die-swell flow. Shear-thinning shrinks the domain of convergence; it also reduces die-swell. A slip boundary condition in the vicinity of the contact line extends the domain of convergence.     

三维弹塑性有限变形问题边界元法中奇异积分的直接分析去奇解法

作为本文作者研究工作的继续,本文提出了处理三维弹塑性有限变形问题边界元法中二次元区域弱奇及Cauchy 主值奇异积分的二次极坐标变换—分析去奇法.该方法先通过适当的二次极坐标变换降低奇异积分的奇异性,然后利用Causs 散度定理去除Cauchy主值积分的奇异性.通过三维弹塑性及三维有限变形问题数值算例说明该方法具有良好的精度及数值稳定性,并且实施较方便.本文方法可直接推广应用于二阶以上高阶元离散模型奇异积分处理.     

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

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

拟协调元研究综述

拟协调元是有限元中十分重要的、具有特色的一种列式体系. 拟协调元列式简单、灵活, 统一了协 调元、非协调元等列式方法. 在列式中, 拟协调元将几何方程和平衡方程同时弱化, 并强调基函数在有限元空 间中的重要作用; 借助对位移和应变离散精度的控制, 拟协调元保障了单元的收敛性, 并可以利用泰勒展开校 核进行简便直接的收敛性分析. 研究者们利用拟协调元已经构造了大量的优秀的单元, 并广泛地应用到结构问 题、流体流动问题、非线性分析、稳定性和破坏分析等方面. 这些工作集中体现了拟协调元的理论价值和工程 应用价值. 对拟协调列式方法、列式理论和已发表文献中的主要拟协调单元进行了总结. 最后对拟协调的研究 工作进行了展望.      

求解固体力学问题的强?弱耦合形式单元微分法

单元微分法是一种新型强形式有限单元法. 与弱形式算法相比, 该算法直接对控制方程进行离散, 不需要用到数值积分. 因此该算法有较简单的形式, 并且其在计算系数矩阵时具有极高的效率. 但作为一种强形式算法, 单元微分法往往需要较多网格或者更高阶单元才能达到满意的计算精度. 与此同时, 对于一些包含奇异点的模型, 如在多材料界面、间断边界条件、裂纹尖端等处, 传统单元微分法往往得不到较精确的计算结果. 为了克服这些缺点, 本文提出了将伽辽金有限元法与单元微分法相结合的强?弱耦合算法, 即整体模型采用单元微分法的同时, 在奇异点附近或某些关键部件采用有限元法. 该策略在保留单元微分法高效率与简洁形式等优点的同时, 确保了求解奇异问题的精度. 在处理大规模问题时, 针对关键部件采用有限元法, 其他部件采用单元微分法, 可以在得到较精确结果的同时, 极大提高整体计算效率. 在本文中, 给出了两个典型算例, 一个是具有切口的二维问题, 一个是复杂的三维发动机问题. 针对这两个问题, 分析了该耦合算法在求二维奇异问题和三维大规模问题时的精度与效率.      

一种有效的广义特征值分析方法

提出了一种适合于自适应有限元分析中求解广义特征值问题的多重网格方法．这种方法充分利用了初始网格下的结果，通过插值或最小二乘拟合技术来得到网格变化后的新的近似特征向量，然后由多重网格迭代过程实现对结构广义特征值问题的求解．在多重网格迭代的光滑步中，选择了收敛梯度法以提高其收敛率；在粗网格校正步中，则导出了一种近似求解特征向量误差的方程．这种方法将网格离散过程和数值求解过程很好地相结合，建立了一个网格细分后广义特征值问题的快速重分析方法，与传统有限元方法相比较，具有计算简便、计算量少等特点，可以作为结构动力问题自适应有限元分析的一种十分有效的工具．     

Analysis of bonded anisotropic wedges with interface crack under anti-plane shear loading

The antiplane stress analysis of two anisotropic finite wedges with arbitrary radii and apex angles that are bonded together along a common edge is investigated. The wedge radial boundaries can be subjected to displacement-displacement boundary condi- tions, and the circular boundary of the wedge is free from any traction. The new finite complex transforms are employed to solve the problem. These finite complex transforms have complex analogies to both kinds of standard finite Mellin transforms. The traction free condition on the crack faces is expressed as a singular integral equation by using the exact analytical method. The explicit terms for the strength of singularity are extracted, showing the dependence of the order of the stress singularity on the wedge angle, material constants, and boundary conditions. A numerical method is used for solving the resul- tant singular integral equations. The displacement boundary condition may be a general term of the Taylor series expansion for the displacement prescribed on the radial edge of the wedge. Thus, the analysis of every kind of displacement boundary conditions can be obtained by the achieved results from the foregoing general displacement boundary condition. The obtained stress intensity factors （SIFs） at the crack tips are plotted and compared with those obtained by the finite element analysis （FEA）.     

基于(α, β)变换和距离变换的三维弱奇异边界积分近奇异性消除方法

通过非线性变换求解三维弱奇异积分时,变换的雅可比消除了被积函数的奇异性。然而,当积分单元形状较差,如顶角过大或者顶角边长比过大时,弱奇异积分中的近奇异性仍然存在,这将导致弱奇异积分计算精度低甚至计算结果完全错误。因此,本文提出了一种基于(α, β)变换和距离变换的弱奇异积分中的近奇异性消除方法,用于精确计算三维弱奇异积分。首先通过(α, β)变换消除弱奇异积分中α方向的奇异性,并分离出β方向的近奇异性;然后针对β方向的积分函数形式,构造对应的距离变换来消除其近奇异性;最后给出具有大顶角和大边长比的弱奇异积分数值算例。结果表明,采用(α, β)变换和β方向距离变换相结合的方案可以精确计算不同单元形状的弱奇异积分。