DYNAMIC STRESS INTENSITY FACTORS AROUND TWO CRACKS NEAR AN INTERFACE OF TWO DISSIMILAR ELASTIC HALF－PLANES UNDER IN－PLANE SHEAR IMPACT LOAD

Transient stresses around two collinear cracks which lie in parallel with theinterface of the two dissimilar half-planes are studied in this article.The surfaces ofthe cracks are sheared suddenly. Application of the Fourier and Laplace transforms technique reduces the problem to that of solving dual integrai equations.To solvethese,the differences of.the crack surface displacements are expanded in a series offunctions which are automatically zero outside of the cracks. The unknown coefficients accompanied in the series are determined by the Schmidt method. The stress intensity .factors are defined in the Laplace transform domain and these are inverted numerically in the physical space .As an example ,the dynamic stress intensity factors around two cracks in a ceramic and steel bonded composite are numerically calculated.

STRESS INTENSITY FACTORS OF A PLATE WITH TWO CRACKS EMANATING FROM AN ARBITRARY HOLE

In this paper, Muskhelishvili complex function theory and boundary collocation method are used to calculate the stress intensity factors (SIF) of a plate with two cracks emanating from an arbitrary hole. The calculated examples include a circular, elliptical, rectangular, or rhombic hole in a plate. The principle and procedure by the method is not only rather simple, but also has good accuracy. The SIF values calculated compare very favorably with the existing solutions. A t the same time,the method can be used far different finite plate with two cracks emanating from a hole with more complex geometrical and loading conditions. It is an effective unified approach for this kind of fracture problems.

A CLOSED FORM SOLUTION OF STRESS INTENSITY FACTORS FOR THREE DIMENSIONAL FINITE BODIES WITH ECCENTRIC CRACKS

In this paper, a new analytical-engineering method of closed form solution about stress intensity factors for three dimensional finite bodies with eccentric cracks is derived by means of energy release rate method. The results of stress intensity factors can be obtained. The results provided ir this method are in nice agreement with those of the famous alternating method by which only special cases can be solved.

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.

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.

Crack propagation mechanism of compression-shear rock under static-dynamic loading and seepage water pressure

To reveal the water inrush mechanics of underground deep rock mass subjected to dynamic disturbance such as blasting, compression-shear rock crack initiation rule and the evolution of crack tip stress intensity factor are analyzed under static-dynamic loading and seepage water pressure on the basis of theoretical deduction and experimental research. It is shown that the major influence factors of the crack tip stress intensity factor are seepage pressure, dynamic load, static stress and crack angle. The existence of seepage water pressure aggravates propagation of branch cracks. With the seepage pressure increasing, the branch crack experiences unstable extension from stable propagation. The dynamic load in the direction of maximum main stress increases type I crack tip stress intensity factor and its influence on type II crack intensity factor is related with crack angle and material property. Crack initiation angle changes with the dynamic load. The initial crack initiation angle of type I dynamic crack fracture is 70.5°. The compression-shear crack initial strength is related to seepage pressure, confining pressure, and dynamic load. Experimental results verify that the initial crack strength increases with the confining pressure increasing, and decreases with the seepage pressure increasing.

Theoretical and numerical studies of crack initiation and propagation in rock masses under freezing pressure and far-field stress

Water-bearing rocks exposed to freezing temperature can be subjected to freezeethaw cycles leading tocrack initiation and propagation, which are the main causes of frost damage to rocks. Based on theGriffith theory of brittle fracture mechanics, the crack initiation criterion, propagation direction, andcrack length under freezing pressure and far-field stress are analyzed. Furthermore, a calculation methodis proposed for the stress intensity factor （SIF） of the crack tip under non-uniformly distributed freezingpressure. The formulae for the crack/fracture propagation direction and length of the wing crack underfreezing pressure are obtained, and the mechanism for coalescence of adjacent cracks is investigated.In addition, the necessary conditions for different coalescence modes of cracks are studied. Using thetopology theory, a new algorithm for frost crack propagation is proposed, which has the capability todefine the crack growth path and identify and update the cracked elements. A model that incorporatesmultiple cracks is built by ANSYS and then imported into FLAC3D. The SIFs are then calculated using aFISH procedure, and the growth path of the freezing cracks after several calculation steps is demonstratedusing the new algorithm. The proposed method can be applied to rocks containing fillings such asdetritus and slurry. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

SIF-based fracture criterion for interface cracks

The complex stress intensity factor K governing the stress field of an interface crack tip may be split into two parts, i.e.,■ and s^（-iε）, so that K = ■ s^（-iε）, s is a characteristic length and ε is the oscillatory index. ■ has the same dimension as the classical stress intensity factor and characterizes the interface crack tip field. That means a criterion for interface cracks may be formulated directly with■, as Irwin（ASME J. Appl. Mech. 24：361–364, 1957） did in 1957 for the classical fracture mechanics. Then, for an interface crack,it is demonstrated that the quasi Mode I and Mode II tip fields can be defined and distinguished from the coupled mode tip fields. Built upon SIF-based fracture criteria for quasi Mode I and Mode II, the stress intensity factor（SIF）-based fracture criterion for mixed mode interface cracks is proposed and validated against existing experimental results.

Characterizing Dessication Cracking of a Remolded Clay (Eutric Vertisol) Using the Fractal Dimension Approach

Fractal dimension fd, was used as one of the parameters to describe dessicationcracking pattern of a remolded Black Cotton soil (Eutric Vertisol). The fractal dimension computed from filtered, thinned and skeletonized binary images of soil cracks using the Fractal3 software provided an insight into temporal variability of fd as well as its relationship with the Crack Intensity Factor (CIF) and Soil Moisture Content (SMC). The results showed that even for single crack, the fd prior to filtering and thinning were higher than after. Cracking patterns were observedfroma chosen soil sample during dessication and the corresponding relationship between fd and CIF compared and monitored. As the critical SMC decreased during drying (45% to 27%), the CIF soil increased (0.023% - 5.75%), so did the fd (1.233 to 1.7193). The fd showed a positive linear correlation with CIF at r2 = 0.247 (P fd with SMC was best described using a polynomial function at r2 = 0.969 (P fd was sensitive to dessication cracking and therefore on SMC changes. Visual observation of dessication cracking showed that CIF increased and attained stability after day 4 while the computed and logarithmic transformed crack area attained stability between days 7 to 10 gradually decreasing to values below 2%. The estimated crack Cover or Brightness of the digitized binary images also gave better approximation of the CIF though this was slightly higher. Our results showed that dessication cracking of the Eutric Vertisol was independent of antecedent critical SMC and was time-constrained. Further soil cracking therefore stopped once maximum CIF was attained and only widening and deepening of pre-existing cracks continued.

A NEW METHOD OF ANALYZING SURFACE CRACKS

The authors have developed a new line-spring boundary element method in the present paper, which combines the advantage of the line-spring model with that of the boundary element method. This method reduces the three-dimension problem of the surface cracks into a quasi-one-dimension problem and can be used to analyze the surface cracked plate under various loading conditions. In this paper theoretical analyses and numerical verifications are carried out. The calculated results are reported, which indicate that the present method is efficient and can be used to analyze the surface crack problem on a personal computer.

Crack patterns corresponding to the residual strength plateau of ceramics subjected to thermal shock

The formation strength plateau of ceramics is addressed. A set of of 99A1203 are conducted, mechanism of the residual subjected to thermal shock thermal shock experiments where the thin specimens of 1 mm× 10 mm×50 mm exhibit parallel through edge cracks, and thus permit quantitative measurements of the crack patterns. The cracks evolve with the severity of ther- mal shock. It is found that there is a correlation between the length and density of the thermal shock cracks. The increase of crack length weakens the residual strength, whereas the increase of crack density improves it. In a considerably wide temperature range, the two contrary effects just counteract each other; consequently a plateau appears in the variation curve of the residual strength. A comparison between the numerical and experimental results of the residual strength is made, and they are found in good agreement. This work is helpful to a deep understanding of the thermal shock failure of ceramics.

ANALYSIS OF COMPLEX STRESS INTENSITIES FOR CRACKED LAMINATES

Classical plate theory has been used to find out interfacial stress intensity factors in composite laminates. By using a well-known relation between the crack-tip energy release rate and the complex stress intensity factor. a closed-form solution for complex. Stress intensity in terms of external loading and a mode mix parameter for fairly. general composite laminates is given. Then a procedure for determining this mode mix. parameter is presented. followed by numerical results for some laminates. Small scale contact condition is expressed in terms of external loading In particular, a symmetric property of interfacial toughness curye is proven. Finally. the accuracy of failure load predicled by elininating oscllation index is discussed. and an example is presented to show the validity and limitation of β=0 approximation.

Determining representative elementary volume size of in-situ expansive soils subjected to drying-wetting cycles through field test

Cracks resulting from cyclic wetting and drying of expansive soils create discontinuities and anisotropy in the soil.The representative elementary volume(REV)defined by the continuous-media theory cannot be applied to cracked expansive soils that are considered discontinuous media.In this study,direct shear tests of three different scales(30 cm^(2),900 cm^(2),1963 cm^(2))and crack image analysis were carried out on undisturbed soil samples subjected to drying-wetting cycles in-situ.The REV size of expansive soil was investigated using the crack intensity factor(CIF)and soil cohesion.The results show that soil cohesion decreased with increasing sample area,and the development of secondary cracks further exacerbated the size effect of sample on cohesion of the soil.As shrinkage cracks developed,the REV size of the soil gradually increased and plateaued after 3−5 cycles.Under the same drying-wetting cycle conditions,the REV size determined using soil cohesion(REV-C)is 1.75 to 2.97 times the REV size determined using CIF(REV-CIF).Under the influence of shrinkage cracks,the average CIF is positively correlated with the REV size determined using different maximum permissible errors,with the coefficient of correlation greater than 0.9.A method for determining the REV-C based on crack image analysis is proposed,and the REV-C of expansive soil in the study area under different exposure times is given.

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.

Calculation of stress intensity factor in two-dimensional cracks by strain energy density factor procedure

In order to calculate the stress intensity factor(SIF) of crack tips in two-dimensional cracks from the viewpoint of strain energy density, a procedure to use the strain energy density factor to calculate the SIF is proposed. In this paper, the procedure is presented to calculate the SIF of crack tips in mode I cracks, mode II cracks and I+II mixed mode cracks. Meanwhile, the results are compared to those calculated by traditional approaches or other approaches based on strain energy density and verified by theoretical solutions. Furthermore, the effect of mesh density near the crack tip is discussed, and the proper location where the strain energy density factor is calculated is also studied. The results show that the SIF calculated by this procedure is close to not only those calculated by other approaches but also the theoretical solutions, thus it is capable of achieving accurate results.Besides, the mesh density around the crack tip should meet such requirements that, in the circular area created, the first layer of singular elements should have a radius about 0.05 mm and each element has a circumferential directional meshing angle to be15°–20°. Furthermore, for a single element around the crack tip, the strain energy density factor is suggested to be calculated in the location where half of the sector element's radius from the crack tip.

ELECTRO-ELASTIC GREEN'S FUNCTIONS FOR APIEZOELECTRIC HALF-SPACE ANDTHEIR APPLICATION

In this paper, as is studied are the electro-elastic solutions for a piezoelectric halfspace subjected Io a line force, a line charge and a line dislocation, i. e.. Green sfunclions on the basis of Stroh formalism and the concept of analytical continuation,explicit expressions for Green's functions are derived. As a direct application of theresults obtained, an infinite piezoelectric solid containing a semi-infinite crack isexammed. Attention iffocused on the stress and electric displacement fields of a cracktip. The stress and electric displacement intensity .factors are given explicitly.

A simpliied fatigue assessment method for transverse illet welded joints

Under the as-welded condition the fatigue crack initiation period was considered nonexistent and Linear Elastic Fracture Mechanics（LEFM） was used to calculate fatigue strength for a range of weld geometries. Fracture mechanics assessment of welded joints requires accurate solutions for stress intensity factor（SIF）. However, the solutions for the SIF of complex welded joints are dificult to determine due to the complicated correction factors. Three methods for SIF prediction are discussed on illet welded specimens containing continuous or semi-elliptical surface cracks, including the traditional correction method Mk, the approximate correction method Kt, and the suggested additional crack size method（ac＋ae）.The new additional crack parameter ae is used to replace the stress concentration effect of weld proile Mk, which simpliies the calculation process. Experimental results are collected to support fatigue strength assessment of the additional crack size method.