On the temperature rise at the tip of a fast running crack2

The high energy concentration at the tip of a running crack leads to irreversible deformations, and a great amount of the deformation energy is set free as heat. Assuming that this moving heat source is of circular shape, the temperature distribution around the crack tip has been calculated. The temperatures are dependent on the radius of the heat source and the crack velocity. Some examples for the material glass are given. The very high temperatures computed lead to the supposition that the observed light emission during fast fracture is of thermal origin.     

Heating during shearing and opening dominated dynamic fracture of polymers

The heat generated from dissipative mechanisms during shearing and opening dominated dynamic fracture of polymethyl methacrylate and polycarbonate was measured with a remote sensing technique that utilizes the detection of infrared radiation. Significant heating was detected for both materials and both modes of fracture. In the shear dominated experiments, the temperature increase at the crack tip in polymethyl methacrylate was 85 K, the approximate increase necessary to reach the glass transition temperature. An adiabatic shear band followed by a dynamically propagating crack were observed during the shear dominated experiments using polycarbonate. The recorded shear band temperature increase was 45 K. This was followed by an additional 100 K temperature increase from the ensuing crack, raising the temperature above glass transition. The maximum temperature increase recorded for the opening mode experiments was 55 K for polymethyl methacrylate and 105 K for polycarbonate. The results of this study show that temperature effects are significant during the dynamic fracture of polymers. The effects are especially important in the shear dominated case where local temperatures approach or exceed the polymer glass transition temperature.     

Dynamic self-similar debonding of interface at very high velocity: Anti-plane case

This paper analyzes the anti-plane problem of dynamic self-similar debonding of interface at very high velocity. The debonding is modeled as an interface crack propagating self-similarly from zero-length. The extending speed is assumed to be transonic or supersonic. We first consider the dynamic debonding under moving concentrated loads. The moving dislocation model of self-similar propagation of an interface crack is used to formulate the problem to a singular integral equation which is solved analytically. The singularity of stresses near the crack tip is discussed and the dynamic stress intensity factors are presented. Finally the solution of dynamic debonding underx ²-type loads is obtained by using the superposition method.     

The deformation and fracture of balloons

Inflation of balloons provides a straightforward way of achieving large biaxial deformations. Previous studies have shown that when a balloon bursts, crack propagation occurs at very high speed – much higher than would be expected from the low strain modulus and elastic wave velocity of the rubber. The present paper is concerned with studies of the deformation and fracture of cylindrical balloons. On inflation, the deformations of such a balloon pass through an unstable region but subsequently increase monotonically with pressure. In this relatively high pressure region, the ratio of the longitudinal and circumferential extension ratios is broadly in accord with expectations from high-strain elasticity theory when the ratio of the corresponding stresses is taken into account. On bursting, crack speeds up to around 300 m/s in this region. It is shown that these speeds are in accord with large increase in incremental moduli for the highly-strained rubber. Marked changes in crack tip profile observed at very high crack speeds are consistent with control of the rate of growth by inertia rather than by the viscoelastic properties of the rubber (as is believed to be the case at lower speeds). Consistent with this, various elastomers having different glass transition temperatures show similar crack growth behaviour in the very high speed region.     

The determination of dynamic fracture toughness of AISI 4340 steel by the shadow spot method

Dynamic crack propagation experiments have been performed using wedge loaded double cantilever beam specimens of an austenitized, quenched and tempered 4340 steel. Measurements of the dynamic stress intensity factor have been made by means of the optical method of caustics. The interpretation of experimental data, obtained from the shadow spot patterns photographed with a Cranz-Schardin high speed camera, is based on an elastodynamic analysis. The instantaneous value of the dynamic stress intensity factor K^d_I is obtained as a function of crack tip velocity. Finally, the interaction of reflected shear and Rayleigh waves with the moving crack tip stress field is considered.     

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.     

Investigation of the dynamics of the development of cracks by the method of photoelasticity

An experimental investigation was made of the process of the development of cracks in two materials: polymethyl methacrylate and a polymerized epoxy resin. For these materials, determinations were made of their optical and mechanical characteristics, such as density, the speed of sound, the critical value of the coefficient of singularity of the stress field at the tip of cracks. As the dynamic characteristic of the process of the development of cracks, an investigation is made of the dependences connecting the coefficient of the singularity of the stress field at the tip of a moving crack and the rate of development of the crack. The question of the interaction between two cracks moving in a previously elongated sample is discussed.     

Crack propagation in wedged double cantilevered beam specimens

A closed form analytical solution of crack propagation in double cantilevered beam specimens opened at a constant rate has been found. Hamilton's principle for non-conservative systems was applied to describe the crack motion, under the assumption of a Bernoulli-Euler beam. The criterion of crack propagation is a critical bending moment at the crack tip. The calculations of beam motion take into account wave effects in the Bernoulli-Euler theory of elastic beams. The beam shape during the crack motion is found with a similarity transformation and expressed by Fresnel integrals. The boundary conditions satisfied are the fixed ones of zero bending moment and constant beam opening rate at the load end of the specimen and the moving ones of zero deflection and zero slope of the deflected beam at the tip of the moving crack. The fracture represents a moving critical bending moment. The analytical results show that the specific fracture surface energy is a unique function of the ratio of the crack length squared to the time subsequent to loading and this is computed from the recorded time-dependence of the crack length.     

Constant moving crack in a magnetoelectroelastic material under anti-plane shear loading

Analytical solutions for an anti-plane Griffith moving crack inside an infinite magnetoelectroelastic medium under the conditions of permeable crack faces are formulated using integral transform method. The far-field anti-plane mechanical shear and in-plane electrical and magnetic loadings are applied to the magnetoelectroelastic material. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived. Field intensity factors for magnetoelectroelastic material are obtained. The stresses, electric displacements and magnetic inductions at the crack tip show inverse square root singularities. The moving speed of the crack have influence on the dynamic electric displacement intensity factor (DEDIF) and the dynamic magnetic induction intensity factor (DMIIF), while the dynamic stress intensity factor (DSIF) does not depend on the velocity of the moving crack. When the crack is moving at very lower or very higher speeds, the crack will propagate along its original plane; while in the range of M_c1 < M < M_c2, the propagation of the crack possibly brings about the branch phenomena in magnetoelectroelastic media.     

On a moving interface crack with a contact zone in a piezoelectric bimaterial

An inplane problem for a crack moving with constant subsonic speed along the interface of two piezoelectric materials is considered. A mechanically frictionless and electrically permeable contact zone is assumed at the right crack tip whilst for the open part of the crack both electrically permeable and electrically insulated conditions are considered. In the first case a moving concentrated loading is prescribed at the crack faces and in the second case an additional electrical charge at the crack faces is prescribed as well. The main attention is devoted to electrically permeable crack faces. Introducing a moving coordinate system at the leading crack tip the corresponding inhomogeneous combined Dirichlet–Riemann problem is formulated and solved exactly for this case. All electromechanical characteristics at the interface are presented in a closed form for arbitrary contact zone lengths, and further, the transcendental equation for the determination of the real contact zone length is derived. As a particular case of the obtained solution a semi-infinite crack with a contact zone is considered. The numerical analysis performed for a certain piezoelectric bimaterial showed an essential increase of the contact zone length and the associated stress intensity factor especially for the near-critical speed region. Similar investigations have been performed for an electrically insulated crack and the same behavior of the above mentioned parameters is observed.     

Determining stress intensity factors in orthotropic composites from far-field measured temperatures

We demonstrate the ability to determine stress intensity factors in orthotropic materials directly from measured temperatures away from the crack and using far-field expressions for the stresses. This is advantageous, recognizing that recorded thermoelastic data can be very unreliable near the tip of a crack. In addition to singular terms that govern in the immediate vicinity of the crack tip, the present series expressions for the stresses contain higher-order finite terms. Little measured input information is needed and data acquisition positions can be selected largely at the user's discretion.     

含裂纹体的单位分解扩展无网格方法

不连续体的数值模拟尤其是动态裂纹的追踪问题一直是工程界研究的热点和难点问题。无网格方法仅仅需要结点信息,非常适合于求解这类问题。基于单位分解思想,在移动最小二乘近似函数(MLS)中根据裂纹面的不连续位移增加一个Heaviside函数,在裂尖则增加四个扩展函数描述渐进裂纹位移场;应用Galerkin方法推导了平衡方程的离散线性方程,并给出了求解裂纹问题应力强度因子的计算公式。与其他类型的扩展无网格相比,在裂尖处近似函数不需要使用可视准则,很容易生成r1/2奇异;另一个优势是影响域并没有因为裂纹的存在而改变,不会降低方程的稀疏性,求解效率较高。数值算例表明,该方法能方便有效地模拟不连续问题,具有十分广阔的应用空间。     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

PERFECTLY PLASTIC FIELDS AT A RAPIDLY PROPAGATING PLANE-STRESS CRACK TIP

Under the condition that all the perfectly plastic stress components at a crack tiP arethe functions ofθonly,making use of the Mises yield condition,steady-state movingequations and elastic perfectly-plastic constitutive equations,we derive the generallyanalytical expressions of perfectly plastic fields at a rapidly propagating plane-stress cracktip.Applying these generally analytical expressions to the concrete crack,we obtain theanalytical expressions of perfectly plastic fields at the rapidly propagating tips of,modesⅠandⅡplane-stress cracks.     

断裂力学判据的评述

从Inglis 和Griffith 的著名论文到Irwin 和Rice 等的奠基性贡献,对断裂力学中的线弹性断裂力学的K判据,界面断裂力学的G判据,和弹塑性断裂力学的J 判据作了扼要的综述. 介绍了在界面断裂力学G判据的基础上提出的界面断裂力学的K判据,以说明断裂力学的判据存在改进的可能性. 在综述中归纳出断裂力学判据中目前还没有较好解决的几个问题. 在总结以往断裂力学研究经验的基础上,指出裂纹端应力奇异性的源是对断裂力学判据存在的问题作进一步研究的切入点. 探讨了裂纹端应变间断的奇点是裂纹端应力奇异性的源的问题,从而对裂纹端应力强度因子的物理意义进行了讨论. 最后,阐述了进行可靠的裂纹端应力场的弹塑性分析是改进弹塑性断裂力学判据的关键,而进行可靠的裂纹端应力场的弹塑性分析的前提是要通过裂纹端应力奇异性的源的研究来获得作用在裂纹端的造成裂纹端应变间断的有限值应力.      

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.     

功能梯度压电板条中电绝缘型运动裂纹的电弹性场

基于三维弹性理论和压电理论，对材料系数按指数函数规律分布的功能梯度压电板条中的反平面运动裂纹问题进行了求解。利用Fourier积分变换方法将电绝缘型运动裂纹问题化为对偶积分方程，并进一步归结为易于求解的第二类Fredholm积分方程。通过渐近分析，获得了裂纹尖端应力、应变、电位移和电场的解析解，给出了裂纹尖端场各个变量的角分布函数，并求得了裂纹尖端场的强度因子，分析了压电材料物性梯度参数、几何尺寸及裂纹运动速度对它们的影响。结果表明，对于电绝缘型裂纹，功能梯度压电板条中运动裂纹尖端附近的各个场变量都具有-1／2阶的奇异性；当裂纹运动速度增大时，裂纹扩展的方向会偏离裂纹面。     

Plastic flow at a crack tip. Energy failure criterion and its relationship to the J integral

Plastic flow at the tip of a crack moving in an elastic body is considered using the theory of an ideally rigid-plastic body. The material at the crack tip is treated as a body consisting of an elastic outer region and a rigid-plastic inner region. It is shown that this representation is energetically justified for small plastic regions. The distribution of the specific dissipation of the work of internal forces and deformations along the particle trajectory at the crack tip is obtained. The relationship between the specific dissipation of the work of internal forces and the J integral under plane-strain conditions is established.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

Theoretical development and experimental validation of a thermally dissipative cohesive zone model for dynamic fracture of amorphous polymers

A thermally dissipative cohesive zone model is developed for predicting the temperature increase at the tip of a crack propagating dynamically in a nominally brittle material exhibiting a cohesive-type failure such as crazing. The model assumes that fracture energy supplied to the crack tip region that is in excess of that needed for the creation of new free surfaces during crack advance is converted to heat within the cohesive zone. Bulk dissipation mechanisms, such as plasticity, are not accounted for. Several cohesive traction laws are examined, and the model is then used to make predictions of crack tip heating at various crack propagation speeds in the nominally brittle amorphous polymer PMMA, observed to fail by a crazing-type mechanism. The heating predictions are compared to experimental data where the temperature field surrounding a high speed crack in PMMA was measured. Measurements are made in real time using a multi-point high speed HgCdTe infrared radiation detector array. At the same time as temperature, simultaneous measurement of fracture energy is made by a strain gauge technique, and crack tip speed is monitored through a resistance ladder method. Material strength can be estimated through uniaxial tension tests, thus minimizing the need for parameter fitting in the stress-opening traction law. Excellent agreement between experiments and theory is found for two of the cohesive traction law temperature predictions, but only for the case where a single craze is active during the dynamic fracture of PMMA, i.e. crack tip speed up to approximately 0.2c_R. For higher speed fracture where subsurface damage becomes prominent, the line dissipation model of a cohesive zone is inadequate, and a distributed damage model is needed.