The computation of the J integral with the traction boundary integral equation

The solutions of the displacement boundary integral equation (BIE) are not uniquely determined in certain types of boundary conditions. Traction boundary integral equations that have unique solutions in these traction and mixed boundary cases are established. For two‐dimensional linear elasticity problems, the divergence‐free property of the traction boundary integral equation is established. By applying Stokes' theorem, unknown tractions or displacements can be reduced to computation of traction integral potential functions at the boundary points. The same is true of the J integral: it is divergence‐free and the evaluation of the J integral can be inverted into the computation of the J integral potential functions at the boundary points of the cracked body. The J integral can be expressed as the linear combination of the tractions and displacements from the traction BIE on the boundary of the cracked body. Numerical integrals are not needed at all. Selected examples are presented to demonstrate the validity of the traction boundary integral and J integral. Copyright © 2004 John Wiley & Sons, Ltd.     

New analytical approximate solution of the generalised temperature integral for kinetic reactions

Simple closed-form expression of the generalised temperature integral in the basic equation to describe kinetic reactions for solid materials in linear heating process is always suitable for use in determining parameters. Many developed solutions only can give high accuracies on the general conditions. A new analytical approximate solution was deduced in this work. The deviations of this solution from the true value are fully analysed. This solution takes advantage in broader application conditions than other known solutions. The application of the new analytical approximate solution in austenite kinetic reaction in this work reinforces that austenite reaction rate is feasible to be given a priori.     

Analytic preconditioners for the electric field integral equation

Since the advent of the fast multipole method, large‐scale electromagnetic scattering problems based on the electric field integral equation (EFIE) formulation are generally solved by a Krylov iterative solver. A well‐known fact is that the dense complex non‐hermitian linear system associated to the EFIE becomes ill‐conditioned especially in the high‐frequency regime. As a consequence, this slows down the convergence rate of Krylov subspace iterative solvers. In this work, a new analytic preconditioner based on the combination of a finite element method with a local absorbing boundary condition is proposed to improve the convergence of the iterative solver for an open boundary. Some numerical tests precise the behaviour of the new preconditioner. Moreover, comparisons are performed with the analytic preconditioner based on the Calderòn's relations for integral equations for several kinds of scatterers. Copyright © 2004 John Wiley & Sons, Ltd.     

Galerkin boundary integral analysis for the axisymmetric Laplace equation

The boundary integral equation for the axisymmetric Laplace equation is solved by employing modified Galerkin weight functions. The alternative weights smooth out the singularity of the Green's function at the symmetry axis, and restore symmetry to the formulation. As a consequence, special treatment of the axis equations is avoided, and a symmetric‐Galerkin formulation would be possible. For the singular integration, the integrals containing a logarithmic singularity are converted to a non‐singular form and evaluated partially analytically and partially numerically. The modified weight functions, together with a boundary limit definition, also result in a simple algorithm for the post‐processing of the surface gradient. Published in 2005 by John Wiley & Sons, Ltd.     

非等温动力学积分-微分法研究聚苯乙烯热解反应机理及动力学参数

首次采用积分-微分相结合的方法研究了聚苯乙烯非等温热解动力学,通过比较同一机理函数的微分和积分法结果,推测了聚苯乙烯热解所遵循的机理及动力学参数。实验结果显示,在氮气中聚苯乙烯有一个失重阶段,温度范围282.9℃~398.2℃,在此范围内约有96%聚苯乙烯发生了分解,其最大失重速率发生在372.4℃。积分-微分法拟合结果表明,聚苯乙烯在氮气中的热解符合一级化学反应,其动力学模型方程式为:dα/dt=A exp(-E/RT)(1-α)。平均表观活化能E为199.2 kJ/m o l,lnA为36.56。     

A new method for the numerical solution of integral equation approximations

A new numerical technique for solving the Ornstein-Zernike equation is described. It is particularly useful in solving the Ornstein-Zernike equation for approximations and pair potentials (such as the Percus-Yevick and mean spherical approximations for finite ranged potentials) which imply a finiteranged direct correlation function since for such approximations the numerical technique is essentially exact. The only approximation involved in such cases is the discretization of direct and total correlation functions over the finite range on which the direct correlation function is nonzero. Thus, the new method avoids truncation of the total correlation function and should permit the critical point and spinodal curve to be mapped out with greater accuracy than is permitted by existing methods. Preliminary explorations on the stability and accuracy of the method are described.Paper presented at the Tenth Symposium on Thermophysical Properties, June 20–23, 1988, Gaithersburg, Maryland, U.S.A.     

Complex hypersingular integral equation for the piece-wise homogeneous half-plane with cracks

New complex hypersingular integral equation (CHSIE) is derived for the half-plane containing the inclusions (which can have the different elastic properties), holes, notches and cracks of the arbitrary shape. This equation is obtained by superposition of the equations for each homogeneous region in a half-plane. The last equations follow from the use of complex analogs of Somigliana's displacement and stress identities (SDI and SSI) and Melan's fundamental solution (FS) written in a complex form. The universal numerical algorithm suggested before for the analogous problem for a piece-wise homogeneous plane is extended on case of a half plane. The unknown functions are approximated by complex Lagrange polynomials of the arbitrary degree. The asymptotics for the displacement discontinuities (DD) at the crack tips are taken into account. Only two types of the boundary elements (straight segments and circular arcs) are used to approximate the boundaries. All the integrals involved in CHSIE are evaluated in a closed form. A wide range of elasticity problems for a half-plane with cracks, openings and inclusions are solved numerically.     

Boundary integral equation formulation for the generalized micro-polar thermoviscoelasticity

A formulation of the boundary integral equation method for generalized linear micro-polar thermoviscoelasticity is given. Fundamental solutions, in Laplace transform domain, of the corresponding differential equations are obtained. The initial, mixed boundary value problem is considered as an example illustrating the BIE formulation. The results are applicable to the generalized thermoelasticity theories: Lord-Shulman with one relaxation time, Green-Lindsay with two relaxation times, Green-Naghdi theories, and Chandrasekharaiah and Tzou with dual-phase lag, as well as to the dynamic coupled theory. The cases of generalized linear micro-polar thermoviscoelasticity of Kelvin-Voigt model, generalized linear thermoviscoelasticity and generalized thermoelasticity can be obtained from the given results.     

Time marching integral equation method for unsteady transonic flows around airfoils

A time marching integral equation method has been proposed here which does not have the limitation of the time linearized integral equation method in that the latter method can not satisfactorily simulate the shock wave motions. Firstly, a model problem–one dimensional initial and boundary value wave problem is treated to clarify the basic idea of the new method. Then the method is implemented for two dimensional unsteady transonic flow problems. The introduction of the concept of a quasi-velocity-potential simplifies the time marching integral equations and the treatment of trailing vortex sheet condition. The numerical calculations show that the method is reasonable and reliable.     

Boundary integral equation methods for two-dimensional models of crack-field interactions

An introduction to the application of surface integral equation methods to the calculation of eddy current-flaw interactions is presented. Two two-dimensional problems are presented which are solved by the boundary integral equation method. Application of collocation methods reduces the problems to systems of linear algebraic equations. The first problem is that of a closed surface crack in a flat slab with an AC magnetic field parallel to the plane of the crack. The second is that of av-groove crack in the AC field of a pair of parallel wires placed parallel to the vertex of the crack. In both cases, maps of the current densities at the surface are displayed, as well as the impedance changes due to the cracks.     

Material constant sensitivity boundary integral equation for anisotropic solids

In this paper, the material constant sensitivity boundary integral equation is presented, and its numerical solution proposed, based on boundary element techniques. The formulation deals with plane problems with general rectilinear anisotropy. Expressions for the computation of sensitivities for displacements, tractions, strains and stresses are derived, both for boundary and interior points. The sensitivities can be computed with respect to the bulk material properties or to the properties of part of the domain (inclusions, coatings, etc.). To assess the accuracy of the proposed approach, the computed results are compared to analytical ones derived from exact solutions obtained by complex potential theory, when possible, or finite difference derivatives otherwise. Copyright © 2005 John Wiley & Sons, Ltd.     

On the Cauchy principal value of the surface integral in the boundary integral equation of 3D elasticity

A simple demonstration of the existence of the Cauchy principal value (CPV) of the strongly singular surface integral in the Somigliana Identity at a non-smooth boundary point is presented. First a regularization of the strongly singular integral by analytical integration of the singular term in the radial direction in pre-image planes of smooth surface patches is carried out. Then it is shown that the sum of the angular integrals of the characteristic of the tractions of the Kelvin fundamental solution is zero, a formula for the transformation of angles between the tangent plane of a suface patch and the pre-image plane at smooth mapping of the surface patch being derived for this purpose.     

Numerical solution of the variation boundary integral equation for inverse problems

In this paper a procedure to solve the identification inverse problems for two‐dimensional potential fields is presented. The procedure relies on a boundary integral equation (BIE) for the variations of the potential, flux, and geometry. This equation is a linearization of the regular BIE for small changes in the geometry. The aim in the identification inverse problems is to find an unknown part of the boundary of the domain, usually an internal flaw, using experimental measurements as additional information. In this paper this problem is solved without resorting to a minimization of a functional, but by an iterative algorithm which alternately solves the regular BIE and the variation BIE. The variation of the geometry of the flaw is modelled by a virtual strainfield, which allows for greater flexibility in the shape of the assumed flaw. Several numerical examples demonstrate the effectiveness and reliability of the proposed approach. Copyright © 2000 John Wiley & Sons, Ltd.     

Singularity extraction technique for integral equation methods with higher order basis functions on plane triangles and tetrahedra

A numerical solution of integral equations typically requires calculation of integrals with singular kernels. The integration of singular terms can be considered either by purely numerical techniques, e.g. Duffy's method, polar co‐ordinate transformation, or by singularity extraction. In the latter method the extracted singular integral is calculated in closed form and the remaining integral is calculated numerically. This method has been well established for linear and constant shape functions. In this paper we extend the method for polynomial shape functions of arbitrary order. We present recursive formulas by which we can extract any number of terms from the singular kernel defined by the fundamental solution of the Helmholtz equation, or its gradient, and integrate the extracted terms times a polynomial shape function in closed form over plane triangles or tetrahedra. The presented formulas generalize the singularity extraction technique for surface and volume integral equation methods with high‐order basis functions. Numerical experiments show that the developed method leads to a more accurate and robust integration scheme, and in many cases also a faster method than, for example, Duffy's transformation. Copyright © 2003 John Wiley & Sons, Ltd.     

Null-field integral equation approach for free vibration analysis of circular plates with multiple circular holes

In this paper, a semi-analytical approach for the eigenproblem of circular plates with multiple circular holes is presented. Natural frequencies and modes are determined by employing the null-field integral formulation in conjunction with degenerate kernels, tensor rotation and Fourier series. In the proposed approach, all kernel functions are expanded into degenerate (separable) forms and all boundary densities are represented by using Fourier series. By uniformly collocating points on the real boundary and taking finite terms of Fourier series, a linear algebraic system can be constructed. The direct searching approach is adopted to determine the natural frequency through the singular value decomposition (SVD). After determining the unknown Fourier coefficients, the corresponding mode shape is obtained by using the boundary integral equations for domain points. The result of the annular plate, as a special case, is compared with the analytical solution to verify the validity of the present method. For the cases of circular plates with an eccentric hole or multiple circular holes, eigensolutions obtained by the present method are compared well with those of the existing approximate analytical method or finite element method (ABAQUS). Besides, the effect of eccentricity of the hole on the natural frequency and mode is also considered. Moreover, the inherent problem of spurious eigenvalue using the integral formulation is investigated and the SVD updating technique is adopted to suppress the occurrence of spurious eigenvalues. Excellent accuracy, fast rate of convergence and high computational efficiency are the main features of the present method thanks to the semi-analytical procedure.     

A new boundary integral equation method for cracked 2-D anisotropic bodies

By using integration by parts to the traditional boundary integral formulation, a traction boundary integral equation for cracked 2-D anisotropic bodies is derived. The new traction integral equation involves only singularity of order 1/r and no hypersingular term appears. The dislocation densities on the crack surface are introduced and the relations between stress intensity factors and dislocation densities near the crack tip are induced to calculate the stress intensity factors. The boundary element method based on the new equation is established and the singular interpolation functions are introduced to model the singularity of the dislocation density (in the order of ) for crack tip elements. The proposed method can be directly used for the 2-D anisotropic body containing cracks of arbitrary geometric shapes. Several numerical examples demonstrate the validity and accuracy of BEM based on the new boundary integral equation.     

Boundary integral equation method for two dissimilar elastic inclusions in an infinite plate

This paper provides a numerical solution for an infinite plate containing two dissimilar elastic inclusions, which is based on complex variable boundary integral equation (CVBIE). The original problem is decomposed into two problems. One is an interior boundary value problem (BVP) for two elastic inclusions, while other is an exterior BVP for the matrix with notches. After performing discretization for the coupled boundary integral equations (BIEs), a system of algebraic equations is formulated. The inverse matrix technique is suggested to solve the relevant algebraic equations, which can avoid using the assembling of some matrices. Several numerical examples are carried out to prove the efficiency of suggested method and the hoop stress along the interface boundary is evaluated.     

Using Lagrange principle for solving two-dimensional integral equation with a positive kernel

This article is devoted to a Lagrange principle application to an inverse problem of a two-dimensional integral equation of the first kind with a positive kernel. To tackle the ill-posedness of this problem, a new numerical method is developed. The optimal and regularization properties of this method are proved. Moreover, a pseudo-optimal error of the proposed method is considered. The efficiency and applicability of this method are demonstrated in a numerical example of an image deblurring problem with noisy data.     

Local integral equation method for viscoelastic Reissner–Mindlin plates

A meshless local Petrov-Galerkin (MLPG) method is applied to solve static and dynamic bending problems of linear viscoelastic plates described by the Reissner–Mindlin theory. To this end, the correspondence principle is applied. A weak formulation for the set of governing equations in the Reissner–Mindlin theory with a unit test function is transformed into local integral equations on local subdomains in the mean surface of the plate. Nodal points are randomly spread on the mean surface of the plate and each node is surrounded by a circular subdomain to which local integral equations are applied. A meshless approximation based on the moving least-squares (MLS) method is employed in the numerical implementation.     

Numerical solutions of hypersingular integral equation for curved cracks in circular regions

In this paper, numerical solutions of a hypersingular integral equation for curved cracks in circular regions are presented. The boundary of the circular regions is assumed to be traction free or fixed. The suggested complex potential is composed of two parts, the principle part and the complementary part. The principle part can model the property of a curved crack in an infinite plate. For the case of the traction free boundary, the complementary part can compensate the traction on the circular boundary caused by the principle part. Physically, the proposed idea is similar to the image method in electrostatics. By using the crack opening displacement (COD) as the unknown function and traction as right hand term in the equation, a hypersingular integral equation for the curved crack problems in the circular regions is obtained. The equation is solved by using the curve length coordinate method. In order to prove that the suggested method can be used to solve more complicated cases of the curved cracks, several numerical examples are given.