Boundary Element Method Analysis of Three-Dimensional Thermoelastic Fracture Problems Using the Energy Domain Integral

[+] Author and Article Information
R. Balderrama, M. Martinez

Mechanical Engineering School, Universidad Central de Venezuela, Caracas, Venezuela

A. P. Cisilino1

Department of Mechanical Engineering, Welding and Fracture Division – INTEMA – CONICET, Universidad Nacional de Mar del Plata, Av. Juan B. Justo 4302, 7600 Mar del Plata, Argentinacisilino@fi.mdp.edu.ar


To whom correspondence should be addressed.

J. Appl. Mech 73(6), 959-969 (Dec 21, 2005) (11 pages) doi:10.1115/1.2173287 History: Received June 07, 2005; Revised December 21, 2005

A boundary element method (BEM) implementation of the energy domain integral (EDI) methodology for the numerical analysis of three-dimensional fracture problems considering thermal effects is presented in this paper. The EDI is evaluated from a domain representation naturally compatible with the BEM, since stresses, strains, temperatures, and derivatives of displacements and temperatures at internal points can be evaluated using the appropriate boundary integral equations. Special emphasis is put on the selection of the auxiliary function that represents the virtual crack advance in the domain integral. This is found to be a key feature to obtain reliable results at the intersection of the crack front with free surfaces. Several examples are analyzed to demonstrate the efficiency and accuracy of the implementation.

Copyright © 2006 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 3

(a) General cracked body with mechanical and thermal boundary conditions. (b) Crack discretization strategy.

Grahic Jump Location
Figure 4

Schematic of the volume cells in the crack front region illustrating the virtual crack extensions for a corner node, a mid-node, and a surface node

Grahic Jump Location
Figure 5

Boundary element discretization and integration cells

Grahic Jump Location
Figure 6

Influence of parameters wp, wnp, and β on the shape of function φ (one-dimensional case)

Grahic Jump Location
Figure 7

Geometry, dimensions, and boundary conditions for the edge and center cracked specimens

Grahic Jump Location
Figure 1

(a) Definition of the local orthogonal Cartesian coordinates at point η on the crack front. (b) Virtual crack front advance.

Grahic Jump Location
Figure 2

Tubular domain surrounding a segment of the crack front

Grahic Jump Location
Figure 10

Normalized mode I stress intensity factor along the crack front for the central crack in a thick panel: (a) results using bi-quadratic φ and (b) results using optimized φ

Grahic Jump Location
Figure 8

Geometry, dimensions, crack discretization, and boundary conditions for the penny-shaped and annular cracks

Grahic Jump Location
Figure 9

Normalized mode I stress intensity factor along the crack front for the edge crack in a thick panel: (a) results using bi-quadratic φ and (b) results using optimized φ



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In