Effect of Blasting Stress Wave on Dynamic Crack Propagation

Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical characteristics of crack propagation is of great significance for engineering blasting.In this study,ANSYS/LS-DYNA was used for blasting numerical simulation,in which the propagation characteristics of blasting stress waves and stress field distribution at the crack tip were closely observed.Moreover,ABAQUS was applied for simulating the crack propagation path and calculating dynamic stress intensity factors(DSIFs).The universal function was calculated by the fractalmethod.The results show that:the compressive wave causes the crack to close and the reflected tensile wave drives the crack to initiate and propagate,and failure mode is mainly tensile failure.The crack propagation velocity varies with time,which increases at first and then decreases,and the crack arrest occurs due to the attenuation of stress waves and dissipation of the blasting energy.In addition,crack arrest toughness is smaller than the crack initiation toughness,applied pressure waveforms(such as the peak pressure,duration,waveforms,wavelengths and loading rates)have a great influence on DSIFs.It is conducive to our deep understanding or the study of blasting stress waves dominated fracture,suggesting a broad reference for the further development of rock blasting in engineering practice.

Dynamic Stress Intensity Factor and Scattering of SH Wave by Circle Canyon and Crack

The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed first, which is the solution of displacement fields for elastic half space with circle canyon under output plane harmonic line loading at horizontal surface. Then the integral equation for determining the unknown forces in the problem can be changed into the algebraic one and solved numerically so that crack DSIF can be determined. Last when the medium parameters are altered, the influence on the crack DSIF is discussed partially with the displacement between circle canyon and crack.

The dynamic stress intensity factor analysis of adhesively bonded material interface crack with damage under shear loading

This paper studies the dynamic stress intensity factor （DSIF） at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the integral transform, the problem is reduced to algebraic equations and can be solved with the collocation dots method in the Laplace domain. Time response of DSIF is calculated with the inverse Laplace integral transform. The results show that the mode Ⅱ DSIF increases with the shear relaxation parameter, shear module and Poisson ratio, while decreases with the swell relaxation parameter. Damage shielding only occurs at the initial stage of crack propagation. The singular index of crack tip is -0.5 and independent on the material parameters, damage conditions of materials, and time. The oscillatory index is controlled by viscoelastic material parameters.

STUDY ON DYNAMIC STRESS INTENSITY FACTORS OF DISK WITH A RADIAL EDGE CRACK SUBJECTED TO EXTERNAL IMPULSIVE PRESSURE

A dynamic weight function method is presented for dynamic stress intensity factors of circular disk with a radial edge crack under external impulsive pressure. The dynamic stresses in a circular disk are solved under abrupt step external pressure using the eigenfunction method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary conditions. By making use of Fourier- Bessel series expansion, the history and distribution of dynamic stresses in the circular disk are derived. Furthermore, the equation for stress intensity factors under uniform pressure is used as the reference case, the weight function equation for the circular disk containing an edge crack is worked out, and the dynamic stress intensity factor equation for the circular disk containing a radial edge crack can be given. The results indicate that the stress intensity factors under sudden step external pressure vary periodically with time, and the ratio of the maximum value of dynamic stress intensity factors to the corresponding static value is about 2.0.

THE ANALYSIS OF DYNAMIC STRESS INTENSITY FACTOR FOR SEMI-CIRCULAR SURFACE CRACK USING TIME-DOMAIN BEM FORMULATION

The time-domain BEM was developed to analyze the dynamic stress intensity factor ( DSIF) of 3-D elastodynamic crack problems. To simulate the stress singularity along the front of a crack, eight-node isoparametric singular elements were used, and the DSIF for a semi-circular surface crack was firstly calculated based on displacement equation using the time-domain BEM formulation. The new scheme to determine the time step was brought forward. By the dynamic analysis program of time-domain BEM compiled by its, several numerical examples are presented, which demonstrate the unconditional stability and high accuracy of time-domain BEM applied to 3-D elastodynamic crack problems.

Analysis of dynamic stress intensity factors of thick-walled cylinder under internal impulsive pressure

Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is used as the reference case, and then the weight function for a thick-walled cylinder containing a radial edge crack can be worked out. Secondly, the dynamic stresses in uncracked thick-walled cylinders are solved under internal impulsive pressure by using mode shape function method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary condi- tions, and the history and distribution of dynamic stresses in thick-walled cylinders are derived in terms of Fourier-Bessel series. Finally, the dynamic stress intensity factor equations for thick-walled cylinder containing a radial edge crack sub- jected to internal impulsive pressure are given by dynamic weight function method. The finite element method is utilized to verify the results of numerical examples, showing the validity and feasibility of the proposed method.

Analysis of dynamic stress intensity factors of three-point bend specimen containing crack

A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally~ the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method （FEM）. Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.

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.

Effect of crack face contact on dynamic stress intensity factors for a hole-edge crack

In order to determine the dynamic stress intensity factors(DSIFs)for a single edge crack at the center hole of a finite plate under a compressive step loading parallel to the crack,the finite element method was employed to solve the cracked plate problem.The square-root stress singularity around the crack tip was simulated by quarter point singular elements collapsed by 8-node two-dimensional isoparametric elements.The DSIFs with and without considering crack face contact situations were evaluated by using the displacement correlation technique,and the influence of contact interaction between crack surfaces on DSIFs was investigated.The numerical results show that if the contact interaction between crack surfaces is ignored,the negative mode I DSIFs may be obtained and a physically impossible interpenetration or overlap of the crack surfaces will occur.Thus the crack face contact has a significant influence on the mode I DSIFs.

THE THREE-DIMENSIONAL STRESS INTENSITY FACTOR UNDER MOVING LOADS ON THE CRACK FACES

The dynamic stress intensity factor history for a half plane crack in an otherwise unbounded elastic body,with the crack faces subjected to a traction distribution consisting of two pairs of combined mode point loads that move in a direction perpendicular to the crack edge is considered.The analytic expression for the combined mode stress intensity factors as a function of time for any point along the crack edge is obtained.The method of solution is based on the application of integral transform together with the Wiener-Hopf technique and the Cagniard-de Hoop method. Some features of the solution are discussed and graphical results for various point load speeds are presented.

Experimental study on dynamic fracture behavior of three-point bending beam with double deformity inclusions

The dynamic fracture behavior of the three-point bending beam with double deformity inclusions under impact loading is studied by using digital high-speed photography in combination with the transmission-type dynamic caustic method. The experimental results indicate that the fluctuation of crack propagation velocity v first increases and then decreases in the crack propagation process. During the process of crack propagating into the inclusion area, the fracture resistance effect of the circular inclusion is the most significant and the effects of triangular and square inclusions are less obvious. The stress intensity factor near the crack tip increases during the propagation process and reaches its maximum value when the crack tip is close to the inclusions. The crack tip’s dynamic stress intensity factor （ DSIF） decreases when the crack exceeds the middle area of the double inclusions. These results provide an experimental basis and scientific foundation to strengthen the evaluation and fracture analysis of the structure containing deformity inclusions.

Dynamic caustics experimental study on interaction between propagating crack and deformity inclusions in primary structure

The approach combining the dynamic caustics method with high-speed photography technology is used to study the interaction between propagating cracks and three kinds of deformity inclusions（ cylinder inclusion, quadruple inclusion and triangular inclusion） under lowvelocity impact loading. By recording the caustic spots of crack tips at different moments during the crack propagation, the variation regulations of dynamic stress intensity factors（ DSIF） and crack growth velocity with respect to time are obtained. The experimental results showthat the resistance effects to crack growth are varied with different shapes of inclusions in specimens, and the quadruple inclusion＇s effect is more apparent. The distortion degree of caustic spots is affected by the shapes of inclusions as well, and the situation is more serious for cylinder and quadruple inclusions. The overall values of DSIFs of triangular inclusion specimen are greater than the others, and the crack growth velocities, characteristic sizes and DSIFs showprocesses of fluctuations because of the disturbance of reflection waves in specimens. The results provide an experimental basis for the analysis of strength and impact-resistance ability in structures with deformity inclusions.

Caustics study of the effect of glass fibres on dynamic fracture of hardened cement paste

The reflected optical caustics method is applied to study dynamic fracture problems in hardened cement paste. First both the unreinforced cement paste and the glass fibres reinforced cement paste specimens were fabricated and the reflective coating on the surface of the specimen was prepared. Secondly the crack path and the shadow spot patterns during the crack propagation process for the two specimens were recorded by using a multi-spark high speed camera.Thirdly some dynamic parameters of two cement paste specimens including crack onset time the dynamic stress intensity factor and crack growth velocity were determined and analyzed comparatively.This indicates that the glass fibres can improve the fracture resistance and delay fracture time.These results will play an important role in evaluating the dynamic fracture properties of cement paste.

Dynamic caustics test of blast load impact on neighboring different cross-section roadways

Using digital laser dynamic caustics experimental system and conducting simulation experiment researched the influence rule of blasting excavation of a new roadway on neighboring existed different cross-section roadways. The experimental results show that the influence of blast load on adjacent roadway has a good relationship with the cross-section of roadway. The expansion distance of precrack existed in circular, arch-wall, rectangular roadway is respectively 1.76, 1.61 and 0 cm under blast load.At the same time, the direct-blast side of rectangular roadway has more obvious damage compared with circular and arch-wall roadway. It explains that plane reflects more stress wave than arc, so that it exerts more tensile failure in the direct-blast side, which leads to less stress wave diffracting to the precrack in the back-blast side. When the precrack extends, higher value dynamic stress intensity factor in circular roadway works longer than that of arch-wall roadway. Indirectly, it explains that plane's weakening function on stress wave is significantly stronger than arc. Stress wave brings about self-evident influence on the upper and bottom endpoints of the rectangular roadway, and it respectively extends 1.03, 2.06 cm along the line link direction of the center of the blasthole and the upper and bottom endpoints on the right wall.

Dynamic fracture behavior of rock under impact load using the caustics method

We studied the dynamic fracture mechanical behavior of rock under different impact rates. The fracture experiment was a three-point bending beam subjected to different impact loads monitored using the reflected caustics method. The mechanical parameters for fracture of the three-poim bending beam specimen under impact load are analyzed. The mechanism of crack propagation is discussed. Experimental results show that the dynamic stress intensity factor increases before crack initiation. When the dynamic stress intensity factor reaches its maximum value the crack starts to develop. After crack initiation the dynamic stress intensity factor decreases rapidly and oscillates. As the impact rate increases the cracks initiate earlier, the maximum value of crack growth velocity becomes smaller and the values of dynamic stress intensity factor also vary less during crack propagation. The results provide a theoretical basis for the study of rock dynamic fracture.

The transient fracture behavior for a functionally graded layered structure subjected to an in-plane impact load

The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally gradedinterlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors （DSIFs） in time domain are studied and some meaningful conclusions are obtained.

SHEAR WAVE SCATTERING FROM A PARTIALLY DEBONDED PIEZOELECTRIC CYLINDRICAL INCLUSION

The scattering of SH wave by a cylindrical piezoelectric inclusion partially debonded from its surrounding piezoelectric material is investigated using the wave function expansion method and singular integral equation technique. The debonding regions are modeled as mul- tiple arc-shaped interface cracks with non-contacting faces. By expressing the scattered ?elds as wave function expansions with unknown coe?cients, the mixed boundary value problem is ?rstly reduced to a set of simultaneous dual series equations. Then dislocation density functions are introduced as unknowns to transform these dual series equations into Cauchy singular integral equations of the ?rst type, which can be numerically solved easily. The solution is valid for arbi- trary number and size of the debonds. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond and two debonds. The e?ects of incidence direc- tion, crack con?guration and various material parameters on the dynamic stress intensity factors are respectively discussed. The solution of this problem is expected to ?nd applications in the investigation of dynamic fracture properties of piezoelectric materials with cracks.

Dynamic anti-plane analysis for two symmetrically interfacial cracks near circular cavity in piezoelectric bi-materials

The present paper is exposed theoretically to the influence on the dynamic stress intensity factor （DSIF） in the piezoelectric bi-materials model with two symmet- rically permeable interracial cracks near the edges of a circular cavity, subjected to the dynamic incident anti-plane shearing wave （SH-wave）. An available theoretical method to dynamic analysis in the related research field is provided. The formulations are based on Green＇s function method. The DSIFs at the inner and outer tips of the left crack are obtained by solving the boundary value problems with the conjunction and crack- simulation technique. The numerical results are obtained by the FORTRAN language program and plotted to show the influence of the variations of the physical parameters, the structural geometry, and the wave frequencies of incident wave on the dimensionless DSIFs. Comparisons with previous work and between the inner and outer tips are con- cluded.

Effects of electric/magnetic impact on the transient fracture of interface crack in piezoelectric-piezomagnetic sandwich structure: anti-plane case

Due to the incompatibility of the interlaminar deformations,the interface debonding or cracking usually happens in a layered magnetoelectric(ME)structure under an applied load.In this paper,the transient responses of the anti-plane interface cracks in piezoelectric(PE)-piezomagnetic(PM)sandwich structures are studied by the standard methods of the integral transform and singular integral equation.Discussion on the numerical examples indicates that the PE-PM-PE structure under electric impact is more likely to fracture than the PM-PE-PM structure under a magnetic impact.The dynamic stress intensity factors(DSIFs)are more sensitive to the variation of the active layer thickness.The effects of the material constants on the DSIFs are dependent on the roles played by PE and PM media during the deformation process.

Normal compression wave scattering by a permeable crack in a fluid-saturated poroelastic solid

A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In particular, the effect of the wave-induced fluid flow due to the presence of a liquid-saturated crack on the dynamic stress intensity factor is analyzed. Fourier sine and cosine integral transforms in conjunction with Helmholtz potential theory are used to formulate the mixed boundary-value problem as dual integral equations in the frequency domain. The dual integral equations are reduced to a Fredholm integral equation of the second kind. It is found that the stress intensity factor monotonically decreases with increasing frequency, decreasing the fastest when the crack width and the slow wave wavelength are of the same order. The characteristic frequency at which the stress intensity factor decays the fastest shifts to higher frequency values when the crack width decreases.