三维有限体中平片裂纹的超奇异积分方程方法——Ⅰ、理论分析

本文利用单个平片裂纹的基本解,将三维有限体中的平片裂纹问题,归为解一组超奇异积分方程,然后使用主部分析方法,对这组方程的求解作了理论分析,其结果在本文的第Ⅰ部分给出,关于这组方程的数值法求解,则给出于本文的第Ⅱ部分。     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

Characterization of the variable exponent Bessel potential spaces via the Poisson semigroup

Under the standard assumptions on the variable exponent p(x) (log- and decay conditions), we give a characterization of the variable exponent Bessel potential space Bα[Lp(⋅)(Rn)] in terms of the rate of convergence of the Poisson semigroup Pt. We show that the existence of the Riesz fractional derivative Dαf in the space Lp(⋅)(Rn) is equivalent to the existence of the limit . In the pre-limiting case we show that the Bessel potential space is characterized by the condition ‖α(I−Pε)f_‖p(⋅)≦Cεα.     

Adaptive-Additive Multilevel Methods for Hypersingular Integral Equation

In this article we present the combined adaptive-additive multilevel methods for the Galerkin approximation of hypersingular integral equation on the interval o = ( m 1,1). We also derive an efficient and reliable a posteriori error estimate for the error between the exact solution u and the approximated multilevel solution $ \tilde u_ {\cal M} $ , measuring locally the quality of $ \tilde u_ {\cal M} $ . The algorithm is carefully designed to obtain minimal complexity. A limitation of our analysis approach is that the meshes must be assumed to be quasi-uniform.     

Analysis of Iterative Sub-structuring Techniques for Boundary Element Approximation of the Hypersingular Operator in Three Dimensions

The article deals with the analysis of Additive Schwarz preconditioners for the h -version of the boundary element method for the hypersingular integral equation on surfaces in three dimensions. The first preconditioner consists of decomposing into local spaces associated with the subdomain interiors, supplemented with a wirebasket space associated with the subdomain interfaces. The wirebasket correction only involves the inversion of a diagonal matrix, while the interior correction consists of inverting the sub-blocks of the stiffness matrix corresponding to the interior degrees of freedom on each subdomain. It is shown that the condition number of the preconditioned system grows at most as max K H m 1 (1 + log H / h K ) 2 where H is the size of the quasi-uniform subdomains and h K is the size of the elements in subdomain K . A second preconditioner is given that incorporates a coarse space associated with the subdomains. This improves the robustness of the method with respect to the number of subdomains: theoretical analysis shows that growth of the condition number of the preconditioned system is now bounded by max K (1 + log H / h K ) 2 .     

A NUMERICAL BOUNDARY EIGENVALUE PROBLEM FOR ELASTIC CRACKS IN FREE AND HALF SPACE

We present in this paper a numerical method for hypersingular boundary integral equations.This method was developed for planar crack problems:additional edge si...     

Water wave scattering by a submerged circular-arc-shaped plate

The problem of water wave scattering by a thin circular-arc-shaped plate submerged in infinitely deep water is investigated by linear theory. The circular-arc is not necessarily symmetric about the vertical through its center. The problem is formulated in terms of a hypersingular integral equation for a discontinuity of the potential function across the plate. The integral equation is solved approximately using a finite series involving Chebyshev polynomials of the second kind. The unknown constants in the finite series are determined numerically by using the collocation and the Galerkin methods. Both the methods ultimately produce very accurate numerical estimates for the reflection coefficient. The numerical results are depicted graphically against the wave number for a variety of configurations of the arc. Some results are compared with known results available in the literature and good agreement is achieved. The suitability of using a circular-arc-shaped plate as an element of a water wave lens has also been discussed on the basis of the present numerical results.     

Adaptive meshless Galerkin boundary node methods for hypersingular integral equations

Adaptive refinement techniques are developed in this paper for the meshless Galerkin boundary node method for hypersingular boundary integral equations. Two types of error estimators are derived. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two consecutive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the numerical solution itself and its projection. These error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization scheme is presented to accomodate the non-local property of hypersingular integral operators for the needed computable local error indicators. The convergence of the adaptive meshless techniques is verified theoretically. To confirm the theoretical results and to show the efficiency of the adaptive techniques, numerical examples in 2D and 3D with high singularities are provided.     

THE APPROXIMATE COMPUTATION OF HYPERSINGULAR INTEGRALS ON INTERVAL

In this paper some numerical methods for computing hypersingular integrals on interval are given. Using geometric meshes, these methods lead 10 an exponential convergence in the range of engineering compulation. A numerical example shows their effeclivity and accuracy.