Method of descriptive regularization and quality of approximate solutions

A method of solving Fredholm integral equations of the first kind is described, which is based on the a priori knowledge of the arrangement of extrema and inflection points of the desired solution and permits taking account of the fundamental qualitative regularities inherent in the exact solution of the problem.     

Single edge-cracked orthotropic strip under three-point load

The problem of an orthotropic strip having a crack of unit length normal to one edge and subjected to a bending moment resulting from three-point loading is solved using integral transform method. The mixed boundary conditions lead to dual integral equations which are ultimately reduced to a Fredholm integral equation of second kind. The integral equation thus obtained is solved by the method developed by Fox and Goodwin. Numerical solutions for a fibre-reinforced composite material have been carried out to determine the stress intensity factor of an orthotropic medium. The same has been compared with the isotropic case.     

On crack problem in an elastic ponderable layer

The paper deals with the plane problem of stress distribution in an elastic ponderable layer with a stationary edge crack normal to the boundary plane. The layer is situated and fixed on a rigid foundation. The stresses are caused by action of body forces. By using the method of Fourier transforms the problem is reduced to a system of dual integral equations and next, to a Fredholm integral equation of the second kind. The numerical analysis of the Fredholm equation permitted to determine the stress intensity factor and the crack opening displacement. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mutually constrained partial differential and integral equations for an exterior field problem

A new implementation of the mutually constrained partial differential and integral equation method for the exterior 2-dimensional field problem is described. It is shown, that the method is applicable to exterior problems in an inhomogeneous medium. The inhomogeneity is considered in the finite element procedure and in boundary element method, where an adequate Green's function is applied. The temperature distribution around a three-cable system is then computed as an illustration. The eddy-current losses in the cable sheaths are calculated using the Fredholm integral equation of the second kind.     

Numerical expansion-iterative method for analysis of integral equation models arising in one- and two-dimensional electromagnetic scattering

Most integral equations of the first kind are ill-posed, and obtaining their numerical solution needs often to solve a linear system of algebraic equations of large condition number. So, solving this system may be difficult or impossible. Since many problems in one- and two-dimensional scattering from perfectly conducting bodies can be modeled by Fredholm integral equations of the first kind, this paper presents an effective numerical expansion-iterative method for solving them. This method is based on vector forms of block-pulse functions. By using this approach, solving the first kind integral equation reduces to solve a recurrence relation. The approximate solution is most easily produced iteratively via the recurrence relation. Therefore, computing the numerical solution does not need to directly solve any linear system of algebraic equations and to use any matrix inversion. Also, the method practically transforms solving of the first kind Fredholm integral equation which is inherently ill-posed into solving second kind Fredholm integral equation. Another advantage is low cost of setting up the equations without applying any projection method such as collocation, Galerkin, etc. To show convergence and stability of the method, some computable error bounds are obtained. Test problems are provided to illustrate its accuracy and computational efficiency, and some practical one- and two-dimensional scatterers are analyzed by it.     

Permeance computation using indirect boundary integral equation method

An approach using indirect boundary integral equation method is proposed to determine the permeance between ferromagnetic poles in axisymmetric and three-dimensional magnetic systems. A generalised mathematical model is given for both types of magnetic systems. It consists of Fredholm integral equations of the first kind with respect to fictitious magnetic charge density sought in the form of simple layer potential. The system of boundary integral equations is solved using the method of mechanical quadratures. The approach is implemented in its own computer code. Results are presented for axisymmetric poles of electromagnets (cylinders, cones and frustum cones) and for a three-dimensional clapper-type system. Comparisons with known formulas are made and their accuracies are estimated. The approach presented is useful at the stage of preliminary design of magnetic systems. It is also applicable to computation of capacitances and electrical conductances     

An initial value method for dual integral equations

An initial value method is derived for a set of dual integral equations encountered in solving mixed boundary value problems in mathematical physics with a circular line of separation of boundary conditions. It is shown that the solution itself, not just a transform of the solution, of the dual integral equations satisfies a Fredholm integral equation. The initial value problem is derived from this Fredholm equation.     

Circulant preconditioners for ill-conditioned boundary integral equations from potential equations

In this paper, we consider solving potential equations by the boundary integral equation approach. The equations so derived are Fredholm integral equations of the first kind and are known to be ill-conditioned. Their discretized matrices are dense and have condition numbers growing like O(n) where n is the matrix size. We propose to solve the equations by the preconditioned conjugate gradient method with circulant integral operators as preconditioners. These are convolution operators with periodic kernels and hence can be inverted efficiently by using fast Fourier transforms. We prove that the preconditioned systems are well conditioned, and hence the convergence rate of the method is linear. Numerical results for two types of regions are given to illustrate the fast convergence. © 1998 John Wiley & Sons, Ltd.     

Computational methods for solving static field and eddy current problems via Fredholm integral equations

Two-dimensional static field problems can be solved by a method based on Fredholm integral equations (equations of the second kind). This has numerical advantages over the mote commonly used integral equation of the first kind. The method is applicable to both magnetostatic and electrostatic problems formulated in terms of either vector or scalar potentials. It has been extended to the solution of eddy current problems with sinusoidal driving functions. The application of the classical Fredholm equation has been extended to problems containing boundary conditions: 1) potential value, 2) normal derivative value, and 3) an interface condition, all in the same problem. The solutions to the Fredholm equations are single or double (dipole) layers of sources on the problem boundaries and interfaces. This method has been developed into computer codes which use piecewise quadratic approximations to the solutions to the integral equations. Exact integrations are used to replace the integral equations by a matrix equation. The solution to this matrix equation can then be used to directly calculate the field anywhere.     

Transient analysis of a piezoelectric strip with a permeable crack under anti-plane impact loads

The problem of a through permeable crack situated in the mid-plane of a piezoelectric strip is considered under anti-plane impact loads for two cases. The first is that the strip boundaries are free of stresses and of electric displacements, and the second is that the strip boundaries are clamped rigid electrodes. The method adopted is to reduce the mixed initial-boundary value problem, by using integral transform techniques, to dual integral equations, which are further transformed into a Fredholm integral equation of the second kind by introducing an auxiliary function. The dynamic stress intensity factor and energy release rate in the Laplace transform domain are obtained in explicit form in terms of the auxiliary function. Some numerical results for the dynamic stress intensity factor are presented graphically in the physical space by using numerical techniques for solving the resulting Fredholm integral equation and inverting Laplace transform.     

第二类Fredholm积分方程的小波快速算法

介绍了一种求解含有对数核的第二类Fredholm积分方程的有效的方法，该方法先用Nystrom法将积分方程离散，然后用小波矩阵变换方法稀疏系数矩阵，对系数矩阵预处理后再对线性方程组迭代求解。数值结果证明了该方法的有效性。     

Edge crack in an elastic layer resting on Winkler foundation

The paper deals with the stress analysis near a crack tip in an elastic layer resting on Winkler foundation. The edge crack is assumed to be normal to the lower boundary plane. The upper surface of the layer is loaded by given forces normal to the boundary. The considered problem is solved by using the method of Fourier transforms and dual integral equations, which are reduced to a Fredholm integral equation of the second kind. The stress intensity factor is given in the term of solution of the Fredholm integral equation and some numerical results are presented.     

An Analytical Method for Computing the One-Dimensional Backward Wave Problem

The present paper reveals a new computational method for the illposed backward wave problem. The Fourier series is used to formulate a first-kind Fredholm integral equation for the unknown initial data of velocity. Then, we consider a direct regularization to obtain a second-kind Fredholm integral equation. The termwise separable property of kernel function allows us to obtain an analytical solution of regularization type. The sufficient condition of the data for the existence and uniqueness of solution is derived. The error estimate of the regularization solution is provided. Some numerical results illustrate the performance of the new method.     

A note on the Gauss-Jacobi quadrature formulae for singular integral equations of the second kind

A fast and efficient numerical method based on the Gauss-Jacobi quadrature is described that is suitable for solving Fredholm singular integral equations of the second kind that are frequently encountered in fracture and contact mechanics. Here we concentrate on the case when the unknown function is singular at both ends of the interval. Quadrature formulae involve fixed nodal points and provide exact results for polynomials of degree 2n − 1, where n is the number of nodes. Finally, an application of the method to a plane problem involving complete contact is presented.     

Problem of an infinite solid containing a flat annular crack under torsion

An analytic method is developed to find the axisymmetric stress distribution in an infinite elastic solid containing a flat annular crack under axial torsion. By use of Hankel transforms, the solution to the problem is reduced to triple-integral equations involving Bessel functions of order 1. Modifying the method discussed by Cooke[Quart. J. Mech. Appl. Math. 16, 193–203 (1963).], the solution of the triple-integral equations is reduced to a pair of Fredholm integral equations of the second kind. Finally, the approximate expressions for the stress intensity factors are obtained by finding the iterative solution of the pair of Fredholm integral equations.     

Multiple crack problems of antiplane elasticity in an infinite body

Twe elementary solutions are presented for case of a pair of normal or tangential concentrated unit forces acting at a point of both edges of a single crack in an infinite plane isotropic elastic medium. Using these two elementary solutions and the principle of superposition, we found that the multiple crack problems can be easily converted into a system of Fredholm integral equations. Finally, the system obtained is solved numerically and the values of the stress intensity factors at the crack tips can be easily calculated. Two numerical examples are given in this paper. A system of Fredholm integral equations is complex form is also presented. We found that the system of Fredholm integral equations can be easily reduced from the system of singular integral equations given by Panasyuk[1]     

Thermal stresses near a penny-shaped crack in an elastic sphere embedded in an infinite elastic space

This paper gives an analysis of the distribution of thermal stresses in a sphere which is bonded to an infinite elastic medium. The thermal and the elastic properties of the sphere and the elastic infinite medium are assumed to be different. The penny-shaped crack lies on the diametral plane of the sphere and the centre of the crack is the centre of the sphere. By making a suitable representation of the temperature function, the heat conduction problem is reduced to the solution of a Fredholm integral equation of the second kind. Using suitable solution of the thermoelastic displacement differential equation, the problem is then reduced to the solution of a Fredholm integral equation, in which the solution of the earlier integral equation arising from heat conduction problem occurs as a known function. Numerical solutions of these two Fredholm integral equations are obtained. These solutions are used to evaluate numerical values for the stress intensity factors. These values are displayed graphically.     

Penny-shaped interface crack between dissimilar nonhomogeneous elastic layers under axially symmetric torsion

Summary The problem of axially symmetric torsion for dissimilar nonhomogeneous bonded elastic layers containing a penny-shaped interface crack is considered. The mixed boundary value problem is reduced to solving a Fredholm integral equation of the second kind. The Fredholm integral equation is solved numerically by reducing it to a system of simultaneous algebraic equations. Numerical results for the stress intensity factor are presented in the form of graphs.     

Direct second kind boundary integral formulation for Stokes flow problems

A direct boundary element method is formulated for the Stokes flow problem based on an integral equation representation for the components of traction. For problems in which the components of velocity are prescribed on the boundary of the domain, this new formulation results in a hypersingular Fredholm integral equation of the second kind. A method of regularization to evaluate the hypersingular integral is discussed. For certain problems involving flows about particles, the integral equation representation for the tractions is not unique because of the existence of rigid body eigenmodes. A method to constrain out these rigid body modes is also discussed. Several example problems are considered in which this new formulation is compared to more traditional boundary element formulations.     

Synthesis of magnetic fields

The paper deals with some problems of magnetic fields synthesis, depending on determination of the current density distribution, which generates the required magnetic field in the investigated region. Such problems can be reduced to the linear, or nonlinear Fredholm integral equations of the first kind, or to the set of these equations. Fredholm integral equation of the first kind belongs to the class of the ill-posed problems, and for its solving the method of regularisation has been used. In the paper there are given some useful results of synthesis of magnetic fields in few practical configurations.