首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 125 毫秒
1.
李聪  胡斌  胡宗军  牛忠荣 《力学学报》2021,53(4):1038-1048
研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法. 在快速多极边界元法中, 源点对于远场区域的积分采用快速多极展开式计算, 而对于近场区域的积分则直接进行计算. 高阶单元的使用使得近场积分, 尤其是奇异积分和几乎奇异积分的计算更加复杂. 通过引入复数表达对其进行简化, 若边界采用线性单元插值, 近场积分可直接解析计算; 若采用二次单元插值, 则给出一个半解析算法计算近场积分. 高阶单元奇异积分和几乎奇异积分计算难题的解决, 使得高阶单元快速多极边界元法不仅能够计算一般结构, 也能被应用于超薄体结构, 拓宽了高阶单元快速多极边界元法的适用范围. 数值算例表明, 若计算精度一定, 高阶单元快速多极边界元法较常值单元快速多极边界元法使用的单元数量显著减少, 且高阶单元快速多极边界元法计算时间与自由度数量成线性关系, 其计算效率仍处于$O(N)$量级, 因此高阶单元快速多极边界元法可更加高效求解大规模问题.   相似文献   

2.
高阶边界元方法在求解波浪对海上建筑物作用问题中具有诸多优势,但由于它所形成的矩阵是一个满阵,计算量和存储量均为未知量的平方量级,很难满足如多体水动力分析等大尺度多未知量问题的计算需要。本文采用预修正快速傅里叶变换高阶边界元方法(pFFT-HOBEM),将计算量和存储量均降低到未知量的线性量级。通过对不同未知量时该方法与传统边界元法的计算量与存储量的对比,以及该方法自身各步骤计算时间的对比,研究了不同pFFT网格方案对计算量和存储量的影响,并提出了基于计算时间最小化原则的pFFT网格优化方法。采用本文方法研究了四柱结构在不同频率波浪作用下的作用力及波面分布,并对近场干涉发生时的物理现象进行了分析和讨论。  相似文献   

3.
针对边界元法中高阶单元中几乎奇异积分计算难题,解剖了二维边界元法高阶单元的几何特征,定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近,从而分离出积分核中主导的奇异函数部分,其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分,再对奇异核函数积分导出解析公式,从而建立了一种新的半解析法,用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例,采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力,并与线性边界元正则化算法结果作了比较,结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构,采用正则化算法的线性边界元分析比有限元有显著优势,而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。  相似文献   

4.
采用边界元法(BEM )求解实际工程问题时,很大一部分误差来自于离散误差。为此,本文基于Lagrange插值原理,提出了一种三维等参管单元边界元算法,该单元能很好地模拟管状结构的几何外形并对物理量进行高阶插值,大大地消除了离散误差。另外,当在边界元法中使用等参管单元时,提出了一种在等参平面内消除积分奇异性的方法。算例表明,本文算法具有划分网格少,求解精度高的优点。  相似文献   

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

6.
大型边界元方程组的并行直接分块求解算法   总被引:5,自引:0,他引:5  
针对大型边界元方程组和网络微机机群环境提出了一种并行直接分块求解算法,算法基于分块高斯-若当消去法的原理,采用内外存交互技术,并行分块消去方法,节点超行的卷帘存储方案和并行环状循环逐次修正策略,增大了解题规模,提高了计算速度。算例计算结果表明该算法具有较高的并行加速比和并行效率,适用于大型问题的边界元法求解。  相似文献   

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

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

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

10.
高效伟  胡金秀  崔苗 《力学学报》2012,44(2):361-368
提出了一种用多域边界元技术求解大型工程问题的新算法. 首先, 采用三步变量凝聚技术, 将由内部点、边界点和公共结点表述的每一子域的基本边界元代数方程表述成只有公共结点变量为未知量的代数方程, 然后, 根据公共结点的平衡方程和协调条件组集具有稀疏系数特征的总体系统方程组. 为了有效求解该系统方程组, 首次在边界元法中引进一种能有效求解大型非对称稀疏系数矩阵方程组的行消元回代法(REBSM), 该方法可在方程的每一行组形成时进行消元和回代, 当方程组组集完毕后即可得到方程的解, 不需要最后的回代过程. 因为一些项的重复计算在每一行的处理中合并掉, 因此REBSM要比传统的高斯消元法需要较少的内存, 而且计算速度具有数量级的提高, 可为边界元法求解大型工程问题提供有力的方程求解器.  相似文献   

11.
In this article, we discuss how the fast multipole method (FMM) can be implemented on modern parallel computers, ranging from computer clusters to multicore processors and graphics cards (GPU). The FMM is a somewhat difficult application for parallel computing because of its tree structure and the fact that it requires many complex operations which are not regularly structured. Computational linear algebra with dense matrices for example allows many optimizations that leverage the regular computation pattern. FMM can be similarly optimized but we will see that the complexity of the optimization steps is greater. The discussion will start with a general presentation of FMMs. We briefly discuss parallel methods for the FMM, such as building the FMM tree in parallel, and reducing communication during the FMM procedure. Finally, we will focus on porting and optimizing the FMM on GPUs.  相似文献   

12.
The singular boundary method(SBM) is a recent meshless boundary collocation method that remedies the perplexing drawback of fictitious boundary in the method of fundamental solutions(MFS). The basic idea is to use the origin intensity factor to eliminate singularity of the fundamental solution at source. The method has so far been applied successfully to the potential and elasticity problems. However, the SBM solution for large-scale problems has been hindered by the operation count of O(N~3) with direct solvers or O(N~2) with iterative solvers, as well as the memory requirement of O(N~2). In this study, the first attempt was made to combine the fast multipole method(FMM) and the SBM to significantly reduce CPU time and memory requirement by one degree of magnitude, namely, O(N). Based on the complex variable representation of fundamental solutions, the FMM-SBM formulations for both displacement and traction were presented. Numerical examples with up to hundreds of thousands of unknowns have successfully been tested on a desktop computer. These results clearly illustrated that the proposed FMM-SBM was very efficient and promising in solving large-scale plane elasticity problems.  相似文献   

13.
针对虚边界元法,引入快速多极展开和广义极小残值法(GMRES)的思想,以形成快速多极虚边界元法的求解思想,并将此方法用于含圆孔薄板有效弹性模量的模拟分析.由于本文方法采用了"源点"多极展开和"场点"局部展开的组合处理方案,从而使得原问题方程组求解的计算耗时量和储存量降至与所求问题的计算自由度数成线性比例.本文工作的研究目的在于:提高虚边界元法在普通台式机上的运算能力和拓宽虚边界元法对大规模复杂问题的求解(或数值模拟).文中给出了均布圆孔的正方形薄板和之字形分布圆孔薄板二个算例,以验证该方法的可行性,计算精度和计算效率.  相似文献   

14.
The time integration method with four-order accuracy, self-starting and implicit for the diffuse chemical reaction kinetics equation or the transient instantaneous temperature filed equation was presented. The examples show that both accuracy and stability are better than Runge-Kutta method with four-order. The coefficients of the equation are stored with sparse matrix pattern, so an algorithm is presented which combines a compact storage scheme with reduced computation cost. The computation of the competitive and consecutive reaction in the rotating packed bed, taken as examples, shows that the method is effective.  相似文献   

15.
From the formulas of the conjugate gradient, a similarity between a symmetric positive definite (SPD) matrix A and a tridiagonal matrix B is obtained. The elements of the matrix B are determinedby the parameters of the conjugate gradient. The computation of eigenvalues of A is then reduced to the case of the tridiagonal matrix B. The approximation of extreme eigenvalues of A can be obtained as a by-product in the computation of the conjugate gradient ifa computational cost of O(s) arithmetic operations is added, where s is the number of iterations This computational cost is negligible compared with the conjugate gradient. If the matrix A is not SPD, the approximation of the condition number of A can be obtained from the computation of the conjugate gradient on A~T A. Numerical results show that this is a convenient and highly efficient method for computing extreme eigenvalues and the condition number of nonsingular matrices.  相似文献   

16.
规划-迭代型弹塑性摩擦接触多极边界元法   总被引:5,自引:2,他引:3  
提出基于多极展开法的规划-迭代型的高度非线性方程的IGMRES(m)高效求解法,并建立其收敛性理论。本法适用于三维弹塑性摩擦接触多极边界元法,有效处理弹塑性摩擦接触迭代的繁杂和费时问题。数值试验证明,本求解法在确保数值精度的前提下,可减少迭代次数,显著提高计算效率。  相似文献   

17.
The method of regularized Stokeslets (MRS) uses a radially symmetric blob function of infinite support to smooth point forces and allows for evaluation of the resulting flow field. This is a common method to study swimmers at zero Reynolds number where the Stokeslet is the fundamental solution corresponding to the kernel of the single layer potential. Simulating the collective motion of N micro-swimmers using the MRS results in at least N2 pair-wise interactions. Efficient simulation of a large number of swimmers in free space is observed with the implementation of the kernel-independent fast multipole method (FMM) for radial basis functions. We illustrate the complexity of the algorithm on a simple test case where we study regularized point forces, showing that the method is of order N. Additionally, we explore accuracy in time for the MRS where the swimmers are modeled as Kirchhoff rods and the kernel-independent FMM is compared to the direct calculation using the standard MRS. Optimal hydrodynamic efficiency is also explored for different configurations of swimmers.  相似文献   

18.
The finite analytic method (FA) developed in the last decade is an effective numerical method for solving fluid flow problems. However, because of the limitation in the present computer, large round-off errors are found in calculating FA coefficients when Reynolds number is large. This paper investigates the cause of this difficulty and presents a special programming technique in making an accurate computation of FA coefficients. Then a fundamental function known as "P_e" is tabulated by the accurate computation. In practical application the interpolation technique is employed so that the FA coefficients can be obtained reliably and quickly.  相似文献   

19.
In this paper, the dual integral formulation is derived for the modified Helmholtz equation in the propagation of oblique incident wave passing a thin barrier (zero thickness) by employing the concept of fast multipole method (FMM) to accelerate the construction of an influence matrix. By adopting the addition theorem, the four kernels in the dual formulation are expanded into degenerate kernels that separate the field point and the source point. The source point matrices decomposed in the four influence matrices are similar to each other or only to some combinations. There are many zeros or the same influence coefficients in the field point matrices decomposed in the four influence matrices, which can avoid calculating the same terms repeatedly. The separable technique reduces the number of floating‐point operations from O((N)2) to O(N loga(N)), where N is the number of elements and a is a small constant independent of N. Finally, the FMM is shown to reduce the CPU time and memory requirement, thus enabling us to apply boundary element method (BEM) to solve water scattering problems efficiently. Two‐moment FMM formulation was found to be sufficient for convergence in the singular equation. The results are compared well with those of conventional BEM and analytical solutions and show the accuracy and efficiency of the FMM. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号