边裂纹柱扭转的强奇性积分方程解法*

本文利用单裂纹扭转的位错型解答,使用有限部积分的概念和方法,最后将含有单根水平裂纹的柱体扭转问题归为解一个强奇性积分方程,并为其建立了数值求解方法,文末作了若干数值例子的计算,结果令人满意.     

三维横观各向同性介质界面裂纹的边界积分方程方法

基于两相三维横观各向同性介质的基本解和Somigliana恒等式,对三维横观各向同性介质中的任意形状的平片界面裂纹,以裂纹面上的不连续位移为待求参量建立了超奇异积分_微分方程,界面平行于横观各向同性面.根据发散积分的有限部积分理论,应用积分方程方法研究得到裂纹前沿的位移和应力场的表达式、奇性指数以及应力强度因子的不连续位移表达式.在非震荡情形下,超奇异积分_微分方程退化为超奇异积分方程,与均匀介质的超奇异积分方程形式完全相同.     

三维界面裂纹的理论分析

使用界面裂纹的基本解及有限部积分的方法 ,将三维界面裂纹问题 ,归为解一组以裂纹面位移间断为未知函数的超奇异积分微分方程 ,然后为此组方程的求解作了系统的理论分析     

三维无限接合体双材料多界面裂纹的超奇异微积分方程法

利用有限部积分的概念,导出了三维无限接合体中多个界面裂纹,在任意载荷作用下的超奇异微积分方程组.数值分析中,未知的位移间断采用基本分布函数和多项式乘积的形式来近似,其中基本分布函数是根据界而裂纹应力的振荡奇异性来选取的.作为典型算例,研究了存在两个矩形界面裂纹时,裂纹之间距离、裂纹形状及双材料弹性常数对应力强度因子的影响.计算表明,应力强度因子随裂纹间的距离的增大而减小.     

一类有限部积分方程数值解的误差分析

本文研究了一类有限部积分方程数值解的误差.利用离散的极值原理和复合中点公式的超收敛性,获得了积分方程配置格式的误差分析理论,改进了有限部积分方程数值解的相关研究成果.并给出数值实验验证了理论分析的正确性.     

构成单材料裂纹和双材料界面裂纹有限应力集中的一般解析函数

对构成裂纹尖端附近有限应力集中解析函数的方法进行了综述.含裂纹平面问题的应力函数可以用无理函数和指数函数两种型式表示.对单材料裂纹,将裂纹长度作为参数,对无理函数型解析函数采用直接加权积分可以消除裂纹尖端应力的奇异性,构造有限连续的应力函数和尖劈型的张开位移函数.对指数函数型解析函数的间接积分适用于界面裂纹问题,但会使积分区间的应力分布出现正负反转和不合理的张开位移形状;结合选择不同权函数的叠加可以得到满足精度要求的有限应力集中解析函数.给出了中心裂纹和对称边裂纹在面内拉伸、剪切和弯曲等6种受力状态下的基本解.阐述了作为解析函数何以回避裂纹尖端应力奇异性的理由.     

多裂纹问题计算分析的本征COD边界积分方程方法

针对多裂纹问题，若采用常规的数值求解技术，计算效率较低.为实现多裂纹问题的大规模数值模拟，建立了本征裂纹张开位移（crack opening displacement, COD）边界积分方程及其迭代算法，并引入Eshelby矩阵的定义，将多裂纹分为近场裂纹和远场裂纹来处理裂纹间的相互影响.以采用常单元作为离散单元的快速多极边界元法为参照，对提出的计算模型和迭代算法进行了数值验证.结果表明，本征COD边界积分方程方法在处理多裂纹问题时取得较大的改进，其计算效率显著高于传统的边界元法和快速多极边界元法.     

含矩形片状裂纹的三维弹性体分析

本文分析了含有矩形片状裂纹(裂纹上下表面受均匀压力作用)的三维弹性体.借助Fourier积分变换.将问题化归为二个变数的对偶积分方程,并获得裂纹面位移和裂纹前缘应力强度因子的解析表达式.     

半无穷大裂纹端部粘聚力分析

准脆性材料裂纹端部断裂过程区粘聚力是导致非线性断裂特性的重要原因,根据准脆性材料的断裂特性,对存在粘聚力分布的半无穷大裂纹力学分析模型,由变形叠加原理得到以该粘聚应力分布为未知函数的积分方程,通过对积分方程的分析推证,得到了该分布函数解的数学结构和级数型表达式;提出了由实际裂纹张开位移,确定裂纹端部粘聚力分布函数的两种方法:其一由连续的裂纹张开位移通过积分变换求解未知函数级数展开项的系数,其二是由离散的裂纹张开位移数据通过最小二乘法确定该函数;推导出了相应方法求解未知量的代数方程,并且给出了适当的算例和讨论。     

SH波对一维六方准晶中快速传播裂纹的散射

利用积分变换技术,结合Copson方法,研究了SH波对一维六方准晶中快速传播裂纹的散射问题.通过对偶积分方程求解,得到声子场和相位子场应力、位移的解析表达式,同时给出了裂纹尖端标准动应力强度因子的定义.数值算例表明随着裂纹传播速度、裂纹长度、入射角、入射波频率增大,裂纹尖端标准动应力强度因子也在增大.当裂纹长度的增大到一定值时,裂纹的传播速度对标准动应力强度因子影响逐渐减小.研究结果在工程材料应用中有一定的参考价值.     

Theoretical analysis of three-dimensional interface crack

Using the fundamental solution of interface crack and the method of finite-part integral, the problem of three-dimensional interface crack is reduced to solve a set of two-dimensional hypersingular integrodifferential equations with unknown displacement discontinuities of crack surface. Then a systematically theoretical analysis for solving these equations is presented.     

The finite element method for fractional diffusion with spectral fractional Laplacian

This paper focuses on the finite element method for Caputo-type parabolic equation with spectral fractional Laplacian, where the time derivative is in the sense of Caputo with order in (0,1) and the spatial derivative is the spectral fractional Laplacian. The time discretization is based on the Hadamard finite-part integral (or the finite-part integral in the sense of Hadamard), where the piecewise linear interpolation polynomials are used. The spatial fractional Laplacian is lifted to the local spacial derivative by using the Caffarelli–Silvestre extension, where the finite element method is used. Full-discretization scheme is constructed. The convergence and error estimates are obtained. Finally, numerical experiments are presented which support the theoretical results.     

Legendre小波求解超奇异积分

超奇异积分的数值算法一直是近些年来研究的重要课题. 基于超奇异积分的 Hadamard 有限部分积分定义, 本文给出了利用 Legendre 小波计算超奇异积分的方法. 当奇异点位于区间内时, 由于 Legendre 小波具有很好的正交性、显式表达式以及小波函数的可计算性, 将区间内的奇异点变换到区间端点处, 再利用区间端点处 Hadamard 有限部分积分的定义,进而可以计算 p+1(p∈N⁺) 阶超奇异积分. 文中最后给出的算例表明了该方法的可行性和有效性.     

The superconvergence of the composite midpoint rule for the finite-part integral

The composite midpoint rule is probably the simplest one among the Newton-Cotes rules for Riemann integral. However, this rule is divergent in general for Hadamard finite-part integral. In this paper, we turn this rule to a useful one and, apply it to evaluate Hadamard finite-part integral as well as to solve the relevant integral equation. The key point is based on the investigation of its pointwise superconvergence phenomenon, i.e., when the singular point coincides with some a priori known point, the convergence rate of the midpoint rule is higher than what is globally possible. We show that the superconvergence rate of the composite midpoint rule occurs at the midpoint of each subinterval and obtain the corresponding superconvergence error estimate. By applying the midpoint rule to approximate the finite-part integral and by choosing the superconvergence points as the collocation points, we obtain a collocation scheme for solving the finite-part integral equation. More interesting is that the inverse of the coefficient matrix of the resulting linear system has an explicit expression, by which an optimal error estimate is established. Some numerical examples are provided to validate the theoretical analysis.     

Superconvergence of the composite Simpson’s rule for a certain finite-part integral and its applications

The superconvergence phenomenon of the composite Simpson’s rule for the finite-part integral with a third-order singularity is studied. The superconvergence points are located and the superconvergence estimate is obtained. Some applications of the superconvergence result, including the evaluation of the finite-part integrals and the solution of a certain finite-part integral equation, are also discussed and two algorithms are suggested. Numerical experiments are presented to confirm the superconvergence analysis and to show the efficiency of the algorithms.     

The numerical evaluation of Hadamard finite-part integrals

Summary A quadrature rule is described for the numerical evaluation of Hadamard finite-part integrals with a double pole singularity within the range of integration. The rule is based upon the observation that such an integral is the derivative of a Cauchy principal value integral.     

Sigmoidal transformations and the Euler-Maclaurin expansion for evaluating certain Hadamard finite-part integrals

Summary. Starting with some results of Lyness concerning the Euler-Maclaurin expansion of Cauchy principal value integrals over it is shown how, by the use of sigmoidal transformations, good approximations can be found for the Hadamard finite-part integral where The analysis is illustrated by some numerical examples. Received March 13, 1996     

The practical Gauss type rules for Hadamard finite-part integrals using Puiseux expansions

A general framework is constructed for efficiently and stably evaluating the Hadamard finite-part integrals by composite quadrature rules. Firstly, the integrands are assumed to have the Puiseux expansions at the endpoints with arbitrary algebraic and logarithmic singularities. Secondly, the Euler-Maclaurin expansion of a general composite quadrature rule is obtained directly by using the asymptotic expansions of the partial sums of the Hurwitz zeta function and the generalized Stieltjes constant, which shows that the standard numerical integration formula is not convergent for computing the Hadamard finite-part integrals. Thirdly, the standard quadrature formula is recast in two steps. In step one, the singular part of the integrand is integrated analytically and in step two, the regular integral of the remaining part is evaluated using the standard composite quadrature rule. In this stage, a threshold is introduced such that the function evaluations in the vicinity of the singularity are intentionally excluded, where the threshold is determined by analyzing the roundoff errors caused by the singular nature of the integrand. Fourthly, two practical algorithms are designed for evaluating the Hadamard finite-part integrals by applying the Gauss-Legendre and Gauss-Kronrod rules to the proposed framework. Practical error indicator and implementation involved in the Gauss-Legendre rule are addressed. Finally, some typical examples are provided to show that the algorithms can be used to effectively evaluate the Hadamard finite-part integrals over finite or infinite intervals.     

热载荷作用下平面裂纹问题的奇异积分方程与边界无法

从边界积分方程出发,导出了二维裂纹体热传导问题及热弹性问题的积分方程组,继而使用奇异积分方程与边界元相结合的方法,为其建立了相应的数值求解方法。此外,利用奇异积分方程的主部分析法,严格地证明了裂纹尖端温度梯度场的1/√r 奇异性,并且给出了奇性温度梯度场的精确解。最后。对一些典型例子,做了数值计算。     

The numerical evaluation of principal value integrals by finite-part integration

Summary We give a numerical formula for the evaluation of finite-part integrals of the form This method is very convenient for computational purposes since mere scalar products of certain weights and function values have to be calculated. Iff ^(2m-1) (s)=0,m=1,2, ..., [k/2],k>1 the above integral reduces to a generalized principal value integral.