A DISCUSSION ABOUT THE EFFECTIVENESS ON USING A VARIATE-TRANSFORMATION TO FIND OUT SOLUTIONS OF MODE Ⅲ CRACK PROBLEMS IN POWER HARDENING MEDIA

By the use of the transformations of physical plane to strain plane and physical planeto stress plane,an analytic expression of the asymptotic solution near a modeⅢcrack tipin a power hardening medium can be obtained.In this paper the effectiveness of thetransformation is discussed.Analytical results show that the transformation is effectiveexcept for a special limit case of power hardening media—the ideal plastic materials.

Boundary Element Analysis forModeⅢCrack Problems of Thin-Walled Structures from Micro-to Nano-Scales

This paper develops a new numerical framework for modeⅢcrack problems of thin-walled structures by integrating multiple advanced techniques in the boundary element literature.The details of special crack-tip elements for displacement and stress are derived.An exponential transformation technique is introduced to accurately calculate the nearly singular integral,which is the key task of the boundary element simulation of thin-walled structures.Three numerical experiments with different types of cracks are provided to verify the performance of the present numerical framework.Numerical results demonstrate that the present scheme is valid for modeⅢcrack problems of thin-walled structures with the thickness-to-length ratio in the microscale,even nanoscale,regime.

Field structure at mode Ⅲ dynamically propagating crack tip in elastic-viscoplastic materials

An elastic-viscoplastic mechanics model is used to investigate asymptotically the mode Ⅲ dynamically propagating crack tip field in elastic-viscoplastic materials. The stress and strain fields at the crack tip possess the same power-law singularity under a linear-hardening condition. The singularity exponent is uniquely determined by the viscosity coefficient of the material. Numerical results indicate that the motion parameter of the crack propagating speed has little effect on the zone structure at the crack tip. The hardening coefficient dominates the structure of the crack-tip field. However, the secondary plastic zone has little influence on the field. The viscosity of the material dominates the strength of stress and strain fields at the crack tip while it does have certain influence on the crack-tip field structure. The dynamic crack-tip field degenerates into the relevant quasi-static solution when the crack moving speed is zero. The corresponding perfectly-plastic solution is recovered from the linear-hardening solution when the hardening coefficient becomes zero.

A PENNY-SHAPED CRACK IN A TRANSVERSELY ISOTROPIC PIEZOELECTRIC SOLID:MODES Ⅱ AND Ⅲ PROBLEMS

An exact analysis of the modes Ⅱ and Ⅲ problems of a penny- shaped crack in a transversely isotropic piezoelectric medium is performed in this paper.The potential theory method is employed based on the general solution of three-dimensional piezoelasticity and the four harmonics involved are represented by one complex potential.Previous results in potential theory are then utilized to obtain the exact solution that is expressed in terms of elementary functions.Comparison is made between the current results with those published and good agreement is obtained.

General forms of elastic-plastic matching equations for mode-Ⅲ cracks near crack line

Crack line analysis is an effective way to solve elastic-plastic crack problems. Application of the method does not need the traditional small-scale yielding conditions and can obtain sufficiently accurate solutions near the crack line. To address mode- Ⅲ crack problems under the perfect elastic-plastic condition, matching procedures of the crack line analysis method axe summarized and refined to give general forms and formulation steps of plastic field, elastic-plastic boundary, and elastic-plastic matching equations near the crack line. The research unifies mode-III crack problems under different conditions into a problem of determining four integral constants with four matching equations. An example is given to verify correctness, conciseness, and generality of the procedure.

Dislocation distribution model of mode Ⅲ propagation crack

Dislocation distribution functions of the edges of mode Ⅲ propagation crack subjected to three types of loads were studied by the methods of the theory of complex variable functions,by which,the problems researched were facilely transformed into Riemann-Hilbert problems and Keldysh-Sedov mixed boundary value problems. Analytical solutions of stresses,displacement and dynamic stress intensity factor were obtained by the measures of the theory of self-similar functions and corresponding differential and integral operation. In terms of the relationship between dislocation distribution functions and displacements,analytical solutions of dislocation distribution functions were obttained,and variational rules of dislocation distribution functions were depicted.

Detection of internal crack growth in polyethylene pipe using guided wave ultrasonic testing

Despite the success of guided wave ultrasonic inspection for internal defect detection in steel pipes,its application on polyethylene(PE)pipe remains relatively unexplored.The growth of internal cracks in PE pipe severely affects its pressure-holding capacity,hence the early detection of internal cracks is crucial for effective pipeline maintenance strategies.This study extends the scope of guided wave-based ultrasonic testing to detect the growth of internal cracks in a natural gas distribution PE pipe.Laboratory experiments and a finite element model were planned to study the wave-crack interaction at different stages of axially oriented internal crack growth with a piezoceramic transducer-based setup arranged in a pitch-catch configuration.Mode dispersion analysis supplemented with preliminary experiments was performed to isolate the optimal inspection frequency,leading to the selection of the T(0,1)mode at 50-kHz for the investigation.A transmission index based on the energy of the T(0,1)mode was developed to trace the extent of simulated crack growth.The findings revealed an inverse linear correlation between the transmission index and the crack depth for crack growth beyond 20%crack depth.

Dynamic propagation problems concerning surfaces of asymmetrical mode III crack subjected to moving loads

With the theory of complex functions, dynamic propagation problems concerning surfaces of asymmetrical mode Ⅲ crack subjected to moving loads are investigated. General representations of analytical solutions are obtained with self-similar functions. The problems can be easily converted into Riemann-Hilbert problems using this technique. Analytical solutions to stress, displacement and dynamic stress intensity factor under constant and unit-step moving loads on the surfaces of asymmetrical extension crack, respectively, are obtained. By applying these solutions, together with the superposition principle, solutions of discretionarily intricate problems can be found.

In-plane shear(ModeⅡ) crack sub-critical propagation of rock at high temperature

In-plane shear crack sub-critical propagation of rock at high temperature was studied by finite element method and shear-box(i.e.compression-shear) test with newly designed electrically conductive adhesive method.Numerical and experimental results show that the normalized shear(Mode Ⅱ) stress intensity factors,K ⅡT/KT0 is decreased as the temperature increases because high temperature can improve stress distribution at crack tip and reduce the Mode Ⅱ stress intensity factor.Microscopic features of fractured surface are of little pits and secondary micro-cracks in the vicinity(1.5-4.0 mm) of the crack tip.The chevron-shape secondary cracks gradually merge in the length of about 4-5 mm and disappear along the direction of crack propagation.Stable shear crack propagation time is increased with the increasing temperature while the stable shear crack propagation rate is decreased with the increasing temperature,since high temperature can increase the shear(Mode Ⅱ) fracture toughness and prevent the crack growth.It is necessary to ensure the ligament of specimen long enough to measure the maximum unstable crack propagation rate of rock.

Effect of holes on in-plane shear(Mode Ⅱ) crack sub-critical propagation of rock

Shear-box(i.e.compression-shear) test and newly designed electrically conductive adhesive method were used to measure shear crack sub-critical propagation time and rate of sandstone specimen.Different cubic specimens with and without holes were tested to study the effect of holes on the shear crack sub-critical propagation.Numerical and experimental results show that three independent variables of hole,the interval distance S,the distance between the center of hole and the crack tip L,and hole radius R,have different contribution to the ratio of stress intensity factor of the specimen with holes to that of the specimen without hole,KⅡ/KⅡ0.Increasing S and decreasing L and R will result in the decrease of KⅡ/KⅡ0 and help crack arrest.The weight relation of the independent variables for KⅡ/KⅡ0 is S>L>R.The specimen DH3 with the largest value of S and the smallest values of L and R has the longest sub-critical crack propagation time and the smallest sub-critical crack propagation rate.Adding two suitable holes symmetrically to the original crack plane in rock specimen is considered to be a potential method for crack arrest of rock.

EXACT SOLUTIONS OF NEAR CRACK LINE FIELDS FOR MODE I CRACK UNDER PLANE STRESS CONDITION IN AN ELASTIC－PERFECTLY PLASTIC SOLID

The near crack line analysis method has been used in the present paper,The classical small scale yielding conditions have been completely abandoned in the analyses and one inappropriate matching condition used to be used at the elasticplastic boundary has been corrected.The reasonable solution of the plastic stresses near the crack line region has been established.By matching the plastic stresses with the exact elastic stresses at the elastic-plastic boundary,the plastic stresses the length of the plastic zone and the unit normal vector of the elastic-plastic boundary near the crock line region have been obtained for a mode I crack under uniaxial tension,as well as a mode I crack under biaxial tension,which shows that for both conditions the plastic stress componentsσy, and σsy.he length of the plastic zone and the unit normal vector of the elastic-plastic boundary are quite the same while the plastic stress σs is different.

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.

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.

A THEORETICAL MODEL OF MATRIX CRACKING IN COMPOSITE LAMINATES UNDER THERMOMECHANICAL LOADING

A theoretical approach is presented for analyzing the ply crackingin general symmetric lami- nates subjected to any combination ofin-plane mechanical loading and uniform temperature changes. Theequivalent constraint model proposed by the authors in a previouswork is used to account for the cracking in- teraction betweenlaminae in the laminates. By using a superposition schemce and thestress field solutions the energy release rate for a ply cracking isexplicitly as a function of stiffness reduction parameters of thelaminates. The ratio of mode Ⅰ to mode Ⅱ is introduced formconstruction of the fracture criterion. The effects of the laminateparameters and the crack spacing on the energy release rate and themode mixity are illustrated. Finally, the model is used to predictthe thermomechanical load for the first-ply-cracking.

Rock fracture under anti-plane shear (Mode Ⅲ) loading

Anti-plane punch-through shear test and anti-planefour-point bending test are used to study the crack initiation and propagation under anti-plane shear (Mode Ⅲ) loading. The tensile and shear stresses at the crack tip are calculated by finite element method. The results show that under Mode Ⅲ loading the maximum principal stress σ1 at crack tip is smaller or a little larger than the maximum shear stress τmax. Since the tensile strength of brittle rock is much lower than its shear strength, σ1 is easy to reach its critical value before τmax reaches its critical value and thus results in Mode I fracture. The fracture trajectory is helicoid and the normal direction of tangential plane with the fractured helicoid is along the predicted direction of the maximum principal stress at the notch tip. It is further proved that Mode Ⅰ instead of Mode Ⅲ fracture occurs in brittle rock under Mode Ⅲ loading. The fracture mode depending on the fracture mechanism must be distinguished from the loading form.

Asymmetrical dynamic propagation problems on mode III interface crack

By the application of the theory of complex functions, asymmetrical dynamic propagation problems on mode Ⅲ interface crack are studied. The universal representations of analytical solutions are obtained by the approaches of self-similar function. The problems researched can be facilely transformed into Riemann-Hilbert problems and analytical solution to an asymmetrical propagation crack under the condition of point loads and unit-step loads, respectively, is acquired. After those solutions were used by superposition theorem, the solutions of arbitrarily complex problems could be attained.

MESHLESS METHOD FOR 2D MIXED-MODE CRACK PROPAGATION BASED ON VORONOI CELL

A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path.For each increment of the crack extension,the meshless method is applied to carry out a stress analysis of the cracked structure.The J-integral,which can be decomposed into mode Ⅰ and mode Ⅱ for mixed-mode crack,is used for the evaluation of the stress intensity factors(SIFs).The crack-propagation direction,predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method.The meshless results agree well with the experimental ones,which validates the accuracy and efficiency of the method.

FATIGUE GROWTH MODELING OF MIXED-MODE CRACK IN PLANE ELASTIC MEDIA

This paper presents an extension of a displacement discontinuity method with cracktip elements （a boundary element method） proposed by the author for fatigue crack growth analysis in plane elastic media under mixed-mode conditions. The boundary element method consists of the non-singular displacement discontinuity elements presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or right crack-tip element is placed locally at the corresponding left or right crack tip on top of the non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the maximum circumferential stress criterion. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the numerical approach. Crack growth is modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characteristics of some related elements are adjusted according to the manner in which the boundary element method is implemented. As an example, the fatigue growth process of cracks emanating from a circular hole in a plane elastic plate is simulated using the numerical simulation approach.

The effect of fractional thermoelasticity on a two-dimensional problem of a mode I crack in a rotating fiber-reinforced thermoelastic medium

This article is concerned with the effect of rotation on the general model of the equations of the generalized thermoe- lasticity for a homogeneous isotropic elastic half-space solid, whose surface is subjected to a Mode-I crack problem. The fractional order theory of thermoelasticity is used to obtain the analytical solutions for displacement components, force stresses, and temperature. The boundary of the crack is subjected to a prescribed stress distribution and temperature. The normal mode analysis technique is used to solve the resulting non-dimensional coupled governing equations of the problem. The variations of the considered variables with the horizontal distance are illustrated graphically. Some particular cases are also discussed in the context of the problem. Effects of the fractional parameter, reinforcement, and rotation on the varia- tions of different field quantities inside the elastic medium are analyzed graphically. Comparisons are made between the results in the presence and those in the absence of fiber-reinforcing, rotating and fractional parameters.

Development of X-FEM methodology and study on mixed-mode crack propagation

The extended finite element method（X-FEM） is a novel numerical methodology with a great potential for using in multi-scale computation and multi-phase coupling problems. The algorithm is discussed and a program is developed based on X-FEM for simulating mixed-mode crack propagation. The maximum circumferential stress criterion and interaction integral are deduced. Some numerical results are compared with the experimental data to prove the capability and efficiency of the algorithm and the program. Numerical analyses of sub-interfacial crack growth in bi-materials give a clear description of the effiect on fracture made by interface and loading condition.