奇异边界法模拟二维弹性力学问题

奇异边界法是一种新的边界型无网格数值离散方法.该方法使用基本解作为插值基函数,在继承传统边界型方法优点的同时,不需要费时费力的网格划分和奇异积分,数学简单,编程容易,是一个真正的无网格方法.为避免配置点与插值源点重合时带来的基本解源点奇异性,该方法提出了源点强度因子的概念,从而将边界型强格式方法的核心归结为求解源点强度因子.论文首次将该方法应用于求解平面弹性力学问题.数值算例表明,本文算法稳定,效率高,并可达到很高的计算精度.     

Singular boundary method for solving plane strain elastostatic problems

This study documents the first attempt to apply the singular boundary method (SBM), a novel boundary only collocation method, to two-dimensional (2D) elasticity problems. Unlike the method of fundamental solutions (MFS), the source points coincide with the collocation points on the physical boundary by using an inverse interpolation technique to regularize the singularity of the fundamental solution of the equation governing the problems of interest. Three benchmark elasticity problems are tested to demonstrate the feasibility and accuracy of the proposed method through detailed comparisons with the MFS, boundary element method (BEM), and finite element method (FEM).     

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

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

Singular boundary method for 2D dynamic poroelastic problems

This paper extends a strong-form meshless boundary collocation method, named the singular boundary method (SBM), for the solution of dynamic poroelastic problems in the frequency domain, which is governed by Biot equations in the form of mixed displacement–pressure formulation. The solutions to problems are represented by using the fundamental solutions of the governing equations in the SBM formulations. To isolate the singularities of the fundamental solutions, the SBM uses the concept of the origin intensity factors to allow the source points to be placed on the physical boundary coinciding with collocation points, which avoids the auxiliary boundary issue of the method of fundamental solutions (MFS). Combining with the origin intensity factors of Laplace and plane strain elastostatic problems, this study derives the SBM formulations for poroelastic problems. Five examples for 2D poroelastic problems are examined to demonstrate the efficiency and accuracy of the present method. In particular, we test the SBM to the multiply connected domain problem, the multilayer problem and the poroelastic problem with corner stress singularities, which are all under varied ranges of frequencies.     

双调和方程奇异边界法的改进及在Kirchhoff板弯曲中的应用

奇异边界法(SBM)是一种基于边界离散的无网格数值方法,在很多科学计算和工程领域中得到广泛的应用．该方法在处理复杂几何区域或者多连通区域时比基本解方法(MFS)数值计算更为稳定,具有易于实施、精度高等优点．SBM数值计算的关键之处在于源强度因子的计算,特别是相对于Laplace方程更为复杂的双调和方程的边界条件下源强度因子的计算．在高阶导数边界条件下,采用反插或者“加减项”原理计算源强度因子相对繁琐．本文对双调和方程的SBM进行了改进,将其中一个插值基函数改进为非奇异基函数形式,避免计算该基函数的源强度因子,极大简化了SBM的数值计算．本文改进对MFS同样有效,可以作为对传统MFS数值算法的补充．数值算例结果表明,本文提出的改进均能得到误差很小的数值解,且算法稳定,计算效率较高．     

压电材料平面问题的虚边界元-等额配点解法

利用压电材料平面问题的基本解和弹性力学虚边界元方法的基本思想,提出了压电材料平面问题的虚边界元-等额配点解法。该解法继承了传统边界元方法的优点,而避免了传统边界元方法遇到的边界积分奇异性问题。最后给出了压电材料平面问题的一些具体算例,并与解析解作了比较。结果表明本文的方法有很高的精度,是该问题一个十分有效的数值求解方法。     

基于时间依赖基本解的奇异边界法模拟二维狄利克雷边界标量波方程

奇异边界法是一个半解析边界配点强格式方法,具有无数值积分和无网格、编程容易以及数学简单等优点。本文首次将时间依赖基本解运用于奇异边界法,计算模拟二维标量波方程;结合确定源点强度因子的反插值技术,提出了二维狄利克雷边界标量波方程源点强度因子的一个经验公式;引进了解决波方程基本解G奇异性的一种无奇异积分处理方法。数值实验证明,基于时间依赖基本解的奇异边界法可精确高效地模拟二维狄利克雷边界标量波方程,在计算效率、精度、稳定性和适应性等方面有明显优势。     

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.     

On new symplectic elasticity approach for exact bending solutions of rectangular thin plates with two opposite sides simply supported

This paper presents a bridging research between a modeling methodology in quantum mechanics/relativity and elasticity. Using the symplectic method commonly applied in quantum mechanics and relativity, a new symplectic elasticity approach is developed for deriving exact analytical solutions to some basic problems in solid mechanics and elasticity which have long been bottlenecks in the history of elasticity. In specific, it is applied to bending of rectangular thin plates where exact solutions are hitherto unavailable. It employs the Hamiltonian principle with Legendre’s transformation. Analytical bending solutions could be obtained by eigenvalue analysis and expansion of eigenfunctions. Here, bending analysis requires the solving of an eigenvalue equation unlike in classical mechanics where eigenvalue analysis is only required in vibration and buckling problems. Furthermore, unlike the semi-inverse approaches in classical plate analysis employed by Timoshenko and others such as Navier’s solution, Levy’s solution, Rayleigh–Ritz method, etc. where a trial deflection function is pre-determined, this new symplectic plate analysis is completely rational without any guess functions and yet it renders exact solutions beyond the scope of applicability of the semi-inverse approaches. In short, the symplectic plate analysis developed in this paper presents a breakthrough in analytical mechanics in which an area previously unaccountable by Timoshenko’s plate theory and the likes has been trespassed. Here, examples for plates with selected boundary conditions are solved and the exact solutions discussed. Comparison with the classical solutions shows excellent agreement. As the derivation of this new approach is fundamental, further research can be conducted not only on other types of boundary conditions, but also for thick plates as well as vibration, buckling, wave propagation, etc.     

An elasticity solution of a nonhomogeneous half-plane problem

ＡＮＥＬＡＳＴＩＣＩＴＹＳＯＬＵＴＩＯＮＯＦＡＮＯＮＨＯＭＯＧＥＮＥＯＵＳＨＡＬＦ－ＰＬＡＮＥＰＲＯＢＬＥＭＷｏＧｕｏ－ｗｅｉ（沃国纬）（ＳｈａｎｇｈａｉＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ．Ｓｈａｎｇｈａｉ）（ＲｅｃｅｉｖｅｄＪａｎ．５．１９９４：...     

The analysis of crack problems with non-local elasticity

IntroductionTheclassicalconhnuummechanicshasbeenusedtosolvemanyproblemsinmacrofracturemechanics,butencountersdifficulheswhentheeffectofITilcrocharacteristicdimensionshouldbetakenintoaccount.Thestressfieldverynearthecracktipisstillnotclear.Somephenomenaofshortcrackscannotbeexplained["']andsomemechanismoffracturehasnotbeensolvedyet.Thenon-localelashcitytheoryseemsattractivetotheseproblems.Thetheoryofnon-localelasticity,establishedanddevelopedbyEringenetal[3),connectstheclassicalcontinuummechan…     

Relaxation procedures for an iterative MFS algorithm for the stable reconstruction of elastic fields from Cauchy data in two-dimensional isotropic linear elasticity

We investigate two numerical procedures for the Cauchy problem in linear elasticity, involving the relaxation of either the given boundary displacements (Dirichlet data) or the prescribed boundary tractions (Neumann data) on the over-specified boundary, in the alternating iterative algorithm of Kozlov et al. (1991). The two mixed direct (well-posed) problems associated with each iteration are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method, while the optimal value of the regularization parameter is chosen via the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point where the accumulation of noise becomes dominant and the errors in predicting the exact solutions increase, is also presented. The MFS-based iterative algorithms with relaxation are tested for Cauchy problems for isotropic linear elastic materials in various geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the proposed method.     

三维高频声波的矩阵压缩边界节点法模拟

本文采用边界节点法(Boundary Knot Method, BKM)求解三维高频声场．由于高频赫姆霍兹方程的解是振荡的，极大影响了数值求解的精确度，需要在计算区域增加离散点，这会增加计算量．同时对于大规模声学问题，依靠边界节点法形成的插值矩阵为满秩，导致计算量过高和存储量过大．所以本文采用矩阵压缩技术(Matrix Compression, MC)，在有效继承边界节点法高精确度的基础上减少计算内存需求和时间，从而提高计算效率．数值实验表明，MC-BKM 求解精度高、收敛速度快、计算时间少，在高频大规模声波问题中应用前景广泛．     

Dimension Reduction in the Context of Structured Deformations

In the present paper the basic boundary value problems (BVPs) of the full coupled linear theory of elasticity for triple porosity materials are investigated by means of the potential method (boundary integral equation method) and some basic results of the classical theory of elasticity are generalized. In particular, the Green’s identities and the formula of Somigliana type integral representation of regular vector and regular (classical) solutions are presented. The representation of Galerkin type solution is obtained and the completeness of this solution is established. The uniqueness theorems for classical solutions of the internal and external BVPs are proved. The surface (single-layer and double-layer) and volume potentials are constructed and their basic properties are established. Finally, the existence theorems for classical solutions of the BVPs are proved by means of the potential method and the theory of singular integral equations.     

The method of analysis of crack problem in three-dimensional non-local elasticity

ｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｈｅｏｒｅｔｉｃａｌｆｒａｍｅｏｆｎｏｎ＿ｌｏｃａｌｅｌａｓｔｉｃｉｔｙｈａｓｂｅｎｓｅｔｕｐ［１,２］．Ｍａｎｙｉｍｐｏｒｔａｎｔｒｅｓｕｌｔｓｈａｖｅｂｅｎｇｏｔｂｙｕｓｉｎｇｎｏｎ＿ｌｏｃａｌｔｈｅｏｒｙｔｏｓｏｌ...     

A refined theory of elastic thick plates for extensional deformation

Based on elasticity theory, various two-dimensional (2D) equations and solutions for extensional deformation have been deduced systematically and directly from the three-dimensional (3D) theory of thick rectangular plates by using the Papkovich–Neuber solution and the Lur’e method without ad hoc assumptions. These equations and solutions can be used to construct a refined theory of thick plates for extensional deformation. It is shown that the displacements and stresses of the plate can be represented by the displacements and transverse normal strain of the midplane. In the case of homogeneous boundary conditions, the exact solutions for the plate are derived, and the exact equations consist of three governing differential equations: the biharmonic equation, the shear equation, and the transcendental equation. With the present theory a solution of these can satisfy all the fundamental equations of 3D elasticity. Moreover, the refined theory of thick plate for bending deformation constructed by Cheng is improved, and some physical or mathematical explanations and proof are provided to support our justification. It is important to note that the refined theory is consistent with the decomposition theorem by Gregory. In the case of nonhomogeneous boundary conditions, the approximate governing differential equations and solutions for the plate are accurate up to the second-order terms with respect to plate thickness. The correctness of the stress assumptions in the classic plane-stress problems is revised. In an example it is shown that the exact or accurate solutions may be obtained by applying the refined theory deduced herein.     

Two-sided estimates in unilateral elasticity

The paper presents an extension to unilateral problems of the classical method of bounding (above and below) the solutions of linear self-adjoint boundary value problems. Using this extension the solution of the general unilateral problem in linear elasticity is bounded in energy by two suitable defined admissible states belonging to two complementary convex sets.     

Potential Method in the Linear Theory of Viscoelastic Materials with Voids

In the present paper the linear theory of viscoelasticity for Kelvin–Voigt materials with voids is considered and some basic results of the classical theory of elasticity are generalized. Indeed, the basic properties of plane harmonic waves are established. The explicit expression of fundamental solution of the system of equations of steady vibrations is constructed by means of elementary functions. The Green’s formulas in the considered theory are obtained. The uniqueness theorems of the internal and external basic boundary value problems (BVPs) are proved. The representation of Galerkin type solution is obtained and the completeness of this solution is established. The formulas of integral representations of Somigliana type of regular vector and regular (classical) solution are obtained. The Sommerfeld-Kupradze type radiation conditions are established. The basic properties of elastopotentials and singular integral operators are given. Finally, the existence theorems for classical solutions of the internal and external basic BVPs of steady vibrations are proved by using of the potential method (boundary integral method) and the theory of singular integral equations.     

COMPLEX FUNDAMENTAL SOLUTIONS FOR SEMI-INFINITE PLANE AND INFINITE PLANE WITH HOLE UNDER VARIOUS BOUNDARY CONDITIONS

I.IntroductionBoundaryelementmethod(BEM)hasbeenwidelyusedforlinearandnonlinearanalysisofengineeringstructures.ThetriditionalboundaryelementdiscretizationforthecouplingproblemoftbullditliollandpipeinFig.IisshowninFig.2.ConlparcdwithFEM,BEMIlasmoreadvantages.Inthispaper,thegeneralmethodoffindingthecomplexfundamentalsolutions(CFSs)I'orsemi-illfiniteplanealldintllliteplanewithholeundervariousboundaryconditions(aCs)hasbeenestablishedbyLlsillgRiemann-Schwarzsymmetricprincipleandsuperpositio…     

An asymptotic analysis of composite beams with kinematically corrected end effects

A finite element-based beam analysis for anisotropic beams with arbitrary-shaped cross-sections is developed with the aid of a formal asymptotic expansion method. From the equilibrium equations of the linear three-dimensional (3D) elasticity, a set of the microscopic 2D and macroscopic 1D equations are systematically derived by introducing the virtual work concept. Displacements at each order are split into two parts, such as fundamental and warping solutions. First we seek the warping solutions via the microscopic 2D cross-sectional analyses that will be smeared into the macroscopic 1D beam equations. The variations of fundamental solutions enable us to formulate the macroscopic 1D beam problems. By introducing the orthogonality of asymptotic displacements to six beam fundamental solutions, the end effects of a clamped boundary are kinematically corrected without applying the sophisticated decay analysis method. The boundary conditions obtained herein are applied to composite beams with solid and thin-walled cross-sections in order to demonstrate the efficiency and accuracy of the formal asymptotic method-based beam analysis (FAMBA) presented in this paper. The numerical results are compared to those reported in literature as well as 3D FEM solutions.