Waves generated by an inclined-plate wave generator

This paper describes the characteristics of small-amplitude waves generated by a sinusoidally oscillating, inclined paddle-type wavemaker operating in a constant-depth channel. Two-dimensional, linearized potential flow is assumed. A semi-analytical method, the boundary collocation method, is used to establish the relationship between wave amplitude and paddle stroke. The numerical results are compared with the numerical results of the boundary integral equation method. It is found that the boundary collocation method is simpler and more flexible to implement and faster to compute. In addition, the numerical results are in reasonably good agreement with the laboratory experimental data. For the vertical wavemaker, which is a special case of the inclined wavemaker, an analytical series solution can be found. By using the boundary collocation method and the boundary integral equation method to solve the vertical wavemaker problem and comparing the results with the analytical series solution, it is found that the boundary collocation method yields a solution which is much more accurate than that from the boundary integral equation method. Finally, the relationships between wave amplitude and paddle stroke are established for different inclinations of the paddle-type wavemaker, based on the boundary collocation method.     

弹性力学平面问题的无奇异边界积分方程

将弹性力学平面问题归化成无奇异边界积分方程,避免了传统的边界元法中的柯西主值（CPV）积分和Hadamard-Finite-Parts(HFP)积分的计算,建立完整的数值求解体系。     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程，结合广义Betti-Rayleigh互易等式与时域下的边界积分方程，推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解，将经过主部分析的积分核函数分解为静态和动态部分，其中动态积分核不具有奇异性。在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换，采用配置点法计算超奇异积分，获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序，得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律，数值结果稳定且收敛速度快。     

加权残数配点法解正交各向异性板的积分方程

本文推导了一般各向异性板弯曲的积分方程,运用加权残数配点法求解了正交各向异性板弯曲的积发方程,本文将部分配点取在边界上,另一部分配点取在域外,只用关于找度的基本积分方程,而不用关于转角的补充积分方程,简化了方程求解和计算程序,由于正交各向异性板没有争析形式的、实用的基本解,本文提出了两种新的近似基本解;加权双三角级数;广义各向同性板解析形式的基本解和加权双三角级数的叠加,算例表明,本文提出的解法和近似基本解适用于各类边界条件的正交各向异性板,具有简单、可靠、精度高等优点。     

Stress analysis for a layer of finite length bonded to a half-space of identical material

The problem in the plane theory of elasticity of an elastic layer bonded to an elastic half-space of the same material is considered. The formulation is achieved by means of integral transforms and the problem is reduced to the solution of a system of singular integral equations. A numerical solution is accomplished for the half-space and layer by use of the collocation scheme developed by Erdogan and Gupta.     

Error analysis and adaptivity in three-dimensional linear elasticity by the usual and hypersingular boundary contour method

Two related topics are addressed in this article. The first part of the article proves that, for a certain admissible class of problems in linear elasticity, the hypersingular boundary contour method (HBCM) can be collocated at all boundary points on the surface of a three-dimensional (3-D) body, including those on boundary contours, edges and corners, because the HBCM-shape-functions satisfy, a priori, all the smoothness requirements for collocation at these points. In contrast, the hypersingular boundary element method needs, in general, relaxation of some of these smoothness requirements for its shape functions, even for collocation at regular points that lie on the boundaries of boundary elements.A hypersingular residual, obtained from the standard and hypersingular boundary integral equations (HBIEs), has been recently proposed as a local error estimator for a boundary element, for the boundary integral equation. The second part in the present article is concerned with a definition of an analogous local error estimator for the boundary contour method, for 3-D linear elasticity. This error estimator is then used to drive an h-adaptive meshing procedure. Numerical results are presented to demonstrate adaptive meshing for selected example problems.     

Boundary integral equations of unique solutions in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＢＥＭ）ｐｒｏｖｉｄｅｓａｎａｔｔｒａｃｔｉｖｅａｌｔｅｒｎａｔｉｖｅｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｉｔｓｍａｉｎａｄｖａｎｔａｇｅｓａｒｅｅｃｏｎｏｍｉｃａｌａｎｄｐａｒｔｉｃｕｌａｒｌｙｃｏｎｖｅｎｉｅｎｔｆｏｒｕｎｂｏｕｎｄｅｄｄｏｍａｉｎａｎｄｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｐｒｏｂｌｅｍｓ．Ｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎ（ＢＩＥ）ｉｓｔｈｅ…     

非均匀介质散射问题的体积分方程数值解法

将非均匀介质视为某一均匀背景介质中的扰动，可建立用均匀背景介质格林函数作基本解的体积分方程．给出了配置法求解体积分方程的数值方法，首先解得扰动域内各点以速度扰动为权的波场函数，然后回代计算得到观测面上各接收点的散射波场．与边界元法和Ｂｏｒｎ近似法计算结果比较表明该方法具有很高的精度，可得到穿过非     

边界节点法的计算稳定性研究

边界节点法利用满足控制方程的非奇异通解作为基函数,半解析边界数值离散偏微分方程,具有精度高、收敛快、易编程等优点,是一种纯无网格配点方法.但是在求解具体问题时,随着节点数的增加,边界节点法经常得到严重病态的插值矩阵.本文利用有效条件数评价边界节点法求解Helmholtz问题线性方程组的计算稳定性;然后利用三种正则化方法处理其病态的线性方程组,并与高斯消元法比较计算精度和收敛性.通过数值实验,本文研究了有效条件数、误差和正则化方法之间的关系.     

调和函数边界积分方程的充要条件

本文以调和函数的边值问题为例,探讨了边界积分方程的充要条件.文中首次提出了超定问题的概念,并建立了超定问题有解的一个充要条件,它也就是直接变量边界积分方程的一个充要条件.文中首次阐明了边界积分方程与变分原理的内在的联系,还指出了间接变量与直接变量两类边界积分方程之间存在着一一对应的关系.文中的慨念、思路和论点不难用于其它有变分原理的问题的边界积分方程.     

Hermite梯度重构核近似配点法及其在功能梯度材料板的静力分析中的应用

由于直接配点法在求解边值问题时边界上的求解精度较低,本文提出了Hermite梯度重构核近似配点法（HGCM）来改进边界求解精度。重构核近似是无网格法中一种常用的近似函数,但是其在求解高阶导数时格式复杂且非常耗时。HGCM采用梯度重构核近似构建形函数的任意高阶导数,提高了计算效率;通过Hermite配点法构建离散方程,提高了边界求解精度。这种方法在求解对应变系数四阶偏微分方程的功能梯度材料板的静力问题时精度高,计算效率高,并可进一步推广应用于高阶偏微分方程描述的边值问题。     

Quasi-Green＇s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green＇s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi- Green＇s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob- lem. The mode shape differential equations of the free vibration problem of a simply- supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa- tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green＇s function method.     

A boundary collocation method for the solution of a flow problem in a complex three-dimensional porous medium

A flow problem in a complex three-dimensional domain with a free surface and mixed-type boundary conditions is solved by the boundary collocation method. The solution is expressed as a combination of source functions distributed all around the domain close to the boundary, plus a special basis function to take care of a corner singularity. The resulting procedure is compared with the boundary integral elements method and is found to be simpler and more flexible to implement and faster to compute.     

Necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions

A necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions is established in the present paper. A super-determined problem of harmonic functions is proposed for the first time. Then a necessary and sufficient condition for the existence of solution of the super-determined problem is proved. At the same time, it is a necessary and sufficient condition for the correct formulation of boundary integral equations with direct unknown quantities. A relation between boundary integral equations and variational principles is discovered for the first time. And a one-to-one correspondence between boundary integral equations with direct and indirect unknown quantities is indicated. The concept and route of the present paper can be applied to other boundary value problems possessing variational principles.     

混合边界各向异性矩形板的静力分析

采用一般解析解和配点法相结合的方法,求解混合边界各向异性矩形板的弯曲问题.先由弯曲挠度的微分方程求出各种类型的齐次解和特解,然后组成一般解析解,再将板的每个边等分 为很多微小的段,仅对每一微段的中点建立应满足的边界条件,由全部边界条件方程式即可求得全部积分常数.以每边一半边界为平夹、另一半边界为简支或自由的方板为例进行了计算,并与四边均为简支的方板进行了对比,表明理论简单,结果实用.     

奇异边界法:一个新的、简单、无网格、边界配点数值方法

物理力学控制方程的基本解有源点奇异性.因而,传统的观点认为奇异基本解一般不能用做控制方程数值解的基函数;除非源点布置在物理域以外的虚假边界上,与物理边界上的配点分离.后者就是近年来受到广泛关注的基本解方法的基本思路.与这些传统方法不同,文中直接使用基本解做为插值基函数,且源点和配点是同一组物理边界上的离散点.本项工作的一个基本假设是源点奇异时的源点强度因子的存在性.利用待求问题控制方程的已知简单解,提出了一个计算源点强度因子的数值方法,并发现源点强度因子确实存在,且是一个有限值,其大小依赖于边界离散点的分布和边界条件类型.由此,文中提出了一个计算微分方程问题的新数值方法,称为奇异边界方法.该方法数学简单,编程容易,是一个真正的无网格方法.初步数值试验显示该方法精度高,收敛速度快.但有关该方法的数学物理基础还不是十分清楚.     

GREEN QUASIFUNCTION METHOD FOR FREE VIBRATION OF SIMPLY-SUPPORTED TRAPEZOIDAL SHALLOW SPHERICAL SHELL ON WINKLER FOUNDATION

The idea of Green quasifunction method is clarified in detail by considering a free vibration problem of simply-supported trapezoidal shallow spherical shell on Winkler foundation.A Green quasifunction is established by using the fundamental solution and boundary equation of the problem.This function satisfies the homogeneous boundary condition of the problem.The mode shape differential equation of the free vibration problem of simply-supported trapezoidal shallow spherical shell on Winkler foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula.There are multiple choices for the normalized boundary equation.Based on a chosen normalized boundary equation, the irregularity of the kernel of integral equations is avoided.Finally, natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations derived from the integral equations.Numerical results show high accuracy of the Green quasifunction method.     

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

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

Equilibrium of an internal transverse crack in a semiinfinite elastic body with thin coating

Static elasticity problems for a half-plane and a strip weakened by a rectilinear transverse crack are studied. In each case, the upper boundary of the body is reinforced by a flexible patch. Various versions of conditions on the lower boundary are considered in the case of the strip. The crack is maintained in the open state by distributed normal forces. The method of generalized integral transforms reduces solving the problem for the equations of equilibriumto solving a singular integral equation of the first kind with the Cauchy kernel with respect to the derivative of the crack opening function. The solutions of the integral equation are constructed by the small parameter and collocation methods for various combinations of the geometric and physical parameters of the problem, and the structure of the solutions is analyzed. The values of the stress intensity factor (SIF) near the crack vertex are obtained.     

边界积分方程中近奇异积分计算的一种变量替换法

准确估计近奇异边界积分是边界元分析中一项很重要的课题,其重要性仅次于对奇异积分的处理. 近年来已发展了许多方法,都取得了一定程度的成功,但这个问题至今仍未得到彻 底的解决. 基于一种新的变量变换的思想和观点,提交了一种通用的积分变换法, 它非常有效地改善了被积函数的震荡特性,从而消除了积分的近奇异性,在不增加计算量的情况 下, 极大地改进了近奇异积分计算的精度. 数值算例表明,其算法稳定,效率高, 并可达到很高的计算精度,即使区域内点非常地靠近边界,仍可取得很理想的结果.