粘弹性材料界面裂纹的瞬态响应和散射研究

本文研究粘弹性材料界面裂纹对冲击载荷的瞬态响应和对广义平面波的稳态散射。相对于已有广泛研究的弹性材料裂纹瞬态响应和稳态散射问题，本文的研究有三个突出特点：1）粘弹性材料；2）界面裂纹；3）广义平面波入射。粘弹性材料界面裂纹对冲击载荷的瞬态响应和对广义平面波的散射尚无开展研究，本文在弹性材料相应问题的研究基础上，首先开展了这一问题的研究。对于冲击载荷下粘弹性界面裂纹的瞬态响应问题，利用Laplace积分变换方法，将粘弹性材料卷积型本构方程转化为Laplace变换域内的代数型本构方程，从而可以在Laplace变换域内象处理弹性材料的冲击响应一样，将相应的混合边值问题归结为关于裂纹张开位移COD的对偶积分方程，并进一步引入裂纹位错密度函数CDD （Crack Dislocation Density），将对偶积分方程化成关于CDD的奇异积分方程（SIE）。用数值方法求解奇异积分方程得到变换域内的动应力强度因子数值解，最后利用Laplace积分逆变换数值方法得到时间域内的动应力强度因子的时间响应。理论分析考虑了两种裂纹模型，即Griffith界面裂纹和柱面圆弧型界面裂纹。考虑的载荷包括反平面冲击载荷和平面冲击载荷。对于平面冲击载荷，通过对裂尖应力场的奇性分析，首次发现粘弹性界面裂纹裂尖动应力场奇性指数不是常数0.5，而是与震荡指数一样依赖材料参数。针对反平面冲击载荷给出了一个算例，计算了裂尖动应力强度因子的时间响应，并与弹性材料的结果作了比较，发现粘弹性效应的影响不仅使过冲击峰值降低，而且使峰值点后移。粘性效应较大时，过冲击现象甚至不出现。关于粘弹性界面裂纹对广东省义平面波的散射问题，首先研究广义平面波在无裂纹存在的理想界面的反射和透射，再研究由于界面裂纹的存在而产生的附加散射场。利用粘弹性材料的复模量理论，可将粘弹性材料的卷积型相构方程化成频率域内的代数型本构方程。类似弹性平面波的处理，在频率域内将问题最终归结为关于裂纹位错密度CDD的奇异积分方程。数值方法求解奇异积分方程即可得到频率域内的散射场，并进而得到裂尖动应力强度因子和远场位移型函数和散射截面。理论分析考虑了两种裂纹模型：Griffith界面裂纹和柱面圆弧型界面裂纹。研究的入射波有广义的SH波和P波。对于广义平面P波入射的情况，通过对裂尖应力场的奇性分析，同样发现粘弹性界面裂纹裂尖动应力场奇性指数不地常数0.5，而是与震荡指数一样依赖于材料参数。对柱面裂纹散射远场的渐近分析，发现远场位移和应力除含有几何衰减因子外，都含有一个材料衰减速因子。散射截面由于材料衰减因子的存在也成为依赖散射半径的量。为了使散射截面仍有意义，文中提出一种修正办法。对Griffith界面裂纹，给出了一个广义平面SH波入射的算例；对柱面界面裂纹，给出了一个广义平面P波入射的算例。计算了不同入射角和入射频率下裂纹的张开位移和动就应力强度因子，并分析了其依赖关系。求解奇异积分方程的数值方法和Laplace积分逆变换数值方法是本文的基本数值方法。本文对这两种方法作了大量的调研和系统的研究。在对比分析的基础上，对现有的各种方法从原理，适用范围，计数效率，优势及特点进行了归纳总结。并尝试了奇异积分方程的最新数值方法--分片连续函数法，证实了其适用性和方便性．

An alternating boundary integral based method for a Cauchy problem for the Laplace equation in semi-infinite regions

We consider a Cauchy problem for the Laplace equation in a two-dimensional semi-infinite region with a bounded inclusion, i.e. the region is the intersection between a half-plane and the exterior of a bounded closed curve contained in the half-plane. The Cauchy data are given on the unbounded part of the boundary of the region and the aim is to construct the solution on the boundary of the inclusion. In 1989, Kozlov and Maz'ya [10] proposed an alternating iterative method for solving Cauchy problems for general strongly elliptic and formally self-adjoint systems in bounded domains. We extend their approach to our setting and in each iteration step mixed boundary value problems for the Laplace equation in the semi-infinite region are solved. Well-posedness of these mixed problems are investigated and convergence of the alternating procedure is examined. For the numerical implementation an efficient boundary integral equation method is proposed, based on the indirect variant of the boundary integral equation approach. The mixed problems are reduced to integral equations over the (bounded) boundary of the inclusion. Numerical examples are included showing the feasibility of the proposed method.

Fast reconstruction of harmonic functions from Cauchy data using the Dirichlet-to-Neumann map and integral equations

We propose and investigate a method for the stable determination of a harmonic function from knowledge of its value and its normal derivative on a part of the boundary of the (bounded) solution domain (Cauchy problem). We reformulate the Cauchy problem as an operator equation on the boundary using the Dirichlet-to-Neumann map. To discretize the obtained operator, we modify and employ a method denoted as Classic II given in [J. Helsing, Faster convergence and higher accuracy for the Dirichlet–Neumann map, J. Comput. Phys. 228 (2009), pp. 2578–2576, Section 3], which is based on Fredholm integral equations and Nyström discretization schemes. Then, for stability reasons, to solve the discretized integral equation we use the method of smoothing projection introduced in [J. Helsing and B.T. Johansson, Fast reconstruction of harmonic functions from Cauchy data using integral equation techniques, Inverse Probl. Sci. Eng. 18 (2010), pp. 381–399, Section 7], which makes it possible to solve the discretized operator equation in a stable way with minor computational cost and high accuracy. With this approach, for sufficiently smooth Cauchy data, the normal derivative can also be accurately computed on the part of the boundary where no data is initially given.

Cauchy Problem for Differential Equation with Caputo Derivative

The paper is devoted to the study of the Cauchy problem for a nonlinear differential equation of complex order with the Caputo fractional derivative. The equivalence of this problem and a nonlinear Volterra integral equation in the space of continuously differentiable functions is established. On the basis of this result, the existence and uniqueness of the solution of the considered Cauchy problem is proved. The approximate-iterative method by Dzjadyk is used to obtain the approximate solution of this problem. Two numerical examples are given.

Boundary value problems for analytic functions in the class of Cauchy-type integrals with density in

We study the Riemann boundary value problem , for analytic functions in the class of analytic functions represented by the Cauchy-type integrals with density in the spaces with variable exponent. We consider both the case when the coefficient is piecewise continuous and it may be of a more general nature, admitting its oscillation. The explicit formulas for solutions in the variable exponent setting are given. The related singular integral equations in the same setting are also investigated. As an application there is derived some extension of the Szegö-Helson theorem to the case of variable exponents.

A Cauchy Integral Related to a Robot-safety Device System

We introduce a robot-safety device system attended by two different repairmen. The twin system is characterized by the natural feature of cold standby and by an admissible “risky” state. In order to analyse the random behaviour of the entire system (robot, safety device, repair facility) we employ a stochastic process endowed with probability measures satisfying general Hokstad-type differential equations. The solution procedure is based on the theory of sectionally holomorphic functions, characterized by a Cauchy-type integral defined as a Cauchy principal value in double sense. An application of the Sokhotski-Plemelj formulae determines the long-run availability of the robot-safety device. Finally, we consider the particular but important case of deterministic repair.

On a Class of Fractional Type Integral Equations in Variable Exponent Spaces

2000 Mathematics Subject Classification: 45A05, 45B05, 45E05,45P05, 46E30

On a Hypersingular Equation of a Problem, for a Crack in Elastic Media

Mathematics Subject Classification 2010: 45DB05, 45E05, 78A45.

On a 3D-Hypersingular Equation of a Problem for a Crack

MSC 2010: 45DB05, 45E05, 78A45

Numerical solution of the Dirichlet initial boundary value problem for the heat equation in exterior 3-dimensional domains using integral equations

A numerical method for the Dirichlet initial boundary value problem for the heat equation in the exterior and unbounded region of a smooth closed simply connected 3-dimensional domain is proposed and investigated. This method is based on a combination of a Laguerre transformation with respect to the time variable and an integral equation approach in the spatial variables. Using the Laguerre transformation in time reduces the parabolic problem to a sequence of stationary elliptic problems which are solved by a boundary layer approach giving a sequence of boundary integral equations of the first kind to solve. Under the assumption that the boundary surface of the solution domain has a one-to-one mapping onto the unit sphere, these integral equations are transformed and rewritten over this sphere. The numerical discretisation and solution are obtained by a discrete projection method involving spherical harmonic functions. Numerical results are included.

Numerical solution of stress intensity factors of multiple cracks in great number with eigen COD boundary integral equations

Based on the eigen crack opening displacement (COD) boundary integral equations, a newly developed computational approach is proposed for the analysis of multiple crack problems. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix. The interactions among cracks are dealt with by two parts according to the distances of cracks to the current crack. The strong effects of cracks in adjacent group are treated with the aid of the local Eshelby matrix derived from the traction BIEs in discrete form. While the relatively week effects of cracks in far-field group are treated in the iteration procedures. Numerical examples are provided for the stress intensity factors of multiple cracks, up to several thousands in number, with the proposed approach. By comparing with the analytical solutions in the literature as well as solutions of the dual boundary integral equations, the effectiveness and the efficiencies of the proposed approach are verified.

Efficient solution of multiple cracks in great number using eigen COD boundary integral equations with iteration procedure

A newly developed computational approach is proposed in the paper for the analysis of multiple crack problems based on the eigen crack opening displacement (COD) boundary integral equations. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix to determine all the unknown CODs step by step. To deal with the interactions among cracks for multiple crack problems, all cracks in the problem are divided into two groups, namely the adjacent group and the far-field group, according to the distance to the current crack in consideration. The adjacent group contains cracks with relatively small distances but strong effects to the current crack, while the others, the cracks of far-field group are composed of those with relatively large distances. Correspondingly, the eigen COD of the current crack is computed in two parts. The first part is computed by using the fictitious tractions of adjacent cracks via the local Eshelby matrix derived from the traction boundary integral equations in discretized form, while the second part is computed by using those of far-field cracks so that the high computational efficiency can be achieved in the proposed approach. The numerical results of the proposed approach are compared not only with those using the dual boundary integral equations (D-BIE) and the BIE with numerical Green's functions (NGF) but also with those of the analytical solutions in literature. The effectiveness and the efficiency of the proposed approach is verified. Numerical examples are provided for the stress intensity factors of cracks, up to several thousands in number, in both the finite and infinite plates.

Eigenstrain boundary integral equations with local Eshelby matrix for stress analysis of ellipsoidal particles

Aiming at the large scale numerical simulation of particle reinforced materials, the concept of local Eshelby matrix has been introduced into the computational model of the eigenstrain boundary integral equation (BIE) to solve the problem of interactions among particles. The local Eshelby matrix can be considered as an extension of the concepts of Eshelby tensor and the equivalent inclusion in numerical form. Taking the subdomain boundary element method as the control, three-dimensional stress analyses are carried out for some ellipsoidal particles in full space with the proposed computational model. Through the numerical examples, it is verified not only the correctness and feasibility but also the high efficiency of the present model with the corresponding solution procedure, showing the potential of solving the problem of large scale numerical simulation of particle reinforced materials.

Approximate solutions of Fredholm integral equations of the second kind

This note is concerned with the problem of determining approximate solutions of Fredholm integral equations of the second kind. Approximating the solution of a given integral equation by means of a polynomial, an over-determined system of linear algebraic equations is obtained involving the unknown coefficients, which is finally solved by using the least-squares method. Several examples are examined in detail. (c) 2009 Elsevier Inc. All rights reserved.