Stress intensity factor solutions for two-dimensional elastostatic problems by the hypersingular boundary integral equation

Identities for the fundamental solution of thin plate bending problems and the nonuniqueness of the hypersingular BIE solution for multi-connected domains

Boundary Element Analysis of Cracks in Shear Deformable Plates and Shells. Topics in Engineering, Vol 43

This paper analyzed the interaction between two parallel vertical cracks in a transversely isotropic half space by using the boundary element method (DBEM). The crack surface is perpendicular to the free surface and isotropic plane subjected to normal and tangential uniform distributed load respectively. The stress intensity factor (SIF) values of the crack are calculated from the crack opening displacements. And the interaction effect factor of SIF values is employed to quantitatively describe the interaction between two cracks. The variations of interaction effect factor are investigated with the distances between cracks, the side ratios and free surface. Results show that the existence of the free surface exerts more obvious influence on the SIF values of the crack close to free surface in half space than those in infinite medium. But the free surface has almost no influence on the interaction between two cracks.

Adaptive Gaussian process regression as an alternative to FEM for prediction of stress intensity factor to assess fatigue degradation in offshore pipeline

The spurious eigenvalues of an annular domain have been verified for the singular and hypersingular boundary‐element methods (BEMs) and circumvented by using the Burton‐Miller approach. Do they also occur in other formulations: continuous formulations such as the singular and hypersingular boundary integral equations (BIEs), the null‐field BIEs and the fictitious BIEs, or such discrete formulations as the null‐field BEMs and the fictitious BEMs? For the ten formulations of the multiply connected problem the study of otherwise the same issues is continued in the present paper. By using the degenerate kernels and the Fourier series, it is demonstrated analytically for the six continuous formulations of BIEs that spurious eigensolutions depend on the geometry of the inner boundary but not on that of the outer boundary. This conclusion can be extended to the six discrete formulations of BEMs. To filter out the spurious eigenvalues, the CHIEF (combined Helmholtz integral equation formulation) method is used here instead of the Burton‐Miller approach. The optimum number and appropriate positions of the CHIEF points are also addressed. It is then shown that, in the null‐field and fictitious BEMs, the spurious and true eigenvalues can be detected and distinguished by using the singular‐value‐decomposition‐updating techniques in conjunction with the Fredholm alternative theorem. Illustrative examples show the validity of the proposed methodologies.

Thermoelastic fracture mechanics with regularized hypersingular boundary integral equations

Analysis of cracked magnetoelectroelastic composites under time-harmonic loading

Abstract: In this paper, a half space adaptive fast multipole boundary facemethod (FMBFM) is presented for solving the three-dimensional half space ex-terior acoustic problems. In the presented method, the Burton-Miller equationbased on the conventional boundary integral equation (CBIE) and its hyper-singularboundary integral equation (HBIE) is used to deal with the fictitious eigenfrequen-cies problem. The half space Green’s function is employed, thus the tree structurein the fast multipole method can be used only for the real domain. The higher orderelements and an adaptive tree structure are used to improve the efficiency of theFMBFM. This half space adaptive FMBFM for half space acoustic problems is anextension of the adaptive FMBFM for full space acoustic problems developed bythe authors. Numerical examples for half space acoustic problems in this paperdemonstrate the efficiency and validity of this method.Keywords: fast multipole boundary face method, Burton-Miller equation, acous-tic problems, half space, modified hyper-singular boundary integral equation1 IntroductionAcoustic problems in half space are a major subject in engineering. Many researchworks have been done on this aspect (Seybert et al., 1988; Secgin et al., 2010; Caoet al., 2010). The half space Green’s function (Seybert et al., 1988) can be appliedby adding the source solution at the image point to the original full-space Green’sfunction for dealing with half space problems. Sound source localization and iden-tification of a refrigerator were accomplished by using the half space boundaryelement method (BEM) (Secgin et al., 2010). Using the half space Green’s func-tion, the discretization of the rigid infinite plane is removed and only the boundaries

Linear elastic fracture mechanics (LEFM) concept is generally used to study the influence of crack on the performance of structures. In order to study the LEFM concept on damaged structure, the usage of finite element analysis software is implemented to do the simulation of the structure. Mesh generation is one of the most crucial procedures in finite element method. For the structure that crack or damaged, it is very important to determine the accurate local meshing size at the crack tip of the crack itself in order to get the accurate value of stress intensity factor, KI. Pre crack will be introduced to the lag screw based on the von mises’ stress result that had been performed in previous research. This paper shows the influence of local mesh arrangement on numerical value of the stress intensity factor, KI obtained by the displacement method. This study aims to simulate the effect of local meshing which is the singularity region on stress intensity factor, KI to the critical point of failure in screw. Five different set of wedges meshing size are introduced during the simulation of finite element analysis. The number of wedges used to simulate this research is 8, 10, 14, 16 and 20. There are three set of numerical equations used to validate the results which are brown and srawley, gross and brown and Tada equation. The result obtained from the finite element software (ANSYS APDL) has a positive agreement with the numerical analysis which is Brown and Srawley compared to other numerical formula. Radius of first row size of 0.014 and singularity element with 14 numbers of wedges is proved to be the best local meshing for this study.

The modified complex variable function method with the continuity conditions of the resultant force and displacement function are used to formulate the hypersingular integral equations (HSIE) for an inclined crack and a circular arc crack lies in the upper part of bonded dissimilar materials subjected to various remote stresses.The curve length coordinate method and appropriate quadrature formulas are used to solve numerically the unknown crack opening displacement (COD) function and the traction along the crack as the right hand term of HSIE.The obtained COD is then used to compute the stress intensity factors (SIF), which control the stability behavior of bodies or materials containing cracks or flaws.Numerical results showed the behavior of the nondimensional SIF at the crack tips.It is observed that the nondimensional SIF at the crack tips depend on the various remote stresses, the elastic constants ratio, the crack geometries and the distance between the crack and the boundary.

An explosion clad plate composed of copper and mild steel is dealt with as a typical example of the bonded dissimilar plate. Tensile tests are carried out using rectangular plate specimens extracted from the clad plate. An artificial through-the-thickness edge crack or parallel cracks are made in each specimen. The cracks are close and perpendicular to the explosive interface. When a tensile load is applied perpendicularly to the crack plane, stress field near cracks, stress intensity factor, crack opening displacement and fracture strength are examined by experiment as well as elasto-plastic finite element analysis. The stress field and the stress intensity factor are evaluated by an infrared stress imaging system. The effects of the material inhomogeneity, the residual stress, the change of material characteristics and the interaction of parallel cracks on the stress field, the stress intensity factor, the deformation behavior and the fracture strength are revealed. The lower strength material ahead of the bonded interface increases the stress intensity factor and the crack opening displacement. The higher strength material ahead of the bonded interface decreases the stress intensity factor and the crack opening displacement. The effects of the material inhomogeneity and the interaction of parallel cracks on the crack opening displacement can be explained from the square of stress intensity factor at the low applied stress level. The brittle fracture strength of the inhomogeneous specimen can be evaluated from the stress intensity factor and linear elastic fracture mechanics criterion.

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e. double) integration strategy that is competitive with its counterpart (i.e. singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C1 continuity, a fact that allows the employment of the popular continuous element technology. The behaviour of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation‐based BEA, for large‐scale problems.

Failure of pressure vessels and piping due to high temperature applications occurs due to the formation of fatigue cracks caused by cyclic load. It is well known that, the consequences of collapses of pipes causing enormous disruption of daily life. Thus there is a need to design and manufacture the pipes with precision and care. The major cause of crack nucleation in pipes is due to corrosion and internal fluid pressure. The crack-tip stresses are determined using stress intensity factor (SIF). In the present work an attempt has been made to determine the SIF for multiple cracks in a circular pipe subjected to internal fluid pressure. Two surface cracks of same size were introduced at the inner wall of the tube. The crack depth ratio (a/t) ranging between 0.1 and 0.5 and crack aspect ratio (a/c) of 0.6 and 1.0 was considered. Internal fluid pressure of 100 MPa was applied at the inner surface of the pipe and the corresponding SIF was measured. SIF values were calculated with consideration of mode-II and mode-III fracture in order to predict the exact SIF. As available SIF solutions of cracked pipes are limited to mode-I fracture, present work presents the influence of additional influence of mode-II and mode-III fracture. It is observed that, as crack depth ratio increases, SIF also increases considerably for semi-circular cracks. Higher SIF values were observed at the crack surface region [S/S 0 = ±1] compared to crack middle [S/S 0 = 0] region. A crossover in SIF was noted at a crack depth ratio of 0.3. At higher crack depths, SIF values decrease at the crack surface region due to additional influence of mode-II and mode-III fracture. In contrast to semi-circular cracks, SIF values are higher at the crack surface region for semi-elliptic cracks irrespective of the crack depths.

On the conventional boundary integral equation formulation for piezoelectric solids with defects or of thin shapes

The boundary element method (BEM) is an efficient numerical method that approximates solutions of various boundary value problems. Despite its success little research has been performed on the conditioning of the linear systems that appear in the BEM. For a Laplace equation with Dirichlet boundary conditions a remarkable phenomenon is observed; the corresponding boundary integral equation (BIE) is singular for a certain critical size of the 2D domain. As a consequence the discrete counterpart of the BIE, the linear system, is singular too, or at least ill-conditioned. This is reflected by the condition number of the system matrix, which is infinitely large, or at least very large. When the condition number of the BEM-matrix is large, the linear system is difficult to solve and the solution of the system is very sensible to perturbations in the boundary data. For a Laplace equation with mixed boundary conditions a similar phenomenon is observed. The corresponding BEM-matrix consists of two blocks; one block originates from the BEM-matrix belonging to the Dirichlet problem, the other block originates from the BEM-matrix belonging to the Neumann problem. The composite matrix inherits the solvability problems from the Dirichlet block. In other words, for the Laplace equation with mixed boundary conditions there exists also a critical size of the 2D domain for which the BEM-matrix has an infinitely large condition number. Hence the size and shape of the domain affects the solvability of the BEM problem. The critical size of the domain for which the BIE becomes singular is related to the logarithmic capacity of the domain. The logarithmic capacity is a positive real number that is a function of the size and shape of the domain. If this logarithmic capacity is equal to one, the domain is a critical domain, and for this domain the BIE becomes singular. Thus by computing the logarithmic capacity we can a-priori determine whether the BIE will be singular or not. The logarithmic capacity depends linearly on the scale of the domain, and thus a domain with logarithmic capacity equal to one can always be found by rescaling the domain. Unfortunately the logarithmic capacity can only be computed analytically for a few simple domains; for more involved domains the logarithmic capacity can be estimated though. There are several possibilities to avoid large condition numbers, i.e. singular BIEs that appear at critical domains. The first option is to rescale the domain such that the logarithmic capacity is unequal to one. One can also add a supplementary condition to the BIE and the linear system. A drawback of this option is that the linear system has more equations than unknowns and different techniques are required to solve the system. A third option is to slightly modify the fundamental solution of the Laplace operator. This fundamental solution directly appears in the BIE and it can be shown that a suitable modification yields BIEs that do not become singular. The critical domains for which the BIEs become singular do not restrict to Laplace equations only. Also for BIEs applied to the biharmonic equation or the elastostatic equations and the Stokes equations such critical domains exist. As the last two equations are vectorial equations, also the corresponding BIE consists of two equations. As a consequence two critical domains can be found for which these BIEs become singular. To obtain nonsingular BIEs techniques similar to the Laplace case can be used. Unfortunately we cannot a-priori determine the sizes for which the BIEs becomes singular, and thus do not know to what size we should rescale the domain to obtain nonsingular BIEs. The existence of critical domains is in essence caused by the logarithmic term in the fundamental solutions for the elliptic boundary value problems in 2D. This logarithmic term does not depend linearly on the size of the domain. When a domain is scaled, i.e. multiplied by a scale factor, the argument of the logarithm is also multiplied by this scale factor, but the logarithm turns this into an additive term. Thus the logarithm transforms multiplication into addition. This affects the BIEs in such a way that critical domains can appear. The fundamental solutions of boundary value problems in 3D do not contain a logarithmic term. Hence scaling of the domain does not affect the fundamental solution, and consequently also the BIE is not affected. Hence we may safely rescale 3D domains without the risk to encounter a critical domain. An example in which a domain takes many different sizes and shapes is the blowing problem. In this problem a viscous fluid is blown to a desired shape. Typically the time is discretised into a set of discrete time steps, and at each step the shape of the fluid is computed by solving the Stokes equations. When attempting to simulate this problem in 2D, we meet a large number of 2D domains, and we risk that one of these domains is equal to or approaches a critical domain. In such a case the BEM will have difficulties with solving the Stokes equations for that particular domain. When simulating the blowing problem in 3D, no critical domains are encountered. It turns out that the BEM is a very efficient numerical method for this particular 3D problem with a moving boundary. As we are merely interested in the shape of the fluid, we only need to know the flow of its boundary. The BEM does exactly that; it does not compute the flow at the interior of the fluid. Furthermore it is rather easy to include other effects from the blowing problem in the model, such as gravity, surface tension and friction from the contact of the fluid with a wall. As only the boundary of the fluid is discretised, the system matrices that appear in the BEM are smaller than the system matrices that appear when solving the problem with a finite element method, for example. Though the BEM-matrices are dense, while the finite element matrices are sparse, the computational effort for the BEM is relatively low. In short, the BEMis a very appropriate numerical method when solving blowing problems.