Fracture of two three-dimensional parallel internal cracks in brittle solid under ultrasonic fracturing

Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both hard rock drilling and oil and gas recovery. A three-dimensional internal laser-engraved crack(3D-ILC) method was introduced to prefabricate two parallel internal cracks within the samples without any damage to the surface. The samples were subjected to UF. The mechanism of UF was elucidated by analyzing the characteristics of fracture surfaces. The crack propagation path under different ultrasonic parameters was obtained by numerical simulation based on the Paris fatigue model and compared to the experimental results of UF. The results show that the 3D-ILC method is a powerful tool for UF research.Under the action of an ultrasonic field, the fracture surface shows the characteristics of beach marks and contains powder locally, indicating that the UF mechanism includes high-cycle fatigue fracture, shear and friction, and temperature load. The two internal cracks become close under UF. The numerical result obtained by the Paris fatigue model also shows the attraction of the two cracks, consistent with the test results. The 3D-ILC method provides a new tool for the experimental study of UF. Compared to the conventional numerical methods based on the analysis of stress-strain and plastic zone, numerical simulation can be a good alternative method to obtain the crack path under UF.

Discrete Element Modeling of Asphalt Concrete Cracking Using a User-def ined Three-dimensional Micromechanical Approach

We established a user-defined micromechanical model using discrete element method （DEM） to investigate the cracking behavior of asphalt concrete （AC）. Using the ＂Fish＂ language provided in the particle flow code in 3-Demensions （PFC3D）, the air voids and mastics in asphalt concrete were realistically built as two distinct phases. With the irregular shape of individual aggregate particles modeled using a clump of spheres of different sizes, the three-dimensional （3D） discrete element model was able to account for aggregate gradation and fraction. Laboratory uniaxial complex modulus test and indirect tensile strength test were performed to obtain input material parameters for the numerical simulation. A set of the indirect tensile test were simulated to study the cracking behavior of AC at two levels of temperature, i e, -10 ℃ and 15 ℃. The predicted results of the numerical simulation were compared with laboratory experimental measurements. Results show that the 3D DEM model is able to predict accurately the fracture pattern of different asphalt mixtures. Based on the DEM model, the effects of air void content and aggregate volumetric fraction on the cracking behavior of asphalt concrete were evaluated.

THREE-DIMENSIONAL ELLIPTIC CRACK UNDER IMPACT LOADING

The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.

Interactions of penny-shaped cracks in three-dimensional solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.

Simulation of three-dimensional tension-induced cracks based on cracking potential function-incorporated extended finite element method

In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.

Analysis of Three-dimensional Crack Propagation by Using Displacement Discontinuity Method

A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.

A modified maximum tangential tensile stress criterion for three-dimensional crack propagation

The three-dimensional （3D） crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterion, it is supposed that cracks propagate only at crack front in the principal normal plane. The tangential tensile stress at crack front in the principal normal plane in local coordinates is employed to determine crack propagation, which is calculated through coordinate transformation from global to local coordinates. New cracks will propagate when the maximum tangential tensile stress at crack front in the principal normal plane reaches the tensile strength of rock-like materials. Compared with the previous crack propagation criteria, the modified crack propagation criterion is helpful in calculating 3D crack stress intensity factor, and can overcome the limitations of propagation step determined by individual experiences in previous studies. Finally, the 3D crack propagation process is traced by element-free Galerkin method. The numerical results agree well with the experimental ones for a frozen resin sample with prefabricated 3D cracks.

INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.

A STUDY ON STRESS INTENSITY FACTORS AND SINGULAR STRESS FIELDS OF THREE-DIMENSIONAL INTERFACE CRACK

By using the finite-part integral concepts and limit technique,the hypersingular inte- grodifferential equations ofthree-dimensional(3D)planar interface crack were obtained; then thedominant-part analysis of 2D hypersingular integral was further usedto investigate the stress fields near the crack front theoretically,and the accurate formulae were obtained for the singular stressfields and the complex stress intensity factors.

CONVOLUTION OF THE IMPACT THREE-DIMENSIONAL ELASTO-DYNAMICS AND DYNAMIC STRESS INTENSITY FACTOR FOR AN ELLIPTIC CRACK

This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method.

Three-dimensional Monte Carlo simulation of bulk fin field effect transistor

In this paper, we investigate the performance of the bulk fin field effect transistor （FinFET） through a three- dimensional （3D） full band Monte Carlo simulator with quantum correction. Several scattering mechanisms, such as the acoustic and optical phonon scattering, the ionized impurity scattering, the impact ionization scattering and the surface roughness scattering are considered in our simulator. The effects of the substrate bias and the surface roughness scattering near the Si/SiO2 interface on the performance of bulk FinFET are mainly discussed in our work. Our results show that the on-current of bulk FinFET is sensitive to the surface roughness and that we can reduce the substrate leakage current by modulating the substrate bias voltage.

Study on Crack Propagation Parameters of Tunnel Lining Structure Based on Peridynamics

The numerical simulation results utilizing the Peridynamics(PD)method reveal that the initial crack and crack propagation of the tunnel concrete lining structure agree with the experimental data compared to the Japanese prototype lining test.The load structure model takes into account the cracking process and distribution of the lining segment under the influence of local bias pressure and lining thickness.In addition,the influence of preset cracks and lining section formon the crack propagation of the concrete lining model is studied.This study evaluates the stability and sustainability of tunnel structure by the Peridynamics method,which provides a reference for the analysis of the causes of lining cracks,and also lays a foundation for the prevention,reinforcement and repair of tunnel lining cracks.

Three dimensional K-Tz stress fields around the embedded center elliptical crack front in elastic plates

Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios （a/c） of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances （r/a） in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.

An extended J-integral for evaluating fluid-driven cracks in hydraulic fracturing

J-integral has served as a powerful tool in characterizing crack tip status. The main feature, i.e. path- independence, makes it one of the foremost fracture parameters. In order to remain the path- independence for fluid-driven cracks, J-integral is revised. In this paper, we present an extended J-in- tegral explicitly for fluid-driven cracks, e.g. hydraulically induced fractures in petroleum reservoirs, for three-dimensional （3D） problems. Particularly, point-wise 3D extended J-integral is proposed to char- acterize the state of a point along crack front. Besides, applications of the extended J-integral to porous media and thermally induced stress conditions are explored. Numerical results show that the extended J- integral is indeed path-independent, and they are in good agreement with those of equivalent domain integral under linear elastic and elastoplastic conditions. In addition, two distance-independent circular integrals in the K-dominance zone are established, which can be used to calculate the stress intensity factor （SIF）.

THE 3-D STRESS INTENSITY FACTOR FOR A HALF PLANE CRACK IN A TRANSVERSELY ISOTROPIC SOLID DUE TO THE MOTION OF LOADS ON THE CRACK FACES

Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results.

Closed form solution of stress intensity factors for cracks emanating from surface semi-spherical cavity in finite body with energy release rate method

In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors （SIFs） of cracks emanating from a surface semi-spherical cavity in a finite body is derived using the energy release rate theory. A mode of crack opening displacements of a normal slice is established, and the normal slice relevant functions are introduced. The proposed method is both effective and accurate for the problem of three-dimensional cracks emanating from a surface cavity. A series of useful results of SIFs are obtained.

Elastodynamic stress intensity factor due to three-dimensional concentrated transient loading on the faces of a crack

The dynamic stress intensity factor for a semi-infinite crack in an otherwise unbounded elastic body is analyzed The crack is subjected to a pair of suddenly applied point loads on its faces at a distance l away from the crack tip The solution of the problem is obtained by superposition of the solutions of two simpler problems. The first of these problems is Lamb' s problem, while the second problem considers a half space with its surface subjected to the negative of the normal displacement induced by Lamb's problem in the range x>0. The latter is solved by means of integral transforms together with the application of Weiner-Hopf technique and Cagniard-de Hoop method. An exact expression is derived for the mode I stress intensity factor as a function of time for any point along the crack edge. Some features of the solution are discussed.

Elastodynamic Stress Intensity Factor History for a Semi-infinite Crack Under Three-dimensional Transient Loading

The dynamic stress intensity factor history for a semi-infinite crack in an otherwise unbounded elastic body is analyzed. The crack is subjected to a pair of suddenly-applied point loadings on its faces at a distance L away from the crack tip. The exact expression for the mode I stress intensity factor as a function of time is obtained. The method of solution is based on the direct application of integral transforms, the Wiener-Hopf technique and the Cagniard-de Hoop method. Due to the existence of the characteristic length in loading, this problem was long believed a knotty problem. Some features of the solutions are discussed and graphical result for numerical computation is presented.

Three-Dimensional Stress and Stress Intensity for Tensioned Flat Plates with Edge Cracks

The stress in the thickness direction is an important factor influencing the fracture behavior of structural members. A stress σy tensioned flat plate with edge cracks is widely used as an analysis model. The stresses σx and σy for the plate model can be acquired from Neuber＇s solution. However, the solution is applicable only for a perfect plane stress or plane strain state. As a consequence of the thickness of the plate a three-dimensional （3-D） stress state will arise near the crack tip, resulting in a variation of the distribution of σx and σy stresses. A full analysis for the 3-D stress fields for a tensioned flat plate with edge cracks has been therefore carried out. The results show that the 3-D stress field near the crack tip is mainly determined by two factors： the thickness of the plate and the curvature radius at the crack tip. A further analysis has been carried out for the stress intensity near the crack tip. In this paper we give some equations matching to the 3-D stress and stress intensity, which describe precisely the stress state near the crack tip, and which can be applied effectively in engineering analysis.

An improved three-dimensional spherical DDA model for simulating rock failure

In this paper, a discontinuous numerical model, namely SDDARF3D(three-dimensional spherical discontinuous deformation analysis for rock failure), is proposed for simulating the whole process of rock failure. Firstly, within the framework of the classical discontinuous deformation analysis(DDA) method, the formulation of three-dimensional spherical DDA(3D SDDA) is deduced; secondly, a bonding and cracking algorithm is constructed and the SDDARF3 D model is proposed; thirdly, corresponding VC++ calculation code is developed and some verification examples are calculated. The simulated results can intuitively reproduce the failure phenomena of rock mass, indicating that the proposed SDDARF3 D numerical model is correct and effective.