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 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.

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.

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

The problem of a Griffith crack in an unbounded orthotropic functionally graded material subjected to antipole shear impact was studied. The shear moduli in two directions of the functionally graded material were assumed to vary proportionately as definite gradient. By using integral transforms and dual integral equations, the local dynamic stress field was obtained. The results of dynamic stress intensity factor show that increasing shear moduli's gradient of FGM or increasing the shear modulus in direction perpendicular to crack surface can restrain the magnitude of dynamic stress intensity factor.

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 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.

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.

Propagation and Coalescence of Blast-Induced Cracks in PMMA Material Containing an Empty Circular Hole Under Delayed Ignition Blasting Load by Using the Dynamic Caustic Method

In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under delayed ignition blasting loads.The following experimental results are obtained.(1)In directional-fracture-controlled blasting,the dynamic stress intensity factors(DSIFs)and the propagation paths of the blast-induced cracks are obviously influenced by the delayed ignition.(2) The circular hole situated between the two boreholes poses a strong guiding effect on the coelesence of the cracks,causing them to propagate towards each other when cracks are reaching the circular hole area.(3)Blast-induced cracks are not initiated preferentially because of the superimposed effect from the explosive stress waves on the cracking area.(4) By using the scanning electron microscopy(SEM)method,it is verified that the roughness of crack surfaces changes along the crack propagation paths.

Prediction of Three Dimensional Crack Path Brittle Fracture in Weldment under Dynamic Load

Three dimensional dynamic stress intensity factors are analyzed for a curved crack with a second order perturbation method. The method is extended to obtain an approximate representation of a three dimensional dynamic stress intensity factors at the tip of a curved crack. Due to three dimensional curved crack growth the dynamic energy release rate can be calculated by using the Irwin＇s formula. A three dimensional curved crack in materials with inhomogeneous fracture toughness are considered. Paths of a brittle three dimensional curved crack propagating along a welded joint are predicted via the present method, where the effects of dynamic applied stresses, residual stresses, and material deterioration due to welding are taken into considerations.

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.

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.

A HALF PLANE CRACK UNDER THREE-DIMENSIONAL COMBINED MODE IMPACT LOADING

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 suddenly-applied shear line loads is consid- ered. The analytic expression for the combined mode stress intensity factors as a function of time is obtained. The method of solution is based on the application of integral transforms and the Wiener-Hopf technique. Some features of the solutions are discussed and graphical numerical results are presented.

SEMI-ELLIPTIC SURFACE CRACK IN AN ELASTIC SOLID WITH FINITE SIZE UNDER IMPACT LOADING

In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a three-dimensional finite element program is developed to compute the dynamic stress intensity factor. The results reveal that the effects of the solid＇s boundary surface, crack surface, material inertia and stress wave interactions play significant roles in dynamic fracture.