轴对称Poisson方程的Trefftz有限元解法

将径向基函数应用到一类轴对称Poisson方程的数值求解中,提出了一种Trefftz有限元计算格式.非0右端项将问题的特解引入Trefftz单元域内场,致使单元刚度方程涉及区域积分.利用径向基函数对特解近似处理,可消除区域积分,从而保持Trefftz有限元法只含边界积分的优势.为获得特解,选取求解域内所有单元的节点和形心作为基本插值点,而在求解域之外构造一个虚拟边界,在其上布置一定数目的虚拟点作为额外插值点.数值算例验证了该方法的有效性和可行性.     

含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．     

非齐次弹性力学问题双互易边界元方法研究北大核心CSCD

基于弹性力学边界元方法理论,将边界元法与双互易法结合,采用指数型基函数对非齐次项进行插值得到双互易边界积分方程.将边界积分方程离散为代数方程组,利用已知边界条件和方程特解求解方程组,得出域内位移和边界面力.指数型基函数的形状参数是由插值点最近距离的最小值决定,采用这种形状参数变化方案,分析径向基函数(RBF)插值精度以及插值稳定性.再次将指数型基函数应用到双互易边界元法中,分析双互易边界元方法下计算精度及稳定性,验证了指数型插值函数作为双互易边界元方法的径向基函数解决弹性力学域内体力项问题的有效性.     

三维稳态对流扩散问题的无网格求解算法研究

无网格法是一种不需要生成网格就可模拟复杂形状流场计算的流体力学问题求解算法.为了提高基于Galerkin弱积分形式的无网格方法求解三维稳态对流扩散问题的计算效率,提出了在空间离散上采用基于凸多面体节点影响域的无网格形函数,并通过选取适当节点影响半径因子避免节点搜索问题,同时减少系统刚度矩阵带宽.计算中当节点影响因子为1.01时,无网格方法的形函数近似具有插值特性且本质边界条件的施加与有限元一样简单.三维立方体区域的稳态对流扩散数值算例表明:在保证计算精度的同时,采用凸多面体节点影响域的无网格方法比传统无网格方法最高可节省计算时间42%.因此从计算效率和精度考虑,在运用无网格方法求解三维问题时建议采用凸多面体节点影响域的无网格方法.     

弹性动力问题的有限元与特解边界元耦合的数值方法

导出了特解边界元法与有限元法的耦合方程。并应用自由度缩减技术,使耦合方程的自由度缩减到有限元域及其和边界元域的耦合边界上。这样得到的耦合方程不增加原有限元方程的带宽和阶数。耦合方程的求解可以引用求解有限元方程的所有方法,易于程序实现。数值算例结果表明,本文所提出的方法是正确的,是一种较为理想的耦合方法。     

基于参数化水平集法的材料非线性子结构拓扑优化

针对利用传统水平集法进行非线性结构拓扑优化计算过程复杂及计算效率低等问题,将参数化水平集方法引入材料非线性结构拓扑优化中。通过全局径向基函数插值初始水平集函数,建立了以插值系数为设计变量、结构的应变能最小为目标函数、材料用量为约束条件的材料非线性结构拓扑优化模型,利用有限元分析对材料非线性结构建立平衡方程,并用迭代法求解。同时,采用子结构法划分设计区域为若干个子区域,将全自由度平衡方程的求解分解为缩减的平衡方程和多个子结构内部位移的求解,减小了计算成本。算例表明,这种处理非线性关系的方法可以在保证数值稳定的同时提高计算效率,得到边界清晰、结构合理的拓扑优化构形。     

插值型无单元Galerkin比例边界法与有限元法的耦合在压电材料断裂分析中的应用

插值型无单元Galerkin比例边界法是一种只需在边界上采用插值型无单元Galerkin法离散且无需基本解的半解析方法,能有效求解压电材料的断裂问题.为进一歩提高这种方法的适用性,该文提出了一种用于压电材料断裂分析的插值型无单元Galerkin比例边界法耦合有限元法(finite element method,FEM)的分析方法.裂纹周边一定范围的计算域采用插值型无单元Galerkin比例边界法离散,其余区域采用FEM离散.插值型无单元Galerkin比例边界法方程和FEM方程的耦合可利用界面两侧广义位移的连续条件方便地实现.最后,给出了两个数值算例验证了该文所提方法的有效性.     

一种基于H-R变分的杂交广义元方法

基于Hellinger Reissner变分原理,通过构造合适的应力场函数使其能更方便和更准确地得到节点上的应力值,同时结合广义有限元构造广义位移插值的方法,在不提高单元节点数目的前提下提高位移场函数的阶次,从而提高其求解精度．这种方法能同时灵活地构造应力场和位移场,在同等精度条件下能占用较少内存和求解更少的方程数目,计算结果也显示了这种方法的有效性和很高的计算精度．     

多裂纹问题计算分析的本征COD边界积分方程方法

针对多裂纹问题，若采用常规的数值求解技术，计算效率较低.为实现多裂纹问题的大规模数值模拟，建立了本征裂纹张开位移（crack opening displacement, COD）边界积分方程及其迭代算法，并引入Eshelby矩阵的定义，将多裂纹分为近场裂纹和远场裂纹来处理裂纹间的相互影响.以采用常单元作为离散单元的快速多极边界元法为参照，对提出的计算模型和迭代算法进行了数值验证.结果表明，本征COD边界积分方程方法在处理多裂纹问题时取得较大的改进，其计算效率显著高于传统的边界元法和快速多极边界元法.     

A double-layer interpolation method for implementation of BEM analysis of problems in potential theory

A double-layer interpolation method (DLIM) is proposed to improve the performance of the boundary element method (BEM). In the DLIM, the nodes of an element are sorted into two groups: (i) nodes inside the element, called source nodes, and (ii) nodes on the vertices and edges of the element, called virtual nodes. With only source nodes, the element becomes a conventional discontinuous element. Taking into account both source and virtual nodes, the element becomes a standard continuous element. The physical variables are interpolated by continuous elements (first-layer interpolation), while the boundary integral equations are collocated at the source nodes only. We further established additional constraint equations between source and virtual nodes using a moving least-squares (MLS) approximation (second-layer interpolation). Using these constraints, a square coefficient matrix of the overall system of linear equations was finally achieved. The DLIM keeps the main advantages of MLS, such as significantly alleviating the meshing task, while providing much better accuracy than the traditional BEM. The method has been used successfully for solving potential problems in two dimensions. Several numerical examples in comparison with other methods have demonstrated the accuracy and efficiency of our method.     

A solution approach for contact problems based on the dual interpolation boundary face method

The recently proposed dual interpolation boundary face method (DiBFM) has been shown to have a much higher accuracy and improved convergence rates compared with the traditional boundary element method. In addition, the DiBFM has the ability to approximate both continuous and discontinuous fields, and this provides a way to approximate the discontinuous pressure at a contact boundary. This paper presents a solution approach for two dimensional frictionless and frictional contact problems based on the DiBFM. The solution approach is divided into outer and inner iterations. In the outer iteration, the size of the contact zone is determined. Then the elements near the contact boundary are updated to approximate the discontinuous pressure. The inner iteration is used to determine the contact state (sticking or sliding), and is only performed for frictional contact problems. To make the system of equations solvable, the contact constraints and some supplementary equations are also given. Several numerical examples demonstrate the validity and high accuracy of the proposed approach. Furthermore, due to the continuity of elements in DiBFM and the detection of the contact boundary, the pressure oscillations near the contact boundary can be treated.     

一种h型自适应有限单元

ｈ型自适应有限单元在网格局部细划时．会产生非常规节点，从而破坏了一般意义上的单元连续性假定．本文利用参照节点对非常规单元进行坐标和位移插值．为保证单元之间坐标和位移的连续性，本文提出了一组修正的形函数，常用的形函数是它的一个特例．本方法应用于有限元程序时，除形函数外无须做任何改动．算例表明水文的方法具有方法简单、精度高、自由度少、计算量小等优点．     

Large Scale Simulation with Scaled Boundary Finite Element Method

Nowadays scientific and engineering applications often require wave propagation in infinite or unbounded domains. In order to model such applications we separate our model into near-field and far-field. The near-field is represented by the well-known finite element method (FEM), whereas the far-field is mapped by a scaled boundary finite element (SBFE) approach. This latter approach allows wave propagation in infinite domains and suppresses the reflection of waves at the boundary, thus being a suitable method to model wave propagation to infinity. It is non-local in time and space. From a computational point of view, those characteristics are a drawback because they lead to storage consuming calculations with high computational time-effort. The non-locality in space causes fully populated unit-impulse acceleration influence matrices for each time step, leading to immense storage consumption for problems with a large number of degrees of freedom. Additionally, a different influence matrix has to be assembled for each time step which yields unacceptable storage requirements for long simulation times. For long slender domains, where many nodes are rather far from each other and where the influence of the degrees of freedom of those distant nodes is neglectable, substructuring represents an efficient method to reduce storage requirements and computational effort. The presented simulation with substructuring still yields satisfactory results. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An overlapping domain decomposition preconditioner for a class of discontinuous Galerkin approximations of advection-diffusion problems

We consider a scalar advection-diffusion problem and a recently proposed discontinuous Galerkin approximation, which employs discontinuous finite element spaces and suitable bilinear forms containing interface terms that ensure consistency. For the corresponding sparse, nonsymmetric linear system, we propose and study an additive, two-level overlapping Schwarz preconditioner, consisting of a coarse problem on a coarse triangulation and local solvers associated to a family of subdomains. This is a generalization of the corresponding overlapping method for approximations on continuous finite element spaces. Related to the lack of continuity of our approximation spaces, some interesting new features arise in our generalization, which have no analog in the conforming case. We prove an upper bound for the number of iterations obtained by using this preconditioner with GMRES, which is independent of the number of degrees of freedom of the original problem and the number of subdomains. The performance of the method is illustrated by several numerical experiments for different test problems using linear finite elements in two dimensions.

     

求解Brinkman方程的修正弱有限元方法

本文针对Brinkman方程引入了一种修正弱Galerkin(MWG)有限元方法.我们通过具有两个离散弱梯度算子的变分形式来逼近模型. 在MWG方法中, 分别用次数为$k$和$k-1$的不连续分段多项式来近似速度函数$u$和压力函数$p$. MWG方法的主要思想是用内部函数的平均值代替边界函数. 因此, 与WG方法相比, MWG方法在不降低准确性的同时, 具有更少的自由度, 对于任意次数不超过$k-1$ 的多项式,MWG方法均可以满足稳定性条件. MWG 方法具有高度的灵活性, 它允许在具有一定形状正则性的任意多边形或多面体上使用不连续函数. 针对$H^1$和$L^22$范数下的速度和压力近似解, 建立了最优阶误差估计. 数值算例表明了该方法的准确性, 收敛性和稳定性.     

Analysis of three-dimensional interior acoustic fields by using the localized method of fundamental solutions

The traditional method of fundamental solutions has a full interpolation matrix, and thus its solution is computationally expensive, especially for large-scale problems with complicated domains. In this paper, we make a first attempt to apply the localized method of fundamental solutions for analysis of 3D interior acoustic fields. The present method first divides the whole computational domain into some overlapping subdomains, and then expresses physical variables as linear combinations of the fundamental solution in each subdomain. Finally, the method forms a sparse and banded system matrix by satisfying both governing equations at interior nodes and boundary conditions at boundary nodes. We provide four numerical experiments to verify the accuracy and the stability of the method. Comparisons of numerical results and computational time are also made between the present method, the method of fundamental solutions, and the COMSOL software.     

Numerical solution of the system of second-order boundary value problems using the local radial basis functions based differential quadrature collocation method

In this research, we propose a numerical scheme to solve the system of second-order boundary value problems. In this way, we use the Local Radial Basis Function Differential Quadrature (LRBFDQ) method for approximating the derivative. The LRBFDQ method approximates the derivatives by Radial Basis Functions (RBFs) interpolation using a small set of nodes in the support domain of any node. So the new scheme needs much less computational work than the globally supported RBFs collocation method. We use two techniques presented by Bayona et al. (2011, 2012) [29], [30] to determine the optimal shape parameter. Some examples are presented to demonstrate the accuracy and easy implementation of the new technique. The results of numerical experiments are compared with the analytical solution, finite difference (FD) method and some published methods to confirm the accuracy and efficiency of the new scheme presented in this paper.     

Boundary element analysis of uncoupled transient thermo-elastic problems with time- and space-dependent heat sources

A boundary element method (BEM) for the analysis of two- and three-dimensional uncoupled transient thermo-elastic problems involving time- and space-dependent heat sources is presented. The domain integrals are efficiently treated using the Cartesian transformation and the radial integration methods without considering any internal cells. Similar to the dual reciprocity method (DRM), some internal points without any connectivity are considered; however, in contrast to the DRM, any arbitrary mesh-free interpolation method can be used in the present formulation. There is no need to find any particular solutions and the shape functions in the mesh-free interpolation method can be arbitrary and sufficiently complicated. Unlike the DRM, the generated system of equations contains the unknowns only on the boundary. After finding the primary unknowns on the boundary, the temperature, displacement, and stress components at all internal points can directly be found without solving any system of equations. Three examples with different forms of heat sources are presented to demonstrate the efficiency and accuracy of the proposed method. Although the proposed BEM is mathematically more complicated than domain methods, such as the finite element method (FEM), it is more efficient from a modelling viewpoint since only the surface mesh has to be generated in the presented method.