Crack propagation in an elastic solid subjected to general loading— IV. Obliquely incident stress pulse

The stress intensity factors of a half-plane crack extending nonuniformly in an isotropic elastic solid subjected to stress wave loading are considered. A plane stress pulse is obliquely incident on the crack, and the wavefront strikes the crack at some initial time. At some arbitrary later time, the crack begins to extend at a nonuniform rate. It is found that the mode I and mode II stress intensity factors each have the form of the product of a universal function of instantaneous cracktip speed with the stress intensity factor for an equivalent stationary crack. An energy-rate balance fracture criterion is applied to obtain an equation of motion for the crack tip and to determine the actual delay time between the arrival of the incident wave and the onset of fracture as a function of angle of incidence of the loading wave.     

Crack propagation in an elastic solid subjected to general loading—III. Stress wave loading

The stress intensity factor of a half-plane crack extending non-uniformly in an isotropic elastic solid subjected to stress wave loading is determined. A plane stress pulse strikes the crack at time t = 0, the wavefront being parallel to the plane of the crack. At some arbitrary later time t = τ, the crack begins to extend at a non-uniform rate. It is found that the stress intensity factor is a universal function of instantaneous crack-tip velocity times the stress intensity factor for an equivalent stationary crack. An energy rate balance fracture criterion is applied to obtain an equation of motion for the crack tip. The delay time between the arrival of the incident pulse and the onset of fracture is also calculated for this fracture criterion.     

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113–2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent ?1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.     

Stress intensity factors for curved and kinked cracks in plane extension

A singular integral equation containing the crack opening displacement (COD) is developed for solving plane elasticity problems. The crack may contain any number of kinks at different intervals and orientations, such as a saw-tooth shape. Cracks in the form of a sine wave can also be treated. The crack tip stress intensity factors are evaluated for a variety of crack shapes and the results are displayed graphically. The distance between the crack tips is found to be a dominant factor on the crack tip stress intensity while the angle between the tangent to the crack tip and load direction determines the proportion of Mode I and II stress intensity factors.     

Interaction of general plane P wave and cylindrical inclusion partially debonded from its viscoelastic matrix

The interaction of a general plane P wave and an elastic cylindrical inclusion of infinite length partially debonded from its surrounding viscoelastic matrix of infinite extension is investigated. The debonded region is modeled as an arc-shaped interface crack between inclusion and matrix with non-contacting faces. With wave functions expansion and singular integral equation technique, the interaction problem is reduced to a set of simultaneous singular integral equations of crack dislocation density function. By analysis of the fundamental solution of the singular integral equation, it is found that dynamic stress field at the crack tip is oscillatory singular, which is related to the frequency of incident wave. The singular integral equations are solved numerically, and the crack open displacement and dynamic stress intensity factor are evaluated for various incident angles and frequencies. The project supported by the National Natural Science Foundation of China (19872002) and Climbing Foundation of Northern Jiaotong University     

The dynamic behaviour of an oblique edge-crack under impact loading

When a stress wave (tensile or compressive) impinges on a crack existing in an elastic medium, reflection, refraction and diffraction-phenomena take place. A result of diffraction is the loading of the crack. While compressive stress-waves do not create any stress concentration at the tip of an existing crack, tensile stress-waves develop stresses at the tip which may cause a propagation of the crack. If the tensile pulse is weak the crack may propagate by steps under the action only of successive tensile stress-pulses, whereas intermediate compressive-stress pulses do not have any influence.A complete study of the phenomena of incubation, initiation and propagation of cracks in thin plates, when they are subjected to a compressive pulse, which is subsequently reflected from the free boundaries of the plate and changed to complicated wave-trains, was undertaken in this paper, based on the method of caustics. Interesting results that were derived from this experimental study are presented.     

不同节理角度和地应力条件下岩石双孔爆破的数值模拟

在采矿工程、地下交通工程和水利水电工程等岩体开挖工程中,爆破法仍是一种主要的破岩方法. 实际天然岩体中存在的裂隙、节理等不连续面和所处的地应力场环境,都会对爆破荷载传播过程和最终的破碎效果产生重要影响. 本文把岩石爆破视为爆炸应力波动态作用和爆生气体压力准静态作用两个独立的先后作用过程,同时考虑初始应力场的静态作用,建立岩石爆破的力学模型. 基于渗流方程描述爆生气体在爆生裂纹中的传播过程,进而基于流固耦合理论实现爆生气体准静态压力对裂纹尖端的力学作用. 爆炸应力波主要在孔周引起压碎区和径向微裂纹区,随后爆生气体压力楔入径向裂纹尖端,促使裂纹进一步扩展,形成径向主裂纹. 数值模型可以再现孔周压碎区、径向微裂纹区、径向主裂纹区萌生、扩展的完整演化过程. 针对不同节理角度和地应力条件下岩石双孔爆破过程的数值模拟结果表明,初始地应力场的压应力作用不利于爆生裂纹的萌生与扩展,但节理的存在对裂纹的扩展具有明显的导向和促进作用,有利于爆生裂纹沿节理面方向的扩展.      

Concentrated impact loading on one face of a semi-infinite crack in an anisotropic elastic material

The response of an unbounded anisotropic elastic body containing a semi-infinite crack subjected to a concentrated impact force on one of the crack faces is studied. An exact solution of the dynamic stress intensity factors is obtained from a linear superposition of the solution of Lamb’s problem and a solution of a dislocation emitting from the crack tip. The stress intensity factors exhibit square-root singularity upon the arrival of the Rayleigh wave at the crack tip. As the Rayleigh wave passes through the crack tip, the stress intensity factors either instantaneously assume the static values or gradually approach to zero. Several numerical examples are given for isotropic, cubic and orthotropic materials.     

Rapid extension of a penny-shaped crack under torsion

The rapid extension of a penny-shaped crack under torsion is investigated. Both dynamic and quasi-static loading is considered. The wave motion is analyzed through a Green's function technique which leads to an integral equation for the stress field around the crack. Asymptotic expansions for the stress intensity and displacement rate intensity functions which are valid for a small time are obtained for the two types of loading. The propagation of the crack is analyzed through the balance of rates of energy criterion.     

各向异性介质中SH波引起的裂纹扩展

本文利用Green函数法,求解各向异性介质中半无限长裂纹在SH波作用下,以任意速度扩展的问题。首先,利用Laplace变换和Cagniard-de Hoop反演法求解各向异性介质中反平面问题的Green函数,并利用它建立了求解裂纹扩展问题的积分方程。因为方程为Abel型的,所以可得到在SH波作用下,半无限长裂纹扩展问题的解析解。还可求得裂纹端点附近的应力和裂纹表面上位移的表达式。并对裂纹端点附近的奇异性进行讨论。最后讨论了裂纹尖端附近任一点的能量关系。并应用Griffith的能量准则,对裂纹扩展规律进行了讨论。     

横观各向同性三维热弹性力学通解及其势理论法

通过引入两个位移函数，对用位移表达的运动平衡方程作了简化．利用算子理论，严格地导出了横观各向同性非耦合热弹性动力学问题的通解．对于静力学问题，通解的形式可进一步简化成用4个准调和函数来表示．具体考察了横观各向同性体内平面裂纹上下表面有对称分布温度作用的问题，推广了势理论方法，导出了一个积分方程和一个微分-积分方程．针对币状裂纹表面受均布温度作用情形，给出了具体的解。     

弹性波与单侧界裂纹相互作用问题的边界元法

措助边界元法设计了一种迭代修正方法来求解单侧界面裂纹模型与弹性波的相互作用问题,作为对算法的检验,用这种方法我们具体地分析了平面简谐弹性波对一个则界面裂纹的入射,给出了裂纹面的接触形态及应力场。     

Dynamic crack growth behaviour during explosive treatment

On the basis of mixed mode fracture as well as dynamic fracture mechanics crack growth behaviour related to the stress wave during the explosive treatment is studied experimentally and theoretically. The effects of the incident wave, the reflected wave and the stress wave duration on the crack growth behaviour are discussed. A method for restricting crack growth by using an energy absorber has been tested and recommended.     

Stress-strain field near crack tip and calculation of critical stress of crack propagation in a pure bending beam of rectangular section with one-sided mode i crack

This paper studies the stress-strain field near crack tip in a pure bending beam of rectangular section with one-sided mode I crack by the analytic method of Ref. [1], then it gives the stress and strain components at the crack tip when the crack propagates and further it obtains the formulas of calculating the elastic deformed area width, the deformed intensity area width and the equation groups of calculating the critical stress of crack propagation, last the equation group of calculating critical stress of crack propagation is verified by calculating instance. The maximum error is 0.18%. First Received May 7, 1994.     

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.     

Mode-III crack kinking under stress-wave loading

A horizontally polarized step-stress wave is incident on a semi-infinite crack in an elastic solid. At the instant that the crack tip is struck, the crack starts to propagate in the forward direction, but under an angle κπ with the plane of the original crack. In this paper a self-similar solution is obtained for the particle velocity of the diffracted cylindrical wave field. The use of Chaplygin's transformation reduces the problem to the solution of Laplace's equation in a semi-infinite strip containing a slit. The Schwarz-Christoffel transformation is employed to map the semi-infinite strip on a half-plane. An analytic function in the half-plane which satisfies appropriate conditions along the real axis, can subsequently be constructed. The Mode-III stress-intensity factor at the tip of the kinked crack has been computed for angles of incidence varying from normal to grazing incidence, for angles of crack kinking defined by -0.5?κ?0.5, and for arbitrary subsonic crack tip speeds.     

Stress intensity factor determination error by the method of caustics

The caustics created at a crack tip under biaxial loading are examined using both the closed form and the approximate solution. Particular attention is given to the caustics and their initial curves created by rays reflected from the rear and front surfaces of the specimen in addition to those in the transverse direction of the cracked specimen. These quantities can differ by more than 30% depending on whether the exact or the approximate solution is used for determining the stress intensity factors which can deviate as much as 80%. The crack length also affects the outcome. Larger errors are committed by the approximate solution for caustics that are closer to the crack tip. The advantage of using the normalized distance from the crack tip is discussed in contrast to using the absolute distance. Incorrect experimental measurements are also discussed in relation to the resulting error for evaluating the crack tip stress intensity factor.     

爆炸处理时裂纹动态扩展行为的研究

本文研究爆炸处理时裂纹在应力波作用下的动态扩展行为,分别讨论了入射波、反射波和应力波持续时间的作用,并根据线弹性断裂力学复合型准则和动态断裂准则对裂纹的扩展进行了分析,作者建议用吸能块的办法控制爆炸处理时裂纹的扩展,并作了实验验证。     

Sudden deceleration or acceleration of an intersonic shear crack

Sudden jumps in the crack tip velocity were revealed by numerical simulation (in both continuum/cohesive element and molecular dynamics approaches) and experiments for rapid shear cracking. The cracking velocity may accelerate from a sub-Rayleigh speed to the intersonic range, or from an intersonic speed to a higher one, when the reflected impact wave reloads the crack tip. On the other hand, the cracking velocity may decelerate from an intersonic speed to a lower one or recede to the sub-Rayleigh range when the fracture driving force declines. The velocity change encountered during intersonic cracking plays a different role from that in the acceleration or deceleration of a subsonic crack. A crack propagating at an intersonic speed would leave a shear wave trailing behind. When the crack decelerates or accelerates, the effect of the trailing wave will lead to a transition period from one steady-state solution of crack tip singularity to another. This investigation aims at quantifying these processes. The full field solution of an intersonic mode II crack whose speed changed suddenly from one velocity (intersonic or subsonic) to another (intersonic or subsonic) is given in closed form. The solution is facilitated via superposing a series of propagating crack problems that are loaded by dislocations to seal the unwanted crack-face sliding or by concentrated forces moving at various speeds to negate the crack-face traction. In contrast to the subsonic solution, the results in the intersonic case indicate that the elastic fields around the crack tip depend on the deceleration or acceleration history that is traced back over a long time. Singularity matching dictates the jump that may actually take place.