Hybrid extended displacement discontinuity-charge simulation method for analysis of cracks in 2D piezoelectric media

The extended displacement discontinuity method (EDDM) and the charge simulation method (CSM) are combined to develop an efficient approach for analysis of cracks in two-dimensional piezoelectric media. In the proposed hybrid EDD–CSM, the solution for an electrically impermeable crack is approximately expressed by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with the sources placed at chosen points outside the domain of the problem under consideration and the extended Crouch fundamental solutions with the extended displacement discontinuities placed on the crack. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy the conditions on the boundary of the domain and on the crack face. Furthermore, the hybrid EDD–CSM is applied to solve the problems of cracks under electrically permeable condition, as well as under semi-permeable conditions by using an iterative approach. Two important crack problems in fracture mechanics, the center cracks and the edge cracks in piezoelectric strips, are analyzed by the proposed method. The stress intensity factor and the electric displacement intensity factor are calculated. Meanwhile the effects of strip size and the electric boundary conditions on these intensity factors are studied.     

An effective method of stress intensity factor calculation for cracks emanating from a triangular or square hole under biaxial loads

This paper is concerned with stress intensity factors for cracks emanating from a triangular or square hole under biaxial loads by means of a new boundary element method. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfied and the crack‐tip displacement discontinuity elements proposed by the author. In the boundary element implementation, the left or the right crack‐tip displacement discontinuity element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. The method is called a Hybrid Displacement Discontinuity Method (HDDM). Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors for plane elastic crack problems. In addition, the present numerical results can reveal the effect of the biaxial loads on stress intensity factors.     

A numerical analysis of cracks emanating from a square hole in a rectangular plate under biaxial loads

This note concerns with stress intensity factors of cracks emanating from a square hole in rectangular plate under biaxial loads by means of the boundary element method which consists of the non-singular displacement discontinuity element presented by Crouch and Starfied and the crack tip displacement discontinuity elements proposed by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundary. The present numerical results illustrate that the present approach is very effective and accurate for calculating stress intensity factors of complicated cracks in a finite plate and can reveal the effect of the biaxial load and the cracked body geometry on stress intensity factors.     

Stress intensity factors for cracks emanating from a triangular or square hole in an infinite plate by boundary elements

This paper concerns stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure calculated by means of a boundary element method, which consists of constant displacement discontinuity element presented by Crouch and Starfied and crack tip displacement discontinuity elements proposed by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors of plane elasticity crack problems. Specifically, the numerical results of stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure are given.     

Three-dimensional BEM for piezoelectric fracture analysis

A boundary element approach with quadratic isoparametric elements, quarter-point elements and singular quarter-point elements for three-dimensional crack problems in piezoelectric solids under mechanical and electrical loading conditions, is presented in this paper for the first time. The procedure is based on Deeg's fundamental solution for anisotropic piezoelectric materials, and the classical extended displacement boundary integral equation. Stress and electric displacement intensity factors are directly evaluated as system unknowns, and also as functions of the computed nodal displacements and electric potentials at crack faces. Special attention is paid to the fundamental solution evaluation. Several three-dimensional crack problems in transversely isotropic bodies under mechanical and electrical loading conditions are analysed. Numerical solutions computed for prismatic cracked 3D plate problems with a plane strain behaviour are in very good agreement with their corresponding 2D BE solutions. Results for a penny shape crack in a piezoelectric cylinder are presented for the first time. The proposed approach is shown to be a simple, robust and useful tool for stress and electric displacement intensity factors evaluation in piezoelectric media.     

Extended displacement discontinuity Green's functions for three-dimensional transversely isotropic magneto-electro-elastic media and applications

A versatile method is presented to derive the extended displacement discontinuity Green's functions or fundamental solutions by using the integral equation method and the Green's functions of the extended point forces. In particular, the three-dimensional (3D) transversely isotropic magneto-electro-elastic problem is used to demonstrate the method. On this condition, the extended displacement discontinuities include the elastic displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, while the extended forces include the point forces, the point electric charge and the point electric current. Based on the obtained Green's functions, the extended Crouch fundamental solutions are derived and an extended displacement discontinuity method is developed for analysis of cracks in 3D magneto-electro-elastic media. The extended intensity factors of two coplanar and parallel rectangular cracks are calculated under impermeable boundary condition to illustrate the application, accuracy and efficiency of the proposed method.     

Analysis of an arbitrarily oriented crack in a finite piezoelectric plane via the hybrid extended displacement discontinuity-fundamental solution method

In this paper, we analyze an arbitrarily oriented crack in a finite two-dimensional piezoelectric medium with the polarization saturation model near the crack tip. We first derive the extended Green’s functions corresponding to the extended point-displacement discontinuities of an arbitrarily oriented crack based on the Green’s functions of the extended point forces and the Somigliana identity. Then, the extended field intensity factors and the local J-integral near the crack tip are expressed in terms of the extended displacement discontinuity on crack faces. Finally, the nonlinear hybrid extended displacement discontinuity-fundamental solution method is proposed to analyze an electrically nonlinear crack in a finite piezoelectric medium. Numerical examples are carried out for both linear and nonlinear fracture models of the crack under electrically impermeable boundary conditions. The influence of the crack orientation and geometric size on the fracture behaviors of the crack is investigated.     

Wave propagation in cracked elastic slabs and half-space domains—TBEM and MFS approaches

In this paper, the traction boundary element method (TBEM) and the method of fundamental solutions (MFS), formulated in the frequency domain, are used to evaluate the 3D scattered wave field generated by 2D empty cracks embedded in an elastic slab and a half-space. Both models overcome the thin-body difficulty posed when the classical BEM is applied.The crack exhibits arbitrary cross section geometry and null thickness. In neither model are the horizontal formation surfaces discretized, since appropriate fundamental solutions are used to take them into consideration.The TBEM models the crack as a single line. The singular and hypersingular integrals that arise during the TBEM model's implementation are computed analytically, which overcomes one of the drawbacks of this formulation. The results provided by the proposed TBEM model are verified against responses provided by the classical BEM models derived for the case of an empty cylindrical circular cavity.The MFS solution is approximated in terms of a linear combination of fundamental solutions, generated by a set of virtual sources simulating the scattered field produced by the crack, using a domain decomposition technique. To avoid singularities, these fictitious sources are not placed close to the crack, and the use of an enriched function to model the displacement jumps across the crack is unnecessary.The performances of the proposed models are compared and their limitations are shown by solving the case of a C-shaped crack embedded in an elastic slab and a half-space domain.The applicability of these formulations is illustrated by presenting snapshots from computer animations in the time domain for an elastic slab containing an S-shaped crack, after applying an inverse Fourier transformation to the frequency domain computations.     

A hypersingular time‐domain BEM for 2D dynamic crack analysis in anisotropic solids

A hypersingular time‐domain boundary element method (BEM) for transient elastodynamic crack analysis in two‐dimensional (2D), homogeneous, anisotropic, and linear elastic solids is presented in this paper. Stationary cracks in both infinite and finite anisotropic solids under impact loading are investigated. On the external boundary of the cracked solid the classical displacement boundary integral equations (BIEs) are used, while the hypersingular traction BIEs are applied to the crack‐faces. The temporal discretization is performed by a collocation method, while a Galerkin method is implemented for the spatial discretization. Both temporal and spatial integrations are carried out analytically. Special analytical techniques are developed to directly compute strongly singular and hypersingular integrals. Only the line integrals over an unit circle arising in the elastodynamic fundamental solutions need to be computed numerically by standard Gaussian quadrature. An explicit time‐stepping scheme is obtained to compute the unknown boundary data including the crack‐opening‐displacements (CODs). Special crack‐tip elements are adopted to ensure a direct and an accurate computation of the elastodynamic stress intensity factors from the CODs. Several numerical examples are given to show the accuracy and the efficiency of the present hypersingular time‐domain BEM. Copyright © 2008 John Wiley & Sons, Ltd.     

基于X-SBFEM的裂纹体非网格重剖分耦合模型研究

扩展有限元法利用了非网格重剖分技术,但需要基于裂尖解析解构造复杂的插值基函数,计算精度受网格疏密和插值基函数等因素影响。比例边界有限元法则在求解无限域和裂尖奇异性问题优势明显,两者衔接于有限元法理论内,可建立一种结合二者优势的断裂耦合数值模型。该文从虚功原理出发,利用位移协调与力平衡机制,提出了一种断裂计算的新方法X-SBFEM,达到了扩展有限元模拟裂纹主体、比例边界有限元模拟裂尖的目的。在数值算例中,通过边裂纹和混合型裂纹的应力强度因子计算,并与理论解对比,验证了该方法的准确性和有效性。     

The exact form and properties of the deformed transverse internal elastic crack

The exact shape and properties of an internal transverse crack in an infinite plate under conditions of plane stress submitted to a biaxial normal loading at infinity are presented. The shape of the deformed Griffith crack as this was defined by the exact theoretical solution of the problem is an ellipse whose properties depend on the mode of loading of the plate and its mechanical properties. This result is different to the result arising from the respective singular and two-term solutions, which define a parabolic shape for the deformed crack. Furthermore, while the singular solution implies that the tip of the crack remains undisplaced during loading, the two-term and the exact solutions define its displacement, which is the same for both cases. In this way the two-term solution constitutes an intermediate case, which at the crack tip coincides with the exact solution, whereas at the middle of the crack coincides with the singular solution. Interesting results and comparisons were also derived for the curvature of the deformed crack at its tip, which was shown to differ considerably according to the different solutions.     

Transient analysis of dynamic crack propagation in piezoelectric materials

Abstract

In this paper, the transient analysis of semi‐infinite propagating cracks in piezoelectric materials subjected to dynamic anti‐plane concentrated body force is investigated. The crack surface is assumed to be covered with an infinitesimally thin, perfectly conducting electrode that is grounded. In analyzing this problem, it has characteristic lengths and a direct attempt towards solving this problem by transform and Wiener‐Hopf techniques (Noble, 1958) is not applicable. In order to solve this problem, a new fundamental solution for propagating cracks in piezoelectric materials is first established and the transient response of the propagating crack is obtained by superposition of the fundamental solution in the Laplace transform domain. The fundamental solution to be used is the responses of applying exponentially distributed traction in the Laplace transform domain on the propagating crack surface. Taking into account the quasi‐static approximation, exact analytical transient solutions for the dynamic stress intensity factor and the dynamic electric displacement intensity factor are obtained by using the Cagniard‐de Hoop method (Cagnard, 1939; de Hoop, 1960) of Laplace inversion and are expressed in explicit forms. Numerical calculations of dynamic intensity factors are evaluated and the results are discussed in detail. The transient solutions for stationary cracks have been shown to approach the corresponding static values after the shear wave of the piezoelectric material has passed the crack tip.     

A numerical Green's function and dual reciprocity BEM method to solve elastodynamic crack problems

A computational model based on the numerical Green's function (NGF) and the dual reciprocity boundary element method (DR-BEM) is presented for the study of elastodynamic fracture mechanics problems. The numerical Green's function, corresponding to an embedded crack within the infinite medium, is introduced into a boundary element formulation, as the fundamental solution, to calculate the unknown external boundary displacements and tractions and in post-processing determine the crack opening displacements (COD). The domain inertial integral present in the elastodynamic equation is transformed into a boundary integral one by the use of the dual reciprocity technique. The dynamic stress intensity factors (SIF), computed through crack opening displacement values, are obtained for several numerical examples, indicating a good agreement with existing solutions.     

Precise elastic-plastic analysis of crack line field for mode II plane strain crack

The near crack line analysis method has been used to investigate the exact elastic-plastic solutions of a mode II crack under plane strain condition in an elastic-perfectly plastic solid. The significance of this paper is that the assumptions of the conventional small scale yielding theory have been completely abandoned. The inappropriateness of matching conditions formerly taken at the elastic-plastic boundary ths been corrected as well. By eatching the general solution of the plastic stress (but not the special solution that was adopted) with the exact elastic stresses (but not the crack tip K-dominant field) at the elastic-plastic boundary near the crack line, the plastic stresses, the length of the plastic zone and the unit normal vector of the elastic-plastic boundary, which are sufficiently precise near the crack line region, have been given. The solutions are suitable not only under the condition that the plastic region is sufficiently small but also under the condition that the plastic region is large.     

Nanoscale anti-plane cracking of materials with consideration of bulk and surface piezoelectricity effects

The influence of surface effect, including surface elasticity and surface piezoelectricity, on thefracture behavior of piezoelectricmaterials with an anti-plane crack is studied.Based on the coupled surface andinterface elasticity model, the solutions to the problem are obtained by applying the singular integral method. Bycomparing the solutions influenced by the surface piezoelectricity with those affected by the surface elasticity,it is found that the influence of the surface piezoelectricity on the crack opening displacement, the crackelectric potential jump across the crack center, the crack tip stress and electric displacement intensity factorscannot be ignored. Under various electrical boundary conditions, the influence of surface piezoelectricity onthe sliding displacement, crack tip stress and electric displacement intensity factors exhibits the same tendency.Besides, the influence of surface piezoelectricity on the electric displacement intensity factor is independentof the electrical boundary conditions, which is different from the results where only the surface elasticity isconsidered.     

Stress intensity factor computation using the method of fundamental solutions: mixed‐mode problems

The method of fundamental solutions is applied to the computation of stress intensity factors in linear elastic fracture mechanics. The displacements are approximated by linear combinations of the fundamental solutions of the Cauchy–Navier equations of elasticity and the leading terms for the displacement near the crack tip. Two algorithms are developed, one using a single domain and one using domain decomposition. Numerical results are given. Copyright © 2006 John Wiley & Sons, Ltd.     

A new form of exact solutions for mode I,II, III crack problems and implications

A new form of an exact linear elastic solution is obtained for the problem of a crack in a semi-infinite plate subjected, at infinity, to antiplane stress (Mode III) loading. The use of the conformal mapping technique results in a convenient stress and displacement solution for each point of the semi-infinite domain. For completeness, the solution technique is extended to Mode I, II problems with center-cracks as well as the Mode III V-notch problem. It is shown that the limiting case of the V-notch collapses to the stress functions independently derived for the edge-cut. At distances equivalent to 10% of the crack length away from the crack tip, the exact solutions give stresses about 7.5% greater than the one-term results. Ramifications of the exact solutions to finite-element solutions, elastic-plastic and diffusion problems are discussed.     

Analytical forms of higher-order asymptotic elastic-plastic crack-tip fields in a linear hardening material under antiplane shear

An analytical study of the higher-order asymptotic solutions of the stress and strain fields near the traction-free crack tip under antiplane shear in a linear hardening material is investigated. The results show that every term of the asymptotic fields is controlled by both elasticity and plasticity and all the higher-order asymptotic fields are governed by linear nonhomogeneous equations. The first four term solutions are presented analytically and the first four terms are described by two independent parameters J and K ₂. The amplitude of the second order term solution is only dependent on the material properties, but independent of loading and geometry. This paper focuses on the case with traction-free crack surface boundary conditions. The effects of different crack surface boundary conditions, such as clamped and mixed surfaces, on the crack-tip fields are also presented. Comparison of multi-term solution with leading term solution, and finite element solution in an infinite strip with semi-infinite crack under constant displacements along the edges is provided.     

Numerical simulations of cracked plate using XIGA under different loads and boundary conditions

This article deals with the numerical simulation of cracked plate using extended isogeometric analysis (XIGA) under different loads and boundary conditions. The plate formulation is done using first-order shear deformation theory. The crack faces are modeled by the Heaviside function, whereas the singularity in stress field at the crack tip is modeled by crack tip enrichment functions. The stress intensity factors for the cracked plate are numerically computed using a domain-based interaction integral. The results obtained by XIGA for the center and edge crack plate are compared with extended finite element method and/or literature results for different types of loads and boundary conditions.