Trefftz有限元法的研究进展

Trefftz有限元法(Trefftz finite element method,TFEM)是一种高效的数值计算方法,兼有传统有限元法和边界元法的诸多优点.基于双独立插值模式,结合杂交泛函和高斯散度定理,推得仅含边界积分的有限元格式.简述了过去10年间(2007—2016)Trefftz有限元法在单元域内插值函数、源项处理、特殊功能单元以及非各向同性材料等方面的研究进展,并对未来的发展趋势给出了几点展望.     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场、稳定温度场、流体绕流、介质中的渗流等各类位势问题。大量算例均获得了满意的结果。     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场，稳定温度场，流体绕流，介质中的渗流等各类位势问题，大量算例均获得了满意的结果。     

一种改进的边界元法在弹性扭转问题中的应用

本文用一种改进的边界元法分析与计算了椭圆截面等直杆的扭转问题,并与正规的边界元法的解进行比较,其结果完全一致.然而,改进边界元法较正规边界元法需要准备的数据大大减少,计算时间更加缩短.因此,本文方法对求解 Poisson 方程问题是一种经济而行之有效的数值计算方法.     

二维结构与声耦合声场问题的一种预测方法

对于二维结构与声耦合声场问题,本文提出了一种基于间接Trefftz方法的波数法,在该方法中,结构的位移和域内的声压响应分别用精确满足控制方程的与所求域大小无关的波函数通解和由外部激励的特解来近似表示.通过采用加权余量公式,强迫近似解在平均意义上满足边界条件,以此获得各个通解的系数,从而得到最终解.算例表明,波数法和边界元法的结果吻合得很好,并且在具有相同精度和收敛性的情况时,波数法比边界元法所需的自由度少.     

三维弹性问题Taylor展开多极边界元法的误差分析

Taylor展开多极边界元法有效的提高了边界元法的求解效率,使之可用于大规模问题的计算。然而,由于计算中对基本解进行了Taylor级数展开,与传统边界元方法相比计算精度有所下降。本文主要针对三维弹性问题Taylor展开多极边界元法的计算精度和误差进行研究。文中对两种方法的计算精度进行了比较;研究了核函数的Taylor展开性质;推导了三维弹性问题基本解的误差估计公式;给出了Taylor展开多极边界元法中远近场的划分原则。通过具体的算例,证明了该方法的正确性和误差估计公式的有效性,说明了影响Taylor展开多极边界元法求解精度的因素。     

二维正交各向异性位势问题的高阶单元快速多极边界元法

研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法. 在快速多极边界元法中, 源点对于远场区域的积分采用快速多极展开式计算, 而对于近场区域的积分则直接进行计算. 高阶单元的使用使得近场积分, 尤其是奇异积分和几乎奇异积分的计算更加复杂. 通过引入复数表达对其进行简化, 若边界采用线性单元插值, 近场积分可直接解析计算; 若采用二次单元插值, 则给出一个半解析算法计算近场积分. 高阶单元奇异积分和几乎奇异积分计算难题的解决, 使得高阶单元快速多极边界元法不仅能够计算一般结构, 也能被应用于超薄体结构, 拓宽了高阶单元快速多极边界元法的适用范围. 数值算例表明, 若计算精度一定, 高阶单元快速多极边界元法较常值单元快速多极边界元法使用的单元数量显著减少, 且高阶单元快速多极边界元法计算时间与自由度数量成线性关系, 其计算效率仍处于$O(N)$量级, 因此高阶单元快速多极边界元法可更加高效求解大规模问题.      

断裂力学的相似边界元法及其应用

首先对弹性力学的相似边界元法进行了研究,推导了相应的计算公式。与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,当边界单元数目较多时大大减少了计算量。在此基础上,将相似边界元法应用于断裂力学,对路面断裂力学问题进行了计算,与有限元法的结果比较,说明了本文方法在减少计算量的情况下仍能较好地保证精度。     

快速多极子展开技术在高阶边界元方法中的实现

高阶边界元法以较常数元方法计算精度高存储低而在工程计算中得到了广泛的应用,但由于其平方存储和计算量的本质,无法应用于大型工程问题中。本文将快速多极子方法（FMM）应用于高阶边界元中从而使其计算量和存储量分别降为O（Nlog N）和O（N）。通过无限区域中水流绕射算例的数值计算,对FMM高阶边界元法与传统高阶边界元法的运算速度和内存消耗进行了分析对比,结果表明对于大型计算问题FMM高阶边界元算法更加有效。     

Trefftz法中的自锁现象及变量减缩法

用Trefftz型边界解法分析了中厚板弯曲问题,发现了一类新的、因Trefftz函数溢机引起的“自锁现象”,并提出了一种消除这类自锁问题的“变量减缩法”。     

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

Axisymmetric viscoplastic deformation by the boundary element method

This paper presents a boundary element formulation and numerical implementation of the problem of small axisymmetric deformation of viscoplastic bodies. While the extension from planar to axisymmetric problems can be carried out fairly simply for the finite element method (FEM), this is far from true for the boundary element method (BEM). The primary reason for this fact is that the axisymmetric kernels in the integral equations of the BEM contain elliptic functions which cannot be integrated analytically even over boundary elements and internal cells of simple shape. Thus, special methods have to be developed for the efficient and accurate numerical integration of these singular and sensitive kernels over discrete elements. The accurate determination of stress rates by differentiation of the displacement rates presents another formidable challenge.A successful numerical implementation of the boundary element method with elementwise (called the Mixed approach) or pointwise (called the pure BEM or BEM approach) determination of stress rates has been carried out. A computer program has been developed for the solution of general axisymmetric viscoplasticity problems. Comparisons of numerical results from the BEM and FEM, for several illustrative problems, are presented and discussed in the paper. It is possible to get direct solutions for the simpler class of problems for cylinders of uniform cross-section, and these solutions are also compared with the BEM and FEM results for such cases.     

基于虚拟激励法的结构随机振动声辐射分析

本文基于有限元法、边界元法和虚拟激励法，对随机激励下结构振动声辐射问题进行研究。提出了一种计算随机激励下结构振动声辐射问题的新方法，其中，有限元法用于计算结构谐振响应，边界元法用于计算结构振动声辐射，虚拟激励法结合有限元和边界元计算随机激励下结构振动声辐射问题。 数值算例表明，本文方法在计算精度上与传统方法等价，且更具高效性。     

多重应力奇异性及其强度系数的数值分析方法

以具有两个应力奇异性次数的平面问题为例,提出了一种利用普通的数值分析结果确定奇异点附近多重应力奇异性的各阶次数以及相应的应力强度系数的数值分析方法,计算实例表明,本方法可以精确地求得各阶应力奇异性的次数,并且可以很方便地应用外插法确定出对应的应力强度系数。     

三维弹性快速多极边界元法

将静电场多极展开法和广义极小残值法结合于三维弹性问题的边界元法，使其求解的计算量及所需内存量同节点的自由度总数成正比，变革计算结构，加快求解速度以适应大规模数值计算。两者结合的关键点在于边界元法基本解的合理分解，并用广义极小残值法(GMRES)求解方程。轧机支承辊变形场大规模数值算例的总自由度数首次达N=34008并获得成功。清晰地描述了支承辊和工作辊接触区的辊型。     

Boundary element method for model analysis of 2-D composite structure

The boundary element method is used for the modal analysis of free vibration of 2-D composite structures in this paper. Since the particular solution method is used to treat the terms of body forces (inertial forces) in the equation of motion, only static fundamental solutions are needed in solving the problem. For an isotropic cantilever beam, the numerical results obtained by using the BEM presented in this paper are in good agreement, with, those of using FEM or other BEM, but this BEM can also be used to analyze problems for anisotropic materials. For simply supported composite laminated beams, the comparisons of the numerical reslts obtained by this method with the analytical results obtained by 1-D laminated beam theory indicate that if the ratio of length/thickness is greater than 20, the results of the two methods are in good agreement, but if the ratio of length/thickness is less than 20, big errors will occur for 1-D laminated beam theory.     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.     

具有中心圆孔的[0°/90°]_s复合材料层合板的边界元法分析

本文在文[1]的基础上,采用子结构法建立了多层复合板的边界元方法,对具有中心园孔[0°/90°]_s的层合板的层间应力作了计算,同有限元法的结果进行了比较,结果表明,应用边界元法处理这类问题,单元划分少,节约了计算机时,而且有较高的计算精度。     

A coupled FEM/BEM approach and its accuracy for solving crack problems in fracture mechanics

The finite element (FEM) and the boundary element methods (BEM) are well known powerful numerical techniques for solving a wide range of problems in applied science and engineering. Each method has its own advantages and disadvantages, so that it is desirable to develop a combined finite element/boundary element method approach, which makes use of their advantages and reduces their disadvantages. Several coupling techniques are proposed in the literature, but until now the incompatibility of the basic variables remains a problem to be solved. To overcome this problem, a special super-element using boundary elements based on the usual finite element technique of total potential energy minimization has been developed in this paper. The application of the most commonly used approaches in finite element method namely quarter-point elements and J-integrals techniques were examined using the proposed coupling FEM–BEM. The accuracy and efficiency of the proposed approach have been assessed for the evaluation of stress intensity factors (SIF). It was found that the FEM–BEM coupling technique gives more accurate values of the stress intensity factors with fewer degrees of freedom.     

A DOMAIN DECOMPOSITION ALGORITHM WITH FINITE ELEMENT-BOUNDARY ELEMENT COUPLING

A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.