二维弹性新型快速多极虚边界元的最小二乘法

在虚边界元最小二乘法的方程求解中采用新型的快速多极展开和广义极小残值法,提出了一种二维弹性新型快速多极虚边界元最小二乘法的求解思想。基于二维弹性问题原有的快速多极虚边界元最小二乘法的展开格式,通过引入对角化的概念,以更新展开传递格式;相对于原有快速多极算法,该方法可进一步提高计算效率且仍能保证具有较高的计算精度。数值算例说明了该方法的可行性、计算效率、计算精度均较高。     

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

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

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

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

薄板稳定问题的边界元法

本文基于小挠度薄板弯曲问题的基本解,建立了求解薄板稳定问题的边界积分方程,并计算了若干算例,结果表明用边界元法求解薄板的稳定问题是行之有效的.     

电磁弹性固体三维问题的虚边界元-等额配点法

依据弹性力学虚边界元法的基本思想和电磁弹性固体的基本解,提出了电磁弹性固体三维问题的虚边界元-等额配点法.该方法继承传统边界元法优点的同时,有效地避免了传统边界元法的边界积分奇异性的问题.算例表明该方法有很高的精度,是求解电磁弹性固体三维问题的一个有效的数值方法.     

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

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

《无奇异边界元法》内容简介

孙焕纯等著《无奇异边界元法》一书共有上下两篇 ,上篇阐述虚边界元法的理论、方法与应用。虚边界法有三种 :一般配点法 ;最小二乘配点法 (超额配点法 ) ;最小二乘二重积分法。*分别对弹性空间、弹性平面、薄板、薄壳问题给出了一个从弹性空间方程出发的统一的数值解法 ,抛弃了板、壳理论关于变形和应力的一切假设 ,又对位势问题、弹性平面问题等给出了边界积分方程离散化求解的系数阵元素的解析计算式。下篇针对传统边界元直接法与间接法的边界积分方程的充要性问题进行了论述 ,并对位势、弹性平面和薄板等问题建立了充要积分方程 ;其次是…     

弹性力学平面问题可靠度分析的蒙特卡罗边界元法

以样条虚边界元法作为样本试验方法,采用蒙特卡罗法进行弹性力学平面问题可靠度分析.为了提高计算效率,引入Taylor展开和Neumann展开技术,避免在大量样本计算中直接生成影响矩阵及对其进行求逆运算,降低了单次样本计算时间;同时引入重要抽样技术,在相同精度情况下减少了蒙特卡罗法的抽取样本数.算例结果表明,该文提出的Taylor-Neumann展开重要抽样蒙特卡罗样条虚边界元法具有良好的计算精度和相当高的计算效率.     

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

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

弯曲波对含多圆孔薄板的散射与动应力集中

采用复变函数法和多极坐标方法,研究了弯曲波对含有多圆孔薄板的散射问题。通过板的弯曲波动方程和内力方程的推导,求出在入射弯曲波条件下该问题的一般解的函数逼近序列和边界条件的表达式。用展开正交函数的方法将待解的问题归结为对一组无穷代数方程组的求解。最后,给出了含3圆孔薄板的孔边动应力集中系数的结果,并分析了孔间距和波数对动应力分布的影响。     

Scattering of flexural wave in a thin plate with multiple circular holes by using the multipole Trefftz method

The scattering of flexural wave by multiple circular holes in an infinite thin plate is analytically solved by using the multipole Trefftz method. The dynamic moment concentration factor (DMCF) along the edge of circular holes is determined. Based on the addition theorem, the solution of the field represented by multiple coordinate systems centered at each circle can be transformed into one coordinate system centered at one circle, where the boundary conditions are given. In this way, a coupled infinite system of simultaneous linear algebraic equations is derived as an analytical model for the scattering of flexural wave by multiple holes in an infinite plate subject to the incident flexural wave. The formulation is general and is easily applicable to dealing with the problem containing multiple circular holes. Although the number of hole is not limited in our proposed method, the numerical results of an infinite plate with three circular holes are presented in the truncated finite system. The effects of both incident wave number and the central distance among circular holes on the DMCF are investigated. Numerical results show that the DMCF of three holes is larger than that of one, when the space among holes is small and meanwhile the specified direction of incident wave is subjected to the plate.     

含有非圆形双孔的无限平板中应力的解析解研究

无限平板中含有任意形状单个孔的问题可以使用复变函数方法获得其应力解析解.对于无限平板中含有两个圆孔或两个椭圆孔的双连通域问题,也可以利用多种方法进行求解,比如双极坐标法、应力函数法、复变函数法以及施瓦茨交替法等.其中复变函数中的保角变换方法是获得应力解析解的一个重要方法.但目前尚未见到用此方法求解无限板中含有一个正方形孔和一个椭圆孔的问题.当板在无穷远处受有均布载荷和孔边作用垂直均布压力时,利用保角变换方法可以求解板中含有两个特定形状孔的问题.该方法将所讨论的区域映射成象平面里的一个圆环,其中最关键的一步是找出相应的映射函数.基于黎曼映射定理,提出了该映射函数一般形式,并利用最优化方法,找到了该问题的具体映射函数,然后通过孔边应力边界条件建立了求解两个解析函数的基本方程,获得了该问题的应力解析解.运用ANSYS有限单元法与结果进行了对比.研究了孔距、椭圆形孔大小和两孔布置方位对边界切向应力的影响,以及不同载荷下两孔中心线上应力分布规律.      

A new analytical-numerical method for calculating interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses

Based on the elementary solutions and new integral equations, a new analytical-numerical method is proposed to calculate the interacting stresses of multiple circular holes in an infinite elastic plate under both remote stresses and arbitrarily distributed stresses applied to the circular boundaries. The validity of this new analytical-numerical method is verified by the analytical solution of the bi-harmonic stress function method, the numerical solution of the finite element method, and the analytical-numerical solutions of the series expansion and Laurent series methods. Some numerical examples are presented to investigate the effects of the hole geometry parameters (radii and relative positions) and loading conditions (remote stresses and surface stresses) on the interacting tangential stresses and interacting stress concentration factors (SCFs). The results show that whether the interference effect is shielding (k <1) or amplifying (k> 1) depends on the relative orientation of holes (α) and remote stresses (σ_x^∞, σ_y^∞). When the maximum principal stress is aligned with the connecting line of two-hole centers and σ_y^∞ <0.5σ_x^∞, the plate containing two circular holes has greater stability than that containing one circular hole, and the smaller circular hole has greater stability than the bigger one. This new method not only has a simple formulation and high accuracy, but also has an advantage of wide applications over common analytical methods and analytical-numerical methods in calculating the interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses.     

Dynamic stress concentrations in thick plates with two holes based on refined theory

Based on complex variables and conformal mapping, the elastic wave scat- tering and dynamic stress concentrations in the plates with two holes are studied by the refined dynamic equation of plate bending. The problem to be solved is changed to a set of infinite algebraic equations by an orthogonM function expansion method. As examples, under free boundary conditions, the numerical results of the dynamic moment concen- tration factors in the plates with two circular holes are computed. The results indicate that the parameters such as the incident wave number, the thickness of plates, and the spacing between holes have great effects on the dynamic stress distributions. The results are accurate because the refined equation is derived without any engineering hypothese.     

Application of a new fast multipole BEM for simulation of 2D elastic solid with large number of inclusions

A fast multipole method (FMM) is applied for BEM to reduce both the operation and memory requirement in dealing with very large scale problems. In this paper, a new version of fast multipole BEM for 2D elastostatics is presented and used for simulation of 2D elastic solid with a large number of randomly distributed inclusions combined with a similar subregion approach. Generalized minimum residual method (GMRES) is used as an iterative solver to solve the equation system formed by BEM iteratively. The numerical results show that the scheme presented is applicable to certain large scale problems. The project supported by the National Nature Science Foundation of China (10172053) and the Ministry of Education     

A wideband fast multipole boundary element method for half-space/plane-symmetric acoustic wave problems

This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.     

Free vibration analysis of an elliptical plate with eccentric elliptical cut-outs

The elaborated collocation multipole method is employed to obtain a semi-analytical solution, involving proper products of angular and radial Mathieu functions, for the free flexural vibrations of a fully clamped thin elastic plate of elliptical planform containing multiple elliptical cutouts of arbitrary size, location, and orientation. The problem boundary conditions are satisfied by uniformly collocating points on the boundaries, and exactly calculating the normal derivative of plate displacement at the collocation points through use of appropriate directional derivative in each coordinate system. The multipole expansion is truncated to yield a coupled algebraic linear system of equations that is then solved for the nontrivial eigensolutions. Extensive numerical simulations present the first three calculated natural frequencies and the associated deformed mode shapes of an elliptical plate with elliptical/circular cutouts, for a wide range of plate/cutout aspect ratios, and cutout location/orientation parameters. The accuracy of solutions is checked through appropriate convergence studies, and the validity of results is established with the aid of a commercial finite element package as well as by comparison with the data available in the existing literature.     

复数矢径虚拟边界谱方法在二维空穴声辐射和散射问题中的应用

基于复数矢径虚拟边界积分法,通过将虚拟积分曲线上的未知源强密度函数用Fourier级数展开,同时借助快速数值Fourier变换计算程序,提出了一种求解二维任意形状空穴声辐射和散射问题的复数矢径虚拟边界谱方法．该方法具有以下特点：(1)不存在奇异积分处理;(2)采用复数矢径虚拟边界积分方法,不仅保证了解的唯一性,而且由于虚拟源强密度函数采用Fourier级数展开,克服了用单元离散方法不能用于较高频率范围的缺点;(3)采用快速数值Fourier变换技术使计算效率大幅度提高．文中给出的计算结果表明：在求解任意形状二维空穴声辐射和散射问题上较通常采用的FEM、BEM和VBEM更为有效．