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.

Solution of stress intensity factors of multiple cracks in plane elasticity with eigen COD formulation of boundary integral equation

The concept of eigen crack opening displacement （COD） can be defined as the COD of a crack in infinite plate under the tractions acting on the crack surface. By introducing this concept, the eigen COD formulation of boundary integral equation is proposed in this paper, together with the solution procedures for multiple crack problems in plane elasticity. With the proposed approach, the multiple crack problems can be solved with the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix as that in the numerical Green’s function （NGF） approach but without the trouble to determine the complementary solutions since the standard boundary element discretization on the crack surface is no longer required with the proposed approach. Some numerical examples computing the stress intensity factors are presented and compared with those in literature to show the accuracy and the effectiveness of the proposed approach.

A WEIGHTED RESIDUAL METHOD FOR ELASTIC-PLASTIC ANALYSIS NEAR A CRACK TIP AND THE CALCULATION OF THE PLASTIC STRESS INTENSITY FACTORS

In this paper, a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample. The elastic-plastic solutions of the crack lip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field, to calculate the plastic stress intensity factors, are also developed. Therefore, a complete analvsis based on the calculation both for the crack tip field and for the whole crack body field is provided.

METHOD TO CALCULATE BENDING CENTER AND STRESS INTENSITY FACTORS OF CRACKED CYLINDER UNDER SAINT_VENANT BENDING

Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross_section is not thin_walled, but of small torsion rigidity is proposed. Some numerical examples are given.

Evaluation of stress intensity factors for bi-material interface cracks using displacement jump methods

The aim of the present work is to investigate the numerical modeling of interfacial cracks that may appear at the interface between two isotropic elastic materials. The extended finite element method is employed to analyze brittle and bi-material interfacial fatigue crack growth by computing the mixed mode stress intensity factors(SIF). Three different approaches are introduced to compute the SIFs. In the first one, mixed mode SIF is deduced from the computation of the contour integral as per the classical J-integral method,whereas a displacement method is used to evaluate the SIF by using either one or two displacement jumps located along the crack path in the second and third approaches. The displacement jump method is rather classical for mono-materials,but has to our knowledge not been used up to now for a bimaterial. Hence, use of displacement jump for characterizing bi-material cracks constitutes the main contribution of the present study. Several benchmark tests including parametric studies are performed to show the effectiveness of these computational methodologies for SIF considering static and fatigue problems of bi-material structures. It is found that results based on the displacement jump methods are in a very good agreement with those of exact solutions, such as for the J-integral method, but with a larger domain of applicability and a better numerical efficiency(less time consuming and less spurious boundary effect).

THE EVALUATION OF STRESS INTENSITY FACTORS OF PLANE CRACK FOR ORTHOTROPIC PLATE WITH EQUAL PARAMETER BY F2LFEM

In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method （F2LFEM）. The general solution of an orthotropic crack problem is obtained by assimilating the problem with isotropic crack problem, and is employed as the global interpolation function in F2LFEM. In the neighborhood of crack tip of the crack plate, the fractal geometry concept is introduced to achieve the similar meshes having similarity ratio less than one and generate an infinitesimal mesh so that the relationship between the stiffness matrices of two adjacent layers is equal. A large number of degrees of freedom around the crack tip are transformed to a small set of generalized coordinates. Numerical examples show that this method is efficient and accurate in evaluating the stress intensity factor （SIF）.

STRESS INTENSITY FACTORS FOR A FINITE PLATE WITH AN INCLINED CRACK BY BOUNDARY COLLOCATION

In this paper, we combine the Muskhelishvili＇s complex variable method and boundary collocation method, and choose a set of new stress function based on the stress boundary condition of crack surface, the higher precision and less computation are reached. This method is applied to calculating the stress intensity factor for a finite plate with an inclined crack. The influence of θ （the obliquity of crack） on the stress intensity factors, as well as the number of summation terms on the stress intensity factor are studied and graphically represented.

Simulation Research on Stress Intensity Factors of Different Crack Aspect Ratios on Hollow Axles

Because of the wicked service environment of the high speed train, it is possible that the hollow axle of the train may encounter the foreign object damage and form a sharp notch. Under the fatigue loading a crack can initiate from the notch and propagate to failure. It is noted that the stress intensity factor is the control parameter of the crack propagating, for the purpose of getting the more exact propagation characteristics, the stress intensity factor is studied mainly. The service loads of hollow axles are defined, and the stress distribution of hollow axles is obtained according to the load spectrum. The semi-ellipse crack configuration is defined with three parameters： the aspect ratio, the relative depth and the relative location along the crack front. Quarter point 20-node isoparametric degenerate singular elements are used for the region near the crack tip. The finite element model of crack extension of hollow axle is created, and the crack front is dispersed which can realize orthogonal extension. Based on this the stress intensity factors of crack front were calculated, and the distribution rules of the stress intensity factors of different initial crack shapes are obtained. The conclusions are compared with that of the analytic method and they agree with each other very well, and the calculating results show that there is a close relationship between the stress intensity factor and the initial crack shape. For a round crack the stress intensity factor at the surface point increases faster than the one at the center point with the crock propagation. However, for a narrow crack, the results are in contrast with that of a round one. So, all the cracks with different shapes propagate toward to a similar shape, and they grow at this shape to end. The study may contribute to the crack propagate characteristics research.

Compensation of stress intensity factors in hollow cylinders containing several cracks under torsion by electro-elastic coating

In this article, a formulation for a hollow cylinder reinforced with an electroelastic layer is investigated. The hollow cylinder and its electro-elastic coating are under the Saint-Venant torsional loading. First, the solution to the problem containing a Volterra-type screw dislocation is obtained by using the Fourier transform. The problem is then reduced to a set of Cauchy singular integral equations by the distributed dislocation method. Finally, several examples are presented to show the effect of the electro-elastic coating on the reduction of the stress intensity factors at the crack tips.

Analysis of stress intensity factor in orthotropic bi-material mixed interface crack

Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With con- sideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial dif- ferential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit, undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship be- tween the bi-material stress intensity factor and the ratio of bi-material Young＇s modulus are given in the numerical analysis.

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.

Comprehensive investigation of stress intensity factors in rotating disks containing three-dimensional semi-elliptical cracks

Initiation and propagation of cracks in rotating disks may cause catastrophic failures. Therefore, determination of fracture parameters under different working con- ditions is an essential issue. In this paper, a comprehensive study of stress intensity factors （SIFs） in rotating disks containing three-dimensional （3D） semi-elliptical cracks subjected to different working conditions is carried out. The effects of mechanical prop- erties, rotational velocity, and orientation of cracks on SIFs in rotating disks under cen- trifugal loading are investigated. Also, the effects of using composite patches to reduce SIFs in rotating disks are studied. The effects of patching design variables such as mechanical properties, thickness, and ply angle are investigated separately. The modeling and analytical procedure are verified in comparison with previously reported results in the literature.

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.

Finite element simulation of stress intensity factors in elastic-plastic crack growth

A finite element program developed elastic-plastic crack propagation simulation using Fortran language. At each propagation step, the adaptive mesh is automatically refined based on a posteriori h-type refinement using norm stress error estimator. A rosette of quarter-point elements is then constructed around the crack tip to facilitate the prediction of crack growth based on the maximum normal stress criterion and to calculate stress intensity factors under plane stress and plane strain conditions. Crack was modelled to propagate through the inter-element in the mesh. Some examples are presented to show the results of the implementation.

SCATTERING OF SH-WAVE BY CRACKS ORIGINATING AT AN ELLIPTIC HOLE AND DYNAMIC STRESS INTENSITY FACTOR

The method of complex function and the method of Green's function are used to investigate the problem of SH-wave scattering by radial cracks of any limited length along the radius originating at the boundary of an elliptical hole, and the solution of dynamic stress intensity factor at the crack tip was given. A Green's function was constructed for the problem, which is a basic solution of displacement field for an elastic half space containing a half elliptical gap impacted by anti-plane harmonic linear source force at any point of its horizontal boundary. With division of a crack technique, a series of integral equations can be established on the conditions of continuity and the solution of dynamic stress intensity factor can be obtained. The influence of an elliptical hole on the dynamic stress intensity factor at the crack tip was discussed.

METHOD TO CALCULATE BENDING CENTER AND STRESS INTENSITY FACTORS OF CRACKED CYLINDER UNDER SAINT-VENANT BENDING

Using the single crack solution and the regular solution elf plane harmonic function, the problem of Saint-Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross-section is not thin-walled, but of small torsion rigidity is proposed. Some numerical examples are given.

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.

Determination of stress intensity factors in half-plane containing several moving cracks

The dynamic stress intensity factors in a half-plane weakened by several finite moving cracks are investigated by employing the Fourier complex transformation. Stress analysis is performed in a half-plane containing a single dislocation and without dislocation. An exact solution in a closed form to the stress fields and displacement is ob- tained. The Galilean transformation is used to transform between coordinates connected to the cracks. The stress components are of the Cauchy singular kind at the location of dislocation and the point of application of the the influence of crack length and crack running force. Numerical examples demonstrate velocity on the stress intensity factor.

Dynamic stress intensity factor K_(III) and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.