Multigrid methods for the computation of singular solutions and stress intensity factors II: Crack singularities

We consider the Poisson equation −Δu=f with homogeneous Dirichlet boundary condition on a two-dimensional polygonal domain Ω with cracks. Multigrid methods for the computation of singular solutions and stress intensity factors using piecewise linear functions are analyzed. The convergence rate for the stress intensity factors is whenfεL ²(Ω) and whenfεH ¹(Ω). The convergence rate in the energy norm is in the first case and in the second case. The costs of these multigrid methods are proportional to the number of elements in the triangulation. The general case wherefεH ^m (Ω) is also discussed. The work of the first author was partially supported by NSF under grant DMS-96-00133     

Flux intensity functions for the Laplacian at axisymmetric edges

We derive explicit representation formulas for the computation of flux intensity functions for mixed boundary value problems for the Poisson equation in axisymmetric domains with edges. We rely on the decomposition of the boundary value problems in three dimensions by means of partial Fourier analysis with respect to the rotational angle into boundary value problems in the two‐dimensional meridian domain of . Utilizing smooth cutoff functions, the solutions of the reduced problems are analyzed semi‐analytically near corners of the plane meridian domain, and the edge flux intensity functions are constructed via Fourier synthesis and convergence analysis. The formulas are also applicable in the case of crack fronts. The constructive nature of the formulas provides in a straightforward way an efficient strategy for the accurate computation of edge flux intensity functions in axisymmetric domains. A demonstration example that illustrates the application of the formulas is presented. Copyright © 2012 John Wiley & Sons, Ltd.     

Removable isolated singularities of solutions to the weighted p‐Laplacian with singular convection

This paper is concerned with the removability of isolated singularities for the weighted p‐Laplacian with singular convection. Critical singular exponents included in diffusion and convection are exactly determined. Copyright © 2012 John Wiley & Sons, Ltd.     

Decomposition in regular and singular parts (in domains with corners) and stability under perturbation of the geometry

We analyze the Dirichlet problem for the Laplacian in a polygonal domain where boundary and angles depend on a parameter. We use the boundary integral equation, localization and Mellin transformation techniques to show that the solution has a decomposition in regular and singular parts which blow up at certain exceptional angles. We derive a modified decomposition which depends continuously on the angle.     

Closed-form solutions for a rectangular layer with central crack under anti-plane shear stress

We consider the problem of determining the stress distributionin a finite rectangular elastic layer containing a Griffithcrack which is opened by internal shear stress acting alongthe length of the crack. The mode III crack is assumed to belocated in the middle plane of the rectangular layer. The followingtwo problems are considered: (A) the central crack is perpendicularto the two fixed lateral surfaces and parallel to the othertwo stress-free surfaces; (B) all the lateral surfaces of therectangular layer are clamped and the central crack is parallelto the two lateral surfaces. By using Fourier transformations,we reduce the solution of each problem to the solution of dualintegral equations with sine kernels and a weight function whichare solved exactly. Finally, we derive closed-form expressionsfor the stress intensity factor at the tip of the crack andthe numerical values for the stress intensity factor at theedges of the cracks are presented in the form of tables.     

Existence of periodic solutions for second-order differential equations with singularities and the strong force condition

We prove the existence of periodic solutions for the equation

(1)     

Prescribing analytic singularities for solutions of a class of vector fields on the torus

We consider the operator acting on distributions on the two-torus where and are real-valued, real analytic functions defined on the unit circle We prove, among other things, that when changes sign, given any subset of the set of the local extrema of the local primitives of there exists a singular solution of such that the projection of its analytic singular support is furthermore, for any and any closed subset of there exists such that and We also provide a microlocal result concerning the trace of at

     

On the numerical solution of Cauchy type singular integral equations and the determination of stress intensity factors in case of complex singularities

Summary A Cauchy type singular integral equation along a finite real interval and with a weight function with complex singularities at the end-points of the integration interval can be numerically solved by reduction to a system of linear equations, by using an appropriate numerical integration rule associated with the Jacobi polynomials, in exactly the same way used for the case of real singularities. For the numerical solution of such an equation arising in plane elasticity crack problems and the evaluation of stress intensity factors at crack tips, the Lobatto-Jacobi numerical integration rule is the most appropriate.

---

Résumé Une équation intégrale du type de Cauchy singulière le long d'un intervalle fini réel et avec une fonction pondérante ayant des singularités complexes aux extremités de l'intervalle d'intégration peut être résolue numériquement par réduction à un système d'équations linéaires, en utilisant une règle appropriée d'intégration numérique associée aux polynômes de Jacobi, exactement de la même manière que dans le cas des singularités réelles. La façon la plus appropriée de trouver la solution numérique de cette équation telle qu'elle se présente dans les problèmes de fissure en élasticité plane et d'évaluer les facteurs de contrainte aux extremités da la fissure est la méthode d'intégration numérique de Lobatto-Jacobi.

     

Computation of mixed mode stress intensity factors in a four-point bend specimen

Four-point bend specimen is one of the most important specimens of the fracture mechanics because it can produce mixed modes I and II. Therefore, computation of stress intensity factors in this specimen is of practical interest. Several relations have been suggested that no one of them has completely considered the effects of the loading point and crack geometry. In this paper, mixed mode stress intensity factors of the bend specimen are computed by finite element method (FEM) and after validating by comparing with the available results in the literature, the results will be assessed to determine the effects of different crack location and loading distances from the middle of the specimen. Finally, two new coefficients comprising these effects are introduced.     

Influence of inclusion on the stress intensity factors under thermomechanical loads in a PM superalloy

The interaction between a round inclusion and a crack under thermomechanical loading is analyzed based on a modified body force method. The traction-free condition on the crack line is mended by adding the resultant force induced by thermal stress to the force equilibrium equations, so that the coupling of mechanical and thermal loads could be taken into account. The series of integral equations can be discretized to a set of linear equations. Stress intensity factors (SIFs) are obtained through solving the linear equations. The calculated results in this paper are compared to those in open references to validate the method and code. The method is applied to a case of FGH95 PM superalloy containing Al₂O₃ inclusions under mechanical and thermal loads. The results show that the thermal load has little effect on SIF, while the mechanical load is the dominant factor.     

Analysis of the singular function boundary integral method for a biharmonic problem with one boundary singularity

In this article, we analyze the singular function boundary integral method (SFBIM) for a two‐dimensional biharmonic problem with one boundary singularity, as a model for the Newtonian stick‐slip flow problem. In the SFBIM, the leading terms of the local asymptotic solution expansion near the singular point are used to approximate the solution, and the Dirichlet boundary conditions are weakly enforced by means of Lagrange multiplier functions. By means of Green's theorem, the resulting discretized equations are posed and solved on the boundary of the domain, away from the point where the singularity arises. We analyze the convergence of the method and prove that the coefficients in the local asymptotic expansion, also referred to as stress intensity factors, are approximated at an exponential rate as the number of the employed expansion terms is increased. Our theoretical results are illustrated through a numerical experiment. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Existence and multiplicity of positive solutions for some Hardy-Sobolev critical elliptic equation with boundary singularities

In this paper, by Ekeland’s variational principle and strong maximum principle, we consider the existence and multiplicity of positive solutions for some semilinear elliptic equation involving critical Hardy-Sobolev exponents and Hardy terms with boundary singularities.     

Existence and maximal Lp-regularity of solutions for the porous medium equation on manifolds with conical singularities

We consider the porous medium equation on manifolds with conical singularities and show existence, uniqueness, and maximal L^p-regularity of a short-time solution. In particular, we obtain information on the short time asymptotics of the solution near the conical point. Our method is based on bounded imaginary powers results for cone differential operators on Mellin–Sobolev spaces and R-sectoriality perturbation techniques.     

Stress intensity factors for the Brazilian disc with a short central crack: Opening versus closing cracks

Closed form expressions are obtained for the stress intensity factors (SIFs) in case of a Brazilian disc with a short central crack, the length of which does not exceed one fifth of the disc radius. The disc is loaded by uniform radial pressure along two finite symmetric arcs of its periphery. The solution is achieved using the method of complex potentials introduced by Kolosov and Muskhelishvili. The advantage of the expressions obtained is that they are valid both for cracks under opening mode as well as for closing cracks. For the first case (opening cracks) the results of the present study are compared with existing approximate solutions and it is concluded that the agreement is excellent as long as the length of the crack remains relatively small compared to the radius of the disc. Regarding the case of a closing mode crack the procedure proposed here (based on a recent alternative approach of the cracked Brazilian disc) leads to a physically acceptable deformed crack shape instead to an unnatural crack with overlapped lips. At the same moment the dependence of the SIFs on the properties of the material is eliminated.     

Stress intensity factors for a finite plate with an inclined crack by boundary collocation

In this paper, we combine the Muskhelishvili＇s complex variable method and boundary collocation method, and choose a set of new stress function based on the stress boundary condition of crack surface, the higher precision and less computation are reached. This method is applied to calculating the stress intensity factor for a finite plate with an inclined crack. The influence of θ （the obliquity of crack） on the stress intensity factors, as well as the number of summation terms on the stress intensity factor are studied and graphically represented.     

Superelements for the finite element solution of two-dimensional elliptic problems with boundary singularities

A novel singular superelement (SSE) formulation has been developed to overcome the loss of accuracy encountered when applying the standard finite element schemes to two-dimensional elliptic problems possessing a singularity on the boundary arising from an abrupt change of boundary conditions or a reentrant corner. The SSE consists of an inner region over which the known analytic form of the solution in the vicinity of the singular point is utilized, and a transition region in which blending functions are used to provide a smooth transition to the usual linear or quadratic isoparametric elements used over the remainder of the domain. Solution of the finite element equations yield directly the coefficients of the asymptotic series, known as the flux/stress intensity factors in linear heat transfer or elasticity theories, respectively. Numerical examples using the SSE for the Laplace equation and for computing the stress intensity factors in the linear theory of elasticity are given, demonstrating that accurate results can be attained for a moderate computational effort.     

Multigrid preconditioned Krylov subspace methods for three-dimensional numerical solutions of the incompressible Navier–Stokes equations

We consider numerical methods for the incompressible Reynolds averaged Navier–Stokes equations discretized by finite difference techniques on non-staggered grids in body-fitted coordinates. A segregated approach is used to solve the pressure–velocity coupling problem. Several iterative pressure linear solvers including Krylov subspace and multigrid methods and their combination have been developed to compare the efficiency of each method and to design a robust solver. Three-dimensional numerical experiments carried out on scalar and vector machines and performed on different fluid flow problems show that a combination of multigrid and Krylov subspace methods is a robust and efficient pressure solver. This revised version was published online in June 2006 with corrections to the Cover Date.     

A boundary element analysis for stress intensity factors of multiple circular arc cracks in a plane elasticity plate

In this paper, a numerical approach for analyzing interacting multiple cracks in infinite linear elastic media is presented. By extending Bueckner’s principle suited for a crack to a general system containing multiple interacting cracks, the original problem is divided into a homogeneous problem (the one without cracks) subjected to remote loads and a multiple crack problem in an unloaded body with applied tractions on the crack surfaces. Thus, the results in terms of the stress intensity factors (SIFs) can be obtained by considering the latter problem, which is analyzed easily by means of the displacement discontinuity method with crack-tip elements proposed recently by the author. Test examples are given to illustrate that the numerical approach is very accurate for analyzing interacting multiple cracks in an infinite linear elastic media under remote uniform stresses. In addition, the displacement discontinuity method with crack-tip elements is used to analyze a multiple crack problem in a finite plate. It is found that the boundary element method is also very accurate for investigating interacting multiple cracks in a finite plate. Specially, a generalization of Bueckner’s principle and the displacement discontinuity method with crack-tip elements are used to analyze multiple circular arc crack problems in infinite plate in tension (including: Two Collinear Circular Arc Cracks, Three Collinear Circular Arc Cracks, Two Parallel Circular Arc Cracks, Three Parallel Circular Arc Cracks and Two Circular Arc Cracks) in a plane elasticity plate. Many results are given.