改进的奇异边界法模拟三维位势问题

奇异边界法是与基本解法相对应的一种边界型无网格数值离散方法. 该方法提出了源点强度因子的概念, 克服了传统基本解方法中最复杂最头疼的虚拟边界问题.基于边界元法中处理奇异积分的数值处理技术, 导出了源点强度因子的解析表达式, 提出了改进的无网格奇异边界法, 并进一步将该方法应用于三维位势问题. 该方法消除了传统方法中样本点的选取, 在不增加计算量的前提下, 极大地提高了奇异边界法的计算精度与稳定性.      

弹性动力学高阶核无关快速多极边界元法

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nystr?m 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nystr?m 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nystr?m 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.      

A HIGH ORDER KERNEL INDEPENDENT FAST MULTIPOLE BOUNDARY ELEMENT METHOD FOR ELASTODYNAMICS

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nyström 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nyström 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nyström 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.     

一种模拟大规模高频声场的双层奇异边界法

大规模高频声场的数值模拟是一项非常有计算挑战性的课题. 为了解决传统边界型离散方法由于全局支撑的满阵限制, 不易应用于大规模高频声场模拟的计算瓶颈, 本文提出了一种用于模拟大规模高频声场的双层奇异边界法. 在该方法中, 通过引入双层结构, 细网格上的全局支撑的满阵被转化为局部支撑的大规模稀疏矩阵, 传统奇异边界法模拟大规模问题时所面临的高计算量以及过度存储需求遂得以解决. 其次, 双层奇异边界法仅通过粗网格评估远场作用, 且独立于特定的插值核函数. 相较于快速多级方法, 该方法具有更强的适应性和灵活性, 且多层结构使该方法具有一定的预调节作用, 非常适合求解具有大规模、高秩、高条件数特点的高频波矩阵. 在其后的散射球模型算例中, 双层奇异边界法配置10万个节点, 成功模拟了无量纲波数高达160的声散射问题. 在对于人头模型的声散射特性分析中, 双层奇异边界法比COMSOL软件计算速度快了约78.13\mathrm{\%}. 当配置8万个节点时, 双层奇异边界法成功模拟了频率高达25 kHz 的工况, 该频率已远远超出了人耳的听力极限.      

PROGRAM-PATTERN MULTIPOLE BOUNDARY ELEMENT METHOD FOR FRICTIONAL CONTACT

A mathematical program is proposed for the highly nonlinear problem involving frictional contact. A program-pattern using the fast multipole boundary element method (FMBEM) is given for 3-D elastic contact with friction to replace the Monte Carlo method. A newoptimized generalized minimal residual (GMRES) algorithm is presented. Numerical examples demonstrate the validity of the program-pattern optimization model for node-to-surface contact with friction. The GMRES algorithm greatly improves the computational efficiency.     

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

研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法.在快速多极边界元法中,源点对于远场区域的积分采用快速多极展开式计算,而对于近场区域的积分则直接进行计算.高阶单元的使用使得近场积分,尤其是奇异积分和几乎奇异积分的计算更加复杂.通过引入复数表达对其进行简化,若边界采用线性单元插值,...     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。     

GENERAL BOUNDARY ELEMENT METHOD FOR NON-LINEAR PROBLEMS

In this paper the well-known non-linear equation f‴+½ff″=0 with boundary conditions f(0)=0, f′(0)=0 and f (∞)=1 is used as an example to describe the basic ideas of a kind of general boundary element method for non-linear problems whose governing equations and boundary conditions may not contain any linear terms at all.     

MESHLESS METHOD OF DUAL RECIPROCITY HYBRID RADIAL BOUNDARY NODE METHOD FOR ELASTICITY

Combining the radial point interpolation method （RPIM）, the dual reciprocity method （DRM） and the hybrid boundary node method （HBNM）, a dual reciprocity hybrid radial boundary node method （DHRBNM） is proposed for linear elasticity. Compared to DHBNM, RPIM is exploited to replace the moving least square （MLS） in DHRBNM, and it gets rid of the deficiency of MLS approximation, in which shape functions lack the delta function property, the boundary condition can not be applied easily and directly and it＇s computational expense is high. Besides, different approximate functions are discussed in DRM to get the interpolation property, in which the accuracy and efficiency for different basis functions are compared. Then RPIM is also applied in DRM to replace the conical function interpolation, which can greatly improve the accuracy of the present method. To demonstrate the effectiveness of the present method, DHBNM is applied for comparison, and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.     

功能梯度材料平面问题的辛弹性力学解法

将辛弹性力学解法推广用于功能梯度材料平面问题的分析,考虑沿长度方向弹性模量为指数函数变化而泊松比为常数的矩形域平面弹性问题,给出了具体的求解步骤. 提出了移位Hamilton矩阵的新概念,建立起相应的辛共轭正交关系;导出了对应特殊本征值的本征解,发现材料的非均匀特性使特殊本征解的形式发生明显的变化.      

一种边界元分析接触问题的快速迭代算法

弹性摩擦接触问题的接触状态和接触应力的分析需要进行迭代，迭代过程是很费机时的，采用一种快速迭代算法进行迭代计算，其迭代的计算机ＣＰＵ时间约为常用迭代方法时间的１／１０左右。     

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

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

THE MESHLESS VIRTUAL BOUNDARY METHOD AND ITS APPLICATIONS TO 2D ELASTICITY PROBLEMS

A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network (IRBFN) constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations (BIE) can be established while the fundamental solutions to the problems are introduced. The singular value decomposition (SVD) method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.     

大规模边界元模态分析的高效数值方法

随着大规模快速边界元计算技术的发展,在复杂结构的动态设计、振动与噪声分析中愈来愈多地采用边界元法,因此求解大规模边界元特征值问题、进行复杂结构和声场模态分析,成为工程应用中一个十分重要,但却极具挑战性的课题,目前国际上还没有十分有效的数值方法.本文针对边界元法中典型的非线性特征值问题,提出了一种通用、高效的数值解法,称为基于预解矩阵采样的Rayleigh-Ritz投影法,记为RSRR.首先,通过求解一系列频域边界元问题来构造特征向量搜索空间,进而可以采用Rayleigh-Ritz投影,将原问题转化为一个可以采用现有方法求解的小规模缩减特征值问题;其次,为了降低Rayleigh-Ritz投影过程的计算量,基于解析函数的Cauchy积分公式,构造了边界元系数矩阵的插值近似方法,以及缩减特征值问题系数矩阵的快速计算方法,给出了插值项数的估计策略;最后,将RSRR与声学快速边界元法结合,应用于大规模吸声结构的复模态分析.数值算例表明,RSRR方法能够可靠地求出给定频段内的全部特征值和特征向量,具有计算效率高、精度高、通用等优点.     

弹性力学平面问题的位移型解答

本文证明了一个线性常系数偏微分方程的通解定理,利用这个通解定理导出了弹性力学平面问题的位移通解。     

计算含大位移动边界非定常流动的无网格算法

研究了可以计算处理包含大位移动边界非定常流动的无网格算法. 创建了无网格方法中由于大位移边界运动造成的不合要求点的判断标准, 采用高效的填充方法实现点云重构, 运用线性插值方法得到新点参数, 实现了局部重构处理大位移动点. 流场计算方面, 在计算域自动布点基础上, 采用曲线逼近计算导数及HLLC格式计算数值通量, 发展了求解基于无网格的ALE方程组的算法. 最后, 对活塞问题进行了模拟, 结果与解析解相吻合, 验证了算法的准确性; 另外, 计算比较了气流流过静止圆柱以及圆柱在静止流场运动流场, 结果表明方法是成功的.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

SINGULAR PERTURBATION OF NONLINEAR BOUNDARY VALUE PROBLEMS

In this paper we consider the boundary value problem where ε.μ are two positive parameters. Under f_y≤-k<0 and other suitable restrictions, there exists a solution and it satisfied where y_(0,0)(x) is solution of reduced problem while y_i-j,j(x)(j=0,1,...,i;i=1,2,...,m) can be obtained successively from certain linear equations.     

域外奇点法在弹性力学中的应用

本文将域外奇点法推广应用于弹性力学问题,求解了扭转问题,平面弹性问题和平面弯曲问题.计算表明,这种方法简单,计算时间短,精度高。