Boundary Element Analysis of Three-Dimensional Interface Cracks in Transversely Isotropic Bimaterials Using the Energy Domain Integral

2010 ◽  
Vol 454 ◽  
pp. 47-77 ◽  
Author(s):  
N.O. Larrosa ◽  
J.E. Ortiz ◽  
A.P. Cisilino

It is presented in this paper a three-dimensional Boundary Element Method (BEM) implementation of the Energy Domain Integral for the fracture mechanical analysis of three-dimensional interface cracks in transversely isotropic bimaterials. The J-integral is evaluated using a domain representation naturally compatible with the BEM, in which the stresses, strains and derivatives of displacements at internal points are evaluated using their appropriate boundary integral equations. Several examples are solved and the results compared with those available in the literature to demonstrate the efficiency and accuracy of the implementation to solve straight and curved crack-front problems.

2005 ◽  
Vol 73 (6) ◽  
pp. 959-969 ◽  
Author(s):  
R. Balderrama ◽  
A. P. Cisilino ◽  
M. Martinez

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.


2018 ◽  
Vol 183 ◽  
pp. 01042 ◽  
Author(s):  
Igor Vorobtsov ◽  
Aleksandr Belov ◽  
Andrey Petrov

The development of time-step boundary-element scheme for the three dimensional boundaryvalue problems of poroelastodynamics is presented. The poroelastic continuum is described using Biot’s mathematical model. Poroelastic material is assumed to consist of a solid phase constituting an elastic formdefining skeleton and carrying most of the loading, and two fluid phases filling the pores. Dynamic equations of the poroelastic medium are written for unknown functions of displacement of the elastic skeleton and pore pressures of the filling materials. Green’s matrices and, based on it, boundary integral equations are written in Laplace domain. Discrete analogue are obtained by applying the collocation method to a regularized boundary integral equation. Boundary element scheme is based on time-step method of numerical inversion of Laplace transform. A modification of the time-step scheme on the nodes of Runge-Kutta methods is considered. The Runge-Kutta scheme is exemplified with 2-and 3-stage Radau schemes. The results of comparing the two schemes in analyzing a numerical example are presented.


2017 ◽  
Vol 29 (6) ◽  
pp. 1255-1271 ◽  
Author(s):  
MingHao Zhao ◽  
Yuan Li ◽  
CuiYing Fan

An arbitrarily shaped planar crack under different thermal and electric boundary conditions on the crack surfaces is studied in three-dimensional transversely isotropic thermopiezoelectric media subjected to thermal–mechanical–electric coupling fields. Using Hankel transformations, Green functions are derived for unit point extended displacement discontinuities in three-dimensional transversely isotropic thermopiezoelectric media, where the extended displacement discontinuities include the conventional displacement discontinuities, electric potential discontinuity, as well as the temperature discontinuity. On the basis of these Green functions, the extended displacement discontinuity boundary integral equations for arbitrarily shaped planar cracks in the isotropic plane of three-dimensional transversely isotropic thermopiezoelectric media are established under different thermal and electric boundary conditions on the crack surfaces, namely, the thermally and electrically impermeable, permeable, and semi-permeable boundary conditions. The singularities of near-crack border fields are analyzed and the extended stress intensity factors are expressed in terms of the extended displacement discontinuities. The effect of different thermal and electric boundary conditions on the extended stress intensity factors is studied via the extended displacement discontinuity boundary element method. Subsequent numerical results of elliptical cracks subjected to combined thermal–mechanical–electric loadings are obtained.


2021 ◽  
Vol 83 (1) ◽  
pp. 76-86
Author(s):  
A.A. Belov ◽  
A.N. Petrov

The application of non-classical approach of the boundary integral equation method in combination with the integral Laplace transform in time to anisotropic elastic wave modeling is considered. In contrast to the classical approach of the boundary integral equation method which is successfully implemented for solving three-dimensional isotropic problems of the dynamic theory of elasticity, viscoelasticity and poroelasticity, the alternative nonclassical formulation of the boundary integral equations method is presented that employs regular Fredholm integral equations of the first kind (integral equations on a plane wave). The construction of such boundary integral equations is based on the structure of the dynamic fundamental solution. The approach employs the explicit boundary integral equations. The inverse Laplace transform is constructed numerically by the Durbin method. A numerical solution of the dynamic problem of anisotropic elasticity theory based on the boundary integral equations method in a nonclassical formulation is presented. The boundary element scheme of the boundary integral equations method is built on the basis of a regular integral equation of the first kind. The problem is solved in anisotropic formulation for the load acting along the normal in the form of the Heaviside function on the cube face weakened by a cubic cavity. The obtained boundary element solutions are compared with finite element solutions. Numerical results prove the efficiency of using boundary integral equations on a single plane wave in solving three-dimensional anisotropic dynamic problems of elasticity theory. The convergence of boundary element solutions is studied on three schemes of surface discretization. The achieved calculation accuracy is not inferior to the accuracy of boundary element schemes for classical boundary integral equations. Boundary element analysis of solutions for a cube with and without a cavity is carried out.


2008 ◽  
Vol 383 ◽  
pp. 109-121 ◽  
Author(s):  
E. Pineda ◽  
M.H. Aliabadi

This paper presents the development of a new boundary element formulation for analysis of fracture problems in creeping materials. For the creep crack analysis the Dual Boundary Element Method (DBEM), which contains two independent integral equations, was formulated. The implementation of creep strain in the formulation is achieved through domain integrals in both boundary integral equations. The domain, where the creep phenomena takes place, is discretized into quadratic quadrilateral continuous and discontinuous cells. The creep analysis is applied to metals with secondary creep behaviour. This is con…ned to standard power law creep equations. Constant applied loads are used to demonstrate time e¤ects. Numerical results are compared with solutions obtained from the Finite Element Method (FEM) and others reported in the literature.


2016 ◽  
Vol 685 ◽  
pp. 172-176 ◽  
Author(s):  
Leonid Igumnov ◽  
S.Yu. Litvinchuk ◽  
A.A. Belov ◽  
A.A. Ipatov

Problems of the poroviscodynamics are considered. Theory of poroviscoelasticity is based on Biot’s equations of fluid saturated porous media under assumption that the skeleton is viscoelastic. Viscoelastic effects of solid skeleton are modeled by mean of elastic-viscoelastic correspondence principle, using such viscoelastic models as a standard linear solid model and model with weakly singular kernel. The fluid is taken as original in Biot’s formulation without viscoelastic effects. Boundary integral equations method is applied to solve three-dimensional boundary-value problems. Boundary-element method with mixed discretization and matched approximation of boundary functions is used. Solution is obtained in Laplace domain, and then Durbin’s algorithm of numerical inversion of Laplace transform is applied to perform solution in time domain. An influence of viscoelastic parameters on dynamic responses is studied. Numerical example of the surface waves modelling is considered.


2017 ◽  
Vol 11 (4) ◽  
pp. 308-312
Author(s):  
Iaroslav Pasternak ◽  
Heorhiy Sulym

AbstractThe paper presents novel boundary element technique for analysis of anisotropic thermomagnetoelectroelastic solids containing cracks and thin shell-like soft inclusions. Dual boundary integral equations of heat conduction and thermomagnetoelectroelasticity are derived, which do not contain volume integrals in the absence of distributed body heat and extended body forces. Models of 3D soft thermomagnetoelectroelastic thin inclusions are adopted. The issues on the boundary element solution of obtained equations are discussed. The efficient techniques for numerical evaluation of kernels and singular and hypersingular integrals are discussed. Nonlin-ear polynomial mappings are adopted for smoothing the integrand at the inclusion’s front, which is advantageous for accurate evaluation of field intensity factors. Special shape functions are introduced, which account for a square-root singularity of extended stress and heat flux at the inclusion’s front. Numerical example is presented.


Sign in / Sign up

Export Citation Format

Share Document